Research article Special Issues

A mathematical study of the impact of cell plasticity on tumour control probability

  • The tumour control probability (TCP) is a treatment planning tool that evaluates the probability of tumour eradication and helps in the assessment of the relative efficacy of different radiotherapy regimens. The response of tumours to radiation differs greatly even between patients with same types of cancers. Tumour heterogeneity or cellular diversity among cancer cells has a pronounced impact on the success of administered radiotherapy protocols. Tumour heterogeneity can be explained using the cancer stem cells (CSCs) hypothesis, which posits that CSCs are responsible for tumour initiation and propagation as well as therapeutic resistance. Moreover, the existence of plasticity or bidirectional transition between CSCs and non-CSCs indicates that, sometimes, non-CSCs appear to mimic CSC phenotypes, resulting in an increase in resistance. Here, we have developed a stochastic model to investigate the impact of plasticity on the efficacy of radiotherapy. The effect of plasticity on TCP is explored by applying the model to standard and hyper-fractionated schedules for a three week period of treatment as well as standard, hyper-fractionated, and accelerated hyper-fractionated schedules with an equal total dose of 30 Gy. Our results confirm that tumour control becomes increasingly difficult in the presence of plasticity as well as for the most resistant tumours. For the case with equal total dose, it is observed that increasing fractionation, at first enhances the probability of CSCs and tumour removal, but ultimately results in lower TCPS+P and TCPS. In addition, the combination of radiotherapy and targeted therapy (with increasing CSC differentiation) improves both the probability of CSC and tumour removal, in the absence of plasticity. However, in the presence of plasticity, the impact of combination therapy is not significant.

    Citation: Farinaz Forouzannia, Sivabal Sivaloganathan, Mohammad Kohandel. A mathematical study of the impact of cell plasticity on tumour control probability[J]. Mathematical Biosciences and Engineering, 2020, 17(5): 5250-5266. doi: 10.3934/mbe.2020284

    Related Papers:

    [1] Ghassen Haddad, Amira Kebir, Nadia Raissi, Amira Bouhali, Slimane Ben Miled . Optimal control model of tumor treatment in the context of cancer stem cell. Mathematical Biosciences and Engineering, 2022, 19(5): 4627-4642. doi: 10.3934/mbe.2022214
    [2] Hasitha N. Weerasinghe, Pamela M. Burrage, Dan V. Nicolau Jr., Kevin Burrage . Agent-based modeling for the tumor microenvironment (TME). Mathematical Biosciences and Engineering, 2024, 21(11): 7621-7647. doi: 10.3934/mbe.2024335
    [3] Qiaojun Situ, Jinzhi Lei . A mathematical model of stem cell regeneration with epigenetic state transitions. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1379-1397. doi: 10.3934/mbe.2017071
    [4] Samantha L Elliott, Emek Kose, Allison L Lewis, Anna E Steinfeld, Elizabeth A Zollinger . Modeling the stem cell hypothesis: Investigating the effects of cancer stem cells and TGF−β on tumor growth. Mathematical Biosciences and Engineering, 2019, 16(6): 7177-7194. doi: 10.3934/mbe.2019360
    [5] Abdulhamed Alsisi, Raluca Eftimie, Dumitru Trucu . Non-local multiscale approach for the impact of go or grow hypothesis on tumour-viruses interactions. Mathematical Biosciences and Engineering, 2021, 18(5): 5252-5284. doi: 10.3934/mbe.2021267
    [6] J. Ignacio Tello . On a mathematical model of tumor growth based on cancer stem cells. Mathematical Biosciences and Engineering, 2013, 10(1): 263-278. doi: 10.3934/mbe.2013.10.263
    [7] Erin N. Bodine, K. Lars Monia . A proton therapy model using discrete difference equations with an example of treating hepatocellular carcinoma. Mathematical Biosciences and Engineering, 2017, 14(4): 881-899. doi: 10.3934/mbe.2017047
    [8] Bingshuo Wang, Wei Li, Junfeng Zhao, Natasa Trisovic . Longtime evolution and stationary response of a stochastic tumor-immune system with resting T cells. Mathematical Biosciences and Engineering, 2024, 21(2): 2813-2834. doi: 10.3934/mbe.2024125
    [9] Donggu Lee, Sunju Oh, Sean Lawler, Yangjin Kim . Bistable dynamics of TAN-NK cells in tumor growth and control of radiotherapy-induced neutropenia in lung cancer treatment. Mathematical Biosciences and Engineering, 2025, 22(4): 744-809. doi: 10.3934/mbe.2025028
    [10] Hsiu-Chuan Wei . Mathematical modeling of tumor growth: the MCF-7 breast cancer cell line. Mathematical Biosciences and Engineering, 2019, 16(6): 6512-6535. doi: 10.3934/mbe.2019325
  • The tumour control probability (TCP) is a treatment planning tool that evaluates the probability of tumour eradication and helps in the assessment of the relative efficacy of different radiotherapy regimens. The response of tumours to radiation differs greatly even between patients with same types of cancers. Tumour heterogeneity or cellular diversity among cancer cells has a pronounced impact on the success of administered radiotherapy protocols. Tumour heterogeneity can be explained using the cancer stem cells (CSCs) hypothesis, which posits that CSCs are responsible for tumour initiation and propagation as well as therapeutic resistance. Moreover, the existence of plasticity or bidirectional transition between CSCs and non-CSCs indicates that, sometimes, non-CSCs appear to mimic CSC phenotypes, resulting in an increase in resistance. Here, we have developed a stochastic model to investigate the impact of plasticity on the efficacy of radiotherapy. The effect of plasticity on TCP is explored by applying the model to standard and hyper-fractionated schedules for a three week period of treatment as well as standard, hyper-fractionated, and accelerated hyper-fractionated schedules with an equal total dose of 30 Gy. Our results confirm that tumour control becomes increasingly difficult in the presence of plasticity as well as for the most resistant tumours. For the case with equal total dose, it is observed that increasing fractionation, at first enhances the probability of CSCs and tumour removal, but ultimately results in lower TCPS+P and TCPS. In addition, the combination of radiotherapy and targeted therapy (with increasing CSC differentiation) improves both the probability of CSC and tumour removal, in the absence of plasticity. However, in the presence of plasticity, the impact of combination therapy is not significant.


    Tumour heterogeneity, which arises due to genetic and phenotypic diversity as well as environmental differences among cancer cells, has a fundamental impact on treatment outcomes. In particular, the resistance of CSCs to radiotherapy can lead to treatment failure and tumour recurrence [1,2]. This cellular diversity has been observed between different patients and within a single tumour. The cancer stem cell (CSC) hypothesis posits that tumour growth is governed by a scarce subpopulation of cancer cells undergoing symmetric and asymmetric proliferation to regenerate themselves and produce other lineages of cancer cells, thus, contributing to intra-tumour heterogeneity. CSCs autofluorescence analysis and cell surface protein expressions have been used to identify CSCs in different types of cancers [3]. However, cell surface protein biomarkers are not perfect and early progenies of CSCs may share the same protein biomarker as CSCs. Furthermore, experimental observations suggest that sometimes non-CSCs display CSC properties, suggesting some degree of plasticity between CSCs and non-CSCs [4,5]. Several mathematical models have been developed to incorporate the effect of CSCs hypothesis and plasticity in cancer cell dynamics and mammosphere formation assays [6,7,8,9,10,11,12,13,14,15]. However, the impact of plasticity on the effectiveness of radiotherapy protocols has not been investigated, and this is the direction we study in this paper.

    A large proportion of patients suffering from cancer, receive radiation therapy as part of their treatment. The aim of radiotherapy is to achieve a high probability of local tumour control while at the same time maintaining a low risk of associated side effects. TCP is a metric used to differentiate among various radiotherapy protocols and is defined as the probability of tumour-cell extinction by the end of the treatment. Consequently, TCP can be used to predict the performance of potential radiotherapy protocols and suggest the one with an optimum outcome. Several mathematical and computational approaches have been developed in the literature to study TCP [16,17,18,19,20,21,22]. Stochastic models based on Poisson statistics have been used to investigate the probability of tumour control. One well-studied model considers a simple birth-death master equation that includes the stochastic effect of cell kill due to radiation. Thus, TCP is defined as the probability of no cancer cells remaining at the end of treatment [22]. The model is based on clonal evolution theory, which claims that carcinogenesis arises due to random mutations that occur in a single cell. But, based on the CSC hypothesis, the elimination of CSCs is fundamental to achieving a cure, since CSCs are capable of initiating and reestablishing tumour growth. Thus, Dhawan et al. (2013) developed a unidirectional hierarchical model (based on the CSC hypothesis) to determine the probability of eliminating CSCs from a heterogeneous cell population [23]. However, the model does not consider the potential plasticity between non-CSCs and CSCs.

    Stochastic models are useful in the study of the behaviour of a small subset of cells. In this paper, the TCP is calculated when a small number of tumour cells remain and cellular plasticity is taken into consideration. In addition, Protein biomarkers are frequently used to identify CSCs. For example, CD34highCD38low and CD44highCD24low are used to identify CSCs for leukaemia and breast tumors respectively. However, it has been suggested that these biomarkers may not be perfect and early generations of progenitors may share the same biomarker [23]. Moreover, experimental studies show that tumour growth occurs starting from either a single biomarker positive cell (cells bearing biomarkers designed for CSCs isolation) or a single negative biomarker cell (cells not bearing biomarkers designed for CSCs isolation) [10]. Here, we present a stochastic model to study the impact of plasticity on tumour control probability. The model considers two sub-populations of positive biomarker cells (CSCs) and negative biomarker cells (non-CSCs). The radiation induced cell kill rate is assumed to be different between CSCs (S) and non-CSCs (P), to account for the various radio-sensitivities among cancer cells. The cancer kill rate is also considered to be a step function such that the kill rate is high during radiation and lower otherwise. Deriving an analytical solution is not feasible when plasticity is incorporated in the model. Thus, a modified Gillespie algorithm for the reactions with rates changing discontinuously is used to solve the time evolution of the stochastic model and calculate the TCP [24]. Consequently, the probability of CSC removal (TCPS) and the probability of tumour eradication (TCS+P) are evaluated for three different radiotherapy schedules: standard, hyper-fractionated, and accelerated hyper-fractionated. In addition, TCPS and TCS+P are calculated for a therapy that combines radiotherapy and targeted therapy in the presence and absence of plasticity. Both TCPS and TCS+P noticeably decrease for the most resistant tumours. Furthermore, applying combination therapy does not substantially improve the control of either the CSC population or the tumour, in the presence of plasticity.

    In this paper we investigate how plasticity between non-CSCs and CSCs affects the TCP. For this purpose, a two compartment model is used to split the total population of cells into sub-populations of stem cells (S) and non-CSCs (P). Stem cells have the potential to go through unlimited cell divisions to replicate themselves and to replace non-CSCs. Furthermore, it is assumed that there is a degree of plasticity in the system and that non-CSCs can revert to stem cells. In addition, cells can die independently of each other at a rate of Γi, for i=S,P, (stem cell, non-CSCs). Therefore, we consider the following division pathways.

    SS+S:ρS,SP:ρSP,PP+P:ρP,PS:ρPS,S0:ΓS(t),P0:ΓP(t), (2.1)

    where ρS and ρP are, respectively, the rates of self renewal for stem cells and proliferation for non-CSCs. In addition, stem cells can replace non-CSCs at rate ρSP, and non-CSCs can replace stem cells at rate ρPS. Cells can die at rate Γi(t), for i=S,P, (stem cell, non-CSCs).

    CSCs are capable of unlimited proliferation to maintain a tumour. Furthermore, they display resistance to radiotherapy through unregulated radiation-induced DNA repair mechanisms after exposure to treatment. Thus, the elimination of CSCs is essential to control a tumour. Consequently, the probability of eradicating CSCs is a key element in measuring the effectiveness of any particular treatment. In this context, the elimination of CSCs when no plasticity exists in the tumour (i.e., ρPS=0) can lead to tumour control. The study of the probability of CSC removal in model 2.1 with ρPS=0 is equivalent to analyzing the extinction probability in a simple birth-death process, and is independent of the dynamic of non-CSCs as seen in [23] and briefly explained in the supplementary material. In the presence of plasticity, however, the removal of CSCs may not result in tumour control.

    The stochastic dynamics of the model (Eq 2.1) is described using the following probability distribution function for a population of nS stem cells and nP non-CSCs at time t with the initial conditionpnS,nP(t0)=δnSnS0δnPnP0. Here, nS0 and nP0 denote the initial number of stem cells and non-CSCs at time t0.

    dpnS,nP(t)dt=ρSpnS1,nP(t)(nS1)+ρSP(nS+1)pnS+1,nP1(t)+ρP(nP1)pnS,nP1(t)+ρPS(nP+1)pnS1,nP+1(t)+ΓS(nS+1)pnS+1,nP(t)+ΓP(nP+1)pnS,nP+1(t)(ρS+ρSP+ΓS)nSpnS,nP(t)(ρP+ρPS+ΓP)nPpnS,nP(t). (2.2)

    Using the probability generating function U(S,P,t)=j=0i=0pi,j(t)SiPj, we obtain

    U(S,P,t)t=[(S1)(ρSSΓS(t))+(PS)ρSP]U(S,P,t)S+[(P1)(ρPPΓP(t))+(SP)ρPS]U(S,P,t)P, (2.3)

    with initial condition U(S,P,0)=Sn0SPn0P.

    The derivation of an analytical solution for Eq (2.3) is not possible. In addition, the cancer cell death rates are defined as step functions in which the radiotherapy-induced cancer kill rate is higher during each exposure duration and lower otherwise. Thus, the modified Gillespie algorithm [24] is employed to evaluate TCP in the presence of discontinuous death rates changes.

    The radiation induced cell kill is assumed to occur directly in time intervals when fractions of radiation are given. It is also suggested that CSCs are less likely to die than non-CSCs, due to the former's resistance to the therapy. Therefore, the cell kill rate at each treatment time interval ΔT and given dose d is defined as

    fi(t,d)=αidΔT+βid2ΔT, (2.4)

    with different radiobiological parameters αi and βi (i{S,P}) for CSCs and non-CSCs. The hazard function fi(t,d) is developed based on a linear quadratic model, and cell survival after each fraction can be represented by exp(αidβid2), (i{S,P}) [25]. Thus, the cell death rates for CSCs and non-CSCs (i{S,P}) can be considered to be the following step function

    Γi(t)={αidjΔT+βid2jΔTt[tj,tj+ΔT],0otherwise, (2.5)

    where tj and dj are the initial time and the given dose of the jth fraction of radiotherapy, respectively. The values of the model (Eq 2.1) and radio-sensitivity parameters are selected from Forouzannia et al. [7]. Thus, ρS=0.2 (day1), ρSP=0.7 (day1), ρP=0.1 (day1), ρPS=0.05 (day1) and radiobiological parameters are αS=0.14 (Gy1), αP=0.41 (Gy1), βS=0.048 (Gy2), and βP=0.17 (Gy2). The duration of treatment at each fraction of radiation ΔT is assumed to be 15 minutes.

    Radiotherapy targets cancer cells and shrinks tumours, but CSCs have the capacity to escape the effects of ionizing radiation due to their resistance to radiation. In general, tumours rich in CSCs are very difficult to control relative to ones with a lower CSC fraction. As a result, study of the efficiency of radiotherapy protocols for eradicating CSCs is of importance for clinical oncologists. However, removing the CSC population may prove to be insufficient in controlling the disease in the presence of plasticity. Therefore, to investigate the impact of plasticity on tumour control, we evaluate the probability of CSC removal, TCPS, and the probability of tumour eradication, TCPS+P, for the model (Eq 2.1) in the presence and absence of plasticity.

    In this direction, TCPS and TCPS+P are computed for two categories: (ⅰ) standard (ST) and hyper-fractionated (HR) over a three week period of treatment; (ⅱ) standard, hyper-fractionated, and accelerated hyper-fractionated (AC) with an equal total dose of 30 Gy. The conventional treatment regimen delivers a dose of 2 Gy per fraction, once each weekday. For the hyper-fractionated and accelerated hyper-fractionated protocols, the respective doses of 1.2 Gy and 1.5 Gy are given twice each weekday. Although the hyper-fractionated schedule delivers a higher total dose than the standard regimen over a three week period of treatment, these two regimens have the same biological effective dose [7].

    Figure 1 represents TCPS and TCPS+P, applying a standard radiotherapy regimen in both the presence and absence of plasticity in the system. The curves are relatively close to each other in Figure 1a due to the high radiotherapy-induced cell kill rates. Consequently, the cell kill rates at each fraction of radiation are reduced by assuming that βS and βP equal zero, to better distinguish the differences between various possible cases (Figure 1b). The results confirm that both TCPS and TCPS+P are reduced when non-CSCs are able to behave like CSCs.

    Figure 1.  TCPS and TCPS+P for a standard radiotherapy schedule in the absence and presence of plasticity with initial numbers of cells n0S=100 and n0P=100. (a) αS=0.14, αP=0.41, βS=0.048, and βP=0.17. (b) βS=0 and βP=0.

    In category (ⅰ), Figure 2 shows TCPS and TCPS+P for standard and hyper-fractionated schedules. In general, the existence of plasticity leads to a smaller probability of CSC and tumour removal for these radiotherapy schedules. The probability of CSC removal is relatively close for these two radiotherapy schedules in the absence of plasticity, but the variation of the calculated TCPS for the two schedules is more distinct in the presence of plasticity (Figure 2a).

    Figure 2.  (a) TCPS and (b) TCPS+P for Standard and hyper-fractionated schedules in the absence and presence of plasticity. The initial numbers of cells are n0S=100 and n0P=100.

    Figure 3 is the graph of (a) the probability of CSC removal and (b) tumour extinction with respect to dose in the presence and absence of plasticity. The results confirm that a larger total dose is required for CSC and tumour extinction in the presence of plasticity. Furthermore, Figure 3 shows that although the hyper-fractionated schedule uses a larger total dose of radiation to obtain the same TCPS when there is no plasticity, standard and hyper-fractionated protocols nearly use the same amount of total dose to achieve the same TCPS in the presence of plasticity and the same TCPS+P in the presence and absence of plasticity (Figure 3a and 3b).

    Figure 3.  (a) TCPS and (b) TCPS+P with respect to dose for standard and hyperfractionated schedules in the absence and presence of plasticity. The initial numbers of cells are n0S=100 and n0P=100.

    In addition, increasing the plasticity in the tumour decreases the TCPS and TCPS+P, implying greater challenges in controlling the disease (Figure S.1 in supplementary figures). In Figures 4 and 5, TCPS and TCPS+P are graphed to explain the impact of a highly resistant CSC subpopulation on treatment efficacy in both the presence and absence of plasticity. The results confirm that reducing the sensitivity of CSCs to radiation yields a lower probability of control of both the CSC subpopulation and tumour. Consequently, radiotherapy protocols are much less effective against highly resistant tumours than less-resistant ones (Figures 4 and 5). Furthermore, the probability of CSC eradication is comparatively close for standard and hyper-fractionated schedules when plasticity does not exist in a tumour (Figure 4a). Consequently, increasing the number of fractions does not improve the probability of CSC removal significantly for extremely resistant tumours in the absence of plasticity.

    Figure 4.  (a) TCPS and (b) TCPS+P for standard and hyperfractionated schedules in the absence of plasticity for different radiosensitivities among CSCs.
    Figure 5.  (a) TCPS and (b) TCPS+P for the two radiotherapy protocols in the presence of plasticity for different radio-sensitivities among CSCs. The initial numbers of cells are n0S=100 and n0P=100.

    It is necessary to understand that CSC removal can ultimately contribute to tumour control when there is no plasticity in the system, although expecting this absence may not be biologically realistic. CSC elimination is not enough to attain a cure when plasticity exists, but is still important in decreasing the most-resistant subpopulation of cells. Consequently, it is crucial to consider TCPS+P as a measure for evaluating therapeutic regimen efficacy in the presence of plasticity, while simultaneously monitoring CSC elimination.

    Figures 6 and 7 show TCPS and TCPS+P for standard, hyper-fractionated, and accelerated hyper-fractionated schedules with an equal total dose of 30 Gy (category (ⅱ)). In general, the results are similar to those seen in Figures 2 and 3 except that accelerated hyper-fractionated schedule employs larger total dose of radiation to obtain similar TCPS and TCPS+P in the presence and absence of plasticity.

    Figure 6.  (a) TCPS and (b) TCPS+P for Standard, hyper-fractionated, and accelerated hyper-fractionated schedules in the absence and presence of plasticity. The initial numbers of cells are n0S=100 and n0P=100.
    Figure 7.  (a) TCPS and (b) TCPS+P with respect to dose for standard, hyper-fractionated, and accelerated hyper-fractionated schedules in the absence and presence of plasticity. The initial numbers of cells are n0S=100 and n0P=100.

    Figure S.2 depicts that increasing plasticity, decreases TCPS and TCPS+P, which is consistent with the results of Figure S.1. The results reported in Figures S.3 and S.4 also confirm the same outcomes shown in Figures 4 and 5 in which decreasing the sensitivity of cancer cells to the radiation decreases TCPS and TCPS+P. In addition, standard, hyper-fractionated, and accelerated hyper-fractionated schedules have almost the same probability of CSCs removal in the absence of plasticity. It is noteworthy that increasing fractionation for the case when the total administered dose is equal, initially improves TCPS and TCPS+P, but eventually leads to lower probability of CSCs and tumour elimination (Figures S.3 and S.4).

    Combination of treatments

    Experimental results suggest that the fraction of CSCs increases following ionizing radiation, due to the resistance of this subset of cells to the therapy [26,27]. This frequently results in relapse and treatment failure. In this direction, both in vitro and in vivo experiments have shown that certain bone morphogenetic proteins (BMPs) are capable of inducing positive biomarker cells (stem-like cancer cells) to differentiate into negative biomarker cells (non-CSCs) in brain tumours [28]. Thus, applying this targeted strategy is expected to escalate the differentiation of radioresistant cells into non-CSCs that are more sensitive to radiation and have less tumourigenic potential. Therefore, a reduction in the CSC pool can contribute to better therapeutic outcomes. However, the plastic transition from non-CSCs to CSCs will reverse this process and may diminish the impact of BMPs. Thus, it is of interest to investigate the impact on TCP of increasing CSC differentiation rates in the presence and absence of plasticity, and this will be discussed next.

    Here, we consider three types of combination treatments over 15 days, consisting of 10 fractions of radiation (the dose of 2Gy is administered per fraction once a day) over a period of 5 days of targeted strategy in order to trigger CSCs to differentiate into non-CSCs. The first protocol includes 10 days of radiation, followed by a targeted strategy of increasing the differentiation rate ρSP for 5 days (scheme 1). The second starts with 5 days of increasing CSC differentiation, followed by 10 days of radiation (scheme 2). Finally, the third involves 5 days of radiation, followed by 5 days of targeted therapy, and last, another 5 days of radiation (scheme 3). We consider a base protocol for each of these three schedules, denoted as schemes 1-b, 2-b, and 3-b, in which the ρSP does not change.

    Figure 8 shows TCPS and TCPS+P for schemes 1-b, 2-b, and 3-b in the absence and presence of plasticity. As explained in the previous section, both TCPS and TCPS+P decrease when plasticity exists. Comparing the results for these three schedules indicates that the minimum TCPS and TCPS+P are obtained by scheme 3-b, suggesting that a large gap between fractionations reduces the probability of CSC and tumor removal (solid and dashed black curves) in the absence and presence of plasticity. In addition, scheme 2-b reports the highest TCPS+P, but this schedule does not lead to the best control of the CSC population (solid and dashed red curves). Employing the targeted therapy in combination with radiotherapy (schemes 1–3) shows an improvement in TCPS and TCPS+P for scheme 2 when plasticity does not exist, with the most increase occurring in the former (TCPS). However, the impact of targeted therapy on schemes 1 and 3 is limited, with almost no increase in TCPS+P for either of the schemes (Figures 8a and 9a). On the other hand, combining the targeted therapy with radiation in the presence of plasticity enhances TCPS and TCPS+P somewhat for schemes 1–3, with almost no increase in TCPS+P for scheme 1 (Figures 8b and 9b). Similar trends are also observed when ρSP increases to 2.1. In Figure 9b, TCPS decreases over a short period of time and increases after for scheme 3. This behaviour occurs because reducing ρSP after 10 days in the presence of plasticity decreases CSC elimination. Consequently, the generation of CSCs from non-CSCs can reduce the probability of complete CSC removal.

    Figure 8.  TCPS and TCPS+P for schemes 1-b, 2-b, and 3-b in the absence and presence of plasticity, with initial numbers of cells n0S=100 and n0P=100.
    Figure 9.  TCPS and TCPS+P for schemes 1, 2, and 3 in the absence and presence of plasticity, with initial numbers of cells n0S=100 and n0P=100. ρSP is increased to 1.4 during targeted therapy.

    The above results demonstrate that increasing ρSP before radiotherapy has the most effect in improving the probability of CSC removal in the absence of plasticity. In addition, the probability of tumour control increases in this case. However, in the presence of plasticity the tumour control is complicated, and even triggering CSCs to differentiate does not enhance TCPS and TCPS+P significantly. Nevertheless, increasing CSC differentiation before radiotherapy leads to better tumour removal control. But TCPS is still less than the case when CSC differentiation increases after radiotherapy. In general, these outcomes demonstrate that understanding the heterogeneity of tumours is important in order to develop optimal protocols.

    Here, we have presented a stochastic model to investigate the impact of plasticity on the tumor control probability. The radiation induced cell kill rate has been modelled using a step function in which the cancer cell kill rate is high during radiotherapy and small otherwise. Therefore, the time evolution trajectory of a hierarchical stochastic model consisting of CSCs and non-CSCs has been calculated using a modified Gillespie algorithm for rates that are changing discontinuously. Thus, TCPS and TCS+P are defined as the probability of removing CSCs and the probability of eliminating a tumour, respectively. The response of cancer cells to the radiation varies among cancer cells with CSCs showing more resistance to the treatment. TCPS and TCPS+P have been computed for standard, and hyper-fractionated protocols for three weeks of treatment as well as standard, hyper-fractionated, and accelerated hyper-fractionated schedules with the equal total dose of 30 Gy.

    The existence of plasticity between non-CSCs and CSCs decreases the efficiency of radiotherapy. The lower sensitivity to radiation among CSCs also contributes to a lower probability of CSC and tumour removal. In addition, hyper-fractionated, and accelerated hyper-fractionated protocols initially improve the probability of CSCs and tumour removal, but eventually result in lower TCPS and TCPS+P. The results also confirm that combination therapy can improve TCPS and TCPS+P in the absence of plasticity, with the most increase in the former. However, combination therapy does not increase the probability of CSC and tumour removal appreciably in the presence of plasticity. Here, we have assumed that stem cells replace the non-CSCs at rate ρSP. It may be biologically more relevant if we considered that cells undergo differentiation, however we think the results would be quantitatively similar to those reported in this paper. The results indicate that cellular heterogeneity and the existence of plasticity in a tumour significantly impact the efficacy of treatment. Thus, it is important to consider this cellular diversity when determining an appropriate treatment protocol.

    In this paper, we have proposed a simple computational model to qualitatively study the impact of cellular plasticity on tumour control probability. Parameters are selected from [7], in which mamosphere formation data from [27,29] was used to estimate parameters. In [27] the breast cancer cell line is irradiated with a single dose or daily doses of 2 Gy. Thus, the impact of plasticity on TCP is studied when a standard schedule is employed. Moreover, we have been also interested in studying the effects of hyper-fractionated and accelerated hyper-fractionated schedules on TCP in the presence of plasticity. In addition, the TCP of 100% is achieved faster for cells irradiated in the in vitro experiments. Of course, the 100% TCP after a short time may not be clinically achievable; however, we emphasize that our work is a mathematical study and caution should be exercised in interpreting the results in a clinical context. However, we hope that the current work and results will provide impetus for further clinical and experimental investigations.

    (MK) and (SS) are grateful for financial support provided by the Natural Science and Engineering Research Council of Canada (NSERC).

    The authors declare that they have no conflict of interest.

    The stochastic dynamic of model 2.1 in the absence of plasticity (i.e., ρPS=0) can be explained using the following master equation showing the probability of having a population of nS stem cells and nP non-CSCs with the initial number of cells nS0 and nP0 at time t0.

    dpnS,nP(t)dt=ρS(nS1)pnS1,nP(t)ρSnSpnS,nP(t)+ρSP(nS+1)pnS+1,nP1(t)ρSPnSpnS,nP(t)+ρP(nP1)pnS,nP1(t)ρPnPpnS,nP(t)+ΓS(t)(nS+1)pnS+1,nP(t)ΓS(t)nSpnS,nP(t)+ΓP(t)(nP+1)pnS,nP(t)ΓP(t)nPpnS,nP(t). (4.1)

    The initial condition is given next with δi,j representing the Kronecker delta function:

    pnS,nP(t0)=δnSnS0δnPnP0. (4.2)

    Using the probability distribution function introduced above, we can define the marginal probability distribution of the number of CSCs as

    unS(t)=nPpnS,nP(t),

    with the following master equation

    dunS(t)dt=unS1ρS(nS1)(ρS+ρSP+ΓS(t))nSunS+unS+1(ρSP+ΓS(t))(nS+1). (4.3)
    Figure S1.  (a) TCPS and (b) TCPS+P for standard and hyperfractionated schedules with different dedifferentiation rates ρPS, between non-CSCs and CSCs. The initial numbers of cells are n0S=100 and n0P=100.
    Figure S2.  (a) TCPS and (b) TCPS+P for standard, hyperfractionated, and accelerated hyperfractionated schedules (category (ⅱ)) with different dedifferentiation rates ρPS, between non-CSCs and CSCs. The initial numbers of cells are n0S=100 and n0P=100.
    Figure S3.  (a) TCPS and (b) TCPS+P for standard, hyperfractionated, and accelerated hyperfractionated schedules in the absence of plasticity for different radiosensitivities among CSCs.
    Figure S4.  (a) TCPS and (b) TCPS+P for standard, hyperfractionated, and accelerated hyperfractionated schedules in the presence of plasticity for different radiosensitivities among CSCs. The initial numbers of cells are n0S=100 and n0P=100.


    [1] D. Hanahan, R. A. Weinberg, Hallmarks of cancer: the next generation, Cell, 144 (2011), 646-674.
    [2] F. Pajonk, E. Vlashi, W. H. McBride, Radiation resistance of cancer stem cells: the 4 R's of radiobiology revisited, Stem Cells, 28 (2010), 639-648.
    [3] S. Sheila, C. Hawkins, et al, Identification of human brain tumour initiating cells, Nature, 432 (2004), 396-401.
    [4] N. D. Marjanovic, R. A. Weinberg, C. L. Chaffer, Poised with purpose: cell plasticity enhances tumorigenicity, Cell Cycle, 12 (2013), 2713-2714.
    [5] N. D. Marjanovic, R. A. Weinberg, C. L. Chaffer, Cell plasticity and heterogeneity in cancer, Clin. Chem., 59 (2013), 168-179.
    [6] F. Forouzannia, Studies of tumor heterogeneity, tumor microenvironment, and radiotherapy: A mathematical and computational approach, Ph.D thesis, University of Waterloo, 2019.
    [7] F. Forouzannia, H. Enderling, M. Kohandel, Mathematical Modeling of the Effects of Tumor Heterogeneity on the Efficiency of Radiation Treatment Schedule, Bull. Math. Biol., 80 (2018), 283-293.
    [8] F. Forouzannia, S. Sivaloganathan, Cancer Stem Cells, the Tipping Point: Minority Rules?, Curr. Stem Cell Rep., 3 (2017), 240-247.
    [9] R. Ganguly, I. K. Puri, Mathematical model for the cancer stem cell hypothesis, Cell Proliferation, 39 (2006), 3-14.
    [10] P. B. Gupta, C. M. Fillmore, G. Jiang, S. D. Shapira, K. Tao, C. Kuperwasser, et al., Stochastic State Transitions Give Rise to Phenotypic Equilibrium in Populations of Cancer Cells, Cell, 146 (2011), 633-644.
    [11] F. Michor, T. P. Hughes, Y. Iwasa, S. Branford, N. P. Shah, C. L. Sawyers, et al., Dynamics of chronic myeloid leukaemia, Nature, 435 (2005), 1267-1270.
    [12] C. Turner, M. Kohandel, Investigating the link between epithelial-mesenchymal transition and the cancer stem cell phenotype: a mathematical approach, J. Theor. Biol., 265 (2010), 329-335.
    [13] S. L. Weekes, B. Barker, S. Bober, K. Cisneros, J. Cline, A. Thompson, et al., A multicompartment mathematical model of cancer stem cell-driven tumor growth dynamics, Bull. Math. Biol., 76 (2014), 1762-1782.
    [14] S. Zapperi, C. A. M. La Porta, Do cancer cells undergo phenotypic switching? The case for imperfect cancer stem cell markers, Sci. Rep., 2 (2012), 1-7.
    [15] D. Zhou, Y. Wang, B. Wu, A multi-phenotypic cancer model with cell plasticity, J. Theor. Biol., 357 (2014), 35-45.
    [16] A. Dawson, T. Hillen, Derivation of the Tumour Control Probability (TCP) from a Cell Cycle Model, Comput. Math. Methods Med., 7 (2006), 121-141.
    [17] J. Gong, M. M. Dos Santos, C. Finlay, T. Hillen, Are more complicated tumour control probability models better?, Math. Med. Biol., 30 (2011), 1-19.
    [18] W. S. Kendal, Technical Report A closed-form description of tumour control with fractionated radiotherapy and repopulation, Int J Radiat Biol, 73 (1998), 207-210.
    [19] S. L. Tucker, H. D. Thames, M. G. Taylor, How well is the probability of tumor cure after fractionated irradiation described by Poisson statistics?, J. Theor. Biol., 124 (1990), 273-282.
    [20] A. Y. Yakovlev, Comments on the distribution of clonogens in irradiated tumors, Radiat. Res., 134 (1993), 117-120.
    [21] M. Zaider, L. Hanin, Tumor control probability in radiation treatment, Med. Phys., 38 (2011), 574-583.
    [22] M. Zaider, G. N. Minerbo, Tumour control probability: a formulation applicable to any temporal protocol of dose delivery, Phys. Med. Biol., 45 (2000), 279.
    [23] A. Dhawan, M. Kohandel, R. Hill, S. Sivaloganathan, Tumour control probability in cancer stem cells hypothesis, PloS one, 9 (2014), 121-141.
    [24] V. Shahrezaei, J. F. Ollivier, P. S. Swain, Colored extrinsic fluctuations and stochastic gene expression, Mol. Syst. Biol., 4 (2008).
    [25] E. J. Hall, A. J. Giaccia, Radiobiology for the Radiologist, 7th edition, Lippincott Williams Wilkins, Philadelphia, 2012.
    [26] S. Bao, Q. Wu, R. E. McLendon, Y. Hao, Q. Shi, A. B. Hjelmeland, et al., Glioma Stem Cells Promote Radioresistance by Preferential Activation of the DNA Damage Response, Nature, 444 (2006), 756-760. doi: 10.1038/nature05236
    [27] C. Lagadec, E. Vlashi, L. Della Donna, Y. Meng, C. Dekmezian, K. Kim, et al., Survival and selfrenewing capacity of breast cancer initiating cells during fractionated radiation treatment, Breast Cancer Res., 12 (2010), R13.
    [28] S. G. M. Piccirillo, B. A. Reynolds,N. Zanetti, G. Lamorte, E. Binda, G. Broggi, et al., Bone morphogenetic proteins inhibit the tumorigenic potential of human brain tumour-initiating cells, Nature, 444 (2006), 761-765.
    [29] T. M. Phillips, W. H. Mcbride, F. Pajonk, The response of CD24-/low/CD44+ breast cancerinitiating cells to radiation, J. Natl. Cancer Inst., 24 (2006), 1777-1785.
  • This article has been cited by:

    1. Farinaz Forouzannia, Vahid Shahrezaei, Mohammad Kohandel, The impact of random microenvironmental fluctuations on tumor control probability, 2021, 509, 00225193, 110494, 10.1016/j.jtbi.2020.110494
    2. Anuraag Bukkuri, Kenneth J. Pienta, Robert H. Austin, Emma U. Hammarlund, Sarah R. Amend, Joel S. Brown, A mathematical investigation of polyaneuploid cancer cell memory and cross-resistance in state-structured cancer populations, 2023, 13, 2045-2322, 10.1038/s41598-023-42368-8
    3. Anuraag Bukkuri, Modeling stress-induced responses: plasticity in continuous state space and gradual clonal evolution, 2024, 143, 1431-7613, 63, 10.1007/s12064-023-00410-3
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3988) PDF downloads(220) Cited by(3)

Figures and Tables

Figures(13)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog