A reduction in the in situ rates of oxygen and glucose consumption of cells in EMT6/Ro spheroids during growth

Analysis of the effects of bench-scale cell culture platforms and inoculum cell concentrations on PSC aggregate formation and culture

Introduction: Human pluripotent stem cells (hPSCs) provide many opportunities for application in regenerative medicine due to their ability to differentiate into cells from all three germ layers, proliferate indefinitely, and replace damaged or dysfunctional cells. However, such cell replacement therapies require the economical generation of clinically relevant cell numbers. Whereas culturing hPSCs as a two-dimensional monolayer is widely used and relatively simple to perform, their culture as suspended three-dimensional aggregates may enable more economical production in large-scale stirred tank bioreactors. To be more relevant to this biomanufacturing, bench-scale differentiation studies should be initiated from aggregated hPSC cultures. Methods: We compared five available bench-scale platforms for generating undifferentiated cell aggregates of human embryonic stem cells (hESCs) using AggreWell™ plates, low attachment plates on an orbital shaker, roller bottles, spinner flasks, and vertical-wheel bioreactors (PBS-Minis). Thereafter, we demonstrated the incorporation of an hPSC aggregation step prior to directed differentiation to pancreatic progenitors and endocrine cells. Results and discussion: The AggreWell™ system had the highest aggregation yield. The initial cell concentrations had an impact on the size of aggregates generated when using AggreWell™ plates as well as in roller bottles. However, aggregates made with low attachment plates, spinner flasks and PBS-Minis were similar regardless of the initial cell number. Aggregate morphology was compact and relatively homogenously distributed in all platforms except for the roller bottles. The size of aggregates formed in PBS-Minis was modulated by the agitation rate during the aggregation. In all cell culture platforms, the net growth rate of cells in 3D aggregates was lower (range: -0.01-0.022 h-1) than cells growing as a monolayer (range: 0.039-0.045 h-1). Overall, this study describes operating ranges that yield high-quality undifferentiated hESC aggregates using several of the most commonly used bench-scale cell culture platforms. In all of these systems, methods were identified to obtain PSC aggregates with greater than 70% viability, and mean diameters between 60 and 260 mm. Finally, we showed the capacity of hPSC aggregates formed with PBS-Minis to differentiate into viable pancreatic progenitors and endocrine cell types.

Characterization of the Intracellular Acidity Regulation of Brain Tumor Cells and Consequences for Therapeutic Optimization of Temozolomide

A well-known feature of tumor cells is high glycolytic activity, leading to acidification of the tumor microenvironment through extensive lactate production. This acidosis promotes processes such as metastasis, aggressiveness, and invasiveness, which have been associated with a worse clinical prognosis. Moreover, the function and expression of transporters involved in regulation of intracellular pH might be altered. In this study, the capacity of tumor cells to regulate their intracellular pH when exposed to a range of pH from very acidic to basic was characterized in two glioma cell lines (F98 and U87) using a new recently published method of fluorescence imaging. Our results show that the regulation of acidity in tumors is not the same for the two investigated cell lines; U87 cells are able to reduce their intracellular acidity, whereas F98 cells do not exhibit this property. On the other hand, F98 cells show a higher level of resistance to acidity than U87 cells. Intracellular regulation of acidity appears to be highly cell-dependent, with different mechanisms activated to preserve cell integrity and function. This characterization was performed on 2D monolayer cultures and 3D spheroids. Spatial heterogeneities were exhibited in 3D, suggesting a spatially modulated regulation in this context. Based on the corpus of knowledge available in the literature, we propose plausible mechanisms to interpret our results, together with some new lines of investigation to validate our hypotheses. Our results might have implications on therapy, since the activity of temozolomide is highly pH-dependent. We show that the drug efficiency can be enhanced, depending on the cell type, by manipulating the extracellular pH. Therefore, personalized treatment involving a combination of temozolomide and pH-regulating agents can be considered.

Optimization of Size of Nanosensitizers for Antitumor Radiotherapy Using Mathematical Modeling

The efficacy of antitumor radiotherapy can be enhanced by utilizing nonradioactive nanoparticles that emit secondary radiation when activated by a primary beam. They consist of small volumes of a radiosensitizing substance embedded within a polymer layer, which is coated with tumor-specific antibodies. The efficiency of nanosensitizers relies on their successful delivery to the tumor, which depends on their size. Increasing their size leads to a higher concentration of active substance; however, it hinders the penetration of nanosensitizers through tumor capillaries, slows down their movement through the tissue, and accelerates their clearance. In this study, we present a mathematical model of tumor growth and radiotherapy with the use of intravenously administered tumor-specific nanosensitizers. Our findings indicate that their optimal size for achieving maximum tumor radiosensitization following a single injection of their fixed total volume depends on the permeability of the tumor capillaries. Considering physiologically plausible spectra of capillary pore radii, with a nanoparticle polymer layer width of 7 nm, the optimal radius of nanoparticles falls within the range of 13-17 nm. The upper value is attained when considering an extreme spectrum of capillary pores.

Optimization of antitumor radiotherapy fractionation via mathematical modeling with account of 4 R's of radiobiology

A spatially-distributed continuous mathematical model of solid tumor growth and treatment by fractionated radiotherapy is presented. The model explicitly accounts for the factors, widely referred to as 4 R's of radiobiology, which influence the efficacy of radiotherapy fractionation protocols: tumor cell repopulation, their redistribution in cell cycle, reoxygenation and repair of sublethal damage of both tumor and normal tissues. With the use of special algorithm the fractionation protocols that provide increased tumor control probability, compared to standard clinical protocol, are found for various physiologically-based values of model parameters under the constraints of fixed overall normal tissue damage and maximum admissible fractional dose. In particular, it is shown that significant gain in treatment efficacy can be achieved for tumors of low malignancy by the use of protracted hyperfractionated protocols. The optimized non-uniform protocols are characterized by gradual escalation of fractional doses in their last parts, which start after the levels of oxygen and nutrients significantly elevate throughout the tumor and accelerated tumor proliferation manifests itself, which is a well-known experimental phenomenon.

Hybrid computational models of multicellular tumour growth considering glucose metabolism

Cancer cells metabolize glucose through metabolic pathways that differ from those used by healthy and differentiated cells. In particular, tumours have been shown to consume more glucose than their healthy counterparts and to use anaerobic metabolic pathways, even under aerobic conditions. Nevertheless, scientists have still not been able to explain why cancer cells evolved to present an altered metabolism and what evolutionary advantage this might provide them. Experimental and computational models have been increasingly used in recent years to understand some of these biological questions. Multicellular tumour spheroids are effective experimental models as they replicate the initial stages of avascular solid tumour growth. Furthermore, these experiments generate data which can be used to calibrate and validate computational studies that aim to simulate tumour growth. Hybrid models are of particular relevance in this field of research because they model cells as individual agents while also incorporating continuum representations of the substances present in the surrounding microenvironment that may participate in intracellular metabolic networks as concentration or density distributions. Henceforth, in this review, we explore the potential of computational modelling to reveal the role of metabolic reprogramming in tumour growth.

Modeling and measuring glucose diffusion and consumption by colorectal cancer spheroids in hanging drops using integrated biosensors

