Integrative mathematical oncology

A mathematical framework for modelling 3D cell motility: applications to glioblastoma cell migration

The collection of 3D cell tracking data from live images of micro-tissues is a recent innovation made possible due to advances in imaging techniques. As such there is increased interest in studying cell motility in 3D in vitro model systems but a lack of rigorous methodology for analysing the resulting data sets. One such instance of the use of these in vitro models is in the study of cancerous tumours. Growing multicellular tumour spheroids in vitro allows for modelling of the tumour microenvironment and the study of tumour cell behaviours, such as migration, which improves understanding of these cells and in turn could potentially improve cancer treatments. In this paper, we present a workflow for the rigorous analysis of 3D cell tracking data, based on the persistent random walk model, but adaptable to other biologically informed mathematical models. We use statistical measures to assess the fit of the model to the motility data and to estimate model parameters and provide confidence intervals for those parameters, to allow for parametrization of the model taking correlation in the data into account. We use in silico simulations to validate the workflow in 3D before testing our method on cell tracking data taken from in vitro experiments on glioblastoma tumour cells, a brain cancer with a very poor prognosis. The presented approach is intended to be accessible to both modellers and experimentalists alike in that it provides tools for uncovering features of the data set that may suggest amendments to future experiments or modelling attempts.

A Mathematical Study of the Influence of Hypoxia and Acidity on the Evolutionary Dynamics of Cancer

Hypoxia and acidity act as environmental stressors promoting selection for cancer cells with a more aggressive phenotype. As a result, a deeper theoretical understanding of the spatio-temporal processes that drive the adaptation of tumour cells to hypoxic and acidic microenvironments may open up new avenues of research in oncology and cancer treatment. We present a mathematical model to study the influence of hypoxia and acidity on the evolutionary dynamics of cancer cells in vascularised tumours. The model is formulated as a system of partial integro-differential equations that describe the phenotypic evolution of cancer cells in response to dynamic variations in the spatial distribution of three abiotic factors that are key players in tumour metabolism: oxygen, glucose and lactate. The results of numerical simulations of a calibrated version of the model based on real data recapitulate the eco-evolutionary spatial dynamics of tumour cells and their adaptation to hypoxic and acidic microenvironments. Moreover, such results demonstrate how nonlinear interactions between tumour cells and abiotic factors can lead to the formation of environmental gradients which select for cells with phenotypic characteristics that vary with distance from intra-tumour blood vessels, thus promoting the emergence of intra-tumour phenotypic heterogeneity. Finally, our theoretical findings reconcile the conclusions of earlier studies by showing that the order in which resistance to hypoxia and resistance to acidity arise in tumours depend on the ways in which oxygen and lactate act as environmental stressors in the evolutionary dynamics of cancer cells.

Superlinear growth reveals the Allee effect in tumors

Integrating experimental data into ecological models plays a central role in understanding biological mechanisms that drive tumor progression where such knowledge can be used to develop new therapeutic strategies. While the current studies emphasize the role of competition among tumor cells, they fail to explain recently observed superlinear growth dynamics across human tumors. Here we study tumor growth dynamics by developing a model that incorporates evolutionary dynamics inside tumors with tumor-microenvironment interactions. Our results reveal that tumor cells' ability to manipulate the environment and induce angiogenesis drives superlinear growth-a process compatible with the Allee effect. In light of this understanding, our model suggests that, for high-risk tumors that have a higher growth rate, suppressing angiogenesis can be the appropriate therapeutic intervention.

Identifying key questions in the ecology and evolution of cancer

The application of evolutionary and ecological principles to cancer prevention and treatment, as well as recognizing cancer as a selection force in nature, has gained impetus over the last 50 years. Following the initial theoretical approaches that combined knowledge from interdisciplinary fields, it became clear that using the eco-evolutionary framework is of key importance to understand cancer. We are now at a pivotal point where accumulating evidence starts to steer the future directions of the discipline and allows us to underpin the key challenges that remain to be addressed. Here, we aim to assess current advancements in the field and to suggest future directions for research. First, we summarize cancer research areas that, so far, have assimilated ecological and evolutionary principles into their approaches and illustrate their key importance. Then, we assembled 33 experts and identified 84 key questions, organized around nine major themes, to pave the foundations for research to come. We highlight the urgent need for broadening the portfolio of research directions to stimulate novel approaches at the interface of oncology and ecological and evolutionary sciences. We conclude that progressive and efficient cross-disciplinary collaborations that draw on the expertise of the fields of ecology, evolution and cancer are essential in order to efficiently address current and future questions about cancer.

Computational modeling reveals a key role for polarized myeloid cells in controlling osteoclast activity during bone injury repair

Bone-forming osteoblasts and -resorbing osteoclasts control bone injury repair, and myeloid-derived cells such as monocytes and macrophages are known to influence their behavior. However, precisely how these multiple cell types coordinate and regulate each other over time within the bone marrow to restore bone is difficult to dissect using biological approaches. Conversely, mathematical modeling lends itself well to this challenge. Therefore, we generated an ordinary differential equation (ODE) model powered by experimental data (osteoblast, osteoclast, bone volume, pro- and anti-inflammatory myeloid cells) obtained from intra-tibially injured mice. Initial ODE results using only osteoblast/osteoclast populations demonstrated that bone homeostasis could not be recovered after injury, but this issue was resolved upon integration of pro- and anti-inflammatory myeloid population dynamics. Surprisingly, the ODE revealed temporal disconnects between the peak of total bone mineralization/resorption, and osteoblast/osteoclast numbers. Specifically, the model indicated that osteoclast activity must vary greatly (> 17-fold) to return the bone volume to baseline after injury and suggest that osteoblast/osteoclast number alone is insufficient to predict bone the trajectory of bone repair. Importantly, the values of osteoclast activity fall within those published previously. These data underscore the value of mathematical modeling approaches to understand and reveal new insights into complex biological processes.

