Modeling of self-organized avascular tumor growth with a hybrid cellular automaton

Pattern formation in multicellular spheroids is addressed with a hybrid lattice-gas cellular automaton model. Multicellular spheroids serve as experimental model system for the study of avascular tumor growth. Typically, multicellular spheroids consist of a necrotic core surrounded by rings of quiescent and proliferating tumor cells, respectively. Furthermore, after an initial exponential growth phase further spheroid growth is significantly slowed down even if further nutrient is supplied. The cellular automaton model explicitly takes into account mitosis, apoptosis and necrosis as well as nutrient consumption and a diffusible signal that is emitted by cells becoming necrotic. All cells follow identical interaction rules. The necrotic signal induces a chemotactic migration of tumor cells towards maximal signal concentrations. Starting from a small number of tumor cells automaton simulations exhibit the self-organized formation of a layered structure consisting of a necrotic core, a ring of quiescent tumor cells and a thin outer ring of proliferating tumor cells.     

A single-cell-based model of tumor growth in vitro: monolayers and spheroids

To what extent the growth dynamics of tumors is controlled by nutrients, biomechanical forces and other factors at different stages and in different environments is still largely unknown. Here we present a biophysical model to study the spatio-temporal growth dynamics of two-dimensional tumor monolayers and three-dimensional tumor spheroids as a complementary tool to in vitro experiments. Within our model each cell is represented as an individual object and parametrized by cell-biophysical and cell-kinetic parameters that can all be experimentally determined. Hence our modeling strategy allows us to study which mechanisms on the microscopic level of individual cells may affect the macroscopic properties of a growing tumor. We find the qualitative growth kinetics and patterns at early growth stages to be remarkably robust. Quantitative comparisons between computer simulations using our model and published experimental observations on monolayer cultures suggest a biomechanically-mediated form of growth inhibition during the experimentally observed transition from exponential to sub-exponential growth at sufficiently large tumor sizes. Our simulations show that the same transition during the growth of avascular tumor spheroids can be explained largely by the same mechanism. Glucose (or oxygen) depletion seems to determine mainly the size of the necrotic core but not the size of the tumor. We explore the consequences of the suggested biomechanical form of contact inhibition, in order to permit an experimental test of our model. Based on our findings we propose a phenomenological growth law in early expansion phases in which specific biological small-scale processes are subsumed in a small number of effective parameters.     

An on-lattice agent-based Monte Carlo model simulating the growth kinetics of multicellular tumor spheroids

PurposeTo develop an on-lattice agent-based model describing the growth of multicellular tumor spheroids using simple Monte Carlo tools.MethodsCells are situated on the vertices of a cubic grid. Different cell states (proliferative, hypoxic or dead) and cell evolution rules, driven by 10 parameters, and the effects of the culture medium are included. About twenty spheroids of MCF-7 human breast cancer were cultivated and the experimental data were used for tuning the model parameters.ResultsSimulated spheroids showed adequate sizes of the necrotic nuclei and of the hypoxic and proliferative cell phases as a function of the growth time, mimicking the overall characteristics of the experimental spheroids. The relation between the radii of the necrotic nucleus and the whole spheroid obtained in the simulations was similar to the experimental one and the number of cells, as a function of the spheroid volume, was well reproduced. The statistical variability of the Monte Carlo model described the whole volume range observed for the experimental spheroids. Assuming that the model parameters vary within Gaussian distributions it was obtained a sample of spheroids that reproduced much better the experimental findings.ConclusionsThe model developed allows describing the growth of in vitro multicellular spheroids and the experimental variability can be well reproduced. Its flexibility permits to vary both the agents involved and the rules that govern the spheroid growth. More general situations, such as, e. g., tumor vascularization, radiotherapy effects on solid tumors, or the validity of the tumor growth mathematical models can be studied.     

Formulation and numerical simulations of a continuum model of avascular tumor growth

