Clonal interactions in cancer: Integrating quantitative models with experimental and clinical data.

Tumors consist of different genotypically distinct subpopulations-or subclones-of cells. These subclones can influence neighboring clones in a process called "clonal interaction." Conventionally, research on driver mutations in cancer has focused on their cell-autonomous effects that lead to an increase in fitness of the cells containing the driver. Recently, with the advent of improved experimental and computational technologies for investigating tumor heterogeneity and clonal dynamics, new studies have shown the importance of clonal interactions in cancer initiation, progression, and metastasis. In this review we provide an overview of clonal interactions in cancer, discussing key discoveries from a diverse range of approaches to cancer biology research. We discuss common types of clonal interactions, such as cooperation and competition, its mechanisms, and the overall effect on tumorigenesis, with important implications for tumor heterogeneity, resistance to treatment, and tumor suppression. Quantitative models-in coordination with cell culture and animal model experiments-have played a vital role in investigating the nature of clonal interactions and the complex clonal dynamics they generate. We present mathematical and computational models that can be used to represent clonal interactions and provide examples of the roles they have played in identifying and quantifying the strength of clonal interactions in experimental systems. Clonal interactions have proved difficult to observe in clinical data; however, several very recent quantitative approaches enable their detection. We conclude by discussing ways in which researchers can further integrate quantitative methods with experimental and clinical data to elucidate the critical-and often surprising-roles of clonal interactions in human cancers.

Optimal control of combination immunotherapy for a virtual murine cohort in a glioblastoma-immune dynamics model.

The immune checkpoint inhibitor anti-PD-1, commonly used in cancer immunotherapy, has not been successful as a monotherapy for the highly aggressive brain cancer glioblastoma. However, when used in conjunction with a CC-chemokine receptor-2 (CCR2) antagonist, anti-PD-1 has shown efficacy in preclinical studies. In this paper, we aim to optimize treatment regimens for this combination immunotherapy using optimal control theory. We extend a treatment-free glioblastoma-immune dynamics ODE model to include interventions with anti-PD-1 and the CCR2 antagonist. An optimized regimen increases the survival of an average mouse from 32 days post-tumor implantation without treatment to 111 days with treatment. We scale this approach to a virtual murine cohort to evaluate mortality and quality of life concerns during treatment, and predict survival, tumor recurrence, or death after treatment. A parameter identifiability analysis identifies five parameters suitable for personalizing treatment within the virtual cohort. Sampling from these five practically identifiable parameters for the virtual murine cohort reveals that personalized, optimized regimens enhance survival: 84% of the virtual mice survive to day 100, compared to 60% survival in a previously studied experimental regimen. Subjects with high tumor growth rates and low T cell kill rates are identified as more likely to die during and after treatment due to their compromised immune systems and more aggressive tumors. Notably, the MDSC death rate emerges as a long-term predictor of either disease-free survival or death.

Harnessing Flex Point Symmetry to Estimate Logistic Tumor Population Growth.

The observed time evolution of a population is well approximated by a logistic growth function in many research fields, including oncology, ecology, chemistry, demography, economy, linguistics, and artificial neural networks. Initial growth is exponential, then decelerates as the population approaches its limit size, i.e., the carrying capacity. In mathematical oncology, the tumor carrying capacity has been postulated to be dynamically evolving as the tumor overcomes several evolutionary bottlenecks and, thus, to be patient specific. As the relative tumor-over-carrying capacity ratio may be predictive and prognostic for tumor growth and treatment response dynamics, it is paramount to estimate it from limited clinical data. We show that exploiting the logistic function's rotation symmetry can help estimate the population's growth rate and carry capacity from fewer data points than conventional regression approaches. We test this novel approach against published pan-cancer animal and human breast cancer data, achieving a 30% to 40% reduction in the time at which subsequent data collection is necessary to estimate the logistic growth rate and carrying capacity correctly. These results could improve tumor dynamics forecasting and augment the clinical decision-making process.

Characterising Cancer Cell Responses to Cyclic Hypoxia Using Mathematical Modelling.