Fractional Mathematical Oncology: On the potential of non-integer order calculus applied to interdisciplinary models

Mathematical Oncology investigates cancer-related phenomena through mathematical models as comprehensive as possible. Accordingly, an interdisciplinary approach involving concepts from biology to materials science can provide a deeper understanding of biological systems pertaining the disease. In this context, fractional calculus (also referred to as non-integer order) is a branch in mathematical analysis whose tools can describe complex phenomena comprising different time and space scales. Fractional-order models may allow a better description and understanding of oncological particularities, potentially contributing to decision-making in areas of interest such as tumor evolution, early diagnosis techniques and personalized treatment therapies. By following a phenomenological (i.e. mechanistic) approach, the present study surveys and explores different aspects of Fractional Mathematical Oncology, reviewing and discussing recent developments in view of their prospective applications.

Engineering in Medicine To Address the Challenge of Cancer Drug Resistance: From Micro- and Nanotechnologies to Computational and Mathematical Modeling

Drug resistance has profoundly limited the success of cancer treatment, driving relapse, metastasis, and mortality. Nearly all anticancer drugs and even novel immunotherapies, which recalibrate the immune system for tumor recognition and destruction, have succumbed to resistance development. Engineers have emerged across mechanical, physical, chemical, mathematical, and biological disciplines to address the challenge of drug resistance using a combination of interdisciplinary tools and skill sets. This review explores the developing, complex, and under-recognized role of engineering in medicine to address the multitude of challenges in cancer drug resistance. Looking through the "lens" of intrinsic, extrinsic, and drug-induced resistance (also referred to as "tolerance"), we will discuss three specific areas where active innovation is driving novel treatment paradigms: (1) nanotechnology, which has revolutionized drug delivery in desmoplastic tissues, harnessing physiochemical characteristics to destroy tumors through photothermal therapy and rationally designed nanostructures to circumvent cancer immunotherapy failures, (2) bioengineered tumor models, which have benefitted from microfluidics and mechanical engineering, creating a paradigm shift in physiologically relevant environments to predict clinical refractoriness and enabling platforms for screening drug combinations to thwart resistance at the individual patient level, and (3) computational and mathematical modeling, which blends in silico simulations with molecular and evolutionary principles to map mutational patterns and model interactions between cells that promote resistance. On the basis that engineering in medicine has resulted in discoveries in resistance biology and successfully translated to clinical strategies that improve outcomes, we suggest the proliferation of multidisciplinary science that embraces engineering.

Evidence for hypoxia increasing the tempo of evolution in glioblastoma

BACKGROUND

Tumour hypoxia is associated with metastatic disease, and while there have been many mechanisms proposed for why tumour hypoxia is associated with metastatic disease, it remains unclear whether one precise mechanism is the key reason or several in concert. Somatic evolution drives cancer progression and treatment resistance, fuelled not only by genetic and epigenetic mutation but also by selection from interactions between tumour cells, normal cells and physical micro-environment. Ecological habitats influence evolutionary dynamics, but the impact on tempo of evolution is less clear.

METHODS

We explored this complex dialogue with a combined clinical-theoretical approach by simulating a proliferative hierarchy under heterogeneous oxygen availability with an agent-based model. Predictions were compared against histology samples taken from glioblastoma patients, stained to elucidate areas of necrosis and TP53 expression heterogeneity.

RESULTS

Results indicate that cell division in hypoxic environments is effectively upregulated, with low-oxygen niches providing avenues for tumour cells to spread. Analysis of human data indicates that cell division is not decreased under hypoxia, consistent with our results.

CONCLUSIONS

Our results suggest that hypoxia could be a crucible that effectively warps evolutionary velocity, making key mutations more likely. Thus, key tumour ecological niches such as hypoxic regions may alter the evolutionary tempo, driving mutations fuelling tumour heterogeneity.

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

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

High School Internship Program in Integrated Mathematical Oncology (HIP IMO): Five-Year Experience at Moffitt Cancer Center

Modern cancer research, and the wealth of data across multiple spatial and temporal scales, has created the need for researchers that are well versed in the life sciences (cancer biology, developmental biology, immunology), medical sciences (oncology) and natural sciences (mathematics, physics, engineering, computer sciences). College undergraduate education traditionally occurs in disciplinary silos, which creates a steep learning curve at the graduate and postdoctoral levels that increasingly bridge multiple disciplines. Numerous colleges have begun to embrace interdisciplinary curricula, but students who double major in mathematics (or other quantitative sciences) and biology (or medicine) remain scarce. We identified the need to educate junior and senior high school students about integrating mathematical and biological skills, through the lens of mathematical oncology, to better prepare students for future careers at the interdisciplinary interface. The High school Internship Program in Integrated Mathematical Oncology (HIP IMO) at Moffitt Cancer Center has so far trained 59 students between 2015 and 2019. We report here on the program structure, training deliverables, curriculum and outcomes. We hope to promote interdisciplinary educational activities early in a student's career.

Modeling Cancer Dynamics

Cancer is a complex disease that continues to pose formidable challenges to therapeutic interventions. An increased understanding of cancer complexity and in particular tumor growth dynamics is critical to the development of more effective therapies. In this respect, an agent-based model of tumor growth is explored with the consideration of cancer reprograming of metabolism and the immune response, to seek insight about the coupling between these two key determinants of tumor growth dynamics. Ultimately, this exploration is intended to inform the development of therapies that can induce a more effective immune response despite the metabolic constraints of the tumor microenvironment.

Meeting the Needs of A Changing Landscape: Advances and Challenges in Undergraduate Biology Education

