Front instabilities and invasiveness of simulated avascular tumors

Improving cancer treatments via dynamical biophysical models

Despite significant advances in oncological research, cancer nowadays remains one of the main causes of mortality and morbidity worldwide. New treatment techniques, as a rule, have limited efficacy, target only a narrow range of oncological diseases, and have limited availability to the general public due their high cost. An important goal in oncology is thus the modification of the types of antitumor therapy and their combinations, that are already introduced into clinical practice, with the goal of increasing the overall treatment efficacy. One option to achieve this goal is optimization of the schedules of drugs administration or performing other medical actions. Several factors complicate such tasks: the adverse effects of treatments on healthy cell populations, which must be kept tolerable; the emergence of drug resistance due to the intrinsic plasticity of heterogeneous cancer cell populations; the interplay between different types of therapies administered simultaneously. Mathematical modeling, in which a tumor and its microenvironment are considered as a single complex system, can address this complexity and can indicate potentially effective protocols, that would require experimental verification. In this review, we consider classical methods, current trends and future prospects in the field of mathematical modeling of tumor growth and treatment. In particular, methods of treatment optimization are discussed with several examples of specific problems related to different types of treatment.

Unification of aggregate growth models by emergence from cellular and intracellular mechanisms

Multicellular aggregate growth is regulated by nutrient availability and removal of metabolites, but the specifics of growth dynamics are dependent on cell type and environment. Classical models of growth are based on differential equations. While in some cases these classical models match experimental observations, they can only predict growth of a limited number of cell types and so can only be selectively applied. Currently, no classical model provides a general mathematical representation of growth for any cell type and environment. This discrepancy limits their range of applications, which a general modelling framework can enhance. In this work, a hybrid cellular Potts model is used to explain the discrepancy between classical models as emergent behaviours from the same mathematical system. Intracellular processes are described using probability distributions of local chemical conditions for proliferation and death and simulated. By fitting simulation results to a generalization of the classical models, their emergence is demonstrated. Parameter variations elucidate how aggregate growth may behave like one classical growth model or another. Three classical growth model fits were tested, and emergence of the Gompertz equation was demonstrated. Effects of shape changes are demonstrated, which are significant for final aggregate size and growth rate, and occur stochastically.

A multiscale in silico model of endothelial to mesenchymal transformation in a tumor microenvironment

Endothelial to mesenchymal transformation (EndMT) is a process in which endothelial cells gain a mesenchymal-like phenotype in response to mechanobiological signals that results in the remodeling or repair of underlying tissue. While initially associated with embryonic development, this process has since been shown to occur in adult tissue remodeling including wound healing, fibrosis, and cancer. In an attempt to understand the role of EndMT in cancer progression and metastasis, we present a multiscale, three-dimensional, in silico model. The model couples tissue level phenomena such as extracellular matrix remodeling, cellular level phenomena such as migration and proliferation, and chemical transport in the tumor microenvironment to mimic in vitro tissue models of the cancer microenvironment. The model is used to study the presence of EndMT-derived activated fibroblasts (EDAFs) and varying substrate stiffness on tumor cell migration and proliferation. The simulations accurately model the behavior of tumor cells under given conditions. The presence of EDAFs and/or an increase in substrate stiffness resulted in an increase in tumor cell activity. This model lays the foundation of further studies of EDAFs in a tumor microenvironment on a cellular and subcellular physiological level.

A Mathematical Framework for Modelling the Metastatic Spread of Cancer

Cancer is a complex disease that starts with mutations of key genes in one cell or a small group of cells at a primary site in the body. If these cancer cells continue to grow successfully and, at some later stage, invade the surrounding tissue and acquire a vascular network, they can spread to distant secondary sites in the body. This process, known as metastatic spread, is responsible for around 90% of deaths from cancer and is one of the so-called hallmarks of cancer. To shed light on the metastatic process, we present a mathematical modelling framework that captures for the first time the interconnected processes of invasion and metastatic spread of individual cancer cells in a spatially explicit manner-a multigrid, hybrid, individual-based approach. This framework accounts for the spatiotemporal evolution of mesenchymal- and epithelial-like cancer cells, membrane-type-1 matrix metalloproteinase (MT1-MMP) and the diffusible matrix metalloproteinase-2 (MMP-2), and for their interactions with the extracellular matrix. Using computational simulations, we demonstrate that our model captures all the key steps of the invasion-metastasis cascade, i.e. invasion by both heterogeneous cancer cell clusters and by single mesenchymal-like cancer cells; intravasation of these clusters and single cells both via active mechanisms mediated by matrix-degrading enzymes (MDEs) and via passive shedding; circulation of cancer cell clusters and single cancer cells in the vasculature with the associated risk of cell death and disaggregation of clusters; extravasation of clusters and single cells; and metastatic growth at distant secondary sites in the body. By faithfully reproducing experimental results, our simulations support the evidence-based hypothesis that the membrane-bound MT1-MMP is the main driver of invasive spread rather than diffusible MDEs such as MMP-2.

A multiscale model for heterogeneous tumor spheroid in vitro

