Tumor proliferation and invasion are intrinsically coupled and unraveled through tunable spheroid and physics-based models.

Proliferation and invasion are two key drivers of tumor growth that are traditionally considered independent multicellular processes. However, these processes are intrinsically coupled through a maximum carrying capacity, i.e., the maximum spatial cell concentration supported by the tumor volume, total cell count, nutrient access, and mechanical properties of the tissue stroma. We explored this coupling of proliferation and invasion through in vitro and in silico methods where we modulated the mechanical properties of the tumor and the surrounding extracellular matrix. E-cadherin expression and stromal collagen concentration were manipulated in a tunable breast cancer spheroid to determine the overall impacts of these tumor variables on net tumor proliferation and continuum invasion. We integrated these results into a mixed-constitutive formulation to computationally delineate the influences of cellular and extracellular adhesion, stiffness, and mechanical properties of the extracellular matrix on net proliferation and continuum invasion. This framework integrates biological in vitro data into concise computational models of invasion and proliferation to provide more detailed physical insights into the coupling of these key tumor processes and tumor growth. STATEMENT OF SIGNIFICANCE: Tumor growth involves expansion into the collagen-rich stroma through intrinsic coupling of proliferation and invasion within the tumor continuum. These processes are regulated by a maximum carrying capacity that is determined by the total cell count, tumor volume, nutrient access, and mechanical properties of the surrounding stroma. The influences of biomechanical parameters (i.e., stiffness, cell elongation, net proliferation rate and cell-ECM friction) on tumor proliferation or invasion cannot be unraveled using experimental methods alone. By pairing a tunable spheroid system with computational modeling, we delineated the interdependencies of each system parameter on tumor proliferation and continuum invasion, and established a concise computational framework for studying tumor mechanobiology.