A method to predict space radiation biological effectiveness for non-cancer effects following intense Solar Particle Events

In addition to the continuous exposure to cosmic rays, astronauts in space are occasionally exposed to Solar Particle Events (SPE), which involve less energetic particles but can deliver much higher doses. The latter can exceed several Gy in a few hours for the most intense SPEs, for which non-stochastic effects are thus a major concern. To identify adequate shielding conditions that would allow respecting the dose limits established by the various space agencies, the absorbed dose in the considered organ/tissue must be multiplied by the corresponding Relative Biological Effectiveness (RBE), which is a complex quantity depending on several factors including particle type and energy, considered biological effect, level of effect (and thus absorbed dose), etc. While in several studies only the particle-type dependence of RBE is taken into account, in this work we developed and applied a new approach where, thanks to an interface between the FLUKA Monte Carlo transport code and the BIANCA biophysical model, the RBE dependence on particle energy and absorbed dose was also considered. Furthermore, we included in the considered SPE spectra primary particles heavier than protons, which in many studies are neglected. This approach was then applied to the October 2003 SPE (the most intense SPE of solar cycle 23, also known as "Halloween event") and the January 2005 event, which was characterized by a lower fluence but a harder spectrum, i.e., with higher-energy particles. The calculation outcomes were then discussed and compared with the current dose limits established for skin and blood forming organs in case of 30-days missions. This work showed that the BIANCA model, if interfaced to a radiation transport code, can be used to calculate the RBE values associated to Solar Particle Events. More generally, this work emphasizes the importance of taking into account the RBE dependence on particle energy and dose when calculating equivalent doses.

Evaluation of Helium Ion Radiotherapy in Combination with Gemcitabine in Pancreatic Cancer In Vitro

BACKGROUND

Pancreatic cancer is one of the most aggressive and lethal cancers. New treatment strategies are highly warranted. Particle radiotherapy could offer a way to overcome the radioresistant nature of pancreatic cancer because of its biological and physical characteristics. Within particles, helium ions represent an attractive therapy option to achieve the highest possible conformity while at the same time protecting the surrounding normal tissue. The aim of this study was to evaluate the cytotoxic efficacy of helium ion irradiation in pancreatic cancer in vitro.

METHODS

Human pancreatic cancer cell lines AsPC-1, BxPC-3 and Panc-1 were irradiated with photons and helium ions at various doses and treated with gemcitabine. Photon irradiation was performed with a biological cabin X-ray irradiator, and helium ion irradiation was performed with a spread-out Bragg peak using the raster scanning technique at the Heidelberg Ion Beam Therapy Center (HIT). The cytotoxic effect on pancreatic cancer cells was measured with clonogenic survival. The survival curves were compared to the predicted curves that were calculated via the modified microdosimetric kinetic model (mMKM).

RESULTS

The experimental relative biological effectiveness (RBE) of helium ion irradiation ranged from 1.0 to 1.7. The predicted survival curves obtained via mMKM calculations matched the experimental survival curves. Mainly additive cytotoxic effects were observed for the cell lines AsPC-1, BxPC-3 and Panc-1.

CONCLUSION

Our results demonstrate the cytotoxic efficacy of helium ion radiotherapy in pancreatic cancer in vitro as well as the capability of mMKM calculation and its value for biological plan optimization in helium ion therapy for pancreatic cancer. A combined treatment of helium irradiation and chemotherapy with gemcitabine leads to mainly additive cytotoxic effects in pancreatic cancer cell lines. The data generated in this study may serve as the radiobiological basis for future experimental and clinical works using helium ion radiotherapy in pancreatic cancer treatment.

Proton relative biological effectiveness for the induction of DNA double strand breaks based on Geant4

Uncertainties in the relative biological effectiveness (RBE) of proton remains a major barrier to the biological optimization of proton therapy. A large amount of experimental data suggest that proton RBE is variable. As an evolving Monte Carlo code toolkit, Geant4-DNA is able to simulate the initial DNA damage caused by particle beams through physical and chemical interactions at the nanometer scale over a short period of time. This contributes to evaluating the radiobiological effects induced by ionizing radiation. Based on the Geant4-DNA toolkit, this study constructed a DNA geometric model containing 6.32Gbp, simulated the relationship between radiochemical yields (G-values) and their corresponding chemical constructors, and calculated a detailed calculation of the sources of damage and the complexity of damage in DNA strand breaks. The damage model constructed in this study can simulate the relative biological effectiveness (RBE) in the proton Bragg peak region. The results indicate that: (1) When the electron energy is below 400 keV, the yield of OH·account for 18.1% to 25.3% of the total water radiolysis yields. (2) Under the influence of histone clearance function, the yield of indirect damage account for over 72.93% of the yield of DNA strand breaks (SBs). When linear energy transfer (LET) increased from 29.79 (keV/μm) to 64.29 (keV/μm), the yield of double strand breaks (DSB) increased from 17.27% to 32.65%. (3) By investigating the effect of proton Bragg peak depth on the yield of direct DSB (DSBdirect) and total DSB (DSBtotal), theRBEDSBtotandRBEDSBdirlevels of cells show that the RBE value of protons reaches 2.2 in the Bragg peak region.

Equivalent uniform RBE-weighted dose in eye plaque brachytherapy

BACKGROUND

Brachytherapy for ocular melanoma is based on the application of eye plaques with different spatial dose nonuniformity, time-dependent dose rates and relative biological effectiveness (RBE).

PURPOSE

We propose a parameter called the equivalent uniform RBE-weighted dose (EUDRBE) that can be used for quantitative characterization of integrated cell survival in radiotherapy modalities with the variable RBE, dose nonuniformity and dose rate. The EUDRBE is applied to brachytherapy with 125I eye plaques designed by the Collaborative Ocular Melanoma Study (COMS).

METHODS

The EUDRBE is defined as the uniform dose distribution with RBE = 1 that causes equal cell survival for a given nonuniform dose distribution with the variable RBE > 1. The EUDRBE can be used for comparison of cell survival for nonuniform dose distributions with different RBE, because they are compared to the reference dose with RBE = 1. The EUDRBE is applied to brachytherapy with 125I COMS eye plaques that are characterized by a steep dose gradient in tumor base-apex direction, protracted irradiation during time intervals of 3-8 days, and variable dose-rate dependent RBE with a maximum of about 1.4. The simulations are based on dose of 85 Gy prescribed to the farthest intraocular extent of the tumor (tumor apex). To compute the EUDRBE in eye plaque brachytherapy and correct for protracted irradiation, the distributions of physical dose have been converted to non-uniform distributions of biologically effective dose (BED) to include the biological effects of sublethal cellular repair, Our radiobiological analysis considers the combined effects of different time-dependent dose rates, spatial dose non-uniformity, dose fractionation and different RBE and can be used to derive optimized dose regimens brachytherapy.

RESULTS

Our simulations show that the EUDRBE increases with the prescription depths and the maximum increase may achieve 6% for the tumor height of 12 mm. This effect stems from a steep dose gradient within the tumor that increases with the prescription depth. The simulations also show that the EUDRBE increase may achieve 12% with increasing the dose rate when implant duration decreases. The combined effect of dose nonuniformity and dose rate may change the EUDRBE up to 18% for the same dose prescription of 85 Gy to tumor apex. The absolute dose range of 48-61 Gy (RBE) for the EUDRBE computed using 4 or 5 fractions is comparable to the dose prescriptions used in stereotactic body radiation therapy (SBRT) with megavoltage X-rays (RBE = 1) for different cancers. The tumor control probabilities in SBRT and eye plaque brachytherapy are very similar at the level of 80% or higher that support the hypothesis that the selected approximations for the EUDRBE are valid.

CONCLUSIONS

The computed range of the EUDRBE in 125I COMS eye plaque brachytherapy suggests that the selected models and hypotheses are acceptable. The EUDRBE can be useful for analysis of treatment outcomes and comparison of different dose regimens in eye plaque brachytherapy.

Primary study of the relative and compound biological effectiveness model for boron neutron capture therapy based on nanodosimetry

BACKGROUND

The current radiobiological model employed for boron neutron capture therapy (BNCT) treatment planning, which relies on microdosimetry, fails to provide an accurate representation the biological effects of BNCT. The precision in calculating the relative biological effectiveness (RBE) and compound biological effectiveness (CBE) plays a pivotal role in determining the therapeutic efficacy of BNCT. Therefore, this study focuses on how to improve the accuracy of the biological effects of BNCT.

PURPOSE

The purpose of this study is to propose new radiation biology models based on nanodosimetry to accurately assess RBE and CBE for BNCT.

METHODS

Nanodosimetry, rooted in ionization cluster size distributions (ICSD), introduces a novel approach to characterize radiation quality by effectively delineating RBE through the ion track structure at the nanoscale. In the context of prior research, this study presents a computational model for the nanoscale assessment of RBE and CBE. We establish a simplified model of DNA chromatin fiber using the Monte Carlo code TOPAS-nBio to evaluate the applicability of ICSD to BNCT and compute nanodosimetric parameters.

RESULTS

Our investigation reveals that both homogeneous and heterogeneous nanodosimetric parameters, as well as the corresponding biological model coefficients α and β, along with RBE values, exhibit variations in response to varying intracellular 10B concentrations. Notably, the nanodosimetric parameterM 1 C 2 $M_1^{{{\mathrm{C}}}_2}$ effectively captures the fluctuations in model coefficients α and RBE.

CONCLUSION

Our model facilitates a nanoscale analysis of BNCT, enabling predictions of nanodosimetric quantities for secondary ions as well as RBE, CBE, and other essential biological metrics related to the distribution of boron. This contribution significantly enhances the precision of RBE calculations and holds substantial promise for future applications in treatment planning.

Hypoxia-informed RBE-weighted beam orientation optimization for intensity modulated proton therapy

BACKGROUND

Variable relative biological effectiveness (RBE) models in treatment planning have been proposed to optimize the therapeutic ratio of proton therapy. It has been reported that proton RBE decreases with increasing tumor oxygen level, offering an opportunity to address hypoxia-related radioresistance with RBE-weighted optimization.

PURPOSE

Here, we obtain a voxel-level estimation of partial oxygen pressure to weigh RBE values in a single biologically informed beam orientation optimization (BOO) algorithm.

METHODS