As 3D in vitro tissue models become more pervasive, their built-in nutrient, metabolite, compound, and waste gradients increase biological relevance at the cost of analysis simplicity. Investigating these gradients and the resulting metabolic heterogeneity requires invasive and time-consuming methods. An alternative is using electrochemical biosensors and measuring concentrations around the tissue model to obtain size-dependent metabolism data. With our hanging-drop-integrated enzymatic glucose biosensors, we conducted current measurements within hanging-drop compartments hosting spheroids formed from the human colorectal carcinoma cell line HCT116. We developed a physics-based mathematical model of analyte consumption and transport, considering (1) diffusion and enzymatic conversion of glucose to form hydrogen peroxide (H₂O₂) by the glucose-oxidase-based hydrogel functionalization of our biosensors at the microscale; (2) H₂O₂ oxidation at the electrode surface, leading to amperometric H₂O₂ readout; (3) glucose diffusion and glucose consumption by cancer cells in a spherical tissue model at the microscale; (4) glucose and H₂O₂ transport in our hanging-drop compartments at the macroscale; and (5) solvent evaporation, leading to glucose and H₂O₂ upconcentration. Our model relates the measured currents to the glucose concentrations generating the currents. The low limit of detection of our biosensors (0.4 ± 0.1 μM), combined with our current-fitting method, enabled us to reveal glucose dynamics within our system. By measuring glucose dynamics in hanging-drop compartments populated by cancer spheroids of various sizes, we could infer glucose distributions within the spheroid, which will help translate in vitro 3D tissue model results to in vivo.

Improving cancer treatments via dynamical biophysical models

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

Combined Influence of Nutrient Supply Level and Tissue Mechanical Properties on Benign Tumor Growth as Revealed by Mathematical Modeling

A continuous mathematical model of non-invasive avascular tumor growth in tissue is presented. The model considers tissue as a biphasic material, comprised of a solid matrix and interstitial fluid. The convective motion of tissue elements happens due to the gradients of stress, which change as a result of tumor cells proliferation and death. The model accounts for glucose as the crucial nutrient, supplied from the normal tissue, and can reproduce both diffusion-limited and stress-limited tumor growth. Approximate tumor growth curves are obtained semi-analytically in the limit of infinite tissue hydraulic conductivity, which implies instantaneous equalization of arising stress gradients. These growth curves correspond well to the numerical solutions and represent classical sigmoidal curves with a short initial exponential phase, subsequent almost linear growth phase and a phase with growth deceleration, in which tumor tends to reach its maximum volume. The influence of two model parameters on tumor growth curves is investigated: tissue hydraulic conductivity, which links the values of stress gradient and convective velocity of tissue phases, and tumor nutrient supply level, which corresponds to different permeability and surface area density of capillaries in the normal tissue that surrounds the tumor. In particular, it is demonstrated, that sufficiently low tissue hydraulic conductivity (intrinsic, e.g., to tumors arising from connective tissue) and sufficiently high nutrient supply can lead to formation of giant benign tumors, reaching tens of centimeters in diameter, which are indeed observed clinically.

A New Mathematical Model for Controlling Tumor Growth Based on Microenvironment Acidity and Oxygen Concentration

Hypoxia and the pH level of the tumor microenvironment have a great impact on the treatment of tumors. Here, the tumor growth is controlled by regulating the oxygen concentration and the acidity of the tumor microenvironment by introducing a two-dimensional multiscale cellular automata model of avascular tumor growth. The spatiotemporal evolution of tumor growth and metabolic variations is modeled based on biological assumptions, physical structure, states of cells, and transition rules. Each cell is allocated to one of the following states: proliferating cancer, nonproliferating cancer, necrotic, and normal cells. According to the response of the microenvironmental conditions, each cell consumes/produces metabolic factors and updates its state based on some stochastic rules. The input parameters are compatible with cancer biology using experimental data. The effect of neighborhoods during mitosis and simulating spatial heterogeneity is studied by considering multicellular layer structure of tumor. A simple Darwinist mutation is considered by introducing a critical parameter (Nmm) that affects division probability of the proliferative tumor cells based on the microenvironmental conditions and cancer hallmarks. The results show that Nmm regulation has a significant influence on the dynamics of tumor growth, the growth fraction, necrotic fraction, and the concentration levels of the metabolic factors. The model not only is able to simulate the in vivo tumor growth quantitatively and qualitatively but also can simulate the concentration of metabolic factors, oxygen, and acidity graphically. The results show the spatial heterogeneity effects on the proliferation of cancer cells and the rest of the system. By increasing Nmm, tumor shrinkage and significant increasing in the oxygen concentration and the pH value of the tumor microenvironment are observed. The results demonstrate the model's ability, providing an essential tool for simulating different tumor evolution scenarios of a patient and reliable prediction of spatiotemporal progression of tumors for utilizing in personalized therapy.

Angiogenesis and vessel co-option in a mathematical model of diffusive tumor growth: The role of chemotaxis

This work considers the propagation of a tumor from the stage of a small avascular sphere in a host tissue and the progressive onset of a tumor neovasculature stimulated by a pro-angiogenic factor secreted by hypoxic cells. The way new vessels are formed involves cell sprouting from pre-existing vessels and following a trail via a chemotactic mechanism (CM). Namely, it is first proposed a detailed general family of models of the CM, based on a statistical mechanics approach. The key hypothesis is that the CM is composed by two components: i) the well-known bias induced by the angiogenic factor gradient; ii) the presence of stochastic changes of the velocity direction, thus giving rise to a diffusive component. Then, some further assumptions and simplifications are applied in order to derive a specific model to be used in the simulations. The tumor progression is favored by its acidic aggression towards the healthy cells. The model includes the evolution of many biological and chemical species. Numerical simulations show the onset of a traveling wave eventually replacing the host tissue with a fully vascularized tumor. The results of simulations agree with experimental measures of the vasculature density in tumors, even in the case of particularly hypoxic tumors.

Mutant Evolution in Spatially Structured and Fragmented Expanding Populations

Mutant evolution in spatially structured systems is important for a range of biological systems, but aspects of it still require further elucidation. Adding to previous work, we provide a simple derivation of growth laws that characterize the number of mutants of different relative fitness in expanding populations in spatial models of different dimensionalities. These laws are universal and independent of "microscopic" modeling details. We further study the accumulation of mutants and find that, with advantageous and neutral mutants, more of them are present in spatially structured, compared to well-mixed colonies of the same size. The behavior of disadvantageous mutants is subtle: if they are disadvantageous through a reduction in division rates, the result is the same, and it is the opposite if the disadvantage is due to a death rate increase. Finally, we show that in all cases, the same results are observed in fragmented, nonspatial patch models. This suggests that the patterns observed are the consequence of population fragmentation, and not spatial restrictions per se We provide an intuitive explanation for the complex dependence of disadvantageous mutant evolution on spatial restriction, which relies on desynchronized dynamics in different locations/patches, and plays out differently depending on whether the disadvantage is due to a lower division rate or a higher death rate. Implications for specific biological systems, such as the evolution of drug-resistant cell mutants in cancer or bacterial biofilms, are discussed.

Optimization of Dose Fractionation for Radiotherapy of a Solid Tumor with Account of Oxygen Effect and Proliferative Heterogeneity

A spatially-distributed continuous mathematical model of solid tumor growth and treatment by fractionated radiotherapy is presented. The model explicitly accounts for three time and space-dependent factors that influence the efficiency of radiotherapy fractionation schemes—tumor cell repopulation, reoxygenation and redistribution of proliferative states. A special algorithm is developed, aimed at finding the fractionation schemes that provide increased tumor cure probability under the constraints of maximum normal tissue damage and maximum fractional dose. The optimization procedure is performed for varied radiosensitivity of tumor cells under the values of model parameters, corresponding to different degrees of tumor malignancy. The resulting optimized schemes consist of two stages. The first stages are aimed to increase the radiosensitivity of the tumor cells, remaining after their end, sparing the caused normal tissue damage. This allows to increase the doses during the second stages and thus take advantage of the obtained increased radiosensitivity. Such method leads to significant expansions in the curative ranges of the values of tumor radiosensitivity parameters. Overall, the results of this study represent the theoretical proof of concept that non-uniform radiotherapy fractionation schemes may be considerably more effective that uniform ones, due to the time and space-dependent effects.

Mathematical Modeling Shows That the Response of a Solid Tumor to Antiangiogenic Therapy Depends on the Type of Growth

It has been hypothesized that solid tumors with invasive type of growth should possess intrinsic resistance to antiangiogenic therapy, which is aimed at cessation of the formation of new blood vessels and subsequent shortage of nutrient inflow to the tumor. In order to investigate this effect, a continuous mathematical model of tumor growth is developed, which considers variables of tumor cells, necrotic tissue, capillaries, and glucose as the crucial nutrient. The model accounts for the intrinsic motility of tumor cells and for the convective motion, arising due to their proliferation, thus allowing considering two types of tumor growth—invasive and compact—as well as their combination. Analytical estimations of tumor growth speed are obtained for compact and invasive tumors. They suggest that antiangiogenic therapy may provide a several times decrease of compact tumor growth speed, but the decrease of growth speed for invasive tumors should be only modest. These estimations are confirmed by numerical simulations, which further allow evaluating the effect of antiangiogenic therapy on tumors with mixed growth type and highlight the non-additive character of the two types of growth.

Acid-mediated tumor invasion as a function of nutrient source location

Cancer cells have an altered metabolism that increases acid production driving to an extracellular pH significantly lower than normal. This leads to normal cell death, and extracellular matrix degradation allowing the formation of an interstitial gap between cancer and healthy cells. In this work, we present a mathematical model to study the interstitial gap formation and evolution considering a tissue with a non-uniform nutrient distribution. Our results indicate that the interstitial gap onsets at the region with highest nutrient consumption. Due to the gap formation, cancer cells near the interface have more nutrient and space availability. This induces cancer cell reproduction and migration toward the nutrient source. Our simulations suggest a strong correlation between gap size and the distance to the nutrient source. Although we do not find a correlation between tumor growth speed and gap size, our results indicate a high risk of metastasis for tumors that develop an interstitial gap, emphasizing the importance of gap detection as a hallmark for cancer invasion.

Investigation of solid tumor progression with account of proliferation/migration dichotomy via Darwinian mathematical model

A new continuous spatially-distributed model of solid tumor growth and progression is presented. The model explicitly accounts for mutations/epimutations of tumor cells which take place upon their division. The tumor grows in normal tissue and its progression is driven only by competition between populations of malignant cells for limited nutrient supply. Two reasons for the motion of tumor cells in space are taken into consideration, i.e., their intrinsic motility and convective fluxes, which arise due to proliferation of tumor cells. The model is applied to investigation of solid tumor progression under phenotypic alterations that inversely affect cell proliferation rate and cell motility by increasing the value of one of the parameters at the expense of another.It is demonstrated that the crucial feature that gives evolutionary advantage to a cell population is the speed of its intergrowth into surrounding normal tissue. Of note, increase in tumor intergrowth speed in not always associated with increase in motility of tumor cells. Depending on the parameters of functions, that describe phenotypic alterations, tumor cellular composition may evolve towards: (1) maximization of cell proliferation rate, (2) maximization of cell motility, (3) non-extremum values of cell proliferation rate and motility. Scenarios are found, where after initial tendency for maximization of cell proliferation rate, the direction of tumor progression sharply switches to maximization of cell motility, which is accompanied by decrease in total speed of tumor growth.

Three-Dimensional Modeling of Avascular Tumor Growth in Both Static and Dynamic Culture Platforms

Microfluidic cell culture platforms are ideal candidates for modeling the native tumor microenvironment because they can precisely reconstruct in vivo cellular behavior. Moreover, mathematical modeling of tumor growth can pave the way toward description and prediction of growth pattern as well as improving cancer treatment. In this study, a modified mathematical model based on concentration distribution is applied to tumor growth in both conventional static culture and dynamic microfluidic cell culture systems. Apoptosis and necrosis mechanisms are considered as the main inhibitory factors in the model, while tumor growth rate and nutrient consumption rate are modified in both quiescent and proliferative zones. We show that such modification can better predict the experimental results of tumor growth reported in the literature. Using numerical simulations, the effects of the concentrations of the nutrients as well as the initial tumor radius on the tumor growth are investigated and discussed. Furthermore, tumor growth is simulated by taking into account the dynamic perfusion into the proposed model. Subsequently, tumor growth kinetics in a three-dimensional (3D) microfluidic device containing a U-shaped barrier is numerically studied. For this case, the effect of the flow rate of culture medium on tumor growth is investigated as well. Finally, to evaluate the impact of the trap geometry on the tumor growth, a comparison is made between the tumor growth kinetics in two frequently used traps in microfluidic cell culture systems, i.e., the U-shaped barrier and microwell structures. The proposed model can provide insight into better predicting the growth and development of avascular tumor in both static and dynamic cell culture platforms.

Three-dimensional tumor growth in time-varying chemical fields: a modeling framework and theoretical study

Background

Contemporary biological observations have revealed a large variety of mechanisms acting during the expansion of a tumor. However, there are still many qualitative and quantitative aspects of the phenomenon that remain largely unknown. In this context, mathematical and computational modeling appears as an invaluable tool providing the means for conducting in silico experiments, which are cheaper and less tedious than real laboratory experiments.

Results

This paper aims at developing an extensible and computationally efficient framework for in silico modeling of tumor growth in a 3-dimensional, inhomogeneous and time-varying chemical environment. The resulting model consists of a set of mathematically derived and algorithmically defined operators, each one addressing the effects of a particular biological mechanism on the state of the system. These operators may be extended or re-adjusted, in case a different set of starting assumptions or a different simulation scenario needs to be considered.

Conclusion

In silico modeling provides an alternative means for testing hypotheses and simulating scenarios for which exact biological knowledge remains elusive. However, finer tuning of pertinent methods presupposes qualitative and quantitative enrichment of available biological evidence. Validation in a strict sense would further require comprehensive, case-specific simulations and detailed comparisons with biomedical observations.

Mathematical modeling of the proliferation gradient in multicellular tumor spheroids

MultiCellular Tumor Spheroids are 3D cell cultures that can accurately reproduce the behavior of solid tumors. It has been experimentally observed that large spheroids exhibit a decreasing gradient of proliferation from the periphery to the center of these multicellular 3D models: the proportion of proliferating cells is higher in the periphery while the non-proliferating quiescent cells increase in depth. In this paper, we propose to investigate the key mechanisms involved in the establishment of this gradient with a Partial Differential Equations model that mimics the experimental set-up of growing spheroids under different nutrients supply conditions. The model consists of mass balance equations on the two cell populations observed in the data: the proliferating cells and the quiescent cells. The spherical symmetry is used to rewrite the model in radial and relative coordinates. Thanks to a rigorous data postprocessing the model is then fit and compared quantitatively with the experimental quantification of the percentage of proliferating cells from EdU immunodetection on 2D spheroid cryosection images. The results of this calibration show that the proliferation gradient observed in spheroids can be quantitatively reproduced by our model.

