Predicting late magnetic resonance image changes in glioma patients after proton therapy

Deep learning based linear energy transfer calculation for proton therapy

Objective.This study aims to address the limitations of traditional methods for calculating linear energy transfer (LET), a critical component in assessing relative biological effectiveness (RBE). Currently, Monte Carlo (MC) simulation, the gold-standard for accuracy, is resource-intensive and slow for dose optimization, while the speedier analytical approximation has compromised accuracy. Our objective was to prototype a deep-learning-based model for calculating dose-averaged LET (LETd) using patient anatomy and dose-to-water (DW) data, facilitating real-time biological dose evaluation and LET optimization within proton treatment planning systems.Approach. 275 4-field prostate proton Stereotactic Body Radiotherapy plans were analyzed, rendering a total of 1100 fields. Those were randomly split into 880, 110, and 110 fields for training, validation, and testing. A 3D Cascaded UNet model, along with data processing and inference pipelines, was developed to generate patient-specific LETddistributions from CT images and DW. The accuracy of the LETdof the test dataset was evaluated against MC-generated ground truth through voxel-based mean absolute error (MAE) and gamma analysis.Main results.The proposed model accurately inferred LETddistributions for each proton field in the test dataset. A single-field LETdcalculation took around 100 ms with trained models running on a NVidia A100 GPU. The selected model yielded an average MAE of 0.94 ± 0.14 MeV cm-1and a gamma passing rate of 97.4% ± 1.3% when applied to the test dataset, with the largest discrepancy at the edge of fields where the dose gradient was the largest and counting statistics was the lowest.Significance.This study demonstrates that deep-learning-based models can efficiently calculate LETdwith high accuracy as a fast-forward approach. The model shows great potential to be utilized for optimizing the RBE of proton treatment plans. Future efforts will focus on enhancing the model's performance and evaluating its adaptability to different clinical scenarios.

Dosimetric comparison of advanced radiation techniques for scalp-sparing in low-grade gliomas

BACKGROUND

Alopecia causes significant distress for patients and negatively impacts quality of life for low-grade glioma (LGG) patients. We aimed to compare and evaluate variations in dose distribution for scalp-sparing in LGG patients with proton therapy and photon therapy, namely intensity-modulated proton therapy (IMPT), intensity-modulated radiotherapy (IMRT), volumetric modulated arc therapy (VMAT), and helical tomotherapy (HT).

METHODS

This retrospective study utilized a dataset comprising imaging data from 22 patients with LGG who underwent postoperative radiotherapy. Treatment plans were generated for each patient with scalp-optimized (SO) approaches and scalp-non-optimized (SNO) approaches using proton techniques and photons techniques; all plans adhered to the same dose constraint of delivering a total radiation dose of 54.04 Gy to the target volume. All treatment plans were subsequently analyzed.

RESULTS

All the plans generated in this study met the dose constraints for the target volume and OARs. The SO plans resulted in reduced maximum scalp dose (Dmax), mean scalp dose (Dmean), and volume of the scalp receiving 30 Gy (V30) and 40 Gy (V40) compared with SNO plans in all radiation techniques. Among all radiation techniques, the IMPT plans exhibited superior performance compared to other plans for dose homogeneity as for SO plans. Also, IMPT showed lower values for Dmean and Dmax than all photon radiation techniques.

CONCLUSION

Our study provides evidence that the SO approach is a feasible technique for reducing scalp radiation dose. However, it is imperative to conduct prospective trials to assess the benefits associated with this approach.

Proton therapy induces a local microglial neuroimmune response

BACKGROUND AND PURPOSE

Although proton therapy is increasingly being used in the treatment of paediatric and adult brain tumours, there are still uncertainties surrounding the biological effect of protons on the normal brain. Microglia, the brain-resident macrophages, have been shown to play a role in the development of radiation-induced neurotoxicity. However, their molecular and hence functional response to proton irradiation remains unknown. This study investigates the effect of protons on microglia by comparing the effect of photons and protons as well as the influence of age and different irradiated volumes.

MATERIALS AND METHODS

Rats were irradiated with 14 Gy to the whole brain with photons (X-rays), plateau protons, spread-out Bragg peak (SOBP) protons or to 50 % anterior, or 50 % posterior brain sub-volumes with plateau protons. RNA sequencing, validation of microglial priming gene expression using qPCR and high-content imaging analysis of microglial morphology were performed in the cortex at 12 weeks post irradiation.

RESULTS

Photons and plateau protons induced a shared transcriptomic response associated with neuroinflammation. This response was associated with a similar microglial priming gene expression signature and distribution of microglial morphologies. Expression of the priming gene signature was less pronounced in juvenile rats compared to adults and slightly increased in rats irradiated with SOBP protons. High-precision partial brain irradiation with protons induced a local microglial priming response and morphological changes.

CONCLUSION

Overall, our data indicate that the brain responds in a similar manner to photons and plateau protons with a shared local upregulation of microglial priming-associated genes, potentially enhancing the immune response to subsequent inflammatory challenges.

The dirty and clean dose concept: Towards creating proton therapy treatment plans with a photon-like dose response

BACKGROUND

Applying tolerance doses for organs at risk (OAR) from photon therapy introduces uncertainties in proton therapy when assuming a constant relative biological effectiveness (RBE) of 1.1.

