Next Article in Journal
Effects of Rho-Associated Kinase (Rock) Inhibitors (Alternative to Y-27632) on Primary Human Corneal Endothelial Cells
Previous Article in Journal
Clorfl86/RHEX Is a Negative Regulator of SCF/KIT Signaling in Human Skin Mast Cells
Previous Article in Special Issue
Distinct Dynamics of Migratory Response to PD-1 and CTLA-4 Blockade Reveals New Mechanistic Insights for Potential T-Cell Reinvigoration following Immune Checkpoint Blockade
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mathematical Model for Predicting Patient Responses to Combined Radiotherapy with CTLA-4 Immune Checkpoint Inhibitors

Department of Biomedical Engineering and of Biomedicine & Health Sciences, College of Medicine, The Catholic University of Korea, Seoul 06591, Republic of Korea
*
Authors to whom correspondence should be addressed.
Cells 2023, 12(9), 1305; https://doi.org/10.3390/cells12091305
Submission received: 6 February 2023 / Revised: 27 April 2023 / Accepted: 27 April 2023 / Published: 3 May 2023
(This article belongs to the Special Issue Cell-Cell Interaction Modelling of Cancer Immunotherapy Treatments)

Abstract

:
The purpose of this study was to develop a cell–cell interaction model that could predict a tumor’s response to radiotherapy (RT) combined with CTLA-4 immune checkpoint inhibition (ICI) in patients with hepatocellular carcinoma (HCC). The previously developed model was extended by adding a new term representing tremelimumab, an inhibitor of CTLA-4. The distribution of the new immune activation term was derived from the results of a clinical trial for tremelimumab monotherapy (NCT01008358). The proposed model successfully reproduced longitudinal tumor diameter changes in HCC patients treated with tremelimumab (complete response = 0%, partial response = 17.6%, stable disease = 58.8%, and progressive disease = 23.6%). For the non-irradiated tumor control group, adding ICI to RT increased the clinical benefit rate from 8% to 32%. The simulation predicts that it is beneficial to start CTLA-4 blockade before RT in terms of treatment sequences. We developed a mathematical model that can predict the response of patients to the combined CTLA-4 blockade with radiation therapy. We anticipate that the developed model will be helpful for designing clinical trials with the ultimate aim of maximizing the efficacy of ICI-RT combination therapy.

1. Introduction

Radiation therapy (RT) is traditionally considered a local cancer treatment due to its immune-suppressive characteristics. However, current evidence has demonstrated that radiation also possesses immune-stimulatory characteristics [1,2]. Potential immunogenic cell death via radiation triggered the investigation of radiation combined with immune checkpoint inhibitors (ICIs) for metastatic cancer [3].
A significant concern in adding RT to ICI is the development of radiation-induced immune-suppressive effects [4]. Despite the growing number of trials combining RT and ICI, whether the RT regimen in combination with ICI is the appropriate choice remains an open question, especially considering factors such as the timing of the radiation [5], sequencing [6], and patient selection [7]. Numerous studies [8,9,10,11,12] have investigated the integration of these treatments based on preclinical evidence demonstrating a synergistic interaction between them. However, it remains unclear how to optimally integrate these therapeutic modalities in the treatment of cancer patients. Beyond disease-specific factors, there exist numerous unanswered questions regarding optimal sequencing of radiation and ICI, as well as radiation dosing and target selection [13]. Therefore, it is essential to develop optimal strategies to improve the effectiveness of a combined treatment.
Mathematical models have the potential to help find optimal administration protocols, provide a deeper understanding of the dynamics, and aid in the design of clinical trials [14]. Various mathematical models were developed to simulate the interaction between cancer cells and the immune system. These models aim to predict the outcome of various cancer treatments, such as immunotherapy and radiation therapy. For example, one study investigated the effect of the immunotherapy and radiation therapy dose schedule in terms of immune response time characteristics through a mathematical model and included cyclic cytokine-based immunotherapy treatment in a tumor growth model [15]. Similarly, Sotolongo-Grau et al. proposed a dynamical system model for tumor–immune system interaction together with a method to mimic radiation therapy. A large population of virtual patients was simulated following an ideal radiation treatment. A characteristic parameter, the immune system/tumor efficiency ratio (ISTER), was introduced, and its dependence on treatment success and other features was studied [16]. Other studies have presented a mathematical model for predicting Kaplan–Meier survival curves of combined radiation and chemotherapy in patients with non-small cell lung disease (NSCLC) for use in clinical trial design [17].
In this paper, using a mathematical model, we study the response of a solid tumor to a combined RT and ICI treatment with various treatment regimens. We take hepatocellular carcinoma (HCC) treated with radiation and tremelimumab as a specific example. The first will describe the derivation of a mathematical parameter for CTLA-4 blockade with tremelimumab. The second will describe the patient response to RT and ICI with different treatment regimens, i.e., (1) patient selection and (2) the timing of radiation, in order to achieve maximized tumor eradication for both local and metastatic cancers.

