Next Article in Journal
Quantifying the Displacement of Data Matrix Code Modules: A Comparative Study of Different Approximation Approaches for Predictive Maintenance of Drop-on-Demand Printing Systems
Previous Article in Journal
Segmentation of 4D Flow MRI: Comparison between 3D Deep Learning and Velocity-Based Level Sets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting the Tumour Response to Radiation by Modelling the Five Rs of Radiotherapy Using PET Images

1
INSERM UMR 1101 “LaTIM”, CEDEX 3, 29238 Brest, France
2
CERV, European Center for Virtual Reality, ENIB, CEDEX 3, 29238 Brest, France
3
The D-Lab, Department of Precision Medicine, GROW-School for Oncology, Maastricht University, 6211 LK Maastricht, The Netherlands
4
CHRU BREST, 29200 Brest, France
*
Author to whom correspondence should be addressed.
J. Imaging 2023, 9(6), 124; https://doi.org/10.3390/jimaging9060124
Submission received: 23 May 2023 / Revised: 10 June 2023 / Accepted: 13 June 2023 / Published: 20 June 2023
(This article belongs to the Topic Medical Image Analysis)

Abstract

:
Despite the intensive use of radiotherapy in clinical practice, its effectiveness depends on several factors. Several studies showed that the tumour response to radiation differs from one patient to another. The non-uniform response of the tumour is mainly caused by multiple interactions between the tumour microenvironment and healthy cells. To understand these interactions, five major biologic concepts called the “5 Rs” have emerged. These concepts include reoxygenation, DNA damage repair, cell cycle redistribution, cellular radiosensitivity and cellular repopulation. In this study, we used a multi-scale model, which included the five Rs of radiotherapy, to predict the effects of radiation on tumour growth. In this model, the oxygen level was varied in both time and space. When radiotherapy was given, the sensitivity of cells depending on their location in the cell cycle was taken in account. This model also considered the repair of cells by giving a different probability of survival after radiation for tumour and normal cells. Here, we developed four fractionation protocol schemes. We used simulated and positron emission tomography (PET) imaging with the hypoxia tracer 18F-flortanidazole (18F-HX4) images as input data of our model. In addition, tumour control probability curves were simulated. The result showed the evolution of tumours and normal cells. The increase in the cell number after radiation was seen in both normal and malignant cells, which proves that repopulation was included in this model. The proposed model predicts the tumour response to radiation and forms the basis for a more patient-specific clinical tool where related biological data will be included.

1. Introduction

Cancer is considered the leading cause of death in many developed countries [1], and it is a major public health problem around the world [2]. In the future, this disease will become the leading cause of morbidity and mortality in all regions of the world [3]. Cancer is a family of diseases characterized by abnormal cell growth, which damages the cell’s DNA. Radiotherapy represents one of the main curative approaches and plays a central role in cancer therapy. Indeed, more than 50% of patients receive radiotherapy at some point during their treatment [4]. Radiotherapy aims to deliver enough radiation to kill tumour cells but as low as possible to limit damage to surrounding normal cells [5]. That is why fractionated radiotherapy, in most cases, is the common technique used to treat cancer with radiation [6]. Indeed, dividing the total dose into several smaller ones over a period of several days reduces the toxicity of healthy cells [4].
The effectiveness of radiotherapy depends on several factors and especially on the total dose delivered. However, the response of the tumour to a given treatment can vary from one patient to another. Several underlying biological effects that could explain the variation in the biological response of tissues occurring during fractionated radiotherapy have been identified. Such biological principles can be summarized by a simple viewpoint called the five Rs of radiobiology: repair, repopulation, reoxygenation, redistribution in the cell cycle and radiosensitivity [7].
The heterogeneity of the tumour response in patients with the same type of cancer and the non-uniformity of tumour cells within the same individual tumour remain a key challenge for the cure of cancer [8]. Understanding and predicting the response of the tumour and normal tissues to treatment is therefore essential to provide useful information on the effectiveness of radiotherapy. The development of a clinical decision-support system based on predictive and prognostic data factors could help clinicians in their decision-making [9].
In recent years, there has been a significant amount of interest in the search for an optimal and personalized treatment plan [10,11,12]. Many models have been developed to study tumour growth and the radiotherapy response to radiation [13]. Most of these studies were based on simulations using mathematical models to describe the processes determining the response of cells and tissues to radiation. For example, the model proposed by Lind et al. [14] is based on the interaction of two Poisson processes and took into account two distinct types of cell damage, potentially repairable and conditionally repairable.
The effects of repair and repopulation have been modelled [15]. In this work, population dynamics and the cell number after radiation were simulated by coupling the linear-quadratic (LQ) model with a repopulation model. Furthermore, the effects of damage repair, cellular repopulation and redistribution in the cell cycle have been modelled via an analytical model describing the effect of radiation on cell inactivation [16].
Fakir et al. [17] presented a mathematical niche model for estimating tumour control probability (TCP) characterizing the tumour repopulation due to the proliferation of cells that survive after irradiation. A recent study by Yang et al. [18] focused on the process of tumour repopulation during radiotherapy and its detection methods using 18F-FLT PET. Badri et al. [19] developed a stochastic radiation-scheduling concept and incorporated inter-patient variability. This approach was based on a linear-quadratic model, including tumour proliferation, while tumour repopulation was based on a simple exponential law starting after a “kick-off” time [20].
Several experimental studies have shown that tumour cells in a hypoxic environment have a higher probability of surviving after irradiation than cells in a well-oxygenated situation. Oxygen availability is closely linked with the micro-vascular network. In this context, Grogan et al. [21] presented a hybrid multiscale model to predict the cellular response to radiotherapy. They used a cellular automaton approach of tumour growth and a model for oxygen transport from blood vessels. Other models have been developed to study the dynamics of hypoxia and the effect of reoxygenation in fractionated radiotherapy [22,23,24,25,26]. In general, most simulation methods remain theoretical and are rarely compared with actual clinical data. In addition, to our knowledge, no study on the modelling of the tumour response to radiotherapy has taken into account these five processes together. Proposed methods are rarely confronted with the actual tumour microenvironment of a given patient. Furthermore, as stated above, hypoxia and the partial pressure of oxygen (pO2) [27] should be taken into account when the response to radiotherapy [28] is studied. Although several methods capable of simulating pO2 exist [29], the proposed approaches often neglect the dynamic aspects of oxygen pressure during treatment [30].
A mathematical model for the tumour response that included the cell cycle effect in the context of lung cancer was presented by Jeong et al. [24] for understanding the dynamics of tumour cells during radiotherapy. A series of recent studies also included the effects of the cell cycle [31,32]. Other analytical models focused on the impact of intrinsic radiosensitivity on the tumour cell population [33,34,35].
It appears from these studies that an ideal predictive model should integrate all known radiobiological processes [34] to take into account the personal biological parameters of each patient, especially at the beginning of treatment. In the present study, the main objective was to evaluate the potential of incorporating the five Rs of radiobiology in a simulation process. We developed an analytical model that included the effects of repair, repopulation, redistribution, reoxygenation and radiosensitivity on the tumour response to radiotherapy.
This new general approach is illustrated by the implementation of tumour [18F] HX4 PET images as input data. This tracer is a hydrophilic variant of the 2-nitroimidazole class of radiotracers, which is used as a marker of hypoxia for the preclinical evaluation and validation of PET imaging [36,37].

2. Materials and Methods

2.1. Model of Tumour Response to Radiotherapy

