INTRODUCTION: While radiotherapy has long been recognized for its ability to directly ablate cancer cells through necrosis or apoptosis, radiotherapy-induced abscopal effect suggests that its impact extends beyond local tumor destruction thanks to immune response. Cellular proliferation and necrosis have been extensively studied using mathematical models that simulate tumor growth, such as Gompertz law, and the radiation effects, such as the linear-quadratic model. However, the effectiveness of radiotherapy-induced immune responses may vary among patients due to individual differences in radiation sensitivity and other factors. METHODS: We present a novel macroscopic approach designed to quantitatively analyze the intricate dynamics governing the interactions among the immune system, radiotherapy, and tumor progression. Building upon previous research demonstrating the synergistic effects of radiotherapy and immunotherapy in cancer treatment, we provide a comprehensive mathematical framework for understanding the underlying mechanisms driving these interactions. RESULTS: Our method leverages macroscopic observations and mathematical modeling to capture the overarching dynamics of this interplay, offering valuable insights for optimizing cancer treatment strategies. One shows that Gompertz law can describe therapy effects with two effective parameters. This result permits quantitative data analyses, which give useful indications for the disease progression and clinical decisions. DISCUSSION: Through validation against diverse data sets from the literature, we demonstrate the reliability and versatility of our approach in predicting the time evolution of the disease and assessing the potential efficacy of radiotherapy-immunotherapy combinations. This further supports the promising potential of the abscopal effect, suggesting that in select cases, depending on tumor size, it may confer full efficacy to radiotherapy.

Biotech Research Center Technology Innovation Institute Abu Dhabi United Arab Emirates

Faculty of Mathematics and Physics Charles University Prague Czechia

Genomics and molecular oncology unit Istituto Oncologico del Mediterraneo Viagrande Italy

Institute for Applied Computing National Research Council of Italy Rome Italy

Istituto Nazionale di Fisica Nucleare Sezione di Catania Catania Italy

Radiotherapy Unit REM Radioterapia Viagrande Italy

School of Medicine University Kore of Enna Enna Italy

Zobrazit více v PubMed

Gompertz B. Xxiv. on the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies. in a letter to francis baily, esq. frs &c. Philos Trans R Soc London. (1825) 115:513–83. doi: 10.1098/rstl.1825.0026 PubMed DOI PMC

Norton L. A gompertzian model of human breast cancer growth. Cancer Res. (1988) 48:7067–71. PubMed

Vaghi C, Rodallec A, Fanciullino R, Ciccolini J, Mochel JP, Mastri M, et al. . Population modeling of tumor growth curves and the reduced gompertz model improve prediction of the age of experimental tumors. PLoS Comput Biol. (2020) 16:e1007178. doi: 10.1371/journal.pcbi.1007178 PubMed DOI PMC

Verhulst P-F. Notice sur la loi que la population poursuit dans son accroissement. Correspondance Mathématique Physique. (1838) 10:113–21. doi: 10.4000/msh.2893 DOI

Hiam-Galvez KJ, Allen BM, Spitzer MH. Systemic immunity in cancer. Nat Rev Cancer. (2021) 21:345–59. doi: 10.1038/s41568-021-00347-z PubMed DOI PMC

Lin C-W, Xie J, Zhang D, Han KH, Grande G, Wu NC, et al. . Immunity against cancer cells may promote their proliferation and metastasis. Proc Natl Acad Sci. (2020) 117:426–31. doi: 10.1073/pnas.1916833117 PubMed DOI PMC

Zhao X, Shao C. Radiotherapy-mediated immunomodulation and anti-tumor abscopal effect combining immune checkpoint blockade. Cancers (Basel). (2020) 12:2762. doi: 10.3390/cancers12102762 PubMed DOI PMC

Puglisi C, Giuffrida R, Borzì G, Di Mattia P, Costa A, Colarossi C, et al. . Radiosensitivity of cancer stem cells has potential predictive value for individual responses to radiotherapy in locally advanced rectal cancer. Cancers (Basel). (2020) 12:2072–6694. doi: 10.3390/cancers12123672 PubMed DOI PMC

