A Monte Carlo based radiation response modelling framework to assess variability of clinical RBE in proton therapy

Double scattering and pencil beam scanning Monte Carlo workflows for proton therapy retrospective studies on radiation-induced toxicities

PURPOSE

Monte Carlo (MC) simulations can be used to accurately simulate dose and linear energy transfers (LET) distributions, thereby allowing for the calculation of the relative biological effectiveness (RBE) of protons. We present hereby the validation and implementation of a workflow for the Monte Carlo modelling of the double scattered and pencil beam scanning proton beamlines at our institution.

METHODS

The TOPAS/Geant4 MC model of the clinical nozzle has been comprehensively validated against measurements. The validation also included a comparison between simulated clinical treatment plans for four representative patients and the clinical treatment planning system (TPS). Moreover, an in-house tool implemented in Python was tested to assess the variable RBE-weighted dose in proton plans, which was illustrated for a patient case with a developing radiation-induced toxicity.

RESULTS

The simulated range and modulation width closely matches the measurements. Gamma-indexes (3%/3mm 3D), which compare the TPS and MC computations, showed a passing rate superior to 98%. The calculated RBE-weighted dose presented a slight increase at the necrosis location, within the PTV margins. This indicates the need for reporting on the physical and biological effects of irradiation in high dose regions, especially at the healthy tissues and increased LET distributions location.

CONCLUSION

The results demonstrate that the Monte Carlo method can be used to independently validate a TPS calculation, and to estimate LET distributions. The features of the in-house tool can be used to correlate LET and RBE-weighted dose distributions with the incidence of radiation-induced toxicities following proton therapy treatments.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Technical note: Investigation of dose and LETd effect to rectum and bladder by using non-straight laterals in prostate cancer receiving proton therapy

BACKGROUND

Parallel-opposed lateral beams are the conventional beam arrangements in proton therapy for prostate cancer. However, when considering linear energy transfer (LET) and RBE effects, alternative beam arrangements should be investigated.

PURPOSE

To investigate the dose and dose averaged LET (LET_d ) impact of using new beam arrangements rotating beams 5°-15° posteriorly to the laterals in prostate cancer treated with pencil-beam-scanning (PBS) proton therapy.

METHODS

Twenty patients with localized prostate cancer were included in this study. Four proton treatment plans for each patient were generated utilizing 0°, 5°, 10°, and 15° posterior oblique beam pairs relative to parallel-opposed lateral beams. Dose-volume histograms (DVHs) from posterior oblique beams were analyzed. Dose-LET_d -volume histogram (DLVH) was employed to study the difference in dose and LET_d with each beam arrangement. DLVH indices, V ( d , l ) $V( {d,l} )$ , defined as the cumulative absolute volume that has a dose of at least d (Gy[RBE]) and a LET_d of at least l (keV/µm), were calculated for both the rectum and bladder to the whole group of patients and two-sub groups with and without hydrogel spacer. These metrics were tested using Wilcoxon signed-rank test.

RESULTS

Rotating beam angles from laterals to slightly posterior by 5°-15° reduced high LET_d volumes while it increased the dose volume in the rectum and increased LET_d in bladders. Beam angles rotated five degrees posteriorly from laterals (i.e., gantry in 95° and 265°) are proposed since they achieved the optimal balance of better LET_d sparing and minimal dose increase in the rectum. A reduction of V(50 Gy[RBE], 2.6 keV/µm) from 7.41 to 3.96 cc (p < 0.01), and a slight increase of V(50 Gy[RBE], 0 keV/µm) from 20.1 to 21.6 cc (p < 0.01) were observed for the group without hydrogel spacer. The LET_d sparing was less effective for the group with hydrogel spacer, which achieved the reduction of V(50 Gy[RBE], 2.6 keV/µm) from 4.28 to 2.10 cc (p < 0.01).

CONCLUSIONS

Posterior oblique angle plans improved LET_d sparing of the rectum while sacrificing LET_d sparing in the bladder in the treatment of prostate cancer with PBS. Beam angle modification from laterals to slightly posterior may be a strategy to redistribute LET_d and perhaps reduce rectal toxicity risks in prostate cancer patients treated with PBS. However, the effect is reduced for patients with hydrogel spacer.

A systematic review of clinical studies on variable proton Relative Biological Effectiveness (RBE)

Recently, a number of clinical studies have explored links between possible Relative Biological Effectiveness (RBE) elevations and patient toxicities and/or image changes following proton therapy. Our objective was to perform a systematic review of such studies. We applied a "Problem [RBE], Intervention [Protons], Population [Patients], Outcome [Side effect]" search strategy to the PubMed database. From our search, we retrieved studies which: (a) performed novel voxel-wise analyses of patient effects versus physical dose and LET (n = 13), and (b) compared image changes between proton and photon cohorts with regard to proton RBE (n = 9). For each retrieved study, we extracted data regarding: primary tumour type; size of patient cohort; type of image change studied; image-registration method (deformable or rigid); LET calculation method, and statistical methodology. We compared and contrasted their methods in order to discuss the weight of clinical evidence for variable proton RBE. We concluded that clinical evidence for variable proton RBE remains statistically weak at present. Our principal recommendation is that proton centres and clinical trial teams collaborate to standardize follow-up protocols and statistical analysis methods, so that larger patient cohorts can ultimately be considered for RBE analyses.

Exploratory study of seed spots analysis to characterize dose and linear energy transfer effect in adverse event initialization of pencil beam scanning proton therapy

BACKGROUND

Both dose and linear-energy-transfer (LET) could play a substantial role in adverse event (AE) initialization of cancer patients treated with pencil-beam-scanning proton therapy (PBS). However, not all the voxels within the AE regions are directly induced from the dose and LET effect. It is important to study the synergistic effect of dose and LET in AE initialization by only including a subset of voxels that are dosimetrically important.

PURPOSE

To perform exploratory investigation of the dose and LET effects upon AE initialization in PBS using seed spots analysis.

METHODS

113 head and neck (H&N) cancer patients receiving curative PBS were included. Among them, 20 patients experienced unanticipated CTCAEv4.0 grade≥3 AEs (AE group) and 93 patients did not (control group). Within the AE group, 13 AE patients were included in the seed spot analysis to derive the descriptive features of AE initialization and the remaining 7 mandible osteoradionecrosis patients and 93 control patients were used to derive the feature-based volume constraint of mandible osteoradionecrosis. The AE regions were contoured and the corresponding dose-LET volume histograms (DLVHs) of AE regions were generated for all patients in the AE group. We selected high LET voxels (the highest 5% of each dose bin) with a range of moderate to high dose (≥∼40 Gy[RBE]) as critical voxels. Critical voxels which were contiguous with each other were grouped into clusters. Each cluster was considered as a potential independent seed spot for AE initialization. Seed spots were displayed in a 2D dose-LET plane based on their mean dose and LET to derive the descriptive features of AE initialization. A volume constraint of mandible osteoradionecrosis was then established based on the extracted features using a receiver operating characteristic curve.

RESULTS

The product of dose and LET (xBD) was found to be a descriptive feature of seed spots leading to AE initialization in this preliminary study. The derived xBD volume constraint for mandible osteoradionecrosis showed good performance with an area-under-curve of 0.87 (sensitivity of 0.714 and specificity of 0.807 in the leave-one-out cross validation) for the very limited patient data included in this study.

CONCLUSION

Our exploratory study showed that both dose and LET were observed to be important in AE initializations. The derived xBD volume constraint could predict mandible osteoradionecrosis reasonably well in the very limited H&N cancer patient data treated with PBS included in this study. This article is protected by copyright. All rights reserved.

Towards harmonizing clinical linear energy transfer (LET) reporting in proton radiotherapy: a European multi-centric study

BACKGROUND

Clinical data suggest that the relative biological effectiveness (RBE) in proton therapy (PT) varies with linear energy transfer (LET). However, LET calculations are neither standardized nor available in clinical routine. Here, the status of LET calculations among European PT institutions and their comparability are assessed.

MATERIALS AND METHODS

Eight European PT institutions used suitable treatment planning systems with their center-specific beam model to create treatment plans in a water phantom covering different field arrangements and fulfilling commonly agreed dose objectives. They employed their locally established LET simulation environments and procedures to determine the corresponding LET distributions. Dose distributions D1.1 and DRBE assuming constant and variable RBE, respectively, and LET were compared among the institutions. Inter-center variability was assessed based on dose- and LET-volume-histogram parameters.

RESULTS

Treatment plans from six institutions fulfilled all clinical goals and were eligible for common analysis. D1.1 distributions in the target volume were comparable among PT institutions. However, corresponding LET values varied substantially between institutions for all field arrangements, primarily due to differences in LET averaging technique and considered secondary particle spectra. Consequently, DRBE using non-harmonized LET calculations increased inter-center dose variations substantially compared to D1.1 and significantly in mean dose to the target volume of perpendicular and opposing field arrangements (p < 0.05). Harmonizing LET reporting (dose-averaging, all protons, LET to water or to unit density tissue) reduced the inter-center variability in LET to the order of 10-15% within and outside the target volume for all beam arrangements. Consequentially, inter-institutional variability in DRBE decreased to that observed for D1.1.

CONCLUSION

Harmonizing the reported LET among PT centers is feasible and allows for consistent multi-centric analysis and reporting of tumor control and toxicity in view of a variable RBE. It may serve as basis for harmonized variable RBE dose prescription in PT.

A Monte Carlo study of different LET definitions and calculation parameters for proton beam therapy

The strongin vitroevidence that proton Relative Biological Effectiveness (RBE) varies with Linear Energy Transfer (LET) has led to an interest in applying LET within treatment planning. However, there is a lack of consensus on LET definition, Monte Carlo (MC) parameters or clinical methodology. This work aims to investigate how common variations of LET definition may affect potential clinical applications. MC simulations (GATE/GEANT4) were used to calculate absorbed dose and different types of LET for a simple Spread Out Bragg Peak (SOBP) and for four clinical PBT plans covering a range of tumour sites. Variations in the following LET calculation methods were considered: (i) averaging (dose-averaged LET (LET_d) & track-averaged LET); (ii) scoring (LET_dto water, to medium and to mass density); (iii) particle inclusion (LET_dto all protons, to primary protons and to particles); (iv) MC settings (hit type and Maximum Step Size (MSS)). LET distributions were compared using: qualitative comparison, LET Volume Histograms (LVHs), single value criteria (maximum and mean values) and optimised LET-weighted dose models. Substantial differences were found between LET values in averaging, scoring and particle type. These differences depended on the methodology, but for one patient a difference of ∼100% was observed between the maximum LET_dfor all particles and maximum LET_dfor all protons within the brainstem in the high isodose region (4 keVμm^-1and 8 keVμm^-1respectively). An RBE model using LET_dincluding heavier ions was found to predict substantially different LET-weighted dose compared to those using other LET definitions. In conclusion, the selection of LET definition may affect the results of clinical metrics considered in treatment planning and the results of an RBE model. The authors' advocate for the scoring of dose-averaged LET to water for primary and secondary protons using a random hit type and automated MSS.

Per-voxel constraints to minimize hot spots in linear energy transfer-guided robust optimization for base of skull head and neck cancer patients in IMPT

PURPOSE

Due to the employment of quadratic programming using soft constraints to implement dose volume constraints and the "trial-and-error" procedure needed to achieve a clinically acceptable plan, conventional dose volume constraints (upper limit) are not adequately effective in controlling small and isolated hot spots in the dose/linear energy transfer (LET) distribution. Such hot spots can lead to adverse events. In order to mitigate the risk of brain necrosis, one of the most clinically significant adverse events in patients receiving intensity-modulated proton therapy (IMPT) for base of skull (BOS) cancer, we propose per-voxel constraints to minimize hot spots in LET-guided robust optimization.

METHODS AND MATERIALS

Ten BOS cancer patients treated with IMPT were carefully selected by meeting one of the following conditions: (1) diagnosis of brain necrosis during follow-up; and (2) considered high risk for brain necrosis by not meeting dose constraints to the brain. An optimizing structure (BrainOPT) and an evaluating structure (BrainROI) that both contained the aforementioned hot dose regions in the brain were generated for optimization and evaluation, respectively. Two plans were generated for every patient: one using conventional dose-only robust optimization, the other using LET-guided robust optimization. The impact of LET was integrated into the optimization via a term of extra biological dose (xBD). A novel optimization tool of per-voxel constraints to control small and isolated hot spots in either the dose, LET, or combined (dose/LET) distribution was developed and used to minimize dose/LET hot spots of the selected structures. Indices from dose-volume histogram (DVH) and xBD dose-volume histogram (xBDVH) were used in the plan evaluation. A newly developed tool of the dose-LET-volume histogram (DLVH) was also adopted to illustrate the underlying mechanism. Wilcoxon signed-rank test was used for statistical comparison of the DVH and xBDVH indices between the conventional dose-only and the LET-guided robustly optimized plans.

RESULTS

Per-voxel constraints effectively and efficiently minimized dose hot spots in both dose-only and LET-guided robust optimization and LET hot spots in LET-guided robust optimization. Compared to the conventional dose-only robust optimization, the LET-guided robust optimization could generate plans with statistically lower xBD hot spots in BrainROI (VxBD,50 Gy[RBE], p = 0.009; VxBD,60 Gy[RBE], p = 0.025; xBD1cc, p = 0.017; xBD2cc, p = 0.022) with comparable dose coverage, dose hot spots in the target, and dose hot spots in BrainROI. DLVH analysis indicated that LET-guided robust optimization could either reduce LET at the same dose level or redistribute high LET from high dose regions to low dose regions.

CONCLUSION

Per-voxel constraint is a powerful tool to minimize dose/LET hot spots in IMPT. The LET-guided robustly optimized plans outperformed the conventional dose-only robustly optimized plans in terms of xBD hot spots control.

Models for Translational Proton Radiobiology-From Bench to Bedside and Back

The number of proton therapy centers worldwide are increasing steadily, with more than two million cancer patients treated so far. Despite this development, pending questions on proton radiobiology still call for basic and translational preclinical research. Open issues are the on-going discussion on an energy-dependent varying proton RBE (relative biological effectiveness), a better characterization of normal tissue side effects and combination treatments with drugs originally developed for photon therapy. At the same time, novel possibilities arise, such as radioimmunotherapy, and new proton therapy schemata, such as FLASH irradiation and proton mini-beams. The study of those aspects demands for radiobiological models at different stages along the translational chain, allowing the investigation of mechanisms from the molecular level to whole organisms. Focusing on the challenges and specifics of proton research, this review summarizes the different available models, ranging from in vitro systems to animal studies of increasing complexity as well as complementing in silico approaches.

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.

Mixed Effect Modeling of Dose and Linear Energy Transfer Correlations With Brain Image Changes After Intensity Modulated Proton Therapy for Skull Base Head and Neck Cancer

Purpose:

Intensity modulated proton therapy (IMPT) could yield high linear energy transfer (LET) in critical structures and increased biological effect. For head and neck cancers at the skull base this could potentially result in radiation-associated brain image change (RAIC). The purpose of the current study was to investigate voxel-wise dose and LET correlations with RAIC after IMPT.

Methods and Materials:

For 15 patients with RAIC after IMPT, contrast enhancement observed on T1-weighted magnetic resonance imaging was contoured and coregistered to the planning computed tomography. Monte Carlo calculated dose and dose-averaged LET (LET_d) distributions were extracted at voxel level and associations with RAIC were modelled using uni- and multivariate mixed effect logistic regression. Model performance was evaluated using the area under the receiver operating characteristic curve and precision-recall curve.

Results:

An overall statistically significant RAIC association with dose and LET_d was found in both the uni- and multivariate analysis. Patient heterogeneity was considerable, with standard deviation of the random effects of 1.81 (1.30–2.72) for dose and 2.68 (1.93–4.93) for LET_d, respectively. Area under the receiver operating characteristic curve was 0.93 and 0.95 for the univariate dose-response model and multivariate model, respectively. Analysis of the LET_d effect demonstrated increased risk of RAIC with increasing LET_d for the majority of patients. Estimated probability of RAIC with LET_d = 1 keV/μm was 4% (95% confidence interval, 0%, 0.44%) and 29% (95% confidence interval, 0.01%, 0.92%) for 60 and 70 Gy, respectively. The TD₁₅ were estimated to be 63.6 and 50.1 Gy with LET_d equal to 2 and 5 keV/μm, respectively.