PURPOSE

This work introduces the novel dirty and clean dose concept, which allows for creating treatment plans with a more photon-like dose response for OAR and, thus, less uncertainties when applying photon-based tolerance doses.

METHODS

The concept divides the 1.1-weighted dose distribution into two parts: the clean and the dirty dose. The clean and dirty dose are deposited by protons with a linear energy transfer (LET) below and above a set LET threshold, respectively. For the former, a photon-like dose response is assumed, while for the latter, the RBE might exceed 1.1. To reduce the dirty dose in OAR, a MaxDirtyDose objective was added in treatment plan optimization. It requires setting two parameters: LET threshold and max dirty dose level. A simple geometry consisting of one target volume and one OAR in water was used to study the reduction in dirty dose in the OAR depending on the choice of the two MaxDirtyDose objective parameters during plan optimization. The best performing parameter combinations were used to create multiple dirty dose optimized (DDopt) treatment plans for two cranial patient cases. For each DDopt plan, 1.1-weighted dose, variable RBE-weighted dose using the Wedenberg RBE model and dose-average LETd distributions as well as resulting normal tissue complication probability (NTCP) values were calculated and compared to the reference plan (RefPlan) without MaxDirtyDose objectives.

RESULTS

In the water phantom studies, LET thresholds between 1.5 and 2.5 keV/µm yielded the best plans and were subsequently used. For the patient cases, nearly all DDopt plans led to a reduced Wedenberg dose in critical OAR. This reduction resulted from an LET reduction and translated into an NTCP reduction of up to 19 percentage points compared to the RefPlan. The 1.1-weighted dose in the OARs was slightly increased (patient 1: 0.45 Gy(RBE), patient 2: 0.08 Gy(RBE)), but never exceeded clinical tolerance doses. Additionally, slightly increased 1.1-weighted dose in healthy brain tissue was observed (patient 1: 0.81 Gy(RBE), patient 2: 0.53 Gy(RBE)). The variation of NTCP values due to variation of α/β from 2 to 3 Gy was much smaller for DDopt (2 percentage points (pp)) than for RefPlans (5 pp).

CONCLUSIONS

The novel dirty and clean dose concept allows for creating biologically more robust proton treatment plans with a more photon-like dose response. The reduced uncertainties in RBE can, therefore, mitigate uncertainties introduced by using photon-based tolerance doses for OAR.

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

Introduction

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

Methods

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

Results and discussion

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

Linear energy transfer-inclusive models of brainstem necrosis following proton therapy of paediatric ependymoma

Background and Purpose

Radiation-induced brainstem necrosis after proton therapy is a severe toxicity with potential association to uncertainties in the proton relative biological effectiveness (RBE). A constant RBE of 1.1 is assumed clinically, but the RBE is known to vary with linear energy transfer (LET). LET-inclusive predictive models of toxicity may therefore be beneficial during proton treatment planning. Hence, we aimed to construct models describing the association between brainstem necrosis and LET in the brainstem.

Materials and methods

A matched case-control cohort (n = 28, 1:3 case-control ratio) of symptomatic brainstem necrosis was selected from 954 paediatric ependymoma brain tumour patients treated with passively scattered proton therapy. Dose-averaged LET (LET_d) parameters in restricted volumes (L_50%, L_10% and L_0.1cm³, the cumulative LET_d) within high-dose thresholds were included in linear- and logistic regression normal tissue complication probability (NTCP) models.

Results

A 1 keV/µm increase in L_10% to the brainstem volume receiving dose over 54 Gy(RBE) led to an increased brainstem necrosis risk [95% confidence interval] of 2.5 [0.0, 7.8] percentage points. The corresponding logistic regression model had area under the receiver operating characteristic curve (AUC) of 0.76, increasing to 0.84 with the anterior pons substructure as a second parameter. 19 [7, 350] patients with toxicity were required to associate the L_10% (D > 54 Gy(RBE)) and brainstem necrosis with 80% statistical power.

Conclusion

The established models of brainstem necrosis illustrate a potential impact of high LET regions in patients receiving high doses to the brainstem, and thereby support LET mitigation during clinical treatment planning.

Slice2Volume: Fusion of multimodal medical imaging and light microscopy data of irradiation-injured brain tissue in 3D

Comprehending cellular changes of radiation-induced brain injury is crucial to prevent and treat the pathology. We provide a unique open dataset of proton-irradiated mouse brains consisting of medical imaging, radiation dose simulations, and large-scale microscopy images, all registered into a common coordinate system. This allows dose-dependent analyses on single-cell level.

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.

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

Purpose

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

Methods and Materials

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

Results

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

Conclusions

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

Comparing biological effectiveness guided plan optimization strategies for cranial proton therapy: potential and challenges

BACKGROUND

To introduce and compare multiple biological effectiveness guided (BG) proton plan optimization strategies minimizing variable relative biological effectiveness (RBE) induced dose burden in organs at risk (OAR) while maintaining plan quality with a constant RBE.

METHODS