Transient alleviation of tumor hypoxia during first days of antiangiogenic therapy as a result of therapy-induced alterations in nutrient supply and tumor metabolism - Analysis by mathematical modeling

A number of experiments on mouse tumor models, as well as certain clinical data, have demonstrated, that antiangiogenic therapy can lead to transient improvement in tumor oxygenation, that allows to increase efficiency of following radiotherapy. In the majority of works, this phenomenon has been explained by enhanced tumor perfusion due to normalization of capillaries' structure, that results in elevated oxygen inflow in tumor. However, changes in tumor perfusion often haven't been directly measured in relevant works, moreover, antiangiogenic therapy has been proven to have ambiguous effect on tumor perfusion both in mouse tumor models and in clinics. Herein, we suggest that elevation of blood perfusion may be not the only reason for transient alleviation of tumor hypoxia, and that it may manifest itself even under unchanged tumor blood flow. We propose that it may be as well caused by the decrease in tumor oxygen consumption rate (OCR) due to the reduction of tumor proliferation level, caused by nutrient shortage in result of antiangiogenic treatment. We provide detailed explanation of this hypothesis and visualize it using a specially developed mathematical model, which takes into account basic features of tumor growth and antiangiogenic therapy. We investigate the influence of the model parameters on oxygen dynamics; demonstrate, that transient alleviation of tumor hypoxia occurs in a fairly wide range of physiologically justified values of parameters; and point out the major factors, that determine oxygen dynamics during antiangiogenic therapy.

The Exponential-Gompertzian Tumor Growth Model: Data from Six Tumor Cell Lines in Vitro and in Vivo. Estimate of the Transition point from Exponential to Gompertzian Growth and Potential Cinical Implications

The published growth data were examined for six tumor cell lines (FSA, Line 1, MCA-11, EMT6/RO, MGH-U1, MLS) grown in vivo and in vitro as monolayer cultures and as multicell spheroids cultured under different experimental conditions. Serial estimates of tumor sizes were fitted by Gompertzian equations obtained with a non-linear computerized program. When the growth equations of the same tumor growing in different experimental conditions were compared, the Gompertzian parameters α0 (initial specific growth rate) and β (retardation factor) showed a strong linear correlation in all the examined lines, with no exception. This occurrence supports the exponential-Gompertzian growth model, where an early exponential phase (which is virtually not influenced by exogenous factors) is followed by a Gompertzian phase, the characteristics of which are greatly dependent on environmental conditions. The transition between the two phases was estimated to occur when tumor size reached 102–104 cells, depending on the cell line. This kinetic change in tumor growth may be clinically relevant as regards cytotoxic treatments. It could explain some consequences of delays in adjuvant (postoperative) chemotherapy observed in clinical trials on primary breast cancer.

A multiscale model for heterogeneous tumor spheroid in vitro

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

Microfluidic device flow field characterization around tumor spheroids with tunable necrosis produced in an optimized off-chip process

Tumor spheroids are a 3-D tumor model that holds promise for testing cancer therapies in vitro using microfluidic devices. Tailoring the properties of a tumor spheroid is critical for evaluating therapies over a broad range of possible indications. Using human colon cancer cells (HCT-116), we demonstrate controlled tumor spheroid growth rates by varying the number of cells initially seeded into microwell chambers. The presence of a necrotic core in the spheroids could be controlled by changing the glucose concentration of the incubation medium. This manipulation had no effect on the size of the tumor spheroids or hypoxia in the spheroid core, which has been predicted by a mathematical model in computer simulations of spheroid growth. Control over the presence of a necrotic core while maintaining other physical parameters of the spheroid presents an opportunity to assess the impact of core necrosis on therapy efficacy. Using micro-particle imaging velocimetry (micro-PIV), we characterize the hydrodynamics and mass transport of nanoparticles in tumor spheroids in a microfluidic device. We observe a geometrical dependence on the flow rate experienced by the tumor spheroid in the device, such that the "front" of the spheroid experiences a higher flow velocity than the "back" of the spheroid. Using fluorescent nanoparticles, we demonstrate a heterogeneous accumulation of nanoparticles at the tumor interface that correlates with the observed flow velocities. The penetration depth of these nanoparticles into the tumor spheroid depends on nanoparticle diameter, consistent with reports in the literature.

Fractal dimension and universality in avascular tumor growth

For years, the comprehension of the tumor growth process has been intriguing scientists. New research has been constantly required to better understand the complexity of this phenomenon. In this paper, we propose a mathematical model that describes the properties, already known empirically, of avascular tumor growth. We present, from an individual-level (microscopic) framework, an explanation of some phenomenological (macroscopic) aspects of tumors, such as their spatial form and the way they develop. Our approach is based on competitive interaction between the cells. This simple rule makes the model able to reproduce evidence observed in real tumors, such as exponential growth in their early stage followed by power-law growth. The model also reproduces (i) the fractal-space distribution of tumor cells and (ii) the universal growth behavior observed in both animals and tumors. Our analyses suggest that the universal similarity between tumor and animal growth comes from the fact that both can be described by the same dynamic equation-the Bertalanffy-Richards model-even if they do not necessarily share the same biological properties.

Simulation-assisted design of microfluidic sample traps for optimal trapping and culture of non-adherent single cells, tissues, and spheroids

This work focuses on modelling design and operation of "microfluidic sample traps" (MSTs). MSTs regroup a widely used class of microdevices that incorporate wells, recesses or chambers adjacent to a channel to individually trap, culture and/or release submicroliter 3D tissue samples ranging from simple cell aggregates and spheroids, to ex vivo tissue samples and other submillimetre-scale tissue models. Numerous MST designs employing various trapping mechanisms have been proposed in the literature, spurring the development of 3D tissue models for drug discovery and personalized medicine. Yet, there lacks a general framework to optimize trapping stability, trapping time, shear stress, and sample metabolism. Herein, the effects of hydrodynamics and diffusion-reaction on tissue viability and device operation are investigated using analytical and finite element methods with systematic parametric sweeps over independent design variables chosen to correspond to the four design degrees of freedom. Combining different results, we show that, for a spherical tissue of diameter d < 500 μm, the simplest, closest to optimal trap shape is a cube of dimensions w equal to twice the tissue diameter: w = 2d. Furthermore, to sustain tissues without perfusion, available medium volume per trap needs to be 100× the tissue volume to ensure optimal metabolism for at least 24 hours.

The role of necrosis, acute hypoxia and chronic hypoxia in 18F-FMISO PET image contrast: a computational modelling study