Conclusions:

Our results suggest that the LET_d effect could be of clinical significance for some patients; LET_d assessment in clinical treatment plans should therefore be taken into consideration.

A systematic review on the usage of averaged LET in radiation biology for particle therapy

Linear Energy Transfer (LET) is widely used to express the radiation quality of ion beams, when characterizing the biological effectiveness. However, averaged LET may be defined in multiple ways, and the chosen definition may impact the resulting reported value. We review averaged LET definitions found in the literature, and quantify which impact using these various definitions have for different reference setups. We recorded the averaged LET definitions used in 354 publications quantifying the relative biological effectiveness (RBE) of hadronic beams, and investigated how these various definitions impact the reported averaged LET using a Monte Carlo particle transport code. We find that the kind of averaged LET being applied is, generally, poorly defined. Some definitions of averaged LET may influence the reported averaged LET values up to an order of magnitude. For publications involving protons, most applied dose averaged LET when reporting RBE. The absence of what target medium is used and what secondary particles are included further contributes to an ill-defined averaged LET. We also found evidence of inconsistent usage of averaged LET definitions when deriving LET-based RBE models. To conclude, due to commonly ill-defined averaged LET and to the inherent problems of LET-based RBE models, averaged LET may only be used as a coarse indicator of radiation quality. We propose a more rigorous way of reporting LET values, and suggest that ideally the entire particle fluence spectra should be recorded and provided for future RBE studies, from which any type of averaged LET (or other quantities) may be inferred.

Late Side Effects in Normal Mouse Brain Tissue After Proton Irradiation

Radiation-induced late side effects such as cognitive decline and normal tissue complications can severely affect quality of life and outcome in long-term survivors of brain tumors. Proton therapy offers a favorable depth-dose deposition with the potential to spare tumor-surrounding normal tissue, thus potentially reducing such side effects. In this study, we describe a preclinical model to reveal underlying biological mechanisms caused by precise high-dose proton irradiation of a brain subvolume. We studied the dose- and time-dependent radiation response of mouse brain tissue, using a high-precision image-guided proton irradiation setup for small animals established at the University Proton Therapy Dresden (UPTD). The right hippocampal area of ten C57BL/6 and ten C3H/He mice was irradiated. Both strains contained four groups (n_irradiated = 3, n_control = 1) treated with increasing doses (0 Gy, 45 Gy, 65 Gy or 85 Gy and 0 Gy, 40 Gy, 60 Gy or 80 Gy, respectively). Follow-up examinations were performed for up to six months, including longitudinal monitoring of general health status and regular contrast-enhanced magnetic resonance imaging (MRI) of mouse brains. These findings were related to comprehensive histological analysis. In all mice of the highest dose group, first symptoms of blood-brain barrier (BBB) damage appeared one week after irradiation, while a dose-dependent delay in onset was observed for lower doses. MRI contrast agent leakage occurred in the irradiated brain areas and was progressive in the higher dose groups. Mouse health status and survival corresponded to the extent of contrast agent leakage. Histological analysis revealed tissue changes such as vessel abnormalities, gliosis, and granule cell dispersion, which also partly affected the non-irradiated contralateral hippocampus in the higher dose groups. All observed effects depended strongly on the prescribed radiation dose and the outcome, i.e. survival, image changes, and tissue alterations, were very consistent within an experimental dose cohort. The derived dose–response model will determine endpoint-specific dose levels for future experiments and may support generating clinical hypotheses on brain toxicity after proton therapy.

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.

Late Contrast Enhancing Brain Lesions in Proton-Treated Patients With Low-Grade Glioma: Clinical Evidence for Increased Periventricular Sensitivity and Variable RBE

PURPOSE

Late radiation-induced contrast-enhancing brain lesions (CEBLs) on magnetic resonance imaging (MRI) after proton therapy of brain tumors have been observed to occur frequently in regions of high linear energy transfer (LET) and in proximity to the ventricular system. We analyzed 110 patients with low-grade glioma treated with proton therapy to determine whether the risk for CEBLs is increased in proximity to the ventricular system and if there is a relationship between relative biological effectiveness (RBE) and LET.

METHODS AND MATERIALS