2. Materials and Methods

The mathematical model was developed using the Python programming language (Python Software Foundation, Version 3.9, Wilmington, DE, USA).

2.1. Model Equations

The previously developed tumor–immune system interaction model was modified to account for the effects of anti-CTLA-4 immune checkpoint inhibitors [18]. The model was applied to predict patient responses after the combined radiation treatment with CTLA-4 immune checkpoint inhibitors.
Our model consists of four compartments—TI: irradiated tumor cells, I: dying tumor cells, TNI: non-irradiated (metastatic) tumor cells, and L: circulating lymphocytes.
d T I d t = a T ω 1 T I g + T I + T N I L δ * t R 1 e α T D T β T D T 2 T I
d I d t = δ * t R 1 e α T D T β T D T 2 T I r I
d T N I d t = a T N I ω 1 T N I g + T I + T N I L
d L d t = ω 2 1 + δ t r e m T I + T N I g + T I + T N I L + ω 3 I g + I L + s f L δ * t R 1 e α L D L L
The new formalism has only one additional term (𝛿trem) in Equation (4), describing the effects of CTLA-4 blockade. CTLA-4 blockade demonstrated a mild increase in the circulating lymphocytes, leading to enhanced T cell responses [19]. Thus, the new term (𝛿trem) increases the effects of ω2, which stimulate the lymphocytes to destroy cancer cells. The effect of the new term (𝛿trem) is assumed to decrease exponentially with the decay constant (λ), as in Equation (5). The decay constant is calculated as ln(2) divided by the half-life of tremelimumab.
δ t r e m = δ 0 , t r e m · e λ t
For term ( δ * t R ), which describes radiation cell death, instead of directly solving the differential equation, radiation cell death is separately accounted for by the instantaneous changes at distinct time points, as shown by the following:
T n + 1 = T n · e α T D T β T D T 2
I n + 1 = I n + T n · 1 e α T D T β T D T 2
L n + 1 = L n + 1 · e α L D L
Except for the new terms (𝛿trem and λ), other parameters were taken and fixed from the previously developed RT-only model. Table 1 summarizes the resulting average parameters of the population.
The irradiated tumor fraction (ITF) is defined as below:
I T F = T I T I + T N I

2.2. Patient Cohort

A total of 10,000 virtual patients were implemented into Python programming language for the simulation of a mathematical model. The virtual patient population in this study was assumed to have inoperable HCC. We assumed that each patient had a unique tumor volume, lymphocyte count, and tumor radio-sensitivity. All of the patient population was implemented as a normal distribution using the random syntax method of the Python program. In other words, the patient population had a certain form of distribution in this virtual clinical trial. The same distributions taken from previous studies were implemented in this study. This study assumed the tumor density to be 109 times the number of tumor cells per cubic cm. The averages and standard deviations of each distribution are described in detail in Table 1.

2.3. Model Fitting

Currently, immunotherapeutic agents used clinically show therapeutic effects in only a portion of cancer patients [29]. Therefore, the virtual patient population used in our study adjusted the distribution of the CTLA-4 blockade ICI response using the following parameters. The distribution of an additional new term (𝛿trem) in this study is fitted to describe the HCC patient responses to the tremelimumab—a complete response (CR): 0%, partial response (PR): 17.6%, and stable disease (SD): 58.8%. These percentages are based on the response evaluation criteria in solid tumors (RECIST) set of rules 1.1 [28]. To find the distribution, we first calculate the RECIST 1.1 responses for 10,000 patients with a constant 𝛿trem from 0 to 99. In the calculated RECIST 1.1 responses for a constant 𝛿trem, we find a distribution suitable for the CTLA-4 blockade ICI response mentioned above using the convolution method. Detailed model fitting procedures are described in Figure 1.

3. Results

A mathematical model was developed to predict a patient’s response to the combined radiotherapy with immune checkpoint inhibitors. The model was calibrated based on the responses reported in the CTLA-4 monotherapy clinical trial. We investigated the effects of irradiated tumor burden and treatment sequences under the radiation and CTLA-4 blockade combination regimen.