Dose-optimized (DOSEopt) proton pencil beam scanning reference treatment plans were generated for ten cranial patients with prescription doses ≥ 54 Gy(RBE) and ≥ 1 OAR close to the clinical target volume (CTV). For each patient, four additional BG plans were created. BG objectives minimized either proton track-ends, dose-averaged linear energy transfer (LET_d), energy depositions from high-LET protons or variable RBE-weighted dose (D_RBE) in adjacent serially structured OARs. Plan quality (RBE = 1.1) was assessed by CTV dose coverage and robustness (2 mm setup, 3.5% density), dose homogeneity and conformity in the planning target volumes and adherence to OAR tolerance doses. LET_d, D_RBE (Wedenberg model, α/β_CTV = 10 Gy, α/β_OAR = 2 Gy) and resulting normal tissue complication probabilities (NTCPs) for blindness and brainstem necrosis were derived. Differences between DOSEopt and BG optimized plans were assessed and statistically tested (Wilcoxon signed rank, α = 0.05).

RESULTS

All plans were clinically acceptable. DOSEopt and BG optimized plans were comparable in target volume coverage, homogeneity and conformity. For recalculated D_RBE in all patients, all BG plans significantly reduced near-maximum D_RBE to critical OARs with differences up to 8.2 Gy(RBE) (p < 0.05). Direct D_RBE optimization primarily reduced absorbed dose in OARs (average ΔD_mean = 2.0 Gy; average ΔLET_d,mean = 0.1 keV/µm), while the other strategies reduced LET_d (average ΔD_mean < 0.3 Gy; average ΔLET_d,mean = 0.5 keV/µm). LET-optimizing strategies were more robust against range and setup uncertaintes for high-dose CTVs than D_RBE optimization. All BG strategies reduced NTCP for brainstem necrosis and blindness on average by 47% with average and maximum reductions of 5.4 and 18.4 percentage points, respectively.

CONCLUSIONS

All BG strategies reduced variable RBE-induced NTCPs to OARs. Reducing LET_d in high-dose voxels may be favourable due to its adherence to current dose reporting and maintenance of clinical plan quality and the availability of reported LET_d and dose levels from clinical toxicity reports after cranial proton therapy. These optimization strategies beyond dose may be a first step towards safely translating variable RBE optimization in the clinics.

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.

Iatrogenic Influence on Prognosis of Radiation-Induced Contrast Enhancements in Patients with Glioma WHO 1-3 following Photon and Proton Radiotherapy

BACKGROUND AND PURPOSE

Radiation-induced contrast enhancement (RICE) is a common side effect following radiotherapy for glioma, but both diagnosis and handling are challenging. Due to the potential risks associated with RICE and its challenges in differentiating RICE from tumor progression, it is critical to better understand how RICE prognosis depends on iatrogenic influence.

MATERIALS AND METHODS

We identified 99 patients diagnosed with RICE who were previously treated with either photon or proton therapy for World Health Organization (WHO) grade 1-3 primary gliomas. Post-treatment brain MRI-based volumetric analysis and clinical data collection was performed at multiple time points.

RESULTS

The most common histologic subtypes were astrocytoma (50%) and oligodendroglioma (46%). In 67%, it was graded WHO grade 2 and in 86% an IDH mutation was present. RICE first occurred after 16 months (range: 1 - 160) in median. At initial RICE occurrence, 39% were misinterpreted as tumor progression. A tumor-specific therapy including chemotherapy or re-irradiation led to a RICE size progression in 86% and 92% of cases, respectively and RICE symptom progression in 57% and 65% of cases, respectively. A RICE-specific therapy such as corticosteroids or Bevacizumab for larger or symptomatic RICE led to a RICE size regression in 81% of cases with symptom stability or regression in 62% of cases.

CONCLUSIONS

While with chemotherapy and re-irradiation a RICE progression was frequently observed, anti-edematous or anti-VEGF treatment frequently went along with a RICE regression. For RICE, correct diagnosis and treatment decisions are challenging and critical and should be made interdisciplinarily.

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.

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.

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

BACKGROUND AND PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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.

An ion-independent phenomenological relative biological effectiveness (RBE) model for proton therapy

BACKGROUND

A relative biological effectiveness (RBE) of 1.1 is used for proton therapy though clinical evidence of varying RBE was raised. Clinical studies on RBE variability have been conducted for decades for carbon radiation, which could advance the understanding of the clinical proton RBE given an ion-independent RBE model. In this work, such a model, linear and simple, using the beam quantity Q = Z²/E (Z = ion charge, E = kinetic energy per nucleon) was tested and compared to the commonly used, proton-specific and linear energy transfer (LET) based Wedenberg RBE model.

MATERIAL AND METHODS

The Wedenberg and Q models, both predicting RBE_max and RBE_min (i.e., RBE at vanishing and very high dose, respectively), are compared in terms of ion-dependence and prediction power. An experimental in-vitro data ensemble covering 115 publications for various ions was used as dataset.

RESULTS

The model parameter of the Q model was observed to be similar for different ions (in contrast to LET). The Q model was trained without any prior knowledge of proton data. For proton RBE, the differences between experimental data and corresponding predictions of the Wedenberg or the Q model were highly comparable.

CONCLUSIONS

A simple linear RBE model using Q instead of LET was proposed and tested to be able to predict proton RBE using model parameter trained based on only RBE data of other particles in a clinical proton energy range for a large in-vitro dataset. Adding (pre)clinical knowledge from carbon ion therapy may, therefore, reduce the dominating biological uncertainty in proton RBE modelling. This would translate in reduced RBE related uncertainty in proton therapy treatment planning.