In this paper, a novel multiscale method is proposed for the study of heterogeneous tumor spheroid growth in vitro. The entire tumor spheroid is described by an ellipsoid-based model while nutrient and other environmental factors are treated as continua. The ellipsoid-based discrete component is capable of incorporating mechanical effects and deformability, while keeping a minimum set of free variables to describe complex shape variations. Moreover, our purely cell-based description of tumor avoids the complex mutual conversion between a cell-based model and continuum model within a tumor, such as force and mass transformation. This advantage makes it highly suitable for the study of tumor spheroids in vitro whose size are normally less than 800 μm in diameter. In addition, our numerical scheme provides two computational options depending on tumor size. For a small or medium tumor spheroid, a three-dimensional (3D) numerical model can be directly applied. For a large spheroid, we suggest the use of a 3D-adapted 2D cross section configuration, which has not yet been explored in the literature, as an alternative for the theoretical investigation to bridge the gap between the 2D and 3D models. Our model and its implementations have been validated and applied to various studies given in the paper. The simulation results fit corresponding in vitro experimental observations very well.

The influence of soluble fragments of extracellular matrix (ECM) on tumor growth and morphology

A major challenge in matrix-metalloproteinase (MMP) target validation and MMP-inhibitor-drug development for anti-cancer clinical trials is to better understand their complex roles (often competing with each other) in tumor progression. While there is extensive research on the growth-promoting effects of MMPs, the growth-inhibiting effects of MMPs has not been investigated thoroughly. So we develop a continuum model of tumor growth and invasion including chemotaxis and haptotaxis in order to examine the complex interaction between the tumor and its host microenvironment and to explore the inhibiting influence of the gradients of soluble fragments of extracellular matrix (ECM) density on tumor growth and morphology. Previously, it was shown both computationally (in one spatial dimension) and experimentally that the chemotactic pull due to soluble ECM gradients is anti-invasive, contrary to the traditional view of the role of chemotaxis in malignant invasion [1]. With two-dimensional numerical simulation and using a level set based tumor-host interface capturing method, we examine the effects of chemotaxis on the progression and morphology of a tumor growing in nutrient-rich and nutrient-poor microenvironments which was not investigated before. In particular we examine how the geometry of the growing tumor is affected when placed in different environments. We also investigate the effects of varying ECM degradation rate, the production rate of matrix degrading enzymes (MDE), and the conversion of ECM into soluble ECM. We find that chemotaxis due to ECM-fragment gradients strongly influences tumor growth and morphology, and that the instabilities caused by tumor cell proliferation and haptotactic movements can be prevented if chemotaxis is sufficiently strong. The influence of chemotaxis and the above factors on tumor growth and morphology are found to be more prominent in nutrient-poor environments than in nutrient-rich environments. So we extend our investigations of these antinvasive chemotactic influences by examining the effects of cell-cell and cell-ECM adhesion and low proliferation rate for tumors growing in low-nutrient environments. We find that as the extent of chemotaxis increases, the effects of adhesion on tumor growth and shape become negligible. Under conditions of low cell mitosis, chemotaxis may cause the tumor to shrink, as the extent of chemotaxis increases. Both stable and unstable tumor shrinkage are predicted by our model. Unexpectedly, in some cases chemotaxis may contribute toward developing instability where haptotaxis alone induces stable growth.

Tumour heterogeneity promotes collective invasion and cancer metastatic dissemination

Heterogeneity within tumour cell populations is commonly observed in most cancers. However, its impact on metastatic dissemination, one of the primary determinants of the disease prognosis, remains poorly understood. Working with a simplified numerical model of tumour spheroids, we investigated the impact of mechanical heterogeneity on the onset of tumour invasion into surrounding tissues. Our work establishes a positive link between tumour heterogeneity and metastatic dissemination, and recapitulates a number of invasion patterns identified in vivo, such as multicellular finger-like protrusions. Two complementary mechanisms are at play in heterogeneous tumours. A small proportion of stronger cells are able to initiate and lead the escape of cells, while collective effects in the bulk of the tumour provide the coordination required to sustain the invasive process through multicellular streaming. This suggests that the multicellular dynamics observed during metastasis is a generic feature of mechanically heterogeneous cell populations and might rely on a limited and generic set of attributes.

Blood vessel tortuosity selects against evolution of aggressive tumor cells in confined tissue environments: A modeling approach

Cancer is a disease of cellular regulation, often initiated by genetic mutation within cells, and leading to a heterogeneous cell population within tissues. In the competition for nutrients and growth space within the tumors the phenotype of each cell determines its success. Selection in this process is imposed by both the microenvironment (neighboring cells, extracellular matrix, and diffusing substances), and the whole of the organism through for example the blood supply. In this view, the development of tumor cells is in close interaction with their increasingly changing environment: the more cells can change, the more their environment will change. Furthermore, instabilities are also introduced on the organism level: blood supply can be blocked by increased tissue pressure or the tortuosity of the tumor-neovascular vessels. This coupling between cell, microenvironment, and organism results in behavior that is hard to predict. Here we introduce a cell-based computational model to study the effect of blood flow obstruction on the micro-evolution of cells within a cancerous tissue. We demonstrate that stages of tumor development emerge naturally, without the need for sequential mutation of specific genes. Secondly, we show that instabilities in blood supply can impact the overall development of tumors and lead to the extinction of the dominant aggressive phenotype, showing a clear distinction between the fitness at the cell level and survival of the population. This provides new insights into potential side effects of recent tumor vasculature normalization approaches.

