Next Article in Journal
Non-Coding Micro RNAs and Hypoxia-Inducible Factors Are Selenium Targets for Development of a Mechanism-Based Combination Strategy in Clear-Cell Renal Cell Carcinoma—Bench-to-Bedside Therapy
Next Article in Special Issue
Clinically Relevant Radiation Exposure Differentially Impacts Forms of Cell Death in Human Cells of the Innate and Adaptive Immune System
Previous Article in Journal
Effects of an Interchain Disulfide Bond on Tropomyosin Structure: A Molecular Dynamics Study
Previous Article in Special Issue
Stereotactic Radiosurgery and Immune Checkpoint Inhibitors in the Management of Brain Metastases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Optimal Radiation Dose to Induce Robust Systemic Anti-Tumor Immunity

1
Nalecz Institute of Biocybernetics and Biomedical Engineering, Polish Academy of Sciences, Ks. Trojdena 4, 02-109 Warsaw, Poland
2
Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, 12902 Magnolia Drive, Tampa, FL 33612, USA
3
Department of Radiation Oncology, H. Lee Moffitt Cancer Center and Research Institute, 12902 Magnolia Drive, Tampa, FL 33612, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(11), 3377; https://doi.org/10.3390/ijms19113377
Submission received: 13 September 2018 / Revised: 8 October 2018 / Accepted: 23 October 2018 / Published: 29 October 2018
(This article belongs to the Special Issue Partnership of Radiotherapy and Immunotherapy)

Abstract

:
The synergy of radiation and the immune system is currently receiving significant attention in oncology as numerous studies have shown that cancer irradiation can induce strong anti-tumor immune responses. It remains unclear, however, what are the best radiation fractionation protocols to maximize the therapeutic benefits of this synergy. Here, we present a novel mathematical model that can be used to predict and dissect the complexity of the immune-mediated response at multiple tumor sites after applying focal irradiation and systemic immunotherapy. We successfully calibrate the proposed framework with published experimental data, in which two tumors were grown in mice at two spatially-separated sites from which only one was irradiated using various radiation fractionation protocols with and without concurrent systemic immunotherapy. The proposed model is calibrated to fit the temporal dynamics of tumor volume at both sites and can predict changes in immune infiltration in the non-irradiated tumors. The model was then used to investigate additional radiation fractionation protocols. Model simulations suggest that the optimal radiation doses per fraction to maximize anti-tumor immunity are between 10 and 13 Gy, at least for the experimental setting used for model calibration. This work provides the framework for evaluating radiation fractionation protocols for radiation-induced immune-mediated systemic anti-tumor responses.

1. Introduction

Radiation is a potent cancer treatment to sterilize abnormal or potentially abnormal tissues locally. About 50% of all cancer patients receive radiotherapy as part of their treatment, either as monotherapy or in combination with surgery, chemotherapy or, more recently, immunotherapy [1,2]. Despite decades of clinical success that have seen tremendous physical advances, some questions in radiobiology remain elusive including how to trigger the radiation-induced bystander [3,4,5] and the abscopal [6,7] effects reliably to further increase therapeutic efficacy. The abscopal effect, the systemic regression of metastatic nodules outside the radiation field after localized irradiation to one or more tumor sites, has attracted significant attention in the last several years, due in part to spectacular breakthroughs in immunotherapies [8,9,10]. Experimental studies conclusively demonstrated that abscopal effects are immune mediated [6,11]. The historic rarity of clinical observations of systemic responses to local radiation is conceivably due to the lack of additional boosts to the immune system [12], as well as infrequent radiological follow-up of untreated metastases. In a retrospective analysis of pre-and post-therapy radiology images of 47 metastatic melanoma cases treated with both ipilimumab and 65 courses of radiation, a shrinking of tumor lesions outside the radiation field was observed in 16 patients, the majority of which had exhibited continued growth after immunotherapy alone [13]. In a recent prospective clinical trial combining focal radiation therapy with immune stimulation in metastatic patients, 11 of 41 patients (26.8%) had an objective abscopal response [9]. The increasing body of literature and evidence that radiation can have robust systemic effects might also change the paradigm of radiation and surgery sequencing for some early-stage localized diseases [14]. Consensus has arisen that the radiation-induced abscopal effect is triggered by stimulation of the adaptive part of the immune system [6,11,15,16], and T cell-mediated responses in particular. Radiation induces cell stress and an immunogenic type of cell death, which activates and primes antigen-presenting cells (APCs) with tumor-associated antigens [17,18]. Activated APCs travel to the tumor draining lymph node where naive CD8+ T lymphocytes can be activated against tumor cells presenting these specific antigens [19]. These newly-activated lymphocytes are distributed systemically by the circulatory system, from which they can also extravasate at unirradiated tumor sites; this is one of the rate-limiting steps for triggering a radiation-induced abscopal effect [7]. Concurrent immunotherapies could augment radiation-induced immune system activation [20], such as immune checkpoint blockade with anti-CTLA-4 antibody to inhibit receptor CTLA-4 [21,22,23] that blocks T cell function. Whilst significant inroads have been made into understanding the radiation-immune system synergy, one of the clinically most relevant questions remains unanswered: What is the best radiation dose and dose fractionation to maximize immune activation? Current radiation protocols have not specifically focused on enhancing immune responses. In fact, fewer, larger doses induce significantly stronger anti-tumor immunity [9,11]. Understanding the complex, non-linear cytotoxic and immunologic consequences of RT, as well as the optimal radiation dose to induce robust anti-tumor immunity is of high biological interest and clinical value. Inroads have been made to compare different radiation schemas for specific cell lines to identify the best of the tested protocols Dewan et al. [11]; yet, limited resources are not permissive of exhaustive preclinical evaluation of all possible dose and dose fractionation protocols [24]. Here, we propose a mathematical framework informed with experimental data that may simulate experimentally untested protocols to help identify radiotherapy protocols optimized for immune activation. For demonstration purposes, but without loss of generality, we present the framework trained on published experimental data of a murine breast cancer model with two co-evolving tumor sites in an in vivo bilateral tumor model (Figure 1).
We develop a mathematical model that simulates radiation response, radiation-induced immune activation, immune checkpoint blockade therapy and the inter-exchange of activated T cells between tumor sites. The proposed model departs from earlier theoretical constructs of predator-prey systems [25] and is first compared to the experimental data from [11]. This data fitting procedure allows us to estimate the model parameters, which are then used to predict responses to dosage not considered in the experimental setting.