Positron emission tomography (PET) using 18F-fluoromisonidazole (FMISO) is a promising technique for imaging tumour hypoxia, and a potential target for radiotherapy dose-painting. However, the relationship between FMISO uptake and oxygen partial pressure ([Formula: see text]) is yet to be quantified fully. Tissue oxygenation varies over distances much smaller than clinical PET resolution (<100 μm versus  ∼4 mm), and cyclic variations in tumour perfusion have been observed on timescales shorter than typical FMISO PET studies (∼20 min versus a few hours). Furthermore, tracer uptake may be decreased in voxels containing some degree of necrosis. This work develops a computational model of FMISO uptake in millimetre-scale tumour regions. Coupled partial differential equations govern the evolution of oxygen and FMISO distributions, and a dynamic vascular source map represents temporal variations in perfusion. Local FMISO binding capacity is modulated by the necrotic fraction. Outputs include spatiotemporal maps of [Formula: see text] and tracer accumulation, enabling calculation of tissue-to-blood ratios (TBRs) and time-activity curves (TACs) as a function of mean tissue oxygenation. The model is characterised using experimental data, finding half-maximal FMISO binding at local [Formula: see text] of 1.4 mmHg (95% CI: 0.3-2.6 mmHg) and half-maximal necrosis at 1.2 mmHg (0.1-4.9 mmHg). Simulations predict a non-linear non-monotonic relationship between FMISO activity (4 hr post-injection) and mean tissue [Formula: see text] : tracer uptake rises sharply from negligible levels in avascular tissue, peaking at  ∼5 mmHg and declining towards blood activity in well-oxygenated conditions. Greater temporal variation in perfusion increases peak TBRs (range 2.20-5.27) as a result of smaller predicted necrotic fraction, rather than fundamental differences in FMISO accumulation under acute hypoxia. Identical late FMISO uptake can occur in regions with differing [Formula: see text] and necrotic fraction, but simulated TACs indicate that additional early-phase information may allow discrimination of hypoxic and necrotic signals. We conclude that a robust approach to FMISO interpretation (and dose-painting prescription) is likely to be based on dynamic PET analysis.

Limiting the development of anti-cancer drug resistance in a spatial model of micrometastases

While chemoresistance in primary tumors is well-studied, much less is known about the influence of systemic chemotherapy on the development of drug resistance at metastatic sites. In this work, we use a hybrid spatial model of tumor response to a DNA damaging drug to study how the development of chemoresistance in micrometastases depends on the drug dosing schedule. We separately consider cell populations that harbor pre-existing resistance to the drug, and those that acquire resistance during the course of treatment. For each of these independent scenarios, we consider one hypothetical cell line that is responsive to metronomic chemotherapy, and another that with high probability cannot be eradicated by a metronomic protocol. Motivated by experimental work on ovarian cancer xenografts, we consider all possible combinations of a one week treatment protocol, repeated for three weeks, and constrained by the total weekly drug dose. Simulations reveal a small number of fractionated-dose protocols that are at least as effective as metronomic therapy in eradicating micrometastases with acquired resistance (weak or strong), while also being at least as effective on those that harbor weakly pre-existing resistant cells. Given the responsiveness of very different theoretical cell lines to these few fractionated-dose protocols, these may represent more effective ways to schedule chemotherapy with the goal of limiting metastatic tumor progression.

Modelling tumour cell proliferation from vascular structure using tissue decomposition into avascular elements

Computer models allow the mechanistically detailed study of tumour proliferation and its dependency on nutrients. However, the computational study of large vascular tumours requires detailed information on the 3-dimensional vessel network and rather high computation times due to complex geometries. This study puts forward the idea of partitioning vascularised tissue into connected avascular elements that can exchange cells and nutrients between each other. Our method is able to rapidly calculate the evolution of proliferating as well as dead and quiescent cells, and hence a proliferative index, from a given amount and distribution of vascularisation of arbitrary complexity. Applying our model, we found that a heterogeneous vessel distribution provoked a higher proliferative index, suggesting increased malignancy, and increased the amount of dead cells compared to a more static tumour environment when a homogenous vessel distribution was assumed. We subsequently demonstrated that under certain amounts of vascularisation, cell proliferation may even increase when vessel density decreases, followed by a subsequent decrease of proliferation. This effect was due to a trade-off between an increase in compensatory proliferation for replacing dead cells and a decrease of cell population due to lack of oxygen supply in lowly vascularised tumours. Findings were illustrated by an ectopic colorectal cancer mouse xenograft model. Our presented approach can be in the future applied to study the effect of cytostatic, cytotoxic and anti-angiogenic chemotherapy and is ideally suited for translational systems biology, where rapid interaction between theory and experiment is essential.

Inferring Growth Control Mechanisms in Growing Multi-cellular Spheroids of NSCLC Cells from Spatial-Temporal Image Data

We develop a quantitative single cell-based mathematical model for multi-cellular tumor spheroids (MCTS) of SK-MES-1 cells, a non-small cell lung cancer (NSCLC) cell line, growing under various nutrient conditions: we confront the simulations performed with this model with data on the growth kinetics and spatial labeling patterns for cell proliferation, extracellular matrix (ECM), cell distribution and cell death. We start with a simple model capturing part of the experimental observations. We then show, by performing a sensitivity analysis at each development stage of the model that its complexity needs to be stepwise increased to account for further experimental growth conditions. We thus ultimately arrive at a model that mimics the MCTS growth under multiple conditions to a great extent. Interestingly, the final model, is a minimal model capable of explaining all data simultaneously in the sense, that the number of mechanisms it contains is sufficient to explain the data and missing out any of its mechanisms did not permit fit between all data and the model within physiological parameter ranges. Nevertheless, compared to earlier models it is quite complex i.e., it includes a wide range of mechanisms discussed in biological literature. In this model, the cells lacking oxygen switch from aerobe to anaerobe glycolysis and produce lactate. Too high concentrations of lactate or too low concentrations of ATP promote cell death. Only if the extracellular matrix density overcomes a certain threshold, cells are able to enter the cell cycle. Dying cells produce a diffusive growth inhibitor. Missing out the spatial information would not permit to infer the mechanisms at work. Our findings suggest that this iterative data integration together with intermediate model sensitivity analysis at each model development stage, provide a promising strategy to infer predictive yet minimal (in the above sense) quantitative models of tumor growth, as prospectively of other tissue organization processes. Importantly, calibrating the model with two nutriment-rich growth conditions, the outcome for two nutriment-poor growth conditions could be predicted. As the final model is however quite complex, incorporating many mechanisms, space, time, and stochastic processes, parameter identification is a challenge. This calls for more efficient strategies of imaging and image analysis, as well as of parameter identification in stochastic agent-based simulations.

We here present how to parameterize a mathematical agent-based model of growing MCTS almost completely from experimental data. MCTS show a similar establishment of pathophysiological gradients and concentric arrangement of heterogeneous cell populations as found in avascular tumor nodules. We build a process chain of imaging, image processing and analysis, and mathematical modeling. In this model, each individual cell is represented by an agent populating one site of a three dimensional un-structured lattice. The spatio-temporal multi-cellular behavior, including migration, growth, division, death of each cell, is considered by a stochastic process, simulated numerically by the Gillespie algorithm. Processes on the molecular scale are described by deterministic partial differential equations for molecular concentrations, coupled to intracellular and cellular decision processes. The parameters of the multi-scale model are inferred from comparisons to the growth kinetics and from image analysis of spheroid cryosections stained for cell death, proliferation and collagen IV. Our final model assumes ATP to be the critical resource that cells try to keep constant over a wide range of oxygen and glucose medium concentrations, by switching between aerobic and anaerobic metabolism. Besides ATP, lactate is shown to be a possible explanation for the control of the necrotic core size. Direct confrontation of the model simulation results with image data on the spatial profiles of cell proliferation, ECM distribution and cell death, indicates that in addition, the effects of ECM and waste factors have to be added to explain the data. Hence the model is a tool to identify likely mechanisms at work that may subsequently be studied experimentally, proposing a model-guided experimental strategy.

Micro-dissected tumor tissues on chip: an ex vivo method for drug testing and personalized therapy