Puglisi C, Giuffrida R, Borzì G, Illari S, Caronia F, Di Mattia P, et al. . Ex vivo irradiation of lung cancer stem cells identifies the lowest therapeutic dose needed for tumor growth arrest and mass reduction in vivo . Front Oncol. (2022) 12:837400. doi: 10.3389/fonc.2022.837400 PubMed DOI PMC

Ngwa W, Irabor O, Schoenfeld J, Hesser J, Demaria S, Formenti S. Using immunotherapy to boost the abscopal effect. Nat Rev Cancer. (2018) 18:313–22. doi: 10.1038/nrc.2018.6 PubMed DOI PMC

Mole R. Whole body irradiation; radiobiology or medicine? Br J Radiol. (1953) 26:234–41. doi: 10.1259/0007-1285-26-305-234 PubMed DOI

Demaria S, Ng B, Devitt M, Babb J, Kawashima N, Liebes L, et al. . Ionizing radiation inhibition of distant untreated tumors (abscopal effect) is immune mediated. Int J Radiat Oncology Biology Phys. (2004) 58:862–70. doi: 10.1016/j.ijrobp.2003.09.012 PubMed DOI

Twyman-Saint Victor C, Rech A, Maity A, Rengan R, Pauken K, Stelekati E, et al. . Radiation and dual checkpoint blockade activate non-redundant immune mechanisms in cancer. Nature. (2015) 520:373–7. doi: 10.1038/nature14292 PubMed DOI PMC

Dovedi S, Cheadle E, Popple A, Poon E, Morrow M, Stewart R, et al. . Fractionated radiation therapy stimulates anti-tumor immunity mediated by both resident and infiltrating polyclonal t-cell populations when combined with pd1 blockade. Clin Cancer Res. (2017) 23:5514–26. doi: 10.1158/1078-0432.CCR-16-1673 PubMed DOI

Demaria S, Kawashima N, Yang A, Devitt M, Babb J, Allison J, et al. . Immune-mediated inhibition of metastases after treatment with local radiation and ctla-4 blockade in a mouse model of breast cancer. Clin Cancer Res. (2005) 11:728–34. doi: 10.1158/1078-0432.728.11.2 PubMed DOI

Ji D, Song C, Li Y, Xia J, Wu Y, Jia J, et al. . Combination of radiotherapy and suppression of tregs enhances abscopal antitumor effect and inhibits metastasis in rectal cancer. J ImmunoTherapy Cancer. (2018) 8:313–22. doi: 10.1136/jitc-2020-000826 PubMed DOI PMC

McNamee E, Korns Johnson D, Homann D, Clambey E. Hypoxia and hypoxia-inducible factors as regulators of t cell development, differentiation, and function. Immunol Res. (2013) 55:58–70. doi: 10.1007/s12026-012-8349-8 PubMed DOI PMC

Castorina P, Castorina L, Ferini G. Non-homogeneous tumor growth and its implications for radiotherapy: A phenomenological approach. J Personalized Med. (2021) 11:527. doi: 10.3390/jpm11060527 PubMed DOI PMC

Dewan MZ, Galloway AE, Kawashima N, Dewyngaert JK, Babb JS, Formenti SC, et al. . Fractionated but not single-dose radiotherapy induces an immune-mediated abscopal effect when combined with anti–ctla-4 antibodyfractionated radiation synergizes with immunotherapy. Clin Cancer Res. (2009) 15:5379–88. doi: 10.1158/1078-0432.CCR-09-0265 PubMed DOI PMC

Agur Z, Vuk-Pavlović S. Mathematical modeling in immunotherapy of cancer: personalizing clinical trials. Mol Ther. (2012) 20:1–2. doi: 10.1038/mt.2011.272 PubMed DOI PMC

Walker R, Enderling H. From concept to clinic: Mathematically informed immunotherapy. Curr Probl Cancer. (2016) 40(1):68–83. doi: 10.1016/j.currproblcancer.2015.10.004 PubMed DOI

Ng J, Dai T. Radiation therapy and the abscopal effect: a concept comes of age. Ann Transl Med. (2016) 4(6):118–20. doi: 10.21037/atm.2016.01.32 PubMed DOI PMC

Serre R, Benzekry S, Padovani L, Meille C, André N, Ciccolini J, et al. . Mathematical modeling of cancer immunotherapy and its synergy with radiotherapymodeling of radioimmunotherapy. Cancer Res. (2016) 76:4931–40. doi: 10.1158/0008-5472.CAN-15-3567 PubMed DOI