The proposed model was developed to describe the evolution and tumour response at both macroscopic and microscopic scales. In this work, the macroscopic scale at the tumour level is represented by PET images. These patient-dependent images are used as the input data, and at this scale, the model consists of a series of N = N x × N y × N z voxels, where N x , N y , N z were the numbers of voxels in each direction for the PET image.
The total number of cells inside a voxel was constant, according to a fixed cell density µ (Table 1). Each voxel considered four populations of cells: tumour cells, capillary cells, normal cells and dead cells that may be present after irradiation. The number of capillary cells determined the state of oxygenation according to the vascular fraction vf (relative volume of capillary cells within a voxel). Since the amount of oxygen available may be heterogeneous in the cell microenvironment, a specific oxygen histogram was assigned to each voxel. The shape of the histogram depended on the vascular fraction of the voxel and was build according to the method described by Espinoza et al. [38]. Here, a series of six different oxygenation patterns was built, each one corresponding to a specific histogram, as depicted in Figure 1.
Once each voxel has been described according to these parameters, the proportions of cells evolve at each cell cycle time. At the microscopic scale, the model was based on the cellular division to handle tumour growth, according to classical phases of the cellular cycle:
  • G0 phase, also called quiescence phase where the cell has left the active cycle and has stopped dividing.
  • G1 phase, where the cell increases in size.
  • S phase, or synthesis phase, where DNA is duplicated.
  • G2 phase, preparing mitosis (synthesis of enzymes, etc.).
  • M phase or mitosis, which is the last phase of the cell cycle, when division occurs.
Based on the work of Wille et al. [39], we speculated that the duration of the cell cycle Tc is 24 h. The initial cell cycle phase distribution was assumed to be 14 h in G1, 6 h in S, 3 h in G2 and 1 h in M. In general, there are two checkpoints in the cell cycle, one for controlling transition from G1 to S and another for controlling transition from G2 to M. The transition from S to G2 was assumed to be automatic. A cell that was in none of these stages was considered in the quiescence phase (G0).
The cell cycle was arbitrarily divided into intervals of 1 h, which is also the duration of each step in the simulation. A cell that entered phase G1 automatically passed through the first intervals but underwent a transition test at the end of the phase. In the case of success, the cell continued its cycle until the next checkpoint (at G2/M).
Checkpoint tests were given using a random variable that follows Bernoulli’s law [32]. The total number of cells at a checkpoint was thus governed by a random variable that follows the binomial law. Given the large number of cells in the population concerned, this binomial distribution was then approximated based on a normal distribution [32].
Concerning tumour response modelling, a cell in the cellular cycle can die, survive and proliferate.
Table 1. List of parameters of the simulated tumour model.
Table 1. List of parameters of the simulated tumour model.
ParameterSymbolValue
Cell densityµ106/mm3 [40]
Cell cycle timeTc24 h
Upper asymptoteC1 [41]
Cell growth rateB0.075 [41]
Oxygen partial pressure p O 2 Calculated from histograms
Inflection point (the pO2 value at the point of maximum incline) M 26.3 [41]
Capillary cell proliferation (doubling time) t a 612 h [42]
Radiosensitivity coefficient α 0.273 Gy−1 [43]
Radiosensitivity coefficient β 0.045 Gy−2 [43]
Dose d 2 Gy
Maximum   O E R   value m 3 [44]
Oxygen partial tension at
O E R = m + 1 / 2
k 3 mmHg [45]
Half-life of dead cell resorption t r 168 h [22]

2.2. Model Components

A series of conventional radiobiological processes [32,41] was used as the main components of the model.
  • Proliferation of tumour cells: the probability of cell division depends on the current oxygenation status of the cell. Cell proliferation factor PF was calculated as follows:
    P F = C . exp { exp B . p O 2 M }
    where C is the superior asymptote, B the growth rate and M represents the pO2 value at the curve inflexion point.
  • Angiogenesis: the presence of hypoxia may induce angiogenesis [46]. In this case the fraction of the capillary cells was multiplied by the factor P F a :
    P F a = exp ln 2 t a × t i m e
    where ta is the doubling time for capillary cells and t i m e is the simulation time step (1 h).
  • Cell survival after irradiation SF:
    S F = exp ( α m d . O E R p O 2 β m ² d ² O E R p O 2 ² )
    Here, α and β are the radio-sensitivity coefficients, m is the maximum ratio, d is the dose and the O E R is given by the following:
    O E R p O 2 = m p O 2 + k p O 2 + k
    where k is the pO2 at half of the increase from 1 to m.
  • Resorption of dead tumour cells: tumour cells that have died after irradiation will be resorbed after a few days. This was reflected by the fraction of resorption:
    R F = 1 exp ( ln 2 t r × t i m e )
    Here, tr represents the half-life of dead cell resorption.
  • Cell replacement: if the number of cells per voxel decreased after the resorption of dead cells, it was necessary to redistribute some new cells in the voxel to maintain cell density. For the sake of simplicity, resorbed cells were arbitrarily replaced by normal cells or capillary ones, depending on the oxygenation model.

2.3. Factors Influencing Tumour Response According to the 5 Rs

In this study, the developed model was obtained by adding the 5 Rs of radiobiology to the processes described above. In the literature, the values of these different parameters differ depending on the type of cancer. With regard to our study, the parameter values relate to the cancer “rhabdomyosarcoma”. Therefore, to determine each parameter, as well as to combine the model components, we relied on previous studies as explained below.

2.3.1. Repair

The use of radiotherapy as a treatment modality exposes cells to ionizing radiation, which damages DNA. Cell survival after irradiation depends on the cell’s ability to repair itself and on the type of lesion. In fact, if DNA damage is irreparable, cells activate death programs [47]. In the case of normal tissue cells, if they can repair the damage, it can only be beneficial for the treatment result. In contrast, in the case of tumour cells, post-radiation survival allows tumour cells to proliferate [48]. The success of radiotherapy therefore depends on the extent of damage in the exposed tumour tissue [47].
In fractionated radiotherapy, by dividing the total dose into a set of fractions, repair produces increased cell survival by allowing cells to recover after the individual radiation dose and consequently proliferation between fractions. Generally, the linear-quadratic model is often used in radiobiology to describe the response to radiation [1]. It represents a mathematical formula linking cell survival and the radiation dose and depends on two parameters, α and β. The parameter α determines the initial slope of the cellular survival curve, and β represents the accumulation of sublethal damage.
Dose splitting allows time for healthy tissue to recover and repair itself between sessions, which is not the case for tumour cells as their repair system is disrupted by nature.
In this study, the density of surviving normal cells after irradiation is given by the following:
S F N = ( α d ( 1 + d β α ) )
The density formula used in our model was inspired by the work presented by Joiner et al. [1]. Dead normal cells will be replaced as dead tumour cells as described above in “Cell replacement”. The control of DNA repair is closely linked to the redistribution of the cell in the cellular cycle [49]. In fact, checkpoints during the cell cycle ensure that DNA is intact before DNA replication and cell division.

2.3.2. Redistribution in the Cell Cycle

During irradiation, cells can be anywhere in the cell cycle, at a given phase, among G1, S, G2 or M.
Cells are characterized by various radiosensitivities depending on their location in the cell cycle. Cells in phase S are known to be more radioresistant, while cells in other phases are relatively more sensitive to radiation [50]. Thus, the radiation dose delivered for an asynchronous population of cells will probably induce more deaths for cells that are in the sensitive phase. Consequently, cells that survive are those that were in resistant phases. As a consequence, the surviving population is partially synchronized [51].
Damage to DNA after irradiation will block cell progression at cell cycle checkpoints [52] causing cell accumulation at this point. This behaviour is expected to prevent cells from entering into the division phase with fatal damage [53]. Over time, surviving cells will continue to evolve in the cycle. If a second dose of radiation is delivered sometime later, some of these cells will have left the resistant phase and will be in a more sensitive phase, which will allow them to be killed more easily. In our study, we proposed to add a weighting coefficient described by Joiner et al. [1] to manage the sensitivity of cells to their phases before irradiation.

2.3.3. Repopulation

Conventional fractionation consists of delivering a small dose (2 Gy) daily on a weekday during several weeks. This protocol allows normal tissues to repair potential radiation damage between two consecutive fractions. However, surviving tumour cells can also repopulate between fractions. Repopulation plays an important role in the outcome of radiotherapy, since surviving tumour cells can generate new cells that were not present at the beginning of the treatment.
Repopulation may have a significant role when the overall treatment time is long, as seen in classical fractionation protocols. However, there are strategies that can reduce the effect of repopulation through accelerated fractionation, which reduces the overall treatment time by delivering two fractions per day instead of one. In addition, an accelerated repopulation of malignant cells is often observed in a series of common cancers [54]. These types of tumours show accelerated repopulation after two or three weeks of treatment.
In this work, we chose to vary the doubling time Tp, defined as the time required for the tumour cell population to double without natural loss. As an example, during the first two weeks Tp is set to 1200 h and becomes 120 h until the end of treatment [1,55].