In vivo observations show that oxygen levels in tumours can fluctuate on fast and slow timescales. As a result, cancer cells can be periodically exposed to pathologically low oxygen levels; a phenomenon known as cyclic hypoxia. Yet, little is known about the response and adaptation of cancer cells to cyclic, rather than, constant hypoxia. Further, existing in vitro models of cyclic hypoxia fail to capture the complex and heterogeneous oxygen dynamics of tumours growing in vivo. Mathematical models can help to overcome current experimental limitations and, in so doing, offer new insights into the biology of tumour cyclic hypoxia by predicting cell responses to a wide range of cyclic dynamics. We develop an individual-based model to investigate how cell cycle progression and cell fate determination of cancer cells are altered following exposure to cyclic hypoxia. Our model can simulate standard in vitro experiments, such as clonogenic assays and cell cycle experiments, allowing for efficient screening of cell responses under a wide range of cyclic hypoxia conditions. Simulation results show that the same cell line can exhibit markedly different responses to cyclic hypoxia depending on the dynamics of the oxygen fluctuations. We also use our model to investigate the impact of changes to cell cycle checkpoint activation and damage repair on cell responses to cyclic hypoxia. Our simulations suggest that cyclic hypoxia can promote heterogeneity in cellular damage repair activity within vascular tumours.

Assessing the Role of Patient Generation Techniques in Virtual Clinical Trial Outcomes.

Virtual clinical trials (VCTs) are growing in popularity as a tool for quantitatively predicting heterogeneous treatment responses across a population. In the context of a VCT, a plausible patient is an instance of a mathematical model with parameter (or attribute) values chosen to reflect features of the disease and response to treatment for that particular patient. A number of techniques have been introduced to determine the set of model parametrizations to include in a virtual patient cohort. These methodologies generally start with a prior distribution for each model parameter and utilize some criteria to determine whether a parameter set sampled from the priors should be included or excluded from the plausible population. No standard technique exists, however, for generating these prior distributions and choosing the inclusion/exclusion criteria. In this work, we rigorously quantify the impact that VCT design choices have on VCT predictions. Rather than use real data and a complex mathematical model, a spatial model of radiotherapy is used to generate simulated patient data and the mathematical model used to describe the patient data is a two-parameter ordinary differential equations model. This controlled setup allows us to isolate the impact of both the prior distribution and the inclusion/exclusion criteria on both the heterogeneity of plausible populations and on predicted treatment response. We find that the prior distribution, rather than the inclusion/exclusion criteria, has a larger impact on the heterogeneity of the plausible population. Yet, the percent of treatment responders in the plausible population was more sensitive to the inclusion/exclusion criteria utilized. This foundational understanding of the role of virtual clinical trial design should help inform the development of future VCTs that use more complex models and real data.

Tumor containment: a more general mathematical analysis.

Clinical and pre-clinical data suggest that treating some tumors at a mild, patient-specific dose might delay resistance to treatment and increase survival time. A recent mathematical model with sensitive and resistant tumor cells identified conditions under which a treatment aiming at tumor containment rather than eradication is indeed optimal. This model however neglected mutations from sensitive to resistant cells, and assumed that the growth-rate of sensitive cells is non-increasing in the size of the resistant population. The latter is not true in standard models of chemotherapy. This article shows how to dispense with this assumption and allow for mutations from sensitive to resistant cells. This is achieved by a novel mathematical analysis comparing tumor sizes across treatments not as a function of time, but as a function of the resistant population size.

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.

Deep learning characterization of brain tumours with diffusion weighted imaging.

Glioblastoma multiforme (GBM) is one of the most deadly forms of cancer. Methods of characterizing these tumours are valuable for improving predictions of their progression and response to treatment. A mathematical model called the proliferation-invasion (PI) model has been used extensively in the literature to model the growth of these tumours, though it relies on known values of two key parameters: the tumour cell diffusivity and proliferation rate. Unfortunately, these parameters are difficult to estimate in a patient-specific manner, making personalized tumour forecasting challenging. In this paper, we develop and apply a deep learning model capable of making accurate estimates of these key GBM-characterizing parameters while simultaneously producing a full prediction of the tumour progression curve. Our method uses two sets of multi sequence MRI in order to produce estimations and relies on a preprocessing pipeline which includes brain tumour segmentation and conversion to tumour cellularity. We first apply our deep learning model to synthetic tumours to showcase the model's capabilities and identify situations where prediction errors are likely to occur. We then apply our model to a clinical dataset consisting of five patients diagnosed with GBM. For all patients, we derive evidence-based estimates for each of the PI model parameters and predictions for the future progression of the tumour, along with estimates of the parameter uncertainties. Our work provides a new, easily generalizable method for the estimation of patient-specific tumour parameters, which can be built upon to aid physicians in designing personalized treatments.

Modelling Radiation Cancer Treatment with a Death-Rate Term in Ordinary and Fractional Differential Equations.