Over the last 25 years, reforms in undergraduate biology education have transformed the way biology is taught at many institutions of higher education. This has been fueled in part by a burgeoning discipline-based education research community, which has advocated for evidence-based instructional practices based on findings from research. This perspective will review some of the changes to undergraduate biology education that have gained or are currently gaining momentum, becoming increasingly common in undergraduate biology classrooms. However, there are still areas in need of improvement. Although more underrepresented minority students are enrolling in and graduating from biology programs than in the past, there is a need to understand the experiences and broaden participation of other underserved groups in biology and ensure biology classroom learning environments are inclusive. Additionally, although understanding biology relies on understanding concepts from the physical sciences and mathematics, students still rarely connect the concepts they learn from other STEM disciplines to biology. Integrating concepts and practices across the STEM disciplines will be critical for biology graduates as they tackle the biological problems of the twenty-first century.

Self-organization in brain tumors: How cell morphology and cell density influence glioma pattern formation

Modeling cancer cells is essential to better understand the dynamic nature of brain tumors and glioma cells, including their invasion of normal brain. Our goal is to study how the morphology of the glioma cell influences the formation of patterns of collective behavior such as flocks (cells moving in the same direction) or streams (cells moving in opposite direction) referred to as oncostream. We have observed experimentally that the presence of oncostreams correlates with tumor progression. We propose an original agent-based model that considers each cell as an ellipsoid. We show that stretching cells from round to ellipsoid increases stream formation. A systematic numerical investigation of the model was implemented in [Formula: see text]. We deduce a phase diagram identifying key regimes for the dynamics (e.g. formation of flocks, streams, scattering). Moreover, we study the effect of cellular density and show that, in contrast to classical models of flocking, increasing cellular density reduces the formation of flocks. We observe similar patterns in [Formula: see text] with the noticeable difference that stream formation is more ubiquitous compared to flock formation.

How and why to build a mathematical model: A case study using prion aggregation

Biological systems are inherently complex, and the increasing level of detail with which we are able to experimentally probe such systems continually reveals new complexity. Fortunately, mathematical models are uniquely positioned to provide a tool suitable for rigorous analysis, hypothesis generation, and connecting results from isolated in vitro experiments with results from in vivo and whole-organism studies. However, developing useful mathematical models is challenging because of the often different domains of knowledge required in both math and biology. In this work, we endeavor to provide a useful guide for researchers interested in incorporating mathematical modeling into their scientific process. We advocate for the use of conceptual diagrams as a starting place to anchor researchers from both domains. These diagrams are useful for simplifying the biological process in question and distinguishing the essential components. Not only do they serve as the basis for developing a variety of mathematical models, but they ensure that any mathematical formulation of the biological system is led primarily by scientific questions. We provide a specific example of this process from our own work in studying prion aggregation to show the power of mathematical models to synergistically interact with experiments and push forward biological understanding. Choosing the most suitable model also depends on many different factors, and we consider how to make these choices based on different scales of biological organization and available data. We close by discussing the many opportunities that abound for both experimentalists and modelers to take advantage of collaborative work in this field.

Mathematical oncology and it's application in non melanoma skin cancer - A primer for radiation oncology professionals

Cancers of the skin (the majority of which are basal and squamous cell skin carcinomas, but also include the rarer Merkel cell carcinoma) are overwhelmingly the most common of all types of cancer. Most of these are treated surgically, with radiation reserved for those patients with high risk features or anatomical locations less suitable for surgery. Given the high incidence of both basal and squamous cell carcinomas, as well as the relatively poor outcome for Merkel cell carcinoma, it is useful to investigate the role of other disciplines regarding their diagnosis, staging and treatment. Mathematical modelling is one such area of investigation. The use of mathematical modelling is a relatively recent addition to the armamentarium of cancer treatment. It has long been recognised that tumour growth and treatment response is a complex, non-linear biological phenomenon with many mechanisms yet to be understood. Despite decades of research, including clinical, population and basic science approaches, we continue to be challenged by the complexity, heterogeneity and adaptability of tumours, both in individual patients in the oncology clinic and across wider patient populations. Prospective clinical trials predominantly focus on average outcome, with little understanding as to why individual patients may or may not respond. The use of mathematical models may lead to a greater understanding of tumour initiation, growth dynamics and treatment response.

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

MOTIVATION

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

RESULTS

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

AVAILABILITY AND IMPLEMENTATION

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

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A case for the reuse and adaptation of mechanistic computational models to study transplant immunology

Computational mechanistic models constitute powerful tools for summarizing our knowledge in quantitative terms, providing mechanistic understanding, and generating new hypotheses. The present review emphasizes the advantages of reusing publicly available computational models as a way to capitalize on existing knowledge, reduce the number of parameters that need to be adjusted to experimental data, and facilitate hypothesis generation. Finally, it includes a step-by-step example of the reuse and adaptation of an existing model of immune responses to tuberculosis, tumor growth, and blood pathogens, to study donor-specific antibody (DSA) responses. This review aims to illustrate the benefit of leveraging the currently available computational models in immunology to accelerate the study of alloimmune responses, and to encourage modelers to share their models to further advance our understanding of transplant immunology.

Optimizing circadian drug infusion schedules towards personalized cancer chronotherapy

Precision medicine requires accurate technologies for drug administration and proper systems pharmacology approaches for patient data analysis. Here, plasma pharmacokinetics (PK) data of the OPTILIV trial in which cancer patients received oxaliplatin, 5-fluorouracil and irinotecan via chronomodulated schedules delivered by an infusion pump into the hepatic artery were mathematically investigated. A pump-to-patient model was designed in order to accurately represent the drug solution dynamics from the pump to the patient blood. It was connected to semi-mechanistic PK models to analyse inter-patient variability in PK parameters. Large time delays of up to 1h41 between the actual pump start and the time of drug detection in patient blood was predicted by the model and confirmed by PK data. Sudden delivery spike in the patient artery due to glucose rinse after drug administration accounted for up to 10.7% of the total drug dose. New model-guided delivery profiles were designed to precisely lead to the drug exposure intended by clinicians. Next, the complete mathematical framework achieved a very good fit to individual time-concentration PK profiles and concluded that inter-subject differences in PK parameters was the lowest for irinotecan, intermediate for oxaliplatin and the largest for 5-fluorouracil. Clustering patients according to their PK parameter values revealed patient subgroups for each drug in which inter-patient variability was largely decreased compared to that in the total population. This study provides a complete mathematical framework to optimize drug infusion pumps and inform on inter-patient PK variability, a step towards precise and personalized cancer chronotherapy.

