Next Article in Journal
Treatment and Prognosis of COVID-19 Associated Olfactory and Gustatory Dysfunctions
Next Article in Special Issue
Dynamics-Adapted Radiotherapy Dose (DARD) for Head and Neck Cancer Radiotherapy Dose Personalization
Previous Article in Journal
Validation of a Personalized Hearing Screening Mobile Health Application for Persons with Moderate Hearing Impairment
Previous Article in Special Issue
The Effect of Radiotherapy on Diffuse Low-Grade Gliomas Evolution: Confronting Theory with Clinical Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Combinations of Chemotherapy and Radiotherapy in Low-Grade Gliomas: A Mathematical Approach

by
Luis E. Ayala-Hernández
1,*,
Armando Gallegos
1,
Philippe Schucht
2,
Michael Murek
2,
Luis Pérez-Romasanta
3,
Juan Belmonte-Beitia
4 and
Víctor M. Pérez-García
4
1
Departamento de Ciencias Exactas y Tecnología Centro Universitario de los Lagos, Universidad de Guadalajara, Enrique Díaz de León 1144, Colonia Paseos de la Montaña, Lagos de Moreno C.P. 47460, Jalisco, Mexico
2
Universitätsklinik für Neurochirurgie, Bern University Hospital, CH-3010 Bern, Switzerland
3
Radiation Oncology Service, Salamanca University Hospital, 37007 Salamanca, Spain
4
Department of Mathematics, Mathematical Oncology Laboratory (MOLAB), Universidad de Castilla-La Mancha, Avda. Camilo José Cela, 3, 13071 Ciudad Real, Spain
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2021, 11(10), 1036; https://doi.org/10.3390/jpm11101036
Submission received: 18 August 2021 / Revised: 30 September 2021 / Accepted: 11 October 2021 / Published: 16 October 2021
(This article belongs to the Special Issue Mathematical Applications for Clinical Radiotherapy)

Abstract

:
Low-grade gliomas (LGGs) are brain tumors characterized by their slow growth and infiltrative nature. Treatment options for these tumors are surgery, radiation therapy and chemotherapy. The optimal use of radiation therapy and chemotherapy is still under study. In this paper, we construct a mathematical model of LGG response to combinations of chemotherapy, specifically to the alkylating agent temozolomide and radiation therapy. Patient-specific parameters were obtained from longitudinal imaging data of the response of real LGG patients. Computer simulations showed that concurrent cycles of radiation therapy and temozolomide could provide the best therapeutic efficacy in-silico for the patients included in the study. The patient cohort was extended computationally to a set of 3000 virtual patients. This virtual cohort was subject to an in-silico trial in which matching the doses of radiotherapy to those of temozolomide in the first five days of each cycle improved overall survival over concomitant radio-chemotherapy according to RTOG 0424. Thus, the proposed treatment schedule could be investigated in a clinical setting to improve combination treatments in LGGs with substantial survival benefits.

1. Introduction

Adult supratentorial WHO grade II diffuse low-grade gliomas (LGGs) are slowgrowing primary brain tumors that are, in general, incurable due to their infiltrative nature. However, active surgical and oncological treatment may lead to patient survival exceeding 10 years following diagnosis [1].
Treatment typically consists of surgery followed by observation, radiotherapy, chemotherapy, or chemoradiation [2]. The decision as to the specific combination of therapies to be used on each patient is based on the qualitative consideration of different variables including age, tumor grade, performance status, tumor histology and location [3].
Systemic treatments play an important role in the management of LGGs. Procarbazine, lomustine, and vincristine (PCV) and temozolomide (TMZ) have been shown to be effective against low grade gliomas [4,5,6]. Today, both therapies are widely recommended by recent clinical practice guidelines [7,8]. Due to its favorable toxicity profile, prolonged TMZ treatment is a relevant option either as upfront or as adjuvant treatment [9].
With regard to combination treatments, Radiation Therapy Oncology Group (RTOG) 9802, a prospective clinical trial, showed improvements in both progression-free survival (PFS) and overall survival (OS) with six cycles of adjuvant procarbazine, lomustine, and vincristine (PCV) chemotherapy following radiation, when compared with RT alone [10].
There is evidence that the combination of chemotherapy and radiotherapy could be a beneficial strategy for the treatment of LGGs [10]. However, large-scale prospective studies are lacking concerning the combination of radiotherapy with TMZ. In high-grade gliomas (HGG), TMZ used as concurrent and adjuvant chemotherapy has become the standard of care [11]. However, in LGGs, further studies are necessary to ascertain the role of TMZ when used in addition to RT. RTOG 0424, a single-arm phase II study, combined concurrent and adjuvant TMZ with RT to treat patients with high-risk LGG [12]. The 3-year OS rate was 73.1%, which exceeded historical controls. Randomized clinical trials are under way to confirm the advantage of TMZ chemotherapy added to RT for patients with high-risk LGGs [13].
Thus, it seems that the combination of chemotherapy and radiotherapy improves PFS as against chemotherapy or radiotherapy alone for LGGs. The question of treatment schedule then arises: Is a sequential or concurrent treatment the best post-surgery therapeutic option in high-risk gliomas? Clinical trials studying combination treatments of RT and TMZ have adopted the intensive treatment pattern found to be useful in HGGs, consisting of concurrent and adjuvant chemotherapy [11]. However, the natural history of the disease is very different for HGGs and LGGs.
Because of the high cost, human effort, and the time and ethical issues involved in clinical trials and basic medical research, mathematical modeling and analysis can potentially contribute to the discovery of better cancer treatment delivery regimens. This is specially true for combination cancer therapies. Due to the many possible combinations, obtaining an optimal toxicity–efficacy balance is an increasingly complex task. Consequently, standard empirical approaches to optimizing drug dosing and scheduling in patients are now of limited utility. Mathematical modeling can substantially advance this practice through improved rationalization of therapeutic strategies [14,15].
A number of mathematical studies have constructed in-silico models to determine the optimal delivery scheme for either radiation therapy [16,17,18,19,20,21] or chemotherapy [22,23,24,25,26,27] alone. However, although mathematical models have potential for the study of optimal combination treatments [28], no studies have addressed computationally the question of what might be the best combination scheme of TMZ and RT for LGGs.
In this paper, we construct a mathematical model of LGG response to TMZ and RT, describing the longitudinal tumor volumetric dynamics using patient data. We find the parameter regimes best describing the response to treatments. We then use the mathematical model to carry out “in-silico” experiments to explore different treatment regimes. Our results show that concurrent cycles of TMZ and RT could provide substantial survival improvements over concomitant radio-chemotherapy according to RTOG 0424 [12].

2. Materials and Methods

2.1. Patient Population

Our initial patient dataset consisted of 84 biopsy/surgery confirmed LGG patients, followed by magnetic resonance imaging (MRI) at Bern University Hospital between 1990 and 2013. From that dataset, we chose patients satisfying the following inclusion criteria: (i) tumors treated with TMZ and/or RT; (ii) availability of at least 1 MRI before the start of treatment; (iii) availability of at least one MRI after the end of the treatment; (iv) no treatment other than RT and/or TMZ given in the period of study; (v) no malignant transformation in the period of study; (vi) Ki-67 labeling index 6 % when available. A total of 17 patients satisfied the inclusion criteria and were classified according to the 2016 WHO criteria as having either astrocytoma or oligodendroglioma (Table 1).
As outlined in Table 1, nine of the patients considered were classified as having oligodendroglioma, six as having astrocytoma, one with oligoastrocytoma, and one as histology not available. Treatment received was TMZ for eight of these, RT for seven, concurrent RT and TMZ for one and sequential RT and TMZ for another one. Patient ages ranged from 28 to 58 years and the Ki-67 labeling index on diagnosis ranged from 1% to 6%.

2.2. Image Acquisition and Radiological Measurements of Tumor Size

