Variable RBE in proton therapy: comparison of different model predictions and their influence on clinical-like scenarios

Study of linear energy transfer effect on rib fracture in breast cancer patients receiving pencil-beam-scanning proton therapy

BACKGROUND

In breast cancer patients treated with pencil-beam scanning proton therapy (PBS), the increased linear energy transfer (LET) near the end of the proton range can affect nearby ribs. This may associate with a higher risk of rib fractures.

PURPOSE

To study the effect of LET on rib fracture in breast cancer patients treated with PBS using a novel tool of dose-LET volume histogram (DLVH).

METHODS

From a prospective registry of patients treated with post-mastectomy proton therapy to the chest wall and regional lymph nodes for breast cancer between 2015 and 2020, we retrospectively identified rib fracture cases detected after completing treatment. Contemporaneously treated control patients who did not develop rib fracture were matched to patients 2:1 considering prescription dose, boost location, reconstruction status, laterality, chest wall thickness, and treatment year. The DLVH index, V(d, l), defined as volume(V) of the structure with at least dose(d) and dose-averaged LET (l) (LETd), was calculated. DLVH plots between the fracture and control group were compared. Conditional logistic regression (CLR) model was used to establish the relation of V(d, l) and the observed fracture at each combination of d and l. The p-value derived from CLR model shows the statistical difference between fracture patients and the matched control group. Using the 2D p-value map derived from CLR model, the DLVH features associated with the patient outcomes were extracted.

RESULTS

Seven rib fracture patients were identified, and fourteen matched patients were selected for the control group. The median time from the completion of proton therapy to rib fracture diagnosis was 12 months (range 5-14 months). Two patients had grade 2 symptomatic rib fracture while the remaining 5 were grade 1 incidentally detected on imaging. The derived p-value map demonstrated larger V(0-36 Gy[RBE], 4.0-5.0 keV/µm) in patients experiencing fracture (p < 0.1). For example, the p-value for V(30 Gy[RBE], 4.0 keV/um) was 0.069.

CONCLUSION

In breast cancer patients receiving PBS, a larger volume of chest wall receiving moderate dose and high LETd may result in an increased risk of rib fracture.

Efficient proton transport modelling for proton beam therapy and biological quantification

In this work, we present a fundamental mathematical model for proton transport, tailored to capture the key physical processes underpinning Proton Beam Therapy (PBT). The model provides a robust and computationally efficient framework for exploring various aspects of PBT, including dose delivery, linear energy transfer, treatment planning and the evaluation of relative biological effectiveness. Our findings highlight the potential of this model as a complementary tool to more complex and computationally intensive simulation techniques currently used in clinical practice.

Enabling clinical use of linear energy transfer in proton therapy for head and neck cancer - A review of implications for treatment planning and adverse events study

Proton therapy offers significant advantages due to its unique physical and biological properties, particularly the Bragg peak, enabling precise dose delivery to tumors while sparing healthy tissues. However, the clinical implementation is challenged by the oversimplification of the relative biological effectiveness (RBE) as a fixed value of 1.1, which does not account for the complex interplay between dose, linear energy transfer (LET), and biological endpoints. Lack of heterogeneity control or the understanding of the complex interplay may result in unexpected adverse events and suboptimal patient outcomes. On the other hand, expanding our knowledge of variable tumor RBE and LET optimization may provide a better management strategy for radioresistant tumors. This review examines recent advancements in LET calculation methods, including analytical models and Monte Carlo simulations. The integration of LET into plan evaluation is assessed to enhance plan quality control. LET-guided robust optimization demonstrates promise in minimizing high-LET exposure to organs at risk, thereby reducing the risk of adverse events. Dosimetric seed spot analysis is discussed to show its importance in revealing the true LET-related effect upon the adverse event initialization by finding the lesion origins and eliminating the confounding factors from the biological processes. Dose-LET volume histograms (DLVH) are discussed as effective tools for correlating physical dose and LET with clinical outcomes, enabling the derivation of clinically relevant dose-LET volume constraints without reliance on uncertain RBE models. Based on DLVH, the dose-LET volume constraints (DLVC)-guided robust optimization is introduced to upgrade conventional dose-volume constraints-based robust optimization, which optimizes the joint distribution of dose and LET simultaneously. In conclusion, translating the advances in LET-related research into clinical practice necessitates a better understanding of the LET-related biological mechanisms and the development of clinically relevant LET-related volume constraints directly derived from the clinical outcomes. Future research is needed to refine these models and conduct prospective trials to assess the clinical benefits of LET-guided optimization on patient outcomes.

Linear approximation of variable relative biological effectiveness models for proton therapy

The McNamara (MCN) and Wedenberg (WED) RBE weighted dose (DRBE), dose and dose-weighted average LET (LETd) were calculated in twenty brain cancer patients. A linear approximation was made for each RBE model to give best agreement to clinically relevant dosimetric parameters. Additional evaluations were done on twenty head and neck and twenty breast cancer patients.The R2 of the fits was ≥0.94 and ≥0.91 for MCN and WED respectively for α/β values ≥1.0 Gy. The graphs derived in this work can be used to convert RBE-LET slopes derived from clinical data to α/β values in the MCN or WED models.

NRG Oncology White Paper on the Relative Biological Effectiveness in Proton Therapy

This position paper, led by the NRG Oncology Particle Therapy Work Group, focuses on the concept of relative biologic effect (RBE) in clinical proton therapy (PT), with the goal of providing recommendations for the next-generation clinical trials with PT on the best practice of investigating and using RBE, which could deviate from the current standard proton RBE value of 1.1 relative to photons. In part 1, current clinical utilization and practice are reviewed, giving the context and history of RBE. Evidence for variation in RBE is presented along with the concept of linear energy transfer (LET). The intertwined nature of tumor radiobiology, normal tissue constraints, and treatment planning with LET and RBE considerations is then reviewed. Part 2 summarizes current and past clinical data and then suggests the next steps to explore and employ tools for improved dynamic models for RBE. In part 3, approaches and methods for the next generation of prospective clinical trials are explored, with the goal of optimizing RBE to be both more reflective of clinical reality and also deployable in trials to allow clinical validation and interpatient comparisons. These concepts provide the foundation for personalized biologic treatments reviewed in part 4. Finally, we conclude with a summary including short- and long-term scientific focus points for clinical PT. The practicalities and capacity to use RBE in treatment planning are reviewed and considered with more biological data in hand. The intermediate step of LET optimization is summarized and proposed as a potential bridge to the ultimate goal of case-specific RBE planning that can be achieved as a hypothesis-generating tool in near-term proton trials.

LET-based approximation of the microdosimetric kinetic model for proton radiotherapy

BACKGROUND

Phenomenological relative biological effectiveness (RBE) models for proton therapy, based on the dose-averaged linear energy transfer (LET), have been developed to address the apparent RBE increase towards the end of the proton range. The results of these phenomenological models substantially differ due to varying empirical assumptions and fitting functions. In contrast, more theory-based approaches are used in carbon ion radiotherapy, such as the microdosimetric kinetic model (MKM). However, implementing microdosimetry-based models in LET-based proton therapy treatment planning systems poses challenges.

PURPOSE

This work presents a LET-based version of the MKM that is practical for clinical use in proton radiotherapy.

METHODS

At first, we derived an approximation of the Mayo Clinic Florida (MCF) MKM for relatively-sparsely ionizing radiation such as protons. The mathematical formalism of the proposed model is equivalent to the original MKM, but it maintains some key features of the MCF MKM, such as the determination of model parameters from measurable cell characteristics. Subsequently, we carried out Monte Carlo calculations with PHITS in different simulated scenarios to establish a heuristic correlation between microdosimetric quantities and the dose averaged LET of protons.

RESULTS