Accuracy and safety of infusion pumps remain a critical issue in the clinics and the development of accurate mathematical models to optimize drug administration though such devices has a key part to play in the advancement of precision medicine. Here, PK data from cancer patient receiving irinotecan, oxaliplatin and 5-fluorouracil into the hepatic artery via an infusion pump was mathematically investigated. A pump-to-patient model was designed and revealed significant inconsistencies between intended drug profiles and actual plasma concentrations. This mathematical model was then used to suggest improved profiles in order to minimise error and optimise delivery. Physiologically-based PK models of the three drugs were then linked to the pump-to-patient model. The whole framework achieved a very good fit to data and allowed quantifying inter-patient variability in PK parameters and linking them to potential clinical biomarkers via patient clustering. The developed methodology improves our understanding of patient-specific drug pharmacokinetics towards personalized drug administration.

Close Encounters of the Cell Kind: The Impact of Contact Inhibition on Tumour Growth and Cancer Models

Cancer is a complex phenomenon, and the sheer variation in behaviour across different types renders it difficult to ascertain underlying biological mechanisms. Experimental approaches frequently yield conflicting results for myriad reasons, and mathematical modelling of cancer is a vital tool to explore what we cannot readily measure, and ultimately improve treatment and prognosis. Like experiments, models are underpinned by certain biological assumptions, variation of which can lead to divergent predictions. An outstanding and important question concerns contact inhibition of proliferation (CIP), the observation that proliferation ceases when cells are spatially confined by their neighbours. CIP is a characteristic of many healthy adult tissues, but it remains unclear to which extent it holds in solid tumours, which exhibit regions of hyper-proliferation, and apparent breakdown of CIP. What precisely occurs in tumour tissue remains an open question, which mathematical modelling can help shed light on. In this perspective piece, we explore the implications of different hypotheses and available experimental evidence to elucidate the implications of these scenarios. We also outline how erroneous conclusions about the nature of tumour growth may be arrived at by looking selectively at biological data in isolation, and how this might be circumvented.

In silico disease model: from simple networks to complex diseases

In this chapter we consider in silico modeling of diseases starting from some simple to some complex (and mathematical) concepts. Examples and applications of in silico modeling for some important categories of diseases (such as for cancers, infectious diseases, and neuronal diseases) are also given.

Mathematical Modeling of Oncolytic Virotherapy

Mathematical modeling in biology has a long history as it allows the analysis and simulation of complex dynamic biological systems at little cost. A mathematical model trained on experimental or clinical data can be used to generate and evaluate hypotheses, to ask "what if" questions, and to perform in silico experiments to guide future experimentation and validation. Such models may help identify and provide insights into the mechanisms that drive changes in dynamic systems. While a mathematical model may never replace actual experiments, it can synergize with experiments to save time and resources by identifying experimental conditions that are unlikely to yield favorable outcomes, and by using optimization principles to identify experiments that are most likely to be successful. Over the past decade, numerous models have also been developed for oncolytic virotherapy, ranging from merely theoretic frameworks to fully integrated studies that utilize experimental data to generate actionable hypotheses. Here we describe how to develop such models for specific oncolytic virotherapy experimental setups, and which questions can and cannot be answered using integrated mathematical oncology.

Genetic Hallmarks and Heterogeneity of Glioblastoma in the Single-Cell Omics Era

Glioblastoma multiforme is the most common and aggressive malignant primary brain tumor. As implied by its name, the disease displays impressive intrinsic heterogeneity. Among other complications, inter- and intratumoral diversity hamper glioblastoma research and therapy, typically leaving patients with little hope for long-term survival. Extensive genetic analyses, including omics, characterize several recurrent mutations. However, confounding factors mask crucial aspects of the pathology to conventional bulk approaches. In recent years, single-cell omics have made their first appearance in cancer research, and the methodology is about to reach its full potential for glioblastoma too. Here, recent glioblastoma single-cell omics investigations are reviewed, and most promising routes toward less grim prognoses and more efficient therapeutics are discussed.

Model genotype-phenotype mappings and the algorithmic structure of evolution

Cancers are complex dynamic systems that undergo evolution and selection. Personalized medicine approaches in the clinic increasingly rely on predictions of tumour response to one or more therapies; these predictions are complicated by the inevitable evolution of the tumour. Despite enormous amounts of data on the mutational status of cancers and numerous therapies developed in recent decades to target these mutations, many of these treatments fail after a time due to the development of resistance in the tumour. The emergence of these resistant phenotypes is not easily predicted from genomic data, since the relationship between genotypes and phenotypes, termed the genotype-phenotype (GP) mapping, is neither injective nor functional. We present a review of models of this mapping within a generalized evolutionary framework that takes into account the relation between genotype, phenotype, environment and fitness. Different modelling approaches are described and compared, and many evolutionary results are shown to be conserved across studies despite using different underlying model systems. In addition, several areas for future work that remain understudied are identified, including plasticity and bet-hedging. The GP-mapping provides a pathway for understanding the potential routes of evolution taken by cancers, which will be necessary knowledge for improving personalized therapies.

Mathematical Models of Cancer Cell Plasticity

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

Modeling heterogeneous tumor growth dynamics and cell-cell interactions at single-cell and cell-population resolution