2. Results

2.1. Data Fitting

The optimization procedure provided the same set of estimated model parameters for about 80% of the 1000 randomly-generated initial parameter sets, indicating a robust model structure and that there are only a few local minima. The average volumetric differences between model simulation and experimental data for optimal set of parameters (Table 1) was 16 ± 19 mm3 (average relative error of 19.6 ± 16%).
Considering the inherent uncertainty of caliper measurements used in the considered experiment [11], such a low volumetric deviation indicates excellent agreement between the model and analyzed experiment. Figure 2 shows the graphical comparison of measured and simulated tumor volumes for each set of experiments (see Appendix Figure A1 and Figure A2 for a more detailed temporal dynamics visualization).
Most importantly, the model was able to simulate the synergy of radiation with the immune system correctly. In agreement with the experimental data, only combination therapy yielded a substantial response of the secondary tumor site (Figure 2 right panels). At the same time, there was only a small difference between model predicted and experimentally evaluated changes in immune infiltration of the secondary, non-irradiated tumor (Figure 3A).
Interestingly, model parameters indicated a non-monotonic dependence of the fraction of cells that will undergo immunogenic cell death ( A I ) on the radiation dose (Table 1, Figure 3B). The model fit to the data suggests that immunogenic cell death was most prevalent after three fractions of radiation with D = 8 Gy. With the derived parameter set, the tumor volume radiation survival fraction decreased with increasing radiation dose (Figure 3B).

2.2. Predicted Radiation Response

To investigate the response to various radiation fractionation protocols, we needed to interpolate both the values of survival fraction ( S F D ) and the fraction of cells undergoing immunogenic cell death ( A I D ) for radiation doses not tested experimentally. Following Sachs et al. [26], we assumed that the radiation survival fraction can be expressed as:
S F D = exp α D β G ( λ T ) D 2 ,
where G ( λ T ) = 2 ( λ T + exp ( λ T ) 1 ) / ( λ T ) 2 , λ is the repair rate, T is the delivery time, D is the dose and α and β are linear-quadratic model parameters. The above equation was able to fit model-estimated values of S F D for D = 6, 8, 20 (see Table 1) for parameter values α = 2.8 × 10 8 , β = 0.0132 and λ = 2.0358 (Figure 3B). It is worth mentioning that the parameters of the radiation response model (1) are conventionally estimated using in vitro clonogenic survival data after 10–14 days. The values reported here refer to in vivo volumetric tumor survival, and thus, the absolute values may not be directly comparable. To interpolate the non-monotonic dependence of the fraction of cells undergoing immunogenic cell death on radiation dose, we used the log-normal distribution without the restriction that the integral over the whole domain needs to be equal to one:
A I D = min γ D exp ( log ( D ) μ ) 2 ρ , 1 .
Equation (2) is able to fit estimated values of A I D for D = 6, 8 and 20 Gy for parameter values γ = 14.173, μ = 2.448 and ρ = 0.232 (Figure 3B).

2.3. Optimal Radiation Dose and Dose Fractionation