Clinical use and future requirements of relative biological effectiveness: survey among all european proton therapy centres

BACKGROUND AND PURPOSE

The relative biological effectiveness (RBE) varies along the treatment field. However, in clinical practice, a constant RBE of 1.1 is assumed, which can result in undesirable side effects. This study provides an accurate overview of current clinical practice for considering proton RBE in Europe.

MATERIALS AND METHODS

A survey was devised and sent to all proton therapy centres in Europe that treat patients. The online questionnaire consisted of 39 questions addressing various aspects of RBE consideration in clinical practice, including treatment planning, patient follow-up and future demands.

RESULTS

All 25 proton therapy centres responded. All centres prescribed a constant RBE of 1.1, but also applied measures (except for one eye treatment centre) to counteract variable RBE effects such as avoiding beams stopping inside or in front of an organ at risk and putting restrictions on the minimum number and opening angle of incident beams for certain treatment sites. For the future, most centres (16) asked for more retrospective or prospective outcome studies investigating the potential effect of the effect of a variable RBE. To perform such studies, 18 centres asked for LET and RBE calculation and visualisation tools developed by treatment planning system vendors.

CONCLUSION

All European proton centres are aware of RBE variability but comply with current guidelines of prescribing a constant RBE. However, they actively mitigate uncertainty and risk of side effects resulting from increased RBE by applying measures and restrictions during treatment planning. To change RBE-related clinical guidelines in the future more clinical data on RBE are explicitly demanded.

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.

Mechanisms and Review of Clinical Evidence of Variations in Relative Biological Effectiveness in Proton Therapy

Proton therapy is increasingly being used as a radiation therapy modality. There is uncertainty about the biological effectiveness of protons relative to photon therapies as it depends on several physical and biological parameters. Radiation oncology currently applies a constant and generic value for the relative biological effectiveness (RBE) of 1.1, which was chosen conservatively to ensure tumor coverage. The use of a constant value has been challenged particularly when considering normal tissue constraints. Potential variations in RBE have been assessed in several published reviews but have mostly focused on data from clonogenic cell survival experiments with unclear relevance for clinical proton therapy. The goal of this review is to put in vitro findings in relation to clinical observations. Relevant in vivo pathways determining RBE for tumors and normal tissues are outlined, including not only damage to tumor cells and parenchyma but also vascular damage and immune response. Furthermore, the current clinical evidence of varying RBE is reviewed. The assessment can serve as guidance for treatment planning, personalized dose prescriptions, and outcome analysis.

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.

Applications of nanodosimetry in particle therapy planning and beyond

This topical review summarizes underlying concepts of nanodosimetry. It describes the development and current status of nanodosimetric detector technology. It also gives an overview of Monte Carlo track structure simulations that can provide nanodosimetric parameters for treatment planning of proton and ion therapy. Classical and modern radiobiological assays that can be used to demonstrate the relationship between the frequency and complexity of DNA lesion clusters and nanodosimetric parameters are reviewed. At the end of the review, existing approaches of treatment planning based on relative biological effectiveness (RBE) models or dose-averaged linear energy transfer are contrasted with an RBE-independent approach based on nandosimetric parameters. Beyond treatment planning, nanodosimetry is also expected to have applications and give new insights into radiation protection dosimetry.

The Organ Sparing Potential of Different Biological Optimization Strategies in Proton Therapy

Purpose

Variable relative biological effectiveness (RBE) models allow for differences in linear energy transfer (LET), physical dose, and tissue type to be accounted for when quantifying and optimizing the biological damage of protons. These models are complex and fraught with uncertainties, and therefore, simpler RBE optimization strategies have also been suggested. Our aim was to compare several biological optimization strategies for proton therapy by evaluating their performance in different clinical cases.

Methods and Materials

Two different optimization strategies were compared: full variable RBE optimization and differential RBE optimization, which involve applying fixed RBE for the planning target volume (PTV) and variable RBE in organs at risk (OARs). The optimization strategies were coupled to 2 variable RBE models and 1 LET-weighted dose model, with performance demonstrated on 3 different clinical cases: brain, head and neck, and prostate tumors.

Results

In cases with low (α/β)x in the tumor, the full RBE optimization strategies had a large effect, with up to 10% reduction in RBE-weighted dose to the PTV and OARs compared with the reference plan, whereas smaller variations (<5%) were obtained with differential optimization. For tumors with high (α/β)x, the differential RBE optimization strategy showed a greater reduction in RBE-weighted dose to the OARs compared with the reference plan and the full RBE optimization strategy.

Conclusions

Differences between the optimization strategies varied across the studied cases, influenced by both biological and physical parameters. Whereas full RBE optimization showed greater OAR sparing, awareness of underdosage to the target must be carefully considered.

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

BACKGROUND AND PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

Correlation of LET With MRI Changes in Brain and Potential Implications for Normal Tissue Complication Probability for Patients With Meningioma Treated With Pencil Beam Scanning Proton Therapy

PURPOSE

This study aimed to investigate the correlation between imaging changes in brain normal tissue and the spatial distribution of linear energy transfer (LET) for a cohort of patients with meningioma treated with scanned proton beams. Then, assuming imaging changes are induced by cell lethality, we studied the correlation between normal tissue complication probability and LET.

