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

Development and clinical application of a probabilistic robustness evaluation tool for pencil beam scanning proton therapy treatments

PURPOSE

to implement a probabilistic-Robustness-Evaluation (pRE) tool for proton therapy treatments and to correlate these results with the worst-case approach (wRE) implemented in commercial TPS for clinical applications.

MATERIALS AND METHODS

12 skull base patients were planned with a robust multiple field optimization (MFO) approach. 10 years of machine QA were analysed to derive the uncertainties of our treatment system (beam delivery and patient positioning system). For a large cohort of patients, post-treatment imaging was acquired to determine the intra-fraction uncertainty. The pRE, considered explicitly all these uncertainties, the fractionation and range uncertainty. For each plan a wRE with different combinations of range and setup uncertainties was simulated. wRE results were then compared, in terms of target coverage and OAR dose limits, with pRE results.

RESULTS

43,400 dose distributions were analysed. pRE simulations lasted 18.6 h (±11.5 h). The results showed that the combination of wRE uncertainty parameters that surrogated the best pRE results with a confidence level of 95 % were (1.0 mm/3.5 %). The median OAR's dose indexes difference (D1/D1cc) between pRE and wRE was 1.90 (±1.49) GyRBE, while for target D98 and D95 it was -0.66(±0.95) and -0.67 (±0.52) GyRBE, respectively.

CONCLUSION

A tool able to explicitly simulate the source of treatment uncertainties and the effect of the fractionation was implemented to have a more realistic evaluation of plan robustness. This tool was used to find the best wRE parameters that surrogate the pRE results while maintaining clinically acceptable timing. These results are now used in our clinical workflow.

Linear approximation of variable relative biological effectiveness models for proton therapy

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

Robustness evaluation of pencil beam scanning proton therapy treatment planning: A systematic review

Compared to conventional radiotherapy using X-rays, proton therapy, in principle, allows better conformity of the dose distribution to target volumes, at the cost of greater sensitivity to physical, anatomical, and positioning uncertainties. Robust planning, both in terms of plan optimization and evaluation, has gained high visibility in publications on the subject and is part of clinical practice in many centers. However, there is currently no consensus on the methods and parameters to be used for robust optimization or robustness evaluation. We propose to overcome this deficiency by following the modified Delphi consensus method. This method first requires a systematic review of the literature. We performed this review using the PubMed and Web Of Science databases, via two different experts. Potential conflicts were resolved by a third expert. We then explored the different methods before focusing on clinical studies that evaluate robustness on a significant number of patients. Many robustness assessment methods are proposed in the literature. Some are more successful than others and their implementation varies between centers. Moreover, they are not all statistically or mathematically equivalent. The most sophisticated and rigorous methods have seen more limited application due to the difficulty of their implementation and their lack of widespread availability.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Multicentric and collaborative study of Spanish neurosurgical management of pediatric craniopharyngiomas: S-PedCPG.co

PURPOSE

To present a descriptive analysis of pediatric craniopharyngiomas (PedCPG) treated in various Spanish hospitals, defining factors related to recurrence and performing a critical analysis of the results.

METHODS

We undertook a multicenter retrospective review of PedCPG treated between 2000 and 2017. Data collected included epidemiological variables, clinical and radiological characteristics, goal of first surgery, rate of recurrence and its approach, adjuvant treatment, complications and permanent morbidity. Associations were studied between progression and number of progressions and independent variables.

RESULTS

The study involved 69 children from 8 Spanish hospitals. Most of the tumors invaded several intracranial compartments at diagnosis, with the hypothalamus involved in 41.3% of cases. The first treatment strategy was usually gross total resection (GTR) (71%), with some patients treated with radiotherapy or intracystic chemotherapy. The progression rate after first surgery was 53% in a mean follow-up of 88.2 months (range 7-357). In the GTR group 38.8% of tumors recurred, 40% in the group of subtotal resection or biopsy and 93.3% in the cyst fenestration±Ommaya reservoir group. Mortality was 7.2%. Follow-up period, size of the tumor and goal of first surgery were significantly related with progression.

CONCLUSIONS

Our results in terms of disease control, hormonal or visual impairment and mortality were acceptable, but there are several areas for improvement. Our short-term goals should be to create a national register of PedCPG, reach a consensus about a treatment algorithm, and improve diagnosis of hypothalamic dysfunction to avoid preventable morbidity.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

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.

Microdosimetry in low energy proton beam at therapeutic-equivalent fluence rate with silicon 3D-cylindrical microdetectors