Radiological glioma growth was quantified by measurements of tumor diameters on successive MRIs (T2/FLAIR sequences). We computed the tumor volumes using the ellipsoidal approximation. The three largest tumor diameters ( D 1 , D 2 , D 3 ) were measured according to three orthogonal reference planes (axial, coronal and sagittal), and the tumor volumes were estimated using, as stated, the ellipsoidal approximation: V = ( D 1 · D 2 · D 3 ) / 2 , following standard clinical practice [29,30]. To estimate the error of the methodology, we took a different set of glioma patients from another study [31] and compared their volumes, as computed accurately using a semi-automatic segmentation approach, with those computed using the ellipsoidal approximation. Mean differences were 18%, which was the reference level taken for the error in the volume computations.

2.3. Mathematical Model

To provide a mathematical description of the growth of LGGs and how they respond to RT and TMZ, we assumed that in the absence of treatment the tumor consists of proliferative P ( t ) and quiescent Q ( t ) cells. Before the malignant transformation of the tumor, it will keep a low genetic heterogeneity, proliferative cells would grow exponentially, and there would be a bidirectional flow from Q ( t ) to P ( t ) , accounting for the cells’ ability to switch from a dormant state to the cell cycle and back [32]. These assumptions are not expected to hold for high-grade gliomas, which are known to have a high level of genetic/phenotype heterogeneity and thus different compartments corresponding to subpopulations with biological behavior should be incorporated to the model. We assumed that no additional mutations appear in LGGs for long periods of time, leading to genotypes with different biological behaviors. When TMZ or RT are administered, they damage only proliferative cells [33,34]. TMZ produces early death, accounting for necrosis, autophagy and druginduced apoptosis, with rate α 1 , and delayed death through mitotic catastrophe, with rate α 2 [26], the latter resulting in a population of lethally damaged cells D ( t ) that die with rate β 3 . RT at the doses used for LGG treatment (1.8–2.2 Gy) kills cells mostly through mitotic catastrophe. Thus, it was assumed that each radiation dose d in the RT course would move a fraction 1 S p ( d ) of P ( t ) cells into the compartment of damaged cells D ( t ) .
Recent work has considered the alternative of adding a delayed death term in the model equations for tumor cells, rather than including the compartment D ( t ) explicitly in the model [35]. We expect that approach, which does not reduce the number of free parameters, to be equivalent to the one used in this paper.
There is a key biological process that has been omitted in previous mathematical models describing LGG response to RT: the fact that RT stimulates proliferation of quiescent cells [36,37,38]. To account for this process, we will consider that a fraction 1 S q ( d ) of quiescent cells are stimulated to re-enter the cell cycle for each radiation dose d, thus being transferred into the compartment of proliferative cells P ( t ) . Studies regarding repopulation during chemotherapy are scarce, and the mechanisms behind the process are still not clearly understood [39,40], but it seems that there is no similar process in the action of TMZ, or at least that it is not quantitatively as relevant as with RT.
Finally, C ( t ) describes the concentration of TMZ in brain and tumor tissue. In this paper, we will assume a simple pharmacokinetics and consider that TMZ is cleared at a constant rate λ , leading to an exponential decay of the concentration.
Figure 1 displays a diagram of the biological processes accounted for in our mathematical model.
The equations of the mathematical model are given by the following system of ordinary differential equations:
d P d t = ρ P β 1 P + β 2 Q α 1 P C α 2 P C ,
d Q d t = β 1 P β 2 Q ,
d D d t = α 1 P C β 3 D ,
d C d t = λ C .
We assumed the number of cells to be proportional to the volumes occupied by each tumor subpopulation, and thus P ( t ) , Q ( t ) and D ( t ) were measured in cm 3 . In this way, total tumor mass P + Q + D was considered to be proportional to total tumor volume. Since the Ki-67 labeling index provides a measure of the fraction of cells proliferating in tumor tissue, we assumed it to be given by the dimensionless ratio P / ( P + Q ) .
For the particular (but important) case for which there is no treatment, that is,  C ( t ) 0 , we obtain the system:
d P d t = ρ P β 1 P + β 2 Q ,
d Q d t = β 1 P β 2 Q .
Exact solutions for this system of linear differential equations can be found in the Supplementary Information (SI) Section S1.
The initial time, corresponding to the first MRI study, was denoted by t 0 . Initial conditions for Equation (1) were taken to be P 0 = P ( t 0 ) , Q 0 = Q ( t 0 ) , D ( t 0 ) = C ( t 0 ) = 0 . Let f ( t j ) = lim t t j f ( t ) and let C * be the fraction of each chemotherapy dose reaching brain tumor tissue. Drug administration was described as impulses for the times t j , so that
P ( t j ) = P ( t j ) ,
Q ( t j ) = Q ( t j ) ,
D ( t j ) = D ( t j ) ,
C ( t j ) = C ( t j ) + C * .
Radiation therapy fractions was simulated as impulses as in previous studies [17,19]. We assumed n radiation doses d given at times t 1 , t 2 , , t n , so that after each irradiation event at time t k we get
P ( t k ) = S p ( d ) P ( t k ) + [ 1 S q ( d ) ] Q ( t k ) ,
Q ( t k ) = S q ( d ) Q ( t k ) ,
D ( t k ) = D ( t k ) + [ 1 S p ( d ) ] P ( t k ) ,
C ( t k ) = C ( t k ) .
When both therapies are administered simultaneously at times t k , Equation (4) become
P ( t k ) = S p ( d ) P ( t k ) + [ 1 S q ( d ) ] Q ( t k ) ,
Q ( t k ) = S q ( d ) Q ( t k ) ,
D ( t k ) = D ( t k ) + [ 1 S p ( d ) ] P ( t k ) ,
C ( t k ) = C ( t k ) + C * .
The cell–kill parameter S p accounts for all radiation damage that cannot be repaired and will eventually lead to cell death in a delayed form through the population D ( t ) .

2.4. Parameter Estimation

2.4.1. TMZ Concentration in the Brain

The usual chemotherapy schedule consists of cycles of 28 days, with five TMZ oral doses on days 1 to 5 and a break on days 6 to 28. The standard dose per day is c = 150 mg per m 2 of patient body surface, which is usually around 1.6 m 2 for women and 1.9 m 2 for men [41]. To further allow for drug loss during transport to the brain we calculated the value of the dose reaching the tumor as
C * = γ · c · b ,
where γ is the fraction of TMZ getting to 1 cm 3 of brain interstitial fluid (from a unit dose) and b is the patient’s body surface. γ can be calculated using the value of maximal TMZ concentration C m a x for a dose of 150 mg/m 2 taken from the literature [42]. The time to reach peak drug concentration in the brain is less than two hours and thus negligible in comparison with the other time scales in the model. We chose to set the drug concentration reached by each standard dose c to the value C m a x = 0.6 μ g/cm 3 as in Ref. [25].
The dose of TMZ given according to RTOG 0424 [12] during the concurrent stage is 75 mg/m 2 /day [11]. Thus, in that treatment stage we set the drug concentration in brain tumor tissue equal to C m a x / 2 .
Using the value of TMZ half-life clearance time t 1 / 2 [43], and following the same methodology as in Ref. [25] we determined the rate of drug decay
λ = ln ( 1 / 2 ) t 1 / 2 .
A summary of the parameter values involved in the calculation of the TMZ pharmacokinetics is presented in Table 2.

2.4.2. Parameter Fitting

Following the conventional RT administration schedule, consisting of 30 fractions of 1.8 Gy given from Monday to Friday for six consecutive weeks [21], we fixed the radiation dose d = 1.8 Gy. The parameters ρ , β 1 , β 2 , β 3 , α 1 , α 2 , S p (1.8 Gy) and S q (1.8 Gy) and the initial populations P 0 and Q 0 were fitted for each patient’s volumetric data using the Matlab (R2020b, The MathWorks, Inc., Natick, MA, USA) function fmincon.
Parameter fitting was performed for each patient by minimizing the quadratic error between longitudinal volumetric data ( t j , v j ) and simulation results ( t j , V ( t j ) ) obtained from Equations (1) and either (3), (4) or (5), that is,
e ( ρ , β 1 , β 2 , β 3 , α 1 , α 2 , S p , S q ) = i = 1 r v j V ( t j ) 2 r .
At surgery time T, at which the Ki-67% is obtained, we imposed the condition:
100 P ( T ) P ( T ) + Q ( T ) Ki 67 % = 0 .
Finally, for the initial volume, we required that:
P 0 + Q 0 [ 0.82 v 0 , 1.18 v 0 ] cm 3 .
Equation (8) were used to obtain the parameter values. Table 3 displays the results obtained for each patient in our dataset. Time between RT and TMZ for patient 151 was very long (11 years), and distinct values of the Ki-67 labeling index were obtained at the two time points. Thus, different parameter sets were obtained for each treatment.
The numerical solutions of Equation (1) together with (3), (4) or (5) were computed using the Matlab function ode45.

