Relating the proton relative biological effectiveness to tumor control and normal tissue complication probabilities assuming interpatient variability in α/β

A joint physics and radiobiology DREAM team vision - Towards better response prediction models to advance radiotherapy

Radiotherapy developed empirically through experience balancing tumour control and normal tissue toxicities. Early simple mathematical models formalized this practical knowledge and enabled effective cancer treatment to date. Remarkable advances in technology, computing, and experimental biology now create opportunities to incorporate this knowledge into enhanced computational models. The ESTRO DREAM (Dose Response, Experiment, Analysis, Modelling) workshop brought together experts across disciplines to pursue the vision of personalized radiotherapy for optimal outcomes through advanced modelling. The ultimate vision is leveraging quantitative models dynamically during therapy to ultimately achieve truly adaptive and biologically guided radiotherapy at the population as well as individual patient-based levels. This requires the generation of models that inform response-based adaptations, individually optimized delivery and enable biological monitoring to provide decision support to clinicians. The goal is expanding to models that can drive the realization of personalized therapy for optimal outcomes. This position paper provides their propositions that describe how innovations in biology, physics, mathematics, and data science including AI could inform models and improve predictions. It consolidates the DREAM team's consensus on scientific priorities and organizational requirements. Scientifically, it stresses the need for rigorous, multifaceted model development, comprehensive validation and clinical applicability and significance. Organizationally, it reinforces the prerequisites of interdisciplinary research and collaboration between physicians, medical physicists, radiobiologists, and computational scientists throughout model development. Solely by a shared understanding of clinical needs, biological mechanisms, and computational methods, more informed models can be created. Future research environment and support must facilitate this integrative method of operation across multiple disciplines.

Brainstem toxicity after proton or photon therapy in children and young adults with localized intracranial ependymoma: A French retrospective study

BACKGROUND AND PURPOSE

Ependymoma is the third most frequent childhood braintumor. Standard treatment is surgery followed by radiation therapy including proton therapy (PBT). Retrospective studies have reported higher rates of brainstem injury after PBT than after photon therapy (XRT). We report a national multicenter study of the incidence of brainstem injury after XRT versus PBT, and their correlations with dosimetric data.

MATERIAL AND METHODS

We included all patients aged < 25 years who were treated with PBT or XRT for intracranial ependymoma at five French pediatric oncology reference centers between 2007 and 2020. We reviewed pre-irradiation MRI, follow-up MRIs over the 12 months post-treatment and clinical data.

RESULTS

Of the 83 patients, 42 were treated with PBT, 37 with XRT, and 4 with both (median dose: 59.4 Gy, range: 53‑60). No new or progressive symptomatic brainstem injury was found. Four patients presented asymptomatic radiographic changes (punctiform brainstem enhancement and FLAIR hypersignal), with median onset at 3.5 months (range: 3.0‑9.4) after radiation therapy, and median offset at 7.6 months (range: 3.7‑7.9). Two had been treated with PBT, one with XRT, and one with mixed XRT-PBT. Prescribed doses were 59.4, 55.8, 59.4 and 54 Gy.

CONCLUSION

Asymptomatic radiographic changes occurred in 4.8% of patients with ependymoma in a large national series. There was no correlation with dose or technique. No symptomatic brainstem injury was identified.

Combinatorial approach of immuno-proton therapy in cancer: Rationale and potential impact

Cancer management is an expansive, growing, and evolving field. In the last decade or so, immunotherapy (IT) and particle beam therapy have made a tremendous impact in this domain. IT has already established itself as the fourth pillar of oncology. Recent emphasis has been centred around combination therapy, postulating additive or multiplicative effects of combining IT with one or more of the three conventional "pillars," that is, surgery, chemotherapy, and radiotherapy. Radio-IT is being increasingly explored and has shown promising outcomes in both preclinical and clinical settings. Particle beam therapy such as protons, when used as the radiotherapeutic modality in conjunction with IT, can potentially limit toxicities and improve this synergism further. Modern proton therapy has demonstrated a reduction in integral dose of radiation and radiation-induced lymphopenia in various sites. Protons, by virtue of their inherent clinically desirable physical and biological characteristics, namely, high linear energy transfer, relative biological effectiveness of range 1.1-1.6, and proven anti-metastatic and immunogenic potential in preclinical studies, might have a superior immunogenic profile than photons. Proton-IT combination is being studied currently by various groups in lung , head neck and brain tumors, and should be evaluated further in other subsites to replicate preclinical outcomes in a clinical setting. In this review, we summarize the currently available evidence for combinatorial approaches and feasibility of proton and IT combination, and thereafter highlight the emerging challenges for practical application of the same in clinics, while also proposing plausible solutions.

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.

A systematic approach for calibrating a Monte Carlo code to a treatment planning system for obtaining dose, LET, variable proton RBE and out-of-field dose