Three glioblastoma patients with [18 F]-fluoromisonidazole (FMISO)-PET/CT images were selected from the institutional database. Oxygen values were derived from tracer uptake using a nonlinear least squares curve fitting. McNamara RBE, calculated from proton dose, was then weighed using oxygen enhancement ratios (OER) for each voxel and incorporated into the dose fidelity term of the BOO algorithm. The nonlinear optimization problem was solved using a split-Bregman approach, with FISTA as the solver. The proposed hypoxia informed RBE-weighted method (HypRBE) was compared to dose fidelity terms using the constant RBE of 1.1 (cRBE) and the normoxic McNamara RBE model (RegRBE). Tumor homogeneity index (HI), maximum biological dose (Dmax), and D95%, as well as OAR therapeutic index (TI = gEUDCTV /gEUDOAR ) were evaluated along with worst-case statistics after normalization to normal tissue isotoxicity.

RESULTS

Compared to [cRBE, RegRBE], HypRBE increased tumor HI, Dmax, and D95% across all plans by on average [31.3%, 31.8%], [48.6%, 27.1%], and [50.4%, 23.8%], respectively. In the worst-case scenario, the parameters increase on average by [12.5%, 14.7%], [7.3%,-8.9%], and [22.3%, 2.1%]. Despite increased OAR Dmean and Dmax by [8.0%, 3.0%] and [13.1%, -0.1%], HypRBE increased average TI by [22.0%, 21.1%]. Worst-case OAR Dmean, Dmax, and TI worsened by [17.9%, 4.3%], [24.5%, -1.2%], and [9.6%, 10.5%], but in the best cases, HypRBE escalates tumor coverage significantly without compromising OAR dose, increasing the therapeutic ratio.

CONCLUSIONS

We have developed an optimization algorithm whose dose fidelity term accounts for hypoxia-informed RBE values. We have shown that HypRBE selects bE:\Alok\aaeams better suited to deliver high physical dose to low RBE, hypoxic tumor regions while sparing the radiosensitive normal tissue.

The Effects of Particle LET and Fluence on the Complexity and Frequency of Clustered DNA Damage

Motivation

Clustered DNA-lesions are predominantly induced by ionizing radiation, particularly by high-LET particles, and considered as lethal damage. Quantification of this specific type of damage as a function of radiation parameters such as LET, dose rate, dose, and particle type can be informative for the prediction of biological outcome in radiobiological studies. This study investigated the induction and complexity of clustered DNA damage for three different types of particles at an LET range of 0.5-250 keV/μm.

Methods

Nanometric volumes (36.0 nm3) of 15 base-pair DNA with its hydration shell was modeled. Electron, proton, and alpha particles at various energies were simulated to irradiate the nanometric volumes. The number of ionization events, low-energy electron spectra, and chemical yields for the formation of °OH, H°, e aq - , and H2O2 were calculated for each particle as a function of LET. Single- and double-strand breaks (SSB and DSB), base release, and clustered DNA-lesions were computed from the Monte-Carlo based quantification of the reactive species and measured yields of the species responsible for the DNA lesion formation.

Results

The total amount of DNA damage depends on particle type and LET. The number of ionization events underestimates the quantity of DNA damage at LETs higher than 10 keV/μm. Minimum LETs of 9.4 and 11.5 keV/μm are required to induce clustered damage by a single track of proton and alpha particles, respectively. For a given radiation dose, an increase in LET reduces the number of particle tracks, leading to more complex clustered DNA damage, but a smaller number of separated clustered damage sites.

Conclusions

The dependency of the number and the complexity of clustered DNA damage on LET and fluence suggests that the quantification of this damage can be a useful method for the estimation of the biological effectiveness of radiation. These results also suggest that medium-LET particles are more appropriate for the treatment of bulk targets, whereas high-LET particles can be more effective for small targets.

Integrating microdosimetricin vitroRBE models for particle therapy into TOPAS MC using the MicrOdosimetry-based modeliNg for RBE ASsessment (MONAS) tool

Objective.In this paper, we present MONAS (MicrOdosimetry-based modelliNg for relative biological effectiveness (RBE) ASsessment) toolkit. MONAS is a TOPAS Monte Carlo extension, that combines simulations of microdosimetric distributions with radiobiological microdosimetry-based models for predicting cell survival curves and dose-dependent RBE.Approach.MONAS expands TOPAS microdosimetric extension, by including novel specific energy scorers to calculate the single- and multi-event specific energy microdosimetric distributions at different micrometer scales. These spectra are used as physical input to three different formulations of themicrodosimetric kinetic model, and to thegeneralized stochastic microdosimetric model(GSM2), to predict dose-dependent cell survival fraction and RBE. MONAS predictions are then validated against experimental microdosimetric spectra andin vitrosurvival fraction data. To show the MONAS features, we present two different applications of the code: (i) the depth-RBE curve calculation from a passively scattered proton SOBP and monoenergetic12C-ion beam by using experimentally validated spectra as physical input, and (ii) the calculation of the 3D RBE distribution on a real head and neck patient geometry treated with protons.Main results.MONAS can estimate dose-dependent RBE and cell survival curves from experimentally validated microdosimetric spectra with four clinically relevant radiobiological models. From the radiobiological characterization of a proton SOBP and12C fields, we observe the well-known trend of increasing RBE values at the distal edge of the radiation field. The 3D RBE map calculated confirmed the trend observed in the analysis of the SOBP, with the highest RBE values found in the distal edge of the target.Significance.MONAS extension offers a comprehensive microdosimetry-based framework for assessing the biological effects of particle radiation in both research and clinical environments, pushing closer the experimental physics-based description to the biological damage assessment, contributing to bridging the gap between a microdosimetric description of the radiation field and its application in proton therapy treatment with variable RBE.

Incorporating variable RBE in IMPT optimization for ependymoma

PURPOSE

To study the dosimetric impact of incorporating variable relative biological effectiveness (RBE) of protons in optimizing intensity-modulated proton therapy (IMPT) treatment plans and to compare it with conventional constant RBE optimization and linear energy transfer (LET)-based optimization.

METHODS

This study included 10 pediatric ependymoma patients with challenging anatomical features for treatment planning. Four plans were generated for each patient according to different optimization strategies: (1) constant RBE optimization (ConstRBEopt) considering standard-of-care dose requirements; (2) LET optimization (LETopt) using a composite cost function simultaneously optimizing dose-averaged LET (LETd ) and dose; (3) variable RBE optimization (VarRBEopt) using a recent phenomenological RBE model developed by McNamara et al.; and (4) hybrid RBE optimization (hRBEopt) assuming constant RBE for the target and variable RBE for organs at risk. By normalizing each plan to obtain the same target coverage in either constant or variable RBE, we compared dose, LETd , LET-weighted dose, and equivalent uniform dose between the different optimization approaches.

RESULTS

We found that the LETopt plans consistently achieved increased LET in tumor targets and similar or decreased LET in critical organs compared to other plans. On average, the VarRBEopt plans achieved lower mean and maximum doses with both constant and variable RBE in the brainstem and spinal cord for all 10 patients. To compensate for the underdosing of targets with 1.1 RBE for the VarRBEopt plans, the hRBEopt plans achieved higher physical dose in targets and reduced mean and especially maximum variable RBE doses compared to the ConstRBEopt and LETopt plans.

CONCLUSION

We demonstrated the feasibility of directly incorporating variable RBE models in IMPT optimization. A hybrid RBE optimization strategy showed potential for clinical implementation by maintaining all current dose limits and reducing the incidence of high RBE in critical normal tissues in ependymoma patients.

Most DNA repair defects do not modify the relationship between relative biological effectiveness and linear energy transfer in CRISPR-edited cells

BACKGROUND

Cancer is a highly heterogeneous disease, driven by frequent genetic alterations which have significant effects on radiosensitivity. However, radiotherapy for a given cancer type is typically given with a standard dose determined from population-level trials. As a result, a proportion of patients are under- or over-dosed, reducing the clinical benefit of radiotherapy. Biological optimization would not only allow individual dose prescription but also a more efficient allocation of limited resources, such as proton and carbon ion therapy. Proton and ion radiotherapy offer an advantage over photons due to their elevated Relative Biological Effectiveness (RBE) resulting from their elevated Linear Energy Transfer (LET). Despite significant interest in optimizing LET by tailoring radiotherapy plans, RBE's genetic dependence remains unclear.

PURPOSE

The aim of this study is to better define the RBE/LET relationship in a panel of cell lines with different defects in DSB repair pathways, but otherwise identical biological features and genetic background to isolate these effects.

METHODS

Normal human cells (RPE1), genetically modified to introduce defects in DNA double-strand break (DSB) repair genes, ATM, BRCA1, DCLRE1C, LIG4, PRKDC and TP53, were used to map the RBE-LET relationship. Cell survival was measured with clonogenic assays after exposure to photons, protons (LET 1 and 12 keV/µm) and alpha particles (129 keV/µm). Gene knockout sensitizer enhancement ratio (SER) values were calculated as the ratio of the mean inactivation dose (MID) of wild-type cells to repair-deficient cells, and RBE values were calculated as the ratio of the MID of X-ray and particle irradiated cells. 53BP1 foci were used to quantify radiation-induced DSBs and their repair following irradiation.

RESULTS

Deletion of NHEJ genes had the greatest impact on photon sensitivity (ATM-/- SER = 2.0 and Lig4-/- SER = 1.8), with genes associated with HR having smaller effects (BRCA1-/- SER = 1.2). Wild-type cells showed RBEs of 1.1, 1.3, 5.0 for low- and high-LET protons and alpha particles respectively. SERs for different genes were independent of LET, apart from NHEJ knockouts which proved to be markedly hypersensitive across all tested LETs. Due to this hypersensitivity, the impact of high LET was reduced in cell models lacking the NHEJ repair pathway. HR-defective cells had moderately increased sensitivity across all tested LETs, but, notably, the contribution of HR pathway to survival appeared independent of LET. Analysis of 53BP1 foci shows that NHEJ-defective cells had the least DSB repair capacity after low LET exposure, and no visible repair after high LET exposure. HR-defective cells also had slower repair kinetics, but the impact of HR defects is not as severe as NHEJ defects.

CONCLUSIONS

DSB repair defects, particularly in NHEJ, conferred significant radiosensitivity across all LETs. This sensitization appeared independent of LET, suggesting that the contribution of different DNA repair pathways to survival does not depend on radiation quality.