Polyurethane Microgel Based Microtissue: Interface-Guided Assembly and Spreading

Colloidal gels are three-dimensional networks of microgel particles and can be utilized to design microtissues where the differential adhesive interactions between the particles and cells, guided by their surface energetics, are engineered to spatially assemble the cellular and colloidal components into three-dimensional microtissues. In this work we utilized a colloidal interaction approach to design cell-polyurethane (PU) microgel bimodal microtissues using endothelial cells (ECs) as a normal cell model and a nonmalignant breast cancer cell line (MCF-7) as a cancer cell model. PU microgels were developed from a library of segmental polyurethanes with poly(ethylene glycol) soft segment and aliphatic diisocyanate/l-tyrosine based chain extender as hard segment to modulate the interactions between PU colloidal particles and cells. The surface energies of the microgel particles and cells were estimated using Zisman's critical surface tension and van Oss-Good-Chaudhury theory (vOGCT) from liquid contact angle analysis. Binary interaction potentials between colloidal PU particles and cells and the ternary interaction between colloidal PU particle, cell, and collagen I/Matrigel were calculated to explain the formation of microtissues and their spreading in extraneous biomatrix respectively by using classical and extended DLVO theory (XDLVO). Furthermore, rheological analysis and in silico simulations were used to analyze the assembly and spreading of the PU microgel based microtissues. In vitro experiments showed that ECs and MCF-7 displayed more differentiated (EC spreading/MCF-7 lumen formation) character when mixed with microgel particles that were stable in aqueous medium and more undifferentiated character (EC nonspreading/MCF-7 spreading) when mixed with microgel particles unstable in aqueous medium.

The universal growth rate behavior and regime transition in adherent cell colonies

In this work, we used five cell lineages, cultivated in vitro, to show they follow a common functional form to the growth rate: a sigmoidal curve, suggesting that competition and cooperation (usual mechanisms for systems with this behavior) might be present. Both theoretical and experimental investigations, on the causes of this behavior, are challenging for the research field; since the sigmoidal form to the growth rate seems to absorb important properties of such systems, e.g., cell deformation and statistical interactions. We shed some light on this subject by showing how cell spreading affects the radius behavior of the growing colonies. Doing numerical time derivatives of the experimental data, we obtained the growth rates. Using reduced variables for the time and rates, we obtained the collapse of all colonies growth rates onto one curve with sigmoidal shape. This suggests a universal-type behavior, with regime transition related to a morphological transition of adherent cell colonies.

Emergent Stratification in Solid Tumors Selects for Reduced Cohesion of Tumor Cells: A Multi-Cell, Virtual-Tissue Model of Tumor Evolution Using CompuCell3D

Tumor cells and structure both evolve due to heritable variation of cell behaviors and selection over periods of weeks to years (somatic evolution). Micro-environmental factors exert selection pressures on tumor-cell behaviors, which influence both the rate and direction of evolution of specific behaviors, especially the development of tumor-cell aggression and resistance to chemotherapies. In this paper, we present, step-by-step, the development of a multi-cell, virtual-tissue model of tumor somatic evolution, simulated using the open-source CompuCell3D modeling environment. Our model includes essential cell behaviors, microenvironmental components and their interactions. Our model provides a platform for exploring selection pressures leading to the evolution of tumor-cell aggression, showing that emergent stratification into regions with different cell survival rates drives the evolution of less cohesive cells with lower levels of cadherins and higher levels of integrins. Such reduced cohesivity is a key hallmark in the progression of many types of solid tumors.

Bridging scales in cancer progression: mapping genotype to phenotype using neural networks

In this review we summarise our recent efforts in trying to understand the role of heterogeneity in cancer progression by using neural networks to characterise different aspects of the mapping from a cancer cells genotype and environment to its phenotype. Our central premise is that cancer is an evolving system subject to mutation and selection, and the primary conduit for these processes to occur is the cancer cell whose behaviour is regulated on multiple biological scales. The selection pressure is mainly driven by the microenvironment that the tumour is growing in and this acts directly upon the cell phenotype. In turn, the phenotype is driven by the intracellular pathways that are regulated by the genotype. Integrating all of these processes is a massive undertaking and requires bridging many biological scales (i.e. genotype, pathway, phenotype and environment) that we will only scratch the surface of in this review. We will focus on models that use neural networks as a means of connecting these different biological scales, since they allow us to easily create heterogeneity for selection to act upon and importantly this heterogeneity can be implemented at different biological scales. More specifically, we consider three different neural networks that bridge different aspects of these scales and the dialogue with the micro-environment, (i) the impact of the micro-environment on evolutionary dynamics, (ii) the mapping from genotype to phenotype under drug-induced perturbations and (iii) pathway activity in both normal and cancer cells under different micro-environmental conditions.

Lattice-based model of ductal carcinoma in situ suggests rules for breast cancer progression to an invasive state