Objective. While integration of variable relative biological effectiveness (RBE) has not reached full clinical implementation, the importance of having the ability to recalculate proton treatment plans in a flexible, dedicated Monte Carlo (MC) code cannot be understated . Here we provide a step-wise method for calibrating dose from a MC code to a treatment planning system (TPS), to obtain required parameters for calculating linear energy transfer (LET), variable RBE and in general enabling clinical realistic research studies beyond the capabilities of a TPS.Approach. Initially, Pristine Bragg peaks (PBP) were calculated in both the Eclipse TPS and the FLUKA MC code. A rearranged Bortfeld energy-range relation was applied to the initial energy of the beam to fine-tune the range of the MC code at 80% dose level distal to the PBP. The energy spread was adapted by dividing the TPS range by the MC range for dose level 80%-20% distal to the PBP. Density and relative proton stopping power were adjusted by comparing the TPS and MC for different Hounsfield units. To find the relationship of dose per primary particle from the MC to dose per monitor unit in the TPS, integration was applied to the area of the Bragg curve. The calibration was validated for spread-out Bragg peaks (SOBP) in water and patient treatment plans. Following the validation, variable RBE were calculated using established models.Main results.The PBPs ranges were within ±0.3mm threshold, and a maximum of 5.5% difference for the SOBPs was observed. The patient validation showed excellent dose agreement between the TPS and MC, with the greatest differences for the lung tumor patient.Significance. Aprocedure for calibrating a MC code to a TPS was developed and validated. The procedure enables MC-based calculation of dose, LET, variable RBE, advanced (secondary) particle tracking and more from treatment plans.

Modelling the influence of radiosensitivity on development of second primary cancer in out-of-field organs following proton therapy for paediatric cranial cancer

OBJECTIVE

Radiobiological modelling the risks of second primary cancer (SPC) after proton therapy (PT) for childhood cranial cancer remains largely unknown. Organ-specific dose-response risk factors such as radiosensitivity require exploration. This study compared the influence of radiosensitivity data (slope of βEAR) on children's lifetime attributable risks (LAR) of SPC development in out-of-field organs following cranial scattering and scanning PT.

METHODS

Out-of-field radiosensitivity parameter estimates for organs (α/β and βEAR) were sourced from literature. Physical distances for 13 out-of-field organs were measured and input into Schneider's SPC model. Sensitivity analyses were performed as a function of radiosensitivity (α/β of 1-10 Gy) and initial slope (βEAR) from Japanese/UK data to estimate the influence on the risk of radiation-induced SPC following scattering and scanning PT.

RESULTS

Models showed similar LAR of SPC estimates for age and sex-matched paediatric phantoms, however, for breast there was a significant increase using Japanese βEAR data. For most organs, scattering PT demonstrated a larger risk of LAR for SPC which increased with α/β.

CONCLUSION

Breast tissue exhibited the highest susceptibility in calculated LAR risk, demonstrating the importance for accurate data input when estimating LAR of SPC.

ADVANCES IN KNOWLEDGE

The findings of this study demonstrated younger female patients undergoing cranial proton therapy have a higher risk of developing second primary cancer of the breast tissue. Long-term multicenter registries are important to improve predictive radiobiological modelling studies of side effects.

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.

Proton therapy and limited surgery for paediatric and adolescent patients with craniopharyngioma (RT2CR): a single-arm, phase 2 study

BACKGROUND

Compared with photon therapy, proton therapy reduces exposure of normal brain tissue in patients with craniopharyngioma, which might reduce cognitive deficits associated with radiotherapy. Because there are known physical differences between the two methods of radiotherapy, we aimed to estimate progression-free survival and overall survival distributions for paediatric and adolescent patients with craniopharyngioma treated with limited surgery and proton therapy, while monitoring for excessive CNS toxicity.

METHODS

In this single-arm, phase 2 study, patients with craniopharyngioma at St Jude Children's Research Hospital (Memphis TN, USA) and University of Florida Health Proton Therapy Institute (Jacksonville, FL, USA) were recruited. Patients were eligible if they were aged 0-21 years at the time of enrolment and had not been treated with previous radiotherapeutic or intracystic therapies. Eligible patients were treated using passively scattered proton beams, 54 Gy (relative biological effect), and a 0·5 cm clinical target volume margin. Surgical treatment was individualised before proton therapy and included no surgery, single procedures with catheter and Ommaya reservoir placement through a burr hole or craniotomy, endoscopic resection, trans-sphenoidal resection, craniotomy, or multiple procedure types. After completing treatment, patients were evaluated clinically and by neuroimaging for tumour progression and evidence of necrosis, vasculopathy, permanent neurological deficits, vision loss, and endocrinopathy. Neurocognitive tests were administered at baseline and once a year for 5 years. Outcomes were compared with a historical cohort treated with surgery and photon therapy. The coprimary endpoints were progression-free survival and overall survival. Progression was defined as an increase in tumour dimensions on successive imaging evaluations more than 2 years after treatment. Survival and safety were also assessed in all patients who received photon therapy and limited surgery. This study is registered with ClinicalTrials.gov, NCT01419067.