Variable RBE in proton radiotherapy: a comparative study with the predictive Mayo Clinic Florida microdosimetric kinetic model and phenomenological models of cell survival

Objectives. (1) To examine to what extent the cell- and exposure- specific information neglected in the phenomenological proton relative biological effectiveness (RBE) models could influence the computed RBE in proton therapy. (2) To explore similarities and differences in the formalism and the results between the linear energy transfer (LET)-based phenomenological proton RBE models and the microdosimetry-based Mayo Clinic Florida microdosimetric kinetic model (MCF MKM). (3) To investigate how the relationship between the RBE and the dose-mean proton LET is affected by the proton energy spectrum and the secondary fragments.Approach. We systematically compared six selected phenomenological proton RBE models with the MCF MKM in track-segment simulations, monoenergetic proton beams in a water phantom, and two spread-out Bragg peaks. A representative comparison within vitrodata for human glioblastoma cells (U87 cell line) is also included.Main results. Marked differences were observed between the results of the phenomenological proton RBE models, as reported in previous studies. The dispersion of these models' results was found to be comparable to the spread in the MCF MKM results obtained by varying the cell-specific parameters neglected in the phenomenological models. Furthermore, while single cell-specific correlation between RBE and the dose-mean proton LET seems reasonable above 2 keVμm-1, caution is necessary at lower LET values due to the relevant contribution of secondary fragments. The comparison within vitrodata demonstrates comparable agreement between the MCF MKM predictions and the results of the phenomenological models.Significance. The study highlights the importance of considering cell-specific characteristics and detailed radiation quality information for accurate RBE calculations in proton therapy. Furthermore, these results provide confidence in the use of the MCF MKM for clonogenic survival RBE calculations in proton therapy, offering a more mechanistic approach compared to phenomenological models.

Effects of tissue heterogeneity and comparisons of collapsed cone and Monte Carlo fast neutron patient dosimetry using the University of Washington clinical neutron therapy system (CNTS)

Fast neutron therapy is a high linear energy transfer (LET) radiation treatment modality offering advantages over low LET radiations. Multileaf collimator technology reduces normal-tissue dose (toxicity) and makes neutron therapy more comparable to MV x-ray treatments. Published clinical-trial and other experiences with fast neutron therapy are reported. Early comparative studies failed to consider differences in target-dose spatial conformality between x-ray and neutron treatments, which is especially important for organs-at-risk close to tumor targets. Treatments planning systems (TPS) for high-energy neutrons lag behind TPS tools for MV x-rays, creating challenges for comparative studies of clinical outcomes. A previously published Monte Carlo model of the University of Washington (UW) Clinical Neutron Therapy System (CNTS) is refined and integrated with the RayStation TPS as an external dose planning/verification tool. The collapsed cone (CC) dose calculations in the TPS are based on measured dose profiles and output factors in water, with the absolute dose determined using a tissue-equivalent ionization chamber. For comparison, independent (external) Monte Carlo simulation computes dose on a voxel-by-voxel basis using an atlas that maps Hounsfield Unit (HU) numbers to elemental composition and density. Although the CC algorithm in the TPS accurately computes neutron dose to water compared to Monte Carlo calculations, calculated dose to water differs from bone or tissue depending largely on hydrogen content. Therefore, the elemental composition of tissue and bone, rather than the material or electron density, affects fast neutron dose. While the CC algorithm suffices for reproducible patient dosimetry in fast neutron therapy, adopting methods that consider tissue heterogeneity would enhance patient-specific neutron dose accuracy relative to national standards for other types of ionizing radiation. Corrections for tissue composition have a significant impact on absolute dose and the relative biological effectiveness (RBE) of neutron treatments compared to other radiation types (MV x-rays, protons, and carbon ions).

Anti-tumor Effect of High Doses of Carbon Ions and X-Rays during Irradiation of Ehrlich Ascites Carcinoma Cells Ex Vivo

The effect of carbon ions (12C) with the energy of 400 MeV/nucleon on the dynamics of induction and growth rate of solid tumors in mice under irradiation of Ehrlich ascites carcinoma cells (EAC) ex vivo at doses of 5-30 Gy relative to the action of equally effective doses of X-ray radiation was studied. The dynamics of tumor induction under the action of 12C and X-rays had a similar character and depended on the dose during 3 months of observation. The value of the latent period, both when irradiating cells with 12C and X-ray, increased with increasing dose, and the interval for tumor induction decreased. The rate of tumor growth after ex vivo irradiation of EAC cells was independent of either dose or type of radiation. The dose at which EAC tumors are not induced within 90 days was 30 Gy for carbon ions and 60 Gy for X-rays. The value of the relative biological effectiveness of carbon ions, calculated from an equally effective dose of 50% probability of tumors, was 2.59.

X-ray and DNA Damage: Limitations of the Dose as a Parameter for In Vitro Studies

A century of studies has demonstrated that the magnitude of a radiation dose determines the extent of its biological effect. However, different types of radiation show different levels of effectiveness. Although all types of X-rays are usually considered to be equivalent, several authors have demonstrated an inverse relationship between photon energy and the biological effectiveness of the X-ray. Nonetheless, the differences among 50-100 keV X-rays are usually considered absent. However, comparing different types of X-rays with different energies is not easy since they are often used with different dose rates, and the latter can be a confounding factor. We compared the biological effectiveness of X-rays with different photon energies but with the same dose rate. Moreover, we also studied X-ray with different dose rates but the same photon energy. Biological effectiveness was assessed measuring DNA damage and cell survival. We confirmed that both the dose rate and photon energy influence the effectiveness of an X-ray. Moreover, we observed that differences in the 50-100 keV range are detectable after controlling for dose-rate variations. Our results, confirming those of previous studies in a more consistent way (and accompanied by hypotheses on the importance of the number of incident photons), underline the limitations of using the dose as the sole parameter for in vitro studies.

Microdosimetry-based investigation of biological effectiveness of252Cf brachytherapy source: TOPAS Monte Carlo study

Objective.To investigate biological effectiveness of252Cf brachytherapy source using Monte Carlo-calculated microdosimetric distributions.Approach.252Cf source capsule was placed at the center of the spherical water phantom and phase-space data were scored as a function of radial distance in water (R= 1-5 cm) using TOPAS Monte Carlo code. The phase-space data were used to calculate microdosimetric distributions at 1μm site size. Using these distributions, Relative Biological Effectiveness (RBE), mean quality factor (Q̅) and Oxygen Enhancement Ratio (OER) were calculated as a function ofR.Main results.The overall shapes of the microdosimetric distributions are comparable at all the radial distances in water. However, slight variation in the bin-wise yield is observed withR. RBE,Q̅and OER are insensitive to R over the range 1-5 cm. Microdosimetric kinetic model based RBE values are about 2.3 and 2.8 for HSG tumour cells and V79 cells, respectively, whereas biological weighting function-based RBE is about 2.8. ICRP 60 and ICRU 40 recommendation-basedQ̅values are about 14.5 and 16, respectively. Theory of dual radiation action based RBE is 11.4. The calculated value of OER is 1.6.Significance.This study demonstrates the relative insensitivity of RBE,Q̅and OER radially away from the252Cf source along the distances of 1-5 cm in water.

Plan robustness and RBE influence for proton dose painting by numbers for head and neck cancers

PURPOSE

To investigate the feasibility of dose painting by numbers (DPBN) with respect to robustness for proton therapy for head and neck cancers (HNC), and to study the influence of variable RBE on the TCP and OAR dose burden.

METHODS AND MATERIALS

Data for 19 patients who have been scanned pretreatment with PET-FDG and subsequently treated with photon therapy were used in the study. A dose response model developed for photon therapy was implemented in a TPS, allowing DPBN plans to be created. Conventional homogeneous dose and DPBN plans were created for each patient, optimized with either fixed RBE = 1.1 or a variable RBE model. Robust optimization was used to create clinically acceptable plans. To estimate the maximum potential loss in TCP due to actual SUV variations from the pre-treatment imaging, we applied a test case with randomized SUV distribution.

RESULTS

Regardless of the use of variable RBE for optimization or evaluation, a statistically significant increase (p < 0.001) in TCP was found for DPBN plans as compared to homogeneous dose plans. Randomizing the SUV distribution decreased the TCP for all plans. A correlation between TCP increase and variance of the SUV distribution and target volume was also found.

CONCLUSION

DPBN for protons and HNC is feasible and could lead to a TCP gain. Risks associated with the temporal variation of SUV distributions could be mitigated by imposing minimum doses to targets. The correlation found between TCP increase and SUV variance and target volume may be used for patient selection.

Treatment of primary or recurrent non-resectable pancreatic cancer with proton beam irradiation combined with gemcitabine-based chemotherapy

BACKGROUND

Pancreatic cancer accounts for around 4.6% of cancers deaths worldwide per year. Despite many advances in treatment regimes, the prognosis is still poor. Only 20% of tumors are primarily resectable. Recurrences-both with distant metastasis as well as locoregional-are frequent. For patients with primary nonresectable localized disease or localized recurrences, we offered chemoradiation to achieve local control over a long period of time. We here report our results on combined chemoradiation of pancreatic tumors and local recurrences using proton beam therapy.

MATERIALS AND METHODS

We report on 25 patients with localized nonresectable pancreatic cancer (15 patients) or local recurrent disease (10 patients). All patients were treated with combined proton radiochemotherapy. Overall survival, progression-free survival, local control, and treatment-related toxicity were analyzed using statistically methods.

RESULTS

Median RT dose was 54.0 Gy (RBE) for proton irradiation. The toxicity of treatment was acceptable. Four CTCAE grade III and IV adverse events (bone marrow disfunction, gastrointestinal [GI] disorders, stent dislocation, myocardial infarction) were recorded during or directly after the end of radiotherapy; two of them were related to combined chemoradiation (bone marrow disfunction, GI disorders). Six weeks after radiotherapy, one additional grade IV toxicity was reported (ileus, caused by peritoneal carcinomatosis, not treatment related). The median progression-free survival was 5.9 months and median overall survival was 11.0 months. The pretherapy CA19‑9 level was a statistically significant prognostic factor for enhanced overall survival. Local control at 6 months and 12 months were determined to be 86% and 80%, respectively.

CONCLUSION

Combined proton chemoradiation leads to high local control rates. Unfortunately, PFS and OS are driven by distant metastasis and were not improved compared to historical data and reports. With this in mind, enhanced chemotherapeutical regimes, in combination with local irradiation, should be evaluated.

