Effects of Proton Craniospinal Radiation on Vertebral Body Growth Retardation in Children.

PURPOSE: It is of great interest to physicians and patients/patients' families to be able to predict the amount of growth decrement after craniospinal irradiation (CSI). Little data exist on the effect of proton CSI. Our aim was to determine the effect of proton CSI on vertebral body (VB) growth retardation, and to identify factors associated with growth delay. METHODS AND MATERIALS: We performed a retrospective outcome data analysis of 80 patients <16 years old with central nervous system tumors who received proton radiation therapy (PRT) at the Massachusetts General Hospital between 2002 and 2010 with available spinal magnetic resonance imaging. Forty-eight patients received CSI, and 32 patients with brain tumors who received focal cranial irradiation served as controls. VB height was measured midline using sagittal T1-weighted contrast or noncontrast enhanced magnetic resonance imaging of the spine. Measurements were repeated at multiple levels (C3, C3-C4, T4, T4-T5, C3-T6, T4-T7, L3, L1-L5) on available scans for the duration of follow-up. Data were fitted using a mixed-effects multivariable regression model, including follow-up time, CSI dose, age at CSI, and pretreatment VB percentile as parameters. RESULTS: Median follow-up was 69.6 months for patients treated with proton CSI and 52.9 months for the control group. There was a significant association of CSI dose, follow-up time, age at treatment, and pretreatment VB percentile with VB growth retardation. Growth retardation was shown to be independent of gender or growth hormone deficiency. CONCLUSIONS: Although the current practice of PRT CSI delivery allows for sparing of the organs anterior to the spine, the vertebral column receives radiation therapy because of its close proximity to the targeted spinal canal. In growing children, the whole VB has generally been included so that growth impairment is even across the VB. We present a quantitative model predicting the growth retardation of patients treated with PRT CSI based on age at treatment, CSI dose, follow-up time, and pretreatment growth percentile.

Tattoo-Free Setup for Patients With Breast Cancer Receiving Regional Nodal Irradiation.

PURPOSE: Patients undergoing regional nodal irradiation (RNI) with either 3-dimensional conformal radiation therapy (3DCRT) planning or volumetric modulated arc therapy (VMAT) receive permanent tattoos to assist with daily setup alignment and verification. With the advent of surface imaging, tattoos may not be necessary to ensure setup accuracy. We compared the accuracy of conventional tattoo-based setups to those without reference to tattoos. METHODS AND MATERIALS: Forty-eight patients receiving RNI at our institution from July 2019 to December 2020 were identified. All patients received tattoos per standard of care. Twenty-four patients underwent setup using tattoos for initial positioning followed by surface and x-ray imaging. A subsequent 24 patients underwent positioning using surface imaging followed by x-ray imaging without reference to tattoos. Patient cohorts were balanced by treatment technique and use of deep inspiration breath hold. Treatment (including setup and delivery) time and x-ray-based shifts after surface imaging were recorded. RESULTS: Among patients in the tattoo group receiving 3DCRT RNI, the average treatment time per fraction was 21.35 versus 19.75 minutes in the 3DCRT RNI no-tattoo cohort (P = .03). Mean 3D vector shifts for patients in the tattoo cohort were 5.6 versus 4.4 mm in the no-tattoo cohort. The average treatment time per fraction for the tattoo VMAT RNI cohort was 23.16 versus 20.82 minutes in the no-tattoo VMAT RNI cohort (P = .08). Mean 3D vector shifts for the patients in the tattoo VMAT cohort were 5.5 versus 7.1 mm in the no-tattoo VMAT cohort. Breath hold technique and body mass index did not affect accuracy in a consistent or clinically relevant way. CONCLUSIONS: Using a combination of surface and x-ray imaging, without reference to tattoos, provides excellent accuracy in alignment and setup verification among patients receiving RNI for breast cancer, regardless of treatment technique and with reduced treatment time. Skin-based tattoos are no longer warranted for patients receiving supine RNI.

Radiation Necrosis with Proton Therapy in a Patient with Aarskog-Scott Syndrome and Medulloblastoma.