FINDINGS

Between Aug 22, 2011, and Jan 19, 2016, 94 patients were enrolled and treated with surgery and proton therapy, of whom 49 (52%) were female, 45 (48%) were male, 62 (66%) were White, 16 (17%) were Black, two (2%) were Asian, and 14 (15%) were other races, and median age was 9·39 years (IQR 6·39-13·38) at the time of radiotherapy. As of data cutoff (Feb 2, 2022), median follow-up was 7·52 years (IQR 6·28-8·53) for patients who did not have progression and 7·62 years (IQR 6·48-8·54) for the full cohort of 94 patients. 3-year progression-free survival was 96·8% (95% CI 90·4-99·0; p=0·89), with progression occurring in three of 94 patients. No deaths occurred at 3 years, such that overall survival was 100%. At 5 years, necrosis had occurred in two (2%) of 94 patients, severe vasculopathy in four (4%), and permanent neurological conditions in three (3%); decline in vision from normal to abnormal occurred in four (7%) of 54 patients with normal vision at baseline. The most common grade 3-4 adverse events were headache (six [6%] of 94 patients), seizure (five [5%]), and vascular disorders (six [6%]). No deaths occurred as of data cutoff.

INTERPRETATION

Proton therapy did not improve survival outcomes in paediatric and adolescent patients with craniopharyngioma compared with a historical cohort, and severe complication rates were similar. However, cognitive outcomes with proton therapy were improved over photon therapy. Children and adolescents treated for craniopharyngioma using limited surgery and post-operative proton therapy have a high rate of tumour control and low rate of severe complications. The outcomes achieved with this treatment represent a new benchmark to which other regimens can be compared.

FUNDING

American Lebanese Syrian Associated Charities, American Cancer Society, the US National Cancer Institute, and Research to Prevent Blindness.

Does variable RBE affect toxicity risks for mediastinal lymphoma patients? NTCP-based evaluation after proton therapy treatment

INTRODUCTION

Mediastinal lymphoma (ML) is a solid malignancy affecting young patients. Modern combined treatments allow obtaining good survival probability, together with a long life expectancy, and therefore with the need to minimize treatment-related toxicities. We quantified the expected toxicity risk for different organs and endpoints in ML patients treated with intensity-modulated proton therapy (IMPT) at our centre, accounting also for uncertainties related to variable RBE.

METHODS

Treatment plans for ten ML patients were recalculated with a TOPAS-based Monte Carlo code, thus retrieving information on LET and allowing the estimation of variable RBE. Published NTCP models were adopted to calculate the toxicity risk for hypothyroidism, heart valve defects, coronary heart disease and lung fibrosis. NTCP was calculated assuming both constant (i.e. 1.1) and variable RBE. The uncertainty associated with individual radiosensitivity was estimated by random sampling α/β values before RBE evaluation.

RESULTS

Variable RBE had a minor impact on hypothyroidism risk for 7 patients, while it led to significant increase for the remaining three (+24% risk maximum increase). Lung fibrosis was slightly affected by variable RBE, with a maximum increase of ≅ 1%. This was similar for heart valve dysfunction, with the exception of one patient showing an about 10% risk increase, which could be explained by means of large heart volume and D₁ increase.

DISCUSSION

The use of NTCP models allows for identifying those patients associated with a higher toxicity risk. For those patients, it might be worth including variable RBE in plan evaluation.

The Mayo Clinic Florida Microdosimetric Kinetic Model of Clonogenic Survival: Application to Various Repair-Competent Rodent and Human Cell Lines

The relative biological effectiveness (RBE) calculations used during the planning of ion therapy treatments are generally based on the microdosimetric kinetic model (MKM) and the local effect model (LEM). The Mayo Clinic Florida MKM (MCF MKM) was recently developed to overcome the limitations of previous MKMs in reproducing the biological data and to eliminate the need for ion-exposed in vitro data as input for the model calculations. Since we are considering to implement the MCF MKM in clinic, this article presents (a) an extensive benchmark of the MCF MKM predictions against corresponding in vitro clonogenic survival data for 4 rodent and 10 cell lines exposed to ions from ¹H to ²³⁸U, and (b) a systematic comparison with published results of the latest version of the LEM (LEM IV). Additionally, we introduce a novel approach to derive an approximate value of the MCF MKM model parameters by knowing only the animal species and the mean number of chromosomes. The overall good agreement between MCF MKM predictions and in vitro data suggests the MCF MKM can be reliably used for the RBE calculations. In most cases, a reasonable agreement was found between the MCF MKM and the LEM IV.

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.

