Modelling of cell killing due to sparsely ionizing radiation in normoxic and hypoxic conditions and an extension to high LET radiation

Modeling of ionizing radiation-induced chromosome aberration and tumor prevalence based on two classes of DNA double-strand breaks clustering in chromatin domains

There has been some controversy over the use of radiobiological models when modeling the dose-response curves of ionizing radiation (IR)-induced chromosome aberration and tumor prevalence, as those curves usually show obvious non-targeted effects (NTEs) at low doses of high linear energy transfer (LET) radiation. The lack of understanding the contribution of NTEs to IR-induced carcinogenesis can lead to distinct deviations of relative biological effectiveness (RBE) estimations of carcinogenic potential, which are widely used in radiation risk assessment and radiation protection. In this work, based on the initial pattern of two classes of IR-induced DNA double-strand breaks (DSBs) clustering in chromatin domains and the subsequent incorrect repair processes, we proposed a novel radiobiological model to describe the dose-response curves of two carcinogenic-related endpoints within the same theoretical framework. The representative experimental data was used to verify the consistency and validity of the present model. The fitting results indicated that, compared with targeted effect (TE) and NTE models, the current model has better fitting ability when dealing with the experimental data of chromosome aberration and tumor prevalence induced by multiple types of IR with different LETs. Notably, the present model without introducing an NTE term was adequate to describe the dose-response curves of IR-induced chromosome aberration and tumor prevalence with NTEs in low-dose regions. Based on the fitting parameters, the LET-dependent RBE values were calculated for three given low doses. Our results showed that the RBE values predicted by the current model gradually decrease with the increase of doses for the endpoints of chromosome aberration and tumor prevalence. In addition, the calculated RBE was also compared with those evaluated from other models. These analyses show that the proposed model can be used as an alternative tool to well describe dose-response curves of multiple carcinogenic-related endpoints and effectively estimate RBE in low-dose regions.

The Mayo Clinic Florida Microdosimetric Kinetic Model of Clonogenic Survival: Application to Various Repair-Competent Rodent and Human Cell Lines

The relative biological effectiveness (RBE) calculations used during the planning of ion therapy treatments are generally based on the microdosimetric kinetic model (MKM) and the local effect model (LEM). The Mayo Clinic Florida MKM (MCF MKM) was recently developed to overcome the limitations of previous MKMs in reproducing the biological data and to eliminate the need for ion-exposed in vitro data as input for the model calculations. Since we are considering to implement the MCF MKM in clinic, this article presents (a) an extensive benchmark of the MCF MKM predictions against corresponding in vitro clonogenic survival data for 4 rodent and 10 cell lines exposed to ions from 1H to 238U, and (b) a systematic comparison with published results of the latest version of the LEM (LEM IV). Additionally, we introduce a novel approach to derive an approximate value of the MCF MKM model parameters by knowing only the animal species and the mean number of chromosomes. The overall good agreement between MCF MKM predictions and in vitro data suggests the MCF MKM can be reliably used for the RBE calculations. In most cases, a reasonable agreement was found between the MCF MKM and the LEM IV.

The DNA Double-Strand Break Repair in Glioma: Molecular Players and Therapeutic Strategies

Gliomas are the most frequent type of tumor in the central nervous system, which exhibit properties that make their treatment difficult, such as cellular infiltration, heterogeneity, and the presence of stem-like cells responsible for tumor recurrence. The response of this type of tumor to chemoradiotherapy is poor, possibly due to a higher repair activity of the genetic material, among other causes. The DNA double-strand breaks are an important type of lesion to the genetic material, which have the potential to trigger processes of cell death or cause gene aberrations that could promote tumorigenesis. This review describes how the different cellular elements regulate the formation of DNA double-strand breaks and their repair in gliomas, discussing the therapeutic potential of the induction of this type of lesion and the suppression of its repair as a control mechanism of brain tumorigenesis.

Impact of DNA Repair Kinetics and Dose Rate on RBE Predictions in the UNIVERSE

Accurate knowledge of the relative biological effectiveness (RBE) and its dependencies is crucial to support modern ion beam therapy and its further development. However, the influence of different dose rates of the reference radiation and ion beam are rarely considered. The ion beam RBE-model within our "UNIfied and VERSatile bio response Engine" (UNIVERSE) is extended by including DNA damage repair kinetics to investigate the impact of dose-rate effects on the predicted RBE. It was found that dose-rate effects increase with dose and biological effects saturate at high dose-rates, which is consistent with data- and model-based studies in the literature. In a comparison with RBE measurements from a high dose in-vivo study, the predictions of the presented modification were found to be improved in comparison to the previous version of UNIVERSE and existing clinical approaches that disregard dose-rate effects. Consequently, DNA repair kinetics and the different dose rates applied by the reference and ion beams might need to be considered in biophysical models to accurately predict the RBE. Additionally, this study marks an important step in the further development of UNIVERSE, extending its capabilities in giving theoretical guidance to support progress in ion beam therapy.

