The growth and survival of myeloma cells are greatly affected by their surrounding microenvironment. To understand the molecular mechanism and the impact of stiffness on the fate of myeloma initiating cells (MICs), we develop a systems biological model to reveal the dynamic regulations by integrating reverse phase protein array data and the stiffness associated pathway.