In this paper we present a continuum mathematical model for a multicellular spheroid that mimics the micro-environment within avascular tumor growth. The model consists of a coupled system of non-linear convection-diffusion-reaction equations. This system is solved using a previously developed conservative Galerkin characteristics method. In the model considered, there are three cell types: the proliferative cells, the quiescent non-dividing cells which stay in the G₀ phase of the cell cycle and the necrotic cells. The model includes viable cell diffusion, diffusion of cellular material and the removal of necrotic cells. We assume that the nutrients diffuse passively and are consumed by the proliferative and quiescent tumor cells depending on the availability of resources (oxygen, glucose, etc.). The numerical simulations are performed using different sets of parameters, including biologically realistic ones, to explore the effects of each of these model parameters on reaching the steady state. The present results, taken together with those reported earlier, indicate that the removal of necrotic cells and the diffusion of cellular material have significant effects on the steady state, reflecting growth saturation, the number of viable cells, and the spheroid size.     

Predicting the growth of glioblastoma multiforme spheroids using a multiphase porous media model

Tumor spheroids constitute an effective in vitro tool to investigate the avascular stage of tumor growth. These three-dimensional cell aggregates reproduce the nutrient and proliferation gradients found in the early stages of cancer and can be grown with a strict control of their environmental conditions. In the last years, new experimental techniques have been developed to determine the effect of mechanical stress on the growth of tumor spheroids. These studies report a reduction in cell proliferation as a function of increasingly applied stress on the surface of the spheroids. This work presents a specialization for tumor spheroid growth of a previous more general multiphase model. The equations of the model are derived in the framework of porous media theory, and constitutive relations for the mass transfer terms and the stress are formulated on the basis of experimental observations. A set of experiments is performed, investigating the growth of U-87MG spheroids both freely growing in the culture medium and subjected to an external mechanical pressure induced by a Dextran solution. The growth curves of the model are compared to the experimental data, with good agreement for both the experimental settings. A new mathematical law regulating the inhibitory effect of mechanical compression on cancer cell proliferation is presented at the end of the paper. This new law is validated against experimental data and provides better results compared to other expressions in the literature.     

A cellular automata model of tumor-immune system interactions

We present a hybrid cellular automata-partial differential equation model of moderate complexity to describe the interactions between a growing tumor next to a nutrient source and the immune system of the host organism. The model allows both temporal and two-dimensional spatial evolution of the system under investigation and is comprised of biological cell metabolism rules derived from both the experimental and mathematical modeling literature. We present numerical simulations that display behaviors which are qualitatively similar to those exhibited in tumor-immune system interaction experiments. These include spherical tumor growth, stable and unstable oscillatory tumor growth, satellitosis and tumor infiltration by immune cells. Finally, the relationship between these different growth regimes and key system parameters is discussed.     

Variations in tumor cell growth rates and metabolism with oxygen concentration, glucose concentration, and extracellular pH.

Tumors and multicellular tumor spheroids can develop gradients in oxygen concentration, glucose concentration, and extracellular pH as they grow. In order to calculate these gradients and assess their impact on tumor growth, it is necessary to quantify the effect of these variables on tumor cell metabolism and growth. In this work, the oxygen consumption rates, glucose consumption rates, and growth rates of EMT6/Ro mouse mammary tumor cells were measured at a variety of oxygen concentrations, glucose concentrations, and extracellular pH levels. At an extracellular pH of 7.25, the oxygen consumption rate of EMT6/Ro cells increased by nearly a factor of 2 as the glucose concentration was decreased from 5.5 mM to 0.4 mM. This effect of glucose concentration on oxygen consumption rate, however, was slight at an extracellular pH of 6.95 and disappeared completely at an extracellular pH of 6.60. The glucose consumption rate of EMT6/Ro cells increased by roughly 40% when the oxygen concentration was reduced from 0.21 mM to 0.023 mM and decreased by roughly 60% when the extracellular pH was decreased from 7.25 to 6.95. The growth rate of EMT6/Ro cells decreased with decreasing oxygen concentration and extracellular pH; however, severe conditions were required to stop cell growth (0.0082 mM oxygen and an extracellular pH of 6.60). Empirical correlations were developed from these data to express EMT6/Ro cell growth rates, oxygen consumption rates, and glucose consumption rates, as functions of oxygen concentration, glucose concentration, and extracellular pH. These empirical correlations make it possible to mathematically model the gradients in oxygen concentration, glucose concentration, and extracellular pH in EMT6/Ro multicellular spheroids by solution of the diffusion/reaction equations. Computations such as these, along with oxygen and pH microelectrode measurements in EMT6/Ro multicellular spheroids, indicated that nutrient concentration and pH levels in the inner regions of spheroids were low enough to cause significant changes in nutrient consumption rates and cell growth rates. However, pH and oxygen concentrations measured or calculated in EMT6/Ro spheroids where quiescent cells have been observed were not low enough to cause the cessation of cell growth, indicating that the observed quiescence must have been due to factors other than acidic pH, oxygen depletion, or glucose depletion.     