METHODS AND MATERIALS

Magnetic resonance imaging T2/fluid attenuated inversion recovery acquired at different intervals after proton radiation were coregistered with the planning computed tomography (CT) images from 26 patients with meningioma with abnormalities after proton radiation therapy. For this purpose, the T2/fluid attenuated inversion recovery areas not on the original magnetic resonance images were contoured, and the LET values for each voxel in the patient geometry were calculated to investigate the correlation between the position of imaging changes and the LET at those positions. To separate the effect of the dose as the inductor of these changes, we compared the LET in these areas with a sample of voxels matching the dose distributions across the image change areas. Patients with a higher LET in image change areas were grouped to verify whether they shared common characteristics.

RESULTS

Eleven of the patients showed higher dose-averaged LET (LETd) in imaging change regions than in the group of voxels with the same dose. This group of patients had significantly shallower targets for their treatment than the other 15 and used fewer beams and angles.

CONCLUSIONS

This study points toward the possibility that areas with imaging change are more likely to occur in regions with high dose or in areas with lower dose but increased LETd. The effect of LETd on imaging changes seems to be more relevant when treating superficial lesions with few nonopposed beams. However, most patients did not show a spatial correlation between their image changes and the LETd values, limiting the cases for the possible role of high LET as a toxicity inductor.

cRGDyK-modified procaine liposome inhibits the proliferation and motility of glioma cells via the ERK/p38MAPK pathway

Glioma is a common type of primary tumor in the central nervous system. Glioma has been increasing in incidence yearly and is a serious threat to human life and health. The aim of the present study was to prepare liposomes for enhanced penetration of the blood-brain barrier and targeting of glioma. A procaine-loaded liposome modified with the cyclic pentapeptide cRGDyK (Pro/cRGDyK-L) was designed and developed. The particle size, ζ potential, encapsulation efficiency, release profile, stability and hemolysis of Pro/cRGDyK-L were characterized in vitro. The targeting and antitumor effects of Pro/cRGDyK-L were also investigated in vitro and in vivo. The results suggested that the cRGDyK peptide significantly facilitated the ability of liposomes to transfer procaine across the BBB and improved the cellular uptake of procaine by C6 glioma cells. The results further demonstrated that Pro/cRGDyK-L strongly suppressed cell motility, stimulated apoptosis and induced cell cycle arrest. The findings further confirmed that Pro/cRGDyK-L exhibited superior antitumor effects by targeting the ERK/p38MAPK pathway and thereby suppressed tumor growth in mice. In conclusion, the present study indicated the potential of Pro/cRGDyK-L as a means to provide improved therapeutic effects on glioma through the ERK/p38MAPK pathway.

Brain-Specific Relative Biological Effectiveness of Protons Based on Long-term Outcome of Patients With Nasopharyngeal Carcinoma

PURPOSE

Uncertainties in relative biological effectiveness (RBE) constitute a major pitfall of the use of protons in clinics. An RBE value of 1.1, which is based on cell culture and animal models, is currently used in clinical proton planning. The purpose of this study was to determine RBE for temporal lobe radiographic changes using long-term follow-up data from patients with nasopharyngeal carcinoma.

METHODS AND MATERIALS

Five hundred sixty-six patients with newly diagnosed nasopharyngeal carcinoma received double-scattering proton therapy or intensity modulated radiation therapy at our institutions. The 2 treatment cohorts were well matched. Proton dose distributions were simulated using Monte Carlo and compared with those obtained from the proton clinical treatment planning system. Late treatment effect was defined as development of enhancement of temporal lobe on T1-weighted magnetic resonance imaging, with or without accompanying clinical symptoms. The tolerance dose was calculated with receiving operator characteristic analysis and the Youden index. Tolerance curves, expressed as a cumulative dose-volume histogram, were generated using the cutoff points.

RESULTS

With a median follow-up period >5 years for both cohorts, 10% of proton patients and 4% of patients undergoing intensity modulated radiation therapy developed temporal lobe enhancement in unilateral temporal lobe. There was no significant difference in dose distributions between the Monte Carlo method and treatment planning system. The tolerance dose-volume levels were V10 (26.1%), V20 (21.9%), V30 (14.0%), V40 (7.7%), V50 (4.8%), and V60 (3.3%) for proton therapy (P < .03). Comparison of the two tolerance curves revealed that tolerance doses of proton treatments were lower than that of photon treatments at all dose levels. The dose tolerance at D1% was 58.56 Gy for protons and 69.07 Gy for photons. The RBE for temporal lobe enhancement from proton treatments were calculated to be 1.18.

CONCLUSIONS

Using long-term clinical outcome of patients with nasopharyngeal carcinoma, our data suggest that the RBE for temporal lobe enhancement is 1.18 at D1%. A prospective study in a large cohort would be necessary to confirm these findings.

A framework for voxel-based assessment of biological effect after proton radiotherapy in pediatric brain cancer patients using multi-modal imaging

INTRODUCTION

The exact dependence of biological effect on dose and linear energy transfer (LET) in human tissue when delivering proton therapy is unknown. In this study, we propose a framework for measuring this dependency using multi-modal image-based assays with deformable registrations within imaging sessions and across time.

MATERIALS AND METHODS