In this work we show the first microdosimetry measurements on a low energy proton beam with therapeutic-equivalent fluence rates by using the second generation of 3D-cylindrical microdetectors. The sensors belong to an improved version of a novel silicon-based 3D-microdetector design with electrodes etched inside silicon, which were manufactured at the National Microelectronics Centre (IMB-CNM, CSIC) in Spain. A new microtechnology has been employed using quasi-toroid electrodes of 25μm diameter and a depth of 20μm within the silicon bulk, resulting in a well-defined cylindrical radiation sensitive volume. These detectors were tested at the 18 MeV proton beamline of the cyclotron at the National Accelerator Centre (CNA, Spain). They were assembled into an in-house low-noise readout electronics to assess their performance at a therapeutic-equivalent fluence rate. Microdosimetry spectra of lineal energy were recorded at several proton energies starting from 18 MeV by adding 50μm thick tungsten foils gradually at the exit-window of the cyclotron external beamline, which corresponds to different depths along the Bragg curve. The experimentalyF¯values in silicon cover from (5.7 ± 0.9) to (8.5 ± 0.4) keV μm^-1in the entrance to (27.4 ± 2.3) keV μm^-1in the distal edge. Pulse height energy spectra were crosschecked with Monte Carlo simulations and an excellent agreement was obtained. This work demonstrates the capability of the second generation 3D-microdetectors to assess accurate microdosimetric distributions at fluence rates as high as those used in clinical centers in proton therapy.

Neurocognitive impairment following proton therapy for paediatric brain tumour: a systematic review of post-therapy assessments

BACKGROUND

Proton therapy (PT), frequently utilised to treat paediatric brain tumour (PBT) patients, eliminates exit dose and minimises dose to healthy tissues that theoretically can mitigate treatment-related effects including cognitive deficits. As clinical outcome data are emerging, we aimed to systematically review current evidence of cognitive changes following PT of PBT.

MATERIALS AND METHODS

We searched PubMed and Scopus electronic databases to identify eligible reports on cognitive changes following PT of PBT according to PRISMA guidelines. Reports were extracted for information on demographics and cognitive outcomes. Then, they were systematically reviewed based on three themes: (1) comparison with photon therapy, (2) comparison with baseline cognitive measures, to population normative mean or radiotherapy-naïve PBT patients and (3) effects of dose distribution to cognition.

RESULTS

Thirteen reports (median size (range): 70 (12-144)) were included. Four reports compared the cognitive outcome between PBT patients treated with proton to photon therapy and nine compared with baseline/normative mean/radiotherapy naïve from which two reported the effects of dose distribution. Reports found significantly poorer cognitive outcome among patients treated with photon therapy compared with proton therapy especially in general cognition and working memory. Craniospinal irradiation (CSI) was consistently associated with poorer cognitive outcome while focal therapy was associated with minor cognitive change/difference. In limited reports available, higher doses to the hippocampus and temporal lobes were implicated to larger cognitive change.

CONCLUSION

Available evidence suggests that PT causes less cognitive deficits compared with photon therapy. Children who underwent focal therapy with proton were consistently shown to have low risk of cognitive deficit suggesting the need for future studies to separate them from CSI. Evidence on the effect of dose distribution to cognition in PT is yet to mature.

RBE for proton radiation therapy - a Nordic view in the international perspective

BACKGROUND

This paper presents an insight into the critical discussions and the current strategies of the Nordic countries for handling the variable proton relative biological effectiveness (RBE) as presented at The Nordic Collaborative Workshop for Particle Therapy that took place at the Skandion Clinic on 14th and 15th of November 2019.

MATERIAL AND METHODS

In the current clinical practice at the two proton centres in operation at the date, Skandion Clinic, and the Danish Centre for Particle Therapy, a constant proton RBE of 1.1 is applied. The potentially increased effectiveness at the end of the particle range is however considered at the stage of treatment planning at both places based on empirical observations and knowledge. More elaborated strategies to evaluate the plans and mitigate the problem are intensely investigated internationally as well at the two centres. They involve the calculation of the dose-averaged linear energy transfer (LET_d) values and the assessment of their distributions corroborated with the distribution of the dose and the location of the critical clinical structures.

RESULTS

Methods and tools for LET_d calculations are under different stages of development as well as models to account for the variation of the RBE with LET_d, dose per fraction, and type of tissue. The way they are currently used for evaluation and optimisation of the plans and their robustness are summarised. A critical but not exhaustive discussion of their potential future implementation in the clinical practice is also presented.

CONCLUSIONS

The need for collaboration between the clinical proton centres in establishing common platforms and perspectives for treatment planning evaluation and optimisation is highlighted as well as the need of close interaction with the research academic groups that could offer a complementary perspective and actively help developing methods and tools for clinical implementation of the more complex metrics for considering the variable effectiveness of the proton beams.

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.