We simulated the response of both primary and secondary tumors to a single dose irradiation to the primary one and compare final overall tumor burden ( V 1 (32 days ) + V 2 ( 32 days ) ) to simulated responses to radiotherapy divided into multiple fractions with the same total or biologically effective dose (BED). We assumed that radiation is delivered on a daily basis starting from Day 12, and following [27,28], we used the standard BED formula:
BED = n D ( 1 + D / θ ) .
with θ = 20 Gy. In all cases, we simulated concurrent 9H10 immunotherapy using protocols from the experimental setup that was used to calibrate the model.
The differences in final tumor volumes dependent on radiation fractionation were primarily governed by the response of the secondary tumor as the primary tumor was almost completely eradicated for a total dose of 60 Gy independent of fractionation schedule. Model simulations suggested that the overall tumor response could vary by more than one order of magnitude depending on the radiation protocol. For a total dose of 40 Gy divided into three fractions and immunotherapy administered at Days 12, 15 and 18, the overall tumor burden at Day 32 was 12 mm3, compared to 513 mm3 if the same total dose was delivered in 15 fractions of 2.67 Gy each (Figure 4A).
Similarly, the model predicts that for 40 Gy BED, the overall tumor burden can be reduced by a factor of almost 20 for the best radiation fractionation (30.5 mm3 for 12.36 Gy × 2 compared to 552 mm3 for 2.38 Gy × 15; Figure 4B). Moreover, simulations suggest that optimal number of fractions and dose per fraction may allow up to a three-fold reduction of the total dose whilst providing a similar volumetric outcome (53.8 mm3 for 7.5 Gy × 8 fractions to a total of 60 Gy total dose vs. 65.3 mm3 for 10 Gy × 2 for a total of 20 Gy; Figure 4A), thereby potentially reducing normal tissue toxicity.
Systemic analysis of the optimal number of fractions and dose per fraction, i.e., the number of fractions that leads to the smallest overall tumor burden, for different total and biologically effective doses indicates that there is a certain range of doses per fraction to maximize therapeutic effect. For a clinically relevant total dose (Figure 4C) or biologically effective dose (Figure 4D) of 60 Gy, the model predicted the optimal treatment with 4–5 fractions of 10 Gy–13 Gy each. Model simulations further indicate that timing of immunotherapy within the tested protocols contributed significantly less to overall outcome compared to the radiation fractionation schema. However, the relatively simple model may not be appropriate to investigate timing of immunotherapy thoroughly, as the detailed pharmacokinetics and dynamics of blood absorption from the peritoneal cavity were omitted.

3. Discussion

Since the advent of immunotherapy, the synergy of radiation with the immune system has received unprecedented attention. The possibility of rationally inducing systemic anti-tumor immunity with local therapy has the flavor of the long-sought “magic bullet” for cancer cure. While clinical trials for immunotherapy in combination with radiotherapy are plentiful, very little attention is given to a priori optimization of radiation dose and dose fractionation to induce robust anti-tumor immunity. The number of clinically feasible radiation dose and dose fractionations is plentiful and an exhaustive pre-clinical evaluation impossible. Progress in integrated mathematical oncology may render such analyses feasible [24]. Here, we introduced a novel mathematical model calibrated with experimental data to make inroads into deciphering the complexity of radiation and immune system synergy. This framework estimates the optimal number of radiation fractions and radiation dose per fraction to elicit the most pronounced systemic immune-mediated tumor responses for clinically relevant total and biologically effective radiation doses.
Model simulations trained on a specific experimental study suggest that for the chosen cell line in the chosen mouse model, stereotactic radiation may have the highest likelihood of inducing robust anti-tumor immunity, which may contribute to the eradication of the tumor targeted by radiation, as well as tumors and individual cancer cells outside the radiation field. However, to draw conclusions about immune activation, longitudinal blood samples would need to be collected and analyzed for change in circulating immune populations and compared to conventional radiation doses. Standard of care radiation for many localized tumors is delivered in 2-Gy daily fractions for a total of 50–70 Gy [27]. According to model predictions, conventional radiation schemas may not be able to elicit strong immune-mediated tumor responses. It is conceivable that systemic responses would be especially important for patients with lymph node involvement, circulating tumor cells or subclinical or undiagnosed metastatic deposits [14]. This will add an additional layer of complexity to balance the different biological consequences of radiation therapy: induced maximum cell killing whilst sparing immune cells and maximizing subsequent immunity. Recent studies have shown that lymph node irradiation may hinder the maturation and infiltration of immune cells [29], yet lymph nodes are routinely irradiated to eradicate metastatic deposits. While the presented framework may provide a useful tool in preclinical research, caution is warranted before clinical translation. Nevertheless, our work indicates that the results of the recent proof-of-principle trial of 3.5 Gy × 10 fractions with concurrent immunotherapy with more than 25% of patients exhibiting objective systemic responses to focal radiation [9] would be further enhanced if 3–4 stereotactic radiation fractions were applied.
The presented results, however, may be biased by the experimental data used for model calibration. More complex metastatic disease distributions, i.e., combinations of tumors in different organs, need to be evaluated using the proposed framework [7,30]. To increase confidence in the clinical validity of framework predictions, prospective in vivo experiments should evaluate the radiation dose-dependent probability of inducing immunogenic cell death. Nevertheless, the presented work provides the first inroads and quantitative foundations to evaluating radiation fractionation protocols for inducing immune-mediated systemic anti-tumor responses. With the continuously increasing number of clinical trials combining radiation and various forms of immunotherapy, such a framework could become an invaluable tool to evaluate clinical responses, validate the model, help design subsequent clinical protocols and, eventually, support individual patient treatment planning.

4. Materials and Methods

4.1. Experimental Data

Herein, we use the data reported by Dewan and colleagues in [11], in which TSAmouse breast carcinoma cells were injected subcutaneously into syngeneic mice at two spatially-separated sites (left flank, Day 0, primary tumor; and right flank, Day 2, secondary tumor; Figure 1). Mice received one of four different treatment protocols: Group 1: no treatment; Group 2: irradiation of the primary tumor only with different fractionation protocols: 20 Gy × 1 fraction, 8 Gy × 3 and 6 Gy × 5 (daily doses for each fractionated protocol; radiation delivered at 600 cGy/min); Group 3: three doses of CTLA-4-blocking monoclonal antibody (9H10; immune checkpoint blockade) given at least two days apart, starting on Day 12, 14 or 16; Group 4: combination radiation and immunotherapy (protocol combinations as in Groups 2 and 3). Tumor volumes at both flanks were evaluated every 2–3 days until Day 32. Treatment was initiated when both tumors reached palpable sizes (Day 12; 32 mm3 and 21 mm3 on average for primary and secondary tumors, respectively). At the end of the experiment (Day 35), secondary tumors from all groups, except for those that received radiation using fractionations other than 3 × 8 Gy, were excised and analyzed by fluorescence microscopy for the presence of CD8+ T cells. The reported mean value at each experimental time point was extracted from [11] for this study.