In cancer research and personalized medicine, new tissue culture models are needed to better predict the response of patients to therapies. With a concern for the small volume of tissue typically obtained through a biopsy, we describe a method to reproducibly section live tumor tissue to submillimeter sizes. These micro-dissected tissues (MDTs) share with spheroids the advantages of being easily manipulated on-chip and kept alive for periods extending over one week, while being biologically relevant for numerous assays. At dimensions below ~420 μm in diameter, as suggested by a simple metabolite transport model and confirmed experimentally, continuous perfusion is not required to keep samples alive, considerably simplifying the technical challenges. For the long-term culture of MDTs, we describe a simple microfluidic platform that can reliably trap samples in a low shear stress environment. We report the analysis of MDT viability for eight different types of tissues (four mouse xenografts derived from human cancer cell lines, three from ovarian and prostate cancer patients, and one from a patient with benign prostatic hyperplasia) analyzed by both confocal microscopy and flow cytometry over an 8-day incubation period. Finally, we provide a proof of principle for chemosensitivity testing of human tissue from a cancer patient performed using the described MDT chip method. This technology has the potential to improve treatment success rates by identifying potential responders earlier during the course of treatment and providing opportunities for direct drug testing on patient tissues in early drug development stages.

A Computational Model for Investigating Tumor Apoptosis Induced by Mesenchymal Stem Cell-Derived Secretome

Apoptosis is a programmed cell death that occurs naturally in physiological and pathological conditions. Defective apoptosis can trigger the development and progression of cancer. Experiments suggest the ability of secretome derived from mesenchymal stem cells (MSC) to induce apoptosis in cancer cells. We develop a hybrid discrete-continuous multiscale model to further investigate the effect of MSC-derived secretome in tumor growth. The model encompasses three biological scales. At the molecular scale, a system of ordinary differential equations regulate the expression of proteins involved in apoptosis signaling pathways. At the cellular scale, discrete equations control cellular migration, phenotypic switching, and proliferation. At the extracellular scale, a system of partial differential equations are employed to describe the dynamics of microenvironmental chemicals concentrations. The simulation is able to produce both avascular tumor growth rate and phenotypic patterns as observed in the experiments. In addition, we obtain good quantitative agreements with the experimental data on the apoptosis of HeLa cancer cells treated with MSC-derived secretome. We use this model to predict the growth of avascular tumor under various secretome concentrations over time.

A Glycolysis-Based In Silico Model for the Solid Tumor Growth

Cancer-tumor growth is a complex process depending on several biological factors, such as the chemical microenvironment of the tumor, the cellular metabolic profile, and its proliferation rate. Several mathematical models have been developed for identifying the interactions between tumor cells and tissue microenvironment, since they play an important role in tumor formation and progression. Toward this direction we propose a new continuum model of avascular glioma-tumor growth, which incorporates a new factor, namely, the glycolytic potential of cancer cells, to express the interactions of three different tumor-cell populations (proliferative, hypoxic, and necrotic) with their tissue microenvironment. The glycolytic potential engages three vital nutrients, i.e., oxygen, glucose, and lactate, which provide cells with the necessary energy for their survival and proliferation. Extensive simulations are performed for different evolution times and various proliferation rates, in order to investigate how the tumor growth is affected. According to medical experts, the experimental observations indicate that the model predicts quite satisfactorily the overall tumor growth as well as the expansion of each region separately. Following extensive evaluation, the proposed model may provide an essential tool for patient-specific tumor simulation and reliable prediction of glioma spatiotemporal expansion.

Spatio-Temporal Dynamics of Hypoxia during Radiotherapy

Tumour hypoxia plays a pivotal role in cancer therapy for most therapeutic approaches from radiotherapy to immunotherapy. The detailed and accurate knowledge of the oxygen distribution in a tumour is necessary in order to determine the right treatment strategy. Still, due to the limited spatial and temporal resolution of imaging methods as well as lacking fundamental understanding of internal oxygenation dynamics in tumours, the precise oxygen distribution map is rarely available for treatment planing. We employ an agent-based in silico tumour spheroid model in order to study the complex, localized and fast oxygen dynamics in tumour micro-regions which are induced by radiotherapy. A lattice-free, 3D, agent-based approach for cell representation is coupled with a high-resolution diffusion solver that includes a tissue density-dependent diffusion coefficient. This allows us to assess the space- and time-resolved reoxygenation response of a small subvolume of tumour tissue in response to radiotherapy. In response to irradiation the tumour nodule exhibits characteristic reoxygenation and re-depletion dynamics which we resolve with high spatio-temporal resolution. The reoxygenation follows specific timings, which should be respected in treatment in order to maximise the use of the oxygen enhancement effects. Oxygen dynamics within the tumour create windows of opportunity for the use of adjuvant chemotherapeutica and hypoxia-activated drugs. Overall, we show that by using modelling it is possible to follow the oxygenation dynamics beyond common resolution limits and predict beneficial strategies for therapy and in vitro verification. Models of cell cycle and oxygen dynamics in tumours should in the future be combined with imaging techniques, to allow for a systematic experimental study of possible improved schedules and to ultimately extend the reach of oxygenation monitoring available in clinical treatment.

The effect of intrinsic and acquired resistances on chemotherapy effectiveness

Although chemotherapy is one of the most common treatments for cancer, it can be only partially successful. Drug resistance is the main cause of the failure of chemotherapy. In this work, we present a mathematical model to study the impact of both intrinsic (preexisting) and acquired (induced by the drugs) resistances on chemotherapy effectiveness. Our simulations show that intrinsic resistance could be as dangerous as acquired resistance. In particular, our simulations suggest that tumors composed by even a small fraction of intrinsically resistant cells may lead to an unsuccessful therapy very quickly. Our results emphasize the importance of monitoring both intrinsic and acquired resistances during treatment in order to succeed and the importance of doing more experimental and genetic research in order to develop a pretreatment clinical test to avoid intrinsic resistance.

A non-phenomenological model of competition and cooperation to explain population growth behaviors

This paper is an extension of a previous work which proposes a non-phenomenological model of population growth that is based on the interactions among the individuals of a population. In addition to what had already been studied—that the individuals interact competitively—in the present work it is also considered that the individuals interact cooperatively. As a consequence of this new consideration, a richer dynamics is observed. For instance, besides getting the population models already reached from the original version of the model (as the Malthus, Verhulst, Gompertz, Richards, Bertalanffy and power-law growth models), the new formulation also reaches the von Foerster growth model and also a regime of divergence of the population at a finite time. An agent-based model is also presented in order to give support to the analytical results. Moreover, this new approach of the model explains the Allee effect as an emergent behavior of the cooperative and competitive interactions among the individuals. The Allee effect is the characteristic of some populations of increasing the population growth rate in a small-sized population. Whereas the models presented in the literature explain the Allee effect with phenomenological ideas, the model presented here explains this effect by the interactions between the individuals. The model is tested with empirical data to justify its formulation. Another interesting macroscopic emergent behavior from the model proposed is the observation of a regime of population divergence at a finite time. It is interesting that this characteristic is observed in humanity's global population growth. It is shown that in a regime of cooperation, the model fits very well to the human population growth data from 1000 AD to nowadays.

A matter of timing: identifying significant multi-dose radiotherapy improvements by numerical simulation and genetic algorithm search