Gong C, Milberg O, Wang B, Vicini P, Narwal R, Roskos L, et al. . A computational multiscale agent-based model for simulating spatio-temporal tumour immune response to pd1 and pdl1 inhibition. J R Soc Interface. (2017) 14:20170320. doi: 10.1098/rsif.2017.0320 PubMed DOI PMC

Marconi R, Strolin S, Bossi G, Strigari L. A meta-analysis of the abscopal effect in preclinical models: Is the biologically effective dose a relevant physical trigger? PLoS One. (2017) 12:e0171559. doi: 10.1371/journal.pone.0171559 PubMed DOI PMC

Vanpouille-Box C, Alard A, Aryankalayil MJ, Sarfraz Y, Diamond JM, Schneider RJ, et al. . Dna exonuclease trex1 regulates radiotherapy-induced tumour immunogenicity. Nat Commun. (2017) 8:15618. doi: 10.1038/ncomms15618 PubMed DOI PMC

Chakwizira A, Ahlstedt J, Nittby Redebrandt H, Ceberg C. Mathematical modelling of the synergistic combination of radiotherapy and indoleamine-2, 3-dioxygenase (ido) inhibitory immunotherapy against glioblastoma. Br J Radiol. (2018) 91:20170857. doi: 10.1259/bjr.20170857 PubMed DOI PMC

Kosinsky Y, Dovedi SJ, Peskov K, Voronova V, Chu L, Tomkinson H, et al. . Radiation and pd-(l) 1 treatment combinations: immune response and dose optimization via a predictive systems model. J immunotherapy Cancer. (2018) 6:1–15. doi: 10.1186/s40425-018-0327-9 PubMed DOI PMC

Liu Y, Dong Y, Kong L, Shi F, Zhu H, Yu J. Abscopal effect of radiotherapy combined with immune checkpoint inhibitors. J Hematol Oncol. (2018) 11:1–15. doi: 10.1186/s13045-018-0647-8 PubMed DOI PMC

Valentinuzzi D, Simončič U, Uršič K, Vrankar M, Turk M, Jeraj R. Predicting tumour response to anti-pd-1 immunotherapy with computational modelling. Phys Med Biol. (2019) 64:025017. doi: 10.1088/1361-6560/aaf96c PubMed DOI

Friedrich T, Henthorn N, Durante M. Modeling radioimmune response—current status and perspectives. Front Oncol. (2021) 11:647272. doi: 10.3389/fonc.2021.647272 PubMed DOI PMC

Malinzi J, Basita KB, Padidar S, Adeola HA. Prospect for application of mathematical models in combination cancer treatments. Inf Med Unlocked. (2021) 23:100534. doi: 10.1016/j.imu.2021.100534 DOI

Bekker RA, Kim S, Pilon-Thomas S, Enderling H. Mathematical modeling of radiotherapy and its impact on tumor interactions with the immune system. Neoplasia. (2022) 28:100796. doi: 10.1016/j.neo.2022.100796 PubMed DOI PMC

Altrock PM, Liu LL, Michor F. The mathematics of cancer: integrating quantitative models. Nat Rev Cancer. (2015) 15:730–45. doi: 10.1038/nrc4029 PubMed DOI

Brady R, Enderling H. Mathematical models of cancer: when to predict novel therapies, and when not to. Bull Math Biol. (2019) 81:3722–31. doi: 10.1007/s11538-019-00640-x PubMed DOI PMC

Browning AP, Simpson MJ. Geometric analysis enables biological insight from complex non-identifiable models using simple surrogates. PLoS Comput Biol. (2023) 19:e1010844. doi: 10.1371/journal.pcbi.1010844 PubMed DOI PMC

Wheldon TE. Mathematical models in cancer research. Adam Hilger Publisher. (1988).

Guiot C, Degiorgis PG, Delsanto PP, Gabriele P, Deisboeck TS. Does tumor growth follow a “universal law”? J Theor Biol. (2003) 225:147–51. doi: 10.1016/S0022-5193(03)00221-2 PubMed DOI

Castorina P, Delsanto P, Guiot C. Classification scheme for phenomenological universalities in growth problems in physics and other sciences. Phys Rev Lett. (2006) 96:188701. doi: 10.1103/PhysRevLett.96.188701 PubMed DOI