PURPOSE: Medulloblastoma is known to be associated with multiple cancer-predisposition syndromes. In this article, we explore a possible association among a patient's Aarskog-Scott syndrome, development of medulloblastoma, and subsequent brainstem radiation necrosis. CASE PRESENTATION: A 5-year-old male with Aarskog-Scott syndrome initially presented to his pediatrician with morning emesis, gait instability, and truncal weakness. He was ultimately found to have a posterior fossa tumor with pathology consistent with group 3 medulloblastoma. After receiving a gross total resection and standard proton beam radiation therapy with concurrent vincristine, he was noted to develop brainstem radiation necrosis, for which he underwent therapy with high-dose dexamethasone, bevacizumab, and hyperbaric oxygen therapy with radiographic improvement and clinical stabilization. CONCLUSION: Based on several possible pathologic correlates in the FDG1 pathway, there exists a potential association between this patient's Aarskog-Scott syndrome and medulloblastoma, which needs to be investigated further. In patients with underlying, rare genetic syndromes, further caution should be taken when evaluating chemotherapy and radiation dosimetry planning.

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.

A dual-stream deep convolutional network for reducing metal streak artifacts in CT images.

Machine learning and deep learning are rapidly finding applications in the medical imaging field. In this paper, we address the long-standing problem of metal artifacts in computed tomography (CT) images by training a dual-stream deep convolutional neural network for streak removal. While many metal artifact reduction methods exist, even state-of-the-art algorithms fall short in some clinical applications. Specifically, proton therapy planning requires high image quality with accurate tumor volumes to ensure treatment success. We explore a dual-stream deep network structure with residual learning to correct metal streak artifacts after a first-pass by a state-of-the-art interpolation-based algorithm, NMAR. We provide the network with a mask of the streaks in order to focus attention on those areas. Our experiments compare a mean squared error loss function with a perceptual loss function to emphasize preservation of image features and texture. Both visual and quantitative metrics are used to assess the resulting image quality for metal implant cases. Success may be due to the duality of information processing, with one network stream performing local structure correction, while the other stream provides an attention mechanism to destreak effectively. This study shows that image-domain deep learning can be highly effective for metal artifact reduction (MAR), and highlights the benefits and drawbacks of different loss functions for solving a major CT reconstruction challenge.

Tattoo free setup for partial breast irradiation: A feasibility study.

PURPOSE: Patients undergoing external beam accelerated partial breast irradiation (APBI) receive permanent tattoos to aid with daily setup alignment and verification. With the advent of three-dimensional (3D) body surface imaging and two-dimensional (2D) x-ray imaging-based matching to surgical clips, tattoos may not be necessary to ensure setup accuracy. We compared the accuracy of conventional tattoo-based setups to a patient setup without tattoos. MATERIALS/METHODS: Twenty consecutive patients receiving APBI at our institution from July 10, 2017 to February 13, 2018 were identified. All patients received tattoos per standard of care. Ten patients underwent setup using tattoos for initial positioning followed by surface imaging and 2D matching of surgical clips. The other ten patients underwent positioning using surface imaging followed by 2D matching without reference to tattoos. Overall setup time and orthogonal x-ray-based shifts after surface imaging per fraction were recorded. Shift data were used to calculate systematic and random error. RESULTS: Among ten patients in the "no tattoo" group, the average setup time per fraction was 6.83 min vs 8.03 min in the tattoo cohort (P < 0.01). Mean 3D vector shifts for patients in the "no tattoo" group were 4.6 vs 5.9 mm in the "tattoo" cohort (P = NS). Mean systematic errors in the "no tattoo" group were: 1.2 mm (1.5 mm SD) superior/inferior, 0.5 mm (1.6 mm SD) right/left, and 2.3 mm (1.9 mm SD) anterior/posterior directions. Mean systematic errors in the "tattoo" group were: 0.8 mm (2.2 mm SD) superior/inferior, 0.3 mm (2.5 mm SD) right/left, and 1.4 mm (4.4 mm SD) anterior/posterior directions. The random errors in the "no tattoo" group ranged from 0.6 to 0.7 mm vs 1.2 to 1.7 mm in the "tattoo" group. CONCLUSIONS: Using both surface imaging and 2D matching to surgical clips provides excellent accuracy in APBI patient alignment and setup verification with reduced setup time relative to the tattoo cohort. Skin-based tattoos may no longer be warranted for patients receiving external beam APBI.

Radiation tolerance of the optic pathway in patients treated with proton and photon radiotherapy.