Experimentally determined relative biological effectiveness of cyclotron-based epithermal neutrons designed for clinical BNCT: in vitro study

A neutron beam for boron neutron capture therapy (BNCT) of deep-seated tumours is designed to maintain a high flux of epithermal neutrons, while keeping the thermal and fast neutron component as low as possible. These neutrons (thermal and fast) have a high relative biological effectiveness in comparison with high energy photon beams used for conventional X-ray radiotherapy. In the past, neutrons for the purpose of BNCT were generated using nuclear reactors. However, there are various challenges that arise when installing a reactor in a hospital environment. From 2006, the Kyoto University Research Reactor Institute, in collaboration with Sumitomo Heavy Industries, began the development of an accelerator-based neutron source for clinical BNCT in a bid to overcome the shortcomings of a nuclear reactor-based neutron source. Following installation and beam performance testing, in vitro studies were performed to assess the biological effect of the neutron beam. Four different cell lines were prepared and irradiated using the accelerator-based neutron source. Following neutron and gamma ray irradiation, the survival curve for each cell line was calculated. The biological end point to determine the relative biological effectiveness (RBE) was set to 10% cell survival, and the D10 for each cell line was determined. The RBE of the accelerator-based neutron beam was evaluated to be 2.62.

Estimation of relative biological effectiveness of 225Ac compared to 177Lu during [225Ac]Ac-PSMA and [177Lu]Lu-PSMA radiopharmaceutical therapy using TOPAS/TOPAS-nBio/MEDRAS

AIM

Over recent years, [225Ac]Ac-PSMA and [177Lu]Lu-PSMA radiopharmaceutical therapy have evolved as a promising treatment option for advanced prostate cancer. Especially for alpha particle emitter treatments, there is still a need for improving dosimetry, which requires accurate values of relative biological effectiveness (RBE). To achieve that, consideration of DNA damages in the cell nucleus and knowledge of the energy deposition in the location of the DNA at the nanometer scale are required. Monte Carlo particle track structure simulations provide access to interactions at this level. The aim of this study was to estimate the RBE of 225Ac compared to 177Lu. The initial damage distribution after radionuclide decay and the residual damage after DNA repair were considered.

METHODS

This study employed the TOol for PArtcile Simulation (TOPAS) based on the Geant4 simulation toolkit. Simulation of the nuclear DNA and damage scoring were performed using the TOPAS-nBio extension of TOPAS. DNA repair was modeled utilizing the Python-based program MEDRAS (Mechanistic DNA Repair and Survival). Five different cell geometries of equal volume and two radionuclide internalization assumptions as well as two cell arrangement scenarios were investigated. The radionuclide activity (number of source points) was adopted based on SPECT images of patients undergoing the above-mentioned therapies.

RESULTS

Based on the simulated dose-effect curves, the RBE of 225Ac compared to 177Lu was determined in a wide range of absorbed doses to the nucleus. In the case of spherical geometry, 3D cell arrangement and full radionuclide internalization, the RBE based on the initial damage had a constant value of approximately 2.14. Accounting for damage repair resulted in RBE values ranging between 9.38 and 1.46 for 225Ac absorbed doses to the nucleus between 0 and 50 Gy, respectively.

CONCLUSION

In this work, the consideration of DNA repair of the damage from [225Ac]Ac-PSMA and [177Lu]Lu-PSMA revealed a dose dependency of the RBE. Hence, this work suggested that DNA repair is an important aspect to understand response to different radiation qualities.

Influence of beam pruning techniques on LET and RBE in proton arc therapy

Introduction

Proton arc therapy (PAT) is an emerging treatment modality that holds promise to improve target volume coverage and reduce linear energy transfer (LET) in organs at risk. We aimed to investigate if pruning the highest energy layers in each beam direction could increase the LET in the target and reduce LET in tissue and organs at risk (OAR) surrounding the target volume, thus reducing the relative biological effectiveness (RBE)-weighted dose and sparing healthy tissue.

Methods

PAT plans for a germinoma, an ependymoma and a rhabdomyosarcoma patient were created in the Eclipse treatment planning system with a prescribed dose of 54 Gy(RBE) using a constant RBE of 1.1 (RBE1.1). The PAT plans was pruned for high energy spots, creating several PAT plans with different amounts of pruning while maintaining tumor coverage, denoted PX-PAT plans, where X represents the amount of pruning. All plans were recalculated in the FLUKA Monte Carlo software, and the LET, physical dose, and variable RBE-weighted dose from the phenomenological Rørvik (ROR) model and an LET weighted dose (LWD) model were evaluated.

Results and discussion

For the germinoma case, all plans but the P6-PAT reduced the mean RBE-weighted dose to the surrounding healthy tissue compared to the PAT plan. The LET was increasingly higher within the PTV for each pruning iteration, where the mean LET from the P6-PAT plan was 1.5 keV / μm higher than for the PAT plan, while the P4- and P5-PAT plans provided an increase of 0.4 and 0.7 keV / μm , respectively. The other plans increased the LET by a smaller margin compared to the PAT plan. Likewise, the LET values to the healthy tissue were reduced for each degree of pruning. Similar results were found for the ependymoma and the rhabdomyosarcoma case. We demonstrated a PAT pruning technique that can increase both LET and RBE in the target volume and at the same time decreased values in healthy tissue, without affecting the target volume dose coverage.

Ionization detail parameters and cluster dose: a mathematical model for selection of nanodosimetric quantities for use in treatment planning in charged particle radiotherapy

Objective. To propose a mathematical model for applying ionization detail (ID), the detailed spatial distribution of ionization along a particle track, to proton and ion beam radiotherapy treatment planning (RTP).Approach. Our model provides for selection of preferred ID parameters (Ip) for RTP, that associate closest to biological effects. Cluster dose is proposed to bridge the large gap between nanoscopicIpand macroscopic RTP. Selection ofIpis demonstrated using published cell survival measurements for protons through argon, comparing results for nineteenIp:Nk,k= 2, 3, …, 10, the number of ionizations in clusters ofkor more per particle, andFk,k= 1, 2, …, 10, the number of clusters ofkor more per particle. We then describe application of the model to ID-based RTP and propose a path to clinical translation.Main results. The preferredIpwereN4andF5for aerobic cells,N5andF7for hypoxic cells. Significant differences were found in cell survival for beams having the same LET or the preferredNk. Conversely, there was no significant difference forF5for aerobic cells andF7for hypoxic cells, regardless of ion beam atomic number or energy. Further, cells irradiated with the same cluster dose for theseIphad the same cell survival. Based on these preliminary results and other compelling results in nanodosimetry, it is reasonable to assert thatIpexist that are more closely associated with biological effects than current LET-based approaches and microdosimetric RBE-based models used in particle RTP. However, more biological variables such as cell line and cycle phase, as well as ion beam pulse structure and rate still need investigation.Significance. Our model provides a practical means to select preferredIpfrom radiobiological data, and to convertIpto the macroscopic cluster dose for particle RTP.

GPU-accelerated calculation of proton microdosimetric spectra as a function of target size, proton energy, and bounding volume size

Objective.Shortcomings of dose-averaged linear energy transfer (LETD), the quantity which is most commonly used to quantify proton relative biological effectiveness, have long been recognized. Microdosimetric spectra may overcome the limitations of LETDbut are extremely computationally demanding to calculate. A systematic library of lineal energy spectra for monoenergetic protons could enable rapid determination of microdosimetric spectra in a clinical environment. The objective of this work was to calculate and validate such a library of lineal energy spectra.Approach. SuperTrack, a GPU-accelerated CUDA/C++ based application, was developed to superimpose tracks calculated using Geant4 onto targets of interest and to compute microdosimetric spectra. Lineal energy spectra of protons with energies from 0.1 to 100 MeV were determined in spherical targets of diameters from 1 nm to 10μm and in bounding voxels with side lengths of 5μm and 3 mm.Main results.Compared to an analogous Geant4-based application, SuperTrack is up to 3500 times more computationally efficient if each track is resampled 1000 times. Dose spectra of lineal energy and dose-mean lineal energy calculated with SuperTrack were consistent with values published in the literature and with comparison to a Geant4 simulation. Using SuperTrack, we developed the largest known library of proton microdosimetric spectra as a function of primary proton energy, target size, and bounding volume size.Significance. SuperTrack greatly increases the computational efficiency of the calculation of microdosimetric spectra. The elevated lineal energy observed in a 3 mm side length bounding volume suggests that lineal energy spectra determined experimentally or computed in small bounding volumes may not be representative of the lineal energy spectra in voxels of a dose calculation grid. The library of lineal energy spectra calculated in this work could be integrated with a treatment planning system for rapid determination of lineal energy spectra in patient geometries.

Comparing the DNA-damage RBE of intraoperative and conventional electron beams using a hybrid simulation approach

PURPOSE

Employing electron beam for radiotherapy purposes now has been established as one of the standard cancer treatment modalities. Both dedicated intraoperative and conventional electron beams can be employed in patient irradiation. Due to the differences between accelerating structure and electron beam delivery of dedicated intraoperative radiotherapy (IORT) machines and conventional ones, the initial energy spectra of the produced electron beam by these machines may be different. Accordingly, this study aims to evaluate whether these spectral differences can affect the relevant relative biological effectiveness (RBE) values of intraoperative and conventional electron beams.

MATERIALS AND METHODS

A hybrid Monte Carlo simulation approach was considered. At first, the head LIAC12 machine (as an IORT accelerator) and Varian 2100C/D (as a conventional accelerator) were simulated by MCNPX code and electron energy spectra at different depths and off-axis distances were scored for two nominal electron energies of 6 and 12 MeV at the field sizes of 6 and 10 cm. Then, the calculated spectra were imported to MCDS code to estimate the induced DNA-damage RBE values. Finally, the obtained RBE values for intraoperative and conventional electron beams were compared together.

RESULTS

The results showed that the RBE values of the intraoperative electron beam are superior to those obtained for conventional electron beam at the same energy/field size combination. Variations of the depth can regularly affect the RBE value for both conventional and intraoperative electron beams, while no ordered variation trend was observed for RBE with changing the off-axis distance. Variations of electron energy and field size can also influence the RBE value for both types of studied electron beams.

CONCLUSIONS