2.5. RT-Only Virtual Experiments

We performed simulations for the patients who received only RT, which consisted of shifting the RT starting time to later times to test if the model results were in agreement with the observations in the clinical trial by Van den Bent et al. [44], who found that deferring RT in LGGs did not have an effect on survival (see Section 3.1.2).
We also simulated the tumor behavior when the time interval between RT sessions was increased. To take into account the 18% error in real tumor volumes, we choose randomized tumor volumes within the range defined by error and we computed the value of the parameters best describing the tumor evolution (see Section 3.2).

2.6. In-Silico Patients

The data available in Table 3 provide either parameters for RT or TMZ, but only one patient received both treatments concurrently. Then, in order to construct virtual patients through whom to compare the efficacy of combination treatments, we computed the most representative ranges for each parameter in Table 3, that is, P 0 + Q 0 [ 2.52 , 146 ] cm 3 , Ki-67 0 { 1 , 1.5 , , 6 } % , ρ [ 0.01 , 0.15 ] 1/day, β 1 [ 0.23 , 10.60 ] 1/day, β 3 [ 1 × 10 5 , 0.021 ] 1/day, S q [ 0 , 0.94 ] , S p [ 0 , 0.7 ] , α 1 [ 2.3 , 45 ] cm 3 / ( μ g day) and α 2 [ 0.37 , 20 ] cm 3 / ( μ g day).
The parameters for adding a virtual therapy or creating an in-silico patient were then chosen by randomly sampling the above corresponding intervals with uniform distribution, and β 2 was calculated using
β 2 = Ki - 67 0 % β 1 100 Ki 67 0 % ρ 100 .
This equation was obtained by taking the limit P ( t ) /( P ( t ) + Q ( t ) ) as t goes to infinity in the analytic solution of Equation (1) (see SI Section S1). In this way, we ensure values of Ki-67 in the range characteristic of low-grade gliomas.
For radiation therapy treatment, the number of fractions was taken to be 30, which is in line with usual clinical practice. For chemotherapy treatments, we assumed the number of TMZ cycles to be 12, that is, the average number of cycles received by the real patients included in this study (see Table 3).

2.7. In-Silico Trial Comparing Different Treatment Modalities/Schedules

To find the best therapeutic approach for large patient cohorts, we designed a threearm in-silico trial. A total of 3000 virtual patients were created by choosing parameters as described in Section 2.6.
Patients in each arm were randomized for the following treatments:
I.
Control group without treatment;
II.
RT (1.8 Gy)+TMZ (concurrently 75 mg/m 2 /day and post-radiation 150–200 mg/m 2 / day) according to RTOG 0424. For simplicity, we assumed a fixed-dose equal to 150 mg/m 2 /day in all TMZ cycles;
III.
RT+TMZ matching the doses of RT (1.8 Gy) to those of TMZ (150 mg/m 2 /day) in the first five days of each cycle (of length 28 days).
Using the MatSurv MATLAB package [45], Kaplan–Meier survival curves for overall survival were then calculated. A death event was assumed to happen when tumor volume reached a fatal size (V F ) randomly distributed between 113.1 and 268.08 cm 3 , which is satisfied when V F ( P 0 + Q 0 ) > 40 cm 3 , to avoid P 0 + Q 0 V F . In line with other work, Pérez-García et al. [26], where a fatal volume for LGGs was assumed to be 268.08 cm 3 .
As low-grade gliomas are less aggressive and their size is measured in T2/FLAIR MRI sequences, we assumed a slightly larger volume and added the variability as a surrogate of tumor location, which substantially influences the critical tumor size. Patients who only reached the fatal volume after 15 years were considered censored events.

3. Results

3.1. Validation of the Mathematical Model

3.1.1. The Mathematical Model Given by Equation (1) Described the Tumor Longitudinal Dynamics in Patients

The mathematical model was fitted for each patient as described in Section 2.4, and best parameters obtained as listed in Table 3. The model was flexible enough to describe the dynamics in all patients faithfully.
Figure 2 shows some examples and, in Figure 3, the dynamics of each cell population for patient 36 is displayed. The agreement between the observed volumes and the simulated values was excellent. It can be seen that the tumor evolution was different for each patient, no matter whether RT, TMZ or both were administered. In addition, the dynamics of the virtual Ki-67 reflects the low mitotic activity of LGGs reported in clinical evidence [2]. Best fits for all patients are shown in the SI Section S2.

3.1.2. Delaying Radiation Therapy Did Not Have an Effect on Overall Survival

The treatment starting time was moved to different times after the real time for patients who received only RT, with the exception of patient 107, in whose case it was moved to a time earlier than the real one, because at the beginning of RT the tumor volume had almost reached the fatal limit, which was taken to be 268.08 cm 3 , corresponding to the largest volume defined in Section 2.7; that is, a sphere of 8 cm in diameter. This choice allowed enough simulation time to compare the evolution of the tumor growth curves. Figure 4 shows typical results. It can be seen that the model fully agrees with the results of Van den Bent et al. [44].

3.2. Protracted Radiotherapy Schemes Lead to Increased Tumor Growth

Some mathematical modeling work on LGGs has concluded that it is possible to space the sessions of RT without losing the therapeutic effect [18,19]. However these studies did not account for the repopulation of the compartment of proliferative cells due to the stimulation of quiescent cells after RT. We used our model to investigate the implications of considering this biological mechanism. Figure 5 displays the tumor longitudinal dynamics when the time interval between RT fractions was chosen to be either 5 days ( Θ = 5 days), 2 days ( Θ = 2 days) or the standard interval ( Θ = 1 day with a break at weekends). In addition, we studied if these results are affected by the 18% error of the tumor volumes. Figure 6 shows that qualitatively the main results are not affected by the error. It can be seen that, at least for these specific dose interval choices, the larger the spacing between doses, the larger the final tumor volume was found to be.

3.3. Concurrent Cycles of RT and TMZ Produce the Longest-Lasting Therapeutic Effect

Next, we used our mathematical model as a test bed to carry out in-silico experiments, adding the missing therapy as described in Section 2.6. The parameter values of the corresponding virtual therapy are listed in Table 4. Our goal was to find optimal treatment regimens.
Based on the structure of the model interactions, we hypothesized that it would be possible to space the RT sessions as long as they are administered simultaneously with TMZ, without losing therapeutic effect. The rational for this hypothesis is that radiotherapy stimulates quiescent cells (which are the largest subpopulation in LGGs) to proliferate, and once in their state of proliferation, TMZ is capable of damaging them. To test this, we computed the overall survival gain (OSG) using two different schemes. The first was concomitant radio-chemotherapy according to RTOG 0424, where RT and TMZ were started at the same time and were given concurrently in weeks 1–6. After a four-week break, TMZ cycles were given without any overlap with RT. The second scheme delivered only cycles of chemotherapy, but radiotherapy (until the course was completed) was given only on the days where chemotherapy was administered (see Figure 7). Matching both therapies led to a longer OS (Figure 8). Similar results were obtained for patient 36 and all of the virtual patients (see Table 4).

3.4. In-Silico Simulation of Clinical Trials