Cancer is a complex, dynamic disease that despite recent advances remains mostly incurable. Inter- and intratumoral heterogeneity are generally considered major drivers of therapy resistance, metastasis, and treatment failure. Recent advances in high-throughput experimentation have produced a wealth of data on tumor heterogeneity and researchers are increasingly turning to mathematical modeling to aid in the interpretation of these complex datasets. In this mini-review, we discuss three important classes of approaches for modeling cellular dynamics within heterogeneous tumors: agent-based models, population dynamics, and multiscale models. An important new focus, for which we provide an example, is the role of intratumoral cell-cell interactions.

Toward Personalized Computer Simulation of Breast Cancer Treatment: A Multiscale Pharmacokinetic and Pharmacodynamic Model Informed by Multitype Patient Data

The usefulness of mechanistic models to disentangle complex multiscale cancer processes, such as treatment response, has been widely acknowledged. However, a major barrier for multiscale models to predict treatment outcomes in individual patients lies in their initialization and parametrization, which needs to reflect individual cancer characteristics accurately. In this study, we use multitype measurements acquired routinely on a single breast tumor, including histopathology, MRI, and molecular profiling, to personalize parts of a complex multiscale model of breast cancer treated with chemotherapeutic and antiangiogenic agents. The model accounts for drug pharmacokinetics and pharmacodynamics. We developed an open-source computer program that simulates cross-sections of tumors under 12-week therapy regimens and used it to individually reproduce and elucidate treatment outcomes of 4 patients. Two of the tumors did not respond to therapy, and model simulations were used to suggest alternative regimens with improved outcomes dependent on the tumor's individual characteristics. It was determined that more frequent and lower doses of chemotherapy reduce tumor burden in a low proliferative tumor while lower doses of antiangiogenic agents improve drug penetration in a poorly perfused tumor. Furthermore, using this model, we were able to correctly predict the outcome in another patient after 12 weeks of treatment. In summary, our model bridges multitype clinical data to shed light on individual treatment outcomes. SIGNIFICANCE: Mathematical modeling is used to validate possible mechanisms of tumor growth, resistance, and treatment outcome.

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

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

Integrating Mathematical Modeling into the Roadmap for Personalized Adaptive Radiation Therapy

In current radiation oncology practice, treatment protocols are prescribed based on the average outcomes of large clinical trials, with limited personalization and without adaptations of dose or dose fractionation to individual patients based on their individual clinical responses. Predicting tumor responses to radiation and comparing predictions against observed responses offers an opportunity for novel treatment evaluation. These analyses can lead to protocol adaptation aimed at the improvement of patient outcomes with better therapeutic ratios. We foresee the integration of mathematical models into radiation oncology to simulate individual patient tumor growth and predict treatment response as dynamic biomarkers for personalized adaptive radiation therapy (RT).

Spatially constrained tumour growth affects the patterns of clonal selection and neutral drift in cancer genomic data

Quantification of the effect of spatial tumour sampling on the patterns of mutations detected in next-generation sequencing data is largely lacking. Here we use a spatial stochastic cellular automaton model of tumour growth that accounts for somatic mutations, selection, drift and spatial constraints, to simulate multi-region sequencing data derived from spatial sampling of a neoplasm. We show that the spatial structure of a solid cancer has a major impact on the detection of clonal selection and genetic drift from both bulk and single-cell sequencing data. Our results indicate that spatial constrains can introduce significant sampling biases when performing multi-region bulk sampling and that such bias becomes a major confounding factor for the measurement of the evolutionary dynamics of human tumours. We also propose a statistical inference framework that incorporates spatial effects within a growing tumour and so represents a further step forwards in the inference of evolutionary dynamics from genomic data. Our analysis shows that measuring cancer evolution using next-generation sequencing while accounting for the numerous confounding factors remains challenging. However, mechanistic model-based approaches have the potential to capture the sources of noise and better interpret the data.

The 2019 mathematical oncology roadmap

Whether the nom de guerre is Mathematical Oncology, Computational or Systems Biology, Theoretical Biology, Evolutionary Oncology, Bioinformatics, or simply Basic Science, there is no denying that mathematics continues to play an increasingly prominent role in cancer research. Mathematical Oncology-defined here simply as the use of mathematics in cancer research-complements and overlaps with a number of other fields that rely on mathematics as a core methodology. As a result, Mathematical Oncology has a broad scope, ranging from theoretical studies to clinical trials designed with mathematical models. This Roadmap differentiates Mathematical Oncology from related fields and demonstrates specific areas of focus within this unique field of research. The dominant theme of this Roadmap is the personalization of medicine through mathematics, modelling, and simulation. This is achieved through the use of patient-specific clinical data to: develop individualized screening strategies to detect cancer earlier; make predictions of response to therapy; design adaptive, patient-specific treatment plans to overcome therapy resistance; and establish domain-specific standards to share model predictions and to make models and simulations reproducible. The cover art for this Roadmap was chosen as an apt metaphor for the beautiful, strange, and evolving relationship between mathematics and cancer.

Systems biology: perspectives on multiscale modeling in research on endocrine-related cancers

Drawing on concepts from experimental biology, computer science, informatics, mathematics and statistics, systems biologists integrate data across diverse platforms and scales of time and space to create computational and mathematical models of the integrative, holistic functions of living systems. Endocrine-related cancers are well suited to study from a systems perspective because of the signaling complexities arising from the roles of growth factors, hormones and their receptors as critical regulators of cancer cell biology and from the interactions among cancer cells, normal cells and signaling molecules in the tumor microenvironment. Moreover, growth factors, hormones and their receptors are often effective targets for therapeutic intervention, such as estrogen biosynthesis, estrogen receptors or HER2 in breast cancer and androgen receptors in prostate cancer. Given the complexity underlying the molecular control networks in these cancers, a simple, intuitive understanding of how endocrine-related cancers respond to therapeutic protocols has proved incomplete and unsatisfactory. Systems biology offers an alternative paradigm for understanding these cancers and their treatment. To correctly interpret the results of systems-based studies requires some knowledge of how in silico models are built, and how they are used to describe a system and to predict the effects of perturbations on system function. In this review, we provide a general perspective on the field of cancer systems biology, and we explore some of the advantages, limitations and pitfalls associated with using predictive multiscale modeling to study endocrine-related cancers.