Fractional calculus has recently been applied to the mathematical modelling of tumour growth, but its use introduces complexities that may not be warranted. Mathematical modelling with differential equations is a standard approach to study and predict treatment outcomes for population-level and patient-specific responses. Here, we use patient data of radiation-treated tumours to discuss the benefits and limitations of introducing fractional derivatives into three standard models of tumour growth. The fractional derivative introduces a history-dependence into the growth function, which requires a continuous death-rate term for radiation treatment. This newly proposed radiation-induced death-rate term improves computational efficiency in both ordinary and fractional derivative models. This computational speed-up will benefit common simulation tasks such as model parameterization and the construction and running of virtual clinical trials.

Tumor Evolution Models of Phase-Field Type with Nonlocal Effects and Angiogenesis.

In this survey article, a variety of systems modeling tumor growth are discussed. In accordance with the hallmarks of cancer, the described models incorporate the primary characteristics of cancer evolution. Specifically, we focus on diffusive interface models and follow the phase-field approach that describes the tumor as a collection of cells. Such systems are based on a multiphase approach that employs constitutive laws and balance laws for individual constituents. In mathematical oncology, numerous biological phenomena are involved, including temporal and spatial nonlocal effects, complex nonlinearities, stochasticity, and mixed-dimensional couplings. Using the models, for instance, we can express angiogenesis and cell-to-matrix adhesion effects. Finally, we offer some methods for numerically approximating the models and show simulations of the tumor's evolution in response to various biological effects.

Global stability and parameter analysis reinforce therapeutic targets of PD-L1-PD-1 and MDSCs for glioblastoma.

Glioblastoma (GBM) is an aggressive primary brain cancer that currently has minimally effective treatments. Like other cancers, immunosuppression by the PD-L1-PD-1 immune checkpoint complex is a prominent axis by which glioma cells evade the immune system. Myeloid-derived suppressor cells (MDSCs), which are recruited to the glioma microenviroment, also contribute to the immunosuppressed GBM microenvironment by suppressing T cell functions. In this paper, we propose a GBM-specific tumor-immune ordinary differential equations model of glioma cells, T cells, and MDSCs to provide theoretical insights into the interactions between these cells. Equilibrium and stability analysis indicates that there are unique tumorous and tumor-free equilibria which are locally stable under certain conditions. Further, the tumor-free equilibrium is globally stable when T cell activation and the tumor kill rate by T cells overcome tumor growth, T cell inhibition by PD-L1-PD-1 and MDSCs, and the T cell death rate. Bifurcation analysis suggests that a treatment plan that includes surgical resection and therapeutics targeting immune suppression caused by the PD-L1-PD1 complex and MDSCs results in the system tending to the tumor-free equilibrium. Using a set of preclinical experimental data, we implement the approximate Bayesian computation (ABC) rejection method to construct probability density distributions that estimate model parameters. These distributions inform an appropriate search curve for global sensitivity analysis using the extended fourier amplitude sensitivity test. Sensitivity results combined with the ABC method suggest that parameter interaction is occurring between the drivers of tumor burden, which are the tumor growth rate and carrying capacity as well as the tumor kill rate by T cells, and the two modeled forms of immunosuppression, PD-L1-PD-1 immune checkpoint and MDSC suppression of T cells. Thus, treatment with an immune checkpoint inhibitor in combination with a therapeutic targeting the inhibitory mechanisms of MDSCs should be explored.

Spatial cumulant models enable spatially informed treatment strategies and analysis of local interactions in cancer systems.

Theoretical and applied cancer studies that use individual-based models (IBMs) have been limited by the lack of a mathematical formulation that enables rigorous analysis of these models. However, spatial cumulant models (SCMs), which have arisen from theoretical ecology, describe population dynamics generated by a specific family of IBMs, namely spatio-temporal point processes (STPPs). SCMs are spatially resolved population models formulated by a system of differential equations that approximate the dynamics of two STPP-generated summary statistics: first-order spatial cumulants (densities), and second-order spatial cumulants (spatial covariances). We exemplify how SCMs can be used in mathematical oncology by modelling theoretical cancer cell populations comprising interacting growth factor-producing and non-producing cells. To formulate model equations, we use computational tools that enable the generation of STPPs, SCMs and mean-field population models (MFPMs) from user-defined model descriptions (Cornell et al. Nat Commun 10:4716, 2019). To calculate and compare STPP, SCM and MFPM-generated summary statistics, we develop an application-agnostic computational pipeline. Our results demonstrate that SCMs can capture STPP-generated population density dynamics, even when MFPMs fail to do so. From both MFPM and SCM equations, we derive treatment-induced death rates required to achieve non-growing cell populations. When testing these treatment strategies in STPP-generated cell populations, our results demonstrate that SCM-informed strategies outperform MFPM-informed strategies in terms of inhibiting population growths. We thus demonstrate that SCMs provide a new framework in which to study cell-cell interactions, and can be used to describe and perturb STPP-generated cell population dynamics. We, therefore, argue that SCMs can be used to increase IBMs' applicability in cancer research.

