Physical characteristics at the turnover-points of relative biological effect (RBE) with linear energy transfer (LET)

Modelling tissue specific RBE for different radiation qualities based on a multiscale characterization of energy deposition

PURPOSE

We present the nanoCluE model, which uses nano- and microdosimetric quantities to model RBE for protons and carbon ions. Under the hypothesis that nano- and microdosimetric quantities correlates with the generation of complex DNA double strand breakes, we wish to investigate whether an improved accuracy in predicting LQ parameters may be achieved, compared to some of the published RBE models.

METHODS

The model is based on experimental LQ data for protons and carbon ions. We generated a database of track structure data for a number of proton and carbon ion kinetic energies with the Geant4-DNA Monte Carlo code. These data were used to obtain both a nanodosimetric quantity and a set of microdosimetric quantities. The latter were tested with different parameterizations versus experimental LQ-data to select the variable and parametrization that yielded the best fit.

RESULTS

For protons, the nanoCluE model yielded, for the ratio of the linear LQ term versus the test data, a root mean square error (RMSE) of 1.57 compared to 1.31 and 1.30 for two earlier other published proton models. For carbon ions the RMSE was 2.26 compared to 3.24 and 5.24 for earlier published carbon ion models.

CONCLUSION

These results demonstrate the feasibility of the nanoCluE RBE model for carbon ions and protons. The increased accuracy for carbon ions as compared to two other considered models warrants further investigation.

The influence of hypoxia on LET and RBE relationships with implications for ultra-high dose rates and FLASH modelling

Objective. To investigate relationships between linear energy transfer (LET), fluence rates, changes in radiosensitivity and the oxygen enhancement ratio (OER) in different ion beams and extend these concepts to ultra-high dose rate (UHDR) or FLASH effects. Approach. LET values providing maximum relative biological effect (RBE), designated as LETU, are found for neon, carbon and helium beams. Proton experiments show reduced RBEs with depth in scattered (divergent) beams, but not with scanned beams, suggesting that instantaneous fluence rates (related to track separation distances) can modify RBE, all other RBE-determining factors being equal. Micro-volumetric energy transfer per μm3 (mVET) is defined by LET × fluence. High fluence rates will increase mVET rates, with proportional shifts of LETU to lower values due to more rapid energy transfer. From the relationship between LETU and OER at conventional dose rates, OER reductions in UHDR/FLASH exposures can be estimated and biological effective dose analysis of experimental lung and skin reactions becomes feasible. Main results. The Furusawa et al data show that hypoxic LETU values exceed their oxic counterparts. OER reduces from around 3–1.25 at LETU, although the relative radiosensitivities of the oxic and hypoxic α parameters (the OER(α)) exceed those of the standard OER values. Increased fluence rates are predicted to reduce LETU and OER. Large FLASH single doses will minimise RBE increments due to the β parameter reducing by a factor of 0.5–0.25 consistent with oxygen depletion, causing radioresistance. Similar results will occur for photons. Tissue α/β ratios increase by around 10 in FLASH conditions, agreeing with derived ion-beam dose rate equations. Significance. Increasing dose rates enhance local energy deposition rate per unit volume, probably causing oxygen depletion and radioresistance in pre-existing hypoxic sites during UHDR/FLASH exposures. The modelled equations provide testable hypotheses for further dose rate investigations in photon, proton and ion beams.

Further development of spinal cord retreatment dose estimation: including radiotherapy with protons and light ions

PURPOSE

A graphical user interface (GUI) was developed to aid in the assessment of changes in the radiation tolerance of spinal cord/similar central nervous system tissues with time between two individual treatment courses.

METHODS

The GUI allows any combination of photons, protons (or ions) to be used as the initial, or retreatment, radiotherapy courses. Allowances for clinical circumstances, of reduced tolerance, can also be made. The radiobiological model was published previously and has been incorporated with additional checks and safety features, to be as safe to use as possible. The proton option includes use of a fixed RBE of 1.1 (set as the default), or a variable RBE, the latter depending on the proton linear energy transfer (LET) for organs at risk. This second LET-based approach can also be used for ions, by changing the LET parameters.

RESULTS

GUI screenshots are used to show the input and output parameters for different clinical situations used in worked examples. The results from the GUI are in agreement with manual calculations, but the results are now rapidly available without tedious and error-prone manual computations. The software outputs provide a maximum dose limit boundary, which should not be exceeded. Clinicians may also choose to further lower the number of treatment fractions, whilst using the same dose per fraction (or conversely a lower dose per fraction but with the same number of fractions) in order to achieve the intended clinical benefit as safely as possible.

CONCLUSIONS

The new GUI will allow scientific-based estimations of time related radiation tolerance changes in the spinal cord and similar central nervous tissues (optic chiasm, brainstem), which can be used to guide the choice of retreatment dose fractionation schedules, with either photons, protons or ions.

Update of the particle irradiation data ensemble (PIDE) for cell survival

The particle irradiation data ensemble (PIDE) is the largest database of cell survival data measured after exposure to ion beams and photon reference radiation. We report here on the updated version of the PIDE database and demonstrate how to investigate generic properties of radiation dose response using these sets of raw data. The database now contains information of over 1100 pairs of photon and ion dose response curves. It provides the originally published raw data of cell survival in addition to given linear quadratic (LQ) model parameters. If available, the raw data were used to derive LQ model parameters in the same way for all experiments. To demonstrate the extent of the database and the variability among experiments we focus on the dose response curves after ion and photon radiation separately in a first step. Furthermore, we discuss the capability and the limitations of the database for analyzing properties of low and high linear energy transfer (LET) radiation response based on multiple experiments. PIDE is freely available to the research community under www.gsi.de/bio-pide.

High LET-Like Radiation Tracks at the Distal Side of Accelerated Proton Bragg Peak

Proton therapy is a type of hadron radiotherapy used for treating solid tumors. Unlike heavy charged elements, proton radiation is considered to be low LET (Linear Energy Transfer) radiation, like X-rays. However, the clinical SOBP (Spread Out Bragg Peak) proton radiation is considered to be higher in relative biological effectiveness (RBE) than both X-ray and their own entrance region. The RBE is estimated to be 1.1–1.2, which can be attributed to the higher LET at the SOBP region than at the entrance region. In order to clarify the nature of higher LET near the Bragg peak of proton radiation and its potential cytotoxic effects, we utilized a horizontal irradiation system with CHO cells. Additionally, we examined DNA repair mutants, analyzed cytotoxicity with colony formation, and assessed DNA damage and its repair with γ-H2AX foci assay in a high-resolution microscopic scale analysis along with the Bragg peak. Besides confirming that the most cytotoxic effects occurred at the Bragg peak, extended cytotoxicity was observed a few millimeters after the Bragg peak. γ-H2AX foci numbers reached a maximum at the Bragg peak and reduced dramatically after the Bragg peak. However, in the post-Bragg peak region, particle track-like structures were sporadically observed. This region contains foci that are more difficult to repair. The peak and post-Bragg peak regions contain rare high LET-like radiation tracks and can cause cellular lethality. This may have caused unwanted side effects and complexities of outputs for the proton therapy treatment.

Fast neutron energy based modelling of biological effectiveness with implications for proton and ion beams

A practical neutron energy dependent RBE model has been developed, based on the relationship between a mono-energetic neutron energy and its likely recoil proton energy. Essentially, the linear energy transfer (LET) values of the most appropriate recoil proton energies are then used to modify the linear quadratic model radiosensitivities (α and β) from their reference LET radiation values to provide the RBE estimates. Experimental neutron studies published by Hall (including some mono-energetic beams ranging from 0.2 to 15 MeV), Broerse, Berry, and data from the Clatterbridge and Detroit clinical neutron beams, which all contain some data from a spectrum of neutron energies, are used to derive single effective neutron energies (NE_eff) for each spectral beam. These energies yield a recoil proton spectrum, but with an effective mean proton energy (being around 50% of NE_eff). The fractional increase in LET is given by the recoil proton LET divided by the proton (LET_U) value which provides the highest RBE. This ratio is then used to determine the change in the linear-quadratic model α and β parameters, from those of the reference radiation, to estimate the RBE. The predicted proton recoil RBE is then reasonably close to the experimental neutron RBE values found when taking into account the variation inherent in biological experiments. The work has some important consequences. The data of Hall et al (1975 Radiat. Res. 64 245-55) shows that the highest RBE values are found with neutron energies around 0.3-0.4 MeV, but this energy cannot possibly generate recoil proton energies which are higher, as necessary for a 0.68 MeV proton with a 30.5 keV μm^-1 LET_U (the LET value which provides the maximum obtainable RBE for a specified ion). For 0.4 MeV neutrons with proton recoil energies of around 0.2 MeV, the latter have a LET of around 62.88 keV μm^-1. This could have an impact on proton beam RBE modelling. However, this is compensated by finding that the maximum radiosensitivity for mono-energetic neutrons was around 1.7 times larger than previously suggested from experimental ion beam studies, probably due to the necessary spreading out of Bragg peaks for ion beam experimental purposes, sampling errors and particle range considerations. This semi-empirical model can be used with minimal computer support and could have applications in ionic beams and in radioprotection.

Assessment of RBE-Weighted Dose Models for Carbon Ion Therapy Toward Modernization of Clinical Practice at HIT: In Vitro, in Vivo, and in Patients

PURPOSE

Present-day treatment planning in carbon ion therapy is conducted with assumptions for a limited number of tissue types and models for effective dose. Here, we comprehensively assess relative biological effectiveness (RBE) in carbon ion therapy and associated models toward the modernization of current clinical practice in effective dose calculation.

METHODS

Using 2 human (A549, H460) and 2 mouse (B16, Renca) tumor cell lines, clonogenic cell survival assay was performed for examination of changes in RBE along the full range of clinical-like spread-out Bragg peak (SOBP) fields. Prediction power of the local effect model (LEM1 and LEM4) and the modified microdosimetric kinetic model (mMKM) was assessed. Experimentation and analysis were carried out in the frame of a multidimensional end point study for clinically relevant ranges of physical dose (D), dose-averaged linear energy transfer (LET_d), and base-line photon radio-sensitivity (α/β)_x. Additionally, predictions were compared against previously reported RBE measurements in vivo and surveyed in patient cases.

RESULTS

RBE model prediction performance varied among the investigated perspectives, with mMKM prediction exhibiting superior agreement with measurements both in vitro and in vivo across the 3 investigated end points. LEM1 and LEM4 performed their best in the highest LET conditions but yielded overestimations and underestimations in low/midrange LET conditions, respectively, as demonstrated by comparison with measurements. Additionally, the analysis of patient treatment plans revealed substantial variability across the investigated models (±20%-30% uncertainty), largely dependent on the selected model and absolute values for input tissue parameters α_x and β_x.

CONCLUSION

RBE dependencies in vitro, in vivo, and in silico were investigated with respect to various clinically relevant end points in the context of tumor-specific tissue radio-sensitivity assignment and accurate RBE modeling. Discovered model trends and performances advocate upgrading current treatment planning schemes in carbon ion therapy and call for verification via clinical outcome analysis with large patient cohorts.