4.2. Mathematical Model of Tumors-Immune System Interaction

We develop a mathematical model that simulates radiation response, radiation-induced immune activation, immune checkpoint blockade therapy and inter-exchange of activated T cells between tumor sites. The proposed model departs from earlier theoretical constructs of predator-prey systems [25]. Based on the experimental setting [11], we consider two tumor sites ( i = 1 , 2 ) at two spatially-separated locations each characterized by time-dependent volume V i ( t ) . At both sites, we formalize the temporal dynamics of four co-existing and interacting populations: (1) viable cancer cells (volume C i ( t ) mm3); (2) cancer cells dying in a non-immunogenic manner (volume D i ( t ) mm3); (3) cancer cells dying in an immunogenic manner (volume I i ( t ) mm3); and (4) activated tumor-specific cytotoxic T cells (effector cells; density E i ( t ) cells/mm3). Assuming that immune cells do not contribute significantly to the observed tumor volume, we denote the total measurable volume with:
V i ( t ) = C i ( t ) + D i ( t ) + I i ( t ) .
The volume occupied by viable cancer cells is assumed to follow logistic growth with site-independent tumor carrying capacity, K, and growth rate, r, modulated by the predation of immunocompetent effector cells. Here, carrying capacity and growth rate are assumed to be the same for each tumor site following the experimental setup of cell line-initiated subcutaneous tumors in contralateral flanks of the mouse. The equation governing the growth of each tumor site is:
d C i ( t ) d t = r C i ( t ) 1 V i ( t ) K a C i ( t ) E i ( t ) .
The above equation is analogous to the one considered in [31,32] and in our earlier work [7]. Herein, however, the equation considers population volumes or densities ( C i ( t ) , V i ( t ) and E i ( t ) ) instead of individual cell numbers. A detailed description of all parameters with specific values and units can be found in Table 1.
To limit model complexity, we assume that cancer cell death due to competition for resources or immune cell predation is non-immunogenic:
d D i ( t ) d t = r C i ( t ) V i ( t ) K + a C i ( t ) E i ( t ) d D i ( t ) ,
in which parameter d denotes a fixed clearance rate of dying cells.
After primary tumor ( V 1 ) irradiation, we introduce an additional volume of dying cells. Radiation-induced cell death is described using the standard approach [27,28]; at each irradiation moment t j + , the model system state is instantaneously modified as follows:
C 1 ( t j + ) = C 1 ( t j ) S F D , D 1 ( t j + ) = D 1 ( t j ) + C 1 ( t j ) ( 1 S F D ) ( 1 A I D ) , I 1 ( t j + ) = I 1 ( t j ) + C 1 ( t j ) ( 1 S F D ) A I D ,
where t j and t j + denote the times immediately before and after irradiation, respectively; S F D denotes the fraction of viable cancer cells surviving radiation with dose D; and A I D is the dose-dependent fraction of cancer cells that undergo immunogenic cell death. Consequently, 1 A I D denotes the fraction of non-immunogenic cell death events. Here, irradiation is the only source of cells in the compartment I 1 ( t ) from which they are cleared with rate d:
d I 1 ( t ) d t = d I 1 ( t ) .
We consider APCs that recognize the same tumor-specific antigens at both tumor sites and travel to the tumor-draining lymph node to activate naive T cells. Those activated T cell traffick with the blood system before extravasating at one of the two tumor sites. According to our prior work [7], each trafficking T cell may extravasate at the i-th tumor site with time-dependent probability p i ( t ) , which in the considered setting can be expressed as:
p i ( t ) = V i ( t ) V 1 ( t ) + V 2 ( t ) .
For tumors that are located in different anatomic organs, Equation (8) will include additional terms. However, this could still be simulated with the previously proposed mathematical framework [7]. Without radiation and immunotherapy, effector cells are being recruited proportionally to total volumetric tumor burden V 1 ( t ) + V 2 ( t ) . The recruitment term for each tumor site i can be expressed as w p i ( t ) V 1 ( t ) + V 2 ( t ) where parameter w is the overall recruitment rate. Explicit consideration of immunotherapy effects and immunogenic cell death after radiation yields the following equation describing the number of cytotoxic T cells:
d V i ( t ) E i ( t ) d t = l E i ( t ) V i ( t ) + 1 + u ( t ) p i ( t ) w j = 1 2 1 + w 2 I j ( t ) V j ( t ) V j ( t ) ,
where the first and second terms are spontaneous exhaustion and recruitment of immune cells, respectively. Function u ( t ) describes the effect of immunotherapy and under usual pharmacokinetic assumptions [33] can be expressed as:
u ( t ) = e 0 t c ( s ) exp c l r ( t s ) d s ,
where c ( s ) is drug administration rate at time s, c l r is the clearance rate of the considered inhibitor and parameter e describes the drug impact on the system. Model Equations (3)–(10) have been implemented in MATLAB software and solved using built-in ordinary differential equation solver ode45 with lowered default tolerances (RelToland AbsTolequal to 10 8 ) to increase solution accuracy.