3T MRI scans were prospectively collected from 6 pediatric brain cancer patients before they underwent proton therapy treatment, and every 3 months for a year after treatment. Scans included T1-weighted with contrast enhancement (T1), T2-FLAIR (T2) and fractional anisotropy (FA) images. In addition, the planning CT, dose distributions and Monte Carlo-calculated LET distributions were collected. A multi-modal deformable image registration framework was used to create a dataset of dose, LET and imaging intensities at baseline and follow-up on a voxel-by-voxel basis. We modelled the biological effect of dose and LET from proton therapy using imaging changes over time as a surrogate for biological effect. We investigated various models to show the feasibility of the framework to model imaging changes. To account for interpatient and intrapatient variations, we used a nested generalized linear mixed regression model. The models were applied to predict imaging changes over time as a function of dose and LET for each modality.

RESULTS

Using the nested models to predict imaging changes, we saw a decrease in the FA signal as a function of dose; however, the signal increased with increasing LET. Similarly, we saw an increase in T2 signal as a function of dose, but a decrease in signal with LET. We saw no changes in T1 voxel values as a function of either dose or LET.

CONCLUSIONS

The imaging changes could successfully model biological effect as a function of dose and LET using our proposed framework. Due to the low number of patients, the imaging changes observed for FA and T2 scans were not marked enough to draw any firm conclusions.

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.

Update of the EPTN atlas for CT- and MR-based contouring in Neuro-Oncology

BACKGROUND AND PURPOSE

To update the digital online atlas for organs at risk (OARs) delineation in neuro-oncology based on high-quality computed tomography (CT) and magnetic resonance (MR) imaging with new OARs.

MATERIALS AND METHODS

In this planned update of the neurological contouring atlas published in 2018, ten new clinically relevant OARs were included, after thorough discussion between experienced neuro-radiation oncologists (RTOs) representing 30 European radiotherapy-oncology institutes. Inclusion was based on daily practice and research requirements. Consensus was reached for the delineation after critical review. Contouring was performed on registered CT with intravenous (IV) contrast (soft tissue & bone window setting) and 3 Tesla (T) MRI (T1 with gadolinium & T2 FLAIR) images of one patient (1 mm slices). For illustration purposes, delineation on a 7 T MRI without IV contrast from a healthy volunteer was added. OARs were delineated by three experienced RTOs and a neuroradiologist based on the relevant literature.

RESULTS

The presented update of the neurological contouring atlas was reviewed and approved by 28 experts in the field. The atlas is available online and includes in total 25 OARs relevant to neuro-oncology, contoured on CT and MRI T1 and FLAIR (3 T & 7 T). Three-dimensional (3D) rendered films are also available online.

CONCLUSION

In order to further decrease inter- and intra-observer OAR delineation variability in the field of neuro-oncology, we propose the use of this contouring atlas in photon and particle therapy, in clinical practice and in the research setting. The updated atlas is freely available on www.cancerdata.org.

Clinical Implementation of Proton Therapy Using Pencil-Beam Scanning Delivery Combined With Static Apertures

Proton therapy makes use of the favorable depth-dose distribution with its characteristic Bragg peak to spare normal tissue distal of the target volume. A steep dose gradient would be desired in lateral dimensions, too. The widespread spot scanning delivery technique is based, however, on pencil-beams with in-air spot full-widths-at-half-maximum of typically 1 cm or more. This hampers the sparing of organs-at-risk if small-scale structures adjacent to the target volume are concerned. The trimming of spot scanning fields with collimating apertures constitutes a simple measure to increase the transversal dose gradient. The current study describes the clinical implementation of brass apertures in conjunction with the pencil-beam scanning delivery mode at a horizontal, clinical treatment head based on commercial hardware and software components. Furthermore, clinical cases, which comprised craniopharyngiomas, re-irradiations and ocular tumors, were evaluated. The dosimetric benefits of 31 treatment plans using apertures were compared to the corresponding plans without aperture. Furthermore, an overview of the radiation protection aspects is given. Regarding the results, robust optimization considering range and setup uncertainties was combined with apertures. The treatment plan optimizations followed a single-field uniform dose or a restricted multi-field optimization approach. Robustness evaluation was expanded to account for possible deviations of the center of the pencil-beam delivery and the mechanical center of the aperture holder. Supplementary apertures improved the conformity index on average by 15.3%. The volume of the dose gradient surrounding the PTV (evaluated between 80 and 20% dose levels) was decreased on average by 17.6%. The mean dose of the hippocampi could be reduced on average by 2.9 GyRBE. In particular cases the apertures facilitated a sparing of an organ-at-risk, e.g. the eye lens or the brainstem. For six craniopharyngioma cases the inclusion of apertures led to a reduction of the mean dose of 1.5 GyRBE (13%) for the brain and 3.1 GyRBE (16%) for the hippocampi.

Radiation-induced brain injury in patients with meningioma treated with proton or photon therapy

INTRODUCTION

Radiation therapy is often used to treat meningioma with adverse features or when unresectable. Proton therapy has advantages over photon therapy in reducing integral dose to the brain. This study compared the incidence of radiological and clinical adverse events after photon versus proton therapy in the treatment of meningioma.

METHODS