Leveraging Mathematical Modeling to Quantify Pharmacokinetic and Pharmacodynamic Pathways: Equivalent Dose Metric

Treatment response assays are often summarized by sigmoidal functions comparing cell survival at a single timepoint to applied drug concentration. This approach has a limited biophysical basis, thereby reducing the biological insight gained from such analysis. In particular, drug pharmacokinetic and pharmacodynamic (PK/PD) properties are overlooked in developing treatment response assays, and the accompanying summary statistics conflate these processes. Here, we utilize mathematical modeling to decouple and quantify PK/PD pathways. We experimentally modulate specific pathways with small molecule inhibitors and filter the results with mechanistic mathematical models to obtain quantitative measures of those pathways. Specifically, we investigate the response of cells to time-varying doxorubicin treatments, modulating doxorubicin pharmacology with small molecules that inhibit doxorubicin efflux from cells and DNA repair pathways. We highlight the practical utility of this approach through proposal of the “equivalent dose metric.” This metric, derived from a mechanistic PK/PD model, provides a biophysically-based measure of drug effect. We define equivalent dose as the functional concentration of drug that is bound to the nucleus following therapy. This metric can be used to quantify drivers of treatment response and potentially guide dosing of combination therapies. We leverage the equivalent dose metric to quantify the specific intracellular effects of these small molecule inhibitors using population-scale measurements, and to compare treatment response in cell lines differing in expression of drug efflux pumps. More generally, this approach can be leveraged to quantify the effects of various pharmaceutical and biologic perturbations on treatment response.

An Integrative Systems Biology and Experimental Approach Identifies Convergence of Epithelial Plasticity, Metabolism, and Autophagy to Promote Chemoresistance

The evolution of therapeutic resistance is a major cause of death for cancer patients. The development of therapy resistance is shaped by the ecological dynamics within the tumor microenvironment and the selective pressure of the host immune system. These selective forces often lead to evolutionary convergence on pathways or hallmarks that drive progression. Thus, a deeper understanding of the evolutionary convergences that occur could reveal vulnerabilities to treat therapy-resistant cancer. To this end, we combined phylogenetic clustering, systems biology analyses, and molecular experimentation to identify convergences in gene expression data onto common signaling pathways. We applied these methods to derive new insights about the networks at play during transforming growth factor-β (TGF-β)-mediated epithelial–mesenchymal transition in lung cancer. Phylogenetic analyses of gene expression data from TGF-β-treated cells revealed convergence of cells toward amine metabolic pathways and autophagy during TGF-β treatment. Knockdown of the autophagy regulatory, ATG16L1, re-sensitized lung cancer cells to cancer therapies following TGF-β-induced resistance, implicating autophagy as a TGF-β-mediated chemoresistance mechanism. In addition, high ATG16L expression was found to be a poor prognostic marker in multiple cancer types. These analyses reveal the usefulness of combining evolutionary and systems biology methods with experimental validation to illuminate new therapeutic vulnerabilities for cancer.

Prediction of Optimal Drug Schedules for Controlling Autophagy

The effects of molecularly targeted drug perturbations on cellular activities and fates are difficult to predict using intuition alone because of the complex behaviors of cellular regulatory networks. An approach to overcoming this problem is to develop mathematical models for predicting drug effects. Such an approach beckons for co-development of computational methods for extracting insights useful for guiding therapy selection and optimizing drug scheduling. Here, we present and evaluate a generalizable strategy for identifying drug dosing schedules that minimize the amount of drug needed to achieve sustained suppression or elevation of an important cellular activity/process, the recycling of cytoplasmic contents through (macro)autophagy. Therapeutic targeting of autophagy is currently being evaluated in diverse clinical trials but without the benefit of a control engineering perspective. Using a nonlinear ordinary differential equation (ODE) model that accounts for activating and inhibiting influences among protein and lipid kinases that regulate autophagy (MTORC1, ULK1, AMPK and VPS34) and methods guaranteed to find locally optimal control strategies, we find optimal drug dosing schedules (open-loop controllers) for each of six classes of drugs and drug pairs. Our approach is generalizable to designing monotherapy and multi therapy drug schedules that affect different cell signaling networks of interest.

Mechanism-Based Modeling of Tumor Growth and Treatment Response Constrained by Multiparametric Imaging Data

Multiparametric imaging is a critical tool in the noninvasive study and assessment of cancer. Imaging methods have evolved over the past several decades to provide quantitative measures of tumor and healthy tissue characteristics related to, for example, cell number, blood volume fraction, blood flow, hypoxia, and metabolism. Mechanistic models of tumor growth also have matured to a point where the incorporation of patient-specific measures could provide clinically relevant predictions of tumor growth and response. In this review, we identify and discuss approaches that use multiparametric imaging data, including diffusion-weighted magnetic resonance imaging, dynamic contrast-enhanced magnetic resonance imaging, diffusion tensor imaging, contrast-enhanced computed tomography, [18F]fluorodeoxyglucose positron emission tomography, and [18F]fluoromisonidazole positron emission tomography to initialize and calibrate mechanistic models of tumor growth and response. We focus the discussion on brain and breast cancers; however, we also identify three emerging areas of application in kidney, pancreatic, and lung cancers. We conclude with a discussion of the future directions for incorporating multiparametric imaging data and mechanistic modeling into clinical decision making for patients with cancer.

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.

Systems Biology, Systems Medicine, Systems Pharmacology: The What and The Why