We contoured CEBLs identified on follow-up T1-MRI scans and computed dose and dose-averaged LET (LET_d) distributions for all patients using the Monte Carlo method. We then performed cross-validated voxel-level logistic regression to predict local risks for image change and to extract model parameters, such as the RBE. From the voxel-level model, we derived a model for patient-level risk prediction based on the treatment plan.

RESULTS

Of 110 patients, 23 exhibited 1 or several CEBLs on follow-up MRI scans. The voxel-level logistic model has an accuracy as follows: area under the curve of 0.94 and Brier score of 2.6 × 10^-5. Model predictions are a 3-fold increased risk in the 4 mm region around the ventricular system and an LET_d-dependent RBE of, for example, 1.20 for LET_d = 2 keV/μm and 1.50 for LET_d = 5 keV/μm. The patient-level risk model has an accuracy as follows: area under the curve of 0.78 and Brier score of 0.13.

CONCLUSIONS

Our findings present clinical evidence for an increased risk in ventricular proximity and for a proton RBE that increases significantly with increasing LET. We present a voxel-level model that accurately predicts the localization of late MRI contrast change and extrapolate a patient-level model that allows treatment plan-based risk prediction.

High-precision image-guided proton irradiation of mouse brain sub-volumes

BACKGROUND AND PURPOSE

Proton radiotherapy offers the potential to reduce normal tissue toxicity. However, clinical safety margins, range uncertainties, and varying relative biological effectiveness (RBE) may result in a critical dose in tumor-surrounding normal tissue. To assess potential adverse effects in preclinical studies, image-guided proton mouse brain irradiation and analysis of DNA damage repair was established.

MATERIAL AND METHODS

We designed and characterized a setup to shape proton beams with 7 mm range in water and 3 mm in diameter and commissioned a Monte Carlo model for in vivo dose simulation. Cone-beam computed tomography and orthogonal X-ray imaging were used to delineate the right hippocampus and position the mice. The brains of three C3H/HeNRj mice were irradiated with 8 Gy and excised 30 min later. Initial DNA double-strand breaks were visualized by staining brain sections for cell nuclei and γH2AX. Imaged sections were analyzed with an automated and validated processing pipeline to provide a quantitative, spatially resolved radiation damage indicator.

RESULTS

The analyzed DNA damage pattern clearly visualized the radiation effect in the mouse brains and could be mapped to the simulated dose distribution. The proton beam passed the right hippocampus and stopped in the central brain region for all evaluated mice.

CONCLUSION

We established image-guided proton irradiation of mouse brains. The clinically oriented workflow facilitates (back-) translational studies. Geometric accuracy, detailed Monte Carlo dose simulations, and cell-based assessment enable a biologically and spatially resolved analysis of radiation response and RBE.

Reduced radiation-induced toxicity by using proton therapy for the treatment of oropharyngeal cancer

Patients with squamous cell carcinoma of the oropharynx are generally treated with (chemo) radiation. Patients with oropharyngeal cancer have better survival than patients with squamous cell carcinoma of other head and neck subsites, especially when related to human papillomavirus. However, radiotherapy results in a substantial percentage of survivors suffering from significant treatment-related side-effects. Late radiation-induced side-effects are mostly irreversible and may even be progressive, and particularly xerostomia and dysphagia affect health-related quality of life. As the risk of radiation-induced side-effects highly depends on dose to healthy normal tissues, prevention of radiation-induced xerostomia and dysphagia and subsequent improvement of health-relatedquality of life can be obtained by applying proton therapy, which offers the opportunity to reduce the dose to both the salivary glands and anatomic structures involved in swallowing.This review describes the results of the first cohort studies demonstrating that proton therapy results in lower dose levels in multiple organs at risk, which translates into reduced acute toxicity (i.e. up to 3 months after radiotherapy), while preserving tumour control. Next to reducing mucositis, tube feeding, xerostomia and distortion of the sense of taste, protons can improve general well-being by decreasing fatigue and nausea. Proton therapy results in decreased rates of tube feeding dependency and severe weight loss up to 1 year after radiotherapy, and may decrease the risk of radionecrosis of the mandible. Also, the model-based approach for selecting patients for proton therapy in the Netherlands is described in this review and future perspectives are discussed.