The use of a constant RBE=1.1 for proton radiotherapy is no longer appropriate

Monte Carlo simulations of cell survival in proton SOBP

Objective. The objective of this study is to develop a multi-scale modeling approach that accurately predicts radiation-induced DNA damage and survival fraction in specific cell lines.Approach. A Monte Carlo based simulation framework was employed to make the predictions. The FLUKA Monte Carlo code was utilized to estimate absorbed doses and fluence energy spectra, which were then used in the Monte Carlo Damage Simulation code to compute DNA damage yields in Chinese hamster V79 cell lines. The outputs were converted into cell survival fractions using a previously published theoretical model. To reduce the uncertainties of the predictions, new values for the parameters of the theoretical model were computed, expanding the database of experimental points considered in the previous estimation. Simulated results were validated against experimental data, confirming the applicability of the framework for proton beams up to 230 MeV. Additionally, the impact of secondary particles on cell survival was estimated.Main results. The simulated survival fraction versus depth in a glycerol phantom is reported for eighteen different configurations. Two proton spread out Bragg peaks at several doses were simulated and compared with experimental data. In all cases, the simulations follow the experimental trends, demonstrating the accuracy of the predictions up to 230 MeV.Significance. This study holds significant importance as it contributes to the advancement of models for predicting biological responses to radiation, ultimately contributing to more effective cancer treatment in proton therapy.

Linear Energy Transfer and Relative Biological Effectiveness Investigation of Various Structures for a Cohort of Proton Patients With Brain Tumors

Purpose

The current knowledge on biological effects associated with proton therapy is limited. Therefore, we investigated the distributions of dose, dose-averaged linear energy transfer (LET_d), and the product between dose and LET_d (DLET_d) for a patient cohort treated with proton therapy. Different treatment planning system features and visualization tools were explored.

Methods and Materials

For a cohort of 24 patients with brain tumors, the LET_d, DLET_d, and dose was calculated for a fixed relative biological effectiveness value and 2 variable models: plan-based and phenomenological. Dose threshold levels of 0, 5, and 20 Gy were imposed for LET_d visualization. The relationship between physical dose and LET_d and the frequency of LET_d hotspots were investigated.

Results

The phenomenological relative biological effectiveness model presented consistently higher dose values. For lower dose thresholds, the LET_d distribution was steered toward higher values related to low treatment doses. Differences up to 26.0% were found according to the threshold. Maximum LET_d values were identified in the brain, periventricular space, and ventricles. An inverse relationship between LET_d and dose was observed. Frequency information to the domain of dose and LET_d allowed for the identification of clusters, which steer the mean LET_d values, and the identification of higher, but sparse, LET_d values.

Conclusions

Identifying, quantifying, and recording LET distributions in a standardized fashion is necessary, because concern exists over a link between toxicity and LET hotspots. Visualizing DLET_d or dose × LET_d during treatment planning could allow for clinicians to make informed decisions.

Application of a simple DNA damage model developed for electrons to proton irradiation

Proton beam therapy allows irradiating tumor volumes with reduced side effects on normal tissues with respect to conventional x-ray radiotherapy. Biological effects such as cell killing after proton beam irradiations depend on the proton kinetic energy, which is intrinsically related to early DNA damage induction. As such, DNA damage estimation based on Monte Carlo simulations is a research topic of worldwide interest. Such simulation is a mean of investigating the mechanisms of DNA strand break formations. However, past modellings considering chemical processes and DNA structures require long calculation times. Particle and heavy ion transport system (PHITS) is one of the general-purpose Monte Carlo codes that can simulate track structure of protons, meanwhile cannot handle radical dynamics simulation in liquid water. It also includes a simple model enabling the efficient estimation of DNA damage yields only from the spatial distribution of ionizations and excitations without DNA geometry, which was originally developed for electron track-structure simulations. In this study, we investigated the potential application of the model to protons without any modification. The yields of single-strand breaks, double-strand breaks (DSBs) and the complex DSBs were assessed as functions of the proton kinetic energy. The PHITS-based estimation showed that the DSB yields increased as the linear energy transfer (LET) increased, and reproduced the experimental and simulated yields of various DNA damage types induced by protons with LET up to about 30 keV μm−1. These results suggest that the current DNA damage model implemented in PHITS is sufficient for estimating DNA lesion yields induced after protons irradiation except at very low energies (below 1 MeV). This model contributes to evaluating early biological impacts in radiation therapy.

Treatment-planning approaches to intensity modulated proton therapy and the impact on dose-weighted linear energy transfer

PURPOSE

We quantified the effect of various forward-based treatment-planning strategies in proton therapy on dose-weighted linear energy transfer (LETd). By maintaining the dosimetric quality at a clinically acceptable level, we aimed to evaluate the differences in LETd among various treatment-planning approaches and their practicality in minimizing biologic uncertainties associated with LETd.

METHOD

Eight treatment-planning strategies that are achievable in commercial treatment-planning systems were applied on a cylindrical water phantom and four pediatric brain tumor cases. Each planning strategy was compared to either an opposed lateral plan (phantom study) or original clinical plan (patient study). Deviations in mean and maximum LETd from clinically acceptable dose distributions were compared.

RESULTS

In the phantom study, using a range shifter and altering the robust scenarios during optimization had the largest effect on the mean clinical target volume LETd, which was reduced from 4.5 to 3.9 keV/μm in both cases. Variations in the intersection angle between beams had the largest effect on LETd in a ring defined 3 to 5 mm outside the target. When beam intersection angles were reduced from opposed laterals (180°) to 120°, 90°, and 60°, corresponding maximum LETd increased from 7.9 to 8.9, 10.9, and 12.2 keV/μm, respectively. A clear trend in mean and maximum LETd variations in the clinical cases could not be established, though spatial distribution of LETd suggested a strong dependence on patient anatomy and treatment geometry.

CONCLUSION

Changes in LETd from treatment-plan setup follow intuitive trends in a controlled phantom experiment. Anatomical and other patient-specific considerations, however, can preclude generalizable strategies in clinical cases. For pediatric cranial radiation therapy, we recommend using opposed lateral treatment fields to treat midline targets.

Study of relationship between dose, LET and the risk of brain necrosis after proton therapy for skull base tumors

PURPOSE

We investigated the relationship between RBE-weighted dose (DRBE) calculated with constant (cRBE) and variable RBE (vRBE), dose-averaged linear energy transfer (LETd) and the risk of radiographic changes in skull base patients treated with protons.

METHODS

Clinical treatment plans of 45 patients were recalculated with Monte Carlo tool FRED. Radiographic changes (i.e. edema and/or necrosis) were identified by MRI. Dosimetric parameters for cRBE and vRBE were computed. Biological margin extension and voxel-based analysis were employed looking for association of DRBE(vRBE) and LETd with brain edema and/or necrosis.

RESULTS

When using vRBE, Dmax in the brain was above the highest dose limits for 38% of patients, while such limit was never exceeded assuming cRBE. Similar values of Dmax were observed in necrotic regions, brain and temporal lobes. Most of the brain necrosis was in proximity to the PTV. The voxel-based analysis did not show evidence of an association with high LETd values.

CONCLUSIONS

When looking at standard dosimetric parameters, the higher dose associated with vRBE seems to be responsible for an enhanced risk of radiographic changes. However, as revealed by a voxel-based analysis, the large inter-patient variability hinders the identification of a clear effect for high LETd.

Clinical Progress in Proton Radiotherapy: Biological Unknowns

Simple Summary

Proton radiation therapy is a more recent type of radiotherapy that uses proton beams instead of classical photon or X-rays beams. The clinical benefit of proton therapy is that it allows to treat tumors more precisely. As a result, proton radiotherapy induces less toxicity to healthy tissue near the tumor site. Despite the experience in the clinical use of protons, the response of cells to proton radiation, the radiobiology, is less understood. In this review, we describe the current knowledge about proton radiobiology.

Abstract

Clinical use of proton radiation has massively increased over the past years. The main reason for this is the beneficial depth-dose distribution of protons that allows to reduce toxicity to normal tissues surrounding the tumor. Despite the experience in the clinical use of protons, the radiobiology after proton irradiation compared to photon irradiation remains to be completely elucidated. Proton radiation may lead to differential damages and activation of biological processes. Here, we will review the current knowledge of proton radiobiology in terms of induction of reactive oxygen species, hypoxia, DNA damage response, as well as cell death after proton irradiation and radioresistance.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Early clinical results of proton spatially fractionated GRID radiation therapy (SFGRT)

OBJECTIVE

Approximately 70 patients with large and bulky tumors refractory to prior treatments were treated with photon spatially fractionated GRID radiation (SFGRT). We identified 10 additional patients who clinically needed GRID but could not be treated with photons due to adjacent critical organs. We developed a proton SFGRT technique, and we report treatment of these 10 patients.

METHODS

Subject data were reviewed for clinical results and dosimetric data. 50% of the patients were metastatic at the time of treatment and five had previous photon radiation to the local site but not via GRID. They were treated with 15-20 cobalt Gray equivalent using a single proton GRID field with an average beamlet count of 22.6 (range 7-51). 80% received an average adjuvant radiation dose to the GRID region of 40.8Gy (range 13.7-63.8Gy). Four received subsequent systemic therapy.

RESULTS

The median follow-up time was 5.9 months (1.1-18.9). At last follow-up, seven patients were alive and three had died. Two patients who had died from metastatic disease had local shrinkage of tumor. Of those alive, four had complete or partial response, two had partial response but later progressed, and one had no response. For all patients, the tumor regression/local symptom improvement rate was 80%. 50% had acute side-effects of grade1/2 only and all were well-tolerated.

CONCLUSION

In circumstances where patients cannot receive photon GRID, proton SFGRT is clinically feasible and effective, with a similar side-effect profile.

ADVANCES IN KNOWLEDGE

Proton GRID should be considered as a treatment option earlier in the disease course for patients who cannot be treated by photon GRID.

Proton RBE dependence on dose in the setting of hypofractionation

Hypofractionated radiotherapy is attractive concerning patient burden and therapy costs, but many aspects play a role when it comes to assess its safety. While exploited for conventional photon therapy and carbon ion therapy, hypofractionation with protons is only rarely applied. One reason for this is uncertainty in the described dose, mainly due to the relative biological effectiveness (RBE), which is small for protons, but not negligible. RBE is generally dose-dependent, and for higher doses as used in hypofractionation, a thorough RBE evaluation is needed. This review article focuses on the RBE variability in protons and associated issues or implications for hypofractionation.

Introducing Proton Track-End Objectives in Intensity Modulated Proton Therapy Optimization to Reduce Linear Energy Transfer and Relative Biological Effectiveness in Critical Structures

PURPOSE

We propose the use of proton track-end objectives in intensity modulated proton therapy (IMPT) optimization to reduce the linear energy transfer (LET) and the relative biological effectiveness (RBE) in critical structures.

METHODS AND MATERIALS

IMPT plans were generated for 3 intracranial patient cases (1.8 Gy (RBE) in 30 fractions) and 3 head-and-neck patient cases (2 Gy (RBE) in 35 fractions), assuming a constant RBE of 1.1. Two plans were generated for each patient: (1) physical dose objectives only (DOSEopt) and (2) same dose objectives as the DOSEopt plan, with additional proton track-end objectives (TEopt). The track-end objectives penalized protons stopping in the risk volume of choice. Dose evaluations were made using a RBE of 1.1 and the LET-dependent Wedenberg RBE model, together with estimates of normal tissue complication probabilities (NTCPs). In addition, the distributions of proton track-ends and dose-average LET (LET_d) were analyzed.

RESULTS

The TEopt plans reduced the mean LET_d in the critical structures studied by an average of 37% and increased the mean LET_d in the primary clinical target volume (CTV) by an average of 23%. This was achieved through a redistribution of the proton track-ends, concurrently keeping the physical dose distribution virtually unchanged compared to the DOSEopt plans. This resulted in substantial RBE-weighted dose (D_RBE) reductions, allowing the TEopt plans to meet all clinical goals for both RBE models and reduce the NTCPs by 0 to 19 percentage points compared to the DOSEopt plans, assuming the Wedenberg RBE model. The DOSEopt plans met all clinical goals assuming a RBE of 1.1 but failed 10 of 19 normal tissue goals assuming the Wedenberg RBE model.

CONCLUSIONS

Proton track-end objectives allow for LET_d reductions in critical structures without compromising the physical target dose. This approach permits the lowering of D_RBE and NTCP in critical structures, independent of the variable RBE model used, and it could be introduced in clinical practice without changing current protocols based on the constant RBE of 1.1.