On the calculation of the relative biological effectiveness of ion radiation therapy using a biological weighting function, the microdosimetric kinetic model (MKM) and subsequent corrections (non-Poisson MKM and modified MKM)

Objective. To investigate similarities and differences in the formalism, processing, and the results of relative biological effectiveness (RBE) calculations with a biological weighting function (BWF), the microdosimetric kinetic model (MKM) and subsequent modifications (non-Poisson MKM, modified MKM). This includes: (a) the extension of the V79-RBE10% BWF to model the RBE for other clonogenic survival levels; (b) a novel implementation of MKMs as weighting functions; (c) a benchmark against Chinese Hamster lung fibroblast (V79) in vitro data; (d) a study on the effect of pre- or post- processing the average biophysical quantities used for the RBE calculations; (e) a possible modification of the modified MKM parameters to improve the model accuracy at high linear energy transfer (LET). Methodology. Lineal energy spectra were simulated for two spherical targets (diameter = 0.464 or 1.0 μm) using PHITS for 1H, 4He, 12C, 20Ne, 40Ar, 56Fe and 132Xe ions. The results of the in silico calculations were compared with published in vitro data. Main results. All models appear to underestimate the RBE α of hydrogen ions. All MKMs generally overestimate the RBE50%, RBE10% and RBE1% for ions with an LET greater than ∼200 keV μm−1. This overestimation is greater for small surviving fractions and is likely due to the assumption of a radiation-independent quadratic term of clonogenic survival (ß). The overall RBE trends seem to be best described by the novel ‘post-processing average’ implementation of the non-Poisson MKM. In case of calculations with the non-Poisson MKM, pre- or post- processing the average biophysical quantities affects the computed RBE values significantly. Significance. This study presents a systematic analysis of the formalism and results of widely used microdosimetric models of clonogenic survival for ions relevant for cancer particle therapy and space radiation protection. Points for improvements were highlighted and will contribute to the development of upgraded biophysical models.

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

Background

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

Methods and materials

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

Results

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

Discussion

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

Conclusions

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

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-022-02022-5.

Normal tissue complication probability modeling to guide individual treatment planning in pediatric cranial proton and photon radiotherapy

PURPOSE

Proton therapy (PT) is broadly accepted as the gold standard of care for pediatric patients with cranial cancer. The superior dose distribution of PT compared to photon radiotherapy reduces normal tissue complication probability (NTCP) for organs at risk. As NTCPs for pediatric organs are not well understood, clinics generally base radiation response on adult data. However, there is evidence that radiation response strongly depends on the age and even sex of a patient. Furthermore, questions surround the influence of individual intrinsic radiosensitivity (α/β ratio) on pediatric NTCP. While the clinical pediatric NTCP data is scarce, radiobiological modeling and sensitivity analyses can be used to investigate the NTCP trends and its dependence on individual modeling parameters. The purpose of this study was to perform sensitivity analyses of NTCP models to ascertain the dependence of radiosensitivity, sex, and age of a child and predict cranial side-effects following intensity-modulated proton therapy (IMPT) and intensity-modulated radiotherapy (IMRT).

METHODS

Previously, six sex-matched pediatric cranial datasets (5, 9, and 13 years old) were planned in Varian Eclipse treatment planning system (13.7). Up to 108 scanning beam IMPT plans and 108 IMRT plans were retrospectively optimized for a range of simulated target volumes and locations. In this work, dose-volume histograms were extracted and imported into BioSuite Software for radiobiological modeling. Relative-Seriality and Lyman-Kutcher-Burman models were used to calculate NTCP values for toxicity endpoints, where TD50, (based on reported adult clinical data) was varied to simulate sex dependence of NTCP. Plausible parameter ranges, based on published literature for adults, were used in modeling. In addition to sensitivity analyses, a 20% difference in TD50 was used to represent the radiosensitivity between the sexes (with females considered more radiosensitive) for ease of data comparison as a function of parameters such as α/β ratio.

RESULTS

IMPT plans resulted in lower NTCP compared to IMRT across all models (p < 0.0001). For medulloblastoma treatment, the risk of brainstem necrosis (> 10%) and cochlea tinnitus (> 20%) among females could potentially be underestimated considering a lower TD50 value for females. Sensitivity analyses show that the difference in NTCP between sexes was significant (p < 0.0001). Similarly, both brainstem necrosis and cochlea tinnitus NTCP varied significantly (p < 0.0001) across tested α/β as a function of TD50 values (assumption being that TD50 values are 20% lower in females). If the true α/β of these pediatric tissues is higher than expected (α/β ∼ 3), the risk of tinnitus for IMRT can significantly increase (p < 0.0001).

CONCLUSION