To test the previous results, in a large enough virtual patient population, we designed an in-silico study following the methodology described in Section 2.7. A total of 3000 virtual patients were grouped into three arms. For the groups that received treatment, the numbers of TMZ cycles and RT sessions were as described in Section 2.6. Patients surviving more than 15 years in silico were considered to be censored events. Figure 9 shows the Kaplan-Meier curves for overall survival. It is clear that the simultaneous RT+TMZ protocol had the best results in terms of survival.
According to the log-rank test, the differences between all the survival curves were statistically significant (p-value < 10 4 ).
One of the most important questions in clinical trials is to estimate the number of patients required to achieve significance.
To obtain that number computationally, we performed an in-silico study comparing the potential advantages of the RT+TMZ simultaneous schedule over the RTOG 0424 protocol. 80 virtual patients were equally distributed in a two-arms (RT+TMZ simultaneous and RTOG 0424 protocols) in-silico trial. The virtual follow-up period was 15 years. The p-value for the difference between the survival curves of each arm was calculated. We repeated the same experiment for another different 80 virtual patients and we followed the same process until completing 20 trials. Then the number of patients was increased by 80 (40 per arm) until a number of 800 was reached and the same methodology explained above was applied. The results are shown in Figure 10. It can be seen in Figure 10 that, as the number of patients increased, more trials were statistically significant, and in fact, for N 360 , only one trial out of 20 resulted in p > 0.05 . For a population of 720 patients (360 per arm) or more, the probability of finding a positive result in the trial would be very high (≥95%).

4. Discussion and Conclusions

In this paper, we have studied in silico what would be the most effective way, in terms of tumor control, of delivering a combination treatment of radiation therapy and temozolomide for LGGs. Radiation therapy was assumed to consist of 30 fractions of a standard 1.8 Gy dose given daily excluding weekends. Chemotherapy cycles consisted of five consecutive daily doses followed by a rest period. To find the best therapeutic combinations we used a biomathematical model. Several key biological aspects were incorporated into the mathematical model, such as the presence of a majority of quiescent cells and a small fraction of proliferative cells that were the target of cytotoxic therapies. Two different types of death process for TMZ were included, the first accounting for ‘early’ cell death and the second leading to ‘late’ cell death, as happens with RT. Finally, and as a distinctive feature of our modeling approach, we accounted for the accelerated repopulation stimulated by radiation.
It is noteworthy that the mathematical model, despite its simplicity, successfully described the tumor longitudinal volumetric dynamics of the LGG patients in our dataset treated with TMZ and/or RT. The model incorporated different parameters: two related to the initial tumor volume and composition P 0 , Q 0 , four related to biological rates ρ , β 1 , β 2 , β 3 , two for each type of treatment, either radiotherapy S q , S p or chemotherapy α 1 , α 2 and one parameter related to the action of chemotherapy λ . The dynamics of the virtual Ki-67 reflected the very low mitotic activity of LGGs reported in clinical practice [2].
Many mathematical models of LGGs have been constructed to understand the response of this disease to RT [16,17,18,19,46] and chemotherapy [16,22,23,25,26]. Several previous studies [17,19] have proposed that enlarging the interval between RT fractions could be a strategy leading to survival advantages. However, with these models it is not possible to simulate the well-quantified mitotic activity of LGGs (Ki-67 labeling index around 4–5%) and also they do not consider the repopulation of proliferating cancer cells due to the stimulatory effect of RT. As seen in Figure 5, our model predicts worse tumor control when RT sessions are more spaced out. The reason is that the time interval between doses allows for quiescent stimulated cells to grow. Thus, according to our more realistic model, the non-standard schemes studied in [19] would not lead to survival advantages in their clinical application.
To our knowledge, this is the first time that the stimulatory effect produced by radiotherapy on quiescent cells has been taken into account in a mathematical model of LGGs. This effect makes quiescent cells proliferate making the chemotherapy treatment delivered concomitantly more effective. This is a key biological effect of radiotherapy that has to be taken into account in mathematical models to avoid obtaining unrealistic results, at least in the range of doses per fraction typically used for LGGs.
We found that our model is in line with the clinical trial by Van den Bent et al. [44], in the sense that deferring RT in LGGs does not alter survival time. A similar scenario was found with the administration of TMZ, where advancing or delaying treatment did not affect overall survival [26].
The main clinical contribution of our virtual simulations was that the simultaneous RT+TMZ protocol produced an overall survival gain. These in-silico gains are substantial when compared with the long-term results of the Therapy Oncology Group 0424 clinical trial carried out on High-Risk LGG patients [5]. Our research opens the door to new RT+TMZ regimens where RT fractions match TMZ treatment, allowing patients to have longer recovery times between cycles of RT/TMZ, and therefore fewer side effects than the RTOG 0424 protocol.
Our computational study predicts that a clinical trial with 700 patients would provide most likely a positive result showing the effectiveness of this therapeutic schedule. These numbers may seem large to be implemented in a real study. However, we included both virtual astrocytoma and oligodendroglioma patients as in our real cohort. Since astrocytoma patients have a shorter survival, a proof of concept study could focus only on those patients were much lower numbers would suffice.
Although the model faithfully reproduced basic aspects of radiotherapy and TMZ treatment, it has several limitations. The first is that we only had longitudinal volumetric data, thus the model was built to describe tumor evolution as a whole. Such models cannot provide information on how tumor cell density is spatially distributed throughout the brain, and do not account for details of the well-known infiltrative nature of gliomas. An extension of this work could be based on reaction-diffusion models accounting for spatial effects. The second is that, even though there are many studies on the repopulation of cancer cells during RT, this mechanism has not yet been proved in LGGs. In our modeling approach, we did not consider the development of resistance to TMZ that is known to arise after long treatments with the drug [47]. This is why we focused our study on treatments with limited number of cycles and in optimizing the outcome from combinations with RT. To develop a mathematical study of long maintenance treatments with TMZ, one should incorporate additional compartments of resistant and probably persister cells and describe the effect of the treatment in inducing transitions between them.
Since the current dose of radiotherapy is not evidence-based, it would be of the highest scientific interest to use a mathematical approach as the one presented here to find the best dose and fractionation scheme to treat LGGs. However, scarcity of longitudinal data of tumors treated with different radiation therapy schedules would limit the relevance of the results. In Ref. [19], the authors studied the problem theoretically using a linear-quadratic model of response to RT and proposed protracted schemes as the ones best controlling LGG growth. However, without the support of real-world data, it is unclear what would be the performance of those schemes when adding RT stimulatory effects and TMZ.
The main objective of this study was to provide a theoretical framework for the synergy of the simultaneous administration of both TMZ and RT in LGGs, which was not affected by the 18% error of real tumor volume as shown in the simulations. Using this model and the methodology developed here a wider range of possibilities, for example, RT+TMZ different to the protocols here studied, only RT for various doses and schedules and only TMZ in diverse regimens could be analyzed in further works. However, additional data supporting the parameter choice would be necessary.
Since it is very difficult to test the effectiveness of long-cycle strategies either in vitro or in animal models, it is a good example of how computational mathematical models can be useful in providing some preclinical evidence.
We hope this computational study will provide a theoretical grounding for the development of clinical and experimental studies and for the definition of standardized optimal treatment protocols combining TMZ and RT for low-grade glioma patients, with improved survival and better quality of life.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jpm11101036/s1.

Author Contributions

Conceptualization, V.M.P.-G.; methodology, V.M.P.-G., A.G.; software, L.E.A.-H., A.G., J.B.-B., V.M.P.-G.; validation, L.E.A.-H., A.G., J.B.-B., V.M.P.-G.; visualization, V.M.P.-G., A.G.; formal analysis, L.E.A.-H., A.G., J.B.-B., V.M.P.-G.; investigation, L.E.A.-H., A.G., P.S., M.M., L.P.-R., J.B.-B., V.M.P.-G.; resources, V.M.P.-G., M.M.; data curation, V.M.P.-G., P.S., M.M.; writing—original draft preparation, L.E.A.-H., A.G., J.B.-B., V.M.P.-G.; writing—review and editing, L.E.A.-H., A.G., L.P.-R., J.B.-B., V.M.P.-G.; supervision, A.G., J.B.-B., V.M.P.-G.; project administration, V.M.P.-G.; funding acquisition, V.M.P.-G., A.G., P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been supported by the James S. Mc. Donnell Foundation (USA) 21st Century Science Initiative in Mathematical and Complex Systems Approaches for Brain Cancer (Collaborative award 220020450, doi:10.37717/220020560). This research has also been supported by grants awarded to V.M.P-G by the Ministerio de Economía y Competitividad/FEDER, Spain (grant number MTM2015-71200-R), Junta de Comunidades de Castilla-La Mancha, Spain (grant number SBPLY/17/180501/000154). This research has also been supported by a grant awarded to J.B-B by the Junta de Comunidades de Castilla-La Mancha, Spain (grant number SBPLY/19/180501/000211).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and was approved by Kantonale Ethikkommission Bern (Bern, Switzerland), with approval number: 07.09.72.

