A hybrid three-scale model of tumor growth
Abstract
Cancer results from a complex interplay of different biological, chemical, and physical phenomena that span a wide range of time and length scales. Computational modeling may help to unfold the role of multiple evolving factors that exist and interact in the tumor microenvironment. Understanding these complex multiscale interactions is a crucial step toward predicting cancer growth and in developing effective therapies. We integrate different modeling approaches in a multiscale, avascular, hybrid tumor growth model encompassing tissue, cell, and sub-cell scales. At the tissue level, we consider the dispersion of nutrients and growth factors in the tumor microenvironment, which are modeled through reaction–diffusion equations. At the cell level, we use an agent-based model (ABM) to describe normal and tumor cell dynamics, with normal cells kept in homeostasis and cancer cells differentiated into quiescent, proliferative, migratory, apoptotic, hypoxic, and necrotic states. Cell movement is driven by the balance of a variety of forces according to Newton’s second law, including those related to growth-induced stresses. Phenotypic transitions are defined by specific rule of behaviors that depend on microenvironment stimuli. We integrate in each cell/agent a branch of the epidermal growth factor receptor (EGFR) pathway. This pathway is modeled by a system of coupled nonlinear differential equations involving the mass laws of 20 molecules. The rates of change in the concentration of some key molecules trigger proliferation or migration advantage response. The bridge between cell and tissue scales is built through the reaction and source terms of the partial differential equations. Our hybrid model is built in a modular way, enabling the investigation of the role of different mechanisms at multiple scales on tumor progression. This strategy allows representing both the collective behavior due to cell assembly as well as microscopic intracellular phenomena described by signal transduction pathways. Here, we investigate the impact of some mechanisms associated with sustained proliferation on cancer progression. Speci- fically, we focus on the intracellular proliferation/migration-advantage-response driven by the EGFR pathway and on proliferation inhibition due to accumulation of growth-induced stresses. Simulations demonstrate that the model can adequately describe some complex mechanisms of tumor dynamics, including growth arrest in avascular tumors. Both the sub-cell model and growth-induced stresses give rise to heterogeneity in the tumor expansion and a rich variety of tumor behaviors.
Communicated by N. Bellomo
