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

Semi-autonomous wound invasion via matrix-deposited, haptotactic cues

Proper wound healing relies on invasion of fibroblasts via directed migration. While the related experimental and mathematical modeling literature has mainly focused on cell migration directed by soluble cues (chemotaxis), there is ample evidence that fibroblast migration is also directed by insoluble, matrix-bound cues (haptotaxis). Furthermore, numerous studies indicate that fibronectin (FN), a haptotactic ligand for fibroblasts, is present and dynamic in the provisional matrix throughout the proliferative phase of wound healing. In the present work, we show the plausibility of a hypothesis that fibroblasts themselves form and maintain haptotactic gradients in a semi-autonomous fashion. As a precursor to this, we examine the positive control scenario where FN is pre-deposited in the wound matrix, and fibroblasts maintain haptotaxis by removing FN at an appropriate rate. After developing conceptual and quantitative understanding of this scenario, we consider two cases in which fibroblasts activate the latent form of a matrix-loaded cytokine, TGFβ, which upregulates the fibroblasts' own secretion of FN. In the first of these, the latent cytokine is pre-patterned and released by the fibroblasts. In the second, fibroblasts in the wound produce the latent TGFβ, with the presence of the wound providing the only instruction. In all cases, wound invasion is more effective than a negative control model with haptotaxis disabled; however, there is a trade-off between the degree of fibroblast autonomy and the rate of invasion.

State-dependent evolutionary models reveal modes of solid tumour growth

Spatial properties of tumour growth have profound implications for cancer progression, therapeutic resistance and metastasis. Yet, how spatial position governs tumour cell division remains difficult to evaluate in clinical tumours. Here, we demonstrate that faster division on the tumour periphery leaves characteristic genetic patterns, which become evident when a phylogenetic tree is reconstructed from spatially sampled cells. Namely, rapidly dividing peripheral lineages branch more extensively and acquire more mutations than slower-dividing centre lineages. We develop a Bayesian state-dependent evolutionary phylodynamic model (SDevo) that quantifies these patterns to infer the differential division rates between peripheral and central cells. We demonstrate that this approach accurately infers spatially varying birth rates of simulated tumours across a range of growth conditions and sampling strategies. We then show that SDevo outperforms state-of-the-art, non-cancer multi-state phylodynamic methods that ignore differential sequence evolution. Finally, we apply SDevo to single-time-point, multi-region sequencing data from clinical hepatocellular carcinomas and find evidence of a three- to six-times-higher division rate on the tumour edge. With the increasing availability of high-resolution, multi-region sequencing, we anticipate that SDevo will be useful in interrogating spatial growth restrictions and could be extended to model non-spatial factors that influence tumour progression.

The role of acoustofluidics and microbubble dynamics for therapeutic applications and drug delivery

Targeted drug delivery is proposed to reduce the toxic effects of conventional therapeutic methods. For that purpose, nanoparticles are loaded with drugs called nanocarriers and directed toward a specific site. However, biological barriers challenge the nanocarriers to convey the drug to the target site effectively. Different targeting strategies and nanoparticle designs are used to overcome these barriers. Ultrasound is a new, safe, and non-invasive drug targeting method, especially when combined with microbubbles. Microbubbles oscillate under the effect of the ultrasound, which increases the permeability of endothelium, hence, the drug uptake to the target site. Consequently, this new technique reduces the dose of the drug and avoids its side effects. This review aims to describe the biological barriers and the targeting types with the critical features of acoustically driven microbubbles focusing on biomedical applications. The theoretical part covers the historical developments in microbubble models for different conditions: microbubbles in an incompressible and compressible medium and bubbles encapsulated by a shell. The current state and the possible future directions are discussed.

Extracellular matrix density regulates the formation of tumour spheroids through cell migration

In this work, we show how the mechanical properties of the cellular microenvironment modulate the growth of tumour spheroids. Based on the composition of the extracellular matrix, its stiffness and architecture can significantly vary, subsequently influencing cell movement and tumour growth. However, it is still unclear exactly how both of these processes are regulated by the matrix composition. Here, we present a centre-based computational model that describes how collagen density, which modulates the steric hindrance properties of the matrix, governs individual cell migration and, consequently, leads to the formation of multicellular clusters of varying size. The model was calibrated using previously published experimental data, replicating a set of experiments in which cells were seeded in collagen matrices of different collagen densities, hence producing distinct mechanical properties. At an initial stage, we tracked individual cell trajectories and speeds. Subsequently, the formation of multicellular clusters was also analysed by quantifying their size. Overall, the results showed that our model could accurately replicate what was previously seen experimentally. Specifically, we showed that cells seeded in matrices with low collagen density tended to migrate more. Accordingly, cells strayed away from their original cluster and thus promoted the formation of small structures. In contrast, we also showed that high collagen densities hindered cell migration and produced multicellular clusters with increased volume. In conclusion, this model not only establishes a relation between matrix density and individual cell migration but also showcases how migration, or its inhibition, modulates tumour growth.

An in-silico study of cancer cell survival and spatial distribution within a 3D microenvironment

3D cell cultures are in-vitro models representing a significant improvement with respect to traditional monolayers. Their diffusion and applicability, however, are hampered by the complexity of 3D systems, that add new physical variables for experimental analyses. In order to account for these additional features and improve the study of 3D cultures, we here present SALSA (ScAffoLd SimulAtor), a general purpose computational tool that can simulate the behavior of a population of cells cultured in a 3D scaffold. This software allows for the complete customization of both the polymeric template structure and the cell population behavior and characteristics. In the following the technical description of SALSA will be presented, together with its validation and an example of how it could be used to optimize the experimental analysis of two breast cancer cell lines cultured in collagen scaffolds. This work contributes to the growing field of integrated in-silico/in-vitro analysis of biological systems, which have great potential for the study of complex cell population behaviours and could lead to improve and facilitate the effectiveness and diffusion of 3D cell culture models.

Cancer systems immunology

Tumor immunology is undergoing a renaissance due to the recent profound clinical successes of tumor immunotherapy. These advances have coincided with an exponential growth in the development of -omics technologies. Armed with these technologies and their associated computational and modeling toolsets, systems biologists have turned their attention to tumor immunology in an effort to understand the precise nature and consequences of interactions between tumors and the immune system. Such interactions are inherently multivariate, spanning multiple time and size scales, cell types, and organ systems, rendering systems biology approaches particularly amenable to their interrogation. While in its infancy, the field of 'Cancer Systems Immunology' has already influenced our understanding of tumor immunology and immunotherapy. As the field matures, studies will move beyond descriptive characterizations toward functional investigations of the emergent behavior that govern tumor-immune responses. Thus, Cancer Systems Immunology holds incredible promise to advance our ability to fight this disease.