Informed Consent Statement

Patient consent was waived as the data were obtained from retrospective observational clinical studies that were approved by the corresponding institutional review boards.

Data Availability Statement

The study was approved by Kantonale Ethikkommission Bern (Bern, Switzerland), with approval number: 07.09.72.

Acknowledgments

L.E.A-H wishes to thank CONACYT for the financial support granted through scholarship 779584.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
LGGLow-Grade Glioma
OOligodendroglioma
AAstrocytoma
OAOligoastrocytoma
RTRadiation Therapy
TMZTemozolomide
CSCycle Size
PCVProcarbazine, Lomustine and Vincristine
MRIMagnetic Resonance Imaging
OSGOverall Survival Gain
NANot Available

References

  1. Lombardi, G.; Barresi, V.; Castellano, A.; Tabouret, E.; Pasqualetti, F.; Salvalaggio, A.; Cerretti, G.; Caccese, M.; Padovan, M.; Zagonel, V.; et al. Clinical management of diffuse low-grade gliomas. Cancers 2020, 12, 3008. [Google Scholar] [CrossRef]
  2. Soffietti, R.; Baumert, B.; Bello, L.; Von Deimling, A.; Duffau, H.; Frénay, M.; Grisold, W.; Grant, R.; Graus, F.; Hoang-Xuan, K.; et al. Guidelines on management of low-grade gliomas: Report of an EFNS–EANO* Task Force. Eur. J. Neurol. 2010, 17, 1124–1133. [Google Scholar] [CrossRef]
  3. Smits, A.; Jakola, A.S. Clinical presentation, natural history, and prognosis of diffuse low-grade gliomas. Neurosurg. Clin. 2019, 30, 35–42. [Google Scholar] [CrossRef]
  4. Buckner, J.C.; Gesme, D., Jr.; O’Fallon, J.R.; Hammack, J.E.; Stafford, S.; Brown, P.D.; Hawkins, R.; Scheithauer, B.W.; Erickson, B.J.; Levitt, R.; et al. Phase II trial of procarbazine, lomustine, and vincristine as initial therapy for patients with low-grade oligodendroglioma or oligoastrocytoma: Efficacy and associations with chromosomal abnormalities. J. Clin. Oncol. 2003, 21, 251–255. [Google Scholar] [CrossRef] [PubMed]
  5. Fisher, B.J.; Pugh, S.L.; Macdonald, D.R.; Chakravatri, A.; Lesser, G.J.; Fox, S.; Rogers, C.L.; Werner-Wasik, M.; Doyle, T.; Bahary, J.P.; et al. Phase 2 Study of a Temozolomide-Based Chemoradiation Therapy Regimen for High-Risk, Low-Grade Gliomas: Long-Term Results of Radiation Therapy Oncology Group 0424. Int. J. Radiat. Oncol. Biol. Phys. 2020, 107, 720–725. [Google Scholar] [CrossRef] [PubMed]
  6. Baumert, B.G.; Hegi, M.E.; van den Bent, M.J.; von Deimling, A.; Gorlia, T.; Hoang-Xuan, K.; Brandes, A.A.; Kantor, G.; Taphoorn, M.J.; Hassel, M.B.; et al. Temozolomide chemotherapy versus radiotherapy in high-risk low-grade glioma (EORTC 22033-26033): A randomised, open-label, phase 3 intergroup study. Lancet Oncol. 2016, 17, 1521–1532. [Google Scholar] [CrossRef]
  7. Weller, M.; van den Bent, M.; Preusser, M.; Le Rhun, E.; Tonn, J.C.; Minniti, G.; Bendszus, M.; Balana, C.; Chinot, O.; Dirven, L.; et al. EANO guidelines on the diagnosis and treatment of diffuse gliomas of adulthood. Nat. Rev. Clin. Oncol. 2021, 18, 170–186. [Google Scholar] [CrossRef] [PubMed]
  8. Jiang, T.; Nam, D.H.; Ram, Z.; Poon, W.S.; Wang, J.; Boldbaatar, D.; Mao, Y.; Ma, W.; Mao, Q.; You, Y.; et al. Clinical practice guidelines for the management of adult diffuse gliomas. Cancer Lett. 2021, 499, 60–72. [Google Scholar] [CrossRef] [PubMed]
  9. Liu, R.; Solheim, K.; Polley, M.Y.; Lamborn, K.R.; Page, M.; Fedoroff, A.; Rabbitt, J.; Butowski, N.; Prados, M.; Chang, S.M. Quality of life in low-grade glioma patients receiving temozolomide. Neuro-Oncol. 2009, 11, 59–68. [Google Scholar] [CrossRef]
  10. Buckner, J.C.; Shaw, E.G.; Pugh, S.L.; Chakravarti, A.; Gilbert, M.R.; Barger, G.R.; Coons, S.; Ricci, P.; Bullard, D.; Brown, P.D.; et al. Radiation plus procarbazine, CCNU, and vincristine in low-grade glioma. N. Engl. J. Med. 2016, 374, 1344–1355. [Google Scholar] [CrossRef]
  11. Stupp, R.; Mason, W.P.; Van Den Bent, M.J.; Weller, M.; Fisher, B.; Taphoorn, M.J.; Belanger, K.; Brandes, A.A.; Marosi, C.; Bogdahn, U.; et al. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N. Engl. J. Med. 2005, 352, 987–996. [Google Scholar] [CrossRef] [PubMed]
  12. Fisher, B.J.; Hu, C.; Macdonald, D.R.; Lesser, G.J.; Coons, S.W.; Brachman, D.G.; Ryu, S.; Werner-Wasik, M.; Bahary, J.P.; Liu, J.; et al. Phase 2 study of temozolomide-based chemoradiation therapy for high-risk low-grade gliomas: Preliminary results of Radiation Therapy Oncology Group 0424. Int. J. Radiat. Oncol. Biol. Phys. 2015, 91, 497–504. [Google Scholar] [CrossRef]
  13. Wang, J.; Wang, Y.; He, Y.; Guan, H.; He, L.; Mu, X.; Peng, X. Radiotherapy versus radiotherapy combined with temozolomide in high-risk low-grade gliomas after surgery: Study protocol for a randomized controlled clinical trial. Trials 2019, 20, 1–7. [Google Scholar] [CrossRef]
  14. Barbolosi, D.; Ciccolini, J.; Lacarelle, B.; Barlési, F.; André, N. Computational oncology—Mathematical modelling of drug regimens for precision medicine. Nat. Rev. Clin. Oncol. 2016, 13, 242–254. [Google Scholar] [CrossRef]
  15. Malinzi, J.; Basita, K.B.; Padidar, S.; Adeola, H.A. Prospect for application of mathematical models in combination cancer treatments. Inform. Med. Unlocked 2021, 23, 100534. [Google Scholar] [CrossRef]
  16. Ribba, B.; Kaloshi, G.; Peyre, M.; Ricard, D.; Calvez, V.; Tod, M.; Čajavec-Bernard, B.; Idbaih, A.; Psimaras, D.; Dainese, L.; et al. A tumor growth inhibition model for low-grade glioma treated with chemotherapy or radiotherapy. Clin. Cancer Res. 2012, 18, 5071–5080. [Google Scholar] [CrossRef] [PubMed]
  17. Pérez-García, V.M.; Bogdanska, M.; Martínez-González, A.; Belmonte-Beitia, J.; Schucht, P.; Pérez-Romasanta, L.A. Delay effects in the response of low-grade gliomas to radiotherapy: A mathematical model and its therapeutical implications. Math. Med. Biol. J. IMA 2015, 32, 307–329. [Google Scholar] [CrossRef] [PubMed]
  18. Pérez-García, V.M.; Pérez-Romasanta, L.A. Extreme protraction for low-grade gliomas: Theoretical proof of concept of a novel therapeutical strategy. Math. Med. Biol. J. IMA 2016, 33, 253–271. [Google Scholar] [CrossRef]
  19. Henares-Molina, A.; Benzekry, S.; Lara, P.C.; García-Rojo, M.; Pérez-García, V.M.; Martínez-González, A. Non-standard radiotherapy fractionations delay the time to malignant transformation of low-grade gliomas. PLoS ONE 2017, 12, e0178552. [Google Scholar] [CrossRef] [PubMed]
  20. Rodríguez, C.R.; Calvo, G.F.; Ramis-Conde, I.; Belmonte-Beitia, J. Stochastic modelling of slow-progressing tumors: Analysis and applications to the cell interplay and control of low grade gliomas. Commun. Nonlinear Sci. Numer. Simul. 2017, 49, 63–80. [Google Scholar] [CrossRef]
  21. Budia, I.; Alvarez-Arenas, A.; Woolley, T.; Calvo, G.; Belmonte-Beitia, J. Radiation protraction schedules for low-grade gliomas: A comparison between different mathematical models. J. R. Soc. Interface 2019, 16, 20190665. [Google Scholar] [CrossRef] [PubMed]
  22. Mazzocco, P.; Barthélémy, C.; Kaloshi, G.; Lavielle, M.; Ricard, D.; Idbaih, A.; Psimaras, D.; Renard, M.A.; Alentorn, A.; Honnorat, J.; et al. Prediction of response to temozolomide in low-grade glioma patients based on tumor size dynamics and genetic characteristics. CPT Pharmacomet. Syst. Pharmacol. 2015, 4, 728–737. [Google Scholar] [CrossRef] [PubMed]
  23. Mazzocco, P.; Honnorat, J.; Ducray, F.; Ribba, B. Increasing the time interval between PCV chemotherapy cycles as a strategy to improve duration of response in low-grade gliomas: Results from a model-based clinical trial simulation. Comput. Math. Methods Med. 2015, 2015, 297903. [Google Scholar] [CrossRef] [PubMed]
  24. Rojas, C.; Belmonte-Beitia, J.; Pérez-García, V.M.; Maurer, H. Dynamics and optimal control of chemotherapy for low grade gliomas: Insights from a mathematical model. Discret. Contin. Dyn. Syst.-B 2016, 21, 1895–1915. [Google Scholar] [CrossRef]
  25. Bogdańska, M.U.; Bodnar, M.; Belmonte-Beitia, J.; Murek, M.; Schucht, P.; Beck, J.; Pérez-García, V.M. A mathematical model of low grade gliomas treated with temozolomide and its therapeutical implications. Math. Biosci. 2017, 288, 1–13. [Google Scholar] [CrossRef]
  26. Pérez-García, V.M.; Ayala-Hernández, L.E.; Belmonte-Beitia, J.; Schucht, P.; Murek, M.; Raabe, A.; Sepúlveda, J. Computational design of improved standardized chemotherapy protocols for grade II oligodendrogliomas. PLoS Comput. Biol. 2019, 15, e1006778. [Google Scholar] [CrossRef] [PubMed]
  27. Bodnar, M.; Pérez, M.V. Mathematical and numerical analysis of low-grade gliomas model and the effects of chemotherapy. Commun. Nonlinear Sci. Numer. Simul. 2019, 72, 552–564. [Google Scholar] [CrossRef]
  28. Bashkirtseva, I.; Ryashko, L.; López, Á.G.; Seoane, J.M.; Sanjuán, M.A. The effect of time ordering and concurrency in a mathematical model of chemoradiotherapy. Commun. Nonlinear Sci. Numer. Simul. 2021, 96, 105693. [Google Scholar] [CrossRef]
  29. Pallud, J.; Mandonnet, E.; Duffau, H.; Kujas, M.; Guillevin, R.; Galanaud, D.; Taillandier, L.; Capelle, L. Prognostic value of initial magnetic resonance imaging growth rates for World Health Organization grade II gliomas. Ann. Neurol. 2006, 60, 380–383. [Google Scholar] [CrossRef]
  30. Mandonnet, E.; Delattre, J.Y.; Tanguy, M.L.; Swanson, K.R.; Carpentier, A.F.; Duffau, H.; Cornu, P.; Van Effenterre, R.; Alvord, E.C., Jr.; Capelle, L. Continuous growth of mean tumor diameter in a subset of grade II gliomas. Ann. Neurol. 2003, 53, 524–528. [Google Scholar] [CrossRef]
  31. Pérez-Beteta, J.; Molina-García, D.; Ortiz-Alhambra, J.A.; Fernández-Romero, A.; Luque, B.; Arregui, E.; Calvo, M.; Borrás, J.M.; Meléndez, B.; Rodríguez de Lope, Á.; et al. Tumor surface regularity at MR imaging predicts survival and response to surgery in patients with glioblastoma. Radiology 2018, 288, 218–225. [Google Scholar] [CrossRef]
  32. Yeh, A.C.; Ramaswamy, S. Mechanisms of cancer cell dormancy—Another hallmark of cancer? Cancer Res. 2015, 75, 5014–5022. [Google Scholar] [CrossRef]
  33. Mills, C.C.; Kolb, E.; Sampson, V.B. Development of chemotherapy with cell-cycle inhibitors for adult and pediatric cancer therapy. Cancer Res. 2018, 78, 320–325. [Google Scholar] [CrossRef]
  34. Pawlik, T.M.; Keyomarsi, K. Role of cell cycle in mediating sensitivity to radiotherapy. Int. J. Radiat. Oncol. Biol. Phys. 2004, 59, 928–942. [Google Scholar] [CrossRef]
  35. Adenis, L.; Plaszczynski, S.; Grammaticos, B.; Pallud, J.; Badoual, M. Radiotherapy Effects on Diffuse Low-Grade Gliomas: Confronting Theory With Clinical Data. arXiv 2021, arXiv:2107.02534. [Google Scholar]
  36. Fogarty, C.E.; Bergmann, A. Killers creating new life: Caspases drive apoptosis-induced proliferation in tissue repair and disease. Cell Death Differ. 2017, 24, 1390–1400. [Google Scholar] [CrossRef] [PubMed]
  37. Zimmerman, M.A.; Huang, Q.; Li, F.; Liu, X.; Li, C.Y. Cell death–stimulated cell proliferation: A tissue regeneration mechanism usurped by tumors during radiotherapy. Semin. Radiat. Oncol. 2013, 23, 288–295. [Google Scholar] [CrossRef] [PubMed]
  38. He, S.; Cheng, J.; Sun, L.; Wang, Y.; Wang, C.; Liu, X.; Zhang, Z.; Zhao, M.; Luo, Y.; Tian, L.; et al. HMGB1 released by irradiated tumor cells promotes living tumor cell proliferation via paracrine effect. Cell Death Dis. 2018, 9, 648. [Google Scholar] [CrossRef]
  39. Davis, A.J.; Tannock, I.F. Repopulation of tumour cells between cycles of chemotherapy: A neglected factor. Lancet Oncol. 2000, 1, 86–93. [Google Scholar] [CrossRef]
  40. Marcu, L. Tumour repopulation and the role of abortive division in squamous cell carcinomas during chemotherapy. Cell Prolif. 2014, 47, 318–325. [Google Scholar] [CrossRef] [PubMed]
  41. Mosteller, R. Simplified calculation of body-surface area. N. Engl. J. Med. 1987, 317, 1098. [Google Scholar]
  42. Portnow, J.; Badie, B.; Chen, M.; Liu, A.; Blanchard, S.; Synold, T.W. The neuropharmacokinetics of temozolomide in patients with resectable brain tumors: Potential implications for the current approach to chemoradiation. Clin. Cancer Res. 2009, 15, 7092–7098. [Google Scholar] [CrossRef]
  43. Hammond, L.A.; Eckardt, J.R.; Baker, S.D.; Eckhardt, S.G.; Dugan, M.; Forral, K.; Reidenberg, P.; Statkevich, P.; Weiss, G.R.; Rinaldi, D.A.; et al. Phase I and pharmacokinetic study of temozolomide on a daily-for-5-days schedule in patients with advanced solid malignancies. J. Clin. Oncol. 1999, 17, 2604–2613. [Google Scholar] [CrossRef]
  44. Van den Bent, M.; Afra, D.; De Witte, O.; Hassel, M.B.; Schraub, S.; Hoang-Xuan, K.; Malmström, P.; Collette, L.; Piérart, M.; Mirimanoff, R.; et al. Long-term efficacy of early versus delayed radiotherapy for low-grade astrocytoma and oligodendroglioma in adults: The EORTC 22845 randomised trial. Lancet 2005, 366, 985–990. [Google Scholar] [CrossRef]
  45. Creed, J.H.; Gerke, T.A.; Berglund, A.E. MatSurv: Survival analysis and visualization in MATLAB. J. Open Source Softw. 2020, 5, 1830. [Google Scholar] [CrossRef]
  46. Badoual, M.; Gerin, C.; Deroulers, C.; Grammaticos, B.; Llitjos, J.F.; Oppenheim, C.; Varlet, P.; Pallud, J. Oedema-based model for diffuse low-grade gliomas: Application to clinical cases under radiotherapy. Cell Prolif. 2014, 47, 369–380. [Google Scholar] [CrossRef] [PubMed]
  47. Ricard, D.; Kaloshi, G.; Amiel-Benouaich, A.; Lejeune, J.; Marie, Y.; Mandonnet, E.; Kujas, M.; Mokhtari, K.; Taillibert, S.; Laigle-Donadey, F.; et al. Dynamic history of low-grade gliomas before and after temozolomide treatment. Ann. Neurol. Off. J. Am. Neurol. Assoc. Child Neurol. Soc. 2007, 61, 484–490. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Diagram illustrating the model interactions. Proliferative cells P ( t ) grow at rate ρ , some enter a quiescent state Q ( t ) at rate β 1 and vice versa at rate β 2 . Radiation damages a fraction 1 S p ( d ) of P ( t ) cells to produce a subpopulation of lethally damaged cells D ( t ) that die at rate β 3 , and also stimulates a fraction 1 S q ( d ) of Q ( t ) to proliferate. Chemotherapy C ( t ) causes some P ( t ) cells to undergo delayed death (as radiation does) at rate α 1 and immediate death at rate α 2 .