Ductal carcinoma in situ (DCIS) is a heterogeneous group of non-invasive lesions of the breast that result from abnormal proliferation of mammary epithelial cells. Pathologists characterize DCIS by four tissue morphologies (micropapillary, cribriform, solid, and comedo), but the underlying mechanisms that distinguish the development and progression of these morphologies are not well understood. Here we explored the conditions leading to the emergence of the different morphologies of DCIS using a two-dimensional multi-cell lattice-based model that incorporates cell proliferation, apoptosis, necrosis, adhesion, and contractility. We found that the relative rates of cell proliferation and apoptosis governed which of the four morphologies emerged. High proliferation and low apoptosis favored the emergence of solid and comedo morphologies. In contrast, low proliferation and high apoptosis led to the micropapillary morphology, whereas high proliferation and high apoptosis led to the cribriform morphology. The natural progression between morphologies cannot be investigated in vivo since lesions are usually surgically removed upon detection; however, our model suggests probable transitions between these morphologies during breast cancer progression. Importantly, cribriform and comedo appear to be the ultimate morphologies of DCIS. Motivated by previous experimental studies demonstrating that tumor cells behave differently depending on where they are located within the mammary duct in vivo or in engineered tissues, we examined the effects of tissue geometry on the progression of DCIS. In agreement with our previous experimental work, we found that cells are more likely to invade from the end of ducts and that this preferential invasion is regulated by cell adhesion and contractility. This model provides additional insight into tumor cell behavior and allows the exploration of phenotypic transitions not easily monitored in vivo.

Breast cancer is a complex disease that affects women worldwide. One heterogeneous group of lesions, ductal carcinoma in situ (DCIS), often begins as a nonmalignant disease but can readily progress if left untreated. The progression of this disease is not well understood because DCIS is typically removed upon detection. Therefore, computational models might help predict whether DCIS will remain nonmalignant or progress towards invasive ductal carcinoma. Here we used a multi-cell lattice-based model to explore the relative effects of cell proliferation, death, division axis, adhesion and contractility on the development and progression of DCIS. We also examined the emergence and progression of DCIS in physiologically relevant geometries of the mammary duct. Our model suggests several plausible progressions between morphologies of DCIS, and predicts that some regions of a duct are preferential for tumor cell invasion.

Fingering in Stochastic Growth Models