A simple allometric function was found able to describe the relationship between the dose-averaged LET of protons and the dose-mean lineal energy, which includes the contributions of secondary particles. The LET-based MKM was used to model the in vitro clonogenic survival RBE of five human and rodent cell lines (A549, AG01522, CHO, T98G, and U87) exposed to pristine and spread-out Bragg peak (SOBP) proton beams. The results of the LET-based MKM agree well with the biological data in a comparable or better way with respect to the other models included in the study. A sensitivity analysis on the model results was also performed.

CONCLUSIONS

The LET-based MKM integrates the predictive theoretical framework of the MCF MKM with a straightforward mathematical description of the RBE based on the dose-averaged LET, a physical quantity readily available in modern treatment planning systems for proton therapy.

Deep-learning-based segmentation using individual patient data on prostate cancer radiation therapy

PURPOSE

Organ-at-risk segmentation is essential in adaptive radiotherapy (ART). Learning-based automatic segmentation can reduce committed labor and accelerate the ART process. In this study, an auto-segmentation model was developed by employing individual patient datasets and a deep-learning-based augmentation method for tailoring radiation therapy according to the changes in the target and organ of interest in patients with prostate cancer.

METHODS

Two computed tomography (CT) datasets with well-defined labels, including contoured prostate, bladder, and rectum, were obtained from 18 patients. The labels of the CT images captured during radiation therapy (CT2nd) were predicted using CT images scanned before radiation therapy (CT1st). From the deformable vector fields (DVFs) created by using the VoxelMorph method, 10 DVFs were extracted when each of the modified CT and CT2nd images were deformed and registered to the fixed CT1st image. Augmented images were acquired by utilizing 110 extracted DVFs and spatially transforming the CT1st images and labels. An nnU-net autosegmentation network was trained by using the augmented images, and the CT2nd label was predicted. A patient-specific model was created for 18 patients, and the performances of the individual models were evaluated. The results were evaluated by employing the Dice similarity coefficient (DSC), average Hausdorff distance, and mean surface distance. The accuracy of the proposed model was compared with those of models trained with large datasets.

RESULTS

Patient-specific models were developed successfully. For the proposed method, the DSC values of the actual and predicted labels for the bladder, prostate, and rectum were 0.94 ± 0.03, 0.84 ± 0.07, and 0.83 ± 0.04, respectively.

CONCLUSION

We demonstrated the feasibility of automatic segmentation by employing individual patient datasets and image augmentation techniques. The proposed method has potential for clinical application in automatic prostate segmentation for ART.

Quantifying the Dosimetric Impact of Proton Range Uncertainties on RBE-Weighted Dose Distributions in Intensity-Modulated Proton Therapy for Bilateral Head and Neck Cancer

BACKGROUND

In current clinical practice, intensity-modulated proton therapy (IMPT) head and neck cancer (HNC) plans are generated using a constant relative biological effectiveness (cRBE) of 1.1. The primary goal of this study was to explore the dosimetric impact of proton range uncertainties on RBE-weighted dose (RWD) distributions using a variable RBE (vRBE) model in the context of bilateral HNC IMPT plans.

METHODS

The current study included the computed tomography (CT) datasets of ten bilateral HNC patients who had undergone photon therapy. Each patient's plan was generated using three IMPT beams to deliver doses to the CTV_High and CTV_Low for doses of 70 Gy(RBE) and 54 Gy(RBE), respectively, in 35 fractions through a simultaneous integrated boost (SIB) technique. Each nominal plan calculated with a cRBE of 1.1 was subjected to the range uncertainties of ±3%. The McNamara vRBE model was used for RWD calculations. For each patient, the differences in dosimetric metrices between the RWD and nominal dose distributions were compared.

RESULTS

The constrictor muscles, oral cavity, parotids, larynx, thyroid, and esophagus showed average differences in mean dose (Dmean) values up to 6.91 Gy(RBE), indicating the impact of proton range uncertainties on RWD distributions. Similarly, the brachial plexus, brain, brainstem, spinal cord, and mandible showed varying degrees of the average differences in maximum dose (Dmax) values (2.78-10.75 Gy(RBE)). The Dmean and Dmax to the CTV from RWD distributions were within ±2% of the dosimetric results in nominal plans.

CONCLUSION

The consistent trend of higher mean and maximum doses to the OARs with the McNamara vRBE model compared to cRBE model highlighted the need for consideration of proton range uncertainties while evaluating OAR doses in bilateral HNC IMPT plans.

Benchmarking proton RBE models

Objective.To biologically optimise proton therapy, models which can accurately predict variations in proton relative biological effectiveness (RBE) are essential. Current phenomenological models show large disagreements in RBE predictions, due to different model assumptions and differences in the data to which they were fit. In this work, thirteen RBE models were benchmarked against a comprehensive proton RBE dataset to evaluate predictions when all models are fit using the same data and fitting techniques, and to assess the statistical robustness of the models.Approach.Model performance was initially evaluated by fitting to the full dataset, and then a cross-validation approach was applied to assess model generalisability and robustness. The impact of weighting the fit and the choice of biological endpoint (either single or multiple survival levels) was also evaluated.Main results.Fitting the models to a common dataset reduced differences between their predictions, however significant disagreements remained due to different underlying assumptions. All models performed poorly under cross-validation in the weighted fits, suggesting that some uncertainties on the experimental data were significantly underestimated, resulting in over-fitting and poor performance on unseen data. The simplest model, which depends linearly on the LET but has no tissue or dose dependence, performed best for a single survival level. However, when fitting to multiple survival levels simultaneously, more complex models with tissue dependence performed better. All models had significant residual uncertainty in their predictions compared to experimental data.Significance.This analysis highlights that poor quality of error estimation on the dose response parameters introduces substantial uncertainty in model fitting. The significant residual error present in all approaches illustrates the challenges inherent in fitting to large, heterogeneous datasets and the importance of robust statistical validation of RBE models.

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.

Calculation of biological effectiveness of SOBP proton beams: a TOPAS Monte Carlo study

Objective.This study aims to investigate the biological effectiveness of Spread-Out Bragg-Peak (SOBP) proton beams with initial kinetic energies 50-250 MeV at different depths in water using TOPAS Monte Carlo code.Approach.The study modelled SOBP proton beams using TOPAS time feature. Various LET-based models and Repair-Misrepair-Fixation model were employed to calculate Relative Biological Effectiveness (RBE) for V79 cell lines at different on-axis depths based on TOPAS. Microdosimetric Kinetic Model and biological weighting function-based models, which utilize microdosimetric distributions, were also used to estimate the RBE. A phase-space-based method was adopted for calculating microdosimetric distributions.Main results.The trend of variation of RBE with depth is similar in all the RBE models, but the absolute RBE values vary based on the calculation models. RBE sharply increases at the distal edge of SOBP proton beams. In the entrance region of all the proton beams, RBE values at 4 Gy i.e. RBE(4 Gy) resulting from different models are in the range of 1.04-1.07, comparable to clinically used generic RBE of 1.1. Moving from the proximal to distal end of the SOBP, RBE(4 Gy) is in the range of 1.15-1.33, 1.13-1.21, 1.11-1.17, 1.13-1.18 and 1.17-1.21, respectively for 50, 100, 150, 200 and 250 MeV SOBP beams, whereas at the distal dose fall-off region, these values are 1.68, 1.53, 1.44, 1.42 and 1.40, respectively.Significance.The study emphasises application of depth-, dose- and energy- dependent RBE values in clinical application of proton beams.

Exploring Helium Ions' Potential for Post-Mastectomy Left-Sided Breast Cancer Radiotherapy