Modelling the formation of necrotic regions in avascular tumours

The mechanisms underlying the formation of necrotic regions within avascular tumours are not well understood. In this paper, we examine the relative roles of nutrient deprivation and of cell death, from both the proliferating phase of the cell cycle via apoptosis and from the quiescent phase via necrosis, in changing the structure within multicellular tumour spheroids and particularly the accumulation of dead cell material in the centre. A mathematical model is presented and studied that accounts for nutrient diffusion, changes in cell cycling rates, the two different routes to cell death as well as active motion of cells and passive motion of the dead cell material. In studying the accumulation of dead cell matter we do not distinguish between the route by which each was formed. The resulting mathematical model is examined for a number of scenarios. Results show that in many cases the size of the necrotic core is closely correlated with low levels in nutrient concentration. However, in certain cases, particularly where the rate of necrosis is large, the resulting necrotic core can lead to regions of non-negligible nutrient concentration-dependent upon the mode of cell death.     

Diffusion regulated growth characteristics of a spherical prevascular carcinoma

Recently a mathematical model of the prevascular phases of tumor growth by diffusion has been investigated (S. A. Maggelakis and J. A. Adam,Math. Comput. Modeling, in press). In this paper we examine in detail the results and implications of that mathematical model, particularly in the light of recent experimental work carried out on multicellular spheroids. The overall growth characteristics are determined in the present model by four parameters:Q, γ, b, andδ, which depend on information about inhibitor production rates, oxygen consumption rates, volume loss and cell proliferation rates, and measures of the degree of non-uniformity of the various diffusion processes that take place. The integro-differential growth equation is solved for the outer spheroid radiusR ₀(t) and three related inner radii subject to the solution of the governing time-independent diffusion equations (under conditions of diffusive equilibrium) and the appropriate boundary conditions. Hopefully, future experimental work will enable reasonable bounds to be placed on parameter values referred to in this model: meanwhile, specific experimentally-provided initial data can be used to predict subsequent growth characteristics ofin vitro multicellular spheroids. This will be one objective of future studies.     

A model for the growth of multicellular spheroids

Abstract. Based on biological observations and the basic physical properties of tri-dimensional structures, a mathematical expression is derived to relate the growth rate of multicellular spheroids to some easily measurable parameters. This model involves properties both of the individual cells and of the spheroid structure, such as the cell doubling time in monolayer, the rate of cell shedding from the spheroid and the depth of the external rim of cycling cells. The derived growth equation predicts a linear expansion of the spheroid diameter with time. The calculated growth rate for a number of spheroid cell types is in good agreement with experimental data. The model provides a simple and practical view of growth control in spheroids, and is further adapted to include parameters presumably responsible for the growth saturation in large spheroids.     

ATP concentrations in multicellular tumor spheroids assessed by single photon imaging and quantitative bioluminescence

To evaluate the interrelationship among the cellular energy status and the development of necrosis in tumor microregions, local ATP concentrations and the extent of necrosis were determined in multicellular tumor spheroids, i.e., in spherical tumor cell aggregates. The spheroids were grown in rotated suspension cultures using EMT6 cells that were derived from a murine mammary sarcoma. The distribution of viable and necrotic cell areas was assessed by histological investigations. The regional distribution of ATP concentrations was measured with a novel technique using quantitative bioluminescence and single photon imaging. This method makes it possible to determine ATP concentrations in absolute terms with a spatial resolution at the level of a single cell. The results show that ATP concentrations in the center of EMT6 spheroids decrease from values of 1.0 to 1.5 mM in small spheroids with 300 microns in diameter to values close to or at the background level in 750 microns spheroids. Necrosis was detectable in spheroids larger than 300 microns, and virtually no spheroid without necrosis was found at sizes larger than 600 microns. Since the emergence of central necrosis precedes the drop in ATP to undetectably low values, the data suggest that energy metabolism is not or not directly involved in the development of necrosis in tumor spheroids under the growth conditions investigated.     

