Inter-patient variations in relative biological efectiveness for craniospinal irradiation with protons

Cranio-spinal irradiation (CSI) using protons has dosimetric advantages compared to photons and is expected to reduce risk of adverse effects. The proton relative biological effectiveness (RBE) varies with linear energy transfer (LET), tissue type and dose, but a variable RBE has not replaced the constant RBE of 1.1 in clinical treatment planning. We examined inter-patient variations in RBE for ten proton CSI patients. Variable RBE models were used to obtain RBE and RBE-weighted doses. RBE was quantified in terms of dose weighted organ-mean RBE (RBEd = mean RBE-weighted dose/mean physical dose) and effective RBE of the near maximum dose (D2%), i.e. RBED2% = D2%,RBE/D2%,phys, where subscripts RBE and phys indicate that the D2% is calculated based on an RBE model and the physical dose, respectively. Compared to the median RBEd of the patient population, differences up to 15% were observed for the individual RBEd values found for the thyroid, while more modest variations were seen for the heart (6%), lungs (2%) and brainstem (<1%). Large inter-patient variation in RBE could be correlated to large spread in LET and dose for these organs at risk (OARs). For OARs with small inter-patient variations, the results show that applying a population based RBE in treatment planning may be a step forward compared to using RBE of 1.1. OARs with large inter-patient RBE variations should ideally be selected for patient-specific biological or RBE robustness analysis if the physical doses are close to known dose thresholds.publishedVersio

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, the RBE varies with factors such as dose level, linear energy transfer (LET) and tissue type. Multiple RBE models have been developed to account for this biological variation. To enable recalculation of patients treated with double scattering (DS) proton therapy, including LET and variable RBE, 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 obtaining LET and variable RBE doses. 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 variable RBE for the two patients, respectively. We presented a method to recalculate DS proton plans in the FLUKA MC code. The implementation was used to obtain LET and variable RBE dose and can be used for investigating variable RBE for previously treated patients.publishedVersio

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 (RBE1.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 RBE1.1 dose characteristics. Differences across cases and controls related to the dose-averaged LET (LETd) and variable RBE-weighted dose from two RBE models were analyzed in the high-dose region. Results LETd metrics were marginally higher for cases vs. controls for the majority of dose levels and brainstem substructures. Considering areas with doses above 54 Gy(RBE1.1), we found a moderate trend of 13% higher median LETd 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 LETd for cases compared to controls were noticeable across structures and LETd metrics for this patient cohort. While case-control differences were minor, an association with the observed symptomatic brainstem toxicity cannot be ruled out.publishedVersio

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.publishedVersio

Combined RBE and OER optimization in proton therapy with FLUKA based on EF5-PET

Introduction Tumor hypoxia is associated with poor treatment outcome. Hypoxic regions are more radioresistant than well-oxygenated regions, as quantified by the oxygen enhancement ratio (OER). In optimization of proton therapy, including OER in addition to the relative biological effectiveness (RBE) could therefore be used to adapt to patient-specific radioresistance governed by intrinsic radiosensitivity and hypoxia. Methods A combined RBE and OER weighted dose (ROWD) calculation method was implemented in a FLUKA Monte Carlo (MC) based treatment planning tool. The method is based on the linear quadratic model, with α and β parameters as a function of the OER, and therefore a function of the linear energy transfer (LET) and partial oxygen pressure (pO2). Proton therapy plans for two head and neck cancer (HNC) patients were optimized with pO2 estimated from [18F]-EF5 positron emission tomography (PET) images. For the ROWD calculations, an RBE of 1.1 (RBE1.1,OER) and two variable RBE models, Rørvik (ROR) and McNamara (MCN), were used, alongside a reference plan without incorporation of OER (RBE1.1). Results For the HNC patients, treatment plans in line with the prescription dose and with acceptable target ROWD could be generated with the established tool. The physical dose was the main factor modulated in the ROWD. The impact of incorporating OER during optimization of HNC patients was demonstrated by the substantial difference found between ROWD and physical dose in the hypoxic tumor region. The largest physical dose differences between the ROWD optimized plans and the reference plan was 12.2 Gy. Conclusion The FLUKA MC based tool was able to optimize proton treatment plans taking the tumor pO2 distribution from hypoxia PET images into account. Independent of RBE-model, both elevated LET and physical dose were found in the hypoxic regions, which shows the potential to increase the tumor control compared to a conventional optimization approach.publishedVersio

The FLUKA Monte Carlo code coupled with an OER model for biologically weighted dose calculations in proton therapy of hypoxic tumors

Introduction The increased radioresistance of hypoxic cells compared to well-oxygenated cells is quantified by the oxygen enhancement ratio (OER). In this study we created a FLUKA Monte Carlo based tool for inclusion of both OER and relative biological effectiveness (RBE) in biologically weighted dose (ROWD) calculations in proton therapy and applied this to explore the impact of hypoxia. Methods The RBE-weighted dose was adapted for hypoxia by making RBE model parameters dependent on the OER, in addition to the linear energy transfer (LET). The OER depends on the partial oxygen pressure (pO2) and LET. To demonstrate model performance, calculations were done with spread-out Bragg peaks (SOBP) in water phantoms with pO2 ranging from strongly hypoxic to normoxic (0.01–30 mmHg) and with a head and neck cancer proton plan optimized with an RBE of 1.1 and pO2 estimated voxel-by-voxel using [18F]-EF5 PET. An RBE of 1.1 and the Rørvik RBE model were used for the ROWD calculations. Results The SOBP in water had decreasing ROWD with decreasing pO2. In the plans accounting for oxygenation, the median target doses were approximately a factor 1.1 lower than the corresponding plans which did not consider the OER. Hypoxia adapted target ROWDs were considerably more heterogeneous than the RBE1.1-weighted doses. Conclusion We realized a Monte Carlo based tool for calculating the ROWD. Read-in of patient pO2 and estimation of ROWD with flexibility in choice of RBE model was achieved, giving a tool that may be useful in future clinical applications of hypoxia-guided particle therapy.publishedVersio

The FLUKA Monte Carlo code coupled with an OER model for biologically weighted dose calculations in proton therapy of hypoxic tumors

Introduction: The increased radioresistance of hypoxic cells compared to well-oxygenated cells is quantified by the oxygen enhancement ratio (OER). In this study we created a FLUKA Monte Carlo based tool for inclusion of both OER and relative biological effectiveness (RBE) in biologically weighted dose (ROWD) calculations in proton therapy and applied this to explore the impact of hypoxia.Methods: The RBE-weighted dose was adapted for hypoxia by making RBE model parameters dependent on the OER, in addition to the linear energy transfer (LET). The OER depends on the partial oxygen pressure (pO(2)) and LET. To demonstrate model performance, calculations were done with spread-out Bragg peaks (SOBP) in water phantoms with pO(2) ranging from strongly hypoxic to normoxic (0.01-30 mmHg) and with a head and neck cancer proton plan optimized with an RBE of 1.1 and pO(2) estimated voxel-by-voxel using [F-18]-EF5 PET. An RBE of 1.1 and the Rorvik RBE model were used for the ROWD calculations.Results: The SOBP in water had decreasing ROWD with decreasing pO(2). In the plans accounting for oxygenation, the median target doses were approximately a factor 1.1 lower than the corresponding plans which did not consider the OER. Hypoxia adapted target ROWDs were considerably more heterogeneous than the RBE1.1-weighted doses.Conclusion: We realized a Monte Carlo based tool for calculating the ROWD. Read-in of patient pO(2) and estimation of ROWD with flexibility in choice of RBE model was achieved, giving a tool that may be useful in future clinical applications of hypoxia-guided particle therapy.</div

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

Purpose Intensity modulated proton therapy (IMPT) could yield high linear energy transfer (LET) in critical structures and increased biological effect. For head and neck cancers at the skull base this could potentially result in radiation-associated brain image change (RAIC). The purpose of the current study was to investigate voxel-wise dose and LET correlations with RAIC after IMPT. Methods and Materials For 15 patients with RAIC after IMPT, contrast enhancement observed on T1-weighted magnetic resonance imaging was contoured and coregistered to the planning computed tomography. Monte Carlo calculated dose and dose-averaged LET (LETd) distributions were extracted at voxel level and associations with RAIC were modelled using uni- and multivariate mixed effect logistic regression. Model performance was evaluated using the area under the receiver operating characteristic curve and precision-recall curve. Results An overall statistically significant RAIC association with dose and LETd was found in both the uni- and multivariate analysis. Patient heterogeneity was considerable, with standard deviation of the random effects of 1.81 (1.30-2.72) for dose and 2.68 (1.93-4.93) for LETd, respectively. Area under the receiver operating characteristic curve was 0.93 and 0.95 for the univariate dose-response model and multivariate model, respectively. Analysis of the LETd effect demonstrated increased risk of RAIC with increasing LETd for the majority of patients. Estimated probability of RAIC with LETd = 1 keV/µm was 4% (95% confidence interval, 0%, 0.44%) and 29% (95% confidence interval, 0.01%, 0.92%) for 60 and 70 Gy, respectively. The TD15 were estimated to be 63.6 and 50.1 Gy with LETd equal to 2 and 5 keV/µm, respectively. Conclusions Our results suggest that the LETd effect could be of clinical significance for some patients; LETd assessment in clinical treatment plans should therefore be taken into consideration.publishedVersio

Measurements of the Photonuclear Neutron Yield During Radiotherapy Using Bubble Detectors and Thermoluminescence Detectors

Medical linear accelerators (linacs) used for radiotherapy with photons produce an undesirable neutron contribution as treatment energies exceed the (γ,n) threshold of the high-Z constituents of the linac. Secondary neutrons in radiotherapy contribute to additional dose outside the treatment volume in the patient. Treatment techniques offering improved confinement of the dose to the target volume are rapidly evolving and will contribute differently to the production of neutron doses outside the volume at target. The intention in this work has been two-fold: (1) Investigation of the characteristics and the responses of the detectors applied. (2) Measurements of the secondary neutron yield outside the treatment volume of a 15 MV photon beam produced by a Varian 23iX medical linac at Haukeland University Hospital (HUS). Neutron yields measured from the treatment techniques Three-Dimensional Conventional Radiation Therapy (3D-CRT) are compared to Intensity Modulated Radiation Therapy (IMRT). The neutron measurements were performed using bubble detectors and Thermoluminescence Detectors (TLDs). The TLD-600 is photon and thermal neutron sensitive, and the TLD-700 is photon sensitive. Pairs of TLDs were applied together in order to discriminate against the photon component of the mixed field. The bubble detectors employed were the Bubble Detector Spectrometer (BDS) used for obtaining information of the neutron energy spectrum, and the Bubble Detector Thermal (BDT) sensitive to thermal neutrons only. The response of the detectors was also measured in a mixed photon thermal-neutron field of the nuclear reactor at Institute for Energy Technology (IFE) at Kjeller (Norway). The BDS detectors were found to be sensitive in this field, either to neutrons with lower energy than the discrimination thresholds given from the vendor or to the photon component of the field. The BDT sensitivity measured in this field was four times higher than the calibrated sensitivity supplied by the manufacturer. The TLDs applied in the reactor field revealed a strong photon background, which made reliable photon discrimination in the TLD signal a challenge. The spatial distribution of the neutrons produced by the medical linac was also measured outside the target volume with TLDs positioned in a plastic phantom. The neutrons were principally detected in the outer layers of the phantom, close to the linac head where the neutrons are produced. Two dose plans were created (in the Eclipse treatment planning system) in order to compare the neutron yield during the IMRT and 3D-CRT treatment techniques. Measurements were performed during delivery of the plans in a solid state plastic phantom, with TLD and BDS detectors situated outside the primary photon fields. The IMRT dose plan resulted in a higher neutron dose than the dose produced during 3D-CRT. The ratio of the neutron dose produced during the delivery of the two plans scales roughly with the ratio of the radiation output of the linac; the Monitor Units (MUs). The neutron fluences measured were 107 n/cm² per treatment Gy delivered to the target volume. For a full treatment of 70 Gy, the additional neutron doses were 0.2±0.1 Sv for 3D-CRT and 0.4±0.3 Sv IMRT. Two characteristic features of the neutron energy spectrum were observed; the low energy continuous distribution from neutron evaporation processes and a rather distinct peak from direct neutron knock-out processes. The sensitivity of the detectors to neutron energies, and potentially to photons, should be further investigated with respect to the correspondence between the response to fluence and dose. To obtain further information of the neutron energy spectrum, the use of TLDs with moderators of different sizes can be an alternative method