2.3.4. Reoxygenation

Local oxygen tension plays an important role in the radiotherapy response [54]. Several studies showed that the poor prognosis of radiotherapy might be correlated with tumour hypoxia [56,57]. Thomlinson and Gray [58] have suggested that human lung cancer contains hypoxic areas, with necrotic areas closely related to the size of the tumour. In fact, when the tumour size is greater than 1 mm, cells located at a distance from the vessels have reduced access to oxygen, and the central part of the tumour becomes necrotic. On the contrary, when the diameter of tumours is less than 1 mm, they were found to be fully oxygenated [1]. In addition, the effect of ionizing radiation targeting DNA is stabilized in the presence of oxygen. In contrast, in hypoxic environments, tumour cells are two-to-three times more resistant than well-oxygenated cells. Moreover, when radiation is separated into multiple fractions delivered over weeks, cells that are relatively resistant due to hypoxia at the beginning of the treatment may become reoxygenated.
Reoxygenation is the phenomenon of hypoxic tumours becoming oxygenated again during radiotherapy treatment. Numerous studies have shown that the reoxygenation of tumours may occur only 24 to 72 h after irradiation [59]. Reoxygenation is caused by several mechanisms, such as reopening temporarily closed vessels or decreasing the respiration of fatally damaged cells. Other mechanisms need longer intervals to occur, such as cell death due to mitotic catastrophe or ischemia leading to tumour shrinkage and a decreased distance between capillaries and tumour cells that allows oxygen to reach the hypoxic zone.
In this work, oxygen levels varied in time and space. The vascular fraction vf varied during the treatment and changed from one voxel to another. In the model, vf depended on the proportion of capillary cells in voxels. Initially, every voxel contained 80% of tumour cells, and in the remaining 20%, 96.4% of normal cells and 3.6% of capillary cells were present [60]. This distribution would remain valid throughout the simulation, with dead cells replaced by capillary and normal cells in the proportions of 3.6% and 96.4%, respectively. The reoxygenation of malignant cells during radiotherapy can increase their radiosensitivity.

2.3.5. Radiosensitivity

Tumours respond differently to radiation therapy. This variable response is correlated with the intrinsic radiosensitivity of cells. Radiosensitivity was added to the four Rs in 1989 by Steel et al. [61]. This new member of the Rs of radiobiology highlights the fact that there is an intrinsic radiosensitivity in different types of cells.
Tumour cell radiosensitivity is a significant prognostic factor for the overall response of tumours and radiotherapy outcomes. Some studies were interested in trying to determine and predict this inter-individual variability. In theory, it would be interesting to be able to determine the intrinsic radiosensitivity of tumour cells, as well as the radiosensitivity of healthy tissues, for each patient. Radiosensitivity, which is the relative susceptibility of cells to ionizing radiation, is usually indicated by the parameters of the linear-quadratic equations. For this reason, we used specific α and β parameters for each cell type, either healthy or tumoural.

2.4. Simulations

In our study, four fractionation schemes were designed and simulated to quantify the impact on the tumour response to radiotherapy (Figure 2).
In Figure 2, vertical arrows indicate the days when the radiotherapy session takes place. The numbers 1 to 5 represent weeks. The first fractionation scheme studied is conventional or standard fractionation (SP: standard protocol), with a 2 Gy dose delivered 5 days per week for 5 weeks. The second fractionation scheme is sometimes referred to as DAHANCA (DP), in which a dose of 2 Gy is administered 6 days a week instead of 5, leading to an overall treatment time of 4 weeks instead of 5. The third fractionation scheme is a personal scheme (PP) in which a dose of 2 Gy is administered for 25 days without any interruption (weekends are included). The fourth fractionation scheme shown in Figure 2 is the CHART protocol. During this protocol, 36 fractions of 1.5 Gy over 12 consecutive days are administered. The treatment consisted of 3 fractions per day every 6 h. The tumour response was simulated for each fractionation scheme using simulated tumour and then PET images.

2.5. Data

In this study, artificially created tumours with simple geometries were used to demonstrate the main features of the model and to visualize the impact of the different biological effects. The virtual tumour had a diameter of approximately 3 cm in accordance with the work of Espinoza et al. ([60]).
Then, for validation, we used [18F] HX4 PET images from an image database generated by our previous study [62], which are considered a reliable tool for the detection of hypoxic tumour regions [37,63,64]. Data (images) from four mice with rhabdomyosarcoma tumours, with two control mice that received NaCl while two others received TH-302 for 4 days, have been used in this study. A [18F] HX4 PET image was acquired before treatment and a second one was acquired on day four of the treatment. Irradiation was delivered on day 3. During the treatment, animals were exposed to modified oxygen concentrations. Two mice were treated under carbogen (95% oxygen, 5% CO2) and the other two were treated under reduced oxygen breathing (7% oxygen):
  • Mouse 1: NaCl + 95% oxygen;
  • Mouse 2: TH-302 + 95% oxygen;
  • Mouse 3: NaCl + 7% oxygen;
  • Mouse 4: TH-302 + 7% oxygen.
The model input data are PET images and a distribution of the partial pressure of oxygen. The model returns simulated PET images. Throughout this study, we hypothesized that the number of tumour cells is proportional to glucose consumption as described elsewhere [28,41].
To calculate the initial distribution of cells in each voxel, we searched for the voxel with the maximum intensity (imax). Then, we divided the values of each voxel by the max. Then, µt = 80% of the tumour cells were assigned to each voxel, and in the remaining 20%, µn = 96.4% were normal cells and µc = 3.6% capillary cells.
In a voxel, we have 106 cells/mm3 [60], and denoting imax as the maximum intensity among all the intensity values of the voxels, the initial number of tumour cells in each voxel is given by the following:
N t u m = i v o x i m a x × V v o x × μ t × 10 6
where Vvox is the volume of a voxel and ivox is the intensity of the voxel.
Therefore, a new PET image can be generated during the simulation using the updated value of the cell number value to calculate the intensities in this new image.

2.6. Tumour Control Probability (TCP)

TCP (Tumour Control Probability) is a model that predicts and quantifies the biologic radiation response of tumours, the survival. The most used model is based on Poisson statistics describing the distribution of the surviving tumour cells. In general, the LQ model is used for calculating cell survival and it allows for the quantification of TCP for tumours.
Here, a TCP curve was generated by simulating the tumour response to radiotherapy. We consider a tumour controlled if all tumour cells were killed.

3. Results

The presented model allows us to simulate the tumour response under radiotherapy. In Figure 3, the temporal evolution of the total number of tumour cells is presented. During the first 72 h, the increase in the number of tumour cells was due to the absence of radiotherapy. After this initial phase, the effect of radiotherapy becomes visible. A daily decrease in tumour cells followed by a recovery was observed due to the repopulation of tumour cells between two radiations.
In Figure 4, tumour growth images are presented. Qualitatively, the visual comparison between real and simulated images shows that there is real similarity between them. The simulated image was obtained from the real images acquired before the beginning of treatment. Mice were under carbogen. Although these images were obtained after a single irradiation, we can see the effectiveness of radiotherapy, which is not the case in Figure 5.
In Figure 5, it can be noted that there is a certain resistance of the cells, which is quite logical because the mice were under limited oxygenation. In the case of mice under breathing 7% conditions, we imposed, at the beginning of the simulation, a hypoxic histogram to reproduce the conditions under which mice were treated. In addition, in the other case (under carbogen), we imposed a well-oxygenated histogram.
In the case of mouse number 3, the effect of radiotherapy is not visible. In the case of mouse 4, we can see a slight change. This result is also visible in the case of mouse 4, despite treatment with TH-302, an anti-tumour agent that is activated in a hypoxic environment.

3.1. Repair