Due to the scarcity of pediatric NTCP data available, sensitivity analyses were performed using plausible ranges based on published adult data. In the clinical scenario where, if female pediatric patients were 20% more radiosensitive (lower TD50 value), they could be up to twice as likely to experience side-effects of brainstem necrosis and cochlea tinnitus compared to males, highlighting the need for considering the sex in NTCP models. Based on our sensitivity analyses, age and sex of a pediatric patient could significantly affect the resultant NTCP from cranial radiotherapy, especially at higher α/β values.

Cognitive Functions of Pediatric Brain Tumor Survivors Treated With Proton Beam Therapy: A Case Series

Pediatric brain tumor survivors who received proton beam therapy at the University of Tsukuba Hospital from 2004 to 2011 were retrospectively evaluated for cognitive function. Five patients were included. The median age of diagnosis was 5.4 years (range: 1.5 to 12.5 y) and the median follow-up time was 5.8 years (range: 3.1 to 8.1 y). IQ scores at follow-up were decreased in 2 of 5 patients; 1 underwent whole-brain irradiation and the other was examined just after surgical removal of recurrent tumors. Local proton beam therapy may preserve cognitive function in survivors of pediatric brain tumors.

Postsurgical geometrical variations of tumor bed and brainstem during photon and proton therapy for pediatric tumors of the posterior fossa: dosimetric impact and predictive factors

PURPOSE

Brainstem radionecrosis is an important issue during the irradiation of tumors of the posterior fossa. The aim of the present study is to analyze postsurgical geometrical variations of tumor bed (TB) and brainstem (BS) and their impact on dosimetry.

METHODS

Retrospective collection of data from pediatric patients treated at a single institution. Availability of presurgical magnetic resonance imaging (MRI) was verified; availability of at least two postsurgical MRIs was considered a further inclusion criterion. The following metrics were analyzed: total volume, Dice similarity coefficient (DSC), and Haudsdorff distances (HD).

RESULTS

Fourteen patients were available for the quantification of major postsurgical geometrical variations of TB. DSC, HD max, and HD average values were 0.47 (range: 0.08;0.76), 11.3 mm (7.7;24.5), and 2.6 mm (0.7;6.7) between the first and the second postoperative MRI, respectively. Postsurgical geometrical variations of the BS were also observed. Coverage to the TB was reduced in one patient (D95: -2.9 Gy), while D2 to the BS was increased for the majority of patients. Overall, predictive factors for significant geometrical changes were presurgical gross tumor volume (GTV) > 33 mL, hydrocephaly at diagnosis, Luschka foramen involvement, and younger age (≤ 8 years).

CONCLUSION

Major volume changes were observed in this cohort, with some dosimetric impact. The use of a recent co-registration MRI is advised. The 2-3 mm HD average observed should be considered in the planning target volume/planning organ at risk volume (PTV/PRV) margin and/or robust optimization planning. Results from wider efforts are needed to verify our findings.

Clinical Progress in Proton Radiotherapy: Biological Unknowns

Simple Summary

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

Abstract

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

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.

Technical challenges for FLASH proton therapy

There is growing interest in the radiotherapy community in the application of FLASH radiotherapy, wherein the dose is delivered to the entire treatment volume in less than a second. Early pre-clinical evidence suggests that these extremely high dose rates provide significant sparing of healthy tissue compared to conventional radiotherapy without reducing the damage to cancerous cells. This interest has been reflected in the proton therapy community, with early tests indicating that the FLASH effect is also present with high dose rate proton irradiation. In order to deliver clinically relevant doses at FLASH dose rates significant technical hurdles must be overcome in the accelerator technology before FLASH proton therapy can be realised. Of these challenges, increasing the average current from the present clinical range of 1-10 nA to in excess of 100 nA is at least feasible with existing technology, while the necessity for rapid energy adjustment on the order of a few milliseconds is much more challenging, particularly for synchrotron-based systems. However, the greatest challenge is to implement full pencil beam scanning, where scanning speeds 2 orders of magnitude faster than the existing state-of-the-art will be necessary, along with similar improvements in the speed and accuracy of associated dosimetry. Hybrid systems utilising 3D-printed patient specific range modulators present the most likely route to clinical delivery. However, to correctly adapt and develop existing technology to meet the challenges of FLASH, more pre-clinical studies are needed to properly establish the beam parameters that are necessary to produce the FLASH effect.

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.

Radiobiological model-based approach to determine the potential of dose-escalated robust intensity-modulated proton radiotherapy in reducing gastrointestinal toxicity in the treatment of locally advanced unresectable pancreatic cancer of the head

Background

The purpose of this study was to determine the potential of escalated dose radiation (EDR) robust intensity-modulated proton radiotherapy (ro-IMPT) in reducing GI toxicity risk in locally advanced unresectable pancreatic cancer (LAUPC) of the head in term of normal tissue complication probability (NTCP) predictive model.