INTRODUCTION: Radiation-induced optic neuropathy (RION) is a complication of radiation therapy (RT) that causes blindness. We aimed to define the tolerance of the anterior optic pathway to fractionated RT and identify risk factors for RION. MATERIALS/METHODS: Patients with chordoma or chondrosarcoma of the skull base treated with proton and photon therapy between 1983 and 2013, who received a minimum of 30 Gy (relative biologic effectiveness [RBE]) to the anterior optic pathway were assessed. Optic neuropathy with radiographic correlation occurring ≥6 months after completion of RT in the absence of tumor recurrence or other probable cause was diagnosed as RION. RESULTS: Of 514 patients, 17 developed RION. With median follow-up of 4.8 years, cumulative incidence of RION was 1% among patients receiving <59 Gy (RBE) and 5.8% among patients receiving ≥60 Gy (RBE) to the optic pathway. Higher maximum point dose to the optic pathway (subhazard ratio [SHR] = 1.2, 95% CI 1.05-1.2, p = 0.001), older age (SHR = 1.1, 95% CI 1.02-1.08, p < 0.0005), and female sex (SHR = 16.3, 95% CI 2.2-122.4, p = 0.007) were statistically significant risk factors for RION in multivariate analysis. CONCLUSION: In our study cohort, rates of RION were very low with conventionally fractionated RT up to 59 Gy. At doses ≥60 Gy, there is an increased risk of RION, with greater risk for women and older patients.

Relative Biological Effectiveness Uncertainties and Implications for Beam Arrangements and Dose Constraints in Proton Therapy.

Current clinical implementation of proton radiation therapy assumes a constant relative biological effectiveness (RBE) value of 1.1 throughout the treatment field, for both the target and organs at risks. Although few in vivo clinical data suggest that this approximation is clinically significant, in vitro studies demonstrate the dependency of RBE on dose, fractionation, proton energy, and linear energy transfer, as well as patient radiosensitivity and definition of endpoint. This article provides a brief review on the principles and individual factors contributing to RBE uncertainties, with emphasis on clinical practice. Clinical considerations regarding the effect of RBE uncertainties and implications for beam arrangements in proton therapy treatment planning are discussed through clinical examples for treatments of prostate cancer and posterior fossa tumors as well as craniospinal irradiation for medulloblastoma. Approaches on biological optimization in proton therapy are presented, including a discussion on linear energy transfer-based optimization as an alternative method for biological optimization and its implementation both in multicriteria optimization and inverse optimization modules.

Can differences in linear energy transfer and thus relative biological effectiveness compromise the dosimetric advantage of intensity-modulated proton therapy as compared to passively scattered proton therapy?

PURPOSE: To investigate the effect of differences in linear energy transfer (LET) and thus the relative biological effectiveness (RBE) between passively scattered proton therapy (PS) and pencil-beam scanning intensity-modulated proton therapy (IMPT). METHODS: IMPT treatment plans were generated for six ependymoma patients, originally treated with PS, using the original plan's computed tomography image sets and beam directions, and its dose-volume values as optimization constraints. Two beam spot sizes and both single-field optimization (SFO) and multi-field optimization (MFO) techniques were used for each patient. Three-dimensional variable-RBE-weighted dose distributions were computed, using Monte Carlo calculated dose and LET distributions, and a linear dose and LET-based RBE model, and were compared between the two delivery methods. RESULTS: Increased target dose coverage and decreased mean and maximum dose to the OARs was achieved with IMPT compared to PS, for constant RBE value of 1.1. Nevertheless, the maximum variable-RBE-weighted dose to the brainstem, was increased up to 6% for the IMPT plans compared to the corresponding PS plans. CONCLUSIONS: IMPT can be dosimetrically superior to PS for ependymoma patients. However, caution should be exercised so that the increased dose conformity is not counteracted by an increase in radiobiological effect in adjacent critical structures.

Brainstem Injury in Pediatric Patients With Posterior Fossa Tumors Treated With Proton Beam Therapy and Associated Dosimetric Factors.