Castorina P, Carcò D, Guiot C, Deisboeck TS. Tumor growth instability and its implications for chemotherapy. Cancer Res. (2009) 69:8507–15. doi: 10.1158/0008-5472.CAN-09-0653 PubMed DOI PMC

Maaß C, Sachs JP, Hardiansyah D, Mottaghy FM, Kletting P, Glatting G. Dependence of treatment planning accuracy in peptide receptor radionuclide therapy on the sampling schedule. EJNMMI Res. (2016) 6:1–9. doi: 10.1186/s13550-016-0185-8 PubMed DOI PMC

Jung H, Krüger H-J, Brammer I, Zywietz F, Beck-Bornholdt H-P. Cell population kinetics of the rhabdomyosarcoma r1h of the rat after single doses of x-rays. Int J Radiat Biol. (1990) 57:567–89. doi: 10.1080/09553009014552701 PubMed DOI

Hansen O, Grau C, Bentzen SM, Overgaard J. Repopulation in the sccvii squamous cell carcinoma assessed by an in vivo-in vitro excision assay. Radiotherapy Oncol. (1996) 39:137–44. doi: 10.1016/0167-8140(96)01728-8 PubMed DOI

O’Donoghue JA. The response of tumours with gompertzian growth characteristics to fractionated radiotherapy. Int J Radiat Biol. (1997) 72:325–39. doi: 10.1080/095530097143329 PubMed DOI

Braakman S, Pathmanathan P, Moore H. Evaluation framework for systems models. CPT Pharmacometrics Syst Pharmacol. (2022) 11:264–289. doi: 10.1002/psp4.12755 PubMed DOI PMC

Castorina P, Blanchard P. Unified approach to growth and aging in biological, technical and biotechnical systems. SpringerPlus. (2012) 1:1–5. doi: 10.1186/2193-1801-1-7 PubMed DOI PMC

Berendt MJ, North RJ. T-cell-mediated suppression of anti-tumor immunity. an explanation for progressive growth of an immunogenic tumor. J Exp Med. (1980) 151:69–80. doi: 10.1084/jem.151.1.69 PubMed DOI PMC

Gonzalez H, Hagerling C, Werb Z. Roles of the immune system in cancer: from tumor initiation to metastatic progression. Genes Dev. (2018) 32:1267–84. doi: 10.1101/gad.314617.118 PubMed DOI PMC

Castorina P, Carco’ D. Nutrient supply, cell spatial correlation and gompertzian tumor growth. Theory Biosci. (2021) 140:197–203. doi: 10.1007/s12064-021-00344-8 PubMed DOI PMC

Khajanchi S, Ghosh D. The combined effects of optimal control in cancer remission. Appl Mathematics Comput. (2015) 271:375–88. doi: 10.1016/j.amc.2015.09.012 DOI

Khajanchi S, Banerjee S. Quantifying the role of immunotherapeutic drug t11 target structure in progression of Malignant gliomas: Mathematical modeling and dynamical perspective. Math Biosci. (2017) 289:69–77. doi: 10.1016/j.mbs.2017.04.006 PubMed DOI

Royama T. Analytical population dynamics Vol. 10. UK: Springer Science & Business Media; (2012).

Van Leeuwen C, Oei A, Crezee J, Bel A, Franken N, Stalpers L, et al. . The alfa and beta of tumours: a review of parameters of the linear-quadratic model, derived from clinical radiotherapy studies. Radiat Oncol. (2018) 13:1–11. doi: 10.1186/s13014-018-1040-z PubMed DOI PMC

Khajanchi S, Banerjee S. Influence of multiple delays in brain tumor and immune system interaction with t11 target structure as a potent stimulator. Math Biosci. (2018) 302:116–30. doi: 10.1016/j.mbs.2018.06.001 PubMed DOI

Khajanchi S. Stability analysis of a mathematical model for glioma-immune interaction under optimal therapy. Int J Nonlinear Sci Numerical Simulation. (2019) 20:269–85. doi: 10.1515/ijnsns-2017-0206 DOI

Khajanchi S, Banerjee S. A strategy of optimal efficacy of t11 target structure in the treatment of brain tumor. J Biol Syst. (2019) 27:225–55. doi: 10.1142/S0218339019500104 DOI