Methods

For 9 patients, intensity-modulated radiotherapy (IMRT) was compared with ro-IMPT. For all plans, the prescription dose was 59.4GyE (Gray equivalent) in 33 fractions with an equivalent organ at risk (OAR) constraints. Physical dose distribution was evaluated. GI toxicity risk for different endpoints was estimated using published NTCP Lyman Kutcher Burman (LKB) models for stomach, duodenum, small bowel, and combine stomach and duodenum (Stoduo). A Wilcoxon signed-rank test was used for dosimetry parameters and NTCP values comparison.

Result

The dosimetric results have shown that, with similar target coverage, ro-IMPT achieves a significant dose-volume reduction in the stomach, small bowel, and stoduo in low to high dose range in comparison to IMRT. NTCP evaluation for the endpoint gastric bleeding of stomach (10.55% vs. 13.97%, P = 0.007), duodenum (1.87% vs. 5.02%, P = 0.004), and stoduo (5.67% vs. 7.81%, P = 0.008) suggest reduced toxicity by ro-IMPT compared to IMRT. ∆NTCP _{IMRT – ro-IMPT} (using parameter from Pan et al. for gastric bleed) of ≥5 to < 10% was seen in 3 patients (33%) for stomach and 2 patients (22%) for stoduo. An overall GI toxicity relative risk (NTCP_ro-IMPT/NTCP_IMRT) reduction was noted (0.16–0.81) for all GI-OARs except for duodenum (> 1) with endpoint grade ≥ 3 GI toxicity (using parameters from Holyoake et al.).

Conclusion

With similar target coverage and better conformity, ro-IMPT has the potential to substantially reduce the risk of GI toxicity compared to IMRT in EDR of LAUPC of the head. This result needs to be further evaluated in future clinical studies.

Investigating the impact of alpha/beta and LETd on relative biological effectiveness in scanned proton beams: An in vitro study based on human cell lines

PURPOSE

A relative biological effectiveness (RBE) of 1.1 is commonly used in clinical proton therapy, irrespective of tissue type and depth. This in vitro study was conducted to quantify the RBE of scanned protons as a function of the dose-averaged linear energy transfer (LET_d ) and the sensitivity factor (α/ß)_X . Additionally, three phenomenological models (McNamara, Rørvik, and Jones) and one mechanistic model (repair-misrepair-fixation, RMF) were applied to the experimentally derived data.

METHODS

Four human cell lines (FaDu, HaCat, Du145, SKMel) with differential (α/ß)_X ratios were irradiated in a custom-designed irradiation setup with doses between 0 and 6 Gy at proximal, central, and distal positions of a 80 mm spread-out Bragg peak (SOBP) centered at 80 mm (setup A: proton energies 66.5-135.6 MeV) and 155 mm (setup B: proton energies 127.2-185.9 MeV) depth, respectively. LET_d values at the respective cell positions were derived from Monte Carlo simulations performed with the treatment planning system (TPS, RayStation). Dosimetric measurements were conducted to verify dose homogeneity and dose delivery accuracy. RBE values were derived for doses that resulted in 90 % (RBE₉₀ ) and 10 % (RBE₁₀ ) of cell survival, and survival after a 0.5 Gy dose (RBE_0.5Gy ), 2 Gy dose (RBE_2Gy ), and 6 Gy dose (RBE_6Gy ).

RESULTS

LET_d values at sample positions were 1.9, 2.1, 2.5, 2.8, 4.1, and 4.5 keV/µm. For the cell lines with high (α/ß)_X ratios (FaDu, HaCat), the LET_d did not impact on the RBE. For low (α/ß)_X cell lines (Du145, SKMel), LQ-derived survival curves indicated a clear correlation of LET_d and RBE. RBE₉₀ values up to 2.9 and RBE₁₀ values between 1.4 and 1.8 were obtained. Model-derived RBE predictions slightly overestimated the RBE for the high (α/ß)_X cell lines, although all models except the Jones model provided RBE values within the experimental uncertainty. For low (α/ß)_X cell lines, no agreement was found between experiments and model predictions, that is, all models underestimated the measured RBE.

CONCLUSIONS

The sensitivity parameter (α/ß)_X was observed to be a major influencing factor for the RBE of protons and its sensitivity toward LET_d changes. RBE prediction models are applicable for high (α/ß)_X cell lines but do not estimate RBE values with sufficient accuracy in low (α/ß)_X cell lines.

Dose-Volume and Radiobiological Model-Based Comparative Evaluation of the Gastrointestinal Toxicity Risk of Photon and Proton Irradiation Plans in Localized Pancreatic Cancer Without Distant Metastasis