4.3. Parameter Estimation

The goal of the parameter estimation procedure was to identify model parameters for which the error of: (1) the simulated tumor volumes V 1 ( t ) and V 2 ( t ) compared to experimental volume measurements; and (2) simulated relative change of T cell infiltration in the secondary tumor compared to experimental data,
E R R ( p ) = e W 1 j = 1 2 t i V j ( t i , p ) V j , d a t a ( t i ) V j , d a t a ( t i ) 2 + e W 2 E 2 ( 35 , p ) E 2 , c t r l ( 35 , p ) R d a t a 2 ,
is minimized. The given set of parameters values is denoted by p; W 1 denotes the set of all experiments; W 2 denotes the set of experiments in which immune infiltration in the secondary tumor was evaluated at Day 35; t i denotes one of the recorded measurement times; c t r l denotes control experiment (no radiation and immunotherapy); and R d a t a denotes the measured change in immune infiltration relative to the control experiment.
Although 9H10 immunotherapy was injected intraperitoneally in the considered experiment [11], we assume for simplicity that c ( s ) = D ( δ ( s t ¯ 1 ) + δ ( s t ¯ 2 ) + . . . ) , where δ is the Dirac delta function and t i is injection day (we take e = e D in Equation (10)). This assumption may in general lead to a decrease of fit quality and introduces inherent bias towards applying immunotherapy at the same time as radiation, but at the same time, we do not need to estimate additional parameters related to drug release from the intraperitoneal cavity.
For parameter estimation, 124 data points were used to calibrate the values of 13 parameters simultaneously. Only values of two parameters were kept fixed during the minimization process: (1) tumor growth rate r, which was calculated directly from the initial growth data of untreated control experiments; (2) T cell exhaustion rate l, which was taken directly from [34]. For numerical error minimization, we used the deterministic derivative-based trust region method [35] implemented in the lsqnonlin function from MATLAB’s Optimization Toolbox. In order to avoid local minima, we generated 1000 random initial parameter values for the optimization procedure.

Author Contributions

Conceptualization, J.P. and H.E.; methodology, J.P.; software, J.P.; validation, J.P. and H.E.; formal analysis, J.P.; investigation, J.P.; resources, H.E.; data curation, J.P.; writing, original draft preparation, J.P.; writing, review and editing, H.E.; visualization, J.P.; supervision, H.E.; project administration, H.E.; funding acquisition, H.E.

Funding

This project was supported (in part) by the American Cancer Society’s Institutional Research Grant.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
RTRadiotherapy
APCAntigen-presenting cell
CTLCytotoxic T lymphocyte
IRIrradiation
BEDBiologically effective dose

Appendix A. Supplementary Data

Figure A1. Comparison of model simulated growth dynamics with experimental data. Solutions to the proposed model (lines; Equations (3)–(10)) were obtained after performing data fitting (estimated parameters are presented in Table 1). Shown are experimental and simulated tumor volumes for tumors grown without treatment or after monotherapy. Experimental data (circles) derived from Dewan et al. [11].
Figure A1. Comparison of model simulated growth dynamics with experimental data. Solutions to the proposed model (lines; Equations (3)–(10)) were obtained after performing data fitting (estimated parameters are presented in Table 1). Shown are experimental and simulated tumor volumes for tumors grown without treatment or after monotherapy. Experimental data (circles) derived from Dewan et al. [11].
Ijms 19 03377 g0a1
Figure A2. Comparison of model simulated growth dynamics with experimental data. Solutions to the proposed model (lines; Equations (3)–(10)) were obtained after performing data fitting (estimated parameters are presented in Table 1). Shown are experimental and simulated tumor volumes after combination treatment. Experimental data (circles) derived from Dewan et al. [11].
Figure A2. Comparison of model simulated growth dynamics with experimental data. Solutions to the proposed model (lines; Equations (3)–(10)) were obtained after performing data fitting (estimated parameters are presented in Table 1). Shown are experimental and simulated tumor volumes after combination treatment. Experimental data (circles) derived from Dewan et al. [11].
Ijms 19 03377 g0a2