In our study, we modelled the response of normal cells to radiotherapy with the linear quadratic model. In Figure 6, the temporal evolution of normal cells is presented. The blue curve represents the evolution of normal cells using the linear quadratic model, which is not the case in the red curve.
During the first 72 h, the decrease in the total number of normal cells is due to the absence of radiotherapy. In fact, when radiotherapy is not administrated, the number of tumour cells per voxel increased. Since there is a fixed cell density in a voxel, tumour cells continued to proliferate at the expense of other cell types. Then, when radiotherapy was administered, the number of dead tumour cells following irradiation gave way to other cell types. Therefore, the number of normal cells increased after 72 h. Then, a slight change in this number was noted between two irradiations caused by cellular sensitivity to irradiation and cellular repopulation.

3.2. Redistribution in the Cell Cycle

The distribution of cells in the cell cycle affects the outcome of radiation therapy. A weighting coefficient was added to manage the sensitivity of the cells according to their location before irradiation (blue curve Figure 7), which is not the case in the orange curve.
According to the work of Joiner et al. [1], we chose to modify the parameters during each irradiation. The aim was to see the impact of this variation on the result of the irradiation after one week of treatment. We observed in Figure 8 that from the third day, there was a difference between the irradiation administered with variations in the parameters (Figure 8a) and without variation (Figure 8b). It was deduced that variation of the parameters allows for the more rapid control of the tumour.
Figure 9 confirms the difference seen in Figure 8. In Figure 9a, parameters were varied during every irradiation (for each irradiation per day, we have a new coefficient), and in Figure 9b, the value of these parameters was constant during the treatment. In this figure, we have drawn an interactive 3D surface plot. These plots display a three-dimensional graph of the intensities of voxels in the image. In our study, we considered the closer the voxel intensity is to the black colour, the more we controlled the tumour. Indeed, the darkest areas are the areas where there is more oxygen. Thus, fewer tumour cells are in these areas.
In Table 2, a percentage of cells in each phase of the cell cycle has been shown. It is observed that the percentage of cells in G1 after 24 h increases slightly and remains stable in phase M. On the other hand, the percentage of cells in phase S decreased but it increased in phase G. This observed evolution corresponds to experimental data presented in the study [65].
In Figure 10, a percentage of cells in each cell cycle phase before and after radiation therapy is shown. It was observed that there was a decrease in the percentage of cells after radiotherapy in phases G1, S and M, which was not the case for phase G2. Compared to the initial condition, for a dose of 2 Gy, there was a large increase in cells in phase G2 after radiotherapy. The blockage of cells in G2 following radiation exposure can be attributed to the activation of cellular repair mechanisms. As a result, the percentage of cells in M phase decreased at the same time. The percentages of G1 and S cells also declined following irradiation [65].
In Figure 11, evolution of the cell cycle phases over time is presented. The irradiation takes place at t = 20 h. It is observed that following radiotherapy, the percentage of cells in G1 and S decreases and increases in G2. After this disturbance (t = 30 h), the percentage in G1 and S increases and that in G2 decreases.

3.3. Repopulation

During the first two weeks of radiotherapy, Tp will be equal to 1200 h and then change to 120 h until the end of the treatment. In Figure 12, the blue curve represents tumour development considering repopulation, but in the red curve, cell repopulation is not considered.

3.4. Reoxygenation

Here, we considered the spatiotemporal evolution of oxygen during radiotherapy treatment.
In Figure 13, we compared the tumour response with the spatiotemporal evolution of oxygenation (blue curve) and constant oxygenation (red curve). We can see that our model is sensitive to changes in the rate of the partial pressure of oxygen.

3.5. Radiosensitivity

In this study, the radiosensitivity parameters were varied to see their impact on the tumour response. The values of these parameters have been reported by three different studies. At the start, we used linear quadratic model parameters α = 0.339 and β = 0.067, which gave a α/β ratio equal to 5.06 Gy [66]. Second, the tumour response was simulated with α = 0.04 and β = 0.0089 (ratio α/β = 4.94 Gy) [32]. Finally, we worked with an α/β ratio = 6.06 Gy [1].
By modifying the α and β radiosensitivity parameters, a different tumour response is obtained (Figure 14). This shows that our model considers intrinsic radiosensitivity.
To demonstrate the main features of the model and visualize the impact of different fractionation patterns on tumour response, an artificially created virtual tumour with simple geometries was used in this study as a first step. Then, to validate the results obtained with the simulated tumour, the tumour response was simulated using a real image of a rhabdomyosarcoma tumour.

3.6. Tumour Response to Different Fractionation Schemes

3.6.1. Simulated Tumour Response

In this study, four fractionation schemes were simulated. The tumour response to radiotherapy was then studied for each treatment regimen. It can be seen in Figure 15 that the personal protocol allows the tumour to be controlled more quickly than the other protocols. It is closely followed by the CHART protocol.
The following figure confirms the result observed in Figure 15. The simulation of the response of the simulated tumour to radiotherapy leads to different responses for the different treatment regimens (see Figure 16a). The PP and CP protocols provide faster tumour control than the conventional regimen and the DAHANCA protocol.
The tumour response was quantified in terms of cell survival and TCP curves (Figure 16b). The TCD50 (the dose required to control 50% of tumours locally) was calculated. The TCD50 values (in Gy) of the TCP curves are 17.7, 17.9, 20.1 and 23.4 for PP, CP, SD and DP, respectively.

3.6.2. Tumour Response to Radiotherapy

In this study, we used [18F] HX4 PET images of mice with rhabdomyosarcoma as input for our simulation model. Figure 17 demonstrates that with the CHART protocol, we have better local control of the tumour.
In Figure 18, the TCP curve is plotted. Three fractionation protocols have been brought together: CHART, staff and DAHANCA. According to this curve, these three protocols are more efficient in terms of local tumour control than the classic or standard model. As with the simulated images, the values of TCD50 were calculated. The values were found to be 13.23, 13.24, 13.72 and 18.18, respectively for CP, PP, DP and SP.
Figure 19 illustrates tumour development for each fractionation scheme. The images obtained on the seventh day after the start of treatment show that the CHART protocol has better tumour control, especially compared to the conventional protocol.

4. Discussion

In the context of the numerical simulation of tumour growth and response to radiotherapy, we proposed to combine two series of biological processes. The first set includes the main processes often studied independently so far, such as tumour cell proliferation, angiogenesis, cell survival after irradiation, resorption of dead tumour cells and cell replacement. The second set includes the main concepts of what is referred to as the five Rs of radiobiology. These five Rs are cell cycle redistribution, reoxygenation, intrinsic radiosensitivity, cell repopulation and repair. This is the first approach where these concepts are integrated into a unified environment, using both the microscopic cell cycle scale and the macroscopic population scale. To our knowledge, no other study on the modelling of the tumour response to radiotherapy considered these processes altogether. Most previous studies indeed simulated these processes separately [67], for example without considering reoxygenation despite its importance in the response to treatment. The approach presented here has several important advantages, such as versatility and the possibility of using it in different contexts with all kinds of input data, for example patient-specific biological data.
More specifically, we have shown that our approach responds correctly to the modification of input parameters, through application to a virtual cell population. We have also shown that it is possible to quantify the impact of fractionation on the response, by studying a theoretical population. An example of a practical evaluation is also proposed, using, as input data, a set of preclinical [18F] HX4 PET images acquired under different oxygen conditions. Images obtained after simulation were consistent with the real post-irradiation images. Moreover, our model has been used to test [18F] FDG PET images (non-hypoxia tracer) of patients with rectal cancer, and the results are promising (data not shown).
In this study, we investigated the impact of four different fractionation schemes on the evolution of a cell population. The objective was to find the protocol that allowed for better control of the tumour while limiting healthy tissue damage. The results showed that the reduction in treatment time has a considerable effect on local tumour control. In the cases of simulated and real tumours, the CHART protocol showed better tumour control, especially when compared with the conventional protocol. In the case of rectal cancer, studies have shown that the CHART protocol, as a preoperative treatment for rectal cancer, is feasible and appears to be associated with low acute and late toxicity. In the study of S. Brooks et al. [68], the authors found that local control was encouraging and therefore justified a more in-depth evaluation of this protocol. Thus, they proposed the incorporation of chemotherapy, which could also be considered as an option of investigation in the future.
Images used in this work were images of mice treated with NaCl or TH-302. TH-302 is an investigational hypoxia-activated prodrug used to treat cancer. TH-302 is activated only in hypoxic regions where the oxygen level is very low. It was developed to target the hypoxic zones existing in tumours. Generally, it is used combination with chemotherapy and shows rather satisfactory results. Recently, in the study of Peeters et al. [62], the authors studied the efficacy of the combination treatment of TH-302 and radiotherapy. In perspective, it would be better to model the TH-302 effect on the tumour and combine its effect with that of radiotherapy and then compare the simulated images with the images of the mice obtained.
As future work, it will be interesting to simulate different dose-painting strategies targeting hypoxic regions in tumours. Indeed, a non-uniform response can be observed within the tumour due to the heterogeneity of oxygenation. This point is consistent with the fact that hypoxic cells are more radioresistant than well-oxygenated ones, which can negatively affect the radiotherapy result. Simulation could provide insights about the non-uniform dose distribution required to destroy an equal proportion of cells within a heterogeneous tumour. In recent work, research demonstrated that radiotherapy induces an immune response. During radiotherapy, some immunosuppressive barriers might increase. This process can lead to tumour resistance. Boustani et al. [69] reviewed the effect of fractionation and the dose on the anti-tumour immune response. They propose a sixth R to the five Rs of radiobiology. This sixth one is reactivation of the immune system. Therefore, on the horizon, it would be interesting to find a way to model this sixth R.

