- Open Access
Modeling drug resistance in a conjoint normal-tumor setting
Theoretical Biology and Medical Modelling volume 12, Article number: 3 (2015)
In this paper, we modify our previously developed conjoint tumor-normal cell model in order to make a distinction between tumor cells that are responsive to chemotherapy and those that may show resistance.
Using this newly developed core model, the evolution of three cell types: normal, tumor, and drug-resistant tumor cells, is studied through a series of numerical simulations. In addition, we illustrate critical factors that cause different dynamical patterns for normal and tumor cells. Among these factors are the co-dependency of the normal and tumor cells, the cells’ response mechanism to a single or multiple chemotherapeutic treatment, the drug administration sequence, and the treatment starting time.
The results provide us with a deeper understanding of the possible evolution of normal, drug-responsive, and drug-resistant tumor cells during the cancer progression, which may contribute to improving the therapeutic strategies.
Assessing the evolution of cancer, in the presence of surrounding normal cells, is the subject of many biomedical studies. Recently reported evidence strongly indicates that the dynamics of tumor cells and the surrounding normal cells are not independent of each other and may be mutually tuned [1–8]. Examination of the coupled population dynamics of tumor and normal cell populations can potentially provide substantial knowledge that may contribute to the identification of more effective therapeutic interventions, particularly in aging populations. Among the variety of research papers in this field, many are based on the analysis of mathematical and computational models. In many of these models, the growth of normal and tumor cells are considered to be independent and are expressed by such functions as the Gompertz, the logistic, and the exponential equations [9–13]. However, the mutual interaction of tumor cells with surrounding normal cells, which was first mathematically introduced in a conjoint model by Witten , could shed light on some of the complex patterns that can be detected during cancer progression [15–17].
The interaction of tumor and normal cells is not the exclusive factor causing different dynamical patterns during cancer progression, The interaction of cells with the host immune system, therapeutic agents such as chemotherapy, immune therapy, or any other therapeutic interactions are additional factors which can influence the evolution patterns of the cell populations [18–33].
While researchers continuously improve cancer treatment strategies, one of the most serious obstacles in cancer treatment are related to drug resistance, where the chemotherapeutic treatments do not lead to the hoped for outcome. The issues related to the drug resistance have been broadly studied from a variety of different perspectives [34–38].
This work aims to contribute to a deeper understanding of drug resistance effects on cancer progression through the analysis of a new mathematical model and its concomitant computational simulation for a coupled tumor-normal cell framework that is more aligned with experimental evidence. To simulate the population evolution of our model, we have used Mathematical V7.0. Model parameter values are estimated based on values previously introduced in the literature and are given in Table 1 of this paper. Additionally, the other parameters in some parts are varied in order to study the system’s evolution.
This paper is structured as follows: in The basic conjoint tumor-normal cell model Section, we briefly review the normal-tumor cell conjoint model. In Conjoint core model in a chemo-resistance setting section, we introduce the drug resistance assumptions and subsequently modify the conjoint model to make a distinction between tumor-responsive and tumor-resistant cells. In Chemo-treatment strategies in a resistance setting section, we include the effects of chemotherapeutic treatment to the modified conjoint model and we examine and discuss the dynamics of the system. We conclude and examine future research directions.
The basic conjoint tumor-normal cell model
Feizabadi & Witten  extended the earlier work of Witten  proposing the following generalized model to describe the inter-connection between normal and tumor cells. The core model equation system is given by:
where T(t), N(t), K T , K N , r T , r N are respectively the total number of tumor cells at time t, the total number of normal cells at time t, the carrying capacity for the tumor cells, the carrying capacity for the normal cells, and the per capita growth rate for the tumor and normal cells, and f T (N), f N (T) are the functional rules relating normal-to-tumor and tumor-to-normal interaction respectively . Witten and Feizabadi  have shown that one possible set of coupled, nonlinear equations for the tumor-normal cell system may be expressed follows:
where T, N, K T , K N , r T , r N are previously defined. In each equation, the second terms represent the interaction between tumor and normal cells. Here, β and κ have the units of 1/time. Also, for consistency, ρ 0 and ρ 1 have units of cells. T* is the critical size of the tumor and as the size of tumor exceeds the critical size, the normal cells growth rate decreases. Figure 1(a) illustrates the time evolution of normal and tumor cells in a hypothetical environment in which they grow independently (uncoupled) from one another, and where each cell population follows a Gompertzian-like behavior. In this figure the growth parameters are considered to be identical for both normal and tumor cells. In Figure 1(b) and 1(c), the conjoint growth is added to the model. As can be seen by the different parameter values, the growth of the normal and tumor cells can be affected as a result of the cellular interactions. The ability of tumor cells to inhibit the normal cell’s growth increases as the population of tumor cells passes the critical value T. In Figure 1(c) the growth of tumor cells is delayed due to the influence of the surrounding normal cells introduced by a higher value for the interaction parameter β. Figure 1(d) illustrates the case in which the normal cells have died out as a result of the strong interaction effect from the tumor cells. These results are not surprising as this is essentially a competitive exclusion model.
In the next portion of the paper we consider a modified conjoint model in which we make a distinction between the group of tumor cells that are responsive to one type of chemotherapeutic agent and those which are resistant to that same chemotherapeutic agent.
Mathematical models and results
Conjoint core model in a chemo-resistance setting
One of the ongoing challenges to maximizing chemotherapeutic success in cancer treatment is the long-standing challenge of tumor cell resistance to single or multiple drug cocktails . This mechanism, known as chemo-resistance, is complex and depends upon many factors including but not limited to the specific drug, specific tumor, or the specific host’s defense mechanism . Coupled with chemo-resistance is the challenge of age-related sensitivity or insensitivity to various drug cocktails. Thus, a dose that might not be lethal in a 20 year-old patient could well be lethal in a 60 year-old patient.
Drug resistance is classified into three major categories. The first category is associated with pharmacologic resistance or when the drug cannot effectively reach the tumor site due to insufficient pharmo-kinetic properties. The second category is rooted in the biochemistry of the tumor cell, for example when the drugs are not active at the tumor cell sites. The third category is when chemo-resistance results from genetic mutation of the tumor cells [42, 43].
In order to overcome drug resistance, we need to improve treatment efficacy by better understanding the resistance mechanisms and their effect on the cancer progression. This is a complex challenge and, so far has remained beyond traditional clinical and experimental examination. The complexity of the problem has led investigators to further develop their understanding using modeling and simulation methods. In fact, this challenge has been the subject of many theoretical and computational studies [44–51]. In the upcoming section of the paper we focus on this problem by introducing a chemo-resistant tumor cell component to our model. We modify the model as follows.
First, we rely on reported evidence indicating that metastatic tumors with larger sizes are more likely to show resistance to chemotherapeutic drugs [52, 53]. Therefore, the effect of normal cells in reducing the large population of tumor cells is not significant. Consequently, we have assumed that the second term is equation 2a is ignorable. For the purposes of discussion, we have assumed that the drug resistant tumor cells are created as a result of tumor cell mutation. They become resistant tumor cells with a mutation rate μ. Further, we assume that the drug resistant tumor cells also grow under the logistic growth law where the population growth rate, rR, is modified by the density dependent term associated to the total number of tumor cells. In our model, K R is the carrying capacity for the drug resistant tumor cells. Lastly, we remember that the population of the normal cells is controlled by the tumor cell population. Combining all of these assumptions along with our earlier model equations, we obtain the following new equation system:
The behavior of the normal, drug sensitive tumor, and drug resistant tumor cells is simulated in Figure 2. Understanding the evolution of each component becomes more critical in connection with the treatment of the system with chemotherapeutic agents that will be discussed in the next section.
In the drug resistance model, as the population of total tumor cells which now includes both the responsive and the resistant tumor cells, passes the critical value T* the normal cell population decreases in number. In Figure 2(a) the growth rate of both the drug responsive and the resistant tumor cells are considered to be identical. In this case the population of the resistant tumor cells is larger due to the fact that mutation of responsive tumor cells continually decreases the population of the responsive tumor cell population and subsequently increases the size of the resistant tumor cell population. In Figure 2(b), the growth rate of responsive tumor cells is higher than that of the resistant tumor cells. Here, we can see that for a period of time, the population of the responsive tumor cell population is larger. However, ultimately the population of the drug resistant tumor cells becomes higher than that of the tumor responsive cells. In Figure 2(c), the tumor-normal cell coupling coefficient is increased slightly. Under this new condition, the population of normal cells has become smaller than the critical value, T* = 5*105. Therefore, the tumor cells overwhelm the normal cells and the normal cells die out of the system faster than before, In addition, the population of the drug resistant tumor cells continues to grow. We next consider the effect of adding a chemotherapeutic agent to our system.
Chemo-treatment strategies in a resistance setting
The conjoint model, in the presence of the chemo-resistant tumor cells, may also be modified to consider the introduction of chemotherapeutic agents. To systematically investigate the evolution of the cells, we have simulated the system’s dynamics under the following conditions. We first assume that due to the drug resistance, the first chemotherapeutic agent introduced to the system has a cytotoxic effect only upon the drug responsive tumor cells, T. Due to the effect of this toxicity, the population of tumor cells decreases following an interaction with this drug. Secondly, we have simulated the dynamics of the system under a combination therapy, where the second chemotherapeutic agent is effective only on the drug resistant tumor cells, T R . Finally, the effects of the drug cocktail are studied when the time of drug administration is varied.
As suggested by Gardner  and used in other studies [29, 33], the drug interaction may be structured as a ϕ (1 − e− MC)ϕ. Here, φ is the cell population number of the three types of cells: T, T R, N. The parameter C is the concentration of the drug at the site, M is the pharmacokinetic factor, and aφ is the response factor. The function F(C) = a ϕ (1 − e− MC) is the fraction cell kill for a given amount of drug “C”. In the presence of the chemotherapeutic agents, the previous mutation model may now be modified as follows:
In the first two simulations, Figure 3(a,b), the conjoint tumor-normal cell population model is simulated where there is no mutation and hence no resistant tumor cell population. In a drug-free system, the coupling effect and the decrease in of normal cells can be observed as the tumor cells exceed the critical value T*. This same system is then simulated when the tumor cells interact with anti-tumor drugs and the toxicity of the drug kills the tumor cells. As the size of tumor cells decrease, the normal cells recover and subsequently return to a higher population number. In Figure 3(a-f), the number of normal cells increases slightly beyond the carrying capacity (K R = 1*106). At this point which, for instance, can be seen around t = 50 in Figure 2(b), the tumor cell population is less than the magnitude of the critical population, T*; therefore, the second term of equation 4 is positive and slightly increases the population of normal cells. In fact, normal cells react to the presence of small groups of tumor cells in the system.
In Figure 3(c), tumor cells we return to the two population tumor cell model; resistant and non-resistant tumor cells. The appearance of the drug resistant tumor cells at t ~ 50 days and their subsequent growth is illustrated in this figure. Given the mutation rate of μ = 10−3, the population of drug responsive tumor cells decreases. However, the population of drug resistant tumor cells increases over time as there are no chemotherapeutic agents that target this population.
Figure 3(d) illustrates the dynamics of the cell populations when the system interacts with an anti-tumor drug which is effective only on the drug responsive tumor cells. The drug is administered at t = 50 days. The drug responsive tumor cells decrease and die out of the system. Due to the chemotherapeutic treatment, the total number of tumor cells falls below the critical size of tumor cells. Therefore, the normal cells maintain higher population for a period of time. However, mutated drug-resistant tumor cells increase and their population will eventually pass the critical value. As a result, the normal cells start to decrease again. This kind of chemotherapeutic intervention can create a delay in possible organ failure by maintaining a higher number of normal cells for a period of time.
In the next simulation, we introduce the combination therapy protocol. A combination therapy is considered a more effective treatment strategy with cancers that show resistance to some of the chemotherapeutic agents. In this mode of intervention, while the tumor is treated by the recommended chemotherapeutic drug protocol, other chemotherapeutic drugs are also used in order to target those tumor cells that have developed defense mechanisms against the first type of chemotherapeutic agent. In Figure 3(e) illustrate the dynamics of the system under a multiple therapeutic protocol. In this simulation, both drugs are administered at the same time, t = 50. A lower dosage and therefore lower toxicity is considered for the anti-resistant tumor cell population. This mode of intervention was chosen due to the fact that, at the start of chemotherapy, drug responsive tumor cells have a higher population. Therefore, a higher drug dosage was considered for the non-resistant population. As can be seen in Figure 3(e), the normal appear, at first, to be stimulated by the tumor cell population growth but eventually return to their carrying capacity value while the drug sensitive tumor cell population dies out of the system as a result of interaction with the anti-cancer drug. In this simulation, due to the toxicity of the second type of the introduced anti-cancer drug, the maximum population of the drug resistant tumor cells is much smaller as compared to a case in which the tumor is treated with only a single chemotheraputic drug (Figure 3(d)).
In Figure 3(f), both treatments are supposed to be started simultaneously at an earlier time; t = 5 days. As can be seen in this figure, not only have the normal cells reached their carrying capacity, but also the growth of tumor cells has been delayed.
This work, a modification of our previous work, focuses on examining the dynamics of interconnected normal and tumor cells treated with chemotherapeutic agents, when some of the tumor cells show chemo-resistance. We examined these dynamics using a collection of different simulation parameters. Simulations demonstrated that in a conjoint system, normal cells enter a phase of diminished growth as the total number of tumor cells passes the magnitude of a critical tumor cell population . To control the population of tumor cells and the decrease of the population of normal cells, which may lead to organ failure, tumor cells can be treated with chemotherapeutic agents. In order to overcome the drug resistance, implementing a combination treatment is recommended. In a combination therapeutic approach, the dosage and the time of chemotherapy introduction play a critical role in minimizing the population of tumor cells, while maintaining the maximum population of the normal cells. According to our simulations, starting the combinatory therapy in the early stage of the cancer progression may lead to better control of the cancer progression as this treatment protocol can minimize the tumor cell population.
In our simulations, the growth rate and mutation rate of the cells are two other important factors that can potentially cause different evolution patterns. Another factor that plays a significant role in the system dynamics is the dosage of the anti-tumor drugs. It is more probable that a better response is achieved by increasing the drug dosage. However, since a majority of the chemotherapeutic drugs are toxic to normal cells and the host immune system, Consequently, the dosage and the level of toxicity must be carefully considered in order to minimize the potential damage to normal cells and to the patient.
In addition, damages that can be produced by chemotherapy are significant in the presence of an impaired immune system. The lack of inclusion of the effects of the immune system in our model is one of its limitations. Therefore, considering how the interaction of tumor cells with the host immune system may affect tumor progression are elements that can potentially be included in our model to achieve outcomes more aligned with clinical and biological observations.
All calculations were executed on an PC using Mathematica v7.0. Code is available from the first author.
Nguyen-Ngoc KV, Cheung KJ, Brenot A, Shamir ER, Gray RS, Hines WC, et al.: ECM microenvironment regulates collective migration and local dissemination in normal and malignant mammary epithelium.Proc Natl Acad Sci 2012,109(39):E2595–604.
Ko SY, Barengo N, Ladanyi A, Lee JS, Marini F, Lengyel E, et al.: HOXA9 promotes ovarian cancer growth by stimulating cancer-associated fibroblasts.J Clin Invest 2012,122(10):3603–17.
Ludford RJ, Barlow H: Sarcomatous transformation of the stroma of mammary carcinomas that stimulated fibroblastic growth in vitro.Cancer Res 1943, 5:257–64.
Ludford RJ, Barlow H: The influence of malignant cells upon the growth of fibroblasts in vitro.Cancer Res 1994, 8:694–703.
Kasuta H, Takaota T, Nagai Y: Interaction in Culture between Normal and Tumor Cells of Rats. In Cancer Cells in Culture. Edited by: Kasuta H. Tokyo, Japan: University of Tokyo Press; 1968.
Ranadive KJ, Bhide SV: Tissue interactions in normal and malignant cells. Little, Brown & Company: Boston, MA; 1962.
Foley JF, Kennedy BJ, Ross JD: A factor from HeLa cells promoting colonial growth of human fibroblast-like cells in culture.Cancer Res 1963, 23:368–71.
Smith AG, Cress HR: Effect of tumor deoxyribonucleic acid upon tissue cultures of fibroblasts.Lab Invest 1951, 10:898–907.
Clare SE, Nakhlis F, Panetta JC: Molecular biology of breast cancer metastasis. The use of mathematical models to determine relapse and to predict response to chemotherapy in breast cancer.Breast Cancer Res 2000, 2:430–5.
Guiot C, Degiorgis PG, Delsanto PP, Gabriele P, Deisboeck TS: Does tumor growth follow a “universal law”?J Theor Biol 2003, 225:147–51.
Hart D, Shochat E, Agur Z: The growth law of primary breast cancer as inferred from mammography screening trials data.Br J Cancer 1998, 78:382–7.
Simon R, Norton L: The Norton-Simon hypothesis: designing more effective and less toxic chemotherapeutic regimens.Nat Clin Pract Oncol 2006, 3:406–7.
Kozusko F, Bajzer Z: Combining Gompertzian growth and cell population dynamics.Math Biosci 2003, 185:153–67.
Witten TM: Modeling cellular aging and tumorigenic transformation.Math Comp Simul 1982, 24:572–84.
Witten TM: Population models of cellular aging: Theoretical and numerical issues. In Advances in Computer Methods for Partial Differential Equations VI. Edited by: Vichnevetsky R, Stepleman RS. New Jersey: IMACS/Rutgers University; 1987.
Witten TM: A mathematical model for the effects of a lymphokine-like ring shaped particle on the dynamics of a conjoint tumor-normal cell culture: Some implications for the aging cancer question. In Mathematics And Computers In Biomedical Applications. Edited by: DiLisi JC. New York: Elsevier/North-Holland; 1985.
Witten TM: Modeling the dynamics of conjoint tumor-normal cell culture systems: II. Further thoughts on the aging-cancer question. In Modeling Of Biomedical Systems. Edited by: Eisenfeld J, Witten TM. North-Holland: Amsterdam; 1986.
O'Byrne KJ, Dalgleish AG, Browning MJ, Steward WP, Harris AL: The relationship between angiogenesis and the immune response in carcinogenesis and the progression of malignant disease.Eur J Cancer 2000, 36:151–69.
Farrar JD, Katz KH, Windsor J, Thrush G, Scheuermann RH, Uhr JW, et al.: Cancer dormancy. VII. A regulatory role for CD8+ T cells and IFN-gamma in establishing and maintaining the tumor-dormant state.J Immunol 1999,162(5):2842–9.
Muller M, Gounari F, Prifti S, Hacker HJ, Schirrmacher V, Khazaie K: Eblac Z tumor dormancy in bone marrow and lymph nodes: active control of proliferating tumor cells by CD8+ immune T cells.Cancer Res 1998,58(23):5439–46.
Prehn RT, Prehn LM: The flip side of immune surveillance: immune dependency.Immunol Rev 2008, 222:341–56.
Prehn RT: The immune reaction as a stimulator of tumor growth.Science 1972, 176:170–1.
Fleisher MS, Loeb L: Immune reactions against tumor growth in animals with spontaneous tumors.J Med Res 1996, 34:1–19.
Bajzer B, Carr T, Josic K, Russell SJ, Dingli D: Modeling of cancer virotherapy with recombinant measles viruses.J Theor Biol 2008, 252:109–22.
Dingli D, Cascino MD, Josic K, Russell SJ, Bajzer Z: Mathematical modeling of cancer radiovirotherapy.Math Biosci 2006, 199:80–103.
Feizabadi MS, Volk C, Hirschbeck S: A two-compartment model interacting with dynamic drugs.Appl Math Lett 2009, 22:1205–9.
Feizabadi MS, Carbonara J: Two-compartment model interacting with proliferating regulatory factor.Appl Math Lett 2010, 23:30–3.
Feizabadi MS, Witten TM: Chemotherapy in cojoint aging-tumor systems: some simple models for addressing coupled aging-cancer dynamics.Theor Biol and Med Model 2010, 7:21.
de Pillis LG, Radunskaya A: The dynamics of an optimally controlled tumor model: A case study.Math Comp Model 2003, 37:1221–44.
Kirschner D, Panetta JC: Modeling immunotherapy of the tumor-immune interaction.J Math Biol 1998,37(3):235–52.
Feizabadi MS, Witten TM: Modeling the effects of a simple immune system and immunodeficiency on the dynamics of conjointly growing tumor and normal cells.Int J Biol Sci 2011,7(6):700–7.
Kirschner D, Tsygvintsev A: On the global dynamics of a model for tumor immunotherapy.Math Biosci and Eng 2009,6(3):573–83.
de Pillis LG, Gu W, Radunskaya AE: Mixed immunotherapy and chemotherapy of tumors: modeling, applications and biological interpretations.J Theor Biol 2006, 238:841–62.
Jasmine F, Franziska M: Evolution of acquired resistance to anti-cancer therapy.J Theor Biol 2014, 355:10–20.
Tomasetti C: On the probability of random genetic mutations for various types of tumor growth.Bull Math Biol 2012, 74:1379–95.
Wu Q, Yang Z, Nie Y, Shi Y, Fan D: Multi-drug resistance in cancer chemotherapeutics: Mechanisms and lab approaches.Cancer Lett 2014, 347:159–66.
Teicher BA: Cancer drug resistance. Totowa: Humana Press; 2006.
Tomasetti C, Levy D: Role of symmetric and asymmetric division of stem cells in developing drug resistance.Proc Natl Acad Sci U S A 2010,107(39):16766–71.
Witten TM: Modeling the dynamics of conjoint tumor-normal cell culture systems: II. Further thoughts on the aging-cancer question. In Modeling of Biomedical Systems. Edited by: Eisenfeld J, Witten TM. North-Holland, Amsterdam: North Holland Publishing Co; 1986:65–70.
Woodcock J, Griffin JP, Behrman RE: Development of novel combination therapies.N Engl J Med 2011, 364:985–7.
Holohan C, Van Schaeybroeck S, Longley DB, Johnston PG: Cancer drug resistance: an evolving paradigm.Nat Rev Cancer 2013,13(10):714–26.
Gottesman MM: Mechanisms of cancer drug resistance.Annu Rev Med 2002, 53:615–27.
Goldie JH, Coldman AJ: Drug Resistance in Cancer: Mechanisms and Models. Cambridge/New York: Cambridge University Press; 1998.
Tomasetti C, Levy D: An elementary approach to modeling drug resistance in cancer.Math Biosci Eng 2010,7(4):905–18.
Lavi O, Gottesman MM, Levy D: The dynamics of drug resistance: A mathematical perspective.Drug Resist Updat 2012, 15:90–7.
Wodarz D, Komarova NL: Emergence and prevention of resistance against small molecule inhibitors.Semin Cancer Biol 2005, 15:506–14.
Tomasetti C, Levy D: Role of symmetric and asymmetric division of stem cells in developing drug resistance.Proc Natl Acad Sci U S A 2010, 107:16766–71.
Panetta JC: A mathematical model of drug resistance: heterogeneous tumors.Math Biosci 1998, 147:41–61.
Panetta JC: A mathematical model of breast and ovarian cancer treated with paclitaxel.Math Biosci 1997, 146:89–113.
Chou TC: Theoretical basis, experimental design, and computerized simulation of synergism and antagonism in drug combination studies.Pharmacol Rev 2006, 58:621–81.
Kim M, Gillies RJ, Rejniak KA: Current advances in mathematical modeling of anti-cancer drug penetration into tumor tissues.Front Oncol 2013,18(3):1–11.
Liang Y, McDonnell S, Clynes M: Examining the relationship between cancer invasion/metastasis and drug resistance.Curr Cancer Drug Targets 2002,2(3):257–77.
Simstein R, Burow M, Parker A, Weldon C, Beckman B: Apoptosis, chemoresistance, and breast cancer: insights from the mcf-7 cell model system.Exp Biol Med 2003, 228:995–1003.
Gardner SN: A mechanistic, predictive model of dose–response curves for cell cycle phase-specific and nonspecific drugs.Cancer Res 2000, 60:1417–25.
MSF Acknowledges the support of Seton Hall University during her sabbatical leave.
Both authors declare that they have no competing interests.
The original conjoint cell culture equations were drawn from earlier work of TMW. Modifications were made by MSF. Computational work was carried out principally by the first author with suggestions from the second. All other work was executed jointly. There was no funding for this research project. Both authors have read and approved the final manuscript.