Dataset Grace . (2008). Available at: https://plasma-gate.weizmann.ac.il/Grace/.

Nesseler JP, Lee M-H, Nguyen C, Kalbasi A, Sayre JW, Romero T, et al. . Tumor size matters—understanding concomitant tumor immunity in the context of hypofractionated radiotherapy with immunotherapy. Cancers. (2020) 12:714. doi: 10.3390/cancers12030714 PubMed DOI PMC

Lim K, Chan P, Dinniwell R, Fyles A, Haider M, Cho Y-B, et al. . Cervical cancer regression measured using weekly magnetic resonance imaging during fractionated radiotherapy: radiobiologic modeling and correlation with tumor hypoxia. Int J Radiat Oncol Biol Phys. (2008) 70:126–33. doi: 10.1016/j.ijrobp.2007.06.033 PubMed DOI

McMahon SJ. The linear quadratic model: usage, interpretation and challenges. Phys Med Biol. (2018) 64:01TR01. doi: 10.1088/1361-6560/aaf26a PubMed DOI

Filatenkov A, Baker J, Mueller A, Kenkel J, Ahn G, Dutt S, et al. . Ablative tumor radiation can change the tumor immune cell microenvironment to induce durable complete remissions. Clin Cancer Res. (2015) 21:3727. doi: 10.1158/1078-0432.CCR-14-2824 PubMed DOI PMC

Muraro E, Furlan C, Avanzo M, Martorelli D, Comaro E, Rizzo A, et al. . Local highdose radiotherapy induces systemic immunomodulating effects of potential therapeutic relevance in oligometastatic breast cancer. Front Immunol. (2017) 8:1476. doi: 10.3389/fimmu.2017.01476 PubMed DOI PMC

Klug F, Prakash H, Huber P, Seibel T, Bender N, Halama N, et al. . Low-dose irradiation programs macrophage differentiation to an inos+/m1 phenotype that orchestrates effective t cell immunotherapy. Cancer Cell. (2013) 24:589. doi: 10.1016/j.ccr.2013.09.014 PubMed DOI

Song C, Glatstein E, Marks L, Emami B, Grimm J, Sperduto P, et al. . Biological principles of stereotactic body radiation therapy (sbrt) and stereotactic radiation surgery (srs): Indirect cell death. Int J Radiat Oncology Biology Phys. (2021) 110:21. doi: 10.1016/j.ijrobp.2019.02.047 PubMed DOI

Ferini G, Valenti V, Tripoli A, Illari S, Molino L, Parisi S, et al. . Lattice or oxygen-guided radiotherapy: What if they converge? possible future directions in the era of immunotherapy. Cancers (Basel). (2021) 13:3290. doi: 10.3390/cancers13133290 PubMed DOI PMC

Tubin S, Vozenin M, Prezado Y, Durante M, Prise K, Lara P, et al. . Novel unconventional radiotherapy techniques: Current status and future perspectives - report from the 2nd international radiation oncology online seminar. Clin Trans Radiat Oncol. (2023) 40:100605. doi: 10.1016/j.ctro.2023.100605 PubMed DOI PMC

Ferini G, Parisi S, Lillo S, Viola A, Minutoli F, Critelli P, et al. . Impressive results after ”metabolism-guided“ lattice irradiation in patients submitted to palliative radiation therapy: Preliminary results of lattice01 multicenter study. Cancers (Basel). (2022) 14:3909. doi: 10.3390/cancers14163909 PubMed DOI PMC

Ferini G, Valenti V, Viola A, Umana G, Illari S, Parisi S, et al. . First-ever clinical experience with magnetic resonance-based lattice radiotherapy for treating bulky gynecological tumors. Anticancer Res. (2022) 42:4641. doi: 10.21873/anticanres.15968 PubMed DOI

Ferini G, Castorina P, Valenti V, Illari S, Sachpazidis I, Castorina L, et al. . A novel radiotherapeutic approach to treat bulky metastases even from cutaneous squamous cell carcinoma: Its rationale and a look at the reliability of the linear-quadratic model to explain its radiobiological effects. Front Oncol. (2022) 12:809279. doi: 10.3389/fonc.2022.809279 PubMed DOI PMC