5. Conclusions

For radiotherapy, the treatment outcome is determined by several factors. Some factors diminish the effect of radiation therapy, for example, the repopulation of tumour cells and repair from sublethal damage. Some factors increase the local control of the cancer, such as processes of reoxygenation and the redistribution of tumour cell into more sensitive phases of the cell cycle. The impact of the individual factor varies across different tissues, but it is important to consider all of these factors when developing a successful radiotherapy model. In the presented work, we tried to develop a model considering all five Rs of radiobiology. The model developed could be expanded in the future to include patient biological data for a truly specific clinical tool.

Author Contributions

Study concept and design, R.H. and N.B.; PET image acquisition, L.J.D. and P.L.; acquisition, analysis, interpretation of data and drafting of the manuscript, R.H.; review of the manuscript, S.A., P.R., L.G., D.V. and N.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

IRB for this study was waived as only images from a database were used. The image database was generated by the study of reference [63].

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be provided if required.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Joiner, M.; Van der Kogel, A. Basic Clinical Radiobiology, 4th ed.; Hodder Arnold: London, UK, 2009; ISBN 978 0 340 929 667. [Google Scholar]
  2. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer Statistics, 2018. CA Cancer J. Clin. 2018, 68, 7–30. [Google Scholar] [CrossRef] [PubMed]
  3. Bray, F.; Jemal, A.; Grey, N.; Ferlay, J.; Forman, D. Global Cancer Transitions According to the Human Development Index (2008–2030): A Population-Based Study. Lancet Oncol. 2012, 13, 790–801. [Google Scholar] [CrossRef] [PubMed]
  4. Baskar, R.; Lee, K.A.; Yeo, R.; Yeoh, K.W. Cancer and Radiation Therapy: Current Advances and Future Directions. Int. J. Med. Sci. 2012, 9, 193–199. [Google Scholar] [CrossRef] [Green Version]
  5. Deloch, L.; Derer, A.; Hartmann, J.; Frey, B.; Fietkau, R.; Gaipl, U.S. Modern Radiotherapy Concepts and the Impact of Radiation on Immune Activation. Front. Oncol. 2016, 6, 141. [Google Scholar] [CrossRef] [Green Version]
  6. Nahum, A.E. The Radiobiology of Hypofractionation. Clin. Oncol. 2015, 27, 260–269. [Google Scholar] [CrossRef]
  7. Klement, R.J. The Influence of Ketogenic Therapy on the 5 R’s of Radiobiology. Int. J. Radiat. Biol. 2017, 95, 394–407. [Google Scholar] [CrossRef]
  8. Peitzsch, C.; Kurth, I.; Ebert, N.; Dubrovska, A.; Baumann, M. Cancer Stem Cells in Radiation Response: Current Views and Future Perspectives in Radiation Oncology. Int. J. Radiat. Biol. 2019, 95, 900–911. [Google Scholar] [CrossRef] [PubMed]
  9. Lambin, P.; Van Stiphout, R.G.P.M.; Starmans, M.H.W.; Rios-Velazquez, E.; Nalbantov, G.; Aerts, H.J.W.L.; Roelofs, E.; Van Elmpt, W.; Boutros, P.C.; Granone, P.; et al. Predicting Outcomes in Radiation Oncology-Multifactorial Decision Support Systems. Nat. Rev. Clin. Oncol. 2012, 10, 27–40. [Google Scholar] [CrossRef] [Green Version]
  10. Mali, S.B. Adaptive Radiotherapy for Head Neck Cancer. J. Maxillofac. Oral Surg. 2016, 15, 549–554. [Google Scholar] [CrossRef] [Green Version]
  11. Hunt, A.; Hansen, V.N.; Oelfke, U.; Nill, S.; Hafeez, S. Adaptive Radiotherapy Enabled by MRI Guidance. Clin. Oncol. 2018, 30, 711–719. [Google Scholar] [CrossRef] [Green Version]
  12. Belshaw, L.; Agnew, C.E.; Irvine, D.M.; Rooney, K.P.; McGarry, C.K. Adaptive Radiotherapy for Head and Neck Cancer Reduces the Requirement for Rescans during Treatment Due to Spinal Cord Dose. Radiat. Oncol. 2019, 14, 189. [Google Scholar] [CrossRef]
  13. Belfatto, A.; White, D.A.; Zhang, Z.; Zhang, Z.; Cerveri, P.; Baroni, G.; Mason, R.P. Mathematical Modeling of Tumor Response to Radiation: Radio-Sensitivity Correlation with BOLD, TOLD, Δr1 and Δr2∗ Investigated in Large Dunning R3327-AT1 Rat Prostate Tumors. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, EMBS, Milan, Italy, 25–29 August 2015; Volume 2015, pp. 3266–3269. [Google Scholar]
  14. Lind, B.K.; Persson, L.M.; Edgren, M.R.; Hedlöf, I.; Brahme, A. Repairable–Conditionally Repairable Damage Model Based on Dual Poisson Processes. Radiat. Res. 2003, 160, 366–375. [Google Scholar] [CrossRef] [PubMed]
  15. Farayola, M.F.; Shafie, S.; Siam, F.M.; Khan, I. Mathematical Modeling of Radiotherapy Cancer Treatment Using Caputo Fractional Derivative. Comput. Methods Programs Biomed. 2020, 188, 105306. [Google Scholar] [CrossRef] [PubMed]
  16. Zaider, M.; Wuu, C.S.; Minerbo, G.N. The Combined Effects of Sublethal Damage Repair, Cellular Repopulation and Redistribution in the Mitotic Cycle. I. Survival Probabilities after Exposure to Radiation. Radiat. Res. 1996, 145, 457. [Google Scholar] [CrossRef]
  17. Fakir, H.; Hlatky, L.; Li, H.; Sachs, R. Repopulation of Interacting Tumor Cells during Fractionated Radiotherapy: Stochastic Modeling of the Tumor Control Probability. Med. Phys. 2013, 40, 121716. [Google Scholar] [CrossRef] [PubMed]
  18. Yang, J.; Yue, J.B.; Liu, J.; Yu, J.M. Repopulation of Tumor Cells during Fractionated Radiotherapy and Detection Methods (Review). Oncol. Lett. 2014, 7, 1755–1760. [Google Scholar] [CrossRef] [Green Version]
  19. Badri, H.; Watanabe, Y.; Leder, K. Optimal Radiotherapy Dose Schedules under Parametric Uncertainty. Phys. Med. Biol. 2016, 61, 338–364. [Google Scholar] [CrossRef]
  20. Saberian, F.; Ghate, A.; Kim, M. Optimal Fractionation in Radiotherapy with Multiple Normal Tissues. Math. Med. Biol. 2015, 33, 211–252. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Grogan, J.A.; Markelc, B.; Connor, A.J.; Muschel, R.J.; Pitt-Francis, J.M.; Maini, P.K.; Byrne, H.M. Predicting the Influence of Microvascular Structure on Tumor Response to Radiotherapy. IEEE Trans. Biomed. Eng. 2017, 64, 504–511. [Google Scholar] [CrossRef]
  22. Harting, C.; Peschke, P.; Borkenstein, K.; Karger, C.P. Single-Cell-Based Computer Simulation of the Oxygen-Dependent Tumour Response to Irradiation. Phys. Med. Biol. 2007, 52, 4775–4789. [Google Scholar] [CrossRef]
  23. Kempf, H.; Bleicher, M.; Meyer-Hermann, M.; Rocha, S. Spatio-Temporal Dynamics of Hypoxia during Radiotherapy. PLoS ONE 2015, 10, e0133357. [Google Scholar] [CrossRef]
  24. Jeho, J.; Jung, H.O.; Jan-Jakob, S.; Jose, B.; Jeffrey, D.B.; Andrew, N.F.; Shyam, S.R.; Joseph, O.D. Modeling the Cellular Response of Lung Cancer to Radiation Therapy for a Broad Range of Fractionation Schedules. Clin. Cancer Res. 2018, 176, 139–148. [Google Scholar] [CrossRef]
  25. Petersson, K.; Adrian, G.; Butterworth, K.; McMahon, S.J. A Quantitative Analysis of the Role of Oxygen Tension in FLASH Radiation Therapy. Int. J. Radiat. Oncol. Biol. Phys. 2020, 107, 539–547. [Google Scholar] [CrossRef] [PubMed]
  26. Gago-Arias, A.; Aguiar, P.; Espinoza, I.; Sánchez-Nieto, B.; Pardo-Montero, J. Modelling Radiation-Induced Cell Death and Tumour Re-Oxygenation: Local versus Global and Instant versus Delayed Cell Death. Phys. Med. Biol. 2016, 61, 1204–1216. [Google Scholar] [CrossRef]
  27. Chen, Y.; Cairns, R.; Papandreou, I.; Koong, A.; Denko, N.C. Oxygen Consumption Can Regulate the Growth of Tumors, a New Perspective on the Warburg Effect. PLoS ONE 2009, 4, e7033. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Mi, H.; Petitjean, C.; Dubray, B.; Vera, P.; Ruan, S. Prediction of Lung Tumor Evolution during Radiotherapy in Individual Patients with Pet. IEEE Trans. Med. Imaging 2014, 33, 995–1003. [Google Scholar] [CrossRef] [PubMed]
  29. Grimes, D.R.D.; Kannan, P.; Warren, D.R.; Markelc, B.; Bates, R.; Muschel, R.; Partridge, M.; Gray, L.; Conger, A.; Ebert, M.; et al. Estimating Oxygen Distribution from Vasculature in Three-Dimensional Tumour Tissue. J. R. Soc. Interface 2016, 13, 20160070. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Powathil, G.G.; Gordon, K.E.; Hill, L.A.; Chaplain, M.A.J. Modelling the Effects of Cell-Cycle Heterogeneity on the Response of a Solid Tumour to Chemotherapy: Biological Insights from a Hybrid Multiscale Cellular Automaton Model. J. Theor. Biol. 2012, 308, 1–19. [Google Scholar] [CrossRef]
  31. Dawson, A.; Hillen, T. Derivation of the Tumour Control Probability (TCP) from a Cell Cycle Model. Comput. Math. Methods Med. 2006, 7, 121–141. [Google Scholar] [CrossRef] [Green Version]
  32. Apeke, S.; Gaubert, L.; Boussion, N.; Lambin, P.; Visvikis, D.; Rodin, V.; Redou, P. Multi-Scale Modelling and Oxygen Impact on Tumor Temporal Evolution: Application on Rectal Cancer during Radiotherapy. IEEE Trans. Med. Imaging 2018, 37, 871. [Google Scholar] [CrossRef]
  33. Alfonso, J.C.L.; Berk, L. Modeling the Effect of Intratumoral Heterogeneity of Radiosensitivity on Tumor Response over the Course of Fractionated Radiation Therapy. Radiat. Oncol. 2019, 14, 88. [Google Scholar] [CrossRef] [Green Version]
  34. Hedman, M.; Bjrk-Eriksson, T.; Mercke, C.; West, C.; Hesselius, P.; Brodin, O. Comparison of Predicted and Clinical Response to Radiotherapy: A Radiobiology Modelling Study. Acta Oncol. 2009, 48, 584–590. [Google Scholar] [CrossRef]
  35. Hedman, M.; Björk-Eriksson, T.; Brodin, O.; Toma-Dasu, I. Predictive Value of Modelled Tumour Control Probability Based on Individual Measurements of in Vitro Radiosensitivity and Potential Doubling Time. Br. J. Radiol. 2013, 86, 20130015. [Google Scholar] [CrossRef] [Green Version]
  36. Betts, H.M.; O’Connor, R.A.; Christian, J.A.; Vinayakamoorthy, V.; Foweraker, K.; Pascoe, A.C.; Perkins, A.C. Hypoxia Imaging with [18 F]HX4 PET in Squamous Cell Head and Neck Cancers: A Pilot Study for Integration into Treatment Planning. Nucl. Med. Commun. 2019, 40, 73–78. [Google Scholar] [CrossRef] [PubMed]
  37. Dubois, L.J.; Lieuwes, N.G.; Janssen, M.H.M.; Peeters, W.J.M.; Windhorst, A.D.; Walsh, J.C.; Kolb, H.C.; Öllers, M.C.; Bussink, J.; Van Dungeon, G.A.M.S.; et al. Preclinical Evaluation and Validation of [ 18F]HX4,a Promising Hypoxia Marker for Imaging. Proc. Natl. Acad. Sci. USA 2011, 108, 14620–14625. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Espinoza, I.; Peschke, P.; Karger, C.P. A Model to Simulate the Oxygen Distribution in Hypoxic Tumors for Different Vascular Architectures. Med. Phys. 2013, 40, 081703. [Google Scholar] [CrossRef] [PubMed]
  39. Wille, J.; Pittelkow, M.; Shipley, G.; Scott, R. Integrated Control of Growth and Differentiation of Normal Human Prokeratinocytes Cultured in Serum-Free Medium: Clonal Analyses, Growth Kinetics, and Cell Cycle Studies. J. Cell. Physiol. 1984, 121, 31–44. [Google Scholar] [CrossRef] [PubMed]
  40. Steel, G.G. Basic Clinical Radiobiology, 3rd ed.; George, G., Ed.; Hodder Arnold: London, UK, 2002; ISBN 9780340807835. [Google Scholar]
  41. Titz, B.; Jeraj, R. An Imaging-Based Tumour Growth and Treatment Response Model: Investigating the Effect of Tumour Oxygenation on Radiation Therapy Response. Phys. Med. Biol. 2008, 53, 4471. [Google Scholar] [CrossRef] [Green Version]
  42. Denekamp, J.; Hobson, B. Endothelial-Cell Proliferation in Experimental Tumours. Br. J. Cancer 1982, 46, 711–720. [Google Scholar] [CrossRef] [Green Version]
  43. Malaise, E.P.; Fertil, B.; Chavaudra, N.; Guichard, M. Distribution of Radiation Sensitivities for Human Tumor Cells of Specific Histological Types: Comparison of In Vivo to In Vitro Data. J. Radrorron Oncol. Biol. Phys 1986, 12, 617–624. [Google Scholar] [CrossRef]
  44. Nash, J.C.; Dalrymple, G.V.; Moss, A.J.; Baker, M.L. Initial Studies with a Line of Radioresistant Rat Tumor Cells. Radiat. Res. 1974, 60, 280–291. [Google Scholar] [CrossRef]
  45. Hall, E.J.; Giaccia, A.J. Radiobiology for the Radiologist, 6th Edition. Int. J. Radiat. Oncol. 2006, 66, 627. [Google Scholar] [CrossRef]
  46. Krock, B.L.; Skuli, N.; Simon, M.C. Hypoxia-Induced Angiogenesis: Good and Evil. Genes Cancer 2011, 2, 1117–1133. [Google Scholar] [CrossRef] [Green Version]
  47. Schulz, A.; Meyer, F.; Dubrovska, A.; Borgmann, K. Cancer Stem Cells and Radioresistance: DNA Repair and Beyond. Cancers 2019, 11, 862. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Brown, J.M.; Carlson, D.J.; Brenner, D.J. The Tumor Radiobiology of SRS and SBRT: Are More than the 5 Rs Involved? Int. J. Radiat. Oncol. Biol. Phys. 2014, 88, 254–262. [Google Scholar] [CrossRef] [Green Version]
  49. Branzei, D.; Foiani, M. Regulation of DNA Repair throughout the Cell Cycle. Nat. Rev. Mol. Cell Biol. 2008, 9, 297–308. [Google Scholar] [CrossRef] [PubMed]
  50. Liu, Y.; Chen, X.; Hu, Q.; Huang, J.; Ling, C. Resistance to Radiotherapy in Lung Cancer. Int. J. Clin. Exp. Med. 2018, 11, 7628–7642. [Google Scholar]
  51. Withers, H.R. The Four R’s of Radiotherapy; Academic Press, Inc.: Cambridge, MA, USA, 1975; Volume 5. [Google Scholar]
  52. Han, W.; Yu, K.N. Response of Cells to Ionizing Radiation. Adv. Biomed. Sci. Eng. 2009, 59, 204–262. [Google Scholar] [CrossRef] [Green Version]
  53. Naderi, S.; Hunton, I.C.; Wang, J.Y.J. Radiation Dose-Dependent Maintenance of G(2) Arrest Requires Retinoblastoma Protein. Cell Cycle 2002, 1, 192–199. [Google Scholar] [CrossRef] [Green Version]
  54. Brown, J.M. The Hypoxic Cell: A Target for Selective Cancer Therapy—Eighteenth Bruce F. Cain Memorial Award Lecture. Cancer Res. 1999, 59, 5863–5870. [Google Scholar]
  55. Nakamura, T.; Matsumine, A.; Matsubara, T.; Asanuma, K.; Uchida, A.; Sudo, A. Clinical Impact of the Tumor Volume Doubling Time on Sarcoma Patients with Lung Metastases. Clin. Exp. Metastasis 2011, 28, 819–825. [Google Scholar] [CrossRef]
  56. Nordsmark, M.; Bentzen, S.M.; Rudat, V.; Brizel, D.; Lartigau, E.; Stadler, P.; Becker, A.; Adam, M.; Molls, M.; Dunst, J.; et al. Prognostic Value of Tumor Oxygenation in 397 Head and Neck Tumors after Primary Radiation Therapy. An International Multi-Center Study. Radiother. Oncol. 2005, 77, 18–24. [Google Scholar] [CrossRef]
  57. Lewin, T.D.; Maini, P.K.; Moros, E.G.; Enderling, H.; Byrne, H.M. The Evolution of Tumour Composition During Fractionated Radiotherapy: Implications for Outcome. Bull. Math. Biol. 2018, 80, 1207–1235. [Google Scholar] [CrossRef] [Green Version]
  58. Thomlinson, R.H.; Gray, L.H. The Histological Structure of Some Human Lung Cancers and the Possible Implications for Radiotherapy. Br. J. Cancer 1955, 9, 539–549. [Google Scholar] [CrossRef] [Green Version]
  59. Crokart, N.; Jordan, B.F.; Baudelet, C.; Ansiaux, R.; Sonveaux, P.; Grégoire, V.; Beghein, N.; DeWever, J.; Bouzin, C.; Feron, O.; et al. Early Reoxygenation in Tumors after Irradiation: Determining Factors and Consequences for Radiotherapy Regimens Using Daily Multiple Fractions. Int. J. Radiat. Oncol. Biol. Phys. 2005, 63, 901–910. [Google Scholar] [CrossRef]
  60. Espinoza, I.; Peschke, P.; Karger, C.P. A Voxel-Based Multiscale Model to Simulate the Radiation Response of Hypoxic Tumors. Med. Phys. 2015, 42, 90–102. [Google Scholar] [CrossRef]
  61. Steel, G.G.; Mcmillan, T.J.; Peacock, J.H. The 5rs of Radiobiology. Int. J. Radiat. Biol. 1989, 56, 1045–1048. [Google Scholar] [CrossRef] [Green Version]
  62. Peeters, S.G.J.A.; Zegers, C.M.L.; Biemans, R.; Lieuwes, N.G.; Van Stiphout, R.G.P.M.; Yaromina, A.; Sun, J.D.; Hart, C.P.; Windhorst, A.D.; Van Elmpt, W.; et al. TH-302 in Combination with Radiotherapy Enhances the Therapeutic Outcome and Is Associated with Pretreatment [18F]HX4 Hypoxia PET Imaging. Clin. Cancer Res. 2015, 21, 2984–2992. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Sanduleanu, S.; van der Wiel, A.M.A.; Lieverse, R.I.Y.; Marcus, D.; Ibrahim, A.; Primakov, S.; Wu, G.; Theys, J.; Yaromina, A.; Dubois, L.J.; et al. Hypoxia PET Imaging with [18F]-HX4—A Promising Next-Generation Tracer. Cancers 2020, 12, 1322. [Google Scholar] [CrossRef] [PubMed]
  64. De Bruycker, S.; Vangestel, C.; Van Den Wyngaert, T.; Pauwels, P.; Wyffels, L.; Staelens, S.; Stroobants, S. 18F-Flortanidazole Hypoxia PET Holds Promise as a Prognostic and Predictive Imaging Biomarker in a Lung Cancer Xenograft Model Treated with Metformin and Radiotherapy. J. Nucl. Med. 2019, 60, 34–40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Lonati, L.; Barbieri, S.; Guardamagna, I.; Ottolenghi, A.; Baiocco, G. Radiation-Induced Cell Cycle Perturbations: A Computational Tool Validated with Flow-Cytometry Data. Sci. Rep. 2021, 11, 925. [Google Scholar] [CrossRef]
  66. Suwinski, R.; Wzietek, I.; Tarnawski, R.; Namysl-Kaletka, A.; Kryj, M.; Chmielarz, A.; Wydmanski, J. Moderately Low Alpha/Beta Ratio for Rectal Cancer May Best Explain the Outcome of Three Fractionation Schedules of Preoperative Radiotherapy. Int. J. Radiat. Oncol. Biol. Phys. 2007, 69, 793–799. [Google Scholar] [CrossRef] [PubMed]
  67. Dirscherl, T.; Rickhey, M.; Bogner, L. Feasibility of TCP-Based Dose Painting by Numbers Applied to a Prostate Case with 18F-Choline PET Imaging. Z. Med. Phys. 2012, 22, 48–57. [Google Scholar] [CrossRef] [PubMed]
  68. Brooks, S.; Glynne-Jones, R.; Novell, R.; Harrison, M.; Brown, K.; Makris, A. Short Course Continuous, Hyperfractionated, Accelerated Radiation Therapy (CHART) as Preoperative Treatment for Rectal Cancer. Acta Oncol. 2006, 45, 1079–1085. [Google Scholar] [CrossRef] [PubMed]
  69. Boustani, J.; Grapin, M.; Laurent, P.A.; Apetoh, L.; Mirjolet, C. The 6th R of Radiobiology: Reactivation of Anti-Tumor Immune Response. Cancers 2019, 11, 860. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Simulated oxygen pressure (pO2) histograms for different vascular fraction values vf (upper: <1%, <3%, <4%; lower: <5%, <8%, ≥8%).