From the results, it can be concluded the structural differences between the dedicated IORT and conventional Linacs can lead to distinct initial electron energy spectra for intraoperative and conventional electron beams. These physical differences can finally lead to different RBE values for intraoperative and conventional electron beams at the same energy and field size.

Relative biological effectiveness of low-energy X-rays (25 kV) in mutant p53 cancer cells

Low-energy X-rays as used in radiation therapy and diagnostics such as mammography are associated with a certain risk of promoting tumour development, especially in patients with mutations in cancer-related genes like TP53. The present study therefore addressed the relative biological effectiveness (RBE) of low-energy X-rays for two human adenocarcinoma cell lines of the breast (MDA-MB-468) and pancreas (BxPC-3) with a mutation in the TP53 gene. Clonogenic survival and cytogenetic changes in terms of micronuclei (MN) formation were determined following irradiation with 25 kV X-rays and 200 kV reference irradiation in the dose range of 1-8 Gy. Except the frequency of MN-containing binucleated cells (BNC) (BNC + MN/BNC) in breast cancer cells yielding an RBE between 0.6 and 0.8, both cell lines displayed dose-dependent variations of RBE values between 1 and 2 for all biological end points (cell survival, (BNC + MN/BNC), MN/BNC, MN/(BNC + MN)) with increased effectiveness of 25 kV irradiation in pancreatic compared to breast cancer cells. The results confirm previous findings indicating increased effectiveness of low-energy X-rays and underline the necessity of careful risk estimation for cancer screening programmes.

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.

Cs-131 as an experimental tool for the investigation and quantification of the radiotoxicity of intracellular Auger decays in vitro

PURPOSE

In this work, we set out to provide an experimental setup, using Cs-131, with associated dosimetry for studying relative biological effectiveness (RBE) of Auger emitters.

MATERIAL AND METHODS

Cs-131 decays by 100% electron capture producing K- (9%) and L- (80%) Auger electrons with mean energies of 26 keV and 3.5 keV, respectively, plus ≈ 9.4 very low energy electrons (<0.5 keV) per decay. Cs-131 accumulates in the cells through the Na⁺/K⁺-ATPase. By this uptake mechanism and the alkali chemistry of Cs⁺, we argue for its intracellular homogeneous distribution. Cs-131 was added to the cell culture medium of HeLa and V79 Cells. The bio-kinetics of Cs-131 (uptake, release, intracellular distribution) was examined by measuring its intracellular activity concentration over time. Taking advantage of the 100% confluent cellular monolayer, we developed a new and robust dosimetry that is entrusted to a quantity called S_C-value.

RESULTS

The S_C-values evaluated in the cell nucleus are almost independent of the nuclear size and geometry. We obtained dose-rate controlled RBE-values for intracellular Cs-131 decay. Using the γH2AX assay, the RBE was 1 for HeLa cells. Using the clonogenic cell survival, it was 3.9 for HeLa cells and 3.2 for V79 cells.

CONCLUSION

This experimental setup and dosimetry provides reliable RBE-values for Auger emitters in various cell lines.

miR-188-5p and Host MALAT1 Regulate RBE Cell Migration, Invasion, and Apoptosis via Up-regulating PSMD10 in Cholangiocarcinoma

The study is designed to explore the regulatory network that MALAT1 competitively binds with miR-188-5p to up-regulate PSMD10 to facilitate cholangiocarcinoma cell migration and invasion and suppress apoptosis. qRT-PCR and fluorescence in situ hybridization (FISH) were used to examine the expression and positive signal of MALAT1 and miR-188-5p in cholangiocarcinoma tissues and HIBEC, HCCC-9810, RBE, and QBC939 cells. Western blot, qRT-PCR, and immunohistochemistry were selected to detect PSMD10 expression in cholangiocarcinoma tissues and cell lines. Dual luciferase reporter gene assay was adopted to verify that miR-188-5p targeted MALAT1 and PSMD10. qRT-PCR, pull down, and western blot were used to examine the regulation of MALAT1-miR-188-5p-PSMD10 axis. Transwell, wound healing assay, and Tunel cell apoptosis were adopted to respectively detect the regulatory abilities of MALAT1-miR-188-5p-PSMD10 axis on cell invasion, migration, and apoptosis. Western blot was used to detect the regulation mechanism of MALAT1 on Bax, Bcl-2, and caspase-3 proteins. Nude mice subcutaneous xenograft model of cholangiocarcinoma was established to examine the impacts of MALAT1 on subcutaneous tumor growth. Immunohistochemistry was adopted to examine the positive indicator of Ki67 antibodies and SMD10 antibodies in each group. MALAT1 and PSMD10 were highly expressed in cholangiocarcinoma tissues and cell lines, while miR-188-5p was lowly expressed. MALAT1 could competitively bind to miR-188-5p, and miR-188-5p could negatively regulate PSMD10. MALAT1, In-miR-188-5p, and PSMD10 could facilitate cell invasion and migration and inhibit apoptosis, while siMALAT1, miR-188-5p, and siPSMD10 produced an opposite result. MALAT1-miR-188-5p-PSMD10 axis could promote RBE cell invasion and migration and inhibit apoptosis, whereas siMALAT1-In-miR-188-5p-siPSMD10 axis showed an opposite result. On the other hand, it was verified that up-regulation/down-regulation of MALAT1 can inhibit/promote Bax and caspase-3 proteins and promote/inhibit the expression of Bcl-2 protein. MALAT1 could facilitate subcutaneous tumor growth and enhance cell proliferation and positive signal of PSMD10, while miR-188-5p worked in an opposite direction. MALAT1 competitively binds to miR-188-5p to up-regulate mRNA translation and protein expression of PSMD10, thereby facilitating cholangiocarcinoma cell invasion and migration and inhibiting its apoptosis. However, interfering MALAT1-miR-188-5p-PSMD10 axis could inhibit the occurrence and development of cholangiocarcinoma.

Increased relative biological effectiveness and periventricular radiosensitivity in proton therapy of glioma patients

PURPOSE

Currently, there is an intense debate on variations in intra-cerebral radiosensitivity and relative biological effectiveness (RBE) in proton therapy of primary brain tumours. Here, both effects were retrospectively investigated using late radiation-induced brain injuries (RIBI) observed in follow-up after proton therapy of patients with diagnosed glioma.

METHODS

In total, 42 WHO grade 2-3 glioma patients out of a consecutive patient cohort having received (adjuvant) proton radio(chemo)therapy between 2014 and 2017 were eligible for analysis. RIBI lesions (symptomatic or clinically asymptomatic) were diagnosed and delineated on contrast-enhanced T1-weighted magnetic resonance imaging scans obtained in the first two years of follow-up. Correlation of RIBI location and occurrence with dose (D), proton dose-averaged linear energy transfer (LET) and variable RBE dose parameters were tested in voxel- and in patient-wise logistic regression analyses. Additionally, anatomical and clinical parameters were considered. Model performance was estimated through cross-validated area-under-the-curve (AUC) values.

RESULTS

In total, 64 RIBI lesions were diagnosed in 21 patients. The median time between start of proton radio(chemo)therapy and RIBI appearance was 10.2 months. Median distances of the RIBI volume centres to the cerebral ventricles and to the clinical target volume border were 2.1 mm and 1.3 mm, respectively. In voxel-wise regression, the multivariable model with D, D × LET and periventricular region (PVR) revealed the highest AUC of 0.90 (95 % confidence interval: 0.89-0.91) while the corresponding model without D × LET revealed a value of 0.84 (0.83-0.86). In patient-level analysis, the equivalent uniform dose (EUD₁₁, a = 11) in the PVR using a variable RBE was the most prominent predictor for RIBI with an AUC of 0.63 (0.32-0.90).

CONCLUSIONS

In this glioma cohort, an increased radiosensitivity within the PVR was observed as well as a spatial correlation of RIBI with an increased RBE. Both need to be considered when delivering radio(chemo)therapy using proton beams.

Calibration and evaluation of the relative biological effectiveness for carbon-ion radiotherapy in a new relative to a clinically applied treatment planning system

BACKGROUND

The study objective was to validate the relative biological effectiveness (RBE) in RayStation for carbon-ion radiotherapy (CIRT) using the Syngo treatment planning system as reference.

METHODS

Local effect model I was established in RayStation (Ray-LEM) with the same parameters as in LEM I in Syngo (Syngo-LEM). Three cube plans covering most of the tumors treated at our center were generated with Syngo-LEM. Ray-LEM re-calculated the Syngo plans and compared the RBEs to the Syngo counterparts. The results showed that RayStation RBE was smaller than Syngo RBE. To ensure that Ray-LEM reproduced Syngo RBE, the observed deviations were used to scale the maximum RBE (RBEmax) in Ray-LEM. After this calibration, we further compared the RayStation RBE to Syngo RBE using additional plans in both homogeneous phantoms and patients, to ensure that the calibrated Ray-LEM reproduced Syngo RBE even with more complex planning features.

RESULTS

The calibration increased the RBEmax by 2.3% to raise the Ray-LEM RBE. The target mean RBE deviations in the phantom evaluation plans were median: 0.0 (minimum: - 1.1 to maximum: 0.7) %, and the target mean RBE deviations of the clinical target volumes of 16 patient cases were - 0.4 (- 1.5 to 0.2) %.

CONCLUSIONS

The residual RBE difference between RayStation and Syngo was found to be ≤ 1.0%. Thus, we can propose to use RayStation for clinical CIRT treatment planning. However, the potential differences due to the absorbed beam model warrants further exploration.

Clinical dose assessment for scanned carbon-ion radiotherapy using linear energy transfer measurements and Monte Carlo simulations