The Impact of Sub-Millisecond Damage Fixation Kinetics on the In Vitro Sparing Effect at Ultra-High Dose Rate in UNIVERSE

The impact of the exact temporal pulse structure on the potential cell and tissue sparing of ultra-high dose-rate irradiation applied in FLASH studies has gained increasing attention. A previous version of our biophysical mechanistic model (UNIVERSE: UNIfied and VERSatile bio response Engine), based on the oxygen depletion hypothesis, has been extended in this work by considering oxygen-dependent damage fixation dynamics on the sub-milliseconds scale and introducing an explicit implementation of the temporal pulse structure. The model successfully reproduces in vitro experimental data on the fast kinetics of the oxygen effect in irradiated mammalian cells. The implemented changes result in a reduction in the assumed amount of oxygen depletion. Furthermore, its increase towards conventional dose-rates is parameterized based on experimental data from the literature. A recalculation of previous benchmarks shows that the model retains its predictive power, while the assumed amount of depleted oxygen approaches measured values. The updated UNIVERSE could be used to investigate the impact of different combinations of pulse structure parameters (e.g., dose per pulse, pulse frequency, number of pulses, etc.), thereby aiding the optimization of potential clinical application and the development of suitable accelerators.

Combined DNA Damage Repair Interference and Ion Beam Therapy: Development, Benchmark, and Clinical Implications of a Mechanistic Biological Model

PURPOSE

Our purpose was to develop a mechanistic model that describes and predicts radiation response after combined DNA damage repair interference (DDRi) and particle radiation therapy.

METHODS AND MATERIALS

The heterogeneous dose distributions of protons and ⁴He ions were implemented into the "UNIfied and VERSatile bio-response Engine" (UNIVERSE). Predictions for monoenergetic and mixed fields over clinically relevant dose and linear energy transfer range were compared with experimental in vitro survival data measured in this work as well as data available in the literature, including different cell lines and DDR interferences. Ultimately, UNIVERSE predictions were investigated in a patient plan.

RESULTS

UNIVERSE accurately predicts survival of cell lines with and without DDRi in clinical settings of ion beam therapy based only on 3 parameters derived from photon data. With increasing dose or linear energy transfer, the radiosensitizing effect of DDRi decreases, resulting in diminished relative biological effect of ion beam radiation for cells subjected to DDRi in comparison to cells that are not. Similar trends were observed in patient plan recalculations; however, this analysis also suggests that DDRi + particle radiation therapy may better preserve the therapeutic window in comparison to DDRi + photon radiation therapy.

CONCLUSIONS

The presented framework represents the first mechanistic model of combined DDRi and particle radiation therapy comprehensively benchmarked in clinically relevant scenarios and a step toward more personalized treatment. It reveals potential differences between DDRi + photon radiation therapy versus DDRi + particle radiation therapy, which have not been described so far. UNIVERSE could aid in appraising the clinical viability of combined administration of radiosensitizing drugs and charged particle therapy, as well as the identification of patients with known DDR deficiencies in the tumor who might benefit from therapy with light ions, freeing limited space at heavy ion therapy centers.

Spot-Scanning Hadron Arc (SHArc) Therapy: A Study With Light and Heavy Ions

PURPOSE

To evaluate the clinical potential of spot-scanning hadron arc (SHArc) therapy with a heavy-ion gantry.

METHODS AND MATERIALS

A series of in silico studies was conducted via treatment plan optimization in FRoG and the RayStation TPS to compare SHArc therapy against reference plans using conventional techniques with single, parallel-opposed, and 3-field configurations for 3 clinical particle beams (protons [p], helium [4He], and carbon [12C] ions). Tests were performed on water-equivalent cylindrical phantoms for simple targets and clinical-like scenarios with an organ-at-risk in proximity of the target. Effective dose and dose-averaged linear energy transfer (LETD) distributions for SHArc were evaluated against conventional planning techniques applying the modified microdosimetric kinetic model for considering bio-effect with (α/β)x = 2 Gy. A model for hypoxia-induced tumor radio-resistance was developed for particle therapy with dependence on oxygen concentration and particle species/energy (Zeff/β)2 to investigate the impact on effective dose.

RESULTS