PURPOSE: Proton radiation therapy is commonly used in young children with brain tumors for its potential to reduce late effects. However, some proton series report higher rates of brainstem injury (0%-16%) than most photon series (2.2%-8.6%). We report the incidence of brainstem injury and a risk factor analysis in pediatric patients with posterior fossa primary tumors treated with proton radiation therapy at our institution. METHODS AND MATERIALS: The study included 216 consecutive patients treated between 2000 and 2015. Dosimetry was available for all but 4 patients. Grade 2 to 5 late brainstem toxicity was assessed by the National Cancer Institute Common Terminology Criteria for Adverse Events version 4.0. RESULTS: The histologies include medulloblastoma (n=154, 71.3%), ependymoma (n=56, 25.9%), and atypical teratoid rhabdoid tumor (n=6, 2.8%). The median age at irradiation was 6.6 years (range, 0.5-23.1 years); median dose, 54 gray relative biological effectiveness (Gy RBE) (range, 46.8-59.4 Gy RBE); and median follow-up period, 4.2 years (range, 0.1-15.3 years) among 198 survivors. Of the patients, 83.3% received chemotherapy; 70.4% achieved gross total resection. The crude rate of injury was 2.3% in all patients, 1.9% in those with medulloblastoma, 3.6% in those with ependymoma, and 0% in those with atypical teratoid rhabdoid tumor. The 5-year cumulative incidence of injury was 2.0% (95% confidence interval, 0.7%-4.8%). The median brainstem dose (minimum dose received by 50% of brainstem) in the whole cohort was 53.6 Gy RBE (range, 16.5-56.8 Gy RBE); maximum point dose within the brainstem (Dmax), 55.2 Gy RBE (range, 48.4-60.5 Gy RBE); and mean dose, 50.4 Gy RBE (range, 21.1-56.7 Gy RBE). In the 5 patients with injury, the median minimum dose received by 50% of the brainstem was 54.6 Gy RBE (range, 50.2-55.1 Gy RBE); Dmax, 56.2 Gy RBE (range, 55.0-57.1 Gy RBE); mean dose, 51.3 Gy RBE (range, 45.4-54.4 Gy RBE); and median volume of the brainstem receiving ≥55 Gy RBE (V55), 27.4% (range, 0%-59.4%). Of the 5 patients with injury, 4 had a brainstem Dmax in the highest quartile (≥55.8 Gy RBE, P = .016) and a V55 in the highest tertile (>6.0%) of the cohort distribution (P = .047). Of the 5 patients with injury, 3 were aged >6 years (age range, 4.1-22.8 years), and 4 of 5 patients received chemotherapy and achieved gross total resection. CONCLUSIONS: The incidence of injury in pediatric patients with posterior fossa tumors is consistent with previous reports in the photon setting. Our data suggest that when Dmax and V55 are kept <55.8 Gy RBE and ≤6.0%, respectively, the 5-year rate of radiation brainstem injury would be <2%.

Evaluating Intensity Modulated Proton Therapy Relative to Passive Scattering Proton Therapy for Increased Vertebral Column Sparing in Craniospinal Irradiation in Growing Pediatric Patients.

PURPOSE: At present, proton craniospinal irradiation (CSI) for growing children is delivered to the whole vertebral body (WVB) to avoid asymmetric growth. We aimed to demonstrate the feasibility and potential clinical benefit of delivering vertebral body sparing (VBS) versus WVB CSI with passively scattered (PS) and intensity modulated proton therapy (IMPT) in growing children treated for medulloblastoma. METHODS AND MATERIALS: Five plans were generated for medulloblastoma patients, who had been previously treated with CSI PS proton radiation therapy: (1) single posteroanterior (PA) PS field covering the WVB (PS-PA-WVB); (2) single PA PS field that included only the thecal sac in the target volume (PS-PA-VBS); (3) single PA IMPT field covering the WVB (IMPT-PA-WVB); (4) single PA IMPT field, target volume including thecal sac only (IMPT-PA-VBS); and (5) 2 posterior-oblique (-35°, +35°) IMPT fields, with the target volume including the thecal sac only (IMPT2F-VBS). For all cases, 23.4 Gy (relative biologic effectiveness [RBE]) was prescribed to 95% of the spinal canal. The dose, linear energy transfer, and variable-RBE-weighted dose distributions were calculated for all plans using the tool for particle simulation, version 2, Monte Carlo system. RESULTS: IMPT VBS techniques efficiently spared the anterior vertebral bodies (AVBs), even when accounting for potential higher variable RBE predicted by linear energy transfer distributions. Assuming an RBE of 1.1, the V10 Gy(RBE) decreased from 100% for the WVB techniques to 59.5% to 76.8% for the cervical, 29.9% to 34.6% for the thoracic, and 20.6% to 25.1% for the lumbar AVBs, and the V20 Gy(RBE) decreased from 99.0% to 17.8% to 20.0% for the cervical, 7.2% to 7.6% for the thoracic, and 4.0% to 4.6% for the lumbar AVBs when IMPT VBS techniques were applied. The corresponding percentages for the PS VBS technique were higher. CONCLUSIONS: Advanced proton techniques can sufficiently reduce the dose to the vertebral body and allow for vertebral column growth for children with central nervous system tumors requiring CSI. This was true even when considering variable RBE values. A clinical trial is planned for VBS to the thoracic and lumbosacral spine in growing children.