A retrospective review was conducted on patients with meningioma treated with proton or photon therapy at two high-volume tertiary cancer centers. Patients with a history of prior radiation therapy (RT) or less than 3 months of follow-up were excluded. Post-RT imaging changes were categorized into abnormal T2 signal intensities (T2 changes) or abnormal T1 post-contrast and T2 signal intensities (T1c+T2 changes) on magnetic resonance imaging (MRI). Clinical outcomes of adverse events and survival were compared between the proton and photon therapies.

RESULTS

Among the total of 77 patients, 38 patients received proton therapy and 39 patients received photon therapy. The median age at diagnosis was 55 years and median follow-up was 2.2 years. No significant differences in symptomatic adverse events were observed between the two groups: grade ≥ 2 adverse events were seen in 4 (10.5%) patients in the proton group and 3 (7.7%) patients in the photon group (p = 0.67). The 2-year cumulative incidences of T2 changes were 38.3% after proton therapy and 47.7% after photon therapy (p = 0.53) and the 2-year cumulative incidences of T1c+T2 changes were 26.8% after proton therapy and 5.3% after photon therapy (p = 0.02). One patient experienced grade ≥ 4 adverse event in each group (p = 0.99). Estimated 2-year progression-free survival was 79.5% (proton therapy 76.0% vs. photon therapy 81.3%, p = 0.66) and 2-year overall survival was 89.7% (proton therapy 86.6% vs. photon therapy 89.3%, p = 0.65).

CONCLUSIONS

Following RT, high rates of T2 changes were seen in meningioma patients regardless of treatment modality. Proton therapy was associated with significantly higher rates of T1c+T2 changes compared with photon therapy, but severe adverse events were uncommon in both groups and survival outcomes were comparable between the two groups. Future studies will aim at correlating the MRI changes with models that can be incorporated into RT planning to avoid toxicity.

Identification of patient benefit from proton beam therapy in brain tumour patients based on dosimetric and NTCP analyses

BACKGROUND

The limited availability of proton beam therapy (PBT) requires individual treatment selection strategies, such as the model-based approach. In this study, we assessed the dosimetric benefit of PBT compared to photon therapy (XRT), analysed the corresponding changes in normal tissue complication probability (NTCP) on a variety of available models, and illustrated model-based patient selection in an in-silico study for patients with brain tumours.

METHODS

For 92 patients treated at two PBT centres, volumetric modulated arc therapy treatment plans were retrospectively created for comparison with the clinically applied PBT plans. Several dosimetric parameters for the brain excluding tumour and margins, cerebellum, brain stem, frontal and temporal lobes, hippocampi, cochleae, chiasm, optic nerves, lacrimal glands, lenses, pituitary gland, and skin were compared between both modalities using Wilcoxon signed-rank tests. NTCP differences (ΔNTCP) were calculated for 11 models predicting brain necrosis, delayed recall, temporal lobe injury, hearing loss, tinnitus, blindness, ocular toxicity, cataract, endocrine dysfunction, alopecia, and erythema. A patient was assumed to be selected for PBT if ΔNTCP exceeded a threshold of 10 percentage points for at least one of the side-effects.

RESULTS

PBT substantially reduced the dose in almost all investigated OARs, especially in the low and intermediate dose ranges and for contralateral organs. In general, NTCP predictions were significantly lower for PBT compared to XRT, in particular in ipsilateral organs. Considering ΔNTCP of all models, 80 patients (87.0%) would have been selected for PBT in this in-silico study, mainly due to predictions of a model on delayed recall (51 patients).

CONCLUSION

In this study, substantial dose reductions for PBT were observed, mainly in contralateral organs. However, due to the sigmoidal dose response, NTCP was particularly reduced in ipsilateral organs. This underlines that physical dose-volume parameters alone may not be sufficient to describe the clinical relevance between different treatment techniques and highlights potential benefits of NTCP models. Further NTCP models for different modern treatment techniques are mandatory and existing models have to be externally validated in order to implement the model-based approach in clinical practice for cranial radiotherapy.

Distinct imaging patterns of pseudoprogression in glioma patients following proton versus photon radiation therapy

PURPOSE

Criteria by the Radiologic Assessment in Neuro-Oncology (RANO) group outline the diagnosis of pseudoprogression (Ps) after photon therapy for gliomas based on timing and location. We noted that patients receiving proton therapy manifested radiographic changes that appear different than Ps after photon therapy, which could be interpreted as tumor progression. In this study, we retrospectively reviewed MR imaging after proton or photon radiation for gliomas. We propose criteria to characterize proton pseudoprogression (ProPs) as distinct from Ps seen after photons.

METHODS

Post-treatment MR imaging, clinical and pathological data of low grade glioma patients were reviewed. Overall, 57 patients receiving protons were reviewed for the presence of ProPs, and 43 patients receiving photons were reviewed for any equivalent imaging changes. Data collected included the location and timing of the new enhancement, tumor grade, molecular subtype, chemotherapy received, and clinical symptoms.

RESULTS

Fourteen patients (24.6%) had new enhancement following radiation therapy that was unique to treatment with protons. The mean time to development of the ProPs was 15.4 months (7-27 months). We established the following criteria to characterize ProPs: located at the distal end of the proton beam; resolves without tumor-directed therapy; and subjectively multifocal, patchy, and small (< 1 cm). In the group receiving photons, none had changes that met our criteria for ProPs.