Live cell division dynamics monitoring in 3D large spheroid tumor models using light sheet microscopy

Background

Multicellular tumor spheroids are models of increasing interest for cancer and cell biology studies. They allow considering cellular interactions in exploring cell cycle and cell division mechanisms. However, 3D imaging of cell division in living spheroids is technically challenging and has never been reported.

Results

Here, we report a major breakthrough based on the engineering of multicellular tumor spheroids expressing an histone H2B fluorescent nuclear reporter protein, and specifically designed sample holders to monitor live cell division dynamics in 3D large spheroids using an home-made selective-plane illumination microscope.

Conclusions

As illustrated using the antimitotic drug, paclitaxel, this technological advance paves the way for studies of the dynamics of cell divion processes in 3D and more generally for the investigation of tumor cell population biology in integrated system as the spheroid model.     

Differential sensitivity to short-chain ceramide analogues of human intestinal carcinoma cells grown in tumor spheroids versus monolayer culture

The cytotoxic activity of short-chain (C(2)) ceramide was evaluated in human intestinal carcinoma cells grown as multicellular tumor spheroids versus the same cells cultured as monolayers under closely comparable conditions. A decrease in cell number was seen in monolayer cultures of HT-29, Caco-2, and HRT-18 cells, with an EC(50) (concentration for half-maximal toxicity) of between 13 and 23 microM. However, when the same cells were grown in the multicellular spheroid format, C(2) was markedly less potent in reducing cell number, with an EC(50) of between 44 and 63 microM, representing a 1.9- to 4.9-fold decrease in its potency. The chemotherapeutic agents 5-fluorouracil and cisplatin were equally potent against spheroids and monolayer cultures, indicating that although drug access is a problem in conventionally grown tumor spheroids it is not a problem for spheroids grown under the conditions used in this study. Our results suggest that although ceramide is capable of inducing cell death in intestinal carcinoma cells grown in spheroid culture, its cellular toxicity is constrained by influences that are independent of drug access and may be the consequence of the altered cellular relationships. Carcinoma cell populations show an intrinsically decreased responsiveness to the effects of ceramide when they are grown in a three-dimensional culture format.     

The migration of cells in multicell tumor spheroids

A mathematical model is proposed to explain the observed internalization of microspheres and 3H-thymidine labelled cells in steady-state multicellular spheroids. The model uses the conventional ideas of nutrient diffusion and consumption by the cells. In addition, a very simple model of the progress of the cells through the cell cycle is considered. Cells are divided into two classes, those proliferating (being in G1, S, G2 or M phases) and those that are quiescent (being in G0). Furthermore, the two categories are presumed to have different chemotactic responses to the nutrient gradient. The model accounts for the spatial and temporal variations in the cell categories together with mitosis, conversion between categories and cell death. Numerical solutions demonstrate that the model predicts the behavior similar to existing models but has some novel effects. It allows for spheroids to approach a steady-state size in a non-monotonic manner, it predicts self-sorting of the cell classes to produce a thin layer of rapidly proliferating cells near the outer surface and significant numbers of cells within the spheroid stalled in a proliferating state. The model predicts that overall tumor growth is not only determined by proliferation rates but also by the ability of cells to convert readily between the classes. Moreover, the steady-state structure of the spheroid indicates that if the outer layers are removed then the tumor grows quickly by recruiting cells stalled in a proliferating state. Questions are raised about the chemotactic response of cells in differing phases and to the dependency of cell cycle rates to nutrient levels.     

Spatio-temporal modeling of nanoparticle delivery to multicellular tumor spheroids

The inefficiency of nanoparticle penetration in tissues limits the therapeutic efficacy of such formulations for cancer applications. Recent work has indicated that modulation of tissue architecture with enzymes such as collagenase significantly increases macromolecule delivery. In this study we developed a mathematical model of nanoparticle penetration into multicellular spheroids that accounts for radially dependent changes in tumor architecture, as represented by the volume fraction of tissue accessible to nanoparticle diffusion. Parameters such as nanoparticle binding, internalization rate constants, and accessible volume fraction were determined experimentally. Unknown parameters of nanoparticle binding sites per cell in the spheroid and pore shape factor were determined by fitting to experimental data. The model was correlated with experimental studies of the penetration of 40 nm nanoparticles in SiHa multicellular spheroids with and without collagenase treatment and was able to accurately predict concentration profiles of nanoparticles within spheroids. The model was also used to investigate the effects of nanoparticle size. This model contributes toward the understanding of the role of tumor architecture on nanoparticle delivery efficiency.     