SHArc plans exhibited similar target coverage with unique treatment attributes and distributions compared with conventional planning, with carbon ions demonstrating the greatest potential for tumor control and normal tissue sparing among the arc techniques. All SHArc plans exhibited a low-dose bath outside the target volume with a reduced maximum dose in normal tissues compared with single, parallel-opposed, and 3-field configuration plans. Moreover, favorable LETD distributions were made possible using the SHArc approach, with maximum LETD in the r = 5 mm tumor core (~8 keVμm-1, ~30 keVμm-1, and ~150 keVμm-1 for p, 4He, and 12C ions, respectively) and reductions of high-LET regions in normal tissues and organs-at-risk compared with static treatment beam delivery.

CONCLUSION

SHArc therapy offers potential treatment benefits such as increased normal tissue sparing. Without explicit consideration of oxygen concentration during treatment planning and optimization, SHArc-C may mitigate tumor hypoxia-induced loss of efficacy. Findings justify further development of robust SHArc treatment planning toward potential clinical translation.

Modeling Direct and Indirect Action on Cell Survival After Photon Irradiation under Normoxia and Hypoxia

The demand for personalized medicine in radiotherapy has been met by a surge of mechanistic models offering predictions of the biological effect of ionizing radiation under consideration of a growing number of parameters. We present an extension of our existing model of cell survival after photon irradiation to explicitly differentiate between the damage inflicted by the direct and indirect (radicals-mediated) action of ionizing radiation. Within our approach, we assume that the oxygenation status affects the indirect action. The effect of different concentrations of dimethyl sulfoxide (DMSO), an effective radical scavenger, has been simulated at different dose levels in normoxic and hypoxic conditions for various cell lines. Our model is found to accurately predict experimental data available in literature, validating the assumptions made in our approach. The presented extension adds further flexibility to our model and could act as basis for further developments of our model.

Generalized Multi-Hit Model of Radiation-Induced Cell Survival with a Closed-Form Solution: An Alternative Method for Determining Isoeffect Doses in Practical Radiotherapy

The standard linear-quadratic (LQ) model is currently the preferred model for describing the ionizing radiation-induced cell survival curves and tissue responses. And the LQ model is also widely used to calculate isoeffect doses for comparing different fractionated schemes in clinical radiotherapy. Despite its ubiquity, because the actual dose-response curve may appear linear at high doses in the semilogarithmic plot, the application of the LQ model is greatly challenged in the high-dose region, while the dose employed in stereotactic body radiotherapy (SBRT) is often in this area. Alternatively, the biophysical models of radiation-induced effects with a linear-quadratic-linear (LQL) characteristic can well fit the dose-survival curve of cells in vitro. However, most of these LQL models are phenomenological and have not fully considered the biophysical mechanism of radiation-induced damage and repair, and the fitting quality decreases in some high-dose ranges. In this work, to provide an alternative model to describe the cell survival curves in high-dose ranges and predict the biologically effective dose (BED) for SBRT, we propose a novel generalized multi-hit model with a closed-form solution by considering an upper bound on the number of lethal damages induced by radiation that can be repaired in a cell. This model has a clear biophysical basis and a simple expression, and also has the LQL characteristic under low- and high-dose approximate conditions. The experimental data fitting indicated that compared to the standard LQ model and our previously generalized target model, the current model can better fit the radiation-induced cell survival curves in the high-dose ranges (P < 0.05). The current model parameters and parameter ratios were determined from the fits in different kinds of cell lines irradiated with various dose rates and linear energy transfer (LET), which indicates that the model parameters significantly depend on the dose rate and LET. Based on the current model, we derived two equivalence formulae for the BED calculations in the low- and high-dose ranges, and then calculated the BED for the clinical data of SBRT from 17 selected studies. The correlation analysis showed that there were significant linear correlations between the BED at isocenter and planning target volume (PTV) edge calculated by this model and the LQ model (R > 0.86, P < 0.001). In conclusion, the generalized multi-hit model proposed in this work can be used as an alternative tool to handle in vitro radiation-induced cell survival curves in high-dose ranges, and calculate the in vivo BED for comparing the dose fractionation schemes in clinical radiotherapy.

Modeling the Effect of Hypoxia and DNA Repair Inhibition on Cell Survival After Photon Irradiation

Mechanistic approaches to modeling the effects of ionizing radiation on cells are on the rise, promising a better understanding of predictions and higher flexibility concerning conditions to be accounted for. In this work we modified and extended a previously published mechanistic model of cell survival after photon irradiation under hypoxia to account for radiosensitization caused by deficiency or inhibition of DNA damage repair enzymes. The model is shown to be capable of describing the survival data of cells with DNA damage repair deficiency, both under norm- and hypoxia. We find that our parameterization of radiosensitization is invariant under change of oxygen status, indicating that the relevant parameters for both mechanisms can be obtained independently and introduced freely to the model to predict their combined effect.

Effects of irradiation conditions on the radiation sensitivity of microorganisms in the presence of OH-radical scavengers

PURPOSE