Figure 1. Diagram illustrating the model interactions. Proliferative cells P ( t ) grow at rate ρ , some enter a quiescent state Q ( t ) at rate β 1 and vice versa at rate β 2 . Radiation damages a fraction 1 S p ( d ) of P ( t ) cells to produce a subpopulation of lethally damaged cells D ( t ) that die at rate β 3 , and also stimulates a fraction 1 S q ( d ) of Q ( t ) to proliferate. Chemotherapy C ( t ) causes some P ( t ) cells to undergo delayed death (as radiation does) at rate α 1 and immediate death at rate α 2 .
Jpm 11 01036 g001
Figure 2. (a,b,e,f) Examples of longitudinal volumetric dynamics (green circles) with their respective 18 % error bars (green vertical lines) and best fits obtained using Equation (1) (solid purple line). Data shown correspond to patients 10 , 36 , 40 and 66 (green circles). (c,d,g,h) Time evolution of 100 P ( t ) / ( P ( t ) + Q ( t ) ) interpreted as virtual (solid purple line) and initial (red circle) Ki-67 labeling index. The vertical dashed lines indicate the onset and end times of the treatments received by each patient (either RT, QT or both).
Figure 2. (a,b,e,f) Examples of longitudinal volumetric dynamics (green circles) with their respective 18 % error bars (green vertical lines) and best fits obtained using Equation (1) (solid purple line). Data shown correspond to patients 10 , 36 , 40 and 66 (green circles). (c,d,g,h) Time evolution of 100 P ( t ) / ( P ( t ) + Q ( t ) ) interpreted as virtual (solid purple line) and initial (red circle) Ki-67 labeling index. The vertical dashed lines indicate the onset and end times of the treatments received by each patient (either RT, QT or both).
Jpm 11 01036 g002aJpm 11 01036 g002b
Figure 3. (ac) Evolution of P ( t ) , Q ( t ) and D ( t ) for patient 36. The sum of each of these cell populations corresponds to the total tumor volume.
Figure 3. (ac) Evolution of P ( t ) , Q ( t ) and D ( t ) for patient 36. The sum of each of these cell populations corresponds to the total tumor volume.
Jpm 11 01036 g003
Figure 4. (ah) Delaying radiotherapy does not have an effect on overall survival. RT was simulated to start one year later (gray solid line) than when actually received by the patient (solid purple line), except for patient 107, in whose case it was simulated to occur one year earlier. The horizontal black dashed line represents the fatal volume and the green circles, with their respective error bars, the longitudinal volumetric data.
Figure 4. (ah) Delaying radiotherapy does not have an effect on overall survival. RT was simulated to start one year later (gray solid line) than when actually received by the patient (solid purple line), except for patient 107, in whose case it was simulated to occur one year earlier. The horizontal black dashed line represents the fatal volume and the green circles, with their respective error bars, the longitudinal volumetric data.
Jpm 11 01036 g004
Figure 5. (ah) Protracted radiotherapy schemes lead to increased tumor growth in silico. Simulations of longitudinal tumor growth dynamics of patients receiving radiation therapy fractionations spaced at 5 days (black solid line), 2 days (red solid line) and the conventional scheme (purple solid line). The vertical red dashed line shows the treatment start time and the green circles the longitudinal volumetric data.
Figure 5. (ah) Protracted radiotherapy schemes lead to increased tumor growth in silico. Simulations of longitudinal tumor growth dynamics of patients receiving radiation therapy fractionations spaced at 5 days (black solid line), 2 days (red solid line) and the conventional scheme (purple solid line). The vertical red dashed line shows the treatment start time and the green circles the longitudinal volumetric data.
Jpm 11 01036 g005aJpm 11 01036 g005b
Figure 6. Same simulations as in Figure 5 for three different sets of parameters taking into account the 18% error in the real tumor volumes. (a) V 0 = 43.56 , ρ = 0.091 , β 1 = 1.08 , β 2 = 0.03 , β 3 = 0.0076 , S q = 0.861 and S p = 0.019 . (b) V 0 = 44.86 , ρ = 0.094 , β 1 = 1.04 , β 2 = 0.029 , β 3 = 0.0077 , S q = 0.866 and S p = 0.003 and (c) V 0 = 44.99 , ρ = 0.101 , β 1 = 1 , β 2 = 0.027 , β 3 = 0.0081 , S q = 0.869 and S p = 0 .
Figure 6. Same simulations as in Figure 5 for three different sets of parameters taking into account the 18% error in the real tumor volumes. (a) V 0 = 43.56 , ρ = 0.091 , β 1 = 1.08 , β 2 = 0.03 , β 3 = 0.0076 , S q = 0.861 and S p = 0.019 . (b) V 0 = 44.86 , ρ = 0.094 , β 1 = 1.04 , β 2 = 0.029 , β 3 = 0.0077 , S q = 0.866 and S p = 0.003 and (c) V 0 = 44.99 , ρ = 0.101 , β 1 = 1 , β 2 = 0.027 , β 3 = 0.0081 , S q = 0.869 and S p = 0 .
Jpm 11 01036 g006
Figure 7. Schematic representation of the combined treatment protocols. In the RTOG 0424 protocol, RT was given concurrently with TMZ (75 mg/m 2 /day) in weeks 1–6 (top of treatment timeline) while in the simultaneous scheme, RT was only administered together with TMZ (150 mg/m 2 /day) (bottom of treatment timeline).
Figure 7. Schematic representation of the combined treatment protocols. In the RTOG 0424 protocol, RT was given concurrently with TMZ (75 mg/m 2 /day) in weeks 1–6 (top of treatment timeline) while in the simultaneous scheme, RT was only administered together with TMZ (150 mg/m 2 /day) (bottom of treatment timeline).
Jpm 11 01036 g007
Figure 8. (ad) Simulated tumor growth curves for patient 36 and for three in-silico patients treated according to RTOG 0424 (purple solid line) and the RT+TMZ simultaneous (red solid line) regimens. Treatment onset is denoted by a vertical green dashed line and the fatal tumor volume by a horizontal black dashed line.
Figure 8. (ad) Simulated tumor growth curves for patient 36 and for three in-silico patients treated according to RTOG 0424 (purple solid line) and the RT+TMZ simultaneous (red solid line) regimens. Treatment onset is denoted by a vertical green dashed line and the fatal tumor volume by a horizontal black dashed line.
Jpm 11 01036 g008
Figure 9. Kaplan -Meier survival curves for the control group (black line) and for the groups treated according to RTOG 0424 (red line) and the RT+TMZ simultaneous (green line) schedules. Differences between curves were statistically significant p-value (p) < 0.05.
Figure 9. Kaplan -Meier survival curves for the control group (black line) and for the groups treated according to RTOG 0424 (red line) and the RT+TMZ simultaneous (green line) schedules. Differences between curves were statistically significant p-value (p) < 0.05.
Jpm 11 01036 g009
Figure 10. In-silico studies comparing the KM survival curves for virtual patients with a followup period of 15 years treated according to RTOG 0424 vs simultaneous RT+TMZ. The pseudocolor plot shows the p-values for each trial as the number of patients per arm is increased. For a patient population greater than 720 (360 per arm) the probability of finding statistically significant differences between curves is 95 % .
Figure 10. In-silico studies comparing the KM survival curves for virtual patients with a followup period of 15 years treated according to RTOG 0424 vs simultaneous RT+TMZ. The pseudocolor plot shows the p-values for each trial as the number of patients per arm is increased. For a patient population greater than 720 (360 per arm) the probability of finding statistically significant differences between curves is 95 % .
Jpm 11 01036 g010
Table 1. Patients included in the study. Patient 36 was the only one treated with concurrent RT and TMZ. Patient 151 received RT and 11 years later TMZ. The meaning of the abbreviations used in the table can be found at the end of this document.
Table 1. Patients included in the study. Patient 36 was the only one treated with concurrent RT and TMZ. Patient 151 received RT and 11 years later TMZ. The meaning of the abbreviations used in the table can be found at the end of this document.
IDAge at DiagnosisSexHistologyTreatmentKi-67 Labeling Index (%) # RT Sessions # TMZ Cycles
1048FOTMZ5-15
2547MOTMZ5-4
3640MOTMZ & RT33513
4030FART134-
5754MOTMZ4-10
6628MARTNA30-
9134MARTNA30-
10533MOTMZ3-20
10734MART330-
10846MOTMZ3-5
12148MART335-
12427MOART335-
15158MORT37-
TMZ5-13
15946MNATMZNA-12
20351MOTMZ6-10
21358FOTMZ3-20
23444MART230-
Table 2. Parameters describing TMZ concentration in brain tumor tissue.
Table 2. Parameters describing TMZ concentration in brain tumor tissue.
ParameterDescriptionValueSources
cStandard dose per day150 mg/m 2 [26]
bPatient body surface1.9 m 2 (men)[41]
1.6 m 2 (women)[41]
γ Fraction of TMZ reaching2.1 · 10 6 / mL (men)[25]
the brain interstitium 2.5 · 10 6 / mL (women)[42]
t 1 / 2 TMZ half-life clearance time≃2 h[43]
λ Rate of TMZ decay0.3466/h[25]
Table 3. Parameter values obtained for the best fits of patients to Equation (1). Patient 151 was considered independently for each therapy. The units for the initial conditions and parameters are: for P 0 and Q 0 , cm 3 ; for Ki-67 0 % ; for ρ , β 1 , β 2 and β 3 , 1/day. Finally, α 1 and α 2 are measured in cm 3 / ( μ g day). * Estimated values.
Table 3. Parameter values obtained for the best fits of patients to Equation (1). Patient 151 was considered independently for each therapy. The units for the initial conditions and parameters are: for P 0 and Q 0 , cm 3 ; for Ki-67 0 % ; for ρ , β 1 , β 2 and β 3 , 1/day. Finally, α 1 and α 2 are measured in cm 3 / ( μ g day). * Estimated values.
ID#RT Sessions#TMZ Cycles P 0 Q 0 Ki-67 0 ρ β 1 β 2 β 3 S q S p α 1 α 2
Radiotherapy + Temozolomide
3635132.787.630.0231.50.0460.00020.450.552.33.9
Radiotherapy
4034-0.022.510.1510.0080.00110.320--
6630-0.6420.93 *0.061.50.040.00010.390--
9130-0.839.52 *0.0332.10.040.00080.160.04--
10730-1.342.830.0921.080.030.00780.860--
12135-0.18630.10.330.002 1 × 10 5 0.10.7--
12435-0.082.830.0130.870.020.0040.010--
1517-0.8627.930.011.30.040.010.940.09--
23430-1.26.920.0321.60.040.00240.780--
Temozolomide
10-152.243.650.013.50.1840.0008--6.014.31
25-41.732.650.0286.420.330.004--2020
57-101.845.140.010.230.0090.006--159.8
105-204.3141.730.0350.70.020.0018--2.86.5
108-50.415.230.0410.60.320.003--453
151-131.733.550.0120.530.020.0013--14.81
159-12266.13 *0.0123.50.10.004--116
203-101.727.260.0591.60.0980.021--8.20.37
213-200.723.130.0114.70.140.0005--15.215.1
Table 4. Real (in black, see Table 3) and virtual (in blue) values of the parameters related to therapies and the numbers of TMZ cycles and RT sessions received for the patients included. α 1 and α 2 are measured in cm 3 / ( μ g day). The final column corresponds to the overall survival gain (OSG) for each patient.
Table 4. Real (in black, see Table 3) and virtual (in blue) values of the parameters related to therapies and the numbers of TMZ cycles and RT sessions received for the patients included. α 1 and α 2 are measured in cm 3 / ( μ g day). The final column corresponds to the overall survival gain (OSG) for each patient.
ID#RT Sessions#TMZ Cycles S q S p α 1 α 2 OSG (Years)
1030150.280.476.014.31>10
253040.10.652020−1.2
3635130.450.552.33.9>10
4034120.32034.614.9>10
5730100.680.69159.88
6630120.39030.318.7
9130120.160.047.718.3>10
10530200.5102.86.54.7
10730120.86030.33.70.4
1083050.20.484532.9
12135120.10.710.310>10
12435120.01014.113.7>10
1517120.940.095.51.40.26
15130130.830.2414.811.9
15930120.180.17116>10
20330100.60.278.20.371.4
21330200.360.4115.215.1>10
23430120.7802.515.63.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ayala-Hernández, L.E.; Gallegos, A.; Schucht, P.; Murek, M.; Pérez-Romasanta, L.; Belmonte-Beitia, J.; Pérez-García, V.M. Optimal Combinations of Chemotherapy and Radiotherapy in Low-Grade Gliomas: A Mathematical Approach. J. Pers. Med. 2021, 11, 1036. https://doi.org/10.3390/jpm11101036