Objective. Dosimetric commissioning of treatment planning systems (TPS) focuses on validating the agreement of the physical dose with experimental data. For carbon-ion radiotherapy, the commissioning of the relative biological effectiveness (RBE) is necessary to predict the clinical outcome based on the radiation quality of the mixed radiation field. In this study, we proposed a approach for RBE commissioning using Monte Carlo (MC) simulations, which was further strengthen by RBE validation based on linear energy transfer (LET) measurements.Approach. First, we tuned the MC simulation based on the results of dosimetric experiments including the beam ranges, beam sizes, and MU calibrations. Furthermore, we compared simulated results to measured depth- and radial-LET distributions of the 430 MeV u^-1carbon-ion spot beam with a 1.5 mm², 36μm thick silicon detector. The measured dose-averaged LET (LET_d) and RBE were compared with the simulated results. The RBE was calculated based on the mixed beam model with linear-quadratic parameters depending on the LET. Finally, TPS-calculated clinical dose profiles were validated through the tuned MC-based calculations.Main results. A 10 keVμm^-1and 0.15 agreement for LET_dand RBE, respectively, were found between simulation and measurement results obtained for a 2σlateral size of 430 MeV u^-1carbon-ion spot beam in water. These results suggested that the tuned MC simulation can be used with acceptable precision for the RBE and LET calculations of carbon-ion spot beam within the clinical energy range. For physical and clinical doses, the TPS- and MC-based calculations showed good agreements within 1.0% at the centre of the spread-out Bragg peaks.Significance. The tuned MC simulation can accurately reproduce the actual carbon-ion beams, and it can be used to validate the physical and clinical dose distributions calculated by TPS. Moreover, the MC simulation can be used for dosimetric commissioning, including clinical doses, without LET measurements.

A case-control study of linear energy transfer and relative biological effectiveness related to symptomatic brainstem toxicity following pediatric proton therapy

BACKGROUND AND PURPOSE

A fixed relative biological effectiveness (RBE) of 1.1 (RBE_1.1) is used clinically in proton therapy even though the RBE varies with properties such as dose level and linear energy transfer (LET). We therefore investigated if symptomatic brainstem toxicity in pediatric brain tumor patients treated with proton therapy could be associated with a variable LET and RBE.

MATERIALS AND METHODS

36 patients treated with passive scattering proton therapy were selected for a case-control study from a cohort of 954 pediatric brain tumor patients. Nine children with symptomatic brainstem toxicity were each matched to three controls based on age, diagnosis, adjuvant therapy, and brainstem RBE_1.1 dose characteristics. Differences across cases and controls related to the dose-averaged LET (LET_d) and variable RBE-weighted dose from two RBE models were analyzed in the high-dose region.

RESULTS

LET_d metrics were marginally higher for cases vs. controls for the majority of dose levels and brainstem substructures. Considering areas with doses above 54 Gy(RBE_1.1), we found a moderate trend of 13% higher median LET_d in the brainstem for cases compared to controls (P = .08), while the difference in the median variable RBE-weighted dose for the same structure was only 2% (P = .6).

CONCLUSION

Trends towards higher LET_d for cases compared to controls were noticeable across structures and LET_d metrics for this patient cohort. While case-control differences were minor, an association with the observed symptomatic brainstem toxicity cannot be ruled out.

Normal tissue sparing by FLASH as a function of single fraction dose: A quantitative analysis

PURPOSE

The FLASH effect designates normal tissue sparing by ultra-high dose rate (UHDR) compared to conventional dose rate (CONV) irradiation without compromising tumor control. Understanding the magnitude of this effect and its dependency on dose are essential requirements for an optimized clinical translation of FLASH radiation therapy. In this context, we evaluated available experimental data on the magnitudes of normal tissue sparing provided by the FLASH effect as a function of dose, and followed a phenomenological data-driven approach for its parameterization.

METHODS

We gathered available in vivo data of the normal tissue sparing of CONV compared to UHDR single fraction doses and converted it to a common scale using isoeffect dose ratios, hereafter referred to as FLASH modifying factors (FMF). We then evaluated the suitability of a piecewise linear function with two pieces to parametrize FMF × D as a function of dose D.

RESULTS

We found that the magnitude of FMF generally decreases (i.e., sparing increases) as function of single fraction dose and that individual data series can be described by the piecewise linear function. The sparing magnitude appears organ specific. Pooled skin reaction data followed a consistent trend as a function of dose. Average FMF values and their standard deviations were 0.95±0.11 for all data below 10 Gy, 0.92±0.06 for mouse gut data between 10-25 Gy, and 0.96±0.07 and 0.71±0.06 for mammalian skin reaction data between 10-25 Gy and >25 Gy, respectively.

CONCLUSIONS

The magnitude of normal tissue sparing by FLASH is increasing with dose and is dependent on the irradiated tissue. A piecewise linear function can parameterize currently available individual data series.

Spot-scanning hadron arc (SHArc) therapy: A proof of concept using single and multi-ion strategies with helium, carbon, oxygen and neon ions

PURPOSE

To present particle arc therapy treatments using single and multi-ion therapy optimization strategies with helium (⁴ He), carbon (¹² C), oxygen (¹⁶ O) and neon (²⁰ Ne) ion beams.

METHODS AND MATERIALS

An optimization procedure and workflow were devised for spot-scanning hadron arc therapy (SHArc) treatment planning in the PRECISE (PaRticle thErapy using single and Combined Ion optimization StratEgies) treatment planning system (TPS). Physical and biological beam models were developed for helium, carbon, oxygen and neon ions via FLUKA MC simulation. SHArc treatments were optimized using both single ion (¹² C, ¹⁶ O, or ²⁰ Ne) and multi-ion therapy (¹⁶ O+⁴ He or ²⁰ Ne+⁴ He) applying variable relative biological effectiveness (RBE) modeling using a modified microdosimetric kinetic model (mMKM) with (α/β)_x values of 2Gy, 5Gy and 3.1Gy respectively, for glioblastoma, pancreatic adenocarcinoma, and prostate adenocarcinoma patient cases. Dose, effective dose, linear energy transfer (LET) and RBE were computed with the GPU-accelerated dose engine FRoG and dosimetric/biophysical attributes were evaluated in the context of conventional particle and photon-based therapies (e.g., volumetric modulated arc therapy [VMAT]).

RESULTS

All SHArc plans met the target optimization goals (3GyRBE) and demonstrated increased target conformity and substantially lower low-dose bath to surrounding normal tissues than VMAT. SHArc plans using a single ion species (¹² C, ¹⁶ O, or ²⁰ Ne) exhibited favorable LET distributions with the highest-LET components centralized in the target volume, with values ranging from ∼80-170keV/μm, ∼130-220keV/μm and ∼180-350keV/μm, for ¹² C, ¹⁶ O, or ²⁰ Ne, respectively, exceeding mean target LET of conventional particle therapy (¹² C:∼60, ¹⁶ O:∼78 ²⁰ Ne:∼100 keV/μm). Multi-ion therapy with SHArc delivery (SHArc_MIT ) provided a similar level of target LET enhancement as SHArc compared to conventional planning, however, with additional benefits of homogenous physical dose and RBE distributions.

CONCLUSION

Here, we demonstrate that arc delivery of light and heavy ion beams, using either a single ion species (¹² C, ¹⁶ O, or ²⁰ Ne) or combining two ions in a single fraction (¹⁶ O+⁴ He or ²⁰ Ne+⁴ He), affords enhanced physical and biological distributions (e.g., LET) compared with conventional delivery with photons or particle beams. SHArc marks the first single and multi-ion arc therapy treatment optimization approach using light and heavy ions. This article is protected by copyright. All rights reserved.

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.

Technical note: Impact of beamline-specific particle energy spectra on clinical plans in carbon ion beam therapy

PURPOSE

The Local Effect Model version one (LEM I) is applied clinically across Europe to quantify the relative biological effectiveness (RBE) of carbon ion beams. It requires the full particle fluence spectrum differential in energy in each voxel as input parameter. Treatment planning systems (TPSs) use beamline-specific look-up tables generated with Monte Carlo (MC) codes. In this study, the changes in RBE weighted dose were quantified using different levels of details in the simulation or different MC codes.

METHODS

The particle fluence differential in energy was simulated with FLUKA and Geant4 at 500 depths in water in 1-mm steps for 58 initial carbon ion energies (between 120.0 and 402.8 MeV/u). A dedicated beam model was applied, including the full description of the Nozzle using GATE-RTionV1.0 (Geant4.10.03p03). In addition, two tables generated with FLUKA were compared. The starting points of the FLUKA simulations were phase space (PhS) files from, firstly, the Geant4 nozzle simulations, and secondly, a clinical beam model where an analytic approach was used to mimic the beamline. Treatment plans (TPs) were generated with RayStation 8B (RaySearch Laboratories AB, Sweden) for cubic targets in water and 10 clinical patient cases using the clinical beam model. Subsequently, the RBE weighted dose was re-computed using the two other fluence tables (FLUKA PhS or Geant4).

RESULTS

The fluence spectra of the primary and secondary particles simulated with Geant4 and FLUKA generally agreed well for the primary particles. Differences were mainly observed for the secondary particles. Interchanging the two energy spectra (FLUKA vs. GEANT4) to calculate the RBE weighted dose distributions resulted in average deviations of less than 1% in the entrance up to the end of the target region, with a maximum local deviation at the distal edge of the target. In the fragment tail, larger discrepancies of up to 5% on average were found for deep-seated targets. The patient and water phantom cases demonstrated similar results.

CONCLUSION

RBE weighted doses agreed well within all tested setups, confirming the clinical beam model provided by the TPS vendor. Furthermore, the results showed that the open source and generally available MC code Geant4 (in particular using GATE or GATE-RTion) can also be used to generate basic beam data required for RBE calculation in carbon ion therapy.

Reformulated McNamara RBE-weighted beam orientation optimization for intensity modulated proton therapy

PURPOSE

Empirical relative biological effectiveness (RBE) models have been used to estimate the biological dose in proton therapy but do not adequately capture the factors influencing RBE values for treatment planning. We reformulate the McNamara RBE model such that it can be added as a linear biological dose fidelity term within our previously developed sensitivity-regularized and heterogeneity-weighted beam orientation optimization (SHBOO) framework.

METHODS

Based on our SHBOO framework, we formulated the biological optimization problem to minimize total McNamara RBE dose to OARs. We solve this problem using two optimization algorithms: FISTA (McNam-FISTA) and Chambolle-Pock (McNam-CP). We compare their performances with a physical dose optimizer assuming RBE = 1.1 in all structures (PHYS-FISTA) and an LET-weighted dose model (LET-FISTA). Three head and neck patients were planned with the four techniques and compared on dosimetry and robustness.

RESULTS