3.1. Calibration

Figure 2 shows the RECIST 1.1 response depending on the fixed 𝛿trem for 10,000 virtual patients. As the 𝛿trem increases, the PD rate decreases, while the PR and CR rates increase at a 𝛿trem value of 35. The SD rate increases up to 65 % at a 𝛿trem value of 40. Based on these results, we are able to set the mean and standard deviation of the 𝛿trem value and achieve RECIST 1.1 responses in patients enrolled in the ICI-only trial (NCT01008358) (Figure 3).

3.2. Simulation

Figure 4 shows the relative change in tumor size for 100 virtual patients over time after fitting the 𝛿trem value. One line represents the tumor growth dynamics of a single virtual patient with its corresponding RECIST1.1 response (colors).
Figure 5 shows the diameter changes in irradiated and non-irradiated tumors for 100 virtual patients. Radiation clearly plays a role in controlling irradiated tumors (waterfall plots on the top in Figure 5). However, radiation led to worse control of the non-irradiated tumors. The CTLA-4-based ICI monotherapy led to the most effective control of non-irradiated tumors. Although adding radiation to ICI achieves better control over irradiated tumors, no complete or partial responses are shown in the non-irradiated tumors.

3.2.1. Irradiated Tumor Burden

Figure 6 shows results for different treatment regimens in a simulated population of patients with HCC. The results show that the quantitative tumor size changes as a function of the irradiated tumor fraction. The irradiated tumor fraction is applied from 0.1 to 120% to obtain the results. For the irradiated tumor, the irradiated tumor fraction is not a significant factor in reducing the tumor size further. However, for combination regimens, non-irradiated tumor growth increases as the irradiated tumor fraction also increases.

3.2.2. Treatment Sequence

We investigated the effects of the sequencing of the RT and ICI combination (Figure 7). The tumor diameter changes were calculated for two simulated patient groups with 99% and 1% irradiated tumor fractions. In both combination cases, the early initiation of CTLA-4 ICI treatment led to the better management of the non-irradiated tumor (plus value of x-axis).

4. Discussion

Cancer immunotherapy is now an established part of the therapeutic options for solid tumors. More recently, radiation therapy has shown promising results in augmenting the anti-tumor effects of immune checkpoint inhibition, thus leading to an increased number of clinical trials [30,31,32]. Finding an optimal patient group is crucial not only to improve the treatment responses but also to save limited resources. In principle, the mathematical model is useful for modeling tumor growth under therapy [33]. The actionable model can assist in designing clinical trials by providing thought experiments that investigate potential treatment choices—a process known as an in silico clinical trial.

4.1. Treatment Modality

To confirm which treatment modalities would be effective for tumor control, changes in the size of the irradiated and non-irradiated tumors were obtained for each treatment modality (no treatment, ICI alone, RT alone, RT and ICI combination therapy), confirmed in Figure 5. It can seem that performing CTLA-4 blockade ICI and RT together would generally lead to a better prognosis than performing ICI solely. However, a slightly different result can be confirmed for the non-irradiated tumor size change. When RT is performed alone, the treatment method targets the local area, so it may be thought that there will be little effect on non-irradiated tumors. However, when CTLA-4 blockade ICI and RT are performed together, a worse prognosis can be confirmed than when ICI is performed solely (see the second and fourth graphs at the bottom of Figure 5). This result is thought to be a phenomenon caused by radiation damage to immune cells.

4.2. Irradiated Tumor Burden

In Figure 6, the tumor diameter change after 1 year for each treatment modality of the irradiated tumors and non-irradiated tumors was confirmed for the irradiated tumor fraction. The different irradiated tumor burdens lead to different amounts of tumor size changes for each treatment modality. The irradiated tumor burden is not a significant factor in controlling the irradiated tumor. This means that tumor control of the locoregional area is sufficient with RT alone. However, in non-irradiated tumor cases (blue and green lines in Figure 6), when adding CTLA-4-based ICI to RT, the model suggests that the smaller irradiated tumor fraction is better for reducing the non-irradiated tumor.

4.3. Treatment Sequence