Modeling mechanical inhomogeneities in small populations of proliferating monolayers and spheroids

Understanding the mechanical behavior of multicellular monolayers and spheroids is fundamental to tissue culture, organism development, and the early stages of tumor growth. Proliferating cells in monolayers and spheroids experience mechanical forces as they grow and divide and local inhomogeneities in the mechanical microenvironment can cause individual cells within the multicellular system to grow and divide at different rates. This differential growth, combined with cell division and reorganization, leads to residual stress. Multiple different modeling approaches have been taken to understand and predict the residual stresses that arise in growing multicellular systems, particularly tumor spheroids. Here, we show that by using a mechanically robust agent-based model constructed with the peridynamic framework, we gain a better understanding of residual stresses in multicellular systems as they grow from a single cell. In particular, we focus on small populations of cells (1–100 s) where population behavior is highly stochastic and prior investigation has been limited. We compare the average strain energy density of cells in monolayers and spheroids using different growth and division rules and find that, on average, cells in spheroids have a higher strain energy density than cells in monolayers. We also find that cells in the interior of a growing spheroid are, on average, in compression. Finally, we demonstrate the importance of accounting for stochastic fluctuations in the mechanical environment, particularly when the cellular response to mechanical cues is nonlinear. The results presented here serve as a starting point for both further investigation with agent-based models, and for the incorporation of major findings from agent-based models into continuum scale models when explicit representation of individual cells is not computationally feasible.     

Nemosis, a novel way of fibroblast activation, in inflammation and cancer

Malignant cells when grown in suspension, as a rule, proliferate and can form spheroids that have been used as a model of tumor nodules, micrometastases and avascular tumors. In contrast, normal adherent cells cannot be stimulated to grow as multicellular aggregates. Now, recent results show that normal fibroblasts if forced to cluster (spheroid formation) do not grow but undergo a new pathway of cell activation (nemosis) leading to a massive proinflammatory, proteolytic and growth factor response. The clustering and activation are initiated by fibronectin-integrin interaction. The activated fibroblasts are able to modulate the behavior of cancer cells and, furthermore malignant cells boost this activation even further. In this model, the activation of fibroblasts terminates in programmed necrosis-like cell death. Activation of the tumor stroma, especially of fibroblasts, is of critical importance for tumor progression, although mechanisms leading to their activation are still largely uncharacterized. In summary, our results suggest that this kind of fibroblast activation (nemosis) may be involved in pathological conditions such as inflammation and cancer.     

A 3-D model of tumor progression based on complex automata driven by particle dynamics

The dynamics of a growing tumor involving mechanical remodeling of healthy tissue and vasculature is neglected in most of the existing tumor models. This is due to the lack of efficient computational framework allowing for simulation of mechanical interactions. Meanwhile, just these interactions trigger critical changes in tumor growth dynamics and are responsible for its volumetric and directional progression. We describe here a novel 3-D model of tumor growth, which combines particle dynamics with cellular automata concept. The particles represent both tissue cells and fragments of the vascular network. They interact with their closest neighbors via semi-harmonic central forces simulating mechanical resistance of the cell walls. The particle dynamics is governed by both the Newtonian laws of motion and the cellular automata rules. These rules can represent cell life-cycle and other biological interactions involving smaller spatio-temporal scales. We show that our complex automata, particle based model can reproduce realistic 3-D dynamics of the entire system consisting of the tumor, normal tissue cells, blood vessels and blood flow. It can explain phenomena such as the inward cell motion in avascular tumor, stabilization of tumor growth by the external pressure, tumor vascularization due to the process of angiogenesis, trapping of healthy cells by invading tumor, and influence of external (boundary) conditions on the direction of tumor progression. We conclude that the particle model can serve as a general framework for designing advanced multiscale models of tumor dynamics and it is very competitive to the modeling approaches presented before.