PhysiBoSS: a multi-scale agent-based modelling framework integrating physical dimension and cell signalling

MOTIVATION

Due to the complexity and heterogeneity of multicellular biological systems, mathematical models that take into account cell signalling, cell population behaviour and the extracellular environment are particularly helpful. We present PhysiBoSS, an open source software which combines intracellular signalling using Boolean modelling (MaBoSS) and multicellular behaviour using agent-based modelling (PhysiCell).

RESULTS

PhysiBoSS provides a flexible and computationally efficient framework to explore the effect of environmental and genetic alterations of individual cells at the population level, bridging the critical gap from single-cell genotype to single-cell phenotype and emergent multicellular behaviour. PhysiBoSS thus becomes very useful when studying heterogeneous population response to treatment, mutation effects, different modes of invasion or isomorphic morphogenesis events. To concretely illustrate a potential use of PhysiBoSS, we studied heterogeneous cell fate decisions in response to TNF treatment. We explored the effect of different treatments and the behaviour of several resistant mutants. We highlighted the importance of spatial information on the population dynamics by considering the effect of competition for resources like oxygen.

AVAILABILITY AND IMPLEMENTATION

PhysiBoSS is freely available on GitHub (https://github.com/sysbio-curie/PhysiBoSS), with a Docker image (https://hub.docker.com/r/gletort/physiboss/). It is distributed as open source under the BSD 3-clause license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Hybrid Model of Cartilage Regeneration Capturing the Interactions Between Cellular Dynamics and Porosity

To accelerate the development of strategies for cartilage tissue engineering, models are necessary to investigate the interactions between cellular dynamics and the local microenvironment. We use a discrete framework to capture the individual behavior of cells, modeling experiments where cells are seeded in a porous scaffold or hydrogel and over the time course of a month, the scaffold slowly degrades while cells divide and synthesize extracellular matrix constituents. The movement of cells and the ability to proliferate is a function of the local porosity, defined as the volume fraction of fluid in the surrounding region. A phenomenological approach is used to capture a continuous profile for the degrading scaffold and accumulating matrix, which will then change the local porosity throughout the construct. We parameterize the model by first matching total cell counts in the construct to chondrocytes seeded in a polyglycolic acid scaffold (Freed et al. in Biotechnol Bioeng 43:597-604, 1994). We investigate the influence of initial scaffold porosity on the total cell count and spatial profiles of cell and ECM in the construct. Cell counts were higher at day 30 in scaffolds of lower initial porosity, and similar cell counts were obtained using different models of scaffold degradation and matrix accumulation (either uniform or cell-specific). Using this modeling framework, we study the interplay between a phenomenological representation of scaffold architecture and porosity as well as the potential continuous application of growth factors. We determine parameter regimes where large cellular aggregates occur, which can hinder matrix accumulation and cellular proliferation. The developed modeling framework can easily be extended and can be used to identify optimal scaffolds and culture conditions that lead to a desired distribution of extracellular matrix and cell counts throughout the construct.

Mathematical Models of Cancer Cell Plasticity

Quantitative modelling is increasingly important in cancer research, helping to integrate myriad diverse experimental data into coherent pictures of the disease and able to discriminate between competing hypotheses or suggest specific experimental lines of enquiry and new approaches to therapy. Here, we review a diverse set of mathematical models of cancer cell plasticity (a process in which, through genetic and epigenetic changes, cancer cells survive in hostile environments and migrate to more favourable environments, respectively), tumour growth, and invasion. Quantitative models can help to elucidate the complex biological mechanisms of cancer cell plasticity. In this review, we discuss models of plasticity, tumour progression, and metastasis under three broadly conceived mathematical modelling techniques: discrete, continuum, and hybrid, each with advantages and disadvantages. An emerging theme from the predictions of many of these models is that cell escape from the tumour microenvironment (TME) is encouraged by a combination of physiological stress locally (e.g., hypoxia), external stresses (e.g., the presence of immune cells), and interactions with the extracellular matrix. We also discuss the value of mathematical modelling for understanding cancer more generally.

Multiscale modeling of solid stress and tumor cell invasion in response to dynamic mechanical microenvironment

Mathematical models can provide a quantitatively sophisticated description of tumor cell (TC) behaviors under mechanical microenvironment and help us better understand the role of specific biophysical factors based on their influences on the TC behaviors. To this end, we propose an off-lattice cell-based multiscale mathematical model to describe the dynamic growth-induced solid stress during tumor progression and investigate the influence of the mechanical microenvironment on TC invasion. At the cellular level, cell-cell and cell-matrix interactive forces depend on the mechanical properties of the cells and the cancer-associated fibroblasts in the stroma, respectively. The constitutive relationship between the interactive forces and cell migrations obeys the Hooke's law and damping effects. At the tissue level, the integrated growth-induced forces caused by proliferating cells within the simulation region are balanced by the external forces applied by the surrounding host tissues. Then, the cell movements are calculated according to the Newton's second law of motion, and the morphology of TC invasion is updated. The simulation results reveal the continuous changes of the macroscopic mechanical forces due to the interactions among the structural components and the microscopic environmental factors. Moreover, the simulation results demonstrate the adverse effect of the stiffness of tumor tissue on tumor growth and invasion. A decrease in the stiffness of tumor and matrix can promote TCs to proliferate at a much faster rate and invade into the surrounding healthy tissue more easily, whereas an increase in the stiffness can lead to an aggressive morphology of tumor invasion. We envision that the proposed model can be served as a quantitative theoretical platform to study the underlying biophysical role of the mechanical microenvironmental factors during tumor invasion and metastasis.

Continuum descriptions of spatial spreading for heterogeneous cell populations: Theory and experiment

Variability in cell populations is frequently observed in both in vitro and in vivo settings. Intrinsic differences within populations of cells, such as differences in cell sizes or differences in rates of cell motility, can be present even within a population of cells from the same cell line. We refer to this variability as cell heterogeneity. Mathematical models of cell migration, for example, in the context of tumour growth and metastatic invasion, often account for both undirected (random) migration and directed migration that is mediated by cell-to-cell contacts and cell-to-cell adhesion. A key feature of standard models is that they often assume that the population is composed of identical cells with constant properties. This leads to relatively simple single-species homogeneous models that neglect the role of heterogeneity. In this work, we use a continuum modelling approach to explore the role of heterogeneity in spatial spreading of cell populations. We employ a three-species heterogeneous model of cell motility that explicitly incorporates different types of experimentally-motivated heterogeneity in cell sizes: (i) monotonically decreasing; (ii) uniform; (iii) non-monotonic; and (iv) monotonically increasing distributions of cell size. Comparing the density profiles generated by the three-species heterogeneous model with density profiles predicted by a more standard single-species homogeneous model reveals that when we are dealing with monotonically decreasing and uniform distributions a simple and computationally efficient single-species homogeneous model can be remarkably accurate in describing the evolution of a heterogeneous cell population. In contrast, we find that the simpler single-species homogeneous model performs relatively poorly when applied to non-monotonic and monotonically increasing distributions of cell sizes. Additional results for heterogeneity in parameters describing both undirected and directed cell migration are also considered, and we find that similar results apply.

Effects of mutations and immunogenicity on outcomes of anti-cancer therapies for secondary lesions

Cancer development is driven by mutations and selective forces, including the action of the immune system and interspecific competition. When administered to patients, anti-cancer therapies affect the development and dynamics of tumours, possibly with various degrees of resistance due to immunoediting and microenvironment. Tumours are able to express a variety of competing phenotypes with different attributes and thus respond differently to various anti-cancer therapies. In this paper, a mathematical framework incorporating a system of delay differential equations for the immune system activation cycle and an agent-based approach for tumour-immune interaction is presented. The focus is on those metastatic, secondary solid lesions that are still undetected and non-vascularised. By using available experimental data, we analyse the effects of combination therapies on these lesions and investigate the role of mutations on the rates of success of common treatments. Findings show that mutations, growth properties and immunoediting influence therapies' outcomes in nonlinear and complex ways, affecting cancer lesion morphologies, phenotypical compositions and overall proliferation patterns. Cascade effects on final outcomes for secondary lesions are also investigated, showing that actions on primary lesions could sometimes result in unexpected clearances of secondary tumours. This outcome is strongly dependent on the clonal composition of the primary and secondary masses and is shown to allow, in some cases, the control of the disease for years.

A 3D tumor spheroid model for the T98G Glioblastoma cell line phenotypic characterization

Major Glioblastoma's hallmarks include proliferation, invasion and heterogeneity. Biological 3D tumor spheroid models can serve as intermediate systems between traditional 2D cell culture and complex in vivo models. Tumor spheroids have been shown to more accurately reproduce the spatial organization and microenvironmental factors of in vivo micro-tumors, such as relevant gradients of nutrients and other molecular agents, while they maintain cell-to-cell and cell-to-matrix interactions. In vitro 3D assays are useful to monitor these properties. Here, we test the suitability of the well-known T98 G Glioblastoma cell line in such a 3D assay. The doubling time and death rate parameters of T98 G are estimated, as well as their spheroidal growth-expansion curves with and without the presence of basement membrane substrate. The T98 G invasive profile is characterized by collective morphology and proliferation-associated invasion. We show that the T98 G secondary GB cell line exhibits both invasive and proliferative capabilities in 3D and thus, can serve as control cell line for the 3D in vitro study of primary GB cell cultures.

Integrating in vitro experiments with in silico approaches for Glioblastoma invasion: the role of cell-to-cell adhesion heterogeneity

Glioblastoma cells adopt migration strategies to invade into the brain parenchyma ranging from individual to collective mechanisms, whose role and dynamics are not yet fully understood. In this work, we explore Glioblastoma heterogeneity and recapitulate its invasive patterns both in vitro, by utilizing primary cells along with the U87MG cell line, and in silico, by adopting discrete, individual cell-based mathematics. Glioblastoma cells are cultured three-dimensionally in an ECM-like substrate. The primary Glioblastoma spheroids adopt a novel cohesive pattern, mimicking perivascular invasion in the brain, while the U87MG adopt a typical, starburst invasive pattern under the same experimental setup. Mathematically, we focus on the role of the intrinsic heterogeneity with respect to cell-to-cell adhesion. Our proposed mathematical approach mimics the invasive morphologies observed in vitro and predicts the dynamics of tumour expansion. The role of the proliferation and migration is also explored showing that their effect on tumour morphology is different per cell type. The proposed model suggests that allowing cell-to-cell adhesive heterogeneity within the tumour population is sufficient for variable invasive morphologies to emerge which remain originally undetectable by conventional imaging, indicating that exploration in pathological samples is needed to improve our understanding and reveal potential patient-specific therapeutic targets.

Dynamical density-functional-theory-based modeling of tissue dynamics: Application to tumor growth

We present a theoretical framework based on an extension of dynamical density-functional theory (DDFT) for describing the structure and dynamics of cells in living tissues and tumors. DDFT is a microscopic statistical mechanical theory for the time evolution of the density distribution of interacting many-particle systems. The theory accounts for cell-pair interactions, different cell types, phenotypes, and cell birth and death processes (including cell division), to provide a biophysically consistent description of processes bridging across the scales, including describing the tissue structure down to the level of the individual cells. Analysis of the model is presented for single-species and two-species cases, the latter aimed at describing competition between tumor and healthy cells. In suitable parameter regimes, model results are consistent with biological observations. Of particular note, divergent tumor growth behavior, mirroring metastatic and benign growth characteristics, are shown to be dependent on the cell-pair-interaction parameters.

Glioma growth modeling based on the effect of vital nutrients and metabolic products

Glioma brain tumors exhibit considerably aggressive behavior leading to high mortality rates. Mathematical modeling of tumor growth aims to explore the interactions between glioma cells and tissue microenvironment, which affect tumor evolution. Leveraging this concept, we present a three-dimensional model of glioma spatio-temporal evolution based on existing continuum approaches, yet incorporating novel factors of the phenomenon. The proposed model involves the interactions between different tumor cell phenotypes and their microenvironment, investigating how tumor growth is affected by complex biological exchanges. It focuses on the separate and combined effect of vital nutrients and cellular wastes on tumor expansion, leading to the formation of cell populations with different metabolic, proliferative, and diffusive profiles. Several simulations were performed on a virtual and a real glioma, using combinations of proliferation and diffusion rates for different evolution times. The model results were validated on a glioma model available in the literature and a real case of tumor progression. The experimental observations indicate that our model estimates quite satisfactorily the expansion of each region and the overall tumor growth. Based on the individual results, the proposed model may provide an important research tool for patient-specific simulation of different tumor evolution scenarios and reliable estimation of glioma evolution. Graphical Abstract Outline of the mathematical model functionality and application to glioma growth with indicative results.

Translational models of tumor angiogenesis: A nexus of in silico and in vitro models

Emerging evidence shows that endothelial cells are not only the building blocks of vascular networks that enable oxygen and nutrient delivery throughout a tissue but also serve as a rich resource of angiocrine factors. Endothelial cells play key roles in determining cancer progression and response to anti-cancer drugs. Furthermore, the endothelium-specific deposition of extracellular matrix is a key modulator of the availability of angiocrine factors to both stromal and cancer cells. Considering tumor vascular network as a decisive factor in cancer pathogenesis and treatment response, these networks need to be an inseparable component of cancer models. Both computational and in vitro experimental models have been extensively developed to model tumor-endothelium interactions. While informative, they have been developed in different communities and do not yet represent a comprehensive platform. In this review, we overview the necessity of incorporating vascular networks for both in vitro and in silico cancer models and discuss recent progresses and challenges of in vitro experimental microfluidic cancer vasculature-on-chip systems and their in silico counterparts. We further highlight how these two approaches can merge together with the aim of presenting a predictive combinatorial platform for studying cancer pathogenesis and testing the efficacy of single or multi-drug therapeutics for cancer treatment.

A Brownian dynamics tumor progression simulator with application to glioblastoma

Tumor progression modeling offers the potential to predict tumor-spreading behavior to improve prognostic accuracy and guide therapy development. Common simulation methods include continuous reaction-diffusion (RD) approaches that capture mean spatio-temporal tumor spreading behavior and discrete agent-based (AB) approaches which capture individual cell events such as proliferation or migration. The brain cancer glioblastoma (GBM) is especially appropriate for such proliferation-migration modeling approaches because tumor cells seldom metastasize outside of the central nervous system and cells are both highly proliferative and migratory. In glioblastoma research, current RD estimates of proliferation and migration parameters are derived from computed tomography or magnetic resonance images. However, these estimates of glioblastoma cell migration rates, modeled as a diffusion coefficient, are approximately 1-2 orders of magnitude larger than single-cell measurements in animal models of this disease. To identify possible sources for this discrepancy, we evaluated the fundamental RD simulation assumptions that cells are point-like structures that can overlap. To give cells physical size (~10 μm), we used a Brownian dynamics approach that simulates individual single-cell diffusive migration, growth, and proliferation activity via a gridless, off-lattice, AB method where cells can be prohibited from overlapping each other. We found that for realistic single-cell parameter growth and migration rates, a non-overlapping model gives rise to a jammed configuration in the center of the tumor and a biased outward diffusion of cells in the tumor periphery, creating a quasi-ballistic advancing tumor front. The simulations demonstrate that a fast-progressing tumor can result from minimally diffusive cells, but at a rate that is still dependent on single-cell diffusive migration rates. Thus, modeling with the assumption of physically-grounded volume conservation can account for the apparent discrepancy between estimated and measured diffusion of GBM cells and provide a new theoretical framework that naturally links single-cell growth and migration dynamics to tumor-level progression.

Spatially extended hybrid methods: a review

Many biological and physical systems exhibit behaviour at multiple spatial, temporal or population scales. Multiscale processes provide challenges when they are to be simulated using numerical techniques. While coarser methods such as partial differential equations are typically fast to simulate, they lack the individual-level detail that may be required in regions of low concentration or small spatial scale. However, to simulate at such an individual level throughout a domain and in regions where concentrations are high can be computationally expensive. Spatially coupled hybrid methods provide a bridge, allowing for multiple representations of the same species in one spatial domain by partitioning space into distinct modelling subdomains. Over the past 20 years, such hybrid methods have risen to prominence, leading to what is now a very active research area across multiple disciplines including chemistry, physics and mathematics. There are three main motivations for undertaking this review. Firstly, we have collated a large number of spatially extended hybrid methods and presented them in a single coherent document, while comparing and contrasting them, so that anyone who requires a multiscale hybrid method will be able to find the most appropriate one for their need. Secondly, we have provided canonical examples with algorithms and accompanying code, serving to demonstrate how these types of methods work in practice. Finally, we have presented papers that employ these methods on real biological and physical problems, demonstrating their utility. We also consider some open research questions in the area of hybrid method development and the future directions for the field.

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.

A model of the effects of cancer cell motility and cellular adhesion properties on tumour-immune dynamics

We present a three-dimensional model simulating the dynamics of an anti-cancer T-cell response against a small, avascular, early-stage tumour. Interactions at the tumour site are accounted for using an agent-based model (ABM), while immune cell dynamics in the lymph node are modelled as a system of delay differential equations (DDEs). We combine these separate approaches into a two-compartment hybrid ABM-DDE system to capture the T-cell response against the tumour. In the ABM at the tumour site, movement of tumour cells is modelled using effective physical forces with a specific focus on cell-to-cell adhesion properties and varying levels of tumour cell motility, thus taking into account the ability of cancer cells to spread and form clusters. We consider the effectiveness of the immune response over a range of parameters pertaining to tumour cell motility, cell-to-cell adhesion strength and growth rate. We also investigate the dependence of outcomes on the distribution of tumour cells. Low tumour cell motility is generally a good indicator for successful tumour eradication before relapse, while high motility leads, almost invariably, to relapse and tumour escape. In general, the effect of cell-to-cell adhesion on prognosis is dependent on the level of tumour cell motility, with an often unpredictable cross influence between adhesion and motility, which can lead to counterintuitive effects. In terms of overall tumour shape and structure, the spatial distribution of cancer cells in clusters of various sizes has shown to be strongly related to the likelihood of extinction.

Continuum approximations for lattice-free multi-species models of collective cell migration

Cell migration within tissues involves the interaction of many cells from distinct subpopulations. In this work, we present a discrete model of collective cell migration where the motion of individual cells is driven by random forces, short range repulsion forces to mimic crowding, and longer range attraction forces to mimic adhesion. This discrete model can be used to simulate a population of cells that is composed of K ≥ 1 distinct subpopulations. To analyse the discrete model we formulate a hierarchy of moment equations that describe the spatial evolution of the density of agents, pairs of agents, triplets of agents, and so forth. To solve the hierarchy of moment equations we introduce two forms of closure: (i) the mean field approximation, which effectively assumes that the distributions of individual agents are independent; and (ii) a moment dynamics description that is based on the Kirkwood superposition approximation. The moment dynamics description provides an approximate way of incorporating spatial patterns, such as agent clustering, into the continuum description. Comparing the performance of the two continuum descriptions confirms that both perform well when adhesive forces are sufficiently weak. In contrast, the moment dynamics description outperforms the mean field model when adhesive forces are sufficiently large. This is a first attempt to provide an accurate continuum description of a lattice-free, multi-species model of collective cell migration.

Cell-alignment patterns in the collective migration of cells with polarized adhesion

Dictyostelium discoideum (Dd) utilizes inhomogeneities in the distribution of cell-cell adhesion molecules on cell membranes for collective cell migration. A simple example of an inhomogeneity is a front-side (leading-edge) polarization in the distribution at the early streaming stage. Experiments have shown that the polarized cell-cell adhesion induces side-by-side contact between cells [Beug et al., Nature (London) 274, 445 (1978)NATUAS0028-083610.1038/274445a0]. This result is counterintuitive, as one would expect cells to align front to front in contact with each other on the basis of front-side polarization. In this work, we theoretically examine whether front-side polarization induces side-by-side contact in collective cell migration. We construct a model for expressing cells with this polarization based on the two-dimensional cellular Potts model. By a numerical simulation with this model, we find cell-cell alignment wherein cells form lateral arrays with side-by-side contacts as observed in the experiments.

Modeling the two-way feedback between contractility and matrix realignment reveals a nonlinear mode of cancer cell invasion

Cancer cell invasion from primary tumors is mediated by a complex interplay between cellular adhesions, actomyosin-driven contractility, and the physical characteristics of the extracellular matrix (ECM). Here, we incorporate a mechanochemical free-energy-based approach to elucidate how the two-way feedback loop between cell contractility (induced by the activity of chemomechanical interactions such as Ca2+ and Rho signaling pathways) and matrix fiber realignment and strain stiffening enables the cells to polarize and develop contractile forces to break free from the tumor spheroids and invade into the ECM. Interestingly, through this computational model, we are able to identify a critical stiffness that is required by the matrix to break intercellular adhesions and initiate cell invasion. Also, by considering the kinetics of the cell movement, our model predicts a biphasic invasiveness with respect to the stiffness of the matrix. These predictions are validated by analyzing the invasion of melanoma cells in collagen matrices of varying concentration. Our model also predicts a positive correlation between the elongated morphology of the invading cells and the alignment of fibers in the matrix, suggesting that cell polarization is directly proportional to the stiffness and alignment of the matrix. In contrast, cells in nonfibrous matrices are found to be rounded and not polarized, underscoring the key role played by the nonlinear mechanics of fibrous matrices. Importantly, our model shows that mechanical principles mediated by the contractility of the cells and the nonlinearity of the ECM behavior play a crucial role in determining the phenotype of the cell invasion.

Simulating Heterogeneous Tumor Cell Populations

Certain tumor phenomena, like metabolic heterogeneity and local stable regions of chronic hypoxia, signify a tumor's resistance to therapy. Although recent research has shed light on the intracellular mechanisms of cancer metabolic reprogramming, little is known about how tumors become metabolically heterogeneous or chronically hypoxic, namely the initial conditions and spatiotemporal dynamics that drive these cell population conditions. To study these aspects, we developed a minimal, spatially-resolved simulation framework for modeling tissue-scale mixed populations of cells based on diffusible particles the cells consume and release, the concentrations of which determine their behavior in arbitrarily complex ways, and on stochastic reproduction. We simulate cell populations that self-sort to facilitate metabolic symbiosis, that grow according to tumor-stroma signaling patterns, and that give rise to stable local regions of chronic hypoxia near blood vessels. We raise two novel questions in the context of these results: (1) How will two metabolically symbiotic cell subpopulations self-sort in the presence of glucose, oxygen, and lactate gradients? We observe a robust pattern of alternating striations. (2) What is the proper time scale to observe stable local regions of chronic hypoxia? We observe the stability is a function of the balance of three factors related to O2-diffusion rate, local vessel release rate, and viable and hypoxic tumor cell consumption rate. We anticipate our simulation framework will help researchers design better experiments and generate novel hypotheses to better understand dynamic, emergent whole-tumor behavior.

"Active" drops as phantom models for living cells: a mesoscopic particle-based approach

Drops and biological cells share some morphological features and visco-elastic properties. The modelling of drops by mesoscopic non-atomistic models has been carried out to a high degree of success in recent years. We extend such treatment and discuss a simple, drop-like model to describe the interactions of the outer layer of cells with the surfaces of materials. Cells are treated as active mechanical objects that are able to generate adhesion forces. They appear with their true size and are made of "parcels of fluids" or beads. The beads are described by (very) few quantities/parameters related to fundamental chemical forces such as hydrophilicity and lipophilicity that represent an average of the properties of a patch of material or an area of the cell(s) surface. The investigation of adhesion dynamics, motion of individual cells, and the collective behavior of clusters of cells on materials is possible. In the simulations, the drops become active soft matter objects and different from regular droplets they do not fuse when in contact, their trajectories are not Brownian, and they can be forced "to secrete" molecules, to name some of the properties targeted by the modeling. The behavior that emerges from the simulations allows ascribing some cell properties to their mechanics, which are related to their biological features.

A cell-centered, agent-based framework that enables flexible environment granularities

BACKGROUND

Mechanistic explanations of cell-level phenomena typically adopt an observer perspective. Explanations developed from a cell's perspective may offer new insights. Agent-based models lend themselves to model from an individual perspective, and existing agent-based models generally utilize a regular lattice-based environment. A framework which utilizes a cell's perspective in an off-lattice environment could improve the overall understanding of biological phenomena.

RESULTS

We present an agent-based, discrete event framework, with a demonstrative focus on biomimetic agents. The framework was first developed in 2-dimensions and then extended, with a subset of behaviors, to 3-dimensions. The framework is expected to facilitate studies of more complex biological phenomena through exploitation of a dynamic Delaunay and Voronoi off-lattice environment. We used the framework to model biological cells and to specifically demonstrate basic biological cell behaviors in two- and three-dimensional space. Potential use cases are highlighted, suggesting the utility of the framework in various scenarios.

CONCLUSIONS

The framework presented in this manuscript expands on existing cell- and agent-centered methods by offering a new perspective in an off-lattice environment. As the demand for biomimetic models grows, the demand for new methods, such as the presented Delaunay and Voronoi framework, is expected to increase.

The Importance of Neighborhood Scheme Selection in Agent-based Tumor Growth Modeling

Modeling tumor growth has proven a very challenging problem, mainly due to the fact that tumors are highly complex systems that involve dynamic interactions spanning multiple scales both in time and space. The desire to describe interactions in various scales has given rise to modeling approaches that use both continuous and discrete variables, known as hybrid approaches. This work refers to a hybrid model on a 2D square lattice focusing on cell movement dynamics as they play an important role in tumor morphology, invasion and metastasis and are considered as indicators for the stage of malignancy used for early prognosis and effective treatment. Considering various distributions of the microenvironment, we explore how Neumann vs. Moore neighborhood schemes affects tumor growth and morphology. The results indicate that the importance of neighborhood selection is critical under specific conditions that include i) increased hapto/chemo-tactic coefficient, ii) a rugged microenvironment and iii) ECM degradation.