Figure 7 shows the effect of sequencing and timing on treatment efficacy. The treatment application period of ICI and RT was confirmed at monthly intervals for up to 6 months. As can be seen in Figure 7a, negative numbers indicate that RT was performed first, and positive numbers indicate that ICI was performed first. We found that the responses are maximized if CTLA4-based ICI is delivered earlier than RT. Our results show that opposing behavior is demonstrated with PD-L1-based ICI [14]. This shows that even if it is not the PD-L1 inhibitor, when another inhibitor is used, the sequencing or timing may be different from when the CTLA-4 inhibitor is used. The model shows that early radiation-induced lymphocyte damage is a significant factor when adding CTLA-4-based ICI to RT. Therefore, we applied more indirect anti-tumor effects of CTLA-4-based ICI, reflected in w3 in the lymphocyte dynamic equations. Thus, the maintenance of the lymphocyte count is more important in inducing anti-tumor effects in CTLA-4-based ICI, compared to the increased direct effect of PD-L1 (w1).

4.4. Further Study

In this study, we built a mathematical model and simulated the responses of patients undergoing treatment with CTLA-4 immune checkpoint inhibitors and radiotherapy. There were limitations in our study, such as the weakness of the mathematical model, the inability to simulate various types of cancer, and the inability to simulate other inhibitors. We would like to address the issues in a further study.

5. Conclusions

We have developed a mathematical model that can predict patient responses when combined CTLA-4 blockade ICI and RT are performed. Through this model, it is possible to confirm which treatment modality would be effective to apply depending on whether an irradiated tumor or non-irradiated tumor needs to be controlled. It was confirmed that radiation therapy alone was sufficient for irradiated tumor control in a locoregional area, while combination therapy was effective for the metastatic tumor regions without irradiation. In addition, when performing ICI and RT together, it was confirmed that it was better to apply ICI therapy first before RT, and based on the results of Figure 7, we confirmed how much time should be left between treatment modalities in order for the treatment to be more effective. The developed mathematical model can assist in the design of clinical trials through thought experiments that investigate potential treatment choices—a process known as an in silico clinical trial. Therefore, we believe that the developed model can be used for patient selection and designing clinical trials in order to maximize the efficacy of CTLA-4 blockade ICI-RT combination therapy.

Author Contributions