The purpose of the paper was to investigate the protective effect of some scavengers of OH radicals (hydroxyl radicals) on the radiation sensitivity of bacteria (in some cases also yeast) under normoxic (in air) or hypoxic (suboxic) conditions and to compare the obtained results with those published earlier for the yeast, all in a wide interval of irradiation conditions. Another aim was to investigate a possible impact of the reaction order of the reaction between the scavengers and the OH radicals on the protection of the cells.

MATERIALS AND METHODS

In order to study the protective effect of OH scavengers we used various concentrations of methanol and potassium formate (in some cases also ethanol) in isotonic salt solutions. These solutions containing living bacteria (Escherichia coli) or yeast (Saccharomyces cerevisiae) were irradiated with ⁶⁰Co radiation using various doses and dose rates. Irradiation was performed in air, in some cases under the hypoxic conditions. The number of surviving cells was determined prior to and after irradiation in suspension with and without scavengers. The surviving fractions after irradiation with and without scavenger were evaluated.

RESULTS

The surviving fraction of bacterial cells increases linearly with increasing concentration of both scavengers. The fraction of surviving cells does not increase with increasing concentration of the scavengers under suboxic conditions. The protective effect Ϭ increases linearly with increasing scavenging efficiency and this dependence is much sharper under normoxic conditions than under suboxic ones. The specific protection k is much higher for the methanol than for the potassium formate.

CONCLUSIONS

The basic characteristics of the impact of scavengers of OH radicals on radiation sensitivity of both bacteria and yeast are the same in a wide interval of doses and dose rates. The specific protection effect is much higher under the normoxic conditions. This protective effect is inversely proportional to the rate constant of the reaction between the scavenger and the OH radicals. It seems to be obvious that the presence of oxygen during irradiation is a necessary condition for the protective action of the scavengers which may be partially controlled by some transport processes and may be connected with the radiation sensitivity of the cells. On the other hand, the change of the reaction order of the reaction of the scavenger with the OH radicals has turned out to be unimportant.

Impact of 18F-FET PET on Target Volume Definition and Tumor Progression of Recurrent High Grade Glioma Treated with Carbon-Ion Radiotherapy

High-precision radiotherapy (HPR) of recurrent high grade glioma (HGG) requires accurate spatial allocation of these infiltrative tumors. We investigated the impact of ¹⁸F-FET PET on tumor delineation and progression of recurrent HGG after HPR with carbon ions. T₁ contrast enhanced MRI and ¹⁸F-FET-PET scans of 26 HGG patients were fused with radiotherapy planning volumes. PET-positive (PET+) tumor volumes using different isocontours (I%) were systematically investigated and compared with MRI-derived gross tumor volumes (GTV). Standardized uptake ratios (SUR) were further correlated with GTV and tumor progression patterns. In grade IV glioma, SUR > 2.92 significantly correlated with poor median overall survival (6.5 vs 13.1 months, p = 0.00016). We found no reliable SUR cut-off criteria for definition of PET+ volumes. Overall conformity between PET and MRI-based contours was low, with maximum conformities between 0.42–0.51 at I40%. The maximum sensitivity and specificity for PET+ volumes outside of GTV predicting tumor progression were 0.16 (I40%) and 0.52 (I50%), respectively. In 75% of cases, FLAIR hyperintense area covered over 80% of PET+ volumes. ¹⁸F-FET-PET derived SUR has a prognostic impact in grade IV glioma. The value of substantial mismatches between MRI-based GTV and PET+ volumes to improve tumor delineation in radiotherapy awaits further validation in randomized prospective trials.

High resistance to X-rays and therapeutic carbon ions in glioblastoma cells bearing dysfunctional ATM associates with intrinsic chromosomal instability

PURPOSE

To investigate chromosomal instability and radiation response mechanisms in glioblastoma cells.

METHODS AND MATERIALS

We undertook a comparative analysis of two patient-derived glioblastoma cell lines. Their resistance to low and high linear energy transfer (LET) radiation was assessed using clonogenic survival assay and their intrinsic chromosome instability status using fluorescence in situ hybridization. DNA damage was analyzed by pulsed-field gel electrophoresis and by γ-H2AX foci quantification. Expression of DNA damage response proteins was assessed by immunoblot.

RESULTS

Increased radioresistance to X-rays as well as carbon ions was observed in glioblastoma cells exhibiting high levels of naturally occurring chromosomal instability and impaired Ataxia-telangiectasia mutated (ATM) signaling, as reflected by lack of phosphorylation of ATM, CHK2 and p53 after double-strand breaks induction.

CONCLUSION

Our results indicate the existence of highly radioresistant glioblastoma cells, characterized by dysfunctional ATM signaling and high levels of intrinsic chromosomal instability.