A continuum approximation to an off-lattice individual-cell based model of cell migration and adhesion

Cell-cell adhesion plays a key role in the collective migration of cells and in determining correlations in the relative cell positions and velocities. Recently, it was demonstrated that off-lattice individual cell based models (IBMs) can accurately capture the correlations observed experimentally in a migrating cell population. However, IBMs are often computationally expensive and difficult to analyse mathematically. Traditional continuum-based models, in contrast, are amenable to mathematical analysis and are computationally less demanding, but typically correspond to a mean-field approximation of cell migration and so ignore cell-cell correlations. In this work, we address this problem by using an off-lattice IBM to derive a continuum approximation which does take into account correlations. We furthermore show that a mean-field approximation of the off-lattice IBM leads to a single partial integro-differential equation of the same form as proposed by Sherratt and co-workers to model cell adhesion. The latter is found to be only effective at approximating the ensemble averaged cell number density when mechanical interactions between cells are weak. In contrast, the predictions of our novel continuum model for the time-evolution of the ensemble cell number density distribution and of the density-density correlation function are in close agreement with those obtained from the IBM for a wide range of mechanical interaction strengths. In particular, we observe 'front-like' propagation of cells in simulations using both our IBM and our continuum model, but not in the continuum model simulations obtained using the mean-field approximation.