Motivated by the widespread use of hybrid-discrete cellular automata in modeling cancer, two simple growth models are studied on the two dimensional lattice that incorporate a nutrient, assumed to be oxygen. In the first model the oxygen concentration u(x, t) is computed based on the geometry of the growing blob, while in the second one u(x, t) satisfies a reaction-diffusion equation. A threshold θ value exists such that cells give birth at rate β(u(x, t) - θ)⁺ and die at rate δ(θ - u(x, t)⁺. In the first model, a phase transition was found between growth as a solid blob and "fingering" at a threshold θ_c = 0.5, while in the second case fingering always occurs, i.e., θ_c = 0.

The formation of tight tumor clusters affects the efficacy of cell cycle inhibitors: a hybrid model study

Cyclin-dependent kinases (CDKs) are vital in regulating cell cycle progression, and, thus, in highly proliferating tumor cells CDK inhibitors are gaining interest as potential anticancer agents. Clonogenic assay experiments are frequently used to determine drug efficacy against the survival and proliferation of cancer cells. While the anticancer mechanisms of drugs are usually described at the intracellular single-cell level, the experimental measurements are sampled from the entire cancer cell population. This approach may lead to discrepancies between the experimental observations and theoretical explanations of anticipated drug mechanisms. To determine how individual cell responses to drugs that inhibit CDKs affect the growth of cancer cell populations, we developed a spatially explicit hybrid agent-based model. In this model, each cell is equipped with internal cell cycle regulation mechanisms, but it is also able to interact physically with its neighbors. We model cell cycle progression, focusing on the G1 and G2/M cell cycle checkpoints, as well as on related essential components, such as CDK1, CDK2, cell size, and DNA damage. We present detailed studies of how the emergent properties (e.g., cluster formation) of an entire cell population depend on altered physical and physiological parameters. We analyze the effects of CDK1 and CKD2 inhibitors on population growth, time-dependent changes in cell cycle distributions, and the dynamic evolution of spatial cell patterns. We show that cell cycle inhibitors that cause cell arrest at different cell cycle phases are not necessarily synergistically super-additive. Finally, we demonstrate that the physical aspects of cell population growth, such as the formation of tight cell clusters versus dispersed colonies, alter the efficacy of cell cycle inhibitors, both in 2D and 3D simulations. This finding may have implications for interpreting the treatment efficacy results of in vitro experiments, in which treatment is applied before the cells can grow to produce clusters, especially because in vivo tumors, in contrast, form large masses before they are detected and treated.

Cell-based multi-parametric model of cleft progression during submandibular salivary gland branching morphogenesis

Cleft formation during submandibular salivary gland branching morphogenesis is the critical step initiating the growth and development of the complex adult organ. Previous experimental studies indicated requirements for several epithelial cellular processes, such as proliferation, migration, cell-cell adhesion, cell-extracellular matrix (matrix) adhesion, and cellular contraction in cleft formation; however, the relative contribution of each of these processes is not fully understood since it is not possible to experimentally manipulate each factor independently. We present here a comprehensive analysis of several cellular parameters regulating cleft progression during branching morphogenesis in the epithelial tissue of an early embryonic salivary gland at a local scale using an on lattice Monte-Carlo simulation model, the Glazier-Graner-Hogeweg model. We utilized measurements from time-lapse images of mouse submandibular gland organ explants to construct a temporally and spatially relevant cell-based 2D model. Our model simulates the effect of cellular proliferation, actomyosin contractility, cell-cell and cell-matrix adhesions on cleft progression, and it was used to test specific hypotheses regarding the function of these parameters in branching morphogenesis. We use innovative features capturing several aspects of cleft morphology and quantitatively analyze clefts formed during functional modification of the cellular parameters. Our simulations predict that a low epithelial mitosis rate and moderate level of actomyosin contractility in the cleft cells promote cleft progression. Raising or lowering levels of contractility and mitosis rate resulted in non-progressive clefts. We also show that lowered cell-cell adhesion in the cleft region and increased cleft cell-matrix adhesions are required for cleft progression. Using a classifier-based analysis, the relative importance of these four contributing cellular factors for effective cleft progression was determined as follows: cleft cell contractility, cleft region cell-cell adhesion strength, epithelial cell mitosis rate, and cell-matrix adhesion strength.

Branching morphogenesis is a complex and dynamic embryonic process that creates the structure of many adult organs, including the salivary gland. During this process, many cellular changes occur in the epithelial cells, including changes in cell-cell adhesions, cell-extracellular matrix (matrix) adhesions, cell proliferation, and cellular contraction, resulting in formation of clefts in the epithelial cells of the organ. A comprehensive understanding of the relative contributions of these cellular processes has crucial therapeutic implications for organ regeneration and functional restoration of organ structure in diseased salivary glands. Here, we have developed a cell-based model of cleft progression and simulated cleft progression under conditions of altered cell-cell adhesions, cellular contractility, cell-matrix adhesion and cell proliferation to identify the optimum cellular conditions that cause clefts to progress. The model predicts that cleft progression requires a moderate level of cleft cell contractility, a low epithelial proliferation rate, reduced cell-cell adhesion strength in the cleft and high cell-matrix adhesion strength also in the cleft region. The results of our classification analysis demonstrate that cellular contractility in the cleft cells has a significant effect on cleft progression, followed by cell-cell adhesion strength, rate of cell proliferation, and strength of cell-matrix adhesion energies.

Cellular potts modeling of tumor growth, tumor invasion, and tumor evolution

Despite a growing wealth of available molecular data, the growth of tumors, invasion of tumors into healthy tissue, and response of tumors to therapies are still poorly understood. Although genetic mutations are in general the first step in the development of a cancer, for the mutated cell to persist in a tissue, it must compete against the other, healthy or diseased cells, for example by becoming more motile, adhesive, or multiplying faster. Thus, the cellular phenotype determines the success of a cancer cell in competition with its neighbors, irrespective of the genetic mutations or physiological alterations that gave rise to the altered phenotype. What phenotypes can make a cell "successful" in an environment of healthy and cancerous cells, and how? A widely used tool for getting more insight into that question is cell-based modeling. Cell-based models constitute a class of computational, agent-based models that mimic biophysical and molecular interactions between cells. One of the most widely used cell-based modeling formalisms is the cellular Potts model (CPM), a lattice-based, multi particle cell-based modeling approach. The CPM has become a popular and accessible method for modeling mechanisms of multicellular processes including cell sorting, gastrulation, or angiogenesis. The CPM accounts for biophysical cellular properties, including cell proliferation, cell motility, and cell adhesion, which play a key role in cancer. Multiscale models are constructed by extending the agents with intracellular processes including metabolism, growth, and signaling. Here we review the use of the CPM for modeling tumor growth, tumor invasion, and tumor progression. We argue that the accessibility and flexibility of the CPM, and its accurate, yet coarse-grained and computationally efficient representation of cell and tissue biophysics, make the CPM the method of choice for modeling cellular processes in tumor development.

Acid-mediated tumour cell invasion: a discrete modelling approach using the extended Potts model

Acidic extracellular pH has been shown to play a crucial part in the invasive and metastatic cascade of some tumours. In this study, we examine the effect of extracellular acidity on tumour invasion focusing, in particular, on cellular adhesion, proteolytic enzyme activity and cellular proliferation. Our numerical simulations using a cellular Potts model show that, under acidic extracellular pH, changes in cell-matrix adhesion strength has a comparable effect on tumour invasiveness as the increase in proteolytic enzyme activity. We also show that tumour cells cultured under physiological pH tend to be large and the tumours develop a "diffuse" morphology compared to those cultured at acidic pH, which display protruding "fingers" at the advancing front. A key model prediction is the observation that the main effect on invasion from culturing cells at low extracellular pH stems from changes in the intercellular and cell-matrix adhesion strengths and proteolytic enzyme secretion rate. However, we show that the effects of proteolysis needs to be significant as low to moderate changes only has nominal effects on cell invasiveness. We find that the low pH e effects on cell size and proliferation rate have much lower influence on cell invasiveness.

An agent-based model for elasto-plastic mechanical interactions between cells, basement membrane and extracellular matrix

The basement membrane (BM) and extracellular matrix (ECM) play critical roles in developmental and cancer biology, and are of great interest in biomathematics. We introduce a model of mechanical cell-BM-ECM interactions that extends current (visco)elastic models (e.g. [8,16]), and connects to recent agent-based cell models (e.g. [2,3,20,26]). We model the BM as a linked series of Hookean springs, each with time-varying length, thickness, and spring constant. Each BM spring node exchanges adhesive and repulsive forces with the cell agents using potential functions. We model elastic BM-ECM interactions with analogous ECM springs. We introduce a new model of plastic BM and ECM reorganization in response to prolonged strains, and new constitutive relations that incorporate molecular-scale effects of plasticity into the spring constants. We find that varying the balance of BM and ECM elasticity alters the node spacing along cell boundaries, yielding a nonuniform BM thickness. Uneven node spacing generates stresses that are relieved by plasticity over long times. We find that elasto-viscoplastic cell shape response is critical to relieving uneven stresses in the BM. Our modeling advances and results highlight the importance of rigorously modeling of cell-BM-ECM interactions in clinically important conditions with significant membrane deformations and time-varying membrane properties, such as aneurysms and progression from in situ to invasive carcinoma.

Migration of adhesive glioma cells: front propagation and fingering

We investigate the migration of glioma cells as a front propagation phenomenon both theoretically (by using both discrete lattice modeling and a continuum approach) and experimentally. For small effective strength of cell-cell adhesion q, the front velocity does not depend on q. When q exceeds a critical threshold, a fingeringlike front propagation is observed due to cluster formation in the invasive zone. We show that the experiments correspond to the transient regime, before the regime of front propagation is established. We performed an additional experiment on cell migration. A detailed comparison with experimental observations showed that the theory correctly predicts the maximal migration distance but underestimates the migration of the main mass of cells.

Adhesion failures determine the pattern of choroidal neovascularization in the eye: a computer simulation study

Choroidal neovascularization (CNV) of the macular area of the retina is the major cause of severe vision loss in adults. In CNV, after choriocapillaries initially penetrate Bruch's membrane (BrM), invading vessels may regress or expand (CNV initiation). Next, during Early and Late CNV, the expanding vasculature usually spreads in one of three distinct patterns: in a layer between BrM and the retinal pigment epithelium (sub-RPE or Type 1 CNV), in a layer between the RPE and the photoreceptors (sub-retinal or Type 2 CNV) or in both loci simultaneously (combined pattern or Type 3 CNV). While most studies hypothesize that CNV primarily results from growth-factor effects or holes in BrM, our three-dimensional simulations of multi-cell model of the normal and pathological maculae recapitulate the three growth patterns, under the hypothesis that CNV results from combinations of impairment of: 1) RPE-RPE epithelial junctional adhesion, 2) Adhesion of the RPE basement membrane complex to BrM (RPE-BrM adhesion), and 3) Adhesion of the RPE to the photoreceptor outer segments (RPE-POS adhesion). Our key findings are that when an endothelial tip cell penetrates BrM: 1) RPE with normal epithelial junctions, basal attachment to BrM and apical attachment to POS resists CNV. 2) Small holes in BrM do not, by themselves, initiate CNV. 3) RPE with normal epithelial junctions and normal apical RPE-POS adhesion, but weak adhesion to BrM (e.g. due to lipid accumulation in BrM) results in Early sub-RPE CNV. 4) Normal adhesion of RBaM to BrM, but reduced apical RPE-POS or epithelial RPE-RPE adhesion (e.g. due to inflammation) results in Early sub-retinal CNV. 5) Simultaneous reduction in RPE-RPE epithelial binding and RPE-BrM adhesion results in either sub-RPE or sub-retinal CNV which often progresses to combined pattern CNV. These findings suggest that defects in adhesion dominate CNV initiation and progression.

Integrating intracellular dynamics using CompuCell3D and Bionetsolver: applications to multiscale modelling of cancer cell growth and invasion

In this paper we present a multiscale, individual-based simulation environment that integrates CompuCell3D for lattice-based modelling on the cellular level and Bionetsolver for intracellular modelling. CompuCell3D or CC3D provides an implementation of the lattice-based Cellular Potts Model or CPM (also known as the Glazier-Graner-Hogeweg or GGH model) and a Monte Carlo method based on the metropolis algorithm for system evolution. The integration of CC3D for cellular systems with Bionetsolver for subcellular systems enables us to develop a multiscale mathematical model and to study the evolution of cell behaviour due to the dynamics inside of the cells, capturing aspects of cell behaviour and interaction that is not possible using continuum approaches. We then apply this multiscale modelling technique to a model of cancer growth and invasion, based on a previously published model of Ramis-Conde et al. (2008) where individual cell behaviour is driven by a molecular network describing the dynamics of E-cadherin and β-catenin. In this model, which we refer to as the centre-based model, an alternative individual-based modelling technique was used, namely, a lattice-free approach. In many respects, the GGH or CPM methodology and the approach of the centre-based model have the same overall goal, that is to mimic behaviours and interactions of biological cells. Although the mathematical foundations and computational implementations of the two approaches are very different, the results of the presented simulations are compatible with each other, suggesting that by using individual-based approaches we can formulate a natural way of describing complex multi-cell, multiscale models. The ability to easily reproduce results of one modelling approach using an alternative approach is also essential from a model cross-validation standpoint and also helps to identify any modelling artefacts specific to a given computational approach.

Contour instabilities in early tumor growth models

Recent tumor growth models are often based on the multiphase mixture framework. Using bifurcation theory techniques, we show that such models can give contour instabilities. Restricting to a simplified but realistic version of such models, with an elastic cell-to-cell interaction and a growth rate dependent on diffusing nutrients, we prove that the tumor cell concentration at the border acts as a control parameter inducing a bifurcation with loss of the circular symmetry. We show that the finite wavelength at threshold has the size of the proliferating peritumoral zone. We apply our predictions to melanoma growth since contour instabilities are crucial for early diagnosis. Given the generality of the equations, other relevant applications can be envisaged for solving problems of tissue growth and remodeling.

Hybrid models of tumor growth

Cancer is a complex, multiscale process in which genetic mutations occurring at a subcellular level manifest themselves as functional changes at the cellular and tissue scale. The multiscale nature of cancer requires mathematical modeling approaches that can handle multiple intracellular and extracellular factors acting on different time and space scales. Hybrid models provide a way to integrate both discrete and continuous variables that are used to represent individual cells and concentration or density fields, respectively. Each discrete cell can also be equipped with submodels that drive cell behavior in response to microenvironmental cues. Moreover, the individual cells can interact with one another to form and act as an integrated tissue. Hybrid models form part of a larger class of individual-based models that can naturally connect with tumor cell biology and allow for the integration of multiple interacting variables both intrinsically and extrinsically and are therefore perfectly suited to a systems biology approach to tumor growth.

Collective behavior of brain tumor cells: the role of hypoxia

We consider emergent collective behavior of a multicellular biological system. Specifically, we investigate the role of hypoxia (lack of oxygen) in migration of brain tumor cells. We performed two series of cell migration experiments. In the first set of experiments, cell migration away from a tumor spheroid was investigated. The second set of experiments was performed in a typical wound-healing geometry: Cells were placed on a substrate, a scratch was made, and cell migration into the gap was investigated. Experiments show a surprising result: Cells under normal and hypoxic conditions have migrated the same distance in the "spheroid" experiment, while in the "scratch" experiment cells under normal conditions migrated much faster than under hypoxic conditions. To explain this paradox, we formulate a discrete stochastic model for cell dynamics. The theoretical model explains our experimental observations and suggests that hypoxia decreases both the motility of cells and the strength of cell-cell adhesion. The theoretical predictions were further verified in independent experiments.

Mathematical modeling of cancer cell invasion of tissue: biological insight from mathematical analysis and computational simulation

The ability of cancer cells to break out of tissue compartments and invade locally gives solid tumours a defining deadly characteristic. One of the first steps of invasion is the remodelling of the surrounding tissue or extracellular matrix (ECM) and a major part of this process is the over-expression of proteolytic enzymes, such as the urokinase-type plasminogen activator (uPA) and matrix metalloproteinases (MMPs), by the cancer cells to break down ECM proteins. Degradation of the matrix enables the cancer cells to migrate through the tissue and subsequently to spread to secondary sites in the body, a process known as metastasis. In this paper we undertake an analysis of a mathematical model of cancer cell invasion of tissue, or ECM, which focuses on the role of the urokinase plasminogen activation system. The model consists of a system of five reaction-diffusion-taxis partial differential equations describing the interactions between cancer cells, uPA, uPA inhibitors, plasmin and the host tissue. Cancer cells react chemotactically and haptotactically to the spatio-temporal effects of the uPA system. The results obtained from computational simulations carried out on the model equations produce dynamic heterogeneous spatio-temporal solutions and using linear stability analysis we show that this is caused by a taxis-driven instability of a spatially homogeneous steady-state. Finally we consider the biological implications of the model results, draw parallels with clinical samples and laboratory based models of cancer cell invasion using three-dimensional invasion assay, and go on to discuss future development of the model.

Front instabilities and invasiveness of simulated 3D avascular tumors

We use the Glazier-Graner-Hogeweg model to simulate three-dimensional (3D), single-phenotype, avascular tumors growing in an homogeneous tissue matrix (TM) supplying a single limiting nutrient. We study the effects of two parameters on tumor morphology: a diffusion-limitation parameter defined as the ratio of the tumor-substrate consumption rate to the substrate-transport rate, and the tumor-TM surface tension. This initial model omits necrosis and oxidative/hypoxic metabolism effects, which can further influence tumor morphology, but our simplified model still shows significant parameter dependencies. The diffusion-limitation parameter determines whether the growing solid tumor develops a smooth (noninvasive) or fingered (invasive) interface, as in our earlier two-dimensional (2D) simulations. The sensitivity of 3D tumor morphology to tumor-TM surface tension increases with the size of the diffusion-limitation parameter, as in 2D. The 3D results are unexpectedly close to those in 2D. Our results therefore may justify using simpler 2D simulations of tumor growth, instead of more realistic but more computationally expensive 3D simulations. While geometrical artifacts mean that 2D sections of connected 3D tumors may be disconnected, the morphologies of 3D simulated tumors nevertheless correlate with the morphologies of their 2D sections, especially for low-surface-tension tumors, allowing the use of 2D sections to partially reconstruct medically-important 3D-tumor structures.

An off-lattice hybrid discrete-continuum model of tumor growth and invasion

We have developed an off-lattice hybrid discrete-continuum (OLHDC) model of tumor growth and invasion. The continuum part of the OLHDC model describes microenvironmental components such as matrix-degrading enzymes, nutrients or oxygen, and extracellular matrix (ECM) concentrations, whereas the discrete portion represents individual cell behavior such as cell cycle, cell-cell, and cell-ECM interactions and cell motility by the often-used persistent random walk, which can be depicted by the Langevin equation. Using this framework of the OLHDC model, we develop a phenomenologically realistic and bio/physically relevant model that encompasses the experimentally observed superdiffusive behavior (at short times) of mammalian cells. When systemic simulations based on the OLHDC model are performed, tumor growth and its morphology are found to be strongly affected by cell-cell adhesion and haptotaxis. There is a combination of the strength of cell-cell adhesion and haptotaxis in which fingerlike shapes, characteristic of invasive tumor, are observed.

Nonlinear modelling of cancer: bridging the gap between cells and tumours

Despite major scientific, medical and technological advances over the last few decades, a cure for cancer remains elusive. The disease initiation is complex, and including initiation and avascular growth, onset of hypoxia and acidosis due to accumulation of cells beyond normal physiological conditions, inducement of angiogenesis from the surrounding vasculature, tumour vascularization and further growth, and invasion of surrounding tissue and metastasis. Although the focus historically has been to study these events through experimental and clinical observations, mathematical modelling and simulation that enable analysis at multiple time and spatial scales have also complemented these efforts. Here, we provide an overview of this multiscale modelling focusing on the growth phase of tumours and bypassing the initial stage of tumourigenesis. While we briefly review discrete modelling, our focus is on the continuum approach. We limit the scope further by considering models of tumour progression that do not distinguish tumour cells by their age. We also do not consider immune system interactions nor do we describe models of therapy. We do discuss hybrid-modelling frameworks, where the tumour tissue is modelled using both discrete (cell-scale) and continuum (tumour-scale) elements, thus connecting the micrometre to the centimetre tumour scale. We review recent examples that incorporate experimental data into model parameters. We show that recent mathematical modelling predicts that transport limitations of cell nutrients, oxygen and growth factors may result in cell death that leads to morphological instability, providing a mechanism for invasion via tumour fingering and fragmentation. These conditions induce selection pressure for cell survivability, and may lead to additional genetic mutations. Mathematical modelling further shows that parameters that control the tumour mass shape also control its ability to invade. Thus, tumour morphology may serve as a predictor of invasiveness and treatment prognosis.

Microenvironmental independence associated with tumor progression

Tumor-microenvironment interactions are increasingly recognized to influence tumor progression. To understand the competitive dynamics of tumor cells in diverse microenvironments, we experimentally parameterized a hybrid discrete-continuum mathematical model with phenotypic trait data from a set of related mammary cell lines with normal, transformed, or tumorigenic properties. Surprisingly, in a resource-rich microenvironment, with few limitations on proliferation or migration, transformed (but not tumorigenic) cells were most successful and outcompeted other cell types in heterogeneous tumor simulations. Conversely, constrained microenvironments with limitations on space and/or growth factors gave a selective advantage to phenotypes derived from tumorigenic cell lines. Analysis of the relative performance of each phenotype in constrained versus unconstrained microenvironments revealed that, although all cell types grew more slowly in resource-constrained microenvironments, the most aggressive cells were least affected by microenvironmental constraints. A game theory model testing the relationship between microenvironment resource availability and competitive cellular dynamics supports the concept that microenvironmental independence is an advantageous cellular trait in resource-limited microenvironments.

3D multi-cell simulation of tumor growth and angiogenesis

We present a 3D multi-cell simulation of a generic simplification of vascular tumor growth which can be easily extended and adapted to describe more specific vascular tumor types and host tissues. Initially, tumor cells proliferate as they take up the oxygen which the pre-existing vasculature supplies. The tumor grows exponentially. When the oxygen level drops below a threshold, the tumor cells become hypoxic and start secreting pro-angiogenic factors. At this stage, the tumor reaches a maximum diameter characteristic of an avascular tumor spheroid. The endothelial cells in the pre-existing vasculature respond to the pro-angiogenic factors both by chemotaxing towards higher concentrations of pro-angiogenic factors and by forming new blood vessels via angiogenesis. The tumor-induced vasculature increases the growth rate of the resulting vascularized solid tumor compared to an avascular tumor, allowing the tumor to grow beyond the spheroid in these linear-growth phases. First, in the linear-spherical phase of growth, the tumor remains spherical while its volume increases. Second, in the linear-cylindrical phase of growth the tumor elongates into a cylinder. Finally, in the linear-sheet phase of growth, tumor growth accelerates as the tumor changes from cylindrical to paddle-shaped. Substantial periods during which the tumor grows slowly or not at all separate the exponential from the linear-spherical and the linear-spherical from the linear-cylindrical growth phases. In contrast to other simulations in which avascular tumors remain spherical, our simulated avascular tumors form cylinders following the blood vessels, leading to a different distribution of hypoxic cells within the tumor. Our simulations cover time periods which are long enough to produce a range of biologically reasonable complex morphologies, allowing us to study how tumor-induced angiogenesis affects the growth rate, size and morphology of simulated tumors.