References

  1. Delaney, G.; Jacob, S.; Featherstone, C.; Barton, M. The role of radiotherapy in cancer treatment: Estimating optimal utilization from a review of evidence-based clinical guidelines. Cancer 2005, 104, 1129–1137. [Google Scholar] [CrossRef] [PubMed]
  2. Erridge, S.C.; Featherstone, C.; Chalmers, R.; Campbell, J.; Stockton, D.; Black, R. What will be the radiotherapy machine capacity required for optimal delivery of radiotherapy in Scotland in 2015? Eur. J. Cancer 2007, 43, 1802–1809. [Google Scholar] [CrossRef] [PubMed]
  3. Wideł, M.; Lalik, A.; Krzywoń, A.; Poleszczuk, J.; Fujarewicz, K.; Rzeszowska-Wolny, J. The different radiation response and radiation-induced bystander effects in colorectal carcinoma cells differing in p53 status. Mutat. Res.-Fundam. Mol. Mech. Mutagen. 2015, 778, 61–70. [Google Scholar] [CrossRef] [PubMed]
  4. Poleszczuk, J.; Krzywoń, A.; Foryś, U.; Wideł, M. Connecting radiation-induced bystander effects and senescence to improve radiation response prediction. Radiat. Res. 2015, 183, 571–577. [Google Scholar] [CrossRef] [PubMed]
  5. Blyth, B.J.; Sykes, P.J. Radiation-Induced Bystander Effects: What Are They, and How Relevant Are They to Human Radiation Exposures? Radiat. Res. 2011, 176, 139–157. [Google Scholar] [CrossRef] [PubMed]
  6. Demaria, S.; Ng, B.; Devitt, M.L.; Babb, J.S.; Kawashima, N.; Liebes, L.; Formenti, S.C. Ionizing radiation inhibition of distant untreated tumors (abscopal effect) is immune mediated. Int. J. Radiat. Oncol. Biol. Phys. 2004, 58, 862–870. [Google Scholar] [CrossRef] [PubMed]
  7. Poleszczuk, J.; Luddy, K.A.; Prokopiou, S.; Robertson-Tessi, M.; Moros, E.G.; Fishman, M.; Djeu, J.Y.; Finkelstein, S.E.; Enderling, H. Abscopal benefits of localized radiotherapy depend on activated T cell trafficing and distribution between metastatic lesions. Cancer Res. 2016, 76, 1009–1018. [Google Scholar] [CrossRef] [PubMed]
  8. Formenti, S.C.; Demaria, S. Combining radiotherapy and cancer immunotherapy: A paradigm shift. J. Natl. Cancer Inst. 2013, 105, 256–265. [Google Scholar] [CrossRef] [PubMed]
  9. Golden, E.B.; Chhabra, A.; Chachoua, A.; Adams, S.; Donach, M.; Fenton-Kerimian, M.; Friedman, K.; Ponzo, F.; Babb, J.S.; Goldberg, J.; et al. Local radiotherapy and granulocyte-macrophage colony-stimulating factor to generate abscopal responses in patients with metastatic solid tumours: A proof-of-principle trial. Lancet Oncol. 2015, 16, 795–803. [Google Scholar] [CrossRef]
  10. Crittenden, M.; Kohrt, H.; Levy, R.; Jones, J.; Camphausen, K.; Dicker, A.; Demaria, S.; Formenti, S. Current clinical trials testing combinations of immunotherapy and radiation. Semin. Radiat. Oncol. 2015, 25, 54–64. [Google Scholar] [CrossRef] [PubMed]
  11. Dewan, M.Z.; Galloway, A.E.; Kawashima, N.; Dewyngaert, J.K.; Babb, J.S.; Formenti, S.C.; Demaria, S. Fractionated but not single-dose radiotherapy induces an immune-mediated abscopal effect when combined with anti-CTLA-4 antibody. Clin. Cancer Res. 2009, 15, 5379–5388. [Google Scholar] [CrossRef] [PubMed]
  12. Kaminski, J.M.; Shinohara, E.; Summers, J.B.; Niermann, K.J.; Morimoto, A.; Brousal, J. The controversial abscopal effect. Cancer Treat. Rev. 2005, 31, 159–172. [Google Scholar] [CrossRef] [PubMed]
  13. Chandra, R.A.; Wilhite, T.J.; Balboni, T.A.; Alexander, B.M.; Spektor, A.; Ott, P.A.; Ng, A.K.; Hodi, F.S.; Schoenfeld, J.D. A systematic evaluation of abscopal responses following radiotherapy in patients with metastatic melanoma treated with ipilimumab. Oncoimmunology 2015, 4, e1046028. [Google Scholar] [CrossRef] [PubMed]
  14. Poleszczuk, J.; Luddy, K.; Chen, L.; Lee, J.K.; Harrison, L.B.; Czerniecki, B.J.; Soliman, H.; Enderling, H. Neoadjuvant radiotherapy of early-stage breast cancer and long-term disease-free survival. Breast Cancer Res. 2017, 19, 75. [Google Scholar] [CrossRef] [PubMed]
  15. Pilones, K.A.; Kawashima, N.; Yang, A.M.; Babb, J.S.; Formenti, S.C.; Demaria, S. Invariant natural killer T cells regulate breast cancer response to radiation and CTLA-4 blockade. Clin. Cancer Res. 2009, 15, 597–606. [Google Scholar] [CrossRef] [PubMed]
  16. Habets, T.H.; Oth, T.; Houben, A.W.; Huijskens, M.J.; Senden-Gijsbers, B.L.; Schnijderberg, M.C.; Brans, B.; Dubois, L.J.; Lambin, P.; De Saint-Hubert, M.; et al. Fractionated radiotherapy with 3 × 8 Gy induces systemic anti-tumour responses and abscopal tumour inhibition without modulating the humoral anti-tumour response. PLoS ONE 2016, 11, e0159515. [Google Scholar] [CrossRef] [PubMed]
  17. Kroemer, G.; Galluzzi, L.; Kepp, O.; Zitvogel, L. Immunogenic cell death in cancer therapy. Annu. Rev. Immunol. 2013, 31, 51–72. [Google Scholar] [CrossRef] [PubMed]
  18. Reits, E.A.; Hodge, J.W.; Herberts, C.A.; Groothuis, T.A.; Chakraborty, M.; Wansley, E.K.; Camphausen, K.; Luiten, R.M.; de Ru, A.H.; Neijssen, J.; et al. Radiation modulates the peptide repertoire, enhances MHC class I expression, and induces successful anti-tumor immunotherapy. J. Exp. Med. 2006, 203, 1259–1271. [Google Scholar] [CrossRef] [PubMed]
  19. Coussens, L.M.; Werb, Z. Inflammation and cancer. Nature 2002, 420, 860–867. [Google Scholar] [CrossRef] [PubMed]
  20. Vatner, R.E.; Cooper, B.T.; Vanpouille-Box, C.; Demaria, S.; Formenti, S.C. Combinations of immunotherapy and radiation in cancer therapy. Front. Oncol. 2014, 4, 325. [Google Scholar] [CrossRef] [PubMed]
  21. Postow, M.A.; Callahan, M.K.; Wolchok, J.D. Immune checkpoint blockade in cancer therapy. J. Clin. Oncol. 2014, 33, 1974–1982. [Google Scholar] [CrossRef] [PubMed]
  22. Callahan, M.K.; Wolchok, J.D.; Allison, J.P. AntiCTLA-4 antibody therapy: Immune monitoring during clinical development of a novel immunotherapy. Semin. Oncol. 2010, 37, 473–484. [Google Scholar] [CrossRef] [PubMed]
  23. Wolchok, J.D.; Saenger, Y. The Mechanism of Anti-CTLA-4 Activity and the Negative Regulation of T-Cell Activation. Oncologist 2008, 13, 2–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Enderling, H.; Kim, S.; Pilon-Thomas, S. The Accelerating Quest For Optimal Radiation And Immunotherapy Combinations For Local And Systemic Tumor Control. Ther. Radiol. Oncol. 2018, 2, 33. [Google Scholar] [CrossRef]
  25. Walker, R.; Enderling, H. From concept to clinic: Mathematically informed immunotherapy. Curr. Probl. Cancer 2016, 40, 68–83. [Google Scholar] [CrossRef] [PubMed]
  26. Sachs, R.K.; Hahnfeld, P.; Brenner, D.J. The link between low-LET dose-response relations and the underlying kinetics of damage production/repair/misrepair. Int. J. Radiat. Biol. 1997, 72, 351–374. [Google Scholar] [PubMed]
  27. Prokopiou, S.; Moros, E.G.; Poleszczuk, J.; Caudell, J.; Torres-Roca, J.F.; Latifi, K.; Lee, J.K.; Myerson, R.; Harrison, L.B.; Enderling, H. A proliferation saturation index to predict radiation response and personalize radiotherapy fractionation. Radiat. Oncol. 2015, 10, 159. [Google Scholar] [CrossRef] [PubMed]
  28. Poleszczuk, J.; Walker, R.; Moros, E.G.; Latifi, K.; Caudell, J.J.; Enderling, H. Predicting Patient-Specific Radiotherapy Protocols Based on Mathematical Model Choice for Proliferation Saturation Index. Bull. Math. Biol. 2018, 80, 1195–1206. [Google Scholar] [CrossRef] [PubMed]
  29. Zhang, Z.; Niedermann, G. Abscopal Effects With Hypofractionated Schedules Extending Into The Effector Phase Of The Tumor-Specific T-Cell Response. Int. J. Radiat. Onclol. Biol. Phys. 2018, 101, 63–73. [Google Scholar] [CrossRef] [PubMed]
  30. Walker, R.; Poleszczuk, J.; Pilon-Thomas, S.; Anderson, A.; Czerniecki, B.; Harrison, L.; Moros, E.; 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] [PubMed]
  31. 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] [PubMed] [Green Version]
  32. Kuznetsov, V.A.; Knott, G.D. Modeling tumor regrowth and immunotherapy. Math. Comput. Model. 2001, 33, 1275–1287. [Google Scholar] [CrossRef]
  33. Bonate, P.L. Pharmacokinetic-Pharmacodynamic Modeling and Simulation; Springer: Boston, MA, USA, 2011. [Google Scholar]
  34. Kirchner, D.; Panetta, J.C. Modeling immunotherapy of the tumor-immune interaction. J. Math. Biol. 1998, 37, 235–252. [Google Scholar] [CrossRef]
  35. Byrd, R.H.; Schnabel, R.B.; Schultz, G.A. A trust region algorithm for nonlinearly constrained optimization. SIAM J. Numer. Anal. 1987, 24, 1151–1170. [Google Scholar] [CrossRef]