CONCLUSION

Patients who receive protons have unique imaging changes after radiation therapy. ProPs could be mistaken for tumor progression, but typically resolves on follow up. Further studies are needed to understand the radiobiology and pathophysiology underlying these imaging changes.

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.

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.

Brain Necrosis in Adult Patients After Proton Therapy: Is There Evidence for Dependency on Linear Energy Transfer?

PURPOSE

To investigate if radiographic imaging changes defined as necrosis correlate with regions in the brain with elevated linear energy transfer (LET) for proton radiation therapy treatments with partial brain involvement in central nervous system and patients with head and neck cancer.

METHODS AND MATERIALS

Fifty patients with head and neck, skull base, or intracranial tumors who underwent proton therapy between 2004 to 2016 with a minimum prescription dose of 59.4 Gy (relative biological effectiveness) and with magnetic resonance imaging changes indicative of brain necrosis after radiation therapy were retrospectively reviewed. Each treatment plan was recalculated using Monte Carlo simulations to provide accurate dose distributions as well as 3-dimensional distributions of LET. To assess the effect of LET on radiographic imaging changes several voxel-based analyses were performed.

RESULTS

In this patient cohort, LET adjusted for dose was not found to be associated with risk of brain necrosis.

CONCLUSIONS

A voxel-based analysis of brain necrosis as an endpoint is difficult owing to uncertainties in the origin of necrosis, timing of imaging, variability in patient specific radiosensitivity, and the simultaneous effect of dose and LET. Even though it is expected that the LET and thus relative biological effectiveness increases at the end of range, effects in patients might be small compared with interpatient variability of radiosensitivity and might be obscured by other confounding factors.

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.

Proton therapy for selected low grade glioma patients in the Netherlands

Proton therapy offers an attractive alternative to conventional photon-based radiotherapy in low grade glioma patients, delivering radiotherapy with equivalent efficacy to the tumour with less radiation exposure to the brain. In the Netherlands, patients with favourable prognosis based on tumour and patient characteristics can be offered proton therapy. Radiation-induced neurocognitive function decline is a major concern in these long surviving patients. Although level 1 evidence of superior clinical outcome with proton therapy is lacking, the Dutch National Health Care Institute concluded that there is scientific evidence to assume that proton therapy can have clinical benefit by reducing radiation-induced brain damage. Based on this decision, proton therapy is standard insured care for selected low grade glioma patients. Patients with other intracranial tumours can also qualify for proton therapy, based on the same criteria. In this paper, the evidence and considerations that led to this decision are summarised. Additionally, the eligibility criteria for proton therapy and the steps taken to obtain high-quality data on treatment outcome are discussed.

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.

Dose-averaged linear energy transfer per se does not correlate with late rectal complications in carbon-ion radiotherapy

BACKGROUND AND PURPOSE

Several studies have focused on increasing the linear energy transfer (LET) within tumours to achieve higher biological effects in carbon-ion radiotherapy (C-ion RT). However, it remains unclear whether LET affects late complications. We assessed whether physical dose and LET distribution can be specific factors for late rectal complications in C-ion RT.

MATERIALS AND METHODS

Overall, 134 patients with uterine carcinomas were registered and retrospectively analysed. Of 134 patients, 132 who were followed up for >6 months were enrolled. The correlations between the relative biological effectiveness (RBE)-weighted dose based on the Kanai model (the ostensible "clinical dose"), dose-averaged LET (LETd), or physical dose and rectal complications were evaluated. Rectal complications were graded according to the Radiation Therapy Oncology Group/European Organization for Research and Treatment of Cancer criteria.

RESULTS

Nine patients developed grade 3 or 4 late rectal complications. Linear regression analysis found that D_2cc in clinical dose was the sole risk factor for ≥grade 3 late rectal complications (p = 0.012). The receiver operating characteristic analysis found that D_2cc of 60.2 Gy (RBE) was a suitable cut-off value for predicting ≥grade 3 late rectal complications. Among 35 patients whose rectal D_2cc was ≥60.2 Gy (RBE), no correlations were found between severe rectal toxicities and LETd alone or physical dose per se.

CONCLUSION

We demonstrated that severe rectal toxicities were related to the rectal D_2cc of the clinical dose in C-ion RT. However, no correlations were found between severe rectal toxicities and LETd alone or physical dose per se.

Exploiting the full potential of proton therapy: An update on the specifics and innovations towards spatial or temporal optimisation of dose delivery

Prescription and delivery of protons are somewhat different compared to photons and may influence outcomes (tumour control and toxicity). These differences should be taken into account to fully exploit the clinical potential of proton therapy. Innovations in proton therapy treatment are also required to widen the therapeutic window and determine appropriate populations of patients that would benefit from new treatments. Therefore, strategies are now being developed to reduce side effects to critical normal tissues using alternative treatment configurations and new spatial or temporal-driven optimisation approaches. Indeed, spatiotemporal optimisation (based on flash, proton minibeam radiation therapy or hypofractionated delivery methods) has been gaining some attention in proton therapy as a mean of improving (biological and physical) dose distribution. In this short review, the main differences in planning and delivery between protons and photons, as well as some of the latest developments and methodological issues (in silico modelling) related to providing scientific evidence for these new techniques will be discussed.