Systems biology is today such a widespread discipline that it becomes difficult to propose a clear definition of what it really is. For some, it remains restricted to the genomic field. For many, it designates the integrated approach or the corpus of computational methods employed to handle the vast amount of biological or medical data and investigate the complexity of the living. Although defining systems biology might be difficult, on the other hand its purpose is clear: systems biology, with its emerging subfields systems medicine and systems pharmacology, clearly aims at making sense of complex observations/experimental and clinical datasets to improve our understanding of diseases and their treatments without putting aside the context in which they appear and develop. In this short review, we aim to specifically focus on these new subfields with the new theoretical tools and approaches that were developed in the context of cancer. Systems pharmacology and medicine now give hope for major improvements in cancer therapy, making personalized medicine closer to reality. As we will see, the current challenge is to be able to improve the clinical practice according to the paradigm shift of systems sciences.

A Bayesian Sequential Learning Framework to Parameterise Continuum Models of Melanoma Invasion into Human Skin

We present a novel framework to parameterise a mathematical model of cell invasion that describes how a population of melanoma cells invades into human skin tissue. Using simple experimental data extracted from complex experimental images, we estimate three model parameters: (i) the melanoma cell proliferation rate, [Formula: see text]; (ii) the melanoma cell diffusivity, D; and (iii) [Formula: see text], a constant that determines the rate that melanoma cells degrade the skin tissue. The Bayesian sequential learning framework involves a sequence of increasingly sophisticated experimental data from: (i) a spatially uniform cell proliferation assay; (ii) a two-dimensional circular barrier assay; and (iii) a three-dimensional invasion assay. The Bayesian sequential learning approach leads to well-defined parameter estimates. In contrast, taking a naive approach that attempts to estimate all parameters from a single set of images from the same experiment fails to produce meaningful results. Overall, our approach to inference is simple-to-implement, computationally efficient, and well suited for many cell biology phenomena that can be described by low-dimensional continuum models using ordinary differential equations and partial differential equations. We anticipate that this Bayesian sequential learning framework will be relevant in other biological contexts where it is challenging to extract detailed, quantitative biological measurements from experimental images and so we must rely on using relatively simple measurements from complex images.

Modeling three-dimensional invasive solid tumor growth in heterogeneous microenvironment under chemotherapy

A systematic understanding of the evolution and growth dynamics of invasive solid tumors in response to different chemotherapy strategies is crucial for the development of individually optimized oncotherapy. Here, we develop a hybrid three-dimensional (3D) computational model that integrates pharmacokinetic model, continuum diffusion-reaction model and discrete cell automaton model to investigate 3D invasive solid tumor growth in heterogeneous microenvironment under chemotherapy. Specifically, we consider the effects of heterogeneous environment on drug diffusion, tumor growth, invasion and the drug-tumor interaction on individual cell level. We employ the hybrid model to investigate the evolution and growth dynamics of avascular invasive solid tumors under different chemotherapy strategies. Our simulations indicate that constant dosing is generally more effective in suppressing primary tumor growth than periodic dosing, due to the resulting continuous high drug concentration. In highly heterogeneous microenvironment, the malignancy of the tumor is significantly enhanced, leading to inefficiency of chemotherapies. The effects of geometrically-confined microenvironment and non-uniform drug dosing are also investigated. Our computational model, when supplemented with sufficient clinical data, could eventually lead to the development of efficient in silico tools for prognosis and treatment strategy optimization.

Longitudinal Liquid Biopsy and Mathematical Modeling of Clonal Evolution Forecast Time to Treatment Failure in the PROSPECT-C Phase II Colorectal Cancer Clinical Trial

Sequential profiling of plasma cell-free DNA (cfDNA) holds immense promise for early detection of patient progression. However, how to exploit the predictive power of cfDNA as a liquid biopsy in the clinic remains unclear. RAS pathway aberrations can be tracked in cfDNA to monitor resistance to anti-EGFR monoclonal antibodies in patients with metastatic colorectal cancer. In this prospective phase II clinical trial of single-agent cetuximab in RAS wild-type patients, we combine genomic profiling of serial cfDNA and matched sequential tissue biopsies with imaging and mathematical modeling of cancer evolution. We show that a significant proportion of patients defined as RAS wild-type based on diagnostic tissue analysis harbor aberrations in the RAS pathway in pretreatment cfDNA and, in fact, do not benefit from EGFR inhibition. We demonstrate that primary and acquired resistance to cetuximab are often of polyclonal nature, and these dynamics can be observed in tissue and plasma. Furthermore, evolutionary modeling combined with frequent serial sampling of cfDNA allows prediction of the expected time to treatment failure in individual patients. This study demonstrates how integrating frequently sampled longitudinal liquid biopsies with a mathematical framework of tumor evolution allows individualized quantitative forecasting of progression, providing novel opportunities for adaptive personalized therapies.Significance: Liquid biopsies capture spatial and temporal heterogeneity underpinning resistance to anti-EGFR monoclonal antibodies in colorectal cancer. Dense serial sampling is needed to predict the time to treatment failure and generate a window of opportunity for intervention. Cancer Discov; 8(10); 1270-85. ©2018 AACR. See related commentary by Siravegna and Corcoran, p. 1213 This article is highlighted in the In This Issue feature, p. 1195.

Nonlinear studies of tumor morphological stability using a two-fluid flow model