Modelling microtube driven invasion of glioma.

Malignant gliomas are notoriously invasive, a major impediment against their successful treatment. This invasive growth has motivated the use of predictive partial differential equation models, formulated at varying levels of detail, and including (i) "proliferation-infiltration" models, (ii) "go-or-grow" models, and (iii) anisotropic diffusion models. Often, these models use macroscopic observations of a diffuse tumour interface to motivate a phenomenological description of invasion, rather than performing a detailed and mechanistic modelling of glioma cell invasion processes. Here we close this gap. Based on experiments that support an important role played by long cellular protrusions, termed tumour microtubes, we formulate a new model for microtube-driven glioma invasion. In particular, we model a population of tumour cells that extend tissue-infiltrating microtubes. Mitosis leads to new nuclei that migrate along the microtubes and settle elsewhere. A combination of steady state analysis and numerical simulation is employed to show that the model can predict an expanding tumour, with travelling wave solutions led by microtube dynamics. A sequence of scaling arguments allows us reduce the detailed model into simpler formulations, including models falling into each of the general classes (i), (ii), and (iii) above. This analysis allows us to clearly identify the assumptions under which these various models can be a posteriori justified in the context of microtube-driven glioma invasion. Numerical simulations are used to compare the various model classes and we discuss their advantages and disadvantages.

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.

Can the Kuznetsov Model Replicate and Predict Cancer Growth in Humans?

Several mathematical models to predict tumor growth over time have been developed in the last decades. A central aspect of such models is the interaction of tumor cells with immune effector cells. The Kuznetsov model (Kuznetsov et al. in Bull Math Biol 56(2):295-321, 1994) is the most prominent of these models and has been used as a basis for many other related models and theoretical studies. However, none of these models have been validated with large-scale real-world data of human patients treated with cancer immunotherapy. In addition, parameter estimation of these models remains a major bottleneck on the way to model-based and data-driven medical treatment. In this study, we quantitatively fit Kuznetsov's model to a large dataset of 1472 patients, of which 210 patients have more than six data points, by estimating the model parameters of each patient individually. We also conduct a global practical identifiability analysis for the estimated parameters. We thus demonstrate that several combinations of parameter values could lead to accurate data fitting. This opens the potential for global parameter estimation of the model, in which the values of all or some parameters are fixed for all patients. Furthermore, by omitting the last two or three data points, we show that the model can be extrapolated and predict future tumor dynamics. This paves the way for a more clinically relevant application of mathematical tumor modeling, in which the treatment strategy could be adjusted in advance according to the model's future predictions.

Mathematical Modeling of Tumor and Cancer Stem Cells Treated with CAR-T Therapy and Inhibition of TGF-[Formula: see text].

The stem cell hypothesis suggests that there is a small group of malignant cells, the cancer stem cells, that initiate the development of tumors, encourage its growth, and may even be the cause of metastases. Traditional treatments, such as chemotherapy and radiation, primarily target the tumor cells leaving the stem cells to potentially cause a recurrence. Chimeric antigen receptor (CAR) T-cell therapy is a form of immunotherapy where the immune cells are genetically modified to fight the tumor cells. Traditionally, the CAR T-cell therapy has been used to treat blood cancers and only recently has shown promising results against solid tumors. We create an ordinary differential equations model which allows for the infusion of trained CAR-T cells to specifically attack the cancer stem cells that are present in the solid tumor microenvironment. Additionally, we incorporate the influence of TGF-[Formula: see text] which inhibits the CAR-T cells and thus promotes the growth of the tumor. We verify the model by comparing it to available data and then examine combinations of CAR-T cell treatment targeting both non-stem and stem cancer cells and a treatment that reduces the effectiveness of TGF-[Formula: see text] to determine the scenarios that eliminate the tumor.

The Tumor Invasion Paradox in Cancer Stem Cell-Driven Solid Tumors.