Proton therapy presents a promising modality for treating left-sided breast cancer due to its unique dose distribution. Helium ions provide increased conformality thanks to a reduced lateral scattering. Consequently, the potential clinical benefit of both techniques was explored. An explorative treatment planning study involving ten patients, previously treated with VMAT (Volumetric Modulated Arc Therapy) for 50 Gy in 25 fractions for locally advanced, node-positive breast cancer, was carried out using proton pencil beam therapy with a fixed relative biological effectiveness (RBE) of 1.1 and helium therapy with a variable RBE described by the mMKM (modified microdosimetric kinetic model). Results indicated that target coverage was improved with particle therapy for both the clinical target volume and especially the internal mammary lymph nodes compared to VMAT. Median dose value analysis revealed that proton and helium plans provided lower dose on the left anterior descending artery (LAD), heart, lungs and right breast than VMAT. Notably, helium therapy exhibited improved ipsilateral lung sparing over protons. Employing NTCP models as available in the literature, helium therapy showed a lower probability of grade ≤ 2 radiation pneumonitis (22% for photons, 5% for protons and 2% for helium ions), while both proton and helium ions reduce the probability of major coronary events with respect to VMAT.

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.

The risk of radiation-induced neurocognitive impairment and the impact of sparing the hippocampus during pediatric proton cranial irradiation

BACKGROUND AND PURPOSE

Hippocampus is a central component for neurocognitive function and memory. We investigated the predicted risk of neurocognitive impairment of craniospinal irradiation (CSI) and the deliverability and effects of hippocampal sparing. The risk estimates were derived from published NTCP models. Specifically, we leveraged the estimated benefit of reduced neurocognitive impairment with the risk of reduced tumor control.

MATERIAL AND METHODS

For this dose planning study, a total of 504 hippocampal sparing intensity modulated proton therapy (HS-IMPT) plans were generated for 24 pediatric patients whom had previously received CSI. Plans were evaluated with respect to target coverage and homogeneity index to target volumes, maximum and mean dose to OARs. Paired t-tests were used to compare hippocampal mean doses and normal tissue complication probability estimates.

RESULTS

The median mean dose to the hippocampus could be reduced from 31.3 Gy_RBE to 7.3 Gy_RBE (p < .001), though 20% of these plans were not considered clinically acceptable as they failed one or more acceptance criterion. Reducing the median mean hippocampus dose to 10.6 Gy_RBE was possible with all plans considered as clinically acceptable treatment plans. By sparing the hippocampus to the lowest dose level, the risk estimation of neurocognitive impairment could be reduced from 89.6%, 62.1% and 51.1% to 41.0% (p < .001), 20.1% (p < .001) and 29.9% (p < .001) for task efficiency, organization and memory, respectively. Estimated tumor control probability was not adversely affected by HS-IMPT, ranging from 78.5 to 80.5% for all plans.

CONCLUSIONS

We present estimates of potential clinical benefit in terms of neurocognitive impairment and demonstrate the possibility of considerably reducing neurocognitive adverse effects, minimally compromising target coverage locally using HS-IMPT.

Clinical applicability of deep learning-based respiratory signal prediction models for four-dimensional radiation therapy

For accurate respiration gated radiation therapy, compensation for the beam latency of the beam control system is necessary. Therefore, we evaluate deep learning models for predicting patient respiration signals and investigate their clinical feasibility. Herein, long short-term memory (LSTM), bidirectional LSTM (Bi-LSTM), and the Transformer are evaluated. Among the 540 respiration signals, 60 signals are used as test data. Each of the remaining 480 signals was spilt into training and validation data in a 7:3 ratio. A total of 1000 ms of the signal sequence (T_s) is entered to the models, and the signal at 500 ms afterward (P_t) is predicted (standard training condition). The accuracy measures are: (1) root mean square error (RMSE) and Pearson correlation coefficient (CC), (2) accuracy dependency on T_s and P_t, (3) respiratory pattern dependency, and (4) error for 30% and 70% of the respiration gating for a 5 mm tumor motion for latencies of 300, 500, and 700 ms. Under standard conditions, the Transformer model exhibits the highest accuracy with an RMSE and CC of 0.1554 and 0.9768, respectively. An increase in T_s improves accuracy, whereas an increase in P_t decreases accuracy. An evaluation of the regularity of the respiratory signals reveals that the lowest predictive accuracy is achieved with irregular amplitude patterns. For 30% and 70% of the phases, the average error of the three models is <1.4 mm for a latency of 500 ms and >2.0 mm for a latency of 700 ms. The prediction accuracy of the Transformer is superior to LSTM and Bi-LSTM. Thus, the three models have clinically applicable accuracies for a latency <500 ms for 10 mm of regular tumor motion. The clinical acceptability of the deep learning models depends on the inherent latency and the strategy for reducing the irregularity of respiration.

CONVERSION OF DOSE DISTRIBUTION TO CELL SURVIVAL FRACTION THROUGH DNA DAMAGE: A MONTE CARLO STUDY

Ionizing radiation plays an important role in cancer treatment. Radiation is able to damage the genetic material of cells, blocking their ability to divide and proliferate further. Since radiation affects both healthy and malignant tissues, for all radiation treatments, the design of an accurate treatment plan is fundamental. Usually, weight factors, such as the relative biological effectiveness, are applied to estimate the impact of the kind of radiation and the irradiated medium on the dose deposition. However, these factors can only provide a partial estimation of the real effect on tissues. In this work, a flexible system that is able to predict cell survival fractions according to the planned dose distribution is presented. Dose deposition and subsequent DNA damage were simulated with a multi-scale modeling approach by first applying the FLUKA Monte Carlo (MC) code to estimate the absorbed doses and fluence energy spectra and then using the MC Damage Simulation code to compute the DNA damage yields. Lastly, the results are converted into cell survival fraction using a theoretical model. The comparisons between the simulated survival fractions with experimental data are reported for a proton spread out Bragg peak at several doses. The presented approach helps to elucidate radiobiological responses along the Bragg curve and has the flexibility to be extended to a wide range of situations of clinical interest.

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.

Roadmap: helium ion therapy

Helium ion beam therapy for the treatment of cancer was one of several developed and studied particle treatments in the 1950s, leading to clinical trials beginning in 1975 at the Lawrence Berkeley National Laboratory. The trial shutdown was followed by decades of research and clinical silence on the topic while proton and carbon ion therapy made debuts at research facilities and academic hospitals worldwide. The lack of progression in understanding the principle facets of helium ion beam therapy in terms of physics, biological and clinical findings persists today, mainly attributable to its highly limited availability. Despite this major setback, there is an increasing focus on evaluating and establishing clinical and research programs using helium ion beams, with both therapy and imaging initiatives to supplement the clinical palette of radiotherapy in the treatment of aggressive disease and sensitive clinical cases. Moreover, due its intermediate physical and radio-biological properties between proton and carbon ion beams, helium ions may provide a streamlined economic steppingstone towards an era of widespread use of different particle species in light and heavy ion therapy. With respect to the clinical proton beams, helium ions exhibit superior physical properties such as reduced lateral scattering and range straggling with higher relative biological effectiveness (RBE) and dose-weighted linear energy transfer (LET_d) ranging from ∼4 keVμm^-1to ∼40 keVμm^-1. In the frame of heavy ion therapy using carbon, oxygen or neon ions, where LET_dincreases beyond 100 keVμm^-1, helium ions exhibit similar physical attributes such as a sharp lateral penumbra, however, with reduced radio-biological uncertainties and without potentially spoiling dose distributions due to excess fragmentation of heavier ion beams, particularly for higher penetration depths. This roadmap presents an overview of the current state-of-the-art and future directions of helium ion therapy: understanding physics and improving modeling, understanding biology and improving modeling, imaging techniques using helium ions and refining and establishing clinical approaches and aims from learned experience with protons. These topics are organized and presented into three main sections, outlining current and future tasks in establishing clinical and research programs using helium ion beams-A. Physics B. Biological and C. Clinical Perspectives.

Range uncertainty reductions in proton therapy may lead to the feasibility of novel beam arrangements which improve organ-at-risk sparing

PURPOSE