Metal artifacts in computed tomography for radiation therapy planning: dosimetric effects and impact of metal artifact reduction.

A significant and increasing number of patients receiving radiation therapy present with metal objects close to, or even within, the treatment area, resulting in artifacts in computed tomography (CT) imaging, which is the most commonly used imaging method for treatment planning in radiation therapy. In the presence of metal implants, such as dental fillings in treatment of head-and-neck tumors, spinal stabilization implants in spinal or paraspinal treatment or hip replacements in prostate cancer treatments, the extreme photon absorption by the metal object leads to prominent image artifacts. Although current CT scanners include a series of correction steps for beam hardening, scattered radiation and noisy measurements, when metal implants exist within or close to the treatment area, these corrections do not suffice. CT metal artifacts affect negatively the treatment planning of radiation therapy either by causing difficulties to delineate the target volume or by reducing the dose calculation accuracy. Various metal artifact reduction (MAR) methods have been explored in terms of improvement of organ delineation and dose calculation in radiation therapy treatment planning, depending on the type of radiation treatment and location of the metal implant and treatment site. Including a brief description of the available CT MAR methods that have been applied in radiation therapy, this article attempts to provide a comprehensive review on the dosimetric effect of the presence of CT metal artifacts in treatment planning, as reported in the literature, and the potential improvement suggested by different MAR approaches. The impact of artifacts on the treatment planning and delivery accuracy is discussed in the context of different modalities, such as photon external beam, brachytherapy and particle therapy, as well as by type and location of metal implants.

Proton Treatment Techniques for Posterior Fossa Tumors: Consequences for Linear Energy Transfer and Dose-Volume Parameters for the Brainstem and Organs at Risk.

PURPOSE: In proton therapy of posterior fossa tumors, at least partial inclusion of the brainstem in the target is necessary because of its proximity to the tumor and required margins. Additionally, the preferred beam geometry results in directing the field distal edge toward this critical structure, raising concerns for brainstem toxicity. Some treatment techniques place the beam's distal edge within the brainstem (dose-sparing techniques), and others avoid elevated linear energy transfer (LET) of the proton field by placing the distal edge beyond it (LET-sparing techniques). Hybrid approaches are also being used. We examine the dosimetric efficacy of these techniques, accounting for LET-dependent and dose-dependent variable relative biologic effectiveness (RBE) distributions. METHODS: Six techniques were applied in ependymoma cases: (a) 3-field dose-sparing; (b) 3-field LET-sparing; (c) 2-field dose-sparing, wide angles; (d) 2-field LET-sparing, wide angles; (e) 2-field LET-sparing, steep angles; and (f) 2-field LET-sparing with feathered distal end. Monte Carlo calculated dose, LET, and RBE-weighted dose distributions were compared. RESULTS: Decreased LET values in the brainstem by LET-sparing techniques were accompanied by higher, not statistically significant, median dose: 53.6 Gy(RBE), 53.4 Gy(RBE), and 54.3 Gy(RBE) for techniques (b), (d), and (e) versus 52.1 Gy(RBE) for technique (a). Accounting for variable RBE distributions, the brainstem volume receiving at least 55 Gy(RBE) increased from 72.5% for technique (a) to 80.3% for (b) (P<.01) and from 70.7% for technique (c) to 77.6% for (d) (P<.01). Less than 2%, but statistically significant, decrease in maximum variable RBE-weighted brainstem dose was observed for the LET-sparing techniques compared with the corresponding dose-sparing (P=.03 and .004). CONCLUSIONS: Extending the proton range beyond the brainstem to reduce LET results in clinically comparable maximum radiobiologic effective dose to this sensitive structure. However this method significantly increasing the brainstem volume receiving RBE-weighted dose higher than 55 Gy(RBE) with possible consequences based on known dose-volume parameters for increased toxicity.