Influence of the microenvironment on cell fate determination and migration

Several critical cell functions are influenced not only by internal cellular machinery but also by external mechanical and biochemical cues from the surrounding microenvironment. Slight changes to the microenvironment can result in dramatic changes to the cell's phenotype; for example, a change in the nutrients or pH of a tumor microenvironment can result in increased tumor metastasis. While cellular fate and the regulators of cell fate have been studied in detail for several decades now, our understanding of the extracellular regulators remains qualitative and far from comprehensive. In this review, we discuss the microenvironment influence on cell fate in terms of adhesion, migration, and differentiation and focus on both developments in experimental and computation tools to analyze cellular fate.

Evolution of intratumoral phenotypic heterogeneity: the role of trait inheritance

A tumour is a heterogeneous population of cells that competes for limited resources. In the clinic, we typically probe the tumour by biopsy, and then characterize it by the dominant genetic clone. But genotypes are only the first link in the chain of hierarchical events that leads to a specific cell phenotype. The relationship between genotype and phenotype is not simple, and the so-called genotype to phenotype map is poorly understood. Many genotypes can produce the same phenotype, so genetic heterogeneity may not translate directly to phenotypic heterogeneity. We therefore choose to focus on the functional endpoint, the phenotype as defined by a collection of cellular traits (e.g. proliferative and migratory ability). Here, we will examine how phenotypic heterogeneity evolves in space and time and how the way in which phenotypes are inherited will drive this evolution. A tumour can be thought of as an ecosystem, which critically means that we cannot just consider it as a collection of mutated cells but more as a complex system of many interacting cellular and microenvironmental elements. At its simplest, a growing tumour with increased proliferation capacity must compete for space as a limited resource. Hypercellularity leads to a contact-inhibited core with a competitive proliferating rim. Evolution and selection occurs, and an individual cell's capacity to survive and propagate is determined by its combination of traits and interaction with the environment. With heterogeneity in phenotypes, the clone that will dominate is not always obvious as there are both local interactions and global pressures. Several combinations of phenotypes can coexist, changing the fitness of the whole. To understand some aspects of heterogeneity in a growing tumour, we build an off-lattice agent-based model consisting of individual cells with assigned trait values for proliferation and migration rates. We represent heterogeneity in these traits with frequency distributions and combinations of traits with density maps. How the distributions change over time is dependent on how traits are passed on to progeny cells, which is our main enquiry. We bypass the translation of genetics to behaviour by focusing on the functional end result of inheritance of the phenotype combined with the environmental influence of limited space.