Compared to Phys-FISTA, McNam-CP was able to match CTV [HI, Dmax, D95%, D98%] by [0.00, 0.05%, 1.4%, 0.8%]. McNam-FISTA and McNam-CP were able to significantly improve overall OAR [Dmean, Dmax] by an average of [36.1%,26.4%] and [29.6%, 20.3%], respectively. Regarding CTV robustness, worst [Dmax, V95%, D95%, D98%] improvement of [-6.6%, 6.2%, 6.0%, 4.8%] was reported for McNam-FISTA and [2.7%, 2.7%, 5.3%, -4.3%] for McNam-CP under combinations of range and setup uncertainties. For OARs, worst [Dmax, Dmean] were improved by McNam-FISTA and McNam-CP by an average of [25.0%, 19.2%] and [29.5%, 36.5%], respectively. McNam-FISTA considerably improved dosimetry and CTV robustness compared to LET-FISTA, which achieved better worst-case OAR doses.

CONCLUSION

The four optimization techniques deliver comparable biological doses for the head and neck cases. Besides modest CTV coverage and robustness improvement, OAR biological dose and robustness were substantially improved with both McNam-FISTA and McNam-CP, showing potential benefit for directly incorporating McNamara RBE in proton treatment planning.

Dose- and LET-dependent changes in mouse skin contracture up to a year after either single dose or fractionated doses of carbon ion or gamma rays

Time dependence of relative biological effectiveness (RBE) of carbon ions for skin damage was investigated to answer the question of whether the flat distribution of biological doses within a Spread-Out Bragg peak (SOBP) which is designed based on in vitro cell kill could also be flat for in vivo late responding tissue. Two spots of Indian ink intracutaneously injected into the legs of C3H mice were measured by calipers. An equieffective dose to produce 30% skin contraction was calculated from a dose-response curve and used to calculate the RBE of carbon ion beams. We discovered skin contraction progressed after irradiation and then reached a stable/slow progression phase. Equieffective doses decreased with time and the decrease was most prominent for gamma rays and least prominent for 100 keV/μm carbon ions. Survival parameter of alpha but not beta in the linear-quadratic model is closely related to the RBE of carbon ions. Biological doses within the SOBP increased with time but their distribution was still flat up to 1 year after irradiation. The outcomes of skin contraction studies suggest that (i) despite the higher RBE for skin contracture after carbon ions compared to gamma rays, gamma rays can result in a more severe late effect of skin contracture. This is due to the carbon effect saturating at a lower dose than gamma rays, and (ii) the biological dose distribution throughout the SOBP remains approximately the same even one year after exposure.

Alpha-Emitter Radiopharmaceuticals and External Beam Radiotherapy: A Radiobiological Model for the Combined Treatment

Previously published studies combined external beam radiotherapy (EBRT) treatments with different activities of ²²³Ra. The data of two-year overall survival (2y-OS) and neutropenia (TOX) incidence when combining EBRT and ²²³Ra are not homogeneous in literature. We adapted the linear-quadratic model (LQ) to ²²³Ra therapy using brachytherapy formalism for a mixture of radionuclides, considering the contribution of all daughter isotopes in the decay chain. A virtual cohort of patients undergoing ²²³Ra therapy was derived using data from the literature. The doses delivered using ²²³Ra and EBRT were converted into biologically equivalent doses. Fixed-effect logistic regression models were derived for both the 2y-OS and TOX and compared with available literature. Based on the literature search, four studies were identified to have reported the ²²³Ra injection activity levels varying from the placebo (0) to 80 kBq/kg, associated or not with EBRT. Logistic regression models revealed a dose-dependent increase in both the 2y-OS (intercept = -1.364; slope = 0.006; p-value ≤ 0.05) and TOX (-5.035; 0.018; ≤0.05) using the EBRT schedule of 8 Gy in 1 fr. Similar results were obtained for other schedules. Discrepancies between our TOX model and those derived for EBRT combined with chemotherapy are discussed. Radiobiological models allow us to estimate dose-dependent relationships, to predict the OS and TOX following combined ²²³Ra + EBRT treatment, which will guide future treatment optimization.

High-Accuracy Relative Biological Effectiveness Values Following Low-Dose Thermal Neutron Exposures Support Bimodal Quality Factor Response with Neutron Energy

Theoretical evaluations indicate the radiation weighting factor for thermal neutrons differs from the current International Commission on Radiological Protection (ICRP) recommended value of 2.5, which has radiation protection implications for high-energy radiotherapy, inside spacecraft, on the lunar or Martian surface, and in nuclear reactor workplaces. We examined the relative biological effectiveness (RBE) of DNA damage generated by thermal neutrons compared to gamma radiation. Whole blood was irradiated by 64 meV thermal neutrons from the National Research Universal reactor. DNA damage and erroneous DNA double-strand break repair was evaluated by dicentric chromosome assay (DCA) and cytokinesis-block micronucleus (CBMN) assay with low doses ranging 6–85 mGy. Linear dose responses were observed. Significant DNA aberration clustering was found indicative of high ionizing density radiation. When the dose contribution of both the ¹⁴N(n,p)¹⁴C and ¹H(n,γ)²H capture reactions were considered, the DCA and the CBMN assays generated similar maximum RBE values of 11.3 ± 1.6 and 9.0 ± 1.1, respectively. Consequently, thermal neutron RBE is approximately four times higher than the current ICRP radiation weighting factor value of 2.5. This lends support to bimodal peaks in the quality factor for RBE neutron energy response, underlining the importance of radiological protection against thermal neutron exposures.

212Pb-conjugated anti-rat HER2/neu antibody against a neu-N derived murine mammary carcinoma cell line: cell kill and RBE in vitro

PURPOSE

In the current work, the RBE of a ²¹²Pb-conjugated anti-HER2/neu antibody construct has been evaluated, in vitro, by colony formation assay. The RBE was estimated by comparing two absorbed dose-survival curves: the first obtained from the conjugated ²¹²Pb experiments (test radiation), the second obtained by parallel experiments of single bolus irradiation of external beam (reference radiation).

MATERIALS AND METHODS

Mammary carcinoma NT2.5 cells were treated with (0-3.70) kBq/ml of radiolabeled antibody. Nonspecific binding was assessed with addition of excess amount of unlabeled antibody. The colony formation curves were converted from activity concentration to cell nucleus absorbed dose by simulating the decay and transport of all daughter and secondary particles of ²¹²Pb, using the Monte Carlo code GEANT 4.

RESULTS

The radiolabeled antibody yielded an RBE of 8.3 at 37% survival and a survival independent RBE (i.e. RBE2) of 9.9. Unbound/untargeted ²¹²Pb-labeled antibody, as obtained in blocking experiments yielded minimal alpha-particle radiation to cells. Conclusions: These results further highlight the importance of specific targeting toward achieving tumor cell kill and low toxicity to normal tissue.

Microdosimetric characterization of a clinical proton therapy beam: comparison between simulated lineal energy distributions in spherical water targets and experimental measurements with a silicon detector

Objective Treatment planning based on computer simulations were proposed to account for the increase in the relative biological effectiveness (RBE) of proton radiotherapy beams near to the edges of the irradiated volume. Since silicon detectors could be used to validate the results of these simulations, it is important to explore the limitations of this comparison. Approach Microdosimetric measurements with a MicroPlus Bridge V2 silicon detector (thickness = 10 µm) were performed along the Bragg peak of a clinical proton beam. The lineal energy distributions, the dose mean values, and the RBE calculated with a biological weighting function were compared with simulations with PHITS (microdosimetric target = 1 µm water sphere), and published clonogenic survival in vitro RBE data for the V79 cell line. The effect of the silicon-to-water conversion was also investigated by comparing three different methodologies (conversion based on a single value, novel bin-to-bin conversions based on SRIM and PSTAR). Main results Mainly due to differences in the microdosimetric targets, the experimental dose-mean lineal energy and RBE values at the distal edge were respectively up to 53% and 28% lower than the simulated ones. Furthermore, the methodology chosen for the silicon-to-water conversion was proven to affect the dose mean lineal energy and the RBE10 up to 32% and 11% respectively. The best methodology to compensate for this underestimation was the bin-to-bin silicon-to-water conversion based on PSTAR. Significance This work represents the first comparison between PHITS-simulated lineal energy distributions in water targets and corresponding experimental spectra measured with silicon detectors. Furthermore, the effect of the silicon-to-water conversion on the RBE was explored for the first time. The proposed methodology based on the PSTAR bin-to-bin conversion appears to provide superior results with respect to commonly used single scaling factors and is recommended for future studies.

Impact of a spatially dependent dose delivery time structure on the biological effectiveness of scanning proton therapy

PURPOSE

In the scanning beam delivery of protons, different portions of the target are irradiated with different linear energy transfer protons with various time intervals and irradiation times. This research aimed to evaluate the spatially dependent biological effectiveness of protracted irradiation in scanning proton therapy.

METHODS

One and two parallel opposed fields plans were created in water phantom with the prescribed dose of 2 Gy. Three scenarios (instantaneous, continuous, and layered scans) were used with the corresponding beam delivery models. The biological dose (physical dose × relative biological effectiveness) was calculated using the linear quadratic model and the theory of dual radiation action to quantitatively evaluate the dose delivery time effect. In addition, simulations using clinical plans (postoperative seminoma and prostate tumor cases) were conducted to assess the impact of the effects on the dose volume histogram parameters and homogeneity coefficient (HC) in targets.

RESULTS

In a single-field plan of water phantom, when the treatment time was 19 min, the layered-scan scenario showed a decrease of <0.2% (almost 3.3%) in the biological dose from the plan on the distal (proximal) side because of the high (low) dose rate. This is in contrast to the continuous scenario, where the biological dose was almost uniformly decreased over the target by approximately 3.3%. The simulation with clinical geometry showed that the decrease rates in D_99% were 0.9% and 1.5% for every 10 min of treatment time prolongation for postoperative seminoma and prostate tumor cases, respectively, whereas the increase rates in HC were 0.7% and 0.2%.

CONCLUSIONS

In protracted irradiation in scanning proton therapy, the spatially dependent dose delivery time structure in scanning beam delivery can be an important factor for accurate evaluation of biological effectiveness.

Comprehensive comparison of local effect model IV predictions with the particle irradiation data ensemble

PURPOSE

The increased relative biological effectiveness (RBE) of ions is one of the key benefits of ion radiotherapy compared to conventional radiotherapy with photons. To account for the increased RBE of ions during the process of ion radiotherapy treatment planning, a robust model for RBE predictions is indispensable. Currently, at several ion therapy centers the local effect model I (LEM I) is applied to predict the RBE, which varies with biological and physical impacting factors. After the introduction of LEM I, several model improvements were implemented, leading to the current version, LEM IV, which is systematically tested in this study.

METHODS