In proton therapy, dose distributions are currently often conformed to organs at risk (OARs) using the less sharp dose fall-off at the lateral beam edge to reduce the effects of uncertainties in the in vivo proton range. However, range uncertainty reductions may make greater use of the sharper dose fall-off at the distal beam edge feasible, potentially improving OAR sparing. We quantified the benefits of such novel beam arrangements.

METHODS

For each of 10 brain or skull base cases, five treatment plans robust to 2 mm setup and 0%-4% range uncertainty were created for the traditional clinical beam arrangement and a novel beam arrangement making greater use of the distal beam edge to conform the dose distribution to the brainstem. Metrics including the brainstem normal tissue complication probability (NTCP) with the endpoint of necrosis were determined for all plans and all setup and range uncertainty scenarios.

RESULTS

For the traditional beam arrangement, reducing the range uncertainty from the current level of approximately 4% to a potentially achievable level of 1% reduced the brainstem NTCP by up to 0.9 percentage points in the nominal and up to 1.5 percentage points in the worst-case scenario. Switching to the novel beam arrangement at 1% range uncertainty improved these values by a factor of 2, that is, to 1.8 percentage points and 3.2 percentage points, respectively. The novel beam arrangement achieved a lower brainstem NTCP in all cases starting at a range uncertainty of 2%.

CONCLUSION

The benefits of novel beam arrangements may be of the same magnitude or even exceed the direct benefits of range uncertainty reductions. Indirect effects may therefore contribute markedly to the benefits of reducing proton range uncertainties.

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.

Ganetespib selectively sensitizes cancer cells for proximal and distal spread-out Bragg peak proton irradiation

Objective

Hypersensitivity towards proton versus photon irradiation was demonstrated in homologous recombination repair (HRR)-deficient cell lines. Hence, combined treatment concepts targeting HRR provide a rational for potential pharmaceutical exploitation. The HSP90 inhibitor ganetespib (STA-9090) downregulates a multitude of HRR-associated proteins and sensitizes for certain chemotherapeutics. Thus, the radiosensitizing effect of HSP90-inhibiting ganetespib was investigated for reference photon irradiation and proton irradiation at a proximal and distal position in a spread-out Bragg peak (SOBP).

Methods

A549 and FaDu cells were treated with low-dose (2 nM resp. 1 nM) ganetespib and irradiated with 200 kV photons. Proton irradiation was performed at a proximal and a distal position within a SOBP, with corresponding dose-averaged linear-energy transfer (LET_D) values of 2.1 and 4.5 keV/µm, respectively. Cellular survival data was fitted to the linear-quadratic model to calculate relative biological effectiveness (RBE) and the dose-modifying factor (DMF). Additionally, A549 cells were treated with increasing doses of ganetespib and investigated by flow cytometry, immunoblotting, and immunofluorescence microscopy to investigate cell cycle distribution, Rad51 protein levels, and γH2AX foci, respectively.

Results

Low-dosed ganetespib significantly sensitized both cancer cell lines exclusively for proton irradiation at both investigated LET_D, resulting in increased RBE values of 10–40%. In comparison to photon irradiation, the fraction of cells in S/G2/M phase was elevated in response to proton irradiation with 10 nM ganetespib consistently reducing this population. No changes in cell cycle distribution were detected in unirradiated cells by ganetespib alone. Protein levels of Rad51 are downregulated in irradiated A549 cells by 10 nM and also 2 nM ganetespib within 24 h. Immunofluorescence staining demonstrated similar induction and removal of γH2AX foci, irrespective of irradiation type or ganetespib administration.

Conclusion

Our findings illustrate a proton-specific sensitizing effect of low-dosed ganetespib in both employed cell lines and at both investigated SOBP positions. We provide additional experimental data on cellular response and a rational for future combinatorial approaches with proton radiotherapy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-022-02036-z.

Quantification of biological range uncertainties in patients treated at the Krakow proton therapy centre

BACKGROUND

Variable relative biological effectiveness (vRBE) in proton therapy might significantly modify the prediction of RBE-weighted dose delivered to a patient during proton therapy. In this study we will present a method to quantify the biological range extension of the proton beam, which results from the application of vRBE approach in RBE-weighted dose calculation.

METHODS AND MATERIALS

The treatment plans of 95 patients (brain and skull base patients) were used for RBE-weighted dose calculation with constant and the McNamara RBE model. For this purpose the Monte Carlo tool FRED was used. The RBE-weighted dose distributions were analysed using indices from dose-volume histograms. We used the volumes receiving at least 95% of the prescribed dose (V95) to estimate the biological range extension resulting from vRBE approach.

RESULTS

The vRBE model shows higher median value of relative deposited dose and D95 in the planning target volume by around 1% for brain patients and 4% for skull base patients. The maximum doses in organs at risk calculated with vRBE was up to 14 Gy above dose limit. The mean biological range extension was greater than 0.4 cm.

DISCUSSION

Our method of estimation of biological range extension is insensitive for dose inhomogeneities and can be easily used for different proton plans with intensity-modulated proton therapy (IMPT) optimization. Using volumes instead of dose profiles, which is the common method, is more universal. However it was tested only for IMPT plans on fields arranged around the tumor area.

CONCLUSIONS

Adopting a vRBE model results in an increase in dose and an extension of the beam range, which is especially disadvantageous in cancers close to organs at risk. Our results support the need to re-optimization of proton treatment plans when considering vRBE.

Quantifying Systematic RBE-Weighted Dose Uncertainty Arising from Multiple Variable RBE Models in Organ at Risk

PURPOSE

Relative biological effectiveness (RBE) uncertainties have been a concern for treatment planning in proton therapy, particularly for treatment sites that are near organs at risk (OARs). In such a clinical situation, the utilization of variable RBE models is preferred over constant RBE model of 1.1. The problem, however, lies in the exact choice of RBE model, especially when current RBE models are plagued with a host of uncertainties. This paper aims to determine the influence of RBE models on treatment planning, specifically to improve the understanding of the influence of the RBE models with regard to the passing and failing of treatment plans. This can be achieved by studying the RBE-weighted dose uncertainties across RBE models for OARs in cases where the target volume overlaps the OARs. Multi-field optimization (MFO) and single-field optimization (SFO) plans were compared in order to recommend which technique was more effective in eliminating the variations between RBE models.

METHODS

Fifteen brain tumor patients were selected based on their profile where their target volume overlaps with both the brain stem and the optic chiasm. In this study, 6 RBE models were analyzed to determine the RBE-weighted dose uncertainties. Both MFO and SFO planning techniques were adopted for the treatment planning of each patient. RBE-weighted dose uncertainties in the OARs are calculated assuming( α β ) x of 3 Gy and 8 Gy. Statistical analysis was used to ascertain the differences in RBE-weighted dose uncertainties between MFO and SFO planning. Additionally, further investigation of the linear energy transfer (LET) distribution was conducted to determine the relationship between LET distribution and RBE-weighted dose uncertainties.

RESULTS

The results showed no strong indication on which planning technique would be the best for achieving treatment planning constraints. MFO and SFO showed significant differences (P <.05) in the RBE-weighted dose uncertainties in the OAR. In both clinical target volume (CTV)-brain stem and CTV-chiasm overlap region, 10 of 15 patients showed a lower median RBE-weighted dose uncertainty in MFO planning compared with SFO planning. In the LET analysis, 8 patients (optic chiasm) and 13 patients (brain stem) showed a lower mean LET in MFO planning compared with SFO planning. It was also observed that lesser RBE-weighted dose uncertainties were present with MFO planning compared with SFO planning technique.

CONCLUSIONS

Calculations of the RBE-weighted dose uncertainties based on 6 RBE models and 2 different( α β ) x revealed that MFO planning is a better option as opposed to SFO planning for cases of overlapping brain tumor with OARs in eliminating RBE-weighted dose uncertainties. Incorporation of RBE models failed to dictate the passing or failing of a treatment plan. To eliminate RBE-weighted dose uncertainties in OARs, the MFO planning technique is recommended for brain tumor when CTV and OARs overlap.