AMA Style

Ayala-Hernández LE, Gallegos A, Schucht P, Murek M, Pérez-Romasanta L, Belmonte-Beitia J, Pérez-García VM. Optimal Combinations of Chemotherapy and Radiotherapy in Low-Grade Gliomas: A Mathematical Approach. Journal of Personalized Medicine. 2021; 11(10):1036. https://doi.org/10.3390/jpm11101036

Chicago/Turabian Style

Ayala-Hernández, Luis E., Armando Gallegos, Philippe Schucht, Michael Murek, Luis Pérez-Romasanta, Juan Belmonte-Beitia, and Víctor M. Pérez-García. 2021. "Optimal Combinations of Chemotherapy and Radiotherapy in Low-Grade Gliomas: A Mathematical Approach" Journal of Personalized Medicine 11, no. 10: 1036. https://doi.org/10.3390/jpm11101036

APA Style

Ayala-Hernández, L. E., Gallegos, A., Schucht, P., Murek, M., Pérez-Romasanta, L., Belmonte-Beitia, J., & Pérez-García, V. M. (2021). Optimal Combinations of Chemotherapy and Radiotherapy in Low-Grade Gliomas: A Mathematical Approach. Journal of Personalized Medicine, 11(10), 1036. https://doi.org/10.3390/jpm11101036

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