Multi-dose radiotherapy protocols (fraction dose and timing) currently used in the clinic are the product of human selection based on habit, received wisdom, physician experience and intra-day patient timetabling. However, due to combinatorial considerations, the potential treatment protocol space for a given total dose or treatment length is enormous, even for relatively coarse search; well beyond the capacity of traditional in-vitro methods. In constrast, high fidelity numerical simulation of tumor development is well suited to the challenge. Building on our previous single-dose numerical simulation model of EMT6/Ro spheroids, a multi-dose irradiation response module is added and calibrated to the effective dose arising from 18 independent multi-dose treatment programs available in the experimental literature. With the developed model a constrained, non-linear, search for better performing cadidate protocols is conducted within the vicinity of two benchmarks by genetic algorithm (GA) techniques. After evaluating less than 0.01% of the potential benchmark protocol space, candidate protocols were identified by the GA which conferred an average of 9.4% (max benefit 16.5%) and 7.1% (13.3%) improvement (reduction) on tumour cell count compared to the two benchmarks, respectively. Noticing that a convergent phenomenon of the top performing protocols was their temporal synchronicity, a further series of numerical experiments was conducted with periodic time-gap protocols (10 h to 23 h), leading to the discovery that the performance of the GA search candidates could be replicated by 17-18 h periodic candidates. Further dynamic irradiation-response cell-phase analysis revealed that such periodicity cohered with latent EMT6/Ro cell-phase temporal patterning. Taken together, this study provides powerful evidence towards the hypothesis that even simple inter-fraction timing variations for a given fractional dose program may present a facile, and highly cost-effecitive means of significantly improving clinical efficacy.

Influence of individual cell motility on the 2D front roughness dynamics of tumour cell colonies

The dynamics of in situ 2D HeLa cell quasi-linear and quasi-radial colony fronts in a standard culture medium is investigated. For quasi-radial colonies, as the cell population increased, a kinetic transition from an exponential to a constant front average velocity regime was observed. Special attention was paid to individual cell motility evolution under constant average colony front velocity looking for its impact on the dynamics of the 2D colony front roughness. From the directionalities and velocity components of cell trajectories in colonies with different cell populations, the influence of both local cell density and cell crowding effects on individual cell motility was determined. The average dynamic behaviour of individual cells in the colony and its dependence on both local spatio-temporal heterogeneities and growth geometry suggested that cell motion undergoes under a concerted cell migration mechanism, in which both a limiting random walk-like and a limiting ballistic-like contribution were involved. These results were interesting to infer how biased cell trajectories influenced both the 2D colony spreading dynamics and the front roughness characteristics by local biased contributions to individual cell motion. These data are consistent with previous experimental and theoretical cell colony spreading data and provide additional evidence of the validity of the Kardar-Parisi-Zhang equation, within a certain range of time and colony front size, for describing the dynamics of 2D colony front roughness.

In silico analysis of cell cycle synchronisation effects in radiotherapy of tumour spheroids

Tumour cells show a varying susceptibility to radiation damage as a function of the current cell cycle phase. While this sensitivity is averaged out in an unperturbed tumour due to unsynchronised cell cycle progression, external stimuli such as radiation or drug doses can induce a resynchronisation of the cell cycle and consequently induce a collective development of radiosensitivity in tumours. Although this effect has been regularly described in experiments it is currently not exploited in clinical practice and thus a large potential for optimisation is missed. We present an agent-based model for three-dimensional tumour spheroid growth which has been combined with an irradiation damage and kinetics model. We predict the dynamic response of the overall tumour radiosensitivity to delivered radiation doses and describe corresponding time windows of increased or decreased radiation sensitivity. The degree of cell cycle resynchronisation in response to radiation delivery was identified as a main determinant of the transient periods of low and high radiosensitivity enhancement. A range of selected clinical fractionation schemes is examined and new triggered schedules are tested which aim to maximise the effect of the radiation-induced sensitivity enhancement. We find that the cell cycle resynchronisation can yield a strong increase in therapy effectiveness, if employed correctly. While the individual timing of sensitive periods will depend on the exact cell and radiation types, enhancement is a universal effect which is present in every tumour and accordingly should be the target of experimental investigation. Experimental observables which can be assessed non-invasively and with high spatio-temporal resolution have to be connected to the radiosensitivity enhancement in order to allow for a possible tumour-specific design of highly efficient treatment schedules based on induced cell cycle synchronisation.

The sensitivity of a cell to a dose of radiation is largely affected by its current position within the cell cycle. While under normal circumstances progression through the cell cycle will be asynchronous in a tumour mass, external influences such as chemo- or radiotherapy can induce a synchronisation. Such a common progression of the inner clock of the cancer cells results in the critical dependence on the effectiveness of any drug or radiation dose on a suitable timing for its administration. We analyse the exact evolution of the radiosensitivity of a sample tumour spheroid in a computer model, which enables us to predict time windows of decreased or increased radiosensitivity. Fractionated radiotherapy schedules can be tailored in order to avoid periods of high resistance and exploit the induced radiosensitivity for an increase in therapy efficiency. We show that the cell cycle effects can drastically alter the outcome of fractionated irradiation schedules in a spheroid cell system. By using the correct observables and continuous monitoring, the cell cycle sensitivity effects have the potential to be integrated into treatment planing of the future and thus to be employed for a better outcome in clinical cancer therapies.

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.

Tumor growth dynamics: insights into evolutionary processes

Identifying the types of event that drive tumor evolution and progression is crucial for understanding cancer. We suggest that the analysis of tumor growth dynamics can provide a window into tumor biology and evolution by connecting them with the types of genetic change that have occurred. Although fundamentally important, the documentation of tumor growth kinetics is more sparse in the literature than is the molecular analysis of cells. Here, we provide a historical summary of tumor growth patterns and argue that they can be classified into five basic categories. We then illustrate how those categories can provide insights into events that drive tumor progression, by discussing a particular evolutionary model as an example and encouraging such analysis in a more general setting.

Ovarian tumor attachment, invasion, and vascularization reflect unique microenvironments in the peritoneum: insights from xenograft and mathematical models

Ovarian cancer relapse is often characterized by metastatic spread throughout the peritoneal cavity with tumors attached to multiple organs. In this study, interaction of ovarian cancer cells with the peritoneal tumor microenvironment was evaluated in a xenograft model based on intraperitoneal injection of fluorescent SKOV3.ip1 ovarian cancer cells. Intra-vital microscopy of mixed GFP-red fluorescent protein (RFP) cell populations injected into the peritoneum demonstrated that cancer cells aggregate and attach as mixed spheroids, emphasizing the importance of homotypic adhesion in tumor formation. Electron microscopy provided high resolution structural information about local attachment sites. Experimental measurements from the mouse model were used to build a three-dimensional cellular Potts ovarian tumor model (OvTM) that examines ovarian cancer cell attachment, chemotaxis, growth, and vascularization. OvTM simulations provide insight into the relative influence of cancer cell-cell adhesion, oxygen availability, and local architecture on tumor growth and morphology. Notably, tumors on the mesentery, omentum, or spleen readily invade the "open" architecture, while tumors attached to the gut encounter barriers that restrict invasion and instead rapidly expand into the peritoneal space. Simulations suggest that rapid neovascularization of SKOV3.ip1 tumors is triggered by constitutive release of angiogenic factors in the absence of hypoxia. This research highlights the importance of cellular adhesion and tumor microenvironment in the seeding of secondary ovarian tumors on diverse organs within the peritoneal cavity. Results of the OvTM simulations indicate that invasion is strongly influenced by features underlying the mesothelial lining at different sites, but is also affected by local production of chemotactic factors. The integrated in vivo mouse model and computer simulations provide a unique platform for evaluating targeted therapies for ovarian cancer relapse.