Figure 1. Simulated oxygen pressure (pO2) histograms for different vascular fraction values vf (upper: <1%, <3%, <4%; lower: <5%, <8%, ≥8%).
Jimaging 09 00124 g001
Figure 2. Different fractionation schedules. Vertical arrows indicate the days when the dose was performed. Four protocols were studied: Standard Protocol (SP), DAHANCA Protocol (DP), Personal Protocol (PP) and CHART Protocol (CP). In the case of CP, three doses were delivered per day (for other details, see the text).
Figure 2. Different fractionation schedules. Vertical arrows indicate the days when the dose was performed. Four protocols were studied: Standard Protocol (SP), DAHANCA Protocol (DP), Personal Protocol (PP) and CHART Protocol (CP). In the case of CP, three doses were delivered per day (for other details, see the text).
Jimaging 09 00124 g002
Figure 3. Evolution of the tumour cell number during the treatment in the case of mouse 1.
Figure 3. Evolution of the tumour cell number during the treatment in the case of mouse 1.
Jimaging 09 00124 g003
Figure 4. Comparison of simulated and real images for two mice under carbogen (95% oxygen, 5% CO2). Input images were acquired before the beginning of radiotherapy. Simulated and real images were obtained after a single irradiation (for other details, see the text).
Figure 4. Comparison of simulated and real images for two mice under carbogen (95% oxygen, 5% CO2). Input images were acquired before the beginning of radiotherapy. Simulated and real images were obtained after a single irradiation (for other details, see the text).
Jimaging 09 00124 g004
Figure 5. Comparison of simulated and real images for two mice under breathing 7% conditions. Input images were acquired before the beginning of radiotherapy. Simulated and real images were obtained after a single irradiation (for other details see text).
Figure 5. Comparison of simulated and real images for two mice under breathing 7% conditions. Input images were acquired before the beginning of radiotherapy. Simulated and real images were obtained after a single irradiation (for other details see text).
Jimaging 09 00124 g005
Figure 6. The evolution of normal cell numbers during the treatment. Comparison of the evolution using the LQ Model (blue curve) and without it (red curve) (for other details, see the text).
Figure 6. The evolution of normal cell numbers during the treatment. Comparison of the evolution using the LQ Model (blue curve) and without it (red curve) (for other details, see the text).
Jimaging 09 00124 g006
Figure 7. Tumour response and cell distribution in the cell cycle. The evolution of the number of tumour cells took into account the distribution in the cell cycle (blue curve) and without (orange curve).
Figure 7. Tumour response and cell distribution in the cell cycle. The evolution of the number of tumour cells took into account the distribution in the cell cycle (blue curve) and without (orange curve).
Jimaging 09 00124 g007
Figure 8. Tumour response in mice in the presence or absence of variations in radiosensitivity during the cell cycle. (a): Weighting coefficients varied, and in (b), coefficients were constant (for other details, see the text).
Figure 8. Tumour response in mice in the presence or absence of variations in radiosensitivity during the cell cycle. (a): Weighting coefficients varied, and in (b), coefficients were constant (for other details, see the text).
Jimaging 09 00124 g008
Figure 9. The interactive 3D surface plots displayed a three-dimensional graph of the intensities of voxels. (a): Parameters of radiosensitivity were varied with every irradiation. Of note, for each irradiation per day, a new coefficient was used, and in (b), parameters were fixed (for other details, see the text).
Figure 9. The interactive 3D surface plots displayed a three-dimensional graph of the intensities of voxels. (a): Parameters of radiosensitivity were varied with every irradiation. Of note, for each irradiation per day, a new coefficient was used, and in (b), parameters were fixed (for other details, see the text).
Jimaging 09 00124 g009
Figure 10. Percentages of cells in each cell cycle phase (before and after radiotherapy).
Figure 10. Percentages of cells in each cell cycle phase (before and after radiotherapy).
Jimaging 09 00124 g010
Figure 11. Evolution of percentage of cells in each cell cycle phase over time. The cells were irradiated at t = 20 h.
Figure 11. Evolution of percentage of cells in each cell cycle phase over time. The cells were irradiated at t = 20 h.
Jimaging 09 00124 g011
Figure 12. Evolution of number of tumour cells with or without tumour repopulation over time.
Figure 12. Evolution of number of tumour cells with or without tumour repopulation over time.
Jimaging 09 00124 g012
Figure 13. The response of tumour cells according to tumour oxygenation.
Figure 13. The response of tumour cells according to tumour oxygenation.
Jimaging 09 00124 g013
Figure 14. Evolution of tumour cell numbers over time according to the α/β ratio.
Figure 14. Evolution of tumour cell numbers over time according to the α/β ratio.
Jimaging 09 00124 g014
Figure 15. Tumour response to radiotherapy for each fractionation protocol using simulated images: (a) Standard Protocol: SP, (b) DAHANCA Protocol: DP, (c) Personal Protocol: PP and (d) CHART Protocol: CP.
Figure 15. Tumour response to radiotherapy for each fractionation protocol using simulated images: (a) Standard Protocol: SP, (b) DAHANCA Protocol: DP, (c) Personal Protocol: PP and (d) CHART Protocol: CP.
Jimaging 09 00124 g015
Figure 16. Tumour response to radiotherapy. (a) Evolution of tumour cell numbers over time, and (b) the TCP curve for the simulated tumour. The TCD50 values (in Gy) of the TCP curves are 17.7, 17.9, 20.1 and 23.4 for PP, CP, SP and DP, respectively.
Figure 16. Tumour response to radiotherapy. (a) Evolution of tumour cell numbers over time, and (b) the TCP curve for the simulated tumour. The TCD50 values (in Gy) of the TCP curves are 17.7, 17.9, 20.1 and 23.4 for PP, CP, SP and DP, respectively.
Jimaging 09 00124 g016
Figure 17. The evolution of tumour cell numbers over time for a rhabdomyosarcoma tumour based on each protocol: Standard Protocol (SP), DAHANCA Protocol (DP), Personal Protocol (PP) and CHART Protocol (CP).
Figure 17. The evolution of tumour cell numbers over time for a rhabdomyosarcoma tumour based on each protocol: Standard Protocol (SP), DAHANCA Protocol (DP), Personal Protocol (PP) and CHART Protocol (CP).
Jimaging 09 00124 g017
Figure 18. TCP curves for each splitting protocol. The values of TCD50 are 13.23, 13.24, 13.72 and 18.18, respectively for CP, PP, DP and SP.
Figure 18. TCP curves for each splitting protocol. The values of TCD50 are 13.23, 13.24, 13.72 and 18.18, respectively for CP, PP, DP and SP.
Jimaging 09 00124 g018
Figure 19. Tumour images corresponding to the tumour on the seventh day after the start of treatment for each protocol.
Figure 19. Tumour images corresponding to the tumour on the seventh day after the start of treatment for each protocol.
Jimaging 09 00124 g019
Table 2. Evolution of phases (without treatment).
Table 2. Evolution of phases (without treatment).
Time (h)Percentage of Cells in:
G1SG2M
068.4420.048.772.75
2468.5816.4412.222.75
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hami, R.; Apeke, S.; Redou, P.; Gaubert, L.; Dubois, L.J.; Lambin, P.; Visvikis, D.; Boussion, N. Predicting the Tumour Response to Radiation by Modelling the Five Rs of Radiotherapy Using PET Images. J. Imaging 2023, 9, 124. https://doi.org/10.3390/jimaging9060124

AMA Style

Hami R, Apeke S, Redou P, Gaubert L, Dubois LJ, Lambin P, Visvikis D, Boussion N. Predicting the Tumour Response to Radiation by Modelling the Five Rs of Radiotherapy Using PET Images. Journal of Imaging. 2023; 9(6):124. https://doi.org/10.3390/jimaging9060124

Chicago/Turabian Style

Hami, Rihab, Sena Apeke, Pascal Redou, Laurent Gaubert, Ludwig J. Dubois, Philippe Lambin, Dimitris Visvikis, and Nicolas Boussion. 2023. "Predicting the Tumour Response to Radiation by Modelling the Five Rs of Radiotherapy Using PET Images" Journal of Imaging 9, no. 6: 124. https://doi.org/10.3390/jimaging9060124

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