We consider the nonlinear dynamics of an avascular tumor at the tissue scale using a two-fluid flow Stokes model, where the viscosity of the tumor and host microenvironment may be different. The viscosities reflect the combined properties of cell and extracellular matrix mixtures. We perform a linear morphological stability analysis of the tumors, and we investigate the role of nonlinearity using boundary-integral simulations in two dimensions. The tumor is non-necrotic, although cell death may occur through apoptosis. We demonstrate that tumor evolution is regulated by a reduced set of nondimensional parameters that characterize apoptosis, cell-cell/cell-extracellular matrix adhesion, vascularization and the ratio of tumor and host viscosities. A novel reformulation of the equations enables the use of standard boundary integral techniques to solve the equations numerically. Nonlinear simulation results are consistent with linear predictions for nearly circular tumors. As perturbations develop and grow, the linear and nonlinear results deviate and linear theory tends to underpredict the growth of perturbations. Simulations reveal two basic types of tumor shapes, depending on the viscosities of the tumor and microenvironment. When the tumor is more viscous than its environment, the tumors tend to develop invasive fingers and a branched-like structure. As the relative ratio of the tumor and host viscosities decreases, the tumors tend to grow with a more compact shape and develop complex invaginations of healthy regions that may become encapsulated in the tumor interior. Although our model utilizes a simplified description of the tumor and host biomechanics, our results are consistent with experiments in a variety of tumor types that suggest that there is a positive correlation between tumor stiffness and tumor aggressiveness.

SpheroidSim-Preliminary evaluation of a new computational tool to predict the influence of cell cycle time and phase fraction on spheroid growth

BACKGROUND

There is a relative paucity of research that integrates materials science and bioengineering with computational simulations to decipher the intricate processes promoting cancer progression. Therefore, a first-generation computational model, SpheroidSim, was developed that includes a biological data set derived from a bioengineered spheroid model to obtain a quantitative description of cell kinetics.

RESULTS

SpheroidSim is a 3D agent-based model simulating the growth of multicellular cancer spheroids. Cell cycle time and phases mathematically motivated the population growth. SpheroidSim simulated the growth dynamics of multiple spheroids by individually defining a collection of specific phenotypic traits and characteristics for each cell. Experimental data derived from a hydrogel-based spheroid model were fit to the predictions providing insight into the influence of cell cycle time (CCT) and cell phase fraction (CPF) on the cell population. A comparison of the number of active cells predicted for each analysis showed that the value and method used to define CCT had a greater effect on the predicted cell population than CPF. The model predictions were similar to the experimental results for the number of cells, with the predicted total number of cells varying by 8% and 12%, respectively, compared to the experimental data.

CONCLUSIONS

SpheroidSim is a first step in developing a biologically based predictive tool capable of revealing fundamental elements in cancer cell physiology. This computational model may be applied to study the effect of the microenvironment on spheroid growth and other cancer cell types that demonstrate a similar multicellular clustering behavior as the population develops. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:1335-1343, 2018.

The biology and mathematical modelling of glioma invasion: a review

Adult gliomas are aggressive brain tumours associated with low patient survival rates and limited life expectancy. The most important hallmark of this type of tumour is its invasive behaviour, characterized by a markedly phenotypic plasticity, infiltrative tumour morphologies and the ability of malignant progression from low- to high-grade tumour types. Indeed, the widespread infiltration of healthy brain tissue by glioma cells is largely responsible for poor prognosis and the difficulty of finding curative therapies. Meanwhile, mathematical models have been established to analyse potential mechanisms of glioma invasion. In this review, we start with a brief introduction to current biological knowledge about glioma invasion, and then critically review and highlight future challenges for mathematical models of glioma invasion.

Effect of cellular de-differentiation on the dynamics and evolution of tissue and tumor cells in mathematical models with feedback regulation

Tissues are maintained by adult stem cells that self-renew and also differentiate into functioning tissue cells. Homeostasis is achieved by a set of complex mechanisms that involve regulatory feedback loops. Similarly, tumors are believed to be maintained by a minority population of cancer stem cells, while the bulk of the tumor is made up of more differentiated cells, and there is indication that some of the feedback loops that operate in tissues continue to be functional in tumors. Mathematical models of such tissue hierarchies, including feedback loops, have been analyzed in a variety of different contexts. Apart from stem cells giving rise to differentiated cells, it has also been observed that more differentiated cells can de-differentiate into stem cells, both in healthy tissue and tumors, aspects of which have also been investigated mathematically. This paper analyses the effect of de-differentiation on the basic and evolutionary dynamics of cells in the context of tissue hierarchy models that include negative feedback regulation of the cell populations. The models predict that in the presence of de-differentiation, the fixation probability of a neutral mutant is lower than in its absence. Therefore, if de-differentiation occurs, a mutant with identical parameters compared to the wild-type cell population behaves like a disadvantageous mutant. Similarly, the process of de-differentiation is found to lower the fixation probability of an advantageous mutant. These results indicate that the presence of de-differentiation can lower the rates of tumor initiation and progression in the context of the models considered here.

Precision Medicine with Imprecise Therapy: Computational Modeling for Chemotherapy in Breast Cancer

Medical oncology is in need of a mathematical modeling toolkit that can leverage clinically-available measurements to optimize treatment selection and schedules for patients. Just as the therapeutic choice has been optimized to match tumor genetics, the delivery of those therapeutics should be optimized based on patient-specific pharmacokinetic/pharmacodynamic properties. Under the current approach to treatment response planning and assessment, there does not exist an efficient method to consolidate biomarker changes into a holistic understanding of treatment response. While the majority of research on chemotherapies focus on cellular and genetic mechanisms of resistance, there are numerous patient-specific and tumor-specific measures that contribute to treatment response. New approaches that consolidate multimodal information into actionable data are needed. Mathematical modeling offers a solution to this problem. In this perspective, we first focus on the particular case of breast cancer to highlight how mathematical models have shaped the current approaches to treatment. Then we compare chemotherapy to radiation therapy. Finally, we identify opportunities to improve chemotherapy treatments using the model of radiation therapy. We posit that mathematical models can improve the application of anticancer therapeutics in the era of precision medicine. By highlighting a number of historical examples of the contributions of mathematical models to cancer therapy, we hope that this contribution serves to engage investigators who may not have previously considered how mathematical modeling can provide real insights into breast cancer therapy.