Reoptimization of Intensity Modulated Proton Therapy Plans Based on Linear Energy Transfer.

PURPOSE: We describe a treatment plan optimization method for intensity modulated proton therapy (IMPT) that avoids high values of linear energy transfer (LET) in critical structures located within or near the target volume while limiting degradation of the best possible physical dose distribution. METHODS AND MATERIALS: To allow fast optimization based on dose and LET, a GPU-based Monte Carlo code was extended to provide dose-averaged LET in addition to dose for all pencil beams. After optimizing an initial IMPT plan based on physical dose, a prioritized optimization scheme is used to modify the LET distribution while constraining the physical dose objectives to values close to the initial plan. The LET optimization step is performed based on objective functions evaluated for the product of LET and physical dose (LET×D). To first approximation, LET×D represents a measure of the additional biological dose that is caused by high LET. RESULTS: The method is effective for treatments where serial critical structures with maximum dose constraints are located within or near the target. We report on 5 patients with intracranial tumors (high-grade meningiomas, base-of-skull chordomas, ependymomas) in whom the target volume overlaps with the brainstem and optic structures. In all cases, high LET×D in critical structures could be avoided while minimally compromising physical dose planning objectives. CONCLUSION: LET-based reoptimization of IMPT plans represents a pragmatic approach to bridge the gap between purely physical dose-based and relative biological effectiveness (RBE)-based planning. The method makes IMPT treatments safer by mitigating a potentially increased risk of side effects resulting from elevated RBE of proton beams near the end of range.

Recent developments and comprehensive evaluations of a GPU-based Monte Carlo package for proton therapy.

Monte Carlo (MC) simulation is commonly considered as the most accurate dose calculation method for proton therapy. Aiming at achieving fast MC dose calculations for clinical applications, we have previously developed a graphics-processing unit (GPU)-based MC tool, gPMC. In this paper, we report our recent updates on gPMC in terms of its accuracy, portability, and functionality, as well as comprehensive tests on this tool. The new version, gPMC v2.0, was developed under the OpenCL environment to enable portability across different computational platforms. Physics models of nuclear interactions were refined to improve calculation accuracy. Scoring functions of gPMC were expanded to enable tallying particle fluence, dose deposited by different particle types, and dose-averaged linear energy transfer (LETd). A multiple counter approach was employed to improve efficiency by reducing the frequency of memory writing conflict at scoring. For dose calculation, accuracy improvements over gPMC v1.0 were observed in both water phantom cases and a patient case. For a prostate cancer case planned using high-energy proton beams, dose discrepancies in beam entrance and target region seen in gPMC v1.0 with respect to the gold standard tool for proton Monte Carlo simulations (TOPAS) results were substantially reduced and gamma test passing rate (1%/1 mm) was improved from 82.7%-93.1%. The average relative difference in LETd between gPMC and TOPAS was 1.7%. The average relative differences in the dose deposited by primary, secondary, and other heavier particles were within 2.3%, 0.4%, and 0.2%. Depending on source proton energy and phantom complexity, it took 8-17 s on an AMD Radeon R9 290x GPU to simulate [Formula: see text] source protons, achieving less than [Formula: see text] average statistical uncertainty. As the beam size was reduced from 10 × 10 cm2 to 1 × 1 cm2, the time on scoring was only increased by 4.8% with eight counters, in contrast to a 40% increase using only one counter. With the OpenCL environment, the portability of gPMC v2.0 was enhanced. It was successfully executed on different CPUs and GPUs and its performance on different devices varied depending on processing power and hardware structure.

Can We Advance Proton Therapy for Prostate? Considering Alternative Beam Angles and Relative Biological Effectiveness Variations When Comparing Against Intensity Modulated Radiation Therapy.