As a comprehensive RBE model should give consistent results for a large variety of ion species and energies, the particle irradiation data ensemble (PIDE) is used to systematically validate the LEM IV. The database covers over 1100 photon and ion survival experiments in form of their linear-quadratic parameters for a wide range of ion types and energies. This makes the database an optimal tool to challenge the systematic dependencies of the RBE model. After appropriate filtering of the database, 571 experiments were identified and used as test data.

RESULTS

The study confirms that the LEM IV reflects the RBE systematics observed in measurements well. It is able to reproduce the dependence of RBE on the linear energy transfer (LET) as well as on the α_γ /β_γ ratio for several ion species in a wide energy range. Additionally, the systematic quantitative analysis revealed precision capabilities and limits of the model. At lower LET values, the LEM IV tends to underestimate the RBE with an increasing underestimation with increasing atomic number of the ion. At higher LET values, the LEM IV overestimates the RBE for protons or helium ions, whereas the predictions for heavier ions match experimental data well.

CONCLUSIONS

The LEM IV is able to predict general RBE characteristics for several ion species in a broad energy range. The accuracy of the predictions is reasonable considering the small number of input parameters needed by the model. The detailed quantification of possible systematic deviations, however, enables to identify not only strengths but also limitations of the model. The gained knowledge can be used to develop model adjustments to further improve the model accuracy, which is on the way.

Comparison between the results of a recently-developed biological weighting function (V79-RBE10BWF) and thein vitroclonogenic survival RBE10of other repair-competent asynchronized normoxic mammalian cell lines and ions not used for the development of the model

728 simulated microdosimetric lineal energy spectra (26 different ions between 1H and 238U, 28 energy points from 1 to 1000 MeV/n) were used in combination with a recently-developed biological weighting function (Parisi et al., 2020) and 571 published in vitro clonogenic survival curves in order to: 1) assess prediction intervals for the in silico results by deriving an empirical indication of the experimental uncertainty from the dispersion in the in vitro hamster lung fibroblast (V79) data used for the development of the biophysical model; 2) explore the possibility of modeling the relative biological effectiveness (RBE) of the 10% clonogenic survival of asynchronized normoxic repair-competent mammalian cell lines other than the one used for the development of the model (V79); 3) investigate the predictive power of the model through a comparison between in silico results and in vitro data for 10 ions not used for the development of the model. At first, different strategies for the assessment of the in silico prediction intervals were compared. The possible sources of uncertainty responsible for the dispersion in the in vitro data were also shortly reviewed. Secondly, also because of the relevant scatter in the in vitro data, no statistically-relevant differences were found between the RBE10 of the investigated different asynchronized normoxic repair-competent mammalian cell lines. The only exception (Chinese Hamster peritoneal fibroblasts, B14FAF28), is likely due to the limited dataset (all in vitro ion data were extracted from a single publication), systematic differences in the linear energy transfer (LET) calculations for the employed very-heavy ions, and the use of reference photon survival curves extracted from a different publication. Finally, the in silico predictions for the 10 ions not used for the model development were in good agreement with the corresponding in vitro data.

Response of Normal Tissues to Boron Neutron Capture Therapy (BNCT) with 10B-Borocaptate Sodium (BSH) and 10B-Paraboronophenylalanine (BPA)

Boron neutron capture therapy (BNCT) is a cancer-selective radiotherapy that utilizes the cancer targeting ¹⁰B-compound. Cancer cells that take up the compound are substantially damaged by the high liner energy transfer (LET) particles emitted mainly from the ¹⁰B(n, α⁷Li reaction. BNCT can minimize the dose to normal tissues, but it must be performed within the tolerable range of normal tissues. Therefore, it is important to evaluate the response of normal tissues to BNCT. Since BNCT yields a mixture of high and low LET radiations that make it difficult to understand the radiobiological basis of BNCT, it is important to evaluate the relative biological effectiveness (RBE) and compound biological effectiveness (CBE) factors for assessing the responses of normal tissues to BNCT. BSH and BPA are the only ¹⁰B-compounds that can be used for clinical BNCT. Their biological behavior and cancer targeting mechanisms are different; therefore, they affect the CBE values differently. In this review, we present the RBE and CBE values of BPA or BSH for normal tissue damage by BNCT irradiation. The skin, brain (spinal cord), mucosa, lung, and liver are included as normal tissues. The CBE values of BPA and BSH for tumor control are also discussed.

Adaptation of stochastic microdosimetric kinetic model to hypoxia for hypo-fractionated multi-ion therapy treatment planning

For hypo-fractionated multi-ion therapy (HFMIT), the stochastic microdosimetric kinetic (SMK) model had been developed to estimate the biological effectiveness of radiation beams with wide linear energy transfer (LET) and dose ranges. The HFMIT will be applied to radioresistant tumors with oxygen-deficient regions. The response of cells to radiation is strongly dependent on the oxygen condition in addition to radiation type, LET and absorbed dose. This study presents an adaptation of the SMK model to account for oxygen-pressure dependent cell responses, and develops the oxygen-effect-incorporated stochastic microdosimetric kinetic (OSMK) model. In the model, following assumptions were made: the numbers of radiation-induced sublethal lesions (double-strand breaks) are reduced due to lack of oxygen, and the numbers of oxygen-mediated lesions are reduced for radiation with high LET. The model parameters were determined by fitting survival data under aerobic and anoxic conditions for human salivary gland tumor cells and V79 cells exposed to helium-, carbon-, and neon-ion beams over the LET range of 18.5-654.0 keVμm^-1. The OSMK model provided good agreement with the experimental survival data of the cells with determination coefficients >0.9. In terms of oxygen enhancement ratio, the OSMK model reproduced the experimental data behavior, including slight dependence on particle type at the same LET. The OSMK model was then implemented into the in-house treatment planning software for the HFMIT to validate its applicability in clinical practice. A treatment plan with helium- and neon-ion beams was made for a pancreatic cancer case assuming an oxygen-deficient region within the tumor. The biological optimization based on the OSMK model preferentially placed the neon-ion beam to the hypoxic region, while it placed both helium- and neon-ion beams to the surrounding normoxic region. The OSMK model offered the accuracy and usability required for hypoxia-based biological optimization in HFMIT treatment planning.

Handling of beam spectra in training and application of proton RBE models

Published data from cell survival experiments are frequently used as training data for models of proton relative biological effectiveness (RBE). The publications rarely provide full information about the primary particle spectrum of the used beam, or its content of heavy secondary particles. The purpose of this paper is to assess to what extent heavy secondary particles may have been present in published cell survival experiments, and to investigate the impact of non-primary protons for RBE calculations in treatment planning. We used the Monte Carlo code Geant4 to calculate the occurrence of non-primary protons and heavier secondary particles for clinical protons beams in water for four incident energies in the [100, 250] MeV interval. We used the resulting spectra together with a conservative RBE parameterization and an RBE model to map both the rise of RBE at the beam entry surface due to heavy secondary particle buildup, and the difference in estimated RBE if non-primary protons are included or not in the beam quality metric. If included, non-primary protons cause a difference of 2% of the RBE in the plateau region of an spread out Bragg peak and 1% in the Bragg peak. Including non-primary protons specifically for RBE calculations will consequently have a negligible impact and can be ignored. A buildup distance in water of one millimeter was sufficient to reach an equilibrium state of RBE for the four incident energies selected. For the investigated experimental data, 83 out of the 86 data points were found to have been determined with at least that amount of buildup material. Hence, RBE model training data should be interpreted to include the contribution of heavy secondaries.

Can the mean linear energy transfer of organs be directly related to patient toxicities for current head and neck cancer intensity-modulated proton therapy practice?

BACKGROUND AND PURPOSE

The relative biological effectiveness (RBE) of proton therapy is predicted to vary with the dose-weighted average linear energy transfer (LET_d). However, RBE values may substantially vary for different clinical endpoints. Therefore, the aim of this study was to assess the feasibility of relating mean D⋅LET_d parameters to patient toxicity for HNC patients treated with proton therapy.

MATERIALS AND METHODS

The delivered physical dose (D) and the voxel-wise product of D and LET_d (D⋅LET_d) distributions were calculated for 100 head and neck cancer (HNC) proton therapy patients using our TPS (Raystation v6R). The means and covariance matrix of the accumulated D and D⋅LET_d of all relevant organs-at-risk (OARs) were used to simulate 2.500 data sets of different sizes. For each dataset, an attempt was made to add mean D⋅LET_d parameters to a multivariable NTCP model based on mean D parameters of the same OAR for xerostomia, tube feeding and dysphagia. The likelihood of creating an NTCP model with statistically significant parameters (i.e. power) was calculated as a function of the simulated sample size for various RBE models.

RESULTS

The sample size required to have a power of at least 80% to show an independent effect of mean D⋅LET_d parameters on toxicity is over 15,000 patients for all toxicities.

CONCLUSION

For current clinical practice, it is not feasible to directly model NTCP with both mean D and mean D⋅LET_d of OARs. These findings should not be interpreted as a contradiction of previous evidence for the relationship between RBE and LET_d.

Does the uncertainty in relative biological effectiveness affect patient treatment in proton therapy?

Clinical treatment with protons uses the concept of relative biological effectiveness (RBE) to convert the absorbed dose into an RBE-weighted dose that equals the dose for radiotherapy with photons causing the same biological effect. Currently, in proton therapy a constant RBE of 1.1 is generically used. However, empirical data indicate that the RBE is not constant, but increases at the distal edge of the proton beam. This increase in RBE is of concern, as the clinical impact is still unresolved, and clinical studies demonstrating a clinical effect of an increased RBE are emerging. Within the European Particle Therapy Network (EPTN) work package 6 on radiobiology and RBE, a workshop was held in February 2020 in Manchester with one day of discussion dedicated to the impact of proton RBE in a clinical context. Current data on RBE effects, patient outcome and modelling from experimental as well as clinical studies were presented and discussed. Furthermore, representatives from European clinical proton therapy centres, who were involved in patient treatment, laid out their current clinical practice on how to consider the risk of a variable RBE in their centres. In line with the workshop, this work considers the actual impact of RBE issues on patient care in proton therapy by reviewing preclinical data on the relation between linear energy transfer (LET) and RBE, current clinical data sets on RBE effects in patients, and applied clinical strategies to manage RBE uncertainties. A better understanding of the variability in RBE would allow development of proton treatments which are safer and more effective.