How changes in extracellular matrix mechanics and gene expression variability might combine to drive cancer progression

Changes in extracellular matrix (ECM) structure or mechanics can actively drive cancer progression; however, the underlying mechanism remains unknown. Here we explore whether this process could be mediated by changes in cell shape that lead to increases in genetic noise, given that both factors have been independently shown to alter gene expression and induce cell fate switching. We do this using a computer simulation model that explores the impact of physical changes in the tissue microenvironment under conditions in which physical deformation of cells increases gene expression variability among genetically identical cells. The model reveals that cancerous tissue growth can be driven by physical changes in the microenvironment: when increases in cell shape variability due to growth-dependent increases in cell packing density enhance gene expression variation, heterogeneous autonomous growth and further structural disorganization can result, thereby driving cancer progression via positive feedback. The model parameters that led to this prediction are consistent with experimental measurements of mammary tissues that spontaneously undergo cancer progression in transgenic C3(1)-SV40Tag female mice, which exhibit enhanced stiffness of mammary ducts, as well as progressive increases in variability of cell-cell relations and associated cell shape changes. These results demonstrate the potential for physical changes in the tissue microenvironment (e.g., altered ECM mechanics) to induce a cancerous phenotype or accelerate cancer progression in a clonal population through local changes in cell geometry and increased phenotypic variability, even in the absence of gene mutation.