PURPOSE: For prostate treatments, robust evidence regarding the superiority of either intensity modulated radiation therapy (IMRT) or proton therapy is currently lacking. In this study we investigated the circumstances under which proton therapy should be expected to outperform IMRT, particularly the proton beam orientations and relative biological effectiveness (RBE) assumptions. METHODS AND MATERIALS: For 8 patients, 4 treatment planning strategies were considered: (A) IMRT; (B) passively scattered standard bilateral (SB) proton beams; (C) passively scattered anterior oblique (AO) proton beams, and (D) AO intensity modulated proton therapy (IMPT). For modalities (B)-(D) the dose and linear energy transfer (LET) distributions were simulated using the TOPAS Monte Carlo platform and RBE was calculated according to 3 different models. RESULTS: Assuming a fixed RBE of 1.1, our implementation of IMRT outperformed SB proton therapy across most normal tissue metrics. For the scattered AO proton plans, application of the variable RBE models resulted in substantial hotspots in rectal RBE weighted dose. For AO IMPT, it was typically not possible to find a plan that simultaneously met the tumor and rectal constraints for both fixed and variable RBE models. CONCLUSION: If either a fixed RBE of 1.1 or a variable RBE model could be validated in vivo, then it would always be possible to use AO IMPT to dose-boost the prostate and improve normal tissue sparing relative to IMRT. For a cohort without rectum spacer gels, this study (1) underlines the importance of resolving the question of proton RBE within the framework of an IMRT versus proton debate for the prostate and (2) highlights that without further LET/RBE model validation, great care must be taken if AO proton fields are to be considered for prostate treatments.

Incidence of CNS Injury for a Cohort of 111 Patients Treated With Proton Therapy for Medulloblastoma: LET and RBE Associations for Areas of Injury.

BACKGROUND: Central nervous system (CNS) injury is a rare complication of radiation therapy for pediatric brain tumors, but its incidence with proton radiation therapy (PRT) is less well defined. Increased linear energy transfer (LET) and relative biological effectiveness (RBE) at the distal end of proton beams may influence this risk. We report the incidence of CNS injury in medulloblastoma patients treated with PRT and investigate correlations with LET and RBE values. METHODS AND MATERIALS: We reviewed 111 consecutive patients treated with PRT for medulloblastoma between 2002 and 2011 and selected patients with clinical symptoms of CNS injury. Magnetic resonance imaging (MRI) findings for all patients were contoured on original planning scans (treatment change areas [TCA]). Dose and LET distributions were calculated for the treated plans using Monte Carlo system. RBE values were estimated based on LET-based published models. RESULTS: At a median follow-up of 4.2 years, the 5-year cumulative incidence of CNS injury was 3.6% for any grade and 2.7% for grade 3+. Three of 4 symptomatic patients were treated with a whole posterior fossa boost. Eight of 10 defined TCAs had higher LET values than the target but statistically nonsignificant differences in RBE values (P=.12). CONCLUSIONS: Central nervous system and brainstem injury incidence for PRT in this series is similar to that reported for photon radiation therapy. The risk of CNS injury was higher for whole posterior fossa boost than for involved field. Although no clear correlation with RBE values was found, numbers were small and additional investigation is warranted to better determine the relationship between injury and LET.

Assessing the Clinical Impact of Approximations in Analytical Dose Calculations for Proton Therapy.

PURPOSE: To assess the impact of approximations in current analytical dose calculation methods (ADCs) on tumor control probability (TCP) in proton therapy. METHODS: Dose distributions planned with ADC were compared with delivered dose distributions as determined by Monte Carlo simulations. A total of 50 patients were investigated in this analysis with 10 patients per site for 5 treatment sites (head and neck, lung, breast, prostate, liver). Differences were evaluated using dosimetric indices based on a dose-volume histogram analysis, a γ-index analysis, and estimations of TCP. RESULTS: We found that ADC overestimated the target doses on average by 1% to 2% for all patients considered. The mean dose, D95, D50, and D02 (the dose value covering 95%, 50% and 2% of the target volume, respectively) were predicted within 5% of the delivered dose. The γ-index passing rate for target volumes was above 96% for a 3%/3 mm criterion. Differences in TCP were up to 2%, 2.5%, 6%, 6.5%, and 11% for liver and breast, prostate, head and neck, and lung patients, respectively. Differences in normal tissue complication probabilities for bladder and anterior rectum of prostate patients were less than 3%. CONCLUSION: Our results indicate that current dose calculation algorithms lead to underdosage of the target by as much as 5%, resulting in differences in TCP of up to 11%. To ensure full target coverage, advanced dose calculation methods like Monte Carlo simulations may be necessary in proton therapy. Monte Carlo simulations may also be required to avoid biases resulting from systematic discrepancies in calculated dose distributions for clinical trials comparing proton therapy with conventional radiation therapy.