Figure 1. Schematic of the experimental setting and model compartments (circles). Arrows indicate transitions between compartments and T bars indicate killing of cancer cells by immune cells. The proposed model (Equations (3)–(10)) formalizes the relevant mechanisms underlying the in vivo experiment [11] in which carcinoma cells were injected subcutaneously into mice at two spatially-separated sites (left and right flank), and only one site (1, primary tumor) was later irradiated (IR). Some mice received also systemic immunotherapy (9H10), which augments activation of tumor-specific cytotoxic T lymphocytes (CTLs).
Figure 1. Schematic of the experimental setting and model compartments (circles). Arrows indicate transitions between compartments and T bars indicate killing of cancer cells by immune cells. The proposed model (Equations (3)–(10)) formalizes the relevant mechanisms underlying the in vivo experiment [11] in which carcinoma cells were injected subcutaneously into mice at two spatially-separated sites (left and right flank), and only one site (1, primary tumor) was later irradiated (IR). Some mice received also systemic immunotherapy (9H10), which augments activation of tumor-specific cytotoxic T lymphocytes (CTLs).
Ijms 19 03377 g001
Figure 2. Comparison of the model simulated growth dynamics with experimental data. Solutions to the proposed model (Equations (3)–(10)) were obtained after performing data fitting (estimated parameters are presented in Table 1). (A) Experimental and simulated tumor volumes for tumors grown without treatment or after monotherapy. (B) Experimental and simulated tumor volumes after combination treatment. Experimental data (red circles) derived from Dewan et al. [11].
Figure 2. Comparison of the model simulated growth dynamics with experimental data. Solutions to the proposed model (Equations (3)–(10)) were obtained after performing data fitting (estimated parameters are presented in Table 1). (A) Experimental and simulated tumor volumes for tumors grown without treatment or after monotherapy. (B) Experimental and simulated tumor volumes after combination treatment. Experimental data (red circles) derived from Dewan et al. [11].
Ijms 19 03377 g002
Figure 3. Comparison of relative T cell density and radiation response curves. (A) Model simulated and experimentally measured changes in immune cell infiltration at the secondary tumor site after the combination of radiotherapy and 9H10 immunotherapy. (B) Proposed interpolation of the model estimated radiation survival fraction ( S F D ) and proportion of cells undergoing immunogenic cell death ( A I D ) for radiation doses D = 6, 8 and 20 Gy (see Table 1) using Equations (1) and (2).
Figure 3. Comparison of relative T cell density and radiation response curves. (A) Model simulated and experimentally measured changes in immune cell infiltration at the secondary tumor site after the combination of radiotherapy and 9H10 immunotherapy. (B) Proposed interpolation of the model estimated radiation survival fraction ( S F D ) and proportion of cells undergoing immunogenic cell death ( A I D ) for radiation doses D = 6, 8 and 20 Gy (see Table 1) using Equations (1) and (2).
Ijms 19 03377 g003
Figure 4. Optimal radiation fractionation and dose per fraction for immune activation. Dependence of the model predicted overall tumor burden at Day 32, i.e., V 1 ( 32 days ) + V 2 ( 32 days ) , for different total (A) and biologically effective doses (B) for various numbers of fractions. Radiation is delivered daily, and three concurrent doses of 9H10 immunotherapy are applied at different times [11]. (C,D) Optimal number of radiation fractions and corresponding doses per fraction depending on the prescribed total and biologically effective doses. We assume that 9H10 immunotherapy is applied on Days 12,15 and 18.
Figure 4. Optimal radiation fractionation and dose per fraction for immune activation. Dependence of the model predicted overall tumor burden at Day 32, i.e., V 1 ( 32 days ) + V 2 ( 32 days ) , for different total (A) and biologically effective doses (B) for various numbers of fractions. Radiation is delivered daily, and three concurrent doses of 9H10 immunotherapy are applied at different times [11]. (C,D) Optimal number of radiation fractions and corresponding doses per fraction depending on the prescribed total and biologically effective doses. We assume that 9H10 immunotherapy is applied on Days 12,15 and 18.
Ijms 19 03377 g004
Table 1. Model (Equations (3)–(10)) parameters and derived values.
Table 1. Model (Equations (3)–(10)) parameters and derived values.
ParameterDescriptionUnitValue
rViable cancer cells volume doubling time1/day0.195
KTumor carrying capacitymm31423.1
aCTLs’ killing ratemm3/(cell·day)0.0177
dClearance rate of dying cells1/day0.264
S F D Fraction of viable cancer cells that survive
after radiation dose
D = 20 Gy-0.265
D = 8 Gy-0.664
D = 6 Gy-0.783
A I D Fraction of cells that will undergo
immunogenic cell death after
radiation dose
D = 20 Gy-0.194
D = 8 Gy-0.984
D = 6 Gy-0.367
lDecay rate of effector cells1/day0.03
wBaseline T cell recruitment ratecell/ (mm3·day)0.135
w 2 Fold change in the baseline T cell recruitment
rate due to immunogenic cell death-15.37
eInitial fold change in recruitment
of cytotoxic T cells caused by the administered
dose of 9H10 immunotherapy-8.495
c l r 9H10 immunotherapy clearance rate1/day0.967

Share and Cite

MDPI and ACS Style

Poleszczuk, J.; Enderling, H. The Optimal Radiation Dose to Induce Robust Systemic Anti-Tumor Immunity. Int. J. Mol. Sci. 2018, 19, 3377. https://doi.org/10.3390/ijms19113377

AMA Style

Poleszczuk J, Enderling H. The Optimal Radiation Dose to Induce Robust Systemic Anti-Tumor Immunity. International Journal of Molecular Sciences. 2018; 19(11):3377. https://doi.org/10.3390/ijms19113377

Chicago/Turabian Style

Poleszczuk, Jan, and Heiko Enderling. 2018. "The Optimal Radiation Dose to Induce Robust Systemic Anti-Tumor Immunity" International Journal of Molecular Sciences 19, no. 11: 3377. https://doi.org/10.3390/ijms19113377

APA Style

Poleszczuk, J., & Enderling, H. (2018). The Optimal Radiation Dose to Induce Robust Systemic Anti-Tumor Immunity. International Journal of Molecular Sciences, 19(11), 3377. https://doi.org/10.3390/ijms19113377

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