A numerical model of EMT6/Ro spheroid dynamics under irradiation: calibration and estimation of the underlying irradiation-induced cell survival probability

We present extensions to our quasi-2D cellular automata spheroid model that add a cellular kinetics module together with an irradiation and repair module. Significantly, our approach is not based on the Linear Quadratic (LQ) model, instead, we propose a simple two-parameter, algorithmic model which captures the essential biological features of irradiation-induced cell death, repair and associated cell cycle delays. This approach allows us to estimate directly the underlying irradiation-induced cell survival probability. We present the calibration of this extended model both with and without the application of single irradiation doses to the commonly studied (in vitro) EMT6/Ro (mammary carcinoma) cell line. A comparison of the estimated underlying cell survival probability with the in vitro survival probability data confirms the expected differences in the measures.

Mechano-transduction in tumour growth modelling

The evolution of biological systems is strongly influenced by physical factors, such as applied forces, geometry or the stiffness of the micro-environment. Mechanical changes are particularly important in solid tumour development, as altered stromal-epithelial interactions can provoke a persistent increase in cytoskeletal tension, driving the gene expression of a malignant phenotype. In this work, we propose a novel multi-scale treatment of mechano-transduction in cancer growth. The avascular tumour is modelled as an expanding elastic spheroid, whilst growth may occur both as a volume increase and as a mass production within a cell rim. Considering the physical constraints of an outer healthy tissue, we derive the thermo-dynamical requirements for coupling growth rate, solid stress and diffusing biomolecules inside a heterogeneous tumour. The theoretical predictions successfully reproduce the stress-dependent growth curves observed by in vitro experiments on multicellular spheroids.

The buckling of capillaries in solid tumours

We develop a model of the buckling (both planar and axial) of capillaries in cancer tumours, using nonlinear solid mechanics. The compressive stress in the tumour interstitium is modelled as a consequence of the rapid proliferation of the tumour cells, using a multiplicative decomposition of the deformation gradient. In turn, the tumour cell proliferation is determined by the oxygen concentration (which is governed by the diffusion equation) and the solid stress. We apply a linear stability analysis to determine the onset of mechanical instability, and the Liapunov–Schmidt reduction to determine the postbuckling behaviour. We find that planar modes usually go unstable before axial modes, so that our model can explain the buckling of capillaries, but not as easily their tortuosity. We also find that the inclusion of anisotropic growth in our model can substantially affect the onset of buckling. Anisotropic growth also results in a feedback effect that substantially affects the magnitude of the buckle.

Dynamics and morphology characteristics of cell colonies with radially spreading growth fronts

The dynamics of two-dimensional (2D) radially spreading growth fronts of Vero cell colonies was investigated utilizing two types of colonies, namely type I starting from clusters with a small number of cells, which initially exhibited arbitrary-shaped rough growth fronts and progressively approached quasicircular ones as the cell population increased; and type II colonies, starting from a relatively large circular three-dimensional (3D) cell cluster. For large cell population colonies, the fractal dimension of the fronts was D(F) = 1.20±0.05. For low cell populations, the mean colony radius increased exponentially with time, but for large ones the constant radial front velocity 0.20±0.02 μm min(-1) was reached. Colony spreading was accompanied by changes in both cell morphology and average size, and by the formation of very large cells, some of them multinuclear. Therefore the heterogeneity of colonies increased and local driving forces that set in began to influence the 2D growth front kinetics. The retardation effect related to the exponential to constant radial front velocity transition was assigned to a number of possible interferences including the cell duplication and 3D growth in the bulk of the colony. The dynamic scaling analysis of overhang-corrected rough colony fronts, after arc-radius coordinate system transformation, resulted in roughness exponent α = 0.50±0.05 and growth exponent β = 0.32±0.04, for arc lengths greater than 100 μm. This set of scaling exponents agreed with that predicted by the Kardar, Parisi, and Zhang continuous equation. For arc lengths shorter than 2-3 cell diameters, the value α = 0.85±0.05 would be related to a cell front roughening caused by temporarily membrane deformations occasionally interfered by cell proliferation.

The role of the microenvironment in tumor growth and invasion

Mathematical modeling and computational analysis are essential for understanding the dynamics of the complex gene networks that control normal development and homeostasis, and can help to understand how circumvention of that control leads to abnormal outcomes such as cancer. Our objectives here are to discuss the different mechanisms by which the local biochemical and mechanical microenvironment, which is comprised of various signaling molecules, cell types and the extracellular matrix (ECM), affects the progression of potentially-cancerous cells, and to present new results on two aspects of these effects. We first deal with the major processes involved in the progression from a normal cell to a cancerous cell at a level accessible to a general scientific readership, and we then outline a number of mathematical and computational issues that arise in cancer modeling. In Section 2 we present results from a model that deals with the effects of the mechanical properties of the environment on tumor growth, and in Section 3 we report results from a model of the signaling pathways and the tumor microenvironment (TME), and how their interactions affect the development of breast cancer. The results emphasize anew the complexities of the interactions within the TME and their effect on tumor growth, and show that tumor progression is not solely determined by the presence of a clone of mutated immortal cells, but rather that it can be 'community-controlled'.

Investigation of tumor hypoxia using a two-enzyme system for in vitro generation of oxygen deficiency

Background

Oxygen deficiency in tumor tissue is associated with a malign phenotype, characterized by high invasiveness, increased metastatic potential and poor prognosis. Hypoxia chambers are the established standard model for in vitro studies on tumor hypoxia. An enzymatic hypoxia system (GOX/CAT) based on the use of glucose oxidase (GOX) and catalase (CAT) that allows induction of stable hypoxia for in vitro approaches more rapidly and with less operating expense has been introduced recently. Aim of this work is to compare the enzymatic system with the established technique of hypoxia chamber in respect of gene expression, glucose metabolism and radioresistance, prior to its application for in vitro investigation of oxygen deficiency.

Methods

Human head and neck squamous cell carcinoma HNO97 cells were incubated under normoxic and hypoxic conditions using both hypoxia chamber and the enzymatic model. Gene expression was investigated using Agilent microarray chips and real time PCR analysis. ¹⁴C-fluoro-deoxy-glucose uptake experiments were performed in order to evaluate cellular metabolism. Cell proliferation after photon irradiation was investigated for evaluation of radioresistance under normoxia and hypoxia using both a hypoxia chamber and the enzymatic system.

Results

The microarray analysis revealed a similar trend in the expression of known HIF-1 target genes between the two hypoxia systems for HNO97 cells. Quantitative RT-PCR demonstrated different kinetic patterns in the expression of carbonic anhydrase IX and lysyl oxidase, which might be due to the faster induction of hypoxia by the enzymatic system. ¹⁴C-fluoro-deoxy-glucose uptake assays showed a higher glucose metabolism under hypoxic conditions, especially for the enzymatic system. Proliferation experiments after photon irradiation revealed increased survival rates for the enzymatic model compared to hypoxia chamber and normoxia, indicating enhanced resistance to irradiation. While the GOX/CAT system allows independent investigation of hypoxia and oxidative stress, care must be taken to prevent acidification during longer incubation.

Conclusion

The results of our study indicate that the enzymatic model can find application for in vitro investigation of tumor hypoxia, despite limitations that need to be considered in the experimental design.