Background: Radiobiological model-based studies of photon-modulated radiotherapy for pancreatic cancer have reported reduced gastrointestinal (GI) toxicity, although the risk is still high. The purpose of this study was to investigate the potential of 3D-passive scattering proton beam therapy (3D-PSPBT) in limiting GI organ at risk (OAR) toxicity in localized pancreatic cancer based on dosimetric data and the normal tissue complication probability (NTCP) model. Methods: The data of 24 pancreatic cancer patients were retrospectively analyzed, and these patients were planned with intensity-modulated radiotherapy (IMRT), volume-modulated arc therapy (VMAT), and 3D-PSPBT. The tumor was targeted without elective nodal coverage. All generated plans consisted of a 50.4-GyE (Gray equivalent) dose in 28 fractions with equivalent OAR constraints, and they were normalized to cover 50% of the planning treatment volume (PTV) with 100% of the prescription dose. Physical dose distributions were evaluated. GI-OAR toxicity risk for different endpoints was estimated by using published NTCP Lyman-Kutcher-Burman (LKB) models. Analysis of variance (ANOVA) was performed to compare the dosimetric data, and ΔNTCP_IMRT-PSPBT and ΔNTCP_VMAT-PSPBT were also computed. Results: Similar homogeneity and conformity for the clinical target volume (CTV) and PTV were exhibited by all three planning techniques (P > 0.05). 3D-PSPBT resulted in a significant dose reduction for GI-OARs in both the low-intermediate dose range (below 30 GyE) and the highest dose region (D _max and V _{50 GyE}) in comparison with IMRT and VMAT (P < 0.05). Based on the NTCP evaluation, the NTCP reduction for GI-OARs by 3D-PSPBT was minimal in comparison with IMRT and VMAT. Conclusion: 3D-PSPBT results in minimal NTCP reduction and has less potential to substantially reduce the toxicity risk of upper GI bleeding, ulceration, obstruction, and perforation endpoints compared to IMRT and VMAT. 3D-PSPBT may have the potential to reduce acute dose-limiting toxicity in the form of nausea, vomiting, and diarrhea by reducing the GI-OAR treated volume in the low-to-intermediate dose range. However, this result needs to be further evaluated in future clinical studies.

Radiobiological and dosimetric impact of RayStation pencil beam and Monte Carlo algorithms on intensity-modulated proton therapy breast cancer plans

PURPOSE

RayStation treatment planning system employs pencil beam (PB) and Monte Carlo (MC) algorithms for proton dose calculations. The purpose of this study is to evaluate the radiobiological and dosimetric impact of RayStation PB and MC algorithms on the intensity-modulated proton therapy (IMPT) breast plans.

METHODS

The current study included ten breast cancer patients, and each patient was treated with 1-2 proton beams to the whole breast/chestwall (CW) and regional lymph nodes in 28 fractions for a total dose of 50.4 Gy relative biological effectiveness (RBE). A total clinical target volume (CTV_Total) was generated by combining individual CTVs: AxI, AxII, AxIII, CW, IMN, and SCVN. All beams in the study were treated with a range shifter (7.5 cm water equivalent thickness). For each patient, three sets of plans were generated: (a) PB optimization followed by PB dose calculation (PB-PB), (b) PB optimization followed by MC dose calculation (PB-MC), and (c) MC optimization followed by MC dose calculation (MC-MC). For a given patient, each plan was robustly optimized on the CTVs with same parameters and objectives. Treatment plans were evaluated using dosimetric and radiobiological indices (equivalent uniform dose (EUD), tumor control probability (TCP), and normal tissue complication probability (NTCP)).

RESULTS

The results are averaged over ten breast cancer patients. In comparison to PB-PB plans, PB-MC plans showed a reduction in CTV target dose by 5.3% for D_99% and 4.1% for D_95% , as well as a reduction in TCP by 1.5-2.1%. Similarly, PB overestimated the EUD of target volumes by 1.8─3.2 Gy(RBE). In contrast, MC-MC plans achieved similar dosimetric and radiobiological (EUD and TCP) results as the ones in PB-PB plans. A selection of one dose calculation algorithm over another did not produce any noticeable differences in the NTCP of the heart, lung, and skin.

CONCLUSION

If MC is more accurate than PB as reported in the literature, dosimetric and radiobiological results from the current study suggest that PB overestimates the target dose, EUD, and TCP for IMPT breast cancer treatment. The overestimation of dosimetric and radiobiological results of the target volume by PB needs to be further interpreted in terms of clinical outcome.

Effect of Irradiation Time on Biological Effectiveness and Tumor Control Probability in Proton Therapy

PURPOSE

The biological effectiveness of proton beams may decrease with irradiation time because of sublethal damage repair (SLDR). The purpose of this study is to systematically evaluate this effect in hypofractionated proton therapy for various target sizes, depths, and prescribed doses per fraction.

METHODS AND MATERIALS