Cancer stem cells (CSCs) are key in understanding tumor growth and tumor progression. A counterintuitive effect of CSCs is the so-called tumor growth paradox: the effect where a tumor with a higher death rate may grow larger than a tumor with a lower death rate. Here we extend the modeling of the tumor growth paradox by including spatial structure and considering cancer invasion. Using agent-based modeling and a corresponding partial differential equation model, we demonstrate and prove mathematically a tumor invasion paradox: a larger cell death rate can lead to a faster invasion speed. We test this result on a generic hypothetical cancer with typical growth rates and typical treatment sensitivities. We find that the tumor invasion paradox may play a role for continuous and intermittent treatments, while it does not seem to be essential in fractionated treatments. It should be noted that no attempt was made to fit the model to a specific cancer, thus, our results are generic and theoretical.

Temporal optimization of radiation therapy to heterogeneous tumour populations and cancer stem cells.

External beam radiation therapy is a key part of modern cancer treatments which uses high doses of radiation to destroy tumour cells. Despite its widespread usage and extensive study in theoretical, experimental, and clinical works, many questions still remain about how best to administer it. Many mathematical studies have examined optimal scheduling of radiotherapy, and most come to similar conclusions. Importantly though, these studies generally assume intratumoral homogeneity. But in recent years, it has become clear that tumours are not homogeneous masses of cancerous cells, but wildly heterogeneous masses with various subpopulations which grow and respond to treatment differently. One subpopulation of particular importance is cancer stem cells (CSCs) which are known to exhibit higher radioresistence compared with non-CSCs. Knowledge of these differences between cell types could theoretically lead to changes in optimal treatment scheduling. Only a few studies have examined this question, and interestingly, they arrive at apparent conflicting results. However, an understanding of their assumptions reveals a key difference which leads to their differing conclusions. In this paper, we generalize the problem of temporal optimization of dose distribution of radiation therapy to a two cell type model. We do so by creating a mathematical model and a numerical optimization algorithm to find the distribution of dose which leads to optimal cell kill. We then create a data set of optimization solutions and use data analysis tools to learn the relationships between model parameters and the qualitative behaviour of optimization results. Analysis of the model and discussion of biological importance are provided throughout. We find that the key factor in predicting the behaviour of the optimal distribution of radiation is the ratio between the radiosensitivities of the present cell types. These results can provide guidance for treatment in cases where clinicians have knowledge of tumour heterogeneity and of the abundance of CSCs.

Computer simulations suggest that prostate enlargement due to benign prostatic hyperplasia mechanically impedes prostate cancer growth.

Prostate cancer and benign prostatic hyperplasia are common genitourinary diseases in aging men. Both pathologies may coexist and share numerous similarities, which have suggested several connections or some interplay between them. However, solid evidence confirming their existence is lacking. Recent studies on extensive series of prostatectomy specimens have shown that tumors originating in larger prostates present favorable pathological features. Hence, large prostates may exert a protective effect against prostate cancer. In this work, we propose a mechanical explanation for this phenomenon. The mechanical stress fields that originate as tumors enlarge have been shown to slow down their dynamics. Benign prostatic hyperplasia contributes to these mechanical stress fields, hence further restraining prostate cancer growth. We derived a tissue-scale, patient-specific mechanically coupled mathematical model to qualitatively investigate the mechanical interaction of prostate cancer and benign prostatic hyperplasia. This model was calibrated by studying the deformation caused by each disease independently. Our simulations show that a history of benign prostatic hyperplasia creates mechanical stress fields in the prostate that impede prostatic tumor growth and limit its invasiveness. The technology presented herein may assist physicians in the clinical management of benign prostate hyperplasia and prostate cancer by predicting pathological outcomes on a tissue-scale, patient-specific basis.

A Century of Fractionated Radiotherapy: How Mathematical Oncology Can Break the Rules.

Radiotherapy is involved in 50% of all cancer treatments and 40% of cancer cures. Most of these treatments are delivered in fractions of equal doses of radiation (Fractional Equivalent Dosing (FED)) in days to weeks. This treatment paradigm has remained unchanged in the past century and does not account for the development of radioresistance during treatment. Even if under-optimized, deviating from a century of successful therapy delivered in FED can be difficult. One way of exploring the infinite space of fraction size and scheduling to identify optimal fractionation schedules is through mathematical oncology simulations that allow for in silico evaluation. This review article explores the evidence that current fractionation promotes the development of radioresistance, summarizes mathematical solutions to account for radioresistance, both in the curative and non-curative setting, and reviews current clinical data investigating non-FED fractionated radiotherapy.