Hybrid model of erythropoiesis

A hybrid model of cell dynamics is presented. It is illustrated by model examples and applied to study erythropoiesis (red blood cell production). In this approach, cells are considered as discrete objects while intra-cellular proteins and extra-cellular biochemical substances are described with continuous models. Spatial organization of erythropoiesis occurring in specific structures of the bone marrow, called erythroblastic island, is investigated.

Mathematical modelling of spatio-temporal glioma evolution

Background

Gliomas are the most common types of brain cancer, well known for their aggressive proliferation and the invasive behavior leading to a high mortality rate. Several mathematical models have been developed for identifying the interactions between glioma cells and tissue microenvironment, which play an important role in the mechanism of the tumor formation and progression.

Methods

Building and expanding on existing approaches, this paper develops a continuous three-dimensional model of avascular glioma spatio-temporal evolution. The proposed spherical model incorporates the interactions between the populations of four different glioma cell phenotypes (proliferative, hypoxic, hypoglychemic and necrotic) and their tissue microenvironment, in order to investigate how they affect tumor growth and invasion in an isotropic and homogeneous medium. The model includes two key variables involved in the proliferation and invasion processes of cancer cells; i.e. the extracellular matrix and the matrix-degradative enzymes concentrations inside the tumor and its surroundings. Additionally, the proposed model focuses on innovative features, such as the separate and independent impact of two vital nutrients, namely oxygen and glucose, in tumor growth, leading to the formation of cell populations with different metabolic profiles. The model implementation takes under consideration the variations of particular factors, such as the local cell proliferation rate, the variable conversion rates of cells from one category to another and the nutrient-dependent thresholds of conversion. All model variables (cell densities, ingredients concentrations) are continuous and described by reaction-diffusion equations.