Validation of a GPU-based Monte Carlo code (gPMC) for proton radiation therapy: clinical cases study.

Monte Carlo (MC) methods are recognized as the gold-standard for dose calculation, however they have not replaced analytical methods up to now due to their lengthy calculation times. GPU-based applications allow MC dose calculations to be performed on time scales comparable to conventional analytical algorithms. This study focuses on validating our GPU-based MC code for proton dose calculation (gPMC) using an experimentally validated multi-purpose MC code (TOPAS) and compare their performance for clinical patient cases. Clinical cases from five treatment sites were selected covering the full range from very homogeneous patient geometries (liver) to patients with high geometrical complexity (air cavities and density heterogeneities in head-and-neck and lung patients) and from short beam range (breast) to large beam range (prostate). Both gPMC and TOPAS were used to calculate 3D dose distributions for all patients. Comparisons were performed based on target coverage indices (mean dose, V95, D98, D50, D02) and gamma index distributions. Dosimetric indices differed less than 2% between TOPAS and gPMC dose distributions for most cases. Gamma index analysis with 1%/1 mm criterion resulted in a passing rate of more than 94% of all patient voxels receiving more than 10% of the mean target dose, for all patients except for prostate cases. Although clinically insignificant, gPMC resulted in systematic underestimation of target dose for prostate cases by 1-2% compared to TOPAS. Correspondingly the gamma index analysis with 1%/1 mm criterion failed for most beams for this site, while for 2%/1 mm criterion passing rates of more than 94.6% of all patient voxels were observed. For the same initial number of simulated particles, calculation time for a single beam for a typical head and neck patient plan decreased from 4 CPU hours per million particles (2.8-2.9 GHz Intel X5600) for TOPAS to 2.4 s per million particles (NVIDIA TESLA C2075) for gPMC. Excellent agreement was demonstrated between our fast GPU-based MC code (gPMC) and a previously extensively validated multi-purpose MC code (TOPAS) for a comprehensive set of clinical patient cases. This shows that MC dose calculations in proton therapy can be performed on time scales comparable to analytical algorithms with accuracy comparable to state-of-the-art CPU-based MC codes.

Patterns of failure after proton therapy in medulloblastoma; linear energy transfer distributions and relative biological effectiveness associations for relapses.

PURPOSE: The pattern of failure in medulloblastoma patients treated with proton radiation therapy is unknown. For this increasingly used modality, it is important to ensure that outcomes are comparable to those in modern photon series. It has been suggested this pattern may differ from photons because of variations in linear energy transfer (LET) and relative biological effectiveness (RBE). In addition, the use of matching fields for delivery of craniospinal irradiation (CSI) may influence patterns of relapse. Here we report the patterns of failure after the use of protons, compare it to that in the available photon literature, and determine the LET and RBE values in areas of recurrence. METHODS AND MATERIALS: Retrospective review of patients with medulloblastoma treated with proton radiation therapy at Massachusetts General Hospital (MGH) between 2002 and 2011. We documented the locations of first relapse. Discrete failures were contoured on the original planning computed tomography scan. Monte Carlo calculation methods were used to estimate the proton LET distribution. Models were used to estimate RBE values based on the LET distributions. RESULTS: A total of 109 patients were followed for a median of 38.8 months (range, 1.4-119.2 months). Of the patients, 16 experienced relapse. Relapse involved the supratentorial compartment (n=8), spinal compartment (n=11), and posterior fossa (n=5). Eleven failures were isolated to a single compartment; 6 failures in the spine, 4 failures in the supratentorium, and 1 failure in the posterior fossa. The remaining patients had multiple sites of disease. One isolated spinal failure occurred at the spinal junction of 2 fields. None of the 70 patients treated with an involved-field-only boost failed in the posterior fossa outside of the tumor bed. We found no correlation between Monte Carlo-calculated LET distribution and regions of recurrence. CONCLUSIONS: The most common site of failure in patients treated with protons for medulloblastoma was outside of the posterior fossa. The most common site for isolated local failure was the spine. We recommend consideration of spinal imaging in follow-up and careful attention to dose distribution in the spinal junction regions. Development of techniques that do not require field matching may be of benefit. We did not identify a direct correlation between lower LET values and recurrence in medulloblastoma patients treated with proton therapy. Patterns of failure do not appear to differ from those in patients treated with photon therapy.