Plans with a single spread-out Bragg peak beam were created using a constant relative biological effectiveness (RBE) of 1.1 to cover targets of 6 different sizes located at 3 different depths in water. Biological doses of 2, 3, 5, 10, and 20 Gy (RBE) were prescribed to the targets. First, to investigate the depth variation of the biological effectiveness, the biological dose in instantaneous irradiation was recalculated based on the microdosimetric kinetic model. SLDR was then taken into account in the microdosimetric kinetic model during treatments to obtain the irradiation time-dependent biological effectiveness for irradiation time T of 5 to 60 minutes and beam interruption time τ of 0 to 60 minutes. The tumor control probabilities were calculated for single-fraction proton therapy fields of different Ts and τs, and the curative doses were evaluated at a tumor control probability of 90%.

RESULTS

The biological effectiveness decreased with longer T and τ and higher prescribed dose. The maximum decrease in the biological effectiveness was 21% with a 20 Gy (RBE) prescribed dose. In single-fraction proton therapy, the curative dose increased linearly by approximately 33% to 35% with the increase of T from 0 to 60 minutes.

CONCLUSIONS

The biological effectiveness varies largely with T and τ because of SLDR during treatments. This effect was pronounced for high prescribed doses per fraction. Thus, the effect of SLDR needs to be considered in hypofractionated and single-fraction proton therapies in relation to size and depth of the target.

Impact of potentially variable RBE in liver proton therapy

Currently, the relative biological effectiveness (RBE) is assumed to be constant with a value of 1.1 in proton therapy. Although trends of RBE variations are well known, absolute values in patients are associated with considerable uncertainties. This study aims to evaluate the impact of a variable proton RBE in proton therapy liver trials using different fractionation schemes. Sixteen liver cancer cases were evaluated assuming two clinical schedules of 40 Gy/5 fractions and 58.05 Gy/15 fractions. The linear energy transfer (LET) and physical dose distribution in patients were simulated using Monte Carlo. The variable RBE distribution was calculated using a phenomenological model, considering the influence of the LET, fraction size and α/β value. Further, models to predict normal tissue complication probability (NTCP) and tumor control probability (TCP) were used to investigate potential RBE effects on outcome predictions. Applying the variable RBE model to the 5 and 15 fractions schedules results in an increase in mean fraction-size equivalent dose (FED) to the normal liver of 5.0% and 9.6% respectively. For patients with a mean FED to the normal liver larger than 29.8 Gy, this results in a non-negligible increase in the predicted NTCP of the normal liver averaging 11.6%, ranging from 2.7% to 25.6%. On the other hand, decrease in TCP was less than 5% for both fractionation regimens for all patients when assuming a variable RBE instead of constant. Consequently, the difference in TCP between the two fractionation schedules did not change significantly assuming a variable RBE while the impact on the NTCP difference was highly case specific. In addition, both the NTCP and TCP decrease with increasing α/β value for both fractionation schemes, with the decreases being more pronounced when using a variable RBE compared to using RBE  =  1.1. Assuming a constant RBE of 1.1 most likely overestimates the therapeutic ratio in proton therapy for liver cancer, predominantly due to underestimation of the RBE-weighted dose to the normal liver. The impact of applying a variable RBE (as compared to RBE  =  1.1) on the NTCP difference of the two fractionation regimens is case dependent. A variable RBE results in a slight increase in TCP difference. Variations in patient radiosensitivity increase when using a variable RBE.

Proton Relative Biological Effectiveness - Uncertainties and Opportunities

Proton therapy treatments are prescribed using a biological effectiveness relative to photon therapy of 1.1, that is, proton beams are considered to be 10% more biologically effective. Debate is ongoing as to whether this practice needs to be revised. This short review summarizes current knowledge on relative biological effectiveness variations and uncertainties in vitro and in vivo. Clinical relevance is discussed and strategies toward biologically guided treatment planning are presented.

Proton relative biological effectiveness (RBE): a multiscale problem

Proton radiotherapy is undergoing rapid expansion both within the UK and internationally, but significant challenges still need to be overcome if maximum benefit is to be realised from this technique. One major limitation is the persistent uncertainty in proton relative biological effectiveness (RBE). While RBE values are needed to link proton radiotherapy to our existing experience with photon radiotherapy, RBE remains poorly understood and is typically incorporated as a constant dose scaling factor of 1.1 in clinical plans. This is in contrast to extensive experimental evidence indicating that RBE is a function of dose, tissue type, and proton linear energy transfer, among other parameters. In this article, we discuss the challenges associated with obtaining clinically relevant values for proton RBE through commonly-used assays, and highlight the wide range of other experimental end points which can inform our understanding of RBE. We propose that accurate and robust optimization of proton radiotherapy ultimately requires a multiscale understanding of RBE, integrating subcellular, cellular, and patient-level processes.