Results

Several simulations were performed using combinations of growth and invasion rates, for different evolution times. The model results were evaluated by medical experts and validated on experimental glioma models available in the literature, revealing high agreement between simulated and experimental results.

Conclusions

Based on the experimental validation, as well as the evaluation by clinical experts, the proposed model may provide an essential tool for the patient-specific simulation of different tumor evolution scenarios and reliable prognosis of glioma spatio-temporal progression.

A computational model for collective cellular motion in three dimensions: general framework and case study for cell pair dynamics

Cell migration in healthy and diseased systems is a combination of single and collective cell motion. While single cell motion has received considerable attention, our understanding of collective cell motion remains elusive. A new computational framework for the migration of groups of cells in three dimensions is presented, which focuses on the forces acting at the microscopic scale and the interactions between cells and their extracellular matrix (ECM) environment. Cell-cell adhesion, resistance due to the ECM and the factors regulating the propulsion of each cell through the matrix are considered. In particular, our approach emphasizes the role of receptors that mediate cell-cell and cell-matrix interactions, and examines how variation in their properties induces changes in cellular motion. As an important case study, we analyze two interacting cells. Our results show that the dynamics of cell pairs depends on the magnitude and the stochastic nature of the forces. Stronger intercellular stability is generally promoted by surface receptors that move. We also demonstrate that matrix resistance, cellular stiffness and intensity of adhesion contribute to migration behaviors in different ways, with memory effects present that can alter pair motility. If adhesion weakens with time, our findings show that cell pair break-up depends strongly on the way cells interact with the matrix. Finally, the motility for cells in a larger cluster (size 50 cells) is examined to illustrate the full capabilities of the model and to stress the role of cellular pairs in complex cellular structures. Overall, our framework shows how properties of cells and their environment influence the stability and motility of cellular assemblies. This is an important step in the advancement of the understanding of collective motility, and can contribute to knowledge of complex biological processes involving migration, aggregation and detachment of cells in healthy and diseased systems.

Tumor cell migration in complex microenvironments

Tumor cell migration is essential for invasion and dissemination from primary solid tumors and for the establishment of lethal secondary metastases at distant organs. In vivo and in vitro models enabled identification of different factors in the tumor microenvironment that regulate tumor progression and metastasis. However, the mechanisms by which tumor cells integrate these chemical and mechanical signals from multiple sources to navigate the complex microenvironment remain poorly understood. In this review, we discuss the factors that influence tumor cell migration with a focus on the migration of transformed carcinoma cells. We provide an overview of the experimental and computational methods that allow the investigation of tumor cell migration, and we highlight the benefits and shortcomings of the various assays. We emphasize that the chemical and mechanical stimulus paradigms are not independent and that crosstalk between them motivates the development of new assays capable of applying multiple, simultaneous stimuli and imaging the cellular migratory response in real-time. These next-generation assays will more closely mimic the in vivo microenvironment to provide new insights into tumor progression, inform techniques to control tumor cell migration, and render cancer more treatable.

A novel method for simulating the extracellular matrix in models of tumour growth

A novel hybrid continuum-discrete model to simulate tumour growth on a cellular scale is proposed. The lattice-based spatiotemporal model consists of reaction-diffusion equations that describe interactions between cancer cells and their microenvironment. The fundamental ingredients that are typically considered are the nutrient concentration, the extracellular matrix (ECM), and matrix degrading enzymes (MDEs). The in vivo processes are very complex and occur on different levels. This in turn leads to huge computational costs. The main contribution of the present work is therefore to describe the processes on the basis of simplified mathematical approaches, which, at the same time, depict realistic results to understand the biological processes. In this work, we discuss if we have to simulate the MDE or if the degraded matrix can be estimated directly with respect to the cancer cell distribution. Additionally, we compare the results for modelling tumour growth using the common and our simplified approach, thereby demonstrating the advantages of the proposed method. Therefore, we introduce variations of the positioning of the nutrient delivering blood vessels and use different initializations of the ECM. We conclude that the novel method, which does not explicitly model the matrix degrading enzymes, provides means for a straightforward and fast implementation for modelling tumour growth.

Selection of Mesenchymal-Like Metastatic Cells in Primary Tumors - An in silico Investigation

