Introduction
Treatment techniques and dose delivery have improved in radiotherapy. Intensity-modulated radiotherapy (IMRT) employs variable intensities across multiple beams. This yields highly conformal dose distributions.
IMRT generally requires multiple beams, which increases the treatment delivery time [1, 2]. Recently, volumetric modulated arc therapy (VMAT) has enabled treatment using one or two arcs [3]. The dose delivery time can be reduced compared with the IMRT technique. Stereotactic radiosurgery (SRS) and stereotactic body radiotherapy (SBRT) involve a large dose per fraction. This requires a longer dose delivery time than that for conventional radiotherapy. Recently, SBRT was combined with a flattening filter-free (FFF) beam, which uses a non-uniform beam and can reduce the treatment delivery time [4, 5].
Dose-calculation algorithms have advanced with the development of computer science. First, a factor-based algorithm was developed that resolves surface curvature and tissue heterogeneity based on effective spatial dose measures such as the path lengths in water, field boundaries, and depth [7]. Subsequently, a model-based algorithm that predicts patient dose distributions from a dose kernel and primary particle fluence was developed [8–10]. Recently, the Monte Carlo algorithm was introduced in clinical settings. The Monte Carlo algorithm addresses the stochasticity of radiation interaction with matter [11–14]. Furuta et al. developed novel medical applications for the particle and heavy-ion transport code system (PHITS) MC package [15]. PHITS can calculate the linear energy distribution to estimate the biological effectiveness. Earlier, we investigated the biological effectiveness of Lipiodol using flattening filter (FF) and FFF beams [16]. However, the biological effectiveness of the dose delivery time has not been investigated.
A prolonged delivery time affects radiobiological damage. Elkind et al. introduced sub-lethal damage repair (SLDR). Herein, cell death tends to decrease with a longer dose delivery time [6]. Nakano et al. reported a difference in dose delivery time between FFF and flattened filter (FF) beams [17]. Moreover, we proposed a dose compensation model for biological effectiveness based on the interruption time [18]. These studies simulated radiobiological effectiveness using a microkinetic model (MKM). However, these did not investigate the impact of tumour hypoxia or a prolonged delivery time. In clinical radiotherapy, the radiosensitivity of tumour cells decreases in hypoxic regions. This plays an important role in the progression of malignancy. Hall et al. reported that at < 20%, pO2 induced redio-resistance [19]. Hypoxia is an important factor limiting tumour prognosis. Tinganelli et al. demonstrated that the cell survival curve changes with the period during which the cells are exposed to hypoxic conditions [20]. Simulation of radiobiological effects including hypoxia and prolonged delivery time could contribute to addressing this limitation.
The current study improved the microdosimetric model by considering the oxygen dose enhancement in hypoxic tumours. Moreover, it has been used to evaluate the biological dose by varying the dose-delivery time in hypoxic tumours.
Materials and methods
Physical dose and lineal energy distribution in PHITS
To calculate the biological dose, the parameters of linear energy and physical dose were obtained from the Monte Carlo calculation code Particle and Heavy Ion Transport Code System (PHITS). The TrueBeam linear accelerator (Varian Medical Systems, Palo Alto, United States) that used phase-space files provided above the secondary jaw with a 6 MV x-ray beam was modelled in PHITS [18]. The dose was calculated in a virtual water phantom (20 × 20 × 20 cm) with a grid size of 2 mm and photon history of 4.0 × 109. The cut-off energies of photon and electron were set to 0.01 MeV and 0.7 MeV, respectively. The validation results showed that the Monte Carlo calculation and the measurement in a 0.04 cm3–volume CC04 (IBA Dosimetry, TN, United States) chamber for a 10 × 10 cm field at an SSD of 100 cm agreed within 1.0% [16]. The linear energy distribution was calculated using the T-SED function in PHITS [22].
Survival fraction in the MKM
The oxygen enhancement ratio (OER) is calculated as the dose in hypoxia divided by that in the specific oxygen concentration to achieve an equal survival level. The OER was estimated using the hypoxia reduction factor (HRF). It is the ratio of the doses for a specific iso-effect under a given oxygenation condition to that under the condition at 21% pO2. The MKM fit of clonogenic survival at a dose corrected for oxygen distribution was applied to predict the radiobiological effects in hypoxic tumours. The equation of HRF is used in the Howard–Flanders mathematical formalism as follows [22]:
(1)
where m is the maximum HRF, O2[c] is the amount of oxygen for any cell of the tumour, and K is the oxygen partial pressure at which the HRF is half the maximum value. m and K were set as 2.7 and 0.002, respectively, and were fitted to the experimental data obtained by Paul-Gilloteaux et al. [23]. To consider the oxygen effect for the response on the tumour after radiotherapy, the HRF was applied to calculate the dose (D’) corrected at 0.01% pO2 (hypoxia causing tumour death), 0.4% (radiobiological hypoxia), 1.0% pO2 (pathological hypoxia), 2.0% pO2 (physiological hypoxia), and 5.0% pO2 (physoxia) [24].
(2)
where D denotes the physical dose. The survival fraction was calculated using the MKM method proposed by Hawkins et al [25]. The surviving fraction of cells after irradiation, as proposed in a previous study, is expressed as follows [18]:
(3)
where α0 and β0 are the proportional factors to D [Gy−1] and D2 [Gy−2], respectively. The dose-mean lineal energy, yD, was set as 2.32 keV/µm [18]. The radius of the domain, rd, was 0.23 μm, and the density of the domain, ρ, was set as 1.0 g/cm3. was derived from a past study [26]:
(4)
where DR is the dose rate and T is the delivery time during irradiation. This is calculated as follows:
(5)
(a + c) is defined as the potentially lethal lesion repair rate obtained from Matsuya et al. [27]. Eq. (3) can be converted using by considering the following hypoxic reduction factor:
(6)
The biological dose (Dbio) proposed by Inaniwa et al. [26] using Eq. 4 was computed as
(7)
α0 and β0 were estimated from the experiment survival data of CHO-K1 cells after irradiation with X-rays under various oxygen conditions of 0%, 0.5%, and 21% pO2 [28, 29]. The CHO-K1 cells were cultured in Ham’s F12 medium supplemented with 1% penicillin/streptomycin (Biochrom AG, Berlin, Germany) and 10% foetal calf serum (FCS). Cells were maintained in a humidified atmosphere containing 5% CO2 at 37°C. These cells were subcultured at a density of 5 × 104 cells in 25 cm2 culture flasks with a 5 ml culture medium two times a week.
Biological dosimetric impact of dose-delivery time for hypoxic tumour
This study evaluated the biological dosimetric impact of the dose-delivery time on hypoxia using the following two metrics: The ratio of the biological dose with 1–24 Gy/min to that with 0.5 Gy/min for a D of 0–20 Gy at 0.01–5.0% pO2 was defined as RDR. The ratio of the biological dose at 0.01–2.0% pO2 to that at 5.0% pO2 was defined as Roxy.
Results
Validation of the MKM in survival fraction with a different oxygen distribution
The model parameters of α0 and β0 were determined by fitting the model to the experiment survival data [28, 29]. Figure 1 shows the fitting curve of the simulated and experimental survival data for the oxygen concentrations of 0%, 0.5%, and 20% pO2. The error bars indicate the standard deviations of the experimental survival data. The fit was in good agreement with the standard deviation of the experimental survival data. The fitted parameters of α0 and β0 were 0.22 Gy−1 and 0.195 Gy−2, respectively.
Biological dose difference with different oxygen concentrations
Figure 2 shows the Dbio for the physical dose owing to the dose rate of 0.5–24 Gy/min at 0.01–5.0% pO2. Dbio increased with an increase in the physical dose. The difference in Dbio of the dose rate was larger with a higher physical dose.
Figure 3 illustrates RDR. It is defined as the ratio of the Dbio at 1–24 Gy/min to that at 0.5 Gy/min. RDR was investigated at 0.01–5.0% pO2. It increased with a higher physical dose. The maximum RDR was 1.09 at 1 Gy/min, 1.12 at 12 Gy/min, and 1.13 at 24 Gy/min, respectively. Although the high dose rate X-ray beam had a larger RDR for each physical dose, the difference in RDR at 12 Gy/min and 24 Gy/min was within 1% pO2.
Biological dose difference with a different dose rate
Figure 4 shows the Dbio for the physical dose at 0.01–5.0% pO2 with a dose rate of 0.5–24 Gy/min. Dbio increased with higher concentrations of oxygen.
Figure 5 shows the ROxy at 0.01–5.0% pO2 for a D of 0–20 Gy with a dose rate of 0.5–24 Gy/min. A higher concentration of oxygen resulted in a higher ROxy. The ROxy for a physical dose of 2–20 Gy was 0.4–0.42, 0.75–0.76, 0.88–0.89, and 0.95–0.96 at 0.01%, 0.4%, 1.0%, and 2.0% pO2 of oxygen concentration, respectively. The difference in ROxy in the range of 2–20 Gy was within 1.0% for all the oxygen concentrations. The difference in ROxy for 0.5–24 Gy/min was also within 1.0%.
Discussion
The dose delivery time depends on the dose rate. That is, a higher dose rate can reduce the dose delivery time. This study demonstrated that a higher dose rate increased the biological dose. Brehwens et al. reported that DNA repair occurred under irradiation with a photon beam of less than 1 Gy/min [30]. In particular, FFF beams at > 12 Gy/min exhibited a larger biological effect. Nakano et al. reported that a short dose-delivery time increased the relative biological effectiveness [17]. A higher dose rate reduces the SLDR. This may cause an increase in the treatment outcome [31]. Recently, FLASH irradiation with an ultrahigh dose-rate of irradiation (> 40 Gy/s) was developed. It enables the sparing of normal tissue while retaining tumour control. Further studies would be performed to investigate the dose-rate effect in normal tissues [31].
This study investigated the effects of oxygen and the dose rates on tumours with an improved MKM. The biological effectiveness increased with a high concentration of oxygen and higher photon-beam dose rate. Hara et al. investigated the radiobiological effects on hypoxic and oxic cells using high- and low-dose-rate beams [32]. They showed that a higher dose rate decreased the surviving fraction. Additionally, the surviving fraction of oxic cells was lower than that of hypoxic cells at an equal photon beam dose rate. These results support the results of the simulations. The neutralisation of radiotherapy-induced reactive oxygen species decreases at higher oxygen concentrations because oxygen helps neutralise hydrated electrons. This causes an increase in the DNA damage [33, 34]. The current study showed that a beam with a higher dose rate, such as an FFF beam, shortens the dose delivery time and increases the biological effectiveness. In particular, it increased exponentially with a higher dose rate and physical dose. Moreover, the oxygen concentration helped increase the biological effectiveness, although the difference in biological effectiveness between the physical doses at each oxygen concentration was marginal.
A limitation of the current study is that other biological effects such as angiogenesis and the cell cycle were not incorporated. Our simulation incorporated a Monte Carlo calculation that can provide an accurate dose and better estimates of the biological dose or tumour control ratio for radiotherapy treatment planning and trial design by incorporating accurate radiobiological models.
Conclusions
Our proposed model could estimate the cell killing and biological dose under hypoxia in a clinical and realistic patient. A shorter dose-delivery time with a higher oxygen distribution increased the radiobiological effect. It was more effective at higher doses per fraction than at lower doses.
Conflict of interest
None declared.
Funding
None declared.