Hadrontherapy for Thymic Epithelial Tumors: Implementation in Clinical Practice

Radiation therapy is part of recommendations in the adjuvant settings for advanced stage or as exclusive treatment in unresectable thymic epithelial tumors (TETs). However, first-generation techniques delivered substantial radiation doses to critical organs at risk (OARs), such as the heart or the lungs, resulting in noticeable radiation-induced toxicity. Treatment techniques have significantly evolved for TET irradiation, and modern techniques efficiently spare normal surrounding tissues without negative impact on tumor coverage and consequently local control or patient survival. Considering its dosimetric advantages, hadrontherapy (which includes proton therapy and carbon ion therapy) has proved to be worthwhile for TET irradiation in particular for challenging clinical situations such as cardiac tumoral involvement. However, clinical experience for hadrontherapy is still limited and mainly relies on small-size proton therapy studies. This critical review aims to analyze the current status of hadrontherapy for TET irradiation to implement it at a larger scale.

Proton Therapy for Breast Cancer: A Consensus Statement From the Particle Therapy Cooperative Group Breast Cancer Subcommittee

Radiation therapy plays an important role in the multidisciplinary management of breast cancer. Recent years have seen improvements in breast cancer survival and a greater appreciation of potential long-term morbidity associated with the dose and volume of irradiated organs. Proton therapy reduces the dose to nontarget structures while optimizing target coverage. However, there remain additional financial costs associated with proton therapy, despite reductions over time, and studies have yet to demonstrate that protons improve upon the treatment outcomes achieved with photon radiation therapy. There remains considerable heterogeneity in proton patient selection and techniques, and the rapid technological advances in the field have the potential to affect evidence evaluation, given the long latency period for breast cancer radiation therapy recurrence and late effects. In this consensus statement, we assess the data available to the radiation oncology community of proton therapy for breast cancer, provide expert consensus recommendations on indications and technique, and highlight ongoing trials' cost-effectiveness analyses and key areas for future research.

Proton range uncertainty reduction benefits for skull base tumors in terms of normal tissue complication probability (NTCP) and healthy tissue doses

PURPOSE

Proton therapy allows for more conformal dose distributions and lower organ at risk and healthy tissue doses than conventional photon-based radiotherapy, but uncertainties in the proton range currently prevent proton therapy from making full use of these advantages. Numerous developments therefore aim to reduce such range uncertainties. In this work, we quantify the benefits of reductions in range uncertainty for treatments of skull base tumors.

METHODS

The study encompassed 10 skull base patients with clival tumors. For every patient, six treatment plans robust to setup errors of 2 mm and range errors from 0% to 5% were created. The determined metrics included the brainstem and optic chiasm normal tissue complication probability (NTCP) with the endpoints of necrosis and blindness, respectively, as well as the healthy tissue volume receiving at least 70% of the prescription dose.

RESULTS

A range uncertainty reduction from the current level of 4% to a potentially achievable level of 1% reduced the probability of brainstem necrosis by up to 1.3 percentage points in the nominal scenario in which neither setup nor range errors occur and by up to 2.9 percentage points in the worst-case scenario. Such a range uncertainty reduction also reduced the optic chiasm NTCP with the endpoint of blindness by up to 0.9 percentage points in the nominal scenario and by up to 2.2 percentage points in the worst-case scenario. The decrease in the healthy tissue volume receiving at least 70% of the prescription dose ranged from -7.8 to 24.1 cc in the nominal scenario and from -3.4 to 38.4 cc in the worst-case scenario.

CONCLUSION

The benefits quantified as part of this study serve as a guideline of the OAR and healthy tissue dose benefits that range monitoring techniques may be able to achieve. Benefits were observed between all levels of range uncertainty. Even smaller range uncertainty reductions may therefore be beneficial.

A Critical Review of LET-Based Intensity-Modulated Proton Therapy Plan Evaluation and Optimization for Head and Neck Cancer Management

In this review article, we review the 3 important aspects of linear-energy-transfer (LET) in intensity-modulated proton therapy (IMPT) for head and neck (H&N) cancer management. Accurate LET calculation methods are essential for LET-guided plan evaluation and optimization, which can be calculated either by analytical methods or by Monte Carlo (MC) simulations. Recently, some new 3D analytical approaches to calculate LET accurately and efficiently have been proposed. On the other hand, several fast MC codes have also been developed to speed up the MC simulation by simplifying nonessential physics models and/or using the graphics processor unit (GPU)–acceleration approach. Some concepts related to LET are also briefly summarized including (1) dose-weighted versus fluence-weighted LET; (2) restricted versus unrestricted LET; and (3) microdosimetry versus macrodosimetry. LET-guided plan evaluation has been clinically done in some proton centers. Recently, more and more studies using patient outcomes as the biological endpoint have shown a positive correlation between high LET and adverse events sites, indicating the importance of LET-guided plan evaluation in proton clinics. Various LET-guided plan optimization methods have been proposed to generate proton plans to achieve biologically optimized IMPT plans. Different optimization frameworks were used, including 2-step optimization, 1-step optimization, and worst-case robust optimization. They either indirectly or directly optimize the LET distribution in patients while trying to maintain the same dose distribution and plan robustness. It is important to consider the impact of uncertainties in LET-guided optimization (ie, LET-guided robust optimization) in IMPT, since IMPT is sensitive to uncertainties including both the dose and LET distributions. We believe that the advancement of the LET-guided plan evaluation and optimization will help us exploit the unique biological characteristics of proton beams to improve the therapeutic ratio of IMPT to treat H&N and other cancers.

Impact of Multiple Beams on Plan Quality, Linear Energy Transfer Distribution, and Plan Robustness of Intensity Modulated Proton Therapy for Lung Cancer

The increase of proton beam number might provide higher degrees of freedom in the optimization of intensity-modulated proton therapy planning. In this study, we aimed to quantitatively explore the potential benefits of the increased beam number, including dose volume histogram (DVH), linear energy transfer volume histogram, and DVH bandwidth metrics. Twelve patients with lung cancer are retrospectively selected. Four plans were created based on internal target volume (ITV) robust optimization for each patient using the RayStation treatment planning system. Four plans were generated using different numbers (three, five, seven, and nine) of evenly separated coplanar beams. The three-beam plan was considered as the reference plan. Biologically equivalent doses were calculated using both constant relative biological effectiveness (RBE) and variable RBE models, respectively. To evaluate plan quality, DVH metrics in the target [ITV: D_2%, CI, HI] and organs-at-risk [Lung: V_5Gy[RBE], V_20Gy[RBE], V_30Gy[RBE]; Heart D_2%; Spinal cord D_2%] were calculated using both RBE models. To evaluate LET distributions, LET volume histogram metrics [ITV LET_mean and LET_2%; Lung LET_mean and LET_2%; Heart LET_2%; Spinal cord LET_2%] were quantified. To evaluate plan robustness, the metrics using DVH bandwidth [ITV: D_2%, D_99%; Lung: V_5Gy[RBE], V_20Gy[RBE], V_30Gy[RBE]; Heart D_2%; Spinal cord D_2%] were also reported. For plan quality, the increase of proton beam number resulted in fewer target hot spots, improved target dose conformity, improved target dose homogeneity, lower median-dose lung volume, and fewer hot spots in spinal cord. As to LET distributions, target mean LET increased significantly as the beam number increased to seven or more. Lung LET hot spots were significantly reduced with the increase of proton beams. With respect to plan robustness, the robustness of target dose coverage, target hot spots, and low-dose lung volume were improved, while the robustness of heart hot spots became worse as the beam number increased to nine. The robustness of cord hot spots became worse using five and seven beams compared to that using three beams. As the proton beam number increased, plan quality and LET distributions were comparable or significantly improved. The robustness of target dose coverage, target dose hot spots, and low-dose lung volume were significantly improved.

Double-strand breaks measured along a 160 MeV proton Bragg curve using a novel FIESTA-DNA probe in a cell-free environment

Proton radiotherapy treatment planning systems use a constant relative biological effectiveness (RBE) = 1.1 to convert proton absorbed dose into biologically equivalent high-energy photon dose. This method ignores linear energy transfer (LET) distributions, and RBE is known to change as a function of LET. Variable RBE approaches have been proposed for proton planning optimization. Experimental validation of models underlying these approaches is a pre-requisite for their clinical implementation. This validation has to probe every level in the evolution of radiation-induced biological damage leading to cell death, starting from DNA double-strand breaks (DSB). Using a novel FIESTA-DNA probe, we measured the probability of double-strand break (P _DSB) along a 160 MeV proton Bragg curve at two dose levels (30 and 60 Gy (RBE)) and compared it to measurements in a 6 MV photon beam. A machined setup that held an Advanced Markus parallel plate chamber for proton dose verification alongside the probes was fabricated. Each sample set consisted of five 10 μl probes suspended inside plastic microcapillary tubes. These were irradiated with protons to 30 Gy (RBE) at depths of 5-17.5 cm and 60 Gy (RBE) at depths of 10-17.2 cm with 1 mm resolution around Bragg peak. Sample sets were also irradiated using 6MV photons to 20, 40, 60, and 80 Gy. For the 30 Gy (RBE) measurements, increases in P _DSB/Gy were observed at 17.0 cm followed by decreases at larger depth. For the 60 Gy (RBE) measurements, no increase in P _DSB/Gy was observed, but there was a decrease after 17.0 cm. Dose-response for P _DSB between 30 and 60 Gy (RBE) showed less than doubling of P _DSB when dose was doubled. Proton RBE effect from DSB, RBE_P,DSB, was <1 except at the Bragg peak. The experiment showed that the novel probe can be used to perform DNA DSB measurements in a proton beam. To establish relevance to clinical environment, further investigation of the probe's chemical scavenging needs to be performed.

Incorporation of the LETd-weighted biological dose in the evaluation of breast intensity-modulated proton therapy plans

PURPOSE

To evaluate the LETd-weighted biological dose to OARs in proton therapy for breast cancer and to study the relationship of the LETd-weighted biological dose relative to the standard dose (RBE = 1.1) and thereby to provide estimations of the biological dose uncertainties with the standard dose calculations (RBE = 1.1) commonly used in clinical practice.

METHOD

This study included 20 patients who received IMPT treatment to the whole breast/chest wall and regional lymph nodes. The LETd distributions were calculated along with the physical dose using an open-source Monte Carlo simulation package, MCsquare. Using the McMahon linear model, the LETd-weighted biological dose was computed from the physical dose and LETd. OAR doses were compared between the Dose (RBE = 1.1) and the LETd-weighted biological dose, on brachial plexus, rib, heart, esophagus, and Ipsilateral lung.

RESULTS

On average, the LETd-weighted biological dose compared to the Dose (RBE = 1.1) was higher by 8% for the brachial plexus D0.1 cc, 13% for the ribs D0.5 cc, 24% for mean heart dose, and 10% for the esophagus D0.1 cc, respectively. The LETd-weighted doses to the Ipsilateral lung V5, V10, and V20 were comparable to the Dose (RBE = 1.1). No statistically significant difference in biological dose enhancement to OARs was observed between the intact breast group and the CW group, with the exception of the ribs: the CW group experienced slightly greater biological dose enhancement (13% vs. 12%, p = 0.04) to the ribs than the intact breast group.

CONCLUSION

Enhanced biological dose was observed compared to standard dose with assumed RBE of 1.1 for the heart, ribs, esophagus, and brachial plexus in breast/CW and regional nodal IMPT plans. Variable RBE models should be considered in the evaluation of the IMPT breast plans, especially for OARs located near the end of range of a proton beam. Clinical outcome studies are needed to validate model predictions for clinical toxicities.

Spatial Agreement of Brainstem Dose Distributions Depending on Biological Model in Proton Therapy for Pediatric Brain Tumors

Purpose

During radiation therapy for pediatric brain tumors, the brainstem is a critical organ at risk, possibly with different radio-sensitivity across its substructures. In proton therapy, treatment planning is currently performed using a constant relative biological effectiveness (RBE) of 1.1 (RBE_1.1), whereas preclinical studies point toward spatial variability of this factor. To shed light on this biological uncertainty, we investigated the spatial agreement between isodose maps produced by different RBE models, with emphasis on (smaller) substructures of the brainstem.

Methods and Materials

Proton plans were recalculated using Monte Carlo simulations in 3 anonymized pediatric patients with brain tumors (a craniopharyngioma, a low-grade glioma, and a posterior fossa ependymoma) to obtain dose and linear energy transfer distributions. Doses and volume metrics for the brainstem and its substructures were calculated using a constant RBE_1.1, 4 phenomenological RBE models with varying (α/β)_x parameters, and with a simpler linear energy transfer-dependent model. The spatial agreement between the dose distributions of constant RBE_1.1 versus the variable RBE models was compared using the Dice similarity coefficient.

Results

The spatial agreement between the variable RBE dose distributions and RBE_1.1 decreased with increasing isodose levels in all patient cases. The patient with ependymoma showed the greatest variation in dose and dose volumes, where V_50Gy(RBE) in the brainstem increased from 32% (RBE_1.1) to 35% to 49% depending on the applied model, corresponding to a spatial agreement (Dice similarity coefficient) between 0.79 and 0.95. The remaining patients showed similar trends, however, with lower absolute values due to lower brainstem doses.

Conclusions

All phenomenological RBE models fully enclosed the isodose volumes of the constant RBE_1.1, and the volumes based on variable RBE spatially agreed. The spatial agreement was dependent on the isodose level, where higher isodose levels showed larger expansions and less agreement between the variable RBE models and RBE_1.1.

Implementation of a double scattering nozzle for Monte Carlo recalculation of proton plans with variable relative biological effectiveness

A constant relative biological effectiveness (RBE) of 1.1 is currently used in clinical proton therapy. However, theRBEvaries with factors such as dose level, linear energy transfer (LET) and tissue type. MultipleRBEmodels have been developed to account for this biological variation. To enable recalculation of patients treated with double scattering (DS) proton therapy, includingLETand variableRBE, we implemented and commissioned a Monte Carlo (MC) model of a DS treatment nozzle. The main components from the IBA nozzle were implemented in the FLUKA MC code. We calibrated and verified the following entities to experimental measurements: range of pristine Bragg peaks (PBPs) and spread-out Bragg peaks (SOBPs), energy spread, lateral profiles, compensator range degradation, and absolute dose. We recalculated two patients with different field setups, comparing FLUKA vs. treatment planning system (TPS) dose, also obtainingLETand variableRBEdoses. We achieved good agreement between FLUKA and measurements. The range differences between FLUKA and measurements were for the PBPs within ±0.9 mm (83% ⩽ 0.5 mm), and for SOBPs ±1.6 mm (82% ⩽ 0.5 mm). The differences in modulation widths were below 5 mm (79% ⩽ 2 mm). The differences in the distal dose fall off (D80%-D20%) were below 0.5 mm for all PBPs and the lateral penumbras diverged from measurements by less than 1 mm. The mean dose difference (RBE= 1.1) in the target between the TPS and FLUKA were below 0.4% in a three-field plan and below 1.4% in a four-field plan. A dose increase of 9.9% and 7.2% occurred when using variableRBEfor the two patients, respectively. We presented a method to recalculate DS proton plans in the FLUKA MC code. The implementation was used to obtainLETand variableRBEdose and can be used for investigating variableRBEfor previously treated patients.

Impact of range uncertainty on clinical distributions of linear energy transfer and biological effectiveness in proton therapy

PURPOSE

Increased radiation response after proton irradiation, such as late radiation-induced toxicity, is determined by high dose and elevated linear energy transfer (LET). Steep dose-averaged LET (LET_d ) gradients and elevated LET_d occur at the end of proton range and might be particularly sensitive to uncertainties in range prediction. Therefore, this study quantified LET_d distributions and the impact of range uncertainty in robust dose-optimized proton treatment plans and assessed the biological effect in normal tissues and tumors of patients.

METHODS

For each of six cancer patients (two brain, head-and-neck, and prostate), two nominal treatment plans were robustly dose optimized using single- and multi-field optimization, respectively. For each plan, two additional scenarios with ±3.5% range deviation relative to the nominal plan were derived by global rescaling of stopping-power ratios. Dose and LET_d distributions were calculated for each scenario using the beam parameters of the corresponding nominal plan. The variability in relative biological effectiveness (RBE) and probability of late radiation-induced brain toxicity (P_IC ) was assessed.

RESULTS

The optimization technique (single- vs multi-field) had a negligible impact on the LET_d distributions in the clinical target volume (CTV) and in most organs at risk (OARs). LET_d distributions in the CTV were rather homogeneous with arithmetic mean of LET_d below 3.2 keV/µm and robust against range deviations. The RBE variability within the CTV induced by range uncertainty was small (≤0.05, 95% confidence interval). In OARs, LET_d hotspots (>7 keV/µm) occurred and LET_d distributions were inhomogeneous and sensitive to range deviations. LET_d hotspots and the impact of range deviations were most prominent in OARs of brain tumor patients which translated in RBE values exceeding 1.1 in all brain OARs. The near-maximum predicted P_IC in healthy brain tissue of brain tumor patients was smaller than 5% and occurred adjacent to the CTV. Range deviations induced absolute differences in P_IC up to 1.2%.

CONCLUSIONS

Robust dose optimization generates LET_d distributions in the target volume robust against range deviations. The current findings support using a constant RBE within the CTV. The impact of range deviations on the considered probability of late radiation-induced toxicity in brain tissue was limited for robust dose-optimized treatment plans. Incorporation of LET_d in robust optimization frameworks may further reduce uncertainty related to the RBE-weighted dose estimation in normal tissues.

Dosimetric uncertainties impact on cell survival curve with low energy proton

There is an increasing number of radiobiological experiments being conducted with low energy protons (less than 5 MeV) for radiobiological studies due to availability of sub-millimetre focused beam. However, low energy proton has broad microdosimetric spectra which can introduce dosimetric uncertainty. In this work, we quantify the impact of this dosimetric uncertainties on the cell survival curve and how it affects the estimation of the alpha and beta parameters in the LQ formalism. Monte Carlo simulation is used to generate the microdosimetric spectra in a micrometer-sized water sphere under proton irradiation. This is modelled using radiobiological experiment set-up at the Centre of Ion Beam Application (CIBA) in National University of Singapore. Our results show that the microdosimetric spectra can introduce both systematic and random shifts in dose and cell survival; this effect is most pronounced with low energy protons. The alpha and beta uncertainties can be up to 10% and above 30%, respectively for low energy protons passing through thin cell target (about 10 microns). These uncertainties are non-negligible and show that care must be taken in using the cell survival curve and its derived parameters for radiobiological models.

A biological effect-guided optimization approach using beam distal-edge avoidance for intensity-modulated proton therapy

PURPOSE

Linear energy transfer (LET)-guided methods have been applied to intensity-modulated proton therapy (IMPT) to improve its biological effect. However, using LET as a surrogate for biological effect ignores the topological relationship of the scanning spot to different structures of interest. In this study, we developed an optimization method that takes advantage of the continuing increase in LET beyond the physical dose Bragg peak. This method avoids placing high biological effect values in critical structures and increases biological effect in the tumor area without compromising target coverage.

METHODS

We selected the cases of two patients with brain tumors and two patients with head and neck tumors who had been treated with proton therapy at our institution. Three plans were created for each case: a plan based on conventional dose-based optimization (DoseOpt), one based on LET-incorporating optimization (LETOpt), and one based on the proposed distal-edge avoidance-guided optimization method (DEAOpt). In DEAOpt, an L₁ -norm sparsity term, in which the penalty of each scanning spot was set according to the topological relationship between the organ positions and the location of the peak scaled LET-weighted dose (c LETxD) was added to a conventional dose-based optimization objective function. All plans were normalized to give the same target dose coverage. Dose (assuming a constant relative biological effectiveness value of 1.1, as in clinical practice), biological effect (c LETxD), and computing time consumption were evaluated and compared among the three optimization approaches for each patient case.

RESULTS

For all four cases, all three optimization methods generated comparable dose coverage in both target and critical structures. The LETOpt plans and DEAOpt plans reduced biological effect hot spots in critical structures and increased biological effect in the target volumes to a similar extent. For the target, the c LETxD_98% and c LETxD_2% in the DEAOpt plans were on average 7.2% and 11.74% higher than in the DoseOpt plans, respectively. For the brainstem, the c LETxD_mean in the DEAOpt plans was on average 33.38% lower than in the DoseOpt plans. In addition, the DEAOpt method saved 30.37% of the computation cost over the LETOpt method.

CONCLUSIONS

DEAOpt is an alternative IMPT optimization approach that correlates the location of scanning spots with biological effect distribution. IMPT could benefit from the use of DEAOpt because this method not only delivers comparable biological effects to LETOpt plans, but also is faster.

Cancer risk after breast proton therapy considering physiological and radiobiological uncertainties

BACKGROUND

The reduced normal tissue dose burden from protons can reduce the risk of second cancer for breast cancer patients. Breathing motion and the impact of variable relative biological effectiveness (RBE) are however concerns for proton dose distributions. This study aimed to quantify the impact of these factors on risk predictions from proton and photon therapy.

MATERIALS AND METHODS

Twelve patients were planned in free breathing with protons and photons to deliver 50 Gy (RBE) in 25 fractions (assuming RBE = 1.1 for protons) to the left breast. Second cancer risk was evaluated with several models for the lungs, contralateral breast, heart and esophagus as organs at risk (OARs). Plans were recalculated on CT-datasets acquired in extreme phases to account for breathing motion. Proton plans were also recalculated assuming variable RBE for a range of radiobiological parameters.

RESULTS

The OARs received substantially lower doses from protons compared to photons. The highest risks were for the lungs (average second cancer risks of 0.31% and 0.12% from photon and proton plans, respectively). The reduced risk with protons was maintained, even when breathing and/or RBE variation were taken into account. Furthermore, while the total risks from the photon plans were seen to increase with the integral dose, no such correlation was observed for the proton plans.

CONCLUSIONS

Protons have an advantage over the photons with respect to the induction of cancer. Uncertainties in physiological movements and radiobiological parameters affected the absolute risk estimates, but not the general trend of lower risk associated with proton therapy.

Validation of linear energy transfer computed in a Monte Carlo dose engine of a commercial treatment planning system

The relative biological effectiveness (RBE) of protons is highly variable and difficult to quantify. However, RBE is related to the local ionization density, which can be related to the physical measurable dose weighted linear energy transfer (LET_D). The aim of this study was to validate the LET_D calculations for proton therapy beams implemented in a commercially available treatment planning system (TPS) using microdosimetry measurements and independent LET_D calculations (Open-MCsquare (MCS)). The TPS (RayStation v6R) was used to generate treatment plans on the CIRS-731-HN anthropomorphic phantom for three anatomical sites (brain, nasopharynx, neck) for a spherical target (Ø  =  5 cm) with uniform target dose to calculate the LET_D distribution. Measurements were performed at the University Medical Center Groningen proton therapy center (Proteus Plus, IBA) using a µ ⁺-probe utilizing silicon on insulator microdosimeters capable of detecting lineal energies as low as 0.15 keV µm^-1 in tissue. Dose averaged mean lineal energy [Formula: see text] depth-profiles were measured for 70 and 130 MeV spots in water and for the three treatment plans in water and an anthropomorphic phantom. The [Formula: see text] measurements were compared to the LET_D calculated in the TPS and MCS independent dose calculation engine. D · [Formula: see text] was compared to D · LET_D in terms of a gamma-index with a distance-to-agreement criteria of 2 mm and increasing dose difference criteria to determine the criteria for which a 90% pass rate was accomplished. Measurements of D · [Formula: see text] were in good agreement with the D · LET_D calculated in the TPS and MCS. The 90% passing rate threshold was reached at different D · LET_D difference criteria for single spots (TPS: 1% MCS: 1%), treatment plans in water (TPS: 3% MCS: 6%) and treatment plans in an anthropomorphic phantom (TPS: 6% MCS: 1%). We conclude that D · LET_D calculations accuracy in the RayStation TPS and open MCSquare are within 6%, and sufficient for clinical D · LET_D evaluation and optimization. These findings remove an important obstacle in the road towards clinical implementation of D · LET_D evaluation and optimization of proton therapy treatment plans. Novelty and significance The dose weighed linear energy transfer (LET_D) distribution can be calculated for proton therapy treatment plans by Monte Carlo dose engines. The relative biological effectiveness (RBE) of protons is known to vary with the LET_D distribution. Therefore, there exists a need for accurate calculation of clinical LET_D distributions. Previous LET_D validations have focused on general purpose Monte Carlo dose engines which are typically not used clinically. We present the first validation of mean lineal energy [Formula: see text] measurements of the LET_D against calculations by the Monte Carlo dose engines of the Raystation treatment planning system and open MCSquare.

The impact of proton LET/RBE modeling and robustness analysis on base-of-skull and pediatric craniopharyngioma proton plans relative to VMAT

Purpose: Pediatric craniopharyngioma, adult base-of-skull sarcoma and chordoma cases are all regarded as priority candidates for proton therapy. In this study, a dosimetric comparison between volumetric modulated arc therapy (VMAT) and intensity modulated proton therapy (IMPT) was first performed. We then investigated the impact of physical and biological uncertainties. We assessed whether IMPT plans remained dosimetrically superior when such uncertainty estimates were considered, especially with regards to sparing organs at risk (OARs).Methodology: We studied 10 cases: four chondrosarcoma, two chordoma and four pediatric craniopharyngioma. VMAT and IMPT plans were created according to modality-specific protocols. For IMPT, we considered (i) variable RBE modeling using the McNamara model for different values of (α/β)x, and (ii) robustness analysis with ±3 mm set-up and 3.5% range uncertainties.Results: When comparing the VMAT and IMPT plans, the dosimetric advantages of IMPT were clear: IMPT led to reduced integral dose and, typically, improved CTV coverage given our OAR constraints. When physical robustness analysis was performed for IMPT, some uncertainty scenarios worsened the CTV coverage but not usually beyond that achieved by VMAT. Certain scenarios caused OAR constraints to be exceeded, particularly for the brainstem and optical chiasm. However, variable RBE modeling predicted even more substantial hotspots, especially for low values of (α/β)x. Variable RBE modeling often prompted dose constraints to be exceeded for critical structures.Conclusion: For base-of-skull and pediatric craniopharyngioma cases, both physical and biological robustness analyses should be considered for IMPT: these analyses can substantially affect the sparing of OARs and comparisons against VMAT. All proton RBE modeling is subject to high levels of uncertainty, but the clinical community should remain cognizant possible RBE effects. Careful clinical and imaging follow-up, plus further research on end-of-range RBE mitigation strategies such as LET optimization, should be prioritized for these cohorts of proton patients.

First benchmarking of the BIANCA model for cell survival prediction in a clinical hadron therapy scenario

In the framework of RBE modelling for hadron therapy, the BIANCA biophysical model was extended to O-ions and was used to construct a radiobiological database describing the survival of V79 cells as a function of ion type (1  ⩽  Z  ⩽  8) and energy. This database allowed performing RBE predictions in very good agreement with experimental data. A method was then developed to construct analogous databases for different cell lines, starting from the V79 database as a reference. Following interface to the FLUKA Monte Carlo radiation transport code, BIANCA was then applied for the first time to predict cell survival in a typical patient treatment scenario, consisting of two opposing fields of range-equivalent protons or C-ions. The model predictions were found to be in good agreement with CHO cell survival data obtained at the Heidelberg ion-beam therapy (HIT) centre, as well as predictions performed by the local effect model (version LEM IV). This work shows that BIANCA can be used to predict cell survival and RBE not only for V79 and AG01522 cells, as shown previously, but also, in principle, for any cell line of interest. Furthermore, following interface to a transport code like FLUKA, BIANCA can provide predictions of 3D biological dose distributions for hadron therapy treatments, thus laying the foundations for future applications in clinics.

Simulation of a radiobiology facility for the Centre for the Clinical Application of Particles

The Centre for the Clinical Application of Particles' Laser-hybrid Accelerator for Radiobiological Applications (LhARA) facility is being studied and requires simulation of novel accelerator components (such as the Gabor lens capture system), detector simulation and simulation of the ion beam interaction with cells. The first stage of LhARA will provide protons up to 15 MeV for in vitro studies. The second stage of LhARA will use a fixed-field accelerator to increase the energy of the particles to allow in vivo studies with protons and in vitro studies with heavier ions. BDSIM, a Geant4 based accelerator simulation tool, has been used to perform particle tracking simulations to verify the beam optics design done by BeamOptics and these show good agreement. Design parameters were defined based on an EPOCH simulation of the laser source and a series of mono-energetic input beams were generated from this by BDSIM. The tracking results show the large angular spread of the input beam (0.2 rad) can be transported with a transmission of almost 100% whilst keeping divergence at the end station very low (<0.1 mrad). The legacy of LhARA will be the demonstration of technologies that could drive a step-change in the provision of proton and light ion therapy (i.e. a laser source coupled to a Gabor lens capture and a fixed-field accelerator), and a system capable of delivering a comprehensive set of experimental data that can be used to enhance the clinical application of proton and light ion therapy.

Biophysical modeling and experimental validation of relative biological effectiveness (RBE) for 4He ion beam therapy

BACKGROUND

Helium (4He) ion beam therapy provides favorable biophysical characteristics compared to currently administered particle therapies, i.e., reduced lateral scattering and enhanced biological damage to deep-seated tumors like heavier ions, while simultaneously lessened particle fragmentation in distal healthy tissues as observed with lighter protons. Despite these biophysical advantages, raster-scanning 4He ion therapy remains poorly explored e.g., clinical translational is hampered by the lack of reliable and robust estimation of physical and radiobiological uncertainties. Therefore, prior to the upcoming 4He ion therapy program at the Heidelberg Ion-beam Therapy Center (HIT), we aimed to characterize the biophysical phenomena of 4He ion beams and various aspects of the associated models for clinical integration.

METHODS

Characterization of biological effect for 4He ion beams was performed in both homogenous and patient-like treatment scenarios using innovative models for estimation of relative biological effectiveness (RBE) in silico and their experimental validation using clonogenic cell survival as the gold-standard surrogate. Towards translation of RBE models in patients, the first GPU-based treatment planning system (non-commercial) for raster-scanning 4He ion beams was devised in-house (FRoG).

RESULTS

Our data indicate clinically relevant uncertainty of ±5-10% across different model simulations, highlighting their distinct biological and computational methodologies. The in vitro surrogate for highly radio-resistant tissues presented large RBE variability and uncertainty within the clinical dose range.

CONCLUSIONS

Existing phenomenological and mechanistic/biophysical models were successfully integrated and validated in both Monte Carlo and GPU-accelerated analytical platforms against in vitro experiments, and tested using pristine peaks and clinical fields in highly radio-resistant tissues where models exhibit the greatest RBE uncertainty. Together, these efforts mark an important step towards clinical translation of raster-scanning 4He ion beam therapy to the clinic.