In order to metastasize, cancer cells must undergo phenotypic transition from an anchorage-dependent form to a motile form via a process referred to as epithelial to mesenchymal transition. It is currently unclear whether metastatic cells emerge late during tumor progression by successive accumulation of mutations, or whether they derive from distinct cell populations already present during the early stages of tumorigenesis. Similarly, the selective pressures that drive metastasis are poorly understood. Selection of cancer cells with increased proliferative capacity and enhanced survival characteristics may explain how some transformations promote a metastatic phenotype. However, it is difficult to explain how cancer cells that disseminate can emerge due to such selective pressure, since these cells usually remain dormant for prolonged periods of time. In the current study, we have used in silico modeling and simulation to investigate the hypothesis that mesenchymal-like cancer cells evolve during the early stages of primary tumor development, and that these cells exhibit survival and proliferative advantages within the tumor microenvironment. In an agent-based tumor microenvironment model, cancer cell agents with distinct sets of attributes governing nutrient consumption, proliferation, apoptosis, random motility, and cell adhesion were allowed to compete for space and nutrients. These simulation data indicated that mesenchymal-like cancer cells displaying high motility and low adhesion proliferate more rapidly and display a survival advantage over epithelial-like cancer cells. Furthermore, the presence of mesenchymal-like cells within the primary tumor influences the macroscopic properties, emergent morphology, and growth rate of tumors.

Invasion from a cell aggregate--the roles of active cell motion and mechanical equilibrium

Cell invasion from an aggregate into a surrounding extracellular matrix (ECM) is an important process during development disease, e.g., vascular network assembly or tumor progression. To describe the behavior emerging from autonomous cell motility, cell-cell adhesion and contact guidance by ECM filaments, we propose a suitably modified cellular Potts model. We consider an active cell motility process in which internal polarity is governed by a positive feedback from cell displacements, a mechanism that can result in highly persistent motion when constrained by an oriented ECM structure. The model allows us to explore the interplay between haptotaxis, matrix degradation and active cell movement. We show that for certain conditions the cells are able to both invade the ECM and follow the ECM tracks. Furthermore, we argue that enforcing mechanical equilibrium within a bulk cell mass is of key importance in multicellular simulations.

Classifying general nonlinear force laws in cell-based models via the continuum limit

Although discrete cell-based frameworks are now commonly used to simulate a whole range of biological phenomena, it is typically not obvious how the numerous different types of model are related to one another, nor which one is most appropriate in a given context. Here we demonstrate how individual cell movement on the discrete scale modeled using nonlinear force laws can be described by nonlinear diffusion coefficients on the continuum scale. A general relationship between nonlinear force laws and their respective diffusion coefficients is derived in one spatial dimension and, subsequently, a range of particular examples is considered. For each case excellent agreement is observed between numerical solutions of the discrete and corresponding continuum models. Three case studies are considered in which we demonstrate how the derived nonlinear diffusion coefficients can be used to (a) relate different discrete models of cell behavior; (b) derive discrete, intercell force laws from previously posed diffusion coefficients, and (c) describe aggregative behavior in discrete simulations.

Normal morphogenesis of epithelial tissues and progression of epithelial tumors

Epithelial cells organize into various tissue architectures that largely maintain their structure throughout the life of an organism. For decades, the morphogenesis of epithelial tissues has fascinated scientists at the interface of cell, developmental, and molecular biology. Systems biology offers ways to combine knowledge from these disciplines by building integrative models that are quantitative and predictive. Can such models be useful for gaining a deeper understanding of epithelial morphogenesis? Here, we take inventory of some recurring themes in epithelial morphogenesis that systems approaches could strive to capture. Predictive understanding of morphogenesis at the systems level would prove especially valuable for diseases such as cancer, where epithelial tissue architecture is profoundly disrupted.

Investigation of bone resorption within a cortical basic multicellular unit using a lattice-based computational model

In this paper we develop a lattice-based computational model focused on bone resorption by osteoclasts in a single cortical basic multicellular unit (BMU). Our model takes into account the interaction of osteoclasts with the bone matrix, the interaction of osteoclasts with each other, the generation of osteoclasts from a growing blood vessel, and the renewal of osteoclast nuclei by cell fusion. All these features are shown to strongly influence the geometrical properties of the developing resorption cavity including its size, shape and progression rate, and are also shown to influence the distribution, resorption pattern and trajectories of individual osteoclasts within the BMU. We demonstrate that for certain parameter combinations, resorption cavity shapes can be recovered from the computational model that closely resemble resorption cavity shapes observed from microCT imaging of human cortical bone.

Modelling acidosis and the cell cycle in multicellular tumour spheroids

A partial differential equation model is developed to understand the effect that nutrient and acidosis have on the distribution of proliferating and quiescent cells and dead cell material (necrotic and apoptotic) within a multicellular tumour spheroid. The rates of cell quiescence and necrosis depend upon the local nutrient and acid concentrations and quiescent cells are assumed to consume less nutrient and produce less acid than proliferating cells. Analysis of the differences in nutrient consumption and acid production by quiescent and proliferating cells shows low nutrient levels do not necessarily lead to increased acid concentration via anaerobic metabolism. Rather, it is the balance between proliferating and quiescent cells within the tumour which is important; decreased nutrient levels lead to more quiescent cells, which produce less acid than proliferating cells. We examine this effect via a sensitivity analysis which also includes a quantification of the effect that nutrient and acid concentrations have on the rates of cell quiescence and necrosis.

Spatiotemporal effects of a possible chemorepulsion of tumor cells by immune system effectors

Tumor-immune system interplay is a landmark of tumor development, and chemotactic attraction of immune system cells towards a tumor is a landmark of immune control. Since tumor cells are capable of chemotactic and chemorepulsive motion, based on a number of analogies between the behavior of tumor cells and that of bacteria, and on the production of potentially chemorepulsive semaphorins by immune systems effectors, we propose here the possibility of chemorepulsion of tumor cells, elicited by chemicals produced by immune system effectors such as macrophages and cytotoxic T lymphocytes. To study the effects of this hypothesized phenomenon within the framework of the interplay of neoplasias with the innate and adaptive immune system, we appropriately extend two well-known models: the tumor-macrophage models by Owen and Sherratt [Owen and Sherratt, J. Theor. Biol., 1998] and the Matzavinos-Chaplain model of tumor-CTL interaction [Matzavinos, Chaplain and Kuznetsov, Math. Med. Biol., 2004]. Our simulations suggest that this mechanism might allow a faster expansion of tumors, and in the concluding remarks we envisage a new possible route of immunoevasion.

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.