- Open Access
Mathematical models for translational and clinical oncology
Journal of Clinical Bioinformaticsvolume 3, Article number: 23 (2013)
In the context of translational and clinical oncology, mathematical models can provide novel insights into tumor-related processes and can support clinical oncologists in the design of the treatment regime, dosage, schedule, toxicity and drug-sensitivity. In this review we present an overview of mathematical models in this field beginning with carcinogenesis and proceeding to the different cancer treatments. By doing so we intended to highlight recent developments and emphasize the power of such theoretical work.
We first highlight mathematical models for translational oncology comprising epidemiologic and statistical models, mechanistic models for carcinogenesis and tumor growth, as well as evolutionary dynamics models which can help to describe and overcome a major problem in the clinic: therapy resistance. Next we review models for clinical oncology with a special emphasis on therapy including chemotherapy, targeted therapy, radiotherapy, immunotherapy and interaction of cancer cells with the immune system.
As evident from the published studies, mathematical modeling and computational simulation provided valuable insights into the molecular mechanisms of cancer, and can help to improve diagnosis and prognosis of the disease, and pinpoint novel therapeutic targets.
Cancer is still one of the leading causes of death in the world and major efforts have been undertaken to improve diagnosis and therapy of common cancer types. Recently developed technologies (i.e. next generation sequencing) give us unprecedented opportunities to study individual cancer samples at the molecular level and to identify genomic variants and rearrangements . This information will build the basis for the stratification of patients, and for personalized or precision medicine. The increasing complexity of the generated data utilizing various high-throughput technologies for characterizing the genome, epigenome, transcriptome, proteome, metabolome, and interactome pose considerable challenges and therefore plethora of bioinformatics methods and tools for the analysis have been developed . However, the real value of the disparate datasets can be truly exploited only if the data is integrated and will then enable one to comprehensively study molecular mechanisms of cancer cells.
One possibility for data integration is the use of mathematical models. Modeling has been successfully applied in physiology for many decades but only recently the quality and the quantity of biomolecular data became available for the development of causative and predictive models. Due to their importance in cancer mathematical models have also been in the focus of theoretical investigators. For example application of theoretical techniques and the postulation of the "two hit" hypothesis in the early 70s led to the identification of tumor-suppressor genes . Later, in a landmark paper it was shown that cancer results from evolutionary processes occurring within the body .
In the context of translational (i.e. from bench to bedside, or in other words: transforming scientific discoveries arising from laboratory to clinical applications) and clinical oncology, mathematical models can provide novel insights into tumor growth and progression, into tumor-related processes such as angiogenesis, the immune response, and the interaction with the tumor microenvironment, and into the development of drug resistance. Furthermore, modeling can support the clinical oncologists in the design of the treatment regime, dosage, schedule, toxicity and drug-sensitivity. Common treatments against the different types of cancers include surgery, radiation therapy, chemotherapy, targeted therapy or combinations of those to limit the progression of malignant disease, eradicate tumor cells and prolong survival. The information gained from mathematical models can also help in the development and efficacy of clinical trials and treatment protocols, and can accelerate the progress of clinical research in fighting cancer.
To the best of our knowledge there is currently no review study on mathematical models focusing on translational and clinical oncology applications except for a similar attempt made by Swierniak et al.  few years ago. We therefore initiated this work to provide an overview of the field and stimulate the discussion and the development of novel models. While mechanistic models have proven extremely valuable and provided novel insights, there are not considered here and we refer the readers to recent reviews [6–8]. Given the wealth of published studies using mathematical models in cancer, we by no means intended to provide a comprehensive picture. Rather, we selected several topics we believe are of relevance for the readers. Wherever possible, we refer to additional reviews in order to guide interested researchers.
We first highlight mathematical models for translational oncology comprising epidemiologic and statistical models, mechanistic models for carcinogenesis and tumor growth, as well as evolutionary dynamics models , which can help to describe and overcome a major problem in the clinic: therapy resistance. Next we review models for clinical oncology. It should be noted that a survey of application of modeling results in clinics was beyond the scope of this review. Rather, we provide an overview of the models with a special emphasis on therapy including chemotherapy, targeted therapy, radiotherapy, immunotherapy and interaction of cancer cells with the immune system. Table 1 shows the specific categories and the publications used in this work.
Mathematical models for translational oncology
Carcinogenesis and tumor-growth models
Early models that aimed to explain the dynamics of cancer progression were based on experimental and epidemiological data, which indicated that the cancer incidence is often rapidly increased with age and simple patterns could be observed at the population level. Fisher and Hollomon  presented a multicellular model in which mutations occur in different cells within the same cell population and only the combination of all mutations leads to cancer development. As an alternative to this theory, Nordling  suggested that mutations must occur sequentially in the same cell for transformation into cancer cell.
Most mathematical models of cancer progression descend from Armitage and Doll’s  multistage theory, which include major concepts for how to think about incidence, carcinogenesis, and progression. The theory states that carcinogenesis progresses through series of genomic alterations in a single cell and the age-specific incidence of cancers is predicted to increase with a power of age that is one less than the number of alterations. Two other studies, using data comparing inherited and noninherited cases in colon cancer  and retinoblastoma , provided additional empirical evidence for the multistage theory. Knudson used a statistical analysis of the incidence of retinoblastoma in children to explain the role of tumor suppressor genes in sporadic and inherited cancers. This work was later extended to a two-stage stochastic model for the process of cancer initiation and progression , which lead to important subsequent work [15, 16] that helped with characterization of other suppressor genes such as APC in colon cancer and TP53, which is mutated in several human tumors.
Even though these models provided accurate descriptions of cancer incidence data, they were unable to relate the data with the functional changes associated with tumor progression. Since then the understanding of the molecular mechanisms underlying tumor initiation and progression has improved  and mechanistic models that use biological knowledge and biophysical laws to quantify and predict cancer progression were developed.
The growth and development of solid tumors occurs in two stages – avascular and vascular. The early spatio-temporal models [18, 19] of avascular tumor growth describe the interactions between tumor cell population and nutrients and calculate the nutrient concentrations as a function of tumor spheroid radius that is changing due to the rate of cell proliferation. Significant progress was made with the development of new models [20, 21] that introduced the interrelated concepts of cell movement and pressure.
Since tumor induced angiogenesis i.e. the growth of a network of blood vessels, is a crucial component of solid tumor growth, the basic models have been expanded to account for tumor growth during angiogenesis and the increase of tumor availability associated with the expanding vasculature. In order to make a transition from avascular to vascular growth, tumors may secrete diffusible substances called tumor angiogenic factors (TAF). The earliest continuum models of tumor angiogenesis  describe the growth of a capillary network in terms of capillary tip densities and capillary sprout densities in response to TAF. The mathematical models in angiogenesis have mostly focused on describing endothelial cell migration and proliferation through the extracellular matrix [23–27]. A comprehensive overview of models in this area can be found in . Mathematical modeling of blood flow in tumor-induced capillary networks has been described in more recent studies [29–31].
The concept that the successful formation of a tumor depends on vascularization has resulted in developing cancer therapies designed to inhibit the tumor vasculature in order to deprive the tumor from oxygen and nutrients. Several models have focused on exploring the efficiency of such antiangiogenic treatments [32–34]. Using methods of optimal control theory to analyze drug dosing and treatment strategies these studies showed that the combination with other forms of therapy would be beneficial.
Clonal evolution models and therapy resistance
An important conceptual breakthrough in understanding cancer lies in Darwinian and ecological theories: cancer progression is an evolutionary process that results from accumulation of genetic and epigenetic variations in somatic cells [35, 36]. Experimental evidences and recent advances in genetic sequencing technologies – that allowed identification of the genetic alterations in a cancer cell - have revealed the complexity and heterogeneity of cancer progression and have stimulated the use of evolutionary-based approaches in the study of cancer.
Several methods of population dynamics and evolutionary game theory were applied to account for the elementary principles of evolution that lead to tumor initiation and progression. In the earliest models, mutations accumulate in a population of constant or variable size, and they consider only one or two mutations [37, 38]. Newer models are now being used to investigate how the sequence and timing of mutations and the environmental conditions influence tumor progression [39, 40]. An in-depth review of models that describe the evolutionary dynamics of cancer can be found in .
Several studies have focused on the waiting time to cancer development, which may be defined as the time from the first presence of neoplasm, until a critical number of hits (driver mutations) are accumulated and initiate the growth of carcinoma. Beerenwinkel et al.  developed one of the first models that was based on genomic studies of colorectal cancer patients. They related the waiting time to the population size, mutation rate, and the advantage of the driver mutations and showed that selective advantage of mutations has the largest effect on the evolutionary dynamics of tumorigenesis. In a similar manner, Bozic et al. , by fitting their model to glioblastoma and pancreatic cancer data, estimated that driver mutations give an average fitness advantage of 0.4%.
Another characteristic of evolutionary processes is the influence of the local cellular environment on the tumor progression. The progress of tumor is characterized not only by the genetic and epigenetic changes accumulating in the cells, but also by the dynamic interactions between cells within the tumor and between the cells and the constantly changing microenvironment. The microenvironment provides a selective fitness landscape that includes competing for limited resources and active intracellular (initiation of cell proliferation and cell death) and extracellular control mechanisms (the immune system) that aim to restore homeostasis.
There are several studies that utilize mathematical modeling to predict and quantify the interactions of the tumor cells with the surrounding environment during tumor progression [44, 45]. Gatenby et al.  developed a model of carcinogenesis according to which the tumor cells have to overcome six microenvironment barriers that appear as tumor cells proliferate. They proposed that the nature and sequence of the alterations during carcinogenesis are determined by the specific microenvironmental properties that prevent proliferation within changing adaptive landscapes.
An important clinical problem in cancer research that can be analyzed using modeling techniques is the development of resistance to targeted therapies. Resistance to drugs may develop as a consequence of genetic events such as point mutations or gene amplifications. The emergence of resistance to therapy as a result from a single mutational effect has been first introduced in a model of Coldman and Goldie . More recent studies have also used point mutations to explain the evolutionary dynamics of drug resistant cancer cells [48–51]. Other models studied gene amplification as one of the mechanisms that has a strong influence on the evolution of drug resistance [52–56].
Foo et al.  designed a methodology that can be used to investigate optimal drug dosing schedules to avoid resistance conferred by one (epi)genetic mutation. In a recent study, Diaz et al.  showed that tumors became resistant to anti-EGFR antibodies as a result of emergence of resistance mutations in KRAS and other genes that were present in clonal subpopulation within the tumors before the initiation of the treatment. Analyzing data from 20 melanoma patients who received targeted therapy, Bozic et al.  found that simultaneous administration of two drugs is much more effective than sequential therapy. The improved understanding of the evolutionary dynamics of cancer provided by these models can have practical implications in the design and administration of new cancer therapeutics.
In summary, using the overwhelming amount of generated knowledge in tumor biology, mathematical modelers have succeeded in formalizing this knowledge and make it usable for simulations. Moreover, the published models represent a unique basis for testing novel hypotheses which are otherwise difficult or even impossible to test. For example, it is very difficult to obtain samples from early cancer stages or longitudinal samples in order to study the development of tumor heterogeneity. The models presented above enable researchers to address questions which were previously not possible and by using iterative cycles of simulations and experimentation ultimately lead to novel knowledge (Figure 1). Moreover, the maturity of the tools and the availability of data in public databases are additionally supporting the translation of this knowledge into clinical practice.
Mathematical models for clinical oncology
Chemotherapy and targeted therapy
Chemotherapy is widely used therapy against cancer. Proliferating cells undergo different phases during the cell cycle including DNA replication and cell division and different chemotherapeutic compounds are affecting cells in different phases. The greatest challenge associated with chemotherapy is that not only cancer cells but also normal dividing cells are affected. In contrast, targeted cancer therapy  is aimed at specific molecular targets and is therefore more effective and has fewer side effects.
Fister et al.  developed a model that try to solve the problem of maximizing the effect on cancer cells but also maximizing the survival of the host cells. The mathematical model shows that if shorter periods of treatment are used it is possible to use higher doses of the drug and destroy more cancer cells without causing excess damage to the host cells. A more complex model is the cellular automaton model of Byrne et al. . It is a multi-scale model that has a vascular layer, a cellular layer and an intracellular layer. With this model it is possible to get detailed spatio-temporal information about the tumor and the healthy tissue. In general three results are possible with the model; the tumor is eliminated, the tumor continues grow, or an oscillation. It is possible to combine different treatments in one model. The model of de Pillis et al.  based on a previous model  combines different treatments and shows that if the chemotherapy is stopped a system with a undetectable tumor can return to a high tumor state.
The combination of different treatments is one possibility to eliminate the patient’s tumor. Jackson et al.  introduced two different types of cells and investigated the tumors response to different chemotherapeutic strategies. It was possible to estimate the largest size of a tumor that can be eradicated by a bolus injection. With only one cell population the results of bolus and continues drug were similar. With two populations, one of them drug resistant, the continuous infusion increased the time to cure. This indicates that it is important to tailor treatment strategies.
Another interesting aspect is the use of growth factors in the model from Panetta and Adam . They showed that the use of growth factors in their model increases the cell killing up to 20%. Based on the model of Anderson and Chaplain  McDougall  developed a model where the blood viscosity, pressure drop and mean capillary radii can be varied of the surrounding vessels of a solid tumor can be analyzed. The model shows that if there are highly interconnected vessels around the tumor there is a low drug delivery to the tumor itself. It shows that it is important to consider the vasculature around the tumor to find the optimal chemotherapy strategy.
The strategy of chemotherapy in combination with other treatments is being increasingly used. Powathil addressed this in recent publications [67–69] and showed ways to simulate and improve protocols of chemotherapy. It was demonstrated that the cytotoxic effect is dependent on many factors like timing of the drug delivery, time delay between the doses, heterogeneities of the cell cycle, the spatial distribution of the tumor and the surrounding microenvironment. It is noteworthy that these issues have been also investigated in older studies, e.g. using models of phase-specificity of chemotherapeutic drugs published in the 1990s . In this section different methods were shown to model chemotherapy and effects that can lead to a better treatment strategy.
Most targeted cancer therapies use monoclonal antibodies directed against tumor-specific surface proteins or small-molecule drugs against intracellular targets (e.g. tyrosin kinases) . Billy et al.  developed a model that simulates a treatment on the angiogenesis of tumors by gene therapy. The gene therapy is delivered by adenoviruses and influences, the antagonist of vascular endothelial growth factor, endostatin. The simulation showed that there is a critical treatment dose which is important to improve the efficacy.
TGF-β is a cytokine that has an immunosuppressive effect. In the model of Kirschner et al.  it was shown that a treatment with initial delivery of double stranded RNA into tumor cells that is cut by the enzyme Dicer into 21–23 segments known as siRNA inhibits TGF-β production and leads to a controlled oscillatory tumor behavior. Using a combination of experimental data and a mathematical model about the resistance against the monoclonal antibody panitumumab based on the Luria-Delbrück model , Diaz et al.  tested the development of mutations conferring resistance to the antibody. The simulation results suggested a combination of therapies where at least two pathways will be required. The use of ex-vivo activated alloreactive cytotoxic-T-lymphocytes (CTL) is another possibility to direct target the tumor. Kronik et al.  developed a mathematical model to investigate the effect of directly administrated CTL to glioblastomas. They showed that most sensitive parameters were the death rate of CTLs, the initial size of the tumor and the maximal growth rate.
Nanda et al.  developed a mathematical model simulating the drug imatinib mesylate that was approved in 2002 by the FDA for use in newly diagnosed cases of chronic myelogenous leukemia. The results show that a high dosing level from the beginning is optimal. Another interesting aspect of targeted therapy is the use of oncolytic viruses. Wein et al.  showed in their model that a single intratumoral injection in a solid tumor is not enough to effectively spread the virus. Also important is the suppression of the immune-mediated clearance of the virus. In the work of Mok et al.  two additional modifications are shown through mathematical modeling of herpes simplex viruses first the decreasing of the binding affinity of the virus and second the effective diffusion coefficient of the virus through degradation of the tumor extracellular matrix.
The aim of radiotherapy is to destroy the tumor cells but not the host cells. This is possible if the tumor cells are more sensitive to irradiation than the host cells. Mathematical modeling can show strategies and improve treatment protocols to obtain an optimal patient treatment. In this sense Rockne et al.  present a model to investigate the response to various schedules and dose distribution on a virtual tumor. The advantage in the mathematical simulation is that the effect of radiation can be observed continuously. The model suggests that a radiation dose on daily basis is more effective than several treatments per day.
Another important aspect is the general response of cells to radiation. Richard et al.  used a cellular automaton model to investigate these mechanisms after low doses of radiation. Enderling et al.  developed a model that simulates the recurrence after radiotherapy. In the 2D simulations it was shown that if pre-malignant cells reside in the breast post-surgery and survive radiotherapy this cells could be the reason for a recurrence.
Tumor immune-cell interaction and immunotherapy
The immune system plays an important role in tumor progression. Immune processes with different components like chemokines, cytokines or different cell types that work together are highly complex and intertwined. Mathematical modeling has already provided deeper insights and helped to get fundamental knowledge and improve patient’s treatment. For example De Boer et al.  developed a detailed model where they were able to show tumor regression and tumor growth dependent on the antigenicity of tumor-immune interaction. Tumour-infiltrating cytotoxic lymphocytes (TICLs) play an important role in tumor-immune interaction. Matzavinos et al.  developed a spatio-temporal model to investigate the interaction of TICLs and tumors. It is possible to simulate the spatio-temporal dynamics of TICLs in a solid tumor. Kirschner et al.  developed a model that includes immunotherapy with cultured immune cells that have anti-tumor reactivity and additionally IL-2. In simulations a total eradiation of the tumor was only possible with the immune therapy. In the model of de Pillis et al.  the cytolytic effectiveness of tumor specific T-cells was the most sensitive parameter. Following the simulation results the efficacy of the CD8+ T cells and the response to immunotherapy was correlating.
One therapy against superficial bladder cancer is the treatment with Bacillus Calmette-Guerin (BCG). Rentsch et al.  showed with their mathematical model that the dose of BCG and the treatment interval have a positive correlation of tumor extension. Wei  investigated this immunotherapy with a mathematical model and showed that the infection rate and the growth rate of the tumor are the most important parameters for a successful treatment. Rihan et al.  investigated the effect of adoptive cellular immunotherapy and found out that only a combination of the treatment with IL-2 can be used to clear the tumor.
In summary, major contributions for clinical oncology have been made by the modeling community. However, although many models were designed and tested for clinical applications, the use in routine setting is sparse. One way to overcome this is to develop models for very specific applications and rigorously test the performance and the predictive power. Furthermore, the use of the available knowledge should be also part of the decision process. We envision a computational decision support system which is using clinical data, molecular data, publicly available data, as well as simulation results of mathematical models to reach a decision for therapeutic strategy.
In this review we presented an overview of mathematical models for translational oncology and clinical oncology beginning with carcinogenesis and proceeding to the different cancer treatments. By doing so we intended to highlight recent developments in the field and emphasize the power of this theoretical work. As demonstrated in a number of studies, mathematical modeling and computational simulation can provide valuable insights into the molecular mechanisms of cancer, can improve diagnosis and prognosis of the disease, and pinpoint novel therapeutic targets. As can be seen in Figure 1, it is often difficult to attribute the generation of new knowledge either to the modeling or to the experimental work. Regardless the origin, the insights obtained from such cycles of experiments and modeling can improve our understanding of the complexity of cancer progression and eventually be used to stop or at least slow down the processes of tumor initiation, evolution and resistance to therapies.
Vogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz LA, Kinzler KW: Cancer genome landscapes. Science. 2013, 339: 1546-1558. 10.1126/science.1235122.
Pabinger S, Dander A, Fischer M, Snajder R, Sperk M, Efremova M, Krabichler B, Speicher MR, Zschocke J, Trajanoski Z: A survey of tools for variant analysis of next-generation genome sequencing data. Brief Bioinform. 2013
Knudson AG: Mutation and cancer: statistical study of retinoblastoma. Proc Natl Acad Sci USA. 1971, 68: 820-823. 10.1073/pnas.68.4.820.
Nowell PC: The clonal evolution of tumor cell populations. Science. 1976, 194: 23-28. 10.1126/science.959840.
Swierniak A, Kimmel M, Smieja J: Mathematical modeling as a tool for planning anticancer therapy. Eur J Pharmacol. 2009, 625: 108-121. 10.1016/j.ejphar.2009.08.041.
Lowengrub JS, Frieboes HB, Jin F, Chuang Y-L, Li X, Macklin P, Wise SM, Cristini V: Nonlinear modelling of cancer: bridging the gap between cells and tumours. Nonlinearity. 2010, 23: R1-R9. 10.1088/0951-7715/23/1/R01.
Byrne HM, Alarcon T, Owen MR, Webb SD, Maini PK: Modelling aspects of cancer dynamics: a review. Philos Transact A Math Phys Eng Sci. 2006, 364: 1563-1578. 10.1098/rsta.2006.1786.
Deisboeck TS, Wang Z, Macklin P, Cristini V: Multiscale cancer modeling. Annu Rev Biomed Eng. 2011, 13: 127-155. 10.1146/annurev-bioeng-071910-124729.
Stransky B, de Souza SJ: Modeling tumor evolutionary dynamics. Front Physiol. 2012, 3: 480-
Fisher JC, Hollomon JH: A hypothesis for the origin of cancer foci. Cancer. 1951, 4: 916-918. 10.1002/1097-0142(195109)4:5<916::AID-CNCR2820040504>3.0.CO;2-7.
Nordling CO: A new theory on cancer-inducing mechanism. Br J Cancer. 1953, 7: 68-72. 10.1038/bjc.1953.8.
Armitage P, Doll R: The age distribution of cancer and a multi-stage theory of carcinogenesis. Br J Cancer. 1954, 8: 1-12. 10.1038/bjc.1954.1.
Ashley DJ: The two "hit" and multiple "hit" theories of carcinogenesis. Br J Cancer. 1969, 23: 313-328. 10.1038/bjc.1969.41.
Moolgavkar SH, Knudson AG: Mutation and cancer: a model for human carcinogenesis. J Natl Cancer Inst. 1981, 66: 1037-1052.
Luebeck EG, Moolgavkar SH: Multistage carcinogenesis and the incidence of colorectal cancer. Proc Natl Acad Sci U S A. 2002, 99: 15095-15100. 10.1073/pnas.222118199.
Gatenby RA, Vincent TL: An evolutionary model of carcinogenesis. Cancer Res. 2003, 63: 6212-6220.
Hanahan D, Weinberg RA: The hallmarks of cancer. Cell. 2000, 100: 57-70. 10.1016/S0092-8674(00)81683-9.
Greenspan H: Models for the growth of a solid tumor by diffusion. Stud Appl Math. 1972, 52: 317-340.
McElwain DLS, Ponzo PJ: A model for the growth of a solid tumor with non-uniform oxygen consumption. Math Biosci. 1977, 35: 267-279. 10.1016/0025-5564(77)90028-1.
Greenspan HP: On the growth and stability of cell cultures and solid tumors. J Theor Biol. 1976, 56: 229-242. 10.1016/S0022-5193(76)80054-9.
Byrne HM, Chaplain MAJ: Modelling the role of cell-cell adhesion in the growth and development of carcinomas. Math Comput Model. 1996, 24: 1-17.
Balding D, McElwain DL: A mathematical model of tumour-induced capillary growth. J Theor Biol. 1985, 114: 53-73. 10.1016/S0022-5193(85)80255-1.
Byrne HM, Chaplain MA: Mathematical models for tumour angiogenesis: numerical simulations and nonlinear wave solutions. Bull Math Biol. 1995, 57: 461-486.
Chaplain MA: Mathematical modelling of angiogenesis. J Neurooncol. 2000, 50: 37-51. 10.1023/A:1006446020377.
Chaplain MA, Stuart AM: A model mechanism for the chemotactic response of endothelial cells to tumour angiogenesis factor. IMA J Math Appl Med Biol. 1993, 10: 149-168. 10.1093/imammb/10.3.149.
Levine HA, Pamuk S, Sleeman BD, Nilsen-Hamilton M: Mathematical modeling of capillary formation and development in tumor angiogenesis: penetration into the stroma. Bull Math Biol. 2001, 63: 801-863. 10.1006/bulm.2001.0240.
Anderson AR, Chaplain MA: Continuous and discrete mathematical models of tumor-induced angiogenesis. Bull Math Biol. 1998, 60: 857-899. 10.1006/bulm.1998.0042.
Mantzaris NV, Webb S, Othmer HG: Mathematical modeling of tumor-induced angiogenesis. J Math Biol. 2004, 49: 111-187.
McDougall SR, Anderson ARA, Chaplain MAJ, Sherratt JA: Mathematical modelling of flow through vascular networks: implications for tumour-induced angiogenesis and chemotherapy strategies. Bull Math Biol. 2002, 64: 673-702. 10.1006/bulm.2002.0293.
Stéphanou A, McDougall SR, Anderson ARA, Chaplain MAJ: Mathematical modelling of the influence of blood rheological properties upon adaptative tumour-induced angiogenesis. Math Comput Model. 2006, 44: 96-123. 10.1016/j.mcm.2004.07.021.
Stéphanou A, McDougall SR, Anderson ARA, Chaplain MAJ: Mathematical modelling of flow in 2D and 3D vascular networks: applications to anti-angiogenic and chemotherapeutic drug strategies. Math Comput Model. 2005, 41: 1137-1156. 10.1016/j.mcm.2005.05.008.
D’ Onofrio A, Gandolfi A: Tumour eradication by antiangiogenic therapy: analysis and extensions of the model by Hahnfeldt et al. (1999). Math Biosci. 2004, 191: 159-184. 10.1016/j.mbs.2004.06.003.
Ledzewicz U, Schaettler H: Anti-angiogenic therapy in cancer treatment as an optimal control problem. SIAM J Control Optim. 2007, 46 (3): 1052-1079. 10.1137/060665294.
Hahnfeldt P, Panigrahy D, Folkman J, Hlatky L: Tumor development under angiogenic signaling: a dynamical theory of tumor growth, treatment response, and postvascular dormancy. Cancer Res. 1999, 59: 4770-4775.
Merlo LMF, Pepper JW, Reid BJ, Maley CC: Cancer as an evolutionary and ecological process. Nat Rev Cancer. 2006, 6: 924-935. 10.1038/nrc2013.
Greaves M, Maley CC: Clonal evolution in cancer. Nature. 2012, 481: 306-313. 10.1038/nature10762.
Haeno H, Iwasa Y, Michor F: The evolution of two mutations during clonal expansion. Genetics. 2007, 177: 2209-2221. 10.1534/genetics.107.078915.
Nowak MA, Michor F, Komarova NL, Iwasa Y: Evolutionary dynamics of tumor suppressor gene inactivation. Proc Natl Acad Sci U S A. 2004, 101: 10635-10638. 10.1073/pnas.0400747101.
Gerstung M, Eriksson N, Lin J, Vogelstein B, Beerenwinkel N: The temporal order of genetic and pathway alterations in tumorigenesis. PLoS One. 2011, 6: e27136-10.1371/journal.pone.0027136.
Durrett R, Schmidt D, Schweinsberg J: A waiting time problem arising from the study of multi-stage carcinogenesis. Ann Appl Probab. 2009, 19: 676-718. 10.1214/08-AAP559.
Attolini CS-O, Michor F: Evolutionary theory of cancer. Ann N Y Acad Sci. 2009, 1168: 23-51. 10.1111/j.1749-6632.2009.04880.x.
Beerenwinkel N, Antal T, Dingli D, Traulsen A, Kinzler KW, Velculescu VE, Vogelstein B, Nowak MA: Genetic progression and the waiting time to cancer. PLoS Comput Biol. 2007, 3: e225-10.1371/journal.pcbi.0030225.
Bozic I, Antal T, Ohtsuki H, Carter H, Kim D, Chen S, Karchin R, Kinzler KW, Vogelstein B, Nowak MA: Accumulation of driver and passenger mutations during tumor progression. Proc Natl Acad Sci U S A. 2010, 107: 18545-18550. 10.1073/pnas.1010978107.
Quaranta V, Rejniak KA, Gerlee P, Anderson ARA: Invasion emerges from cancer cell adaptation to competitive microenvironments: quantitative predictions from multiscale mathematical models. Semin Cancer Biol. 2008, 18: 338-348. 10.1016/j.semcancer.2008.03.018.
Lee H-O, Silva AS, Concilio S, Li Y-S, Slifker M, Gatenby RA, Cheng JD: Evolution of tumor invasiveness: the adaptive tumor microenvironment landscape model. Cancer Res. 2011, 71: 6327-6337. 10.1158/0008-5472.CAN-11-0304.
Gatenby RA, Gillies RJ: A microenvironmental model of carcinogenesis. Nat Rev Cancer. 2008, 8: 56-61. 10.1038/nrc2255.
Goldie JH, Coldman AJ: A mathematic model for relating the drug sensitivity of tumors to their spontaneous mutation rate. Cancer Treat Rep. 1979, 63: 1727-1733.
Michor F, Nowak MA, Iwasa Y: Evolution of resistance to cancer therapy. Curr Pharm Des. 2006, 12: 261-271. 10.2174/138161206775201956.
Komarova N: Stochastic modeling of drug resistance in cancer. J Theor Biol. 2006, 239: 351-366. 10.1016/j.jtbi.2005.08.003.
Komarova NL, Wodarz D: Drug resistance in cancer: principles of emergence and prevention. Proc Natl Acad Sci U S A. 2005, 102: 9714-9719. 10.1073/pnas.0501870102.
Basanta D, Gatenby RA, Anderson ARA: Exploiting evolution to treat drug resistance: combination therapy and the double bind. Mol Pharm. 2012, 9: 914-921. 10.1021/mp200458e.
Swierniak A, Smieja J: Analysis and optimization of drug resistant and phase-specific cancer chemotherapy models. Math Biosci Eng MBE. 2005, 2: 657-670.
Kimmel M, Axelrod DE: Mathematical models of gene amplification with applications to cellular drug resistance and tumorigenicity. Genetics. 1990, 125: 633-644.
Kimmel M, Swierniak A, Polanski A: Infinite-dimensional model of evolution of drug resistance of cancer cells. J Math Syst Estim Contr. 1998, 8 (1): 1-16.
Smieja J, Swierniak A: Different models of chemotherapy taking into account drug resistance stemming from gene amplification. Int. J. Appl Math Comp Sci. 2003, 13.3: 297-306.
Harnevo LE, Agur Z: Drug resistance as a dynamic process in a model for multistep gene amplification under various levels of selection stringency. Cancer Chemother Pharmacol. 1992, 30: 469-476. 10.1007/BF00685599.
Foo J, Michor F: Evolution of resistance to targeted anti-cancer therapies during continuous and pulsed administration strategies. PLoS Comput Biol. 2009, 5: e1000557-10.1371/journal.pcbi.1000557.
Diaz LA, Williams RT, Wu J, Kinde I, Hecht JR, Berlin J, Allen B, Bozic I, Reiter JG, Nowak MA, Kinzler KW, Oliner KS, Vogelstein B: The molecular evolution of acquired resistance to targeted EGFR blockade in colorectal cancers. Nature. 2012, 486: 537-540.
Bozic I, Reiter JG, Allen B, Antal T, Chatterjee K, Shah P, Moon YS, Yaqubie A, Kelly N, Le DT, Lipson EJ, Chapman PB, Diaz LA, Vogelstein B, Nowak MA: Evolutionary dynamics of cancer in response to targeted combination therapy. Elife. 2013, 2: e00747-10.7554/eLife.00747.
Sawyers C: Targeted cancer therapy. Nature. 2004, 432: 294-297. 10.1038/nature03095.
Fister KR, Panetta JC: Optimal control applied to cell-cycle-specific cancer chemotherapy. SIAM J Appl Math. 2000, 60: 1059-1072. 10.1137/S0036139998338509.
Byrne HM, Owen MR, Alarcon T, Murphy J, Maini PK: Modelling the response of vascular tumours to chemotherapy: a multiscale approach. Math Models Methods Appl Sci. 2006, 16: 1219-1241. 10.1142/S0218202506001522.
De Pillis LG, Gu W, Radunskaya AE: Mixed immunotherapy and chemotherapy of tumors: modeling, applications and biological interpretations. J Theor Biol. 2006, 238: 841-862. 10.1016/j.jtbi.2005.06.037.
De Pillis L, Radunskaya A: A mathematical model of immune response to tumor invasion. Proc Second MIT Conf Comput Fluid Solid Mech KJ Bathe Ed Fluid Solid Mech. 2003, Comput: Solid
Jackson TL, Byrne HM: A mathematical model to study the effects of drug resistance and vasculature on the response of solid tumors to chemotherapy. Math Biosci. 2000, 164: 17-38. 10.1016/S0025-5564(99)00062-0.
Panetta JC, Adam J: A mathematical model of cycle-specific chemotherapy. Math Comput Model. 1995, 22: 67-82. 10.1016/0895-7177(95)00112-F.
Powathil GG, Adamson DJA, Chaplain MAJ: Towards predicting the response of a solid tumour to chemotherapy and radiotherapy treatments: clinical insights from a computational model. PLoS Comput Biol. 2013, 9: e1003120-10.1371/journal.pcbi.1003120.
Powathil G, Kohandel M, Sivaloganathan S, Oza A, Milosevic M: Mathematical modeling of brain tumors: effects of radiotherapy and chemotherapy. Phys Med Biol. 2007, 52: 3291-3306. 10.1088/0031-9155/52/11/023.
Powathil GG, Gordon KE, Hill LA, Chaplain MAJ: Modelling the effects of cell-cycle heterogeneity on the response of a solid tumour to chemotherapy: biological insights from a hybrid multiscale cellular automaton model. J Theor Biol. 2012, 308: 1-19.
Cojocaru L, Agur Z: A theoretical analysis of interval drug dosing for cell-cycle-phase-specific drugs. Math Biosci. 1992, 109: 85-97. 10.1016/0025-5564(92)90053-Y.
Bozic I, Allen B, Nowak MA: Dynamics of targeted cancer therapy. Trends Mol Med. 2012, 18: 311-316. 10.1016/j.molmed.2012.04.006.
Billy F, Ribba B, Saut O, Morre-Trouilhet H, Colin T, Bresch D, Boissel J-P, Grenier E, Flandrois J-P: A pharmacologically based multiscale mathematical model of angiogenesis and its use in investigating the efficacy of a new cancer treatment strategy. J Theor Biol. 2009, 260: 545-562. 10.1016/j.jtbi.2009.06.026.
Kirschner DE, Jackson TL, Arciero JC: A mathematical model of tumor-immune evasion and siRNA treatment. Discrete Contin Dyn Syst - Ser B. 2003, 4: 39-58.
Dewanji A, Luebeck EG, Moolgavkar SH: A generalized Luria-Delbrück model. Math Biosci. 2005, 197: 140-152. 10.1016/j.mbs.2005.07.003.
Kronik N, Kogan Y, Vainstein V, Agur Z: Improving alloreactive CTL immunotherapy for malignant gliomas using a simulation model of their interactive dynamics. Cancer Immunol Immunother CII. 2008, 57: 425-439. 10.1007/s00262-007-0387-z.
Nanda S, Moore H, Lenhart S: Optimal control of treatment in a mathematical model of chronic myelogenous leukemia. Math Biosci. 2007, 210: 143-156. 10.1016/j.mbs.2007.05.003.
Wein LM, Wu JT, Kirn DH: Validation and analysis of a mathematical model of a replication-competent oncolytic virus for cancer treatment: implications for virus design and delivery. Cancer Res. 2003, 63: 1317-1324.
Mok W, Stylianopoulos T, Boucher Y, Jain RK: Mathematical modeling of herpes simplex virus distribution in solid tumors: implications for cancer gene therapy. Clin Cancer Res Off J Am Assoc Cancer Res. 2009, 15: 2352-2360. 10.1158/1078-0432.CCR-08-2082.
Rockne R, Alvord EC, Rockhill JK, Swanson KR: A mathematical model for brain tumor response to radiation therapy. J Math Biol. 2008, 58: 561-578.
Richard M, Kirkby KJ, Webb RP, Kirkby NF: A mathematical model of response of cells to radiation. Nucl Instruments Methods Phys Res Sect B Beam Interactions Mater Atoms. 2007, 255: 18-22. 10.1016/j.nimb.2006.11.077.
Enderling H, Anderson ARA, Chaplain MAJ: A model of breast carcinogenesis and recurrence after radiotherapy. PAMM. 2007, 7: 1121701-1121702. 10.1002/pamm.200700362.
De Boer RJ, Hogeweg P, Dullens HF, De Weger RA, Den Otter W: Macrophage T lymphocyte interactions in the anti-tumor immune response: a mathematical mode. J Immunol Baltim Md 1950. 1985, 134: 2748-2758.
Matzavinos A, Chaplain MAJ, Kuznetsov VA: Mathematical modelling of the spatio-temporal response of cytotoxic T-lymphocytes to a solid tumour. Math Med Biol J IMA. 2004, 21: 1-34. 10.1093/imammb/21.1.1.
Kirschner D, Panetta JC: Modeling immunotherapy of the tumor-immune interaction. J Math Biol. 1998, 37: 235-252. 10.1007/s002850050127.
De Pillis LG, Radunskaya AE, Wiseman CL: A validated mathematical model of cell-mediated immune response to tumor growth. Cancer Res. 2005, 65: 7950-7958.
Rentsch CA, Biot C, Gsponer JR, Bachmann A, Albert ML, Breban R: BCG-mediated bladder cancer immunotherapy: identifying determinants of treatment response using a calibrated mathematical model. PLoS One. 2013, 8: e56327-10.1371/journal.pone.0056327.
Wei H-C: A numerical study of a mathematical model of pulsed immunotherapy for superficial bladder cancer. Jpn J Ind Appl Math. 2013, 30: 441-452. 10.1007/s13160-013-0107-3.
Rihan FA, Safan M, Abdeen MA, Abdel Rahman D: Qualitative and computational analysis of a mathematical model for tumor-immune interactions. J Appl Math. 2012, Article ID 475720: 19-
This work was supported by the Austria Science Fund (Projects Doktoratskolleg W11 Molecular Cell Biology and Oncology and SFB F21 Cell Proliferation and Cell Death in Tumors), the Tiroler Standortagentur (Bioinformatics Tyrol), and by the Austrian Research Promotion Agency (FFG), project Oncotyrol. We apologize to the authors of papers not cited in this review due to space constraints.
The authors declare that they have no competing interests.
RG, ME, PC, HH, and ZT carried out literature search and wrote the manuscript. ZT conceived the study. All authors read and approved the final manuscript.