Conceptualization, W.S.; methodology, W.S.; software, Y.K. and W.S.; validation, B.-Y.C., T.S.S. and W.S.; formal analysis, B.-Y.C., T.S.S. and W.S.; investigation, Y.K. and W.S.; resources, W.S.; data curation, Y.K.; writing—original draft preparation, Y.K.; writing—review and editing, B.-Y.C., T.S.S. and W.S.; visualization, Y.K. and W.S.; supervision, B.-Y.C., T.S.S. and W.S.; project administration, B.-Y.C., T.S.S. and W.S.; funding acquisition, T.S.S. and W.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the National Research Foundation of Korea (NRF, No. 2021R1C1C1005930) funded by the Korean Government (MSIT). This work was also supported by the National Research Foundation of Korea (NRF) grant funded by the Korean Government (MSIT) (Grant No. 2018R1A2B2005343).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data can be obtained from corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Antonia, S.J.; Villegas, A.; Daniel, D.; Vicente, D.; Murakami, S.; Hui, R.; Yokoi, T.; Chiappori, A.; Lee, K.H.; de Wit, M. Durvalumab after chemoradiotherapy in stage III non–small-cell lung cancer. N. Engl. J. Med. 2017, 377, 1919–1929. [Google Scholar] [CrossRef] [PubMed]
  2. Lumniczky, K.; Sáfrány, G. The impact of radiation therapy on the antitumor immunity: Local effects and systemic consequences. Cancer Lett. 2015, 356, 114–125. [Google Scholar] [CrossRef]
  3. Pike, L.R.; Bang, A.; Mahal, B.A.; Taylor, A.; Krishnan, M.; Spektor, A.; Cagney, D.N.; Aizer, A.A.; Alexander, B.M.; Rahma, O. The impact of radiation therapy on lymphocyte count and survival in metastatic cancer patients receiving PD-1 immune checkpoint inhibitors. Int. J. Radiat. Oncol. Biol. Phys. 2019, 103, 142–151. [Google Scholar] [CrossRef] [PubMed]
  4. Harrington, N.; Chambers, K.; Ross, W.; Filion, L. Radiation damage and immune suppression in splenic mononuclear cell populations. Clin. Exp. Immunol. 1997, 107, 417–424. [Google Scholar] [CrossRef]
  5. Young, K.H.; Baird, J.R.; Savage, T.; Cottam, B.; Friedman, D.; Bambina, S.; Messenheimer, D.J.; Fox, B.; Newell, P.; Bahjat, K.S. Optimizing timing of immunotherapy improves control of tumors by hypofractionated radiation therapy. PLoS ONE 2016, 11, e0157164. [Google Scholar] [CrossRef] [PubMed]
  6. Vanpouille-Box, C.; Alard, A.; Aryankalayil, M.J.; Sarfraz, Y.; Diamond, J.M.; Schneider, R.J.; Inghirami, G.; Coleman, C.N.; Formenti, S.C.; Demaria, S. DNA exonuclease Trex1 regulates radiotherapy-induced tumour immunogenicity. Nat. Commun. 2017, 8, 15618. [Google Scholar] [CrossRef]
  7. Iivanainen, S.; Koivunen, J.P. Possibilities of improving the clinical value of immune checkpoint inhibitor therapies in cancer care by optimizing patient selection. Int. J. Mol. Sci. 2020, 21, 556. [Google Scholar] [CrossRef] [PubMed]
  8. Karwacki, J.; Kiełbik, A.; Szlasa, W.; Sauer, N.; Kowalczyk, K.; Krajewski, W.; Saczko, J.; Kulbacka, J.; Szydełko, T.; Małkiewicz, B. Boosting the Immune Response—Combining Local and Immune Therapy for Prostate Cancer Treatment. Cells 2022, 11, 2793. [Google Scholar] [CrossRef] [PubMed]
  9. Tao, Z.; Ding, W.; Cheng, Z.; Feng, Y.; Kang, Z.; Qiu, R.; Zhao, S.; Hu, W.; Zhou, F.; Wu, D. Preclinical Study of Plasmodium Immunotherapy Combined with Radiotherapy for Solid Tumors. Cells 2022, 11, 3600. [Google Scholar] [CrossRef] [PubMed]
  10. Keisari, Y.; Kelson, I. The potentiation of anti-tumor immunity by tumor abolition with alpha particles, protons, or carbon ion radiation and its enforcement by combination with immunoadjuvants or inhibitors of immune suppressor cells and checkpoint molecules. Cells 2021, 10, 228. [Google Scholar] [CrossRef]
  11. Mondini, M.; Levy, A.; Meziani, L.; Milliat, F.; Deutsch, E. Radiotherapy–immunotherapy combinations–perspectives and challenges. Mol. Oncol. 2020, 14, 1529–1537. [Google Scholar] [CrossRef]
  12. Kalbasi, A.; June, C.H.; Haas, N.; Vapiwala, N. Radiation and immunotherapy: A synergistic combination. J. Clin. Investig. 2013, 123, 2756–2763. [Google Scholar] [CrossRef] [PubMed]
  13. Rajeev-Kumar, G.; Pitroda, S.P. Synergizing radiotherapy and immunotherapy: Current challenges and strategies for optimization. Neoplasia 2023, 36, 100867. [Google Scholar] [CrossRef] [PubMed]
  14. Sung, W.; Hong, T.S.; Poznansky, M.C.; Paganetti, H.; Grassberger, C. Mathematical modeling to simulate the effect of adding radiotherapy to immunotherapy and application to hepatocellular carcinoma. Int. J. Radiat. Oncol. Biol. Phys. 2021, 112, 1055–1062. [Google Scholar] [CrossRef] [PubMed]
  15. Rodriguez-Perez, D.; Sotolongo-Grau, O.; Riquelme, R.E.; Sotolongo-Costa, O.; Miranda, J.A.S.; Antoranz, J. Assessment of cancer immunotherapy outcome in terms of the immune response time features. Math. Med. Biol. J. IMA 2007, 24, 287–300. [Google Scholar] [CrossRef] [PubMed]
  16. Sotolongo-Grau, O.; Rodríguez-Pérez, D.; Santos-Miranda, J.; Sotolongo-Costa, O.; Antoranz, J. Immune system–tumour efficiency ratio as a new oncological index for radiotherapy treatment optimization. Math. Med. Biol. J. IMA 2009, 26, 297–307. [Google Scholar] [CrossRef]
  17. Geng, C.; Paganetti, H.; Grassberger, C. Prediction of treatment response for combined chemo-and radiation therapy for non-small cell lung cancer patients using a bio-mathematical model. Sci. Rep. 2017, 7, 13542. [Google Scholar] [CrossRef]
  18. Sung, W.; Grassberger, C.; McNamara, A.L.; Basler, L.; Ehrbar, S.; Tanadini-Lang, S.; Hong, T.S.; Paganetti, H. A tumor-immune interaction model for hepatocellular carcinoma based on measured lymphocyte counts in patients undergoing radiotherapy. Radiother. Oncol. 2020, 151, 73–81. [Google Scholar] [CrossRef] [PubMed]
  19. Fecci, P.E.; Ochiai, H.; Mitchell, D.A.; Grossi, P.M.; Sweeney, A.E.; Archer, G.E.; Cummings, T.; Allison, J.P.; Bigner, D.D.; Sampson, J.H. Systemic CTLA-4 blockade ameliorates glioma-induced changes to the CD4+ T cell compartment without affecting regulatory T-cell function. Clin. Cancer Res. 2007, 13, 2158–2167. [Google Scholar] [CrossRef]
  20. Park, Y.; Choi, D.; Lim, H.K.; Rhim, H.; Kim, Y.-s.; Kim, S.H.; Lee, W.J. Growth rate of new hepatocellular carcinoma after percutaneous radiofrequency ablation: Evaluation with multiphase CT. Am. J. Roentgenol. 2008, 191, 215–220. [Google Scholar] [CrossRef] [PubMed]
  21. Tai, A.; Erickson, B.; Khater, K.A.; Li, X.A. Estimate of radiobiologic parameters from clinical data for biologically based treatment planning for liver irradiation. Int. J. Radiat. Oncol. Biol. Phys. 2008, 70, 900–907. [Google Scholar] [CrossRef] [PubMed]
  22. Kuznetsov, V.A.; Makalkin, I.A.; Taylor, M.A.; Perelson, A.S. Nonlinear dynamics of immunogenic tumors: Parameter estimation and global bifurcation analysis. Bull. Math. Biol. 1994, 56, 295–321. [Google Scholar] [CrossRef]
  23. Serre, R.; Barlesi, F.; Muracciole, X.; Barbolosi, D. Immunologically effective dose: A practical model for immuno-radiotherapy. Oncotarget 2018, 9, 31812. [Google Scholar] [CrossRef] [PubMed]
  24. Sotolongo-Costa, O.; Molina, L.M.; Perez, D.R.; Antoranz, J.; Reyes, M.C. Behavior of tumors under nonstationary therapy. Phys. D Nonlinear Phenom. 2003, 178, 242–253. [Google Scholar] [CrossRef]
  25. Byun, H.K.; Kim, N.; Park, S.; Seong, J. Acute severe lymphopenia by radiotherapy is associated with reduced overall survival in hepatocellular carcinoma. Strahlenther. Onkol. 2019, 195, 1007–1017. [Google Scholar] [CrossRef]
  26. Nakamura, N.; Kusunoki, Y.; Akiyama, M. Radiosensitivity of CD4 or CD8 positive human T-lymphocytes by an in vitro colony formation assay. Radiat. Res. 1990, 123, 224–227. [Google Scholar] [CrossRef] [PubMed]
  27. Walker, R.; Poleszczuk, J.; Pilon-Thomas, S.; Kim, S.; Anderson, A.A.; Czerniecki, B.J.; Harrison, L.B.; Moros, E.G.; Enderling, H. Immune interconnectivity of anatomically distant tumors as a potential mediator of systemic responses to local therapy. Sci. Rep. 2018, 8, 9474. [Google Scholar] [CrossRef]
  28. Sangro, B.; Gomez-Martin, C.; de la Mata, M.; Iñarrairaegui, M.; Garralda, E.; Barrera, P.; Riezu-Boj, J.I.; Larrea, E.; Alfaro, C.; Sarobe, P. A clinical trial of CTLA-4 blockade with tremelimumab in patients with hepatocellular carcinoma and chronic hepatitis C. J. Hepatol. 2013, 59, 81–88. [Google Scholar] [CrossRef]
  29. Galluzzi, L.; Chan, T.A.; Kroemer, G.; Wolchok, J.D.; López-Soto, A. The hallmarks of successful anticancer immunotherapy. Sci. Transl. Med. 2018, 10, eaat7807. [Google Scholar] [CrossRef]
  30. Pilones, K.A.; Vanpouille-Box, C.; Demaria, S. Combination of radiotherapy and immune checkpoint inhibitors. Semin Radiat. Oncol. 2015, 25, 28–33. [Google Scholar] [CrossRef]
  31. Kang, J.; Demaria, S.; Formenti, S. Current clinical trials testing the combination of immunotherapy with radiotherapy. J. Immunother. Cancer 2016, 4, 51. [Google Scholar] [CrossRef] [PubMed]
  32. Cushman, T.R.; Caetano, M.S.; Welsh, J.W.; Verma, V. Overview of ongoing clinical trials investigating combined radiotherapy and immunotherapy. Immunotherapy 2018, 10, 851-0. [Google Scholar] [CrossRef] [PubMed]
  33. Watanabe, Y.; Dahlman, E.L.; Leder, K.Z.; Hui, S.K. A mathematical model of tumor growth and its response to single irradiation. Theor. Biol. Med. Model. 2016, 13, 6. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic diagram of model fitting procedure. The yellow highlighted part is a procedure different from the previous study. Abbreviations: RT = radiation therapy; ICI = immune checkpoint inhibitor; RECIST = response evaluation criteria in solid tumors.
Figure 1. Schematic diagram of model fitting procedure. The yellow highlighted part is a procedure different from the previous study. Abbreviations: RT = radiation therapy; ICI = immune checkpoint inhibitor; RECIST = response evaluation criteria in solid tumors.
Cells 12 01305 g001
Figure 2. The effect of the 𝛿trem (tremelimumab term) on RECIST1.1 responses for 10,000 virtual patients. The green, blue, orange, and red lines indicate complete response (CR), partial response (PR), stable disease (SD), and progressive disease (PD), respectively. Abbreviations: RECIST = response evaluation criteria in solid tumors.
Figure 2. The effect of the 𝛿trem (tremelimumab term) on RECIST1.1 responses for 10,000 virtual patients. The green, blue, orange, and red lines indicate complete response (CR), partial response (PR), stable disease (SD), and progressive disease (PD), respectively. Abbreviations: RECIST = response evaluation criteria in solid tumors.
Cells 12 01305 g002
Figure 3. Based on the RECIST 1.1 response to the constant 𝛿trem, a figure was used to find the distribution of tremelimumab suitable for the values for PR (17.6%), SD (58.8%), and PD (23.6%) [28]. (a) A heatmap showing the difference in relation to the reference according to the mean and standard deviation of tremelimumab. (b) Distribution of suitable-for-reference values and RECIST 1.1 responses. Abbreviations: PR = partial response; SD = stable disease; PD = progressive disease; RECIST = response evaluation criteria in solid tumors.
Figure 3. Based on the RECIST 1.1 response to the constant 𝛿trem, a figure was used to find the distribution of tremelimumab suitable for the values for PR (17.6%), SD (58.8%), and PD (23.6%) [28]. (a) A heatmap showing the difference in relation to the reference according to the mean and standard deviation of tremelimumab. (b) Distribution of suitable-for-reference values and RECIST 1.1 responses. Abbreviations: PR = partial response; SD = stable disease; PD = progressive disease; RECIST = response evaluation criteria in solid tumors.
Cells 12 01305 g003
Figure 4. Tumor diameter change over time for 100 virtual patients treated with tremelimumab. Each line represents the tumor dynamics of a single patient. The green, blue, orange, and red lines indicate complete response (CR), partial response (PR), stable disease (SD), and progressive disease (PD), respectively.
Figure 4. Tumor diameter change over time for 100 virtual patients treated with tremelimumab. Each line represents the tumor dynamics of a single patient. The green, blue, orange, and red lines indicate complete response (CR), partial response (PR), stable disease (SD), and progressive disease (PD), respectively.
Cells 12 01305 g004
Figure 5. Tumor size changes for 100 virtual patients after one year of treatment (upper: irradiated tumor, lower: non-irradiated tumor). Each bar represents a single patient. The RECIST response 1.1 is expressed as CR: green, PR: blue, SD: orange, and PD: dark red. Abbreviations: CR = complete response; PR = partial response; SD = stable disease; PD = progressive disease; RECIST = response evaluation criteria in solid tumors; RT = radiation therapy; ICI = immune checkpoint inhibitor.
Figure 5. Tumor size changes for 100 virtual patients after one year of treatment (upper: irradiated tumor, lower: non-irradiated tumor). Each bar represents a single patient. The RECIST response 1.1 is expressed as CR: green, PR: blue, SD: orange, and PD: dark red. Abbreviations: CR = complete response; PR = partial response; SD = stable disease; PD = progressive disease; RECIST = response evaluation criteria in solid tumors; RT = radiation therapy; ICI = immune checkpoint inhibitor.
Cells 12 01305 g005
Figure 6. Diameter change in the irradiated tumor and non-irradiated tumor after 1 year according to the irradiated tumor fraction ((a): irradiated tumor, (b): non-irradiated tumor). The green line indicates CTLA-4 blockade ICI + RT; the blue line indicates RT only; the gold line indicates CTLA-4 blockade ICI only; and the red line indicates none. Abbreviations: RT = radiation therapy; ICI = immune checkpoint inhibitor.
Figure 6. Diameter change in the irradiated tumor and non-irradiated tumor after 1 year according to the irradiated tumor fraction ((a): irradiated tumor, (b): non-irradiated tumor). The green line indicates CTLA-4 blockade ICI + RT; the blue line indicates RT only; the gold line indicates CTLA-4 blockade ICI only; and the red line indicates none. Abbreviations: RT = radiation therapy; ICI = immune checkpoint inhibitor.
Cells 12 01305 g006
Figure 7. Predicted results of tumor diameter change with different modality combination sequences. (a) Schematic to explain (b,c). Blue and green arrows indicate the start of treatment. On the time axis, negative values (− sign with red arrow in (a)) indicate that RT started before the ICI, whereas positive values (+ sign with red arrow in (a)) indicate that RT started after the ICI. (b) Response to treatment sequence when the irradiated tumor fraction was 99%. (c) Response to treatment sequence when the irradiated tumor fraction was 1%. Blue error bars indicate irradiated tumors and light blue error bars indicate non-irradiated tumors. Abbreviations: RT = radiation therapy; ICI = immune checkpoint inhibitor.
Figure 7. Predicted results of tumor diameter change with different modality combination sequences. (a) Schematic to explain (b,c). Blue and green arrows indicate the start of treatment. On the time axis, negative values (− sign with red arrow in (a)) indicate that RT started before the ICI, whereas positive values (+ sign with red arrow in (a)) indicate that RT started after the ICI. (b) Response to treatment sequence when the irradiated tumor fraction was 99%. (c) Response to treatment sequence when the irradiated tumor fraction was 1%. Blue error bars indicate irradiated tumors and light blue error bars indicate non-irradiated tumors. Abbreviations: RT = radiation therapy; ICI = immune checkpoint inhibitor.
Cells 12 01305 g007
Table 1. Summary of the parameters. LQ: linear quadratic.
Table 1. Summary of the parameters. LQ: linear quadratic.
ParameterFunctionValueRefs.
a Tumor growth 0.01   d 1 [20]
α T / β T Tumor—LQ cell death 14.3   G y [21]
f Lymphocyte decay rate 0.033   d 1 [22]
r Inactivated tumor cell decay rate 0.14   d 1 [23]
ω 1 Tumor-directed lymphocyte efficiency 0.119   d 1 [16,24]
ω 2
ω 3
Tumor/inactivated tumor lymphocyte
recruitment constant
0.003   d 1
0.009   d 1
[16,24]
g Half-saturation constant 7.330 × 10 10 [16,24]
s Lymphocyte regeneration 1.470 × 10 8 d 1 [25]
α T Tumor—LQ cell death 0.139   G y 1 [21]
α L Lymphocytes—LQ cell death 0.737   G y 1 [26,27]
δ t r e m Effectiveness of immune checkpoint inhibitor (tremelimumab)Normally distributed
(μ = 37 σ = 4)
[28]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kim, Y.; Choe, B.-Y.; Suh, T.S.; Sung, W. A Mathematical Model for Predicting Patient Responses to Combined Radiotherapy with CTLA-4 Immune Checkpoint Inhibitors. Cells 2023, 12, 1305. https://doi.org/10.3390/cells12091305

AMA Style

Kim Y, Choe B-Y, Suh TS, Sung W. A Mathematical Model for Predicting Patient Responses to Combined Radiotherapy with CTLA-4 Immune Checkpoint Inhibitors. Cells. 2023; 12(9):1305. https://doi.org/10.3390/cells12091305

Chicago/Turabian Style

Kim, Yongjin, Bo-Young Choe, Tae Suk Suh, and Wonmo Sung. 2023. "A Mathematical Model for Predicting Patient Responses to Combined Radiotherapy with CTLA-4 Immune Checkpoint Inhibitors" Cells 12, no. 9: 1305. https://doi.org/10.3390/cells12091305

APA Style

Kim, Y., Choe, B. -Y., Suh, T. S., & Sung, W. (2023). A Mathematical Model for Predicting Patient Responses to Combined Radiotherapy with CTLA-4 Immune Checkpoint Inhibitors. Cells, 12(9), 1305. https://doi.org/10.3390/cells12091305

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop