Dosimetric benefits of robust treatment planning for intensity modulated proton therapy for base-of-skull cancers

Robust IMPT and follow-up toxicity in skull base chordoma and chondrosarcoma-a single-institution clinical experience

BACKGROUND

Chordomas and chondrosarcomas of the skull base are rare, slowly growing malignant bone neoplasms. Despite their radioresistant properties, proton therapy has been successfully used as an adjunct to resection or as a definitive treatment. Herewith, we present our experience with robustly optimized intensity-modulated proton therapy (IMPT) and related toxicities in skull base chordoma and chondrosarcoma patients treated at HollandPTC, Delft, the Netherlands.

METHODS

Clinical data, treatment plans, and acute toxicities of patients treated between July 2019 and August 2021 were reviewed. CT and 3.0T MRI scans for treatment planning were performed in supine position in a thermoplastic mold. In total, 21 dose optimization and 28 dose evaluation scenarios were simulated. Acute toxicity was scored weekly before and during the treatment according to the CTCAE v4.0. Median follow-up was 35 months (range 12-36 months).

RESULTS

Overall, 9 chordoma and 3 chondrosarcoma patients with 1-3 resections prior to IMPT were included; 4 patients had titanium implants. Brainstem core and surface and spinal cord core and surface were used for nominal plan robust optimization in 11, 10, 8, and 7 patients, respectively. Middle ear inflammation, dry mouth, radiation dermatitis, taste disorder, and/or alopecia of grades 1-3 were noted at the end of treatment among 6 patients without similar complaints at inclusion; symptoms disappeared 3 months following the treatment.

CONCLUSION

Robustly optimized IMPT is clinically feasible as a postoperative treatment for skull base chordoma and chondrosarcoma patients. We observed acceptable early toxicities (grade 1-3) that disappeared within the first 3 months after irradiation.

NRG Oncology and Particle Therapy Co-Operative Group Patterns of Practice Survey and Consensus Recommendations on Pencil-Beam Scanning Proton Stereotactic Body Radiation Therapy and Hypofractionated Radiation Therapy for Thoracic Malignancies

Stereotactic body radiation therapy (SBRT) and hypofractionation using pencil-beam scanning (PBS) proton therapy (PBSPT) is an attractive option for thoracic malignancies. Combining the advantages of target coverage conformity and critical organ sparing from both PBSPT and SBRT, this new delivery technique has great potential to improve the therapeutic ratio, particularly for tumors near critical organs. Safe and effective implementation of PBSPT SBRT/hypofractionation to treat thoracic malignancies is more challenging than the conventionally fractionated PBSPT because of concerns of amplified uncertainties at the larger dose per fraction. The NRG Oncology and Particle Therapy Cooperative Group Thoracic Subcommittee surveyed proton centers in the United States to identify practice patterns of thoracic PBSPT SBRT/hypofractionation. From these patterns, we present recommendations for future technical development of proton SBRT/hypofractionation for thoracic treatment. Among other points, the recommendations highlight the need for volumetric image guidance and multiple computed tomography-based robust optimization and robustness tools to minimize further the effect of uncertainties associated with respiratory motion. Advances in direct motion analysis techniques are urgently needed to supplement current motion management techniques.

A novel fast robust optimization algorithm for intensity-modulated proton therapy with minimum monitor unit constraint

BACKGROUND

Intensity-modulated proton therapy (IMPT) optimizes spot intensities and position, providing better conformability. However, the successful application of IMPT is dependent upon addressing the challenges posed by range and setup uncertainties. In order to address the uncertainties in IMPT, robust optimization is essential.

PURPOSE

This study aims to develop a novel fast algorithm for robust optimization of IMPT with minimum monitor unit (MU) constraint.

METHODS AND MATERIALS

The study formulates a robust optimization problem and proposes a novel, fast algorithm based on the alternating direction method of multipliers (ADMM) framework. This algorithm enables distributed computation and parallel processing. Ten clinical cases were used as test scenarios to evaluate the performance of the proposed approach. The robust optimization method (RBO-NEW) was compared with plans that only consider nominal optimization using CTV (NMO-CTV) without handling uncertainties and PTV (NMO-PTV) to handle the uncertainties, as well as with conventional robust-optimized plans (RBO-CONV). Dosimetric metrics, including D95, homogeneity index, and Dmean, were used to evaluate the dose distribution quality. The area under the root-mean-square dose (RMSD)-volume histogram curves (AUC) and dose-volume histogram (DVH) bands were used to evaluate the robustness of the treatment plan. Optimization time cost was also assessed to measure computational efficiency.

RESULTS

The results demonstrated that the RBO plans exhibited better plan quality and robustness than the NMO plans, with RBO-NEW showing superior computational efficiency and plan quality compared to RBO-CONV. Specifically, statistical analysis results indicated that RBO-NEW was able to reduce the computational time from389.70 ± 207.40 $389.70\pm 207.40$ to228.60 ± 123.67 $228.60\pm 123.67$ s (p < 0.01 $p<0.01$ ) and reduce the mean organ-at-risk (OAR) dose from9.38 ± 12.80 $9.38\pm 12.80$ % of the prescription dose to9.07 ± 12.39 $9.07\pm 12.39$ % of the prescription dose (p < 0.05 $p<0.05$ ) compared to RBO-CONV.

CONCLUSION

This study introduces a novel fast robust optimization algorithm for IMPT treatment planning with minimum MU constraint. Such an algorithm is not only able to enhance the plan's robustness and computational efficiency without compromising OAR sparing but also able to improve treatment plan quality and reliability.

Technical Note: Generalizable and Promptable Artificial Intelligence Model to Augment Clinical Delineation in Radiation Oncology

BACKGROUND

Efficient and accurate delineation of organs at risk (OARs) is a critical procedure for treatment planning and dose evaluation. Deep learning-based auto-segmentation of OARs has shown promising results and is increasingly being used in radiation therapy. However, existing deep learning-based auto-segmentation approaches face two challenges in clinical practice: generalizability and human-AI interaction. A generalizable and promptable auto-segmentation model, which segments OARs of multiple disease sites simultaneously and supports on-the-fly human-AI interaction, can significantly enhance the efficiency of radiation therapy treatment planning.

PURPOSE

Meta's segment anything model (SAM) was proposed as a generalizable and promptable model for next-generation natural image segmentation. We further evaluated the performance of SAM in radiotherapy segmentation.

METHODS

Computed tomography (CT) images of clinical cases from four disease sites at our institute were collected: prostate, lung, gastrointestinal, and head & neck. For each case, we selected the OARs important in radiotherapy treatment planning. We then compared both the Dice coefficients and Jaccard indices derived from three distinct methods: manual delineation (ground truth), automatic segmentation using SAM's 'segment anything' mode, and automatic segmentation using SAM's 'box prompt' mode that implements manual interaction via live prompts during segmentation.

RESULTS

Our results indicate that SAM's segment anything mode can achieve clinically acceptable segmentation results in most OARs with Dice scores higher than 0.7. SAM's box prompt mode further improves Dice scores by 0.1∼0.5. Similar results were observed for Jaccard indices. The results show that SAM performs better for prostate and lung, but worse for gastrointestinal and head & neck. When considering the size of organs and the distinctiveness of their boundaries, SAM shows better performance for large organs with distinct boundaries, such as lung and liver, and worse for smaller organs with less distinct boundaries, like parotid and cochlea.

CONCLUSIONS

Our results demonstrate SAM's robust generalizability with consistent accuracy in automatic segmentation for radiotherapy. Furthermore, the advanced box-prompt method enables the users to augment auto-segmentation interactively and dynamically, leading to patient-specific auto-segmentation in radiation therapy. SAM's generalizability across different disease sites and different modalities makes it feasible to develop a generic auto-segmentation model in radiotherapy.

Beam mask and sliding window-facilitated deep learning-based accurate and efficient dose prediction for pencil beam scanning proton therapy

BACKGROUND

Accurate and efficient dose calculation is essential for on-line adaptive planning in proton therapy. Deep learning (DL) has shown promising dose prediction results in photon therapy. However, there is a scarcity of DL-based dose prediction methods specifically designed for proton therapy. Successful dose prediction method for proton therapy should account for more challenging dose prediction problems in pencil beam scanning proton therapy (PBSPT) due to its sensitivity to heterogeneities.

PURPOSE

To develop a DL-based PBSPT dose prediction workflow with high accuracy and balanced complexity to support on-line adaptive proton therapy clinical decision and subsequent replanning.

METHODS

PBSPT plans of 103 prostate cancer patients (93 for training and the other 10 for independent testing) and 83 lung cancer patients (73 for training and the other 10 for independent testing) previously treated at our institution were included in the study, each with computed tomography scans (CTs), structure sets, and plan doses calculated by the in-house developed Monte-Carlo dose engine (considered as the ground truth in the model training and testing). For the ablation study, we designed three experiments corresponding to the following three methods: (1) Experiment 1, the conventional region of interest (ROI) (composed of targets and organs-at-risk [OARs]) method. (2) Experiment 2, the beam mask (generated by raytracing of proton beams) method to improve proton dose prediction. (3) Experiment 3, the sliding window method for the model to focus on local details to further improve proton dose prediction. A fully connected 3D-Unet was adopted as the backbone. Dose volume histogram (DVH) indices, 3D Gamma passing rates with a criterion of 3%/3 mm/10%, and dice coefficients for the structures enclosed by the iso-dose lines between the predicted and the ground truth doses were used as the evaluation metrics. The calculation time for each proton dose prediction was recorded to evaluate the method's efficiency.

RESULTS

Compared to the conventional ROI method, the beam mask method improved the agreement of DVH indices for both targets and OARs and the sliding window method further improved the agreement of the DVH indices (for lung cancer, CTV D98 absolute deviation: 0.74 ± 0.18 vs. 0.57 ± 0.21 vs. 0.54 ± 0.15 Gy[RBE], ROI vs. beam mask vs. sliding window methods, respectively). For the 3D Gamma passing rates in the target, OARs, and BODY (outside target and OARs), the beam mask method improved the passing rates in these regions and the sliding window method further improved them (for prostate cancer, targets: 96.93% ± 0.53% vs. 98.88% ± 0.49% vs. 99.97% ± 0.07%, BODY: 86.88% ± 0.74% vs. 93.21% ± 0.56% vs. 95.17% ± 0.59%). A similar trend was also observed for the dice coefficients. This trend was especially remarkable for relatively low prescription isodose lines (for lung cancer, 10% isodose line dice: 0.871 ± 0.027 vs. 0.911 ± 0.023 vs. 0.927 ± 0.017). The dose predictions for all the testing cases were completed within 0.25 s.

CONCLUSIONS

An accurate and efficient deep learning-augmented proton dose prediction framework has been developed for PBSPT, which can predict accurate dose distributions not only inside but also outside ROI efficiently. The framework can potentially further reduce the initial planning and adaptive replanning workload in PBSPT.

Proton Pencil-Beam Scanning Stereotactic Body Radiation Therapy and Hypofractionated Radiation Therapy for Thoracic Malignancies: Patterns of Practice Survey and Recommendations for Future Development from NRG Oncology and PTCOG

Stereotactic body radiation therapy (SBRT) and hypofractionation using pencil-beam scanning (PBS) proton therapy (PBSPT) is an attractive option for thoracic malignancies. Combining the advantages of target coverage conformity and critical organ sparing from both PBSPT and SBRT, this new delivery technique has great potential to improve the therapeutic ratio, particularly for tumors near critical organs. Safe and effective implementation of PBSPT SBRT/hypofractionation to treat thoracic malignancies is more challenging than the conventionally-fractionated PBSPT due to concerns of amplified uncertainties at the larger dose per fraction. NRG Oncology and Particle Therapy Cooperative Group (PTCOG) Thoracic Subcommittee surveyed US proton centers to identify practice patterns of thoracic PBSPT SBRT/hypofractionation. From these patterns, we present recommendations for future technical development of proton SBRT/hypofractionation for thoracic treatment. Amongst other points, the recommendations highlight the need for volumetric image guidance and multiple CT-based robust optimization and robustness tools to minimize further the impact of uncertainties associated with respiratory motion. Advances in direct motion analysis techniques are urgently needed to supplement current motion management techniques.

Modelling small block aperture in an in-house developed GPU-accelerated Monte Carlo-based dose engine for pencil beam scanning proton therapy

Purpose

To enhance an in-house graphic-processing-unit (GPU) accelerated virtual particle (VP)-based Monte Carlo (MC) proton dose engine (VPMC) to model aperture blocks in both dose calculation and optimization for pencil beam scanning proton therapy (PBSPT)-based stereotactic radiosurgery (SRS).

Methods and Materials

A module to simulate VPs passing through patient-specific aperture blocks was developed and integrated in VPMC based on simulation results of realistic particles (primary protons and their secondaries). To validate the aperture block module, VPMC was first validated by an opensource MC code, MCsquare, in eight water phantom simulations with 3cm thick brass apertures: four were with aperture openings of 1, 2, 3, and 4cm without a range shifter, while the other four were with same aperture opening configurations with a range shifter of 45mm water equivalent thickness. Then, VPMC was benchmarked with MCsquare and RayStation MC for 10 patients with small targets (average volume 8.4 cc with range of 0.4 - 43.3 cc). Finally, 3 typical patients were selected for robust optimization with aperture blocks using VPMC.

Results

In the water phantoms, 3D gamma passing rate (2%/2mm/10%) between VPMC and MCsquare was 99.71±0.23%. In the patient geometries, 3D gamma passing rates (3%/2mm/10%) between VPMC/MCsquare and RayStation MC were 97.79±2.21%/97.78±1.97%, respectively. Meanwhile, the calculation time was drastically decreased from 112.45±114.08 seconds (MCsquare) to 8.20±6.42 seconds (VPMC) with the same statistical uncertainties of ~0.5%. The robustly optimized plans met all the dose-volume-constraints (DVCs) for the targets and OARs per our institutional protocols. The mean calculation time for 13 influence matrices in robust optimization by VPMC was 41.6 seconds and the subsequent on-the-fly "trial-and-error" optimization procedure took only 71.4 seconds on average for the selected three patients.

Conclusion

VPMC has been successfully enhanced to model aperture blocks in dose calculation and optimization for the PBSPT-based SRS.

Beam mask and sliding window-facilitated deep learning-based accurate and efficient dose prediction for pencil beam scanning proton therapy

PURPOSE

To develop a DL-based PBSPT dose prediction workflow with high accuracy and balanced complexity to support on-line adaptive proton therapy clinical decision and subsequent replanning.

METHODS

PBSPT plans of 103 prostate cancer patients and 83 lung cancer patients previously treated at our institution were included in the study, each with CTs, structure sets, and plan doses calculated by the in-house developed Monte-Carlo dose engine. For the ablation study, we designed three experiments corresponding to the following three methods: 1) Experiment 1, the conventional region of interest (ROI) method. 2) Experiment 2, the beam mask (generated by raytracing of proton beams) method to improve proton dose prediction. 3) Experiment 3, the sliding window method for the model to focus on local details to further improve proton dose prediction. A fully connected 3D-Unet was adopted as the backbone. Dose volume histogram (DVH) indices, 3D Gamma passing rates, and dice coefficients for the structures enclosed by the iso-dose lines between the predicted and the ground truth doses were used as the evaluation metrics. The calculation time for each proton dose prediction was recorded to evaluate the method's efficiency.

RESULTS

Compared to the conventional ROI method, the beam mask method improved the agreement of DVH indices for both targets and OARs and the sliding window method further improved the agreement of the DVH indices. For the 3D Gamma passing rates in the target, OARs, and BODY (outside target and OARs), the beam mask method can improve the passing rates in these regions and the sliding window method further improved them. A similar trend was also observed for the dice coefficients. In fact, this trend was especially remarkable for relatively low prescription isodose lines. The dose predictions for all the testing cases were completed within 0.25s.

Technical note: Investigation of dose and LETd effect to rectum and bladder by using non-straight laterals in prostate cancer receiving proton therapy

BACKGROUND

Parallel-opposed lateral beams are the conventional beam arrangements in proton therapy for prostate cancer. However, when considering linear energy transfer (LET) and RBE effects, alternative beam arrangements should be investigated.

PURPOSE

To investigate the dose and dose averaged LET (LET_d ) impact of using new beam arrangements rotating beams 5°-15° posteriorly to the laterals in prostate cancer treated with pencil-beam-scanning (PBS) proton therapy.

METHODS

Twenty patients with localized prostate cancer were included in this study. Four proton treatment plans for each patient were generated utilizing 0°, 5°, 10°, and 15° posterior oblique beam pairs relative to parallel-opposed lateral beams. Dose-volume histograms (DVHs) from posterior oblique beams were analyzed. Dose-LET_d -volume histogram (DLVH) was employed to study the difference in dose and LET_d with each beam arrangement. DLVH indices, V ( d , l ) $V( {d,l} )$ , defined as the cumulative absolute volume that has a dose of at least d (Gy[RBE]) and a LET_d of at least l (keV/µm), were calculated for both the rectum and bladder to the whole group of patients and two-sub groups with and without hydrogel spacer. These metrics were tested using Wilcoxon signed-rank test.

RESULTS

Rotating beam angles from laterals to slightly posterior by 5°-15° reduced high LET_d volumes while it increased the dose volume in the rectum and increased LET_d in bladders. Beam angles rotated five degrees posteriorly from laterals (i.e., gantry in 95° and 265°) are proposed since they achieved the optimal balance of better LET_d sparing and minimal dose increase in the rectum. A reduction of V(50 Gy[RBE], 2.6 keV/µm) from 7.41 to 3.96 cc (p < 0.01), and a slight increase of V(50 Gy[RBE], 0 keV/µm) from 20.1 to 21.6 cc (p < 0.01) were observed for the group without hydrogel spacer. The LET_d sparing was less effective for the group with hydrogel spacer, which achieved the reduction of V(50 Gy[RBE], 2.6 keV/µm) from 4.28 to 2.10 cc (p < 0.01).

CONCLUSIONS

Posterior oblique angle plans improved LET_d sparing of the rectum while sacrificing LET_d sparing in the bladder in the treatment of prostate cancer with PBS. Beam angle modification from laterals to slightly posterior may be a strategy to redistribute LET_d and perhaps reduce rectal toxicity risks in prostate cancer patients treated with PBS. However, the effect is reduced for patients with hydrogel spacer.

Cardiopulmonary Toxicity Following Intensity-Modulated Proton Therapy (IMPT) Versus Intensity-Modulated Radiation Therapy (IMRT) for Stage III Non-Small Cell Lung Cancer

INTRODUCTION

Intensity-modulated proton therapy (IMPT) has the potential to reduce radiation dose to normal organs when compared to intensity-modulated radiation therapy (IMRT). We hypothesized that IMPT is associated with a reduced rate of cardiopulmonary toxicities in patients with Stage III NSCLC when compared with IMRT.

METHODS

We analyzed 163 consecutively treated patients with biopsy-proven, stage III NSCLC who received IMPT (n = 35, 21%) or IMRT (n = 128, 79%). Patient, tumor, and treatment characteristics were analyzed. Overall survival (OS), freedom-from distant metastasis (FFDM), freedom-from locoregional relapse (FFLR), and cardiopulmonary toxicities (CTCAE v5.0) were calculated using the Kaplan-Meier estimate. Univariate cox regressions were conducted for the final model.

RESULTS

Median follow-up of surviving patients was 25.5 (range, 4.6-58.1) months. Median RT dose was 60 (range, 45-72) Gy [RBE]. OS, FFDM, and FFLR were not different based on RT modality. IMPT provided significant dosimetric pulmonary and cardiac sparing when compared to IMRT. IMPT was associated with a reduced rate of grade more than or equal to 3 pneumonitis (HR 0.25, P = .04) and grade more than or equal to 3 cardiac events (HR 0.33, P = .08). Pre-treatment predicted diffusing capacity for carbon monoxide less than equal to 57% (HR 2.8, P = .04) and forced expiratory volume in the first second less than equal to 61% (HR 3.1, P = .03) were associated with an increased rate of grade more than or equal to 3 pneumonitis.

CONCLUSIONS

IMPT is associated with a reduced risk of clinically significant pneumonitis and cardiac events when compared with IMRT without compromising tumor control in stage III NSCLC. IMPT may provide a safer treatment option, particularly for high-risk patients with poor pretreatment pulmonary function.

Collimating individual beamlets in pencil beam scanning proton therapy, a dosimetric investigation

The purpose of this work is to investigate collimating individual proton beamlets from a dosimetric perspective and to introduce a new device concept, the spot scanning aperture (SSA). The SSA consists of a thin aperture with a small cylindrical opening attached to a robotics system, which allows the aperture to follow and align with individual beamlets during spot delivery. Additionally, a range shifter is incorporated (source-side) for treating shallow depths. Since the SSA trims beamlets spot by spot, the patient-facing portion of the device only needs to be large enough to trim a single proton beamlet. The SSA has been modelled in an open-source Monte-Carlo-based dose engine (MCsquare) to characterize its dosimetric properties in water at depths between 0 and 10 cm while varying the following parameters: the aperture material, thickness, distance to the water phantom, distance between the aperture and attached range shifter, and the aperture opening radius. Overall, the SSA greatly reduced spot sizes for all the aperture opening radii that were tested (1 - 4 mm), especially in comparison with the extended range shifter (ranger shifter placed at 30 cm from patient); greater than 50% when placed less than 10 cm away from the patient at depths in water less than 50 mm. The peak to entrance dose ratio and linear energy transfer was found to depend on the thickness of the aperture and therefore the aperture material. Neutron production rates were also investigated and discussed.

Exploratory study of seed spots analysis to characterize dose and linear energy transfer effect in adverse event initialization of pencil beam scanning proton therapy

BACKGROUND

Both dose and linear-energy-transfer (LET) could play a substantial role in adverse event (AE) initialization of cancer patients treated with pencil-beam-scanning proton therapy (PBS). However, not all the voxels within the AE regions are directly induced from the dose and LET effect. It is important to study the synergistic effect of dose and LET in AE initialization by only including a subset of voxels that are dosimetrically important.

PURPOSE

To perform exploratory investigation of the dose and LET effects upon AE initialization in PBS using seed spots analysis.

METHODS

113 head and neck (H&N) cancer patients receiving curative PBS were included. Among them, 20 patients experienced unanticipated CTCAEv4.0 grade≥3 AEs (AE group) and 93 patients did not (control group). Within the AE group, 13 AE patients were included in the seed spot analysis to derive the descriptive features of AE initialization and the remaining 7 mandible osteoradionecrosis patients and 93 control patients were used to derive the feature-based volume constraint of mandible osteoradionecrosis. The AE regions were contoured and the corresponding dose-LET volume histograms (DLVHs) of AE regions were generated for all patients in the AE group. We selected high LET voxels (the highest 5% of each dose bin) with a range of moderate to high dose (≥∼40 Gy[RBE]) as critical voxels. Critical voxels which were contiguous with each other were grouped into clusters. Each cluster was considered as a potential independent seed spot for AE initialization. Seed spots were displayed in a 2D dose-LET plane based on their mean dose and LET to derive the descriptive features of AE initialization. A volume constraint of mandible osteoradionecrosis was then established based on the extracted features using a receiver operating characteristic curve.

RESULTS

The product of dose and LET (xBD) was found to be a descriptive feature of seed spots leading to AE initialization in this preliminary study. The derived xBD volume constraint for mandible osteoradionecrosis showed good performance with an area-under-curve of 0.87 (sensitivity of 0.714 and specificity of 0.807 in the leave-one-out cross validation) for the very limited patient data included in this study.

CONCLUSION

Our exploratory study showed that both dose and LET were observed to be important in AE initializations. The derived xBD volume constraint could predict mandible osteoradionecrosis reasonably well in the very limited H&N cancer patient data treated with PBS included in this study. This article is protected by copyright. All rights reserved.

GPU-accelerated Monte Carlo-based online adaptive proton therapy - a feasibility study

PURPOSE

To develop an online Graphic-Processing-Unit (GPU)-accelerated Monte-Carlo-based adaptive radiation therapy (ART) workflow for pencil beam scanning (PBS) proton therapy to address inter-fraction anatomical changes in patients treated with PBS.

METHODS AND MATERIALS

A four-step workflow was developed using our in-house developed GPU-accelerated Monte-Carlo-based treatment planning system to implement online Monte-Carlo-based ART for PBS. The first step conducts diffeomorphic demon-based deformable image registration (DIR) to propagate contours on the initial planning CT (pCT) to the verification CT (vCT) to form a new structure set. The second step performs forward dose calculation of the initial plan on the vCT with the propagated contours after manual approval (possible modifications involved). The third step triggers a re-optimization of the plan depending on whether the verification dose meets the clinical requirements or not. A robust evaluation will be done for both the verification plan in the second step and the re-opotimized plan in the third step. The fourth step involves a two-stage (before and after delivery) patient specific quality assurance (PSQA) of the re-optimized plan. The before-delivery PSQA is to compare the plan dose to the dose calculated using an independent fast open-source Monte Carlo code, MCsquare. The after-delivery PSQA is to compare the plan dose to the dose re-calculated using the log file (spot MU, spot position, and spot energy) collected during the delivery. Jaccard index (JI), Dice similarity coefficients (DSCs), and Hausdorff distance (HD) were used to assess the quality of the propagated contours in the first step. A commercial plan evaluation software, ClearCheck™, was integrated into the workflow to carry out efficient plan evaluation. 3D Gamma analysis was used during the fourth step to ensure the accuracy of the plan dose from re-optimization. Three patients with three different disease sites were chosen to evaluate the feasibility of the online ART workflow for PBS.

RESULTS

For all three patients, the propagated contours were found to have good volume conformance [JI (lowest-highest: 0.833-0.983) and DSC (0.909-0.992)] but sub-optimal boundary coincidence [HD (2.37-20.76 mm)] for organs at risk (OARs). The verification dose evaluated by ClearCheck™ showed significant degradation of the target coverage due to the inter-fractional anatomical changes. Re-optimization on the vCT resulted in great improvement of the plan quality to a clinically acceptable level. 3D Gamma analyses of PSQA confirmed the accuracy of the plan dose before delivery (mean Gamma index = 98.74% with a threshold of 2%/2 mm/10%), and after delivery based on the log files (mean Gamma index = 99.05% with a threshold of 2%/2 mm/10%). The average time cost for the complete execution of the workflow was around 858 seconds, excluding the time for manual intervention.

CONCLUSION

The proposed online ART workflow for PBS was demonstrated to be efficient and effective by generating a re-optimized plan that significantly improved the plan quality. This article is protected by copyright. All rights reserved.

Empirical Relative Biological Effectiveness (RBE) for Mandible Osteoradionecrosis (ORN) in Head and Neck Cancer Patients Treated With Pencil-Beam-Scanning Proton Therapy (PBSPT): A Retrospective, Case-Matched Cohort Study

Purpose

To retrospectively investigate empirical relative biological effectiveness (RBE) for mandible osteoradionecrosis (ORN) in head and neck (H&N) cancer patients treated with pencil-beam-scanning proton therapy (PBSPT).

Methods

We included 1,266 H&N cancer patients, of which, 931 patients were treated with volumetric-modulated arc therapy (VMAT) and 335 were treated with PBSPT. Among them, 26 VMAT and 9 PBSPT patients experienced mandible ORN (ORN group), while all others were included in the control group. To minimize the impact of the possible imbalance in clinical factors between VMAT and PBSPT patients in the dosimetric comparison between these two modalities and the resulting RBE quantification, we formed a 1:1 case-matched patient cohort (335 VMAT patients and 335 PBSPT patients including both the ORN and control groups) using the greedy nearest neighbor matching of propensity scores. Mandible dosimetric metrics were extracted from the case-matched patient cohort and statistically tested to evaluate the association with mandibular ORN to derive dose volume constraints (DVCs) for VMAT and PBSPT, respectively. We sought the equivalent constraint doses for VMAT so that the critical volumes of VMAT were equal to those of PBSPT at different physical doses. Empirical RBEs of PBSPT for ORN were obtained by calculating the ratio between the derived equivalent constraint doses and physical doses of PBSPT. Bootstrapping was further used to get the confidence intervals.

Results

Clinical variables of age, gender, tumor stage, prescription dose, chemotherapy, hypertension or diabetes, dental extraction, smoking history, or current smoker were not statistically related to the incidence of ORN in the overall patient cohort. Smoking history was found to be significantly associated with the ORN incidence in PBSPT patients only. V40Gy[RBE], V50Gy[RBE], and V60Gy[RBE] were statistically different (p<0.05) between the ORN and control group for VMAT and PBSPT. Empirical RBEs of 1.58(95%CI: 1.34-1.64), 1.34(95%CI: 1.23-1.40), and 1.24(95%: 1.15-1.26) were obtained for proton dose at 40 Gy[RBE=1.1], 50 Gy[RBE=1.1] and 60 Gy[RBE=1.1], respectively.

Conclusions

Our study suggested that RBEs were larger than 1.1 at moderate doses (between 40 and 60 Gy[RBE=1.1]) with high LET for mandible ORN. RBEs are underestimated in current clinical practice in PBSPT. The derived DVCs can be used for PBSPT plan evaluation and optimization to minimize the incidence rate of mandible ORN.

Proton Therapy for Prostate Cancer: Challenges and Opportunities

Simple Summary

Reported clinical outcomes of proton therapy (PT) for localized prostate cancer are similar to photon-based external beam radiotherapy. Apparently, the dosimetric advantages of PT have yet to be translated to clinical benefits. The suboptimal clinical outcomes of PT might be attributable to inadequate dose prescription, as indicated by the ASCENDE-RT trial. Moreover, uncertainties involved in the treatment planning and delivery processes, as well as technological limitations in PT treatment systems, may lead to discrepancies between planned doses and actual doses delivered to patients. In this article, we reviewed the current status of PT for prostate cancer and discussed different clinical implementations that could potentially improve the clinical outcome of PT for prostate cancer. Various technological advancements under which uncertainties in dose calculations can be minimized, including MRI-guided PT, dual-energy photon-counting CT and high-resolution Monte Carlo-based treatment planning systems, are highlighted.

Abstract

The dosimetric advantages of proton therapy (PT) treatment plans are demonstrably superior to photon-based external beam radiotherapy (EBRT) for localized prostate cancer, but the reported clinical outcomes are similar. This may be due to inadequate dose prescription, especially in high-risk disease, as indicated by the ASCENDE-RT trial. Alternatively, the lack of clinical benefits with PT may be attributable to improper dose delivery, mainly due to geometric and dosimetric uncertainties during treatment planning, as well as delivery procedures that compromise the dose conformity of treatments. Advanced high-precision PT technologies, and treatment planning and beam delivery techniques are being developed to address these uncertainties. For instance, external magnetic resonance imaging (MRI)-guided patient setup rooms are being developed to improve the accuracy of patient positioning for treatment. In-room MRI-guided patient positioning systems are also being investigated to improve the geometric accuracy of PT. Soon, high-dose rate beam delivery systems will shorten beam delivery time to within one breath hold, minimizing the effects of organ motion and patient movements. Dual-energy photon-counting computed tomography and high-resolution Monte Carlo-based treatment planning systems are available to minimize uncertainties in dose planning calculations. Advanced in-room treatment verification tools such as prompt gamma detector systems will be used to verify the depth of PT. Clinical implementation of these new technologies is expected to improve the accuracy and dose conformity of PT in the treatment of localized prostate cancers, and lead to better clinical outcomes. Improvement in dose conformity may also facilitate dose escalation, improving local control and implementation of hypofractionation treatment schemes to improve patient throughput and make PT more cost effective.

Per-voxel constraints to minimize hot spots in linear energy transfer-guided robust optimization for base of skull head and neck cancer patients in IMPT

PURPOSE

Due to the employment of quadratic programming using soft constraints to implement dose volume constraints and the "trial-and-error" procedure needed to achieve a clinically acceptable plan, conventional dose volume constraints (upper limit) are not adequately effective in controlling small and isolated hot spots in the dose/linear energy transfer (LET) distribution. Such hot spots can lead to adverse events. In order to mitigate the risk of brain necrosis, one of the most clinically significant adverse events in patients receiving intensity-modulated proton therapy (IMPT) for base of skull (BOS) cancer, we propose per-voxel constraints to minimize hot spots in LET-guided robust optimization.

METHODS AND MATERIALS

Ten BOS cancer patients treated with IMPT were carefully selected by meeting one of the following conditions: (1) diagnosis of brain necrosis during follow-up; and (2) considered high risk for brain necrosis by not meeting dose constraints to the brain. An optimizing structure (BrainOPT) and an evaluating structure (BrainROI) that both contained the aforementioned hot dose regions in the brain were generated for optimization and evaluation, respectively. Two plans were generated for every patient: one using conventional dose-only robust optimization, the other using LET-guided robust optimization. The impact of LET was integrated into the optimization via a term of extra biological dose (xBD). A novel optimization tool of per-voxel constraints to control small and isolated hot spots in either the dose, LET, or combined (dose/LET) distribution was developed and used to minimize dose/LET hot spots of the selected structures. Indices from dose-volume histogram (DVH) and xBD dose-volume histogram (xBDVH) were used in the plan evaluation. A newly developed tool of the dose-LET-volume histogram (DLVH) was also adopted to illustrate the underlying mechanism. Wilcoxon signed-rank test was used for statistical comparison of the DVH and xBDVH indices between the conventional dose-only and the LET-guided robustly optimized plans.

RESULTS

Per-voxel constraints effectively and efficiently minimized dose hot spots in both dose-only and LET-guided robust optimization and LET hot spots in LET-guided robust optimization. Compared to the conventional dose-only robust optimization, the LET-guided robust optimization could generate plans with statistically lower xBD hot spots in BrainROI (VxBD,50 Gy[RBE], p = 0.009; VxBD,60 Gy[RBE], p = 0.025; xBD1cc, p = 0.017; xBD2cc, p = 0.022) with comparable dose coverage, dose hot spots in the target, and dose hot spots in BrainROI. DLVH analysis indicated that LET-guided robust optimization could either reduce LET at the same dose level or redistribute high LET from high dose regions to low dose regions.

CONCLUSION

Per-voxel constraint is a powerful tool to minimize dose/LET hot spots in IMPT. The LET-guided robustly optimized plans outperformed the conventional dose-only robustly optimized plans in terms of xBD hot spots control.

Single-institution clinical experience using robust intensity modulated proton therapy in chordoma and chondrosarcoma of the mobile spine and sacrum: Feasibility and need for plan adaptation

BACKGROUND

Due to its specific physical characteristics, proton irradiation is especially suited for irradiation of chordomas and chondrosarcoma in the axial skeleton. Robust plan optimization renders the proton beam therapy more predictable upon individual setup errors. Reported experience with the planning and delivery of robustly optimized plans in chordoma and chondrosarcoma of the mobile spine and sacrum, is limited. In this study, we report on the clinical use of robustly optimized, intensity modulated proton beam therapy in these patients.

METHODS

We retrospectively reviewed patient, treatment and acute toxicity data of all patients with chordoma and chondrosarcoma of the mobile spine and sacrum, treated between 1 April 2019 and 1 April 2020 at our institute. Anatomy changes during treatment were evaluated by weekly cone-beam CTs (CBCT), supplemented by scheduled control-CTs or ad-hoc control-CTs. Acute toxicity was scored weekly during treatment and at 3 months after therapy according to CTCAE 4.0.

RESULTS

17 chordoma and 3 chondrosarcoma patients were included. Coverage of the high dose clinical target volume was 99.8% (range 56.1-100%) in the nominal and 80.9% (range 14.3-99.6%) in the voxel-wise minimum dose distribution. Treatment plan adaptation was needed in 5 out of 22 (22.7%) plans. Reasons for plan adaptation were either reduced tumor coverage or increased dose to the OAR.

CONCLUSIONS

Robustly optimized intensity modulated proton beam therapy for chordoma and chondrosarcoma of the mobile spine is feasible. Plan adaptations due to anatomical changes were required in approximately 23 percent of treatment courses.

Adaptive proton therapy

Radiation therapy treatments are typically planned based on a single image set, assuming that the patient's anatomy and its position relative to the delivery system remains constant during the course of treatment. Similarly, the prescription dose assumes constant biological dose-response over the treatment course. However, variations can and do occur on multiple time scales. For treatment sites with significant intra-fractional motion, geometric changes happen over seconds or minutes, while biological considerations change over days or weeks. At an intermediate timescale, geometric changes occur between daily treatment fractions. Adaptive radiation therapy is applied to consider changes in patient anatomy during the course of fractionated treatment delivery. While traditionally adaptation has been done off-line with replanning based on new CT images, online treatment adaptation based on on-board imaging has gained momentum in recent years due to advanced imaging techniques combined with treatment delivery systems. Adaptation is particularly important in proton therapy where small changes in patient anatomy can lead to significant dose perturbations due to the dose conformality and finite range of proton beams. This review summarizes the current state-of-the-art of on-line adaptive proton therapy and identifies areas requiring further research.

Investigate the Dosimetric and Potential Clinical Benefits Utilizing Stereotactic Body Radiation Therapy With Simultaneous Integrated Boost Technique for Locally Advanced Pancreatic Cancer: A Comparison Between Photon and Proton Beam Therapy

Purpose

To investigate the potential clinical benefits of using stereotactic body radiation therapy (SBRT) with simultaneous integrated boost (SIB) technique for locally advanced pancreatic cancer (LAPC) among different treatment modalities and planning strategies, including photon and proton.

Method

A total of 19 patients were retrospectively selected in this study: 13 cases with the tumor located in the head of the pancreas and 6 cases with the tumor in the body of the pancreas. SBRT-SIB plans were generated using volumetric modulated arc therapy (VMAT), two-field Intensity Modulated Proton Therapy (IMPT), and three-field IMPT. The IMPT used the robust optimization parameters of ± 3.5% range and 5-mm setup uncertainties. Root-mean-square deviation dose (RMSD) volume histograms were used to evaluate the target coverage robustness quantitatively. Dosimetric metrics based on the dose-volume histogram (DVH), homogeneity index (HI), and normal tissue complication probability (NTCP) were analyzed to evaluate the potential clinical benefits among different planning groups.

Results

With a similar CTV and SIB coverage, two-field IMPT provided a lower maximum dose for the stomach (median: 18.6GyE, p<0.05) and duodenum (median: 32.62GyE, p<0.05) when the target was located in the head of the pancreas compared to VMAT and three-field IMPT. The risks of gastric bleed (3.42%) and grade ≥ 3 GI toxicity (4.55%) were also decreased. However, for the target in the body of the pancreas, VMAT showed a lower maximum dose for the stomach (median 30.93GyE, p<0.05) and toxicity of gastric bleed (median: 8.67%, p<0.05) compared to two-field IMPT and three-field IMPT, while other maximum doses and NTCPs were similar. The RMSD volume histogram (RVH) analysis shows that three-field IMPT provided better robustness for targets but not for OARs. Instead, three-field IMPT increased the Dmean of organs such as the stomach, duodenum, and intestine.

Conclusion

The results indicated that the tumor locations could play a critical role in determining clinical benefits among different treatment modalities. Two-field IMPT could be a better option for LAPC patients whose tumors are located in the head of the pancreas. It provides lower severe toxicity for the stomach and duodenum. Nevertheless, VMAT is preferred for the body with better protection for the possibility of gastric bleed.

A Critical Review of LET-Based Intensity-Modulated Proton Therapy Plan Evaluation and Optimization for Head and Neck Cancer Management

In this review article, we review the 3 important aspects of linear-energy-transfer (LET) in intensity-modulated proton therapy (IMPT) for head and neck (H&N) cancer management. Accurate LET calculation methods are essential for LET-guided plan evaluation and optimization, which can be calculated either by analytical methods or by Monte Carlo (MC) simulations. Recently, some new 3D analytical approaches to calculate LET accurately and efficiently have been proposed. On the other hand, several fast MC codes have also been developed to speed up the MC simulation by simplifying nonessential physics models and/or using the graphics processor unit (GPU)–acceleration approach. Some concepts related to LET are also briefly summarized including (1) dose-weighted versus fluence-weighted LET; (2) restricted versus unrestricted LET; and (3) microdosimetry versus macrodosimetry. LET-guided plan evaluation has been clinically done in some proton centers. Recently, more and more studies using patient outcomes as the biological endpoint have shown a positive correlation between high LET and adverse events sites, indicating the importance of LET-guided plan evaluation in proton clinics. Various LET-guided plan optimization methods have been proposed to generate proton plans to achieve biologically optimized IMPT plans. Different optimization frameworks were used, including 2-step optimization, 1-step optimization, and worst-case robust optimization. They either indirectly or directly optimize the LET distribution in patients while trying to maintain the same dose distribution and plan robustness. It is important to consider the impact of uncertainties in LET-guided optimization (ie, LET-guided robust optimization) in IMPT, since IMPT is sensitive to uncertainties including both the dose and LET distributions. We believe that the advancement of the LET-guided plan evaluation and optimization will help us exploit the unique biological characteristics of proton beams to improve the therapeutic ratio of IMPT to treat H&N and other cancers.

Technical Note: 4D robust optimization in small spot intensity-modulated proton therapy (IMPT) for distal esophageal carcinoma

PURPOSE

To compare the dosimetric performances of small-spot three-dimensional (3D) and four-dimensional (4D) robustly optimized intensity-modulated proton (IMPT) plans in the presence of uncertainties and interplay effect simultaneously for distal esophageal carcinoma.

METHOD AND MATERIALS

Thirteen (13) patients were selected and re-planned with small-spot ( σ  ~ 2-6 mm) 3D and 4D robust optimization in IMPT, respectively. The internal clinical target volumes (CTV_high3d , CTV_low3d ) were used in 3D robust optimization. Different CTVs (CTV_high4d , CTV_low4d ) were generated by subtracting an inner margin of the motion amplitudes in three cardinal directions from the internal CTVs and used in 4D robust optimization. All patients were prescribed the same dose to CTVs (50 Gy[RBE] for CTV_high3d /CTV_high4d and 45 Gy[RBE] for CTV_low3d /CTV_low4d ). Dose-volume histogram (DVH) indices were calculated to assess plan quality. Comprehensive plan robustness evaluations that consisted of 300 perturbed scenarios (10 different motion patterns to consider irregular motion (sampled from a Gaussian distribution) and 30 different uncertainties scenarios (sampled from a 4D uniform distribution) combined), were performed to quantify robustness to uncertainties and interplay effect simultaneously. Wilcoxon signed-rank test was used for statistical analysis.

RESULTS

Compared to 3D robustly optimized plans, 4D robustly optimized plans had statistically improved target coverage and better sparing of lungs and heart (heart D_mean , P = 0.001; heart V_30Gy[RBE] , P = 0.001) in the nominal scenario. 4D robustly optimized plans had better robustness in target dose coverage (CTV_high4d V_100% , P = 0.002) and the protection of lungs and heart (heart D_mean , P = 0.001; heart V_30Gy[RBE] , P = 0.001) when uncertainties and interplay effect were considered simultaneously.

CONCLUSIONS

Even with small spots in IMPT, 4D robust optimization outperformed 3D robust optimization in terms of normal tissue protection and robustness to uncertainties and interplay effect simultaneously. Our findings support the use of 4D robust optimization to treat distal esophageal carcinoma with small spots in IMPT.

Prospective Assessment of Early Proton Therapy-Induced Optic Neuropathy in Patients With Intracranial, Orbital or Sinonasal Tumors: Impact of A Standardized Ophthalmological Follow Up

PURPOSE

Proton therapy (PT) can be a good option to achieve tumor control while reducing the probability of radiation induced toxicities compared to X-ray-based radiotherapy. However, there are still uncertainties about the effects of PT on the organs in direct contact with the irradiated volume. The aim of this prospective series was to report 6-month follow-up of clinical and functional optic neuropathy rates of patients treated by proton therapy using a standardized comprehensive optic examination.

METHODS AND MATERIALS

Standardized ophthalmological examinations were performed to analyze subclinical anomalies in a systematic way before treatment and 6 months after the end of proton therapy with: Automatic visual field, Visual evoked potential (VEP) and optic coherence of tomography (OCT).

RESULTS

From October 2018 to July 2020 we analyzed 81 eyes. No significant differences were found in the analysis of the clinical examination of visual functions by the radiation oncologist. However, considering VEP, the impairment was statistically significant for both fibers explored at 30'angle (p:0.007) and 60'angle (p <0.001). In patients with toxicity, the distance of the target volume from the optical pathways was more important with a p-value for 30'VEP at 0.035 and for 60'VEP at 0.039.

CONCLUSIONS

These results confirm uncertainties concerning relative biological effectiveness of proton therapy, linear energy transfer appears to be more inhomogeneous especially in areas close to the target volumes. The follow-up of patients after proton therapy is not an easy process to set up but it is necessary to improve our knowledges about the biological effects of proton therapy in real life. Our study which will continue during the coming years, suggests that follow-up with in-depth examinations such as VEP as a biomarker could improve the detection of early abnormalities.

Robust treatment planning with 4D intensity modulated carbon ion therapy for multiple targets in stage IV non-small cell lung cancer

Intensity modulated particle therapy (IMPT) with carbon ions can generate highly conformal treatment plans; however, IMPT is limited in robustness against range and positioning uncertainty. This is particularly true for moving targets, even though all motion states of a 4DCT are considered in 4D-IMPT. Here, we expand 4D-IMPT to include robust non-linear RBE-weighted optimization to explore its potential in improving plan robustness and sparing critical organs. In this study, robust 4D-optimization-based on worst-case optimization on 9 scenarios-was compared to conventional 4D-optimization with PTV margins using 4D dose calculation and robustness analysis for 21 uncertainty scenarios. Slice-by-slice rescanning was used for motion mitigation. Both 4D-optimization strategies were tested on a cohort of 8 multi-lesion lung cancer patients with the goal of prioritizing OAR sparing in a hypofractionated treatment plan. Planning objectives were to keep the OAR volume doses below corresponding limits while simultaneously achieve CTV coverage with D95% ≥ 95 %. For the conventional plans, average D95% was at 98.7% which fulfilled the target objective in 83.2% of scenarios. For the robust plans, average D95% was reduced to 97.6% which still fulfilled the target objective in 80.7% of cases, but led to significantly improved overall OAR sparing: Volume doses were below the limits in 96.2% of cases for the conventional and 99.5% for the robust plans. When considering the particularly critical smaller airways only, fulfillment rates could be increased from 76.2% to 96% for the robust plans. This study has shown that plan robustness of 4D-IMPT could be improved by using robust 4D-optimization, offering greater control over uncertainties in the actual delivered dose. In some cases, this required sacrificing target coverage for the benefit of better OAR sparing.

Beam angle comparison for distal esophageal carcinoma patients treated with intensity-modulated proton therapy

Purpose

To compare the dosimetric performances of intensity‐modulated proton therapy (IMPT) plans generated with two different beam angle configurations (the Right–Left oblique posterior beams and the Superior–Inferior oblique posterior beams) for the treatment of distal esophageal carcinoma in the presence of uncertainties and interplay effect.

Methods and Materials

Twenty patients’ IMPT plans were retrospectively selected, with 10 patients treated with the R‐L oblique posterior beams (Group R‐L) and the other 10 patients treated with the S‐I oblique posterior beams (Group S‐I). Patients in both groups were matched by their clinical target volumes (CTVs—high and low dose levels) and respiratory motion amplitudes. Dose‐volume‐histogram (DVH) indices were used to assess plan quality. DVH bandwidth was calculated to evaluate plan robustness. Interplay effect was quantified using four‐dimensional (4D) dynamic dose calculation with random respiratory starting phase of each fraction. Normal tissue complication probability (NTCP) for heart, liver, and lung was calculated, respectively, to estimate the clinical outcomes. Wilcoxon signed‐rank test was used for statistical comparison between the two groups.

Results

Compared with plans in Group R‐L, plans in Group S‐I resulted in significantly lower liver D_mean and lung V_30Gy[RBE] with slightly higher but clinically acceptable spinal cord D_max. Similar plan robustness was observed between the two groups. When interplay effect was considered, plans in Group S‐I performed statistically better for heart D_mean and V_30Gy[RBE], lung Dmean and V_5Gy[RBE], and liver D_mean, with slightly increased but clinically acceptable spinal cord D_max. NTCP for liver was significantly better in Group S‐I.

Conclusions

IMPT plans in Group S‐I have better sparing of liver, heart, and lungs at the slight cost of spinal cord maximum dose protection, and are more interplay‐effect resilient compared to IMPT plans in Group R‐L. Our study supports the routine use of the S‐I oblique posterior beams for the treatments of distal esophageal carcinoma.

Intensity-modulated proton therapy (IMPT) interplay effect evaluation of asymmetric breathing with simultaneous uncertainty considerations in patients with non-small cell lung cancer

PURPOSE

Intensity-modulated proton therapy (IMPT) is sensitive to uncertainties from patient setup and proton beam range, as well as interplay effect. In addition, respiratory motion may vary from cycle to cycle, and also from day to day. These uncertainties can severely degrade the original plan quality and potentially affect patient's outcome. In this work, we developed a new tool to comprehensively consider the impact of all these uncertainties and provide plan robustness evaluation under them.

METHODS

We developed a comprehensive plan robustness evaluation tool that considered both uncertainties from patient setup and proton beam range, as well as respiratory motion simultaneously. To mimic patients' respiratory motion, the time spent in each phase was randomly sampled based on patient-specific breathing pattern parameters as acquired during the four-dimensional (4D)-computed tomography (CT) simulation. Spots were then assigned to one specific phase according to the temporal relationship between spot delivery sequence and patients' respiratory motion. Dose in each phase was calculated by summing contributions from all the spots delivered in that phase. The final 4D dynamic dose was obtained by deforming all doses in each phase to the maximum exhalation phase. Three hundred (300) scenarios (10 different breathing patterns with 30 different setup and range uncertainty scenario combinations) were calculated for each plan. The dose-volume histograms (DVHs) band method was used to assess plan robustness. Benchmarking the tool as an application's example, we compared plan robustness under both three-dimensional (3D) and 4D robustly optimized IMPT plans for 10 nonrandomly selected patients with non-small cell lung cancer.

RESULTS

The developed comprehensive plan robustness tool had been successfully applied to compare the plan robustness between 3D and 4D robustly optimized IMPT plans for 10 lung cancer patients. In the presence of interplay effect with uncertainties considered simultaneously, 4D robustly optimized plans provided significantly better CTV coverage (D_95% , P = 0.002), CTV homogeneity (D_5% -D_95% , P = 0.002) with less target hot spots (D_5% , P = 0.002), and target coverage robustness (CTV D_95% bandwidth, P = 0.004) compared to 3D robustly optimized plans. Superior dose sparing of normal lung (lung D_mean , P = 0.020) favoring 4D plans and comparable normal tissue sparing including esophagus, heart, and spinal cord for both 3D and 4D plans were observed. The calculation time for all patients included in this study was 11.4 ± 2.6 min.

CONCLUSION

A comprehensive plan robustness evaluation tool was successfully developed and benchmarked for plan robustness evaluation in the presence of interplay effect, setup and range uncertainties. The very high efficiency of this tool marks its clinical adaptation, highly practical and versatile nature, including possible real-time intra-fractional interplay effect evaluation as a potential application for future use.

The Role of Plan Robustness Evaluation in Comparing Protons and Photons Plans - An Application on IMPT and IMRT Plans in Skull Base Chordomas

PURPOSE

To analyze robustness of treatment plans optimized using different approaches in intensity modulated proton therapy (IMPT) and investigate the necessity of robust optimization and evaluation in intensity modulated radiotherapy (IMRT) plans for skull base chordomas.

MATERIALS AND METHODS

Two photon plans, standard IMRT and robustly optimized IMRT (RB-IMRT), and two IMPT plans, robustly optimized multi field optimization (MFO) and hybrid-MFO (HB-MFO), were created in RayStation TPS for five patients previously treated using single field uniform optimization (SFO). Both set-up and range uncertainties were incorporated during robust optimization of IMPT plans whereas only set-up uncertainty was used in RB-IMRT. The dosimetric outcomes from the five planning techniques were compared for every patient using standard dose volume indices and integral dose (ID) estimated for target and organs at risk (OARs). Robustness of each treatment plan was assessed by introducing set-up uncertainties of ±3 mm along the three translational axes and, only in protons, an additional range uncertainty of ±3.5%.

RESULTS

All the five nominal plans provided comparable and clinically acceptable target coverage. In comparison to nominal plans, worst case decrease in D95% of clinical target volume-high risk (CTV-HR) were 11.1%, 13.5%, and 13.6% for SFO, MFO, and HB-MFO plans respectively. The corresponding values were 13.7% for standard IMRT which improved to 11.5% for RB-IMRT. The worst case increased in high dose (D1%) to CTV-HR was highest in IMRT (2.1%) and lowest in SFO (0.7%) plans. Moreover, IMRT showed worst case increases in D1% for all neurological OARs and were lowest for SFO plans. The worst case D1% for brainstem, chiasm, spinal cord, optic nerves, and temporal lobes were increased by 29%, 41%, 30%, 41% and 14% for IMRT and 18%, 21%, 21%, 24%, and 7% for SFO plans, respectively. In comparison to IMRT, RB-IMRT improved D1% of all neurological OARs ranging from 5% to 14% in worst case scenarios.

CONCLUSION

Based on the five cases presented in the current study, all proton planning techniques (SFO, MFO and HB-MFO) were robust both for target coverage and OARs sparing. Standard IMRT plans were less robust than proton plans in regards to high doses to neurological OARs. However, robust optimization applied to IMRT resulted in improved robustness in both target coverage and high doses to OARs. Robustness evaluation may be considered as a part of plan evaluation procedure even in IMRT.

Efficacy and toxicity of particle radiotherapy in WHO grade II and grade III meningiomas: a systematic review

OBJECTIVEAdjuvant radiotherapy has become a common addition to the management of high-grade meningiomas, as immediate treatment with radiation following resection has been associated with significantly improved outcomes. Recent investigations into particle therapy have expanded into the management of high-risk meningiomas. Here, the authors systematically review studies on the efficacy and utility of particle-based radiotherapy in the management of high-grade meningioma.METHODSA literature search was developed by first defining the population, intervention, comparison, outcomes, and study design (PICOS). A search strategy was designed for each of three electronic databases: PubMed, Embase, and Scopus. Data extraction was conducted in accordance with the PRISMA guidelines. Outcomes of interest included local disease control, overall survival, and toxicity, which were compared with historical data on photon-based therapies.RESULTSEleven retrospective studies including 240 patients with atypical (WHO grade II) and anaplastic (WHO grade III) meningioma undergoing particle radiation therapy were identified. Five of the 11 studies included in this systematic review focused specifically on WHO grade II and III meningiomas; the others also included WHO grade I meningioma. Across all of the studies, the median follow-up ranged from 6 to 145 months. Local control rates for high-grade meningiomas ranged from 46.7% to 86% by the last follow-up or at 5 years. Overall survival rates ranged from 0% to 100% with better prognoses for atypical than for malignant meningiomas. Radiation necrosis was the most common adverse effect of treatment, occurring in 3.9% of specified cases.CONCLUSIONSDespite the lack of randomized prospective trials, this review of existing retrospective studies suggests that particle therapy, whether an adjuvant or a stand-alone treatment, confers survival benefit with a relatively low risk for severe treatment-derived toxicity compared to standard photon-based therapy. However, additional controlled studies are needed.

A fast robust optimizer for intensity modulated proton therapy using GPU

Robust optimization has been shown to be effective for stabilizing treatment planning in intensity modulated proton therapy (IMPT), but existing algorithms for the optimization process is time‐consuming. This paper describes a fast robust optimization tool that takes advantage of the GPU parallel computing technologies. The new robust optimization model is based on nine boundary dose distributions — two for ±range uncertainties, six for ±set‐up uncertainties along anteroposterior (A‐P), lateral (R‐L) and superior‐inferior (S‐I) directions, and one for nominal situation. The nine boundary influence matrices were calculated using an in‐house finite size pencil beam dose engine, while the conjugate gradient method was applied to minimize the objective function. The proton dose calculation algorithm and the conjugate gradient method were tuned for heterogeneous platforms involving the CPU host and GPU device. Three clinical cases — one head and neck cancer case, one lung cancer case, and one prostate cancer case — were investigated to demonstrate the clinical feasibility of the proposed robust optimizer. Compared with results from Varian Eclipse (version 13.3), the proposed method is found to be conducive to robust treatment planning that is less sensitive to range and setup uncertainties. The three tested cases show that targets can achieve high dose uniformity while organs at risks (OARs) are in better protection against setup and range errors. Based on the CPU + GPU heterogeneous platform, the execution times of the head and neck cancer case and the prostate cancer case are much less than half of Eclipse, while the run time of the lung cancer case is similar to that of Eclipse. The fast robust optimizer developed in this study can improve the reliability of traditional proton treatment planning in a much faster speed, thus making it possible for clinical utility.

Managing treatment-related uncertainties in proton beam radiotherapy for gastrointestinal cancers

In recent years, there has been rapid adaption of proton beam radiotherapy (RT) for treatment of various malignancies in the gastrointestinal (GI) tract, with increasing number of institutions implementing intensity modulated proton therapy (IMPT). We review the progress and existing literature regarding the technical aspects of RT planning for IMPT, and the existing tools that can help with the management of uncertainties which may impact the daily delivery of proton therapy. We provide an in-depth discussion regarding range uncertainties, dose calculations, image guidance requirements, organ and body cavity filling consideration, implanted devices and hardware, use of fiducials, breathing motion evaluations and both active and passive motion management methods, interplay effect, general IMPT treatment planning considerations including robustness plan evaluation and optimization, and finally plan monitoring and adaptation. These advances have improved confidence in delivery of IMPT for patients with GI malignancies under various scenarios.

Robust Optimization for Intensity Modulated Proton Therapy to Redistribute High Linear Energy Transfer from Nearby Critical Organs to Tumors in Head and Neck Cancer

PURPOSE

We propose linear energy transfer (LET)-guided robust optimization in intensity modulated proton therapy for head and neck cancer. This method simultaneously considers LET and physical dose distributions of tumors and organs at risk (OARs) with uncertainties.

METHODS AND MATERIALS

Fourteen patients with head and neck cancer were included in this retrospective study. Cord, brain stem, brain, and oral cavity were considered. Two algorithms, voxel-wise worst-case robust optimization and LET-guided robust optimization (LETRO), were used to generate intensity modulated proton therapy plans for each patient. The latter method directly optimized LET distributions rather than indirectly as in previous methods. LET-volume histograms (LETVHs) were generated, and high LET was redistributed from nearby OARs to tumors in a user-defined way via LET-volume constraints. Dose-volume histogram indices, such as clinical target volume (CTV) D_98% and D_2%-D_98%, cord D_max, brain stem D_max, brain D_max, and oral cavity D_mean, were calculated. Plan robustness was quantified using the worst-case analysis method. LETVH indices analogous to dose-volume histogram indices were used to characterize LET distributions. The Wilcoxon signed rank test was performed to measure statistical significance.

RESULTS

In the nominal scenario, LETRO provided higher LET distributions in the CTV (unit: keV/μm; CTV LET_98%: 1.18 vs 1.08, LETRO vs RO, P = .0031) while preserving comparable physical dose and plan robustness. LETRO achieved significantly reduced LET distributions in the cord, brain stem, and oral cavity compared with RO (cord LET_max: 7.20 vs 8.20, P = .0010; brain stem LET_max: 10.95 vs 12.05, P = .0007; oral cavity LET_mean: 2.11 vs 3.12, P = .0052) and had comparable physical dose and plan robustness in all OARs. In the worst-case scenario, LETRO achieved significantly higher LET_mean in the CTV, reduced LET_max in the brain, and was comparable to other LETVH indices (CTV LET_mean: 3.26 vs 3.35, P = .0012; brain LET_max: 24.80 vs 22.00, P = .0016).

CONCLUSIONS

LETRO robustly optimized LET and physical dose distributions simultaneously. It redistributed high LET from OARs to targets with slightly modified physical dose and plan robustness.

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

Purpose: Pediatric craniopharyngioma, adult base-of-skull sarcoma and chordoma cases are all regarded as priority candidates for proton therapy. In this study, a dosimetric comparison between volumetric modulated arc therapy (VMAT) and intensity modulated proton therapy (IMPT) was first performed. We then investigated the impact of physical and biological uncertainties. We assessed whether IMPT plans remained dosimetrically superior when such uncertainty estimates were considered, especially with regards to sparing organs at risk (OARs).Methodology: We studied 10 cases: four chondrosarcoma, two chordoma and four pediatric craniopharyngioma. VMAT and IMPT plans were created according to modality-specific protocols. For IMPT, we considered (i) variable RBE modeling using the McNamara model for different values of (α/β)x, and (ii) robustness analysis with ±3 mm set-up and 3.5% range uncertainties.Results: When comparing the VMAT and IMPT plans, the dosimetric advantages of IMPT were clear: IMPT led to reduced integral dose and, typically, improved CTV coverage given our OAR constraints. When physical robustness analysis was performed for IMPT, some uncertainty scenarios worsened the CTV coverage but not usually beyond that achieved by VMAT. Certain scenarios caused OAR constraints to be exceeded, particularly for the brainstem and optical chiasm. However, variable RBE modeling predicted even more substantial hotspots, especially for low values of (α/β)x. Variable RBE modeling often prompted dose constraints to be exceeded for critical structures.Conclusion: For base-of-skull and pediatric craniopharyngioma cases, both physical and biological robustness analyses should be considered for IMPT: these analyses can substantially affect the sparing of OARs and comparisons against VMAT. All proton RBE modeling is subject to high levels of uncertainty, but the clinical community should remain cognizant possible RBE effects. Careful clinical and imaging follow-up, plus further research on end-of-range RBE mitigation strategies such as LET optimization, should be prioritized for these cohorts of proton patients.

Technical Note: Treatment planning system (TPS) approximations matter - comparing intensity-modulated proton therapy (IMPT) plan quality and robustness between a commercial and an in-house developed TPS for nonsmall cell lung cancer (NSCLC)

PURPOSE

Approximate dose calculation methods were used in the nominal dose distribution and the perturbed dose distributions due to uncertainties in a commercial treatment planning system (CTPS) for robust optimization in intensity-modulated proton therapy (IMPT). We aimed to investigate whether the approximations influence plan quality, robustness, and interplay effect of the resulting IMPT plans for the treatment of locally advanced lung cancer patients.

MATERIALS AND METHODS

Ten consecutively treated locally advanced nonsmall cell lung cancer (NSCLC) patients were selected. Two IMPT plans were created for each patient using our in-house developed TPS, named "Solo," and also the CTPS, Eclipse^TM (Varian Medical Systems, Palo Alto, CA, USA), respectively. The plans were designed to deliver prescription doses to internal target volumes (ITV) drawn by a physician on averaged four-dimensional computed tomography (4D-CT). Solo plans were imported back to CTPS, and recalculated in CTPS for fair comparison. Both plans were further verified for each patient by recalculating doses in the inhalation and exhalation phases to ensure that all plans met clinical requirements. Plan robustness was quantified on all phases using dose-volume-histograms (DVH) indices in the worst-case scenario. The interplay effect was evaluated for every plan using an in-house developed software, which randomized starting phases of each field per fraction and accumulated dose in the exhalation phase based on the patient's breathing motion pattern and the proton spot delivery in a time-dependent fashion. DVH indices were compared using Wilcoxon rank-sum test.

RESULTS

Compared to the plans generated using CTPS on the averaged CT, Solo plans had significantly better target dose coverage and homogeneity (normalized by the prescription dose) in the worst-case scenario [ITV D_95% : 98.04% vs 96.28%, Solo vs CTPS, P = 0.020; ITV D_5% -D_95% : 7.20% vs 9.03%, P = 0.049] while all DVH indices were comparable in the nominal scenario. On the inhalation phase, Solo plans had better target dose coverage and cord D_max in the nominal scenario [ITV D_95% : 99.36% vs 98.45%, Solo vs CTPS, P = 0.014; cord D_max : 20.07 vs 23.71 Gy(RBE), P = 0.027] with better target coverage and cord D_max in the worst-case scenario [ITV D_95% : 97.89% vs 96.47%, Solo vs CTPS, P = 0.037; cord D_max : 24.57 vs 28.14 Gy(RBE), P = 0.037]. On the exhalation phase, similar phenomena were observed in the nominal scenario [ITV D_95% : 99.63% vs 98.87%, Solo vs CTPS, P = 0.037; cord D_max : 19.67 vs 23.66 Gy(RBE), P = 0.039] and in the worst-case scenario [ITV D_95% : 98.20% vs 96.74%, Solo vs CTPS, P = 0.027; cord D_max : 23.47 vs 27.93 Gy(RBE), P = 0.027]. In terms of interplay effect, plans generated by Solo had significantly better target dose coverage and homogeneity, less hot spots, and lower esophageal D_mean , and cord D_max [ITV D_95% : 101.81% vs 98.68%, Solo vs CTPS, P = 0.002; ITV D_5% -D_95% : 2.94% vs 7.51%, P = 0.002; cord D_max : 18.87 vs 22.29 Gy(RBE), P = 0.014].

CONCLUSIONS

Solo-generated IMPT plans provide improved cord sparing, better target robustness in all considered phases, and reduced interplay effect compared with CTPS. Consequently, the approximation methods currently used in commercial TPS programs may have space for improvement in generating optimal IMPT plans for patient cases with locally advanced lung cancer.

Lower doses to hippocampi and other brain structures for skull-base meningiomas with intensity modulated proton therapy compared to photon therapy

BACKGROUND AND PURPOSE

Radiotherapy of skull-base meningiomas is challenging due to the close proximity of multiple sensitive organs at risk (OARs). This study systematically compared intensity modulated proton therapy (IMPT), non-coplanar volumetric modulated arc therapy (VMAT) and intensity modulated radiotherapy (IMRT) based on automated treatment planning. Differences in OARs sparing, with specific focus on the hippocampi, and low-dose delivery were quantified.

MATERIALS AND METHODS

Twenty patients, target diameter >3 cm, were included. Automated plan generation was used to calculate a VMAT plan with three non-coplanar arcs, an IMRT plan with nine non-coplanar beams with optimized gantry and couch angles, and an IMPT plan with three patient-specific selected non-coplanar beams. A prescription dose of 50.4 GyRBE in 28 fractions was used. The same set of constraints and prioritized objectives was used. All plans were rescaled to the same target coverage. Repeated measures ANOVA was used to assess the statistical significance of differences in OAR dose parameters between planning techniques.

RESULTS

Compared to VMAT and IMRT, IMPT significantly improved dose conformity to the target volume. Consequently, large dose reductions in OARs were observed. With respect to VMAT, the mean dose and D_40% in the bilateral hippocampus were on average reduced by 48% and 74%, respectively (p ≤ 0.005). With IMPT, the mean dose in the normal brain and volumes receiving 20-30 Gy were up to 47% lower (p ≤ 0.01). When comparing IMPT and IMRT, even larger dose differences in those OARs were observed.

CONCLUSION

For skull-base meningiomas IMPT allows for a considerable dose reduction in the hippocampi, normal brain and other OARs compared to both non-coplanar VMAT and IMRT, which may lead to a clinically relevant reduction of late neurocognitive side effects.

Practical robustness evaluation in radiotherapy - A photon and proton-proof alternative to PTV-based plan evaluation

BACKGROUND AND PURPOSE

A planning target volume (PTV) in photon treatments aims to ensure that the clinical target volume (CTV) receives adequate dose despite treatment uncertainties. The underlying static dose cloud approximation (the assumption that the dose distribution is invariant to errors) is problematic in intensity modulated proton treatments where range errors should be taken into account as well. The purpose of this work is to introduce a robustness evaluation method that is applicable to photon and proton treatments and is consistent with (historic) PTV-based treatment plan evaluations.

MATERIALS AND METHODS

The limitation of the static dose cloud approximation was solved in a multi-scenario simulation by explicitly calculating doses for various treatment scenarios that describe possible errors in the treatment course. Setup errors were the same as the CTV-PTV margin and the underlying theory of 3D probability density distributions was extended to 4D to include range errors, maintaining a 90% confidence level. Scenario dose distributions were reduced to voxel-wise minimum and maximum dose distributions; the first to evaluate CTV coverage and the second for hot spots. Acceptance criteria for CTV D98 and D2 were calibrated against PTV-based criteria from historic photon treatment plans.

RESULTS

CTV D98 in worst case scenario dose and voxel-wise minimum dose showed a very strong correlation with scenario average D98 (R² > 0.99). The voxel-wise minimum dose visualised CTV dose conformity and coverage in 3D in agreement with PTV-based evaluation in photon therapy. Criteria for CTV D98 and D2 of the voxel-wise minimum and maximum dose showed very strong correlations to PTV D98 and D2 (R² > 0.99) and on average needed corrections of -0.9% and +2.3%, respectively.

CONCLUSIONS

A practical approach to robustness evaluation was provided and clinically implemented for PTV-less photon and proton treatment planning, consistent with PTV evaluations but without its static dose cloud approximation.

Dosimetric comparison of distal esophageal carcinoma plans for patients treated with small-spot intensity-modulated proton versus volumetric-modulated arc therapies

BACKGROUND

Esophageal carcinoma is the eighth most common cancer in the world. Volumetric-modulated arc therapy (VMAT) is widely used to treat distal esophageal carcinoma due to high conformality to the target and good sparing of organs at risk (OAR). It is not clear if small-spot intensity-modulated proton therapy (IMPT) demonstrates a dosimetric advantage over VMAT. In this study, we compared dosimetric performance of VMAT and small-spot IMPT for distal esophageal carcinoma in terms of plan quality, plan robustness, and interplay effects.

METHODS

35 distal esophageal carcinoma patients were retrospectively reviewed; 19 patients received small-spot IMPT and the remaining 16 of them received VMAT. Both plans were generated by delivering prescription doses to clinical target volumes (CTVs) on phase-averaged 4D-CT's. The dose-volume-histogram (DVH) band method was used to quantify plan robustness. Software was developed to evaluate interplay effects with randomized starting phases for each field per fraction. DVH indices were compared using Wilcoxon rank-sum test. For fair comparison, all the treatment plans were normalized to have the same CTVhigh D95% in the nominal scenario relative to the prescription dose.

RESULTS

In the nominal scenario, small-spot IMPT delivered statistically significantly lower liver Dmean and V30Gy[RBE] , lung Dmean , heart Dmean compared with VMAT. CTVhigh dose homogeneity and protection of other OARs were comparable between the two treatments. In terms of plan robustness, the IMPT and VMAT plans were comparable for kidney V18Gy[RBE] , liver V30Gy[RBE] , stomach V45Gy[RBE] , lung Dmean , V5Gy[RBE] , and V20Gy[RBE] , cord Dmax and D 0.03 c m 3 , liver Dmean , heart V20Gy[RBE] , and V30Gy[RBE] , but IMPT was significantly worse for CTVhigh D95% , D 2 c m 3 , and D5% -D95% , CTVlow D95% , heart Dmean , and V40Gy[RBE] , requiring careful and experienced adjustments during the planning process and robustness considerations. The small-spot IMPT plans still met the standard clinical requirements after interplay effects were considered.

CONCLUSIONS

Small-spot IMPT decreases doses to heart, liver, and total lung compared to VMAT as well as achieves clinically acceptable plan robustness. Our study supports the use of small-spot IMPT for the treatment of distal esophageal carcinoma.

Robustness Analysis for External Beam Radiation Therapy Treatment Plans: Describing Uncertainty Scenarios and Reporting Their Dosimetric Consequences

PURPOSE

With external beam radiation therapy, uncertainties in treatment planning and delivery can result in an undesirable dose distribution delivered to the patient that can compromise the benefit of treatment. Techniques including geometric margins and probabilistic optimization have been used effectively to mitigate the effects of uncertainties. However, their broad application is inconsistent and can compromise the conclusions derived from cross-technique and cross-modality comparisons.

METHODS AND MATERIALS

Conventional methods to deal with treatment planning and delivery uncertainties are described, and robustness analysis is presented as a framework that is applicable across treatment techniques and modalities.

RESULTS

This report identifies elements that are imperative to include when conducting a robustness analysis and describing uncertainties and their dosimetric effects.

CONCLUSION

The robustness analysis approach described here is presented to promote reliable plan evaluation and dose reporting, particularly during clinical trials conducted across institutions and treatment modalities.

Robust beam orientation optimization for intensity-modulated proton therapy

PURPOSE

Dose conformality and robustness are equally important in intensity modulated proton therapy (IMPT). Despite the obvious implication of beam orientation on both dosimetry and robustness, an automated, robust beam orientation optimization algorithm has not been incorporated due to the problem complexity and paramount computational challenge. In this study, we developed a novel IMPT framework that integrates robust beam orientation optimization (BOO) and robust fluence map optimization (FMO) in a unified framework.

METHODS

The unified framework is formulated to include a dose fidelity term, a heterogeneity-weighted group sparsity term, and a sensitivity regularization term. The L2, 1/2-norm group sparsity is used to reduce the number of active beams from the initial 1162 evenly distributed noncoplanar candidate beams, to between two and four. A heterogeneity index, which evaluates the lateral tissue heterogeneity of a beam, is used to weigh the group sparsity term. With this index, beams more resilient to setup uncertainties are encouraged. There is a symbiotic relationship between the heterogeneity index and the sensitivity regularization; the integrated optimization framework further improves beam robustness against both range and setup uncertainties. This Sensitivity regularization and Heterogeneity weighting based BOO and FMO framework (SHBOO-FMO) was tested on two skull-base tumor (SBT) patients and two bilateral head-and-neck (H&N) patients. The conventional CTV-based optimized plans (Conv) with SHBOO-FMO beams (SHBOO-Conv) and manual beams (MAN-Conv) were compared to investigate the beam robustness of the proposed method. The dosimetry and robustness of SHBOO-FMO plan were compared against the manual beam plan with CTV-based voxel-wise worst-case scenario approach (MAN-WC).

RESULTS

With SHBOO-FMO method, the beams with superior range robustness over manual beams were selected while the setup robustness was maintained or improved. On average, the lowest [D95%, V95%, V100%] of CTV were increased from [93.85%, 91.06%, 70.64%] in MAN-Conv plans, to [98.62%, 98.61%, 96.17%] in SHBOO-Conv plans with range uncertainties. With setup uncertainties, the average lowest [D98%, D95%, V95%, V100%] of CTV were increased from [92.06%, 94.83%, 94.31%, 78.93%] in MAN-Conv plans, to [93.54%, 96.61%, 97.01%, 91.98%] in SHBOO-Conv plans. Compared with the MAN-WC plans, the final SHBOO-FMO plans achieved comparable plan robustness and better OAR sparing, with an average reduction of [Dmean, Dmax] of [6.31, 6.55] GyRBE for the SBT cases and [1.89, 5.08] GyRBE for the H&N cases from the MAN-WC plans.

CONCLUSION

We developed a novel method to integrate robust BOO and robust FMO into IMPT optimization for a unified solution of both BOO and FMO, generating plans with superior dosimetry and good robustness.

Improving dosimetric outcome for hippocampus and cochlea sparing whole brain radiotherapy using spot-scanning proton arc therapy

This feasibility study shows that Spot-scanning Proton Arc therapy (SPArc) is able to significantly reduce the dose to the hippocampus and cochlea compared to both Volumetric Modulated Arc Photon Therapy (VMAT) and the robust optimized Intensity Modulated Proton Therapy (ro-IMPT) plans in whole brain radiotherapy. Furthermore, SPArc not only improves plan robustness but could potentially deliver a treatment as efficient as ro-IMPT when proton system's energy layer switch time is less than 1 s.

Robust optimization for intensity-modulated proton therapy with soft spot sensitivity regularization

PURPOSE

Proton dose distribution is sensitive to uncertainties in range estimation and patient positioning. Currently, the proton robustness is managed by worst-case scenario optimization methods, which are computationally inefficient. To overcome these challenges, we develop a novel intensity-modulated proton therapy (IMPT) optimization method that integrates dose fidelity with a sensitivity term that describes dose perturbation as the result of range and positioning uncertainties.

METHODS

In the integrated optimization framework, the optimization cost function is formulated to include two terms: a dose fidelity term and a robustness term penalizing the inner product of the scanning spot sensitivity and intensity. The sensitivity of an IMPT scanning spot to perturbations is defined as the dose distribution variation induced by range and positioning errors. To evaluate the sensitivity, the spatial gradient of the dose distribution of a specific spot is first calculated. The spot sensitivity is then determined by the total absolute value of the directional gradients of all affected voxels. The fast iterative shrinkage-thresholding algorithm is used to solve the optimization problem. This method was tested on three skull base tumor (SBT) patients and three bilateral head-and-neck (H&N) patients. The proposed sensitivity-regularized method (SenR) was implemented on both clinic target volume (CTV) and planning target volume (PTV). They were compared with conventional PTV-based optimization method (Conv) and CTV-based voxel-wise worst-case scenario optimization approach (WC).

RESULTS

Under the nominal condition without uncertainties, the three methods achieved similar CTV dose coverage, while the CTV-based SenR approach better spared organs at risks (OARs) compared with the WC approach, with an average reduction of [Dmean, Dmax] of [4.72, 3.38]  GyRBE for the SBT cases and [2.54, 3.33] GyRBE for the H&N cases. The OAR sparing of the PTV-based SenR method was comparable with the WC method. The WC method, and SenR approaches all improved the plan robustness from the conventional PTV-based method. On average, under range uncertainties, the lowest [D95%, V95%, V100%] of CTV were increased from [93.75%, 88.47%, 47.37%] in the Conv method, to [99.28%, 99.51%, 86.64%] in the WC method, [97.71%, 97.85%, 81.65%] in the SenR-CTV method and [98.77%, 99.30%, 85.12%] in the SenR-PTV method, respectively. Under setup uncertainties, the average lowest [D95%, V95%, V100%] of CTV were increased from [95.35%, 94.92%, 65.12%] in the Conv method, to [99.43%, 99.63%, 87.12%] in the WC method, [96.97%, 97.13%, 77.86%] in the SenR-CTV method, and [98.21%, 98.34%, 83.88%] in the SenR-PTV method, respectively. The runtime of the SenR optimization is eight times shorter than that of the voxel-wise worst-case method.

CONCLUSION

We developed a novel computationally efficient robust optimization method for IMPT. The robustness is calculated as the spot sensitivity to both range and shift perturbations. The dose fidelity term is then regularized by the sensitivity term for the flexibility and trade-off between the dosimetry and the robustness. In the stress test, SenR is more resilient to unexpected uncertainties. These advantages in combination with its fast computation time make it a viable candidate for clinical IMPT planning.

A novel and individualized robust optimization method using normalized dose interval volume constraints (NDIVC) for intensity-modulated proton radiotherapy

PURPOSE

Intensity-modulated proton therapy (IMPT) is known to be sensitive to patient setup and range uncertainty issues. Multiple robust optimization methods have been developed to mitigate the impact of these uncertainties. Here, we propose a new robust optimization method, which provides an alternative way of robust optimization in IMPT, and is clinically practical, which will enable users to control the balance between nominal plan quality and plan robustness in a user-defined fashion.

METHOD

We calculated nine individual dose distributions which corresponded to one nominal and eight extreme scenarios caused by patient setup and proton beam's range uncertainties. For each voxel, the normalized dose interval (NDI) is defined as the full dose range variation divided by the maximum dose in all uncertainty scenarios (NDI = [max - min dose]/max dose), which was then used to calculate the normalized dose interval volume histogram (NDIVH) curves. The areas under the NDIVH curves were used to quantify plan robustness. A normalized dose interval volume constraint (NDIVC) applied to the target was incorporated to specify the desired robustness which was user-defined. Users could then explore the trade-off between nominal plan quality and plan robustness by adjusting the position of the NDIVCs on the NDIVH curves freely. We benchmarked our method using one lung, five head and neck (H&N), and three prostate cases by comparing our results to those derived using the voxel-wise worst-case robust optimization.

RESULTS

Using the benchmark cases, our new method achieved quality IMPT plans comparable to those derived from the voxel-wise worst-case robust optimization for both nominal plan quality and plan robustness in general; even more conformal and more homogeneous target dose distributions in some cases, if proper NDIVCs were applied. The AUC under NDIVH, as a precise quantitative index of plan robustness, was consistent with DVH bandwidths. Additionally, we demonstrated the feasibility of adjusting the position of NDIVCs in the NDIVH curves which allowed users to explore the trade-off between nominal plan quality and plan robustness.

CONCLUSIONS

The NDIVH-based robust optimization method provided a novel and individualized way of robust optimization in IMPT, and enables users to adjust the balance between nominal plan quality and plan robustness in a user-defined fashion. This method is applicable for continued improvement and developing the next generation of IMPT planning algorithms in the future.

Robust Proton Treatment Planning: Physical and Biological Optimization

Accurate prediction of tumor control and toxicities in radiation therapy faces many uncertainties. Besides interpatient variability in the response to radiation, there are also dosimetric uncertainties, that is, differences between the dose displayed in a treatment planning system and the dose actually delivered to the patient. These uncertainties originate from several sources including imperfect knowledge of the patient geometry, approximation in the physics of radiation interaction with tissues, and uncertainties in the biological effectiveness of radiation. Generally, uncertainties are considered in the treatment planning process by applying margins. In intensity-modulated radiotherapy (IMRT), this leads to the planning target volume (PTV) concept. Intensity-modulated proton therapy (IMPT) is widely considered as the future of proton therapy. The treatment planning methods for IMPT and IMRT are similar and based on mathematical optimization techniques for both modalities. However, the PTV concept has fundamental limitations in IMPT. Therefore, researchers have developed robust optimization methods that directly incorporate uncertainties into the IMPT optimization problem. In recent years, vendors of commercial planning systems have started to implement these methods so that robust IMPT planning becomes available in clinical practice. This article summarizes uncertainties in proton therapy and the limitations of the PTV concept to deal with them. Subsequently, robust optimization techniques to overcome these limitations are reviewed.

Small-spot intensity-modulated proton therapy and volumetric-modulated arc therapies for patients with locally advanced non-small-cell lung cancer: A dosimetric comparative study

Purpose

To compare dosimetric performance of volumetric‐modulated arc therapy (VMAT) and small‐spot intensity‐modulated proton therapy for stage III non‐small‐cell lung cancer (NSCLC).

Methods and Materials

A total of 24 NSCLC patients were retrospectively reviewed; 12 patients received intensity‐modulated proton therapy (IMPT) and the remaining 12 received VMAT. Both plans were generated by delivering prescription doses to clinical target volumes (CTV) on averaged 4D‐CTs. The dose‐volume‐histograms (DVH) band method was used to quantify plan robustness. Software was developed to evaluate interplay effects with randomized starting phases of each field per fraction. DVH indices were compared using Wilcoxon rank sum test.

Results

Compared with VMAT, IMPT delivered significantly lower cord D_max, heart D_mean, and lung V_5 Gy[RBE] with comparable CTV dose homogeneity, and protection of other OARs. In terms of plan robustness, the IMPT plans were statistically better than VMAT plans in heart D_mean, but were statistically worse in CTV dose coverage, cord D_max, lung D_mean, and V_5 Gy[RBE]. Other DVH indices were comparable. The IMPT plans still met the standard clinical requirements with interplay effects considered.

Conclusions

Small‐spot IMPT improves cord, heart, and lung sparing compared to VMAT and achieves clinically acceptable plan robustness at least for the patients included in this study with motion amplitude less than 11 mm. Our study supports the usage of IMPT to treat some lung cancer patients.

Impact of bowel gas and body outline variations on total accumulated dose with intensity-modulated proton therapy in locally advanced cervical cancer patients

INTRODUCTION

Density changes occurring during fractionated radiotherapy in the pelvic region may degrade proton dose distributions. The aim of the study was to quantify the dosimetric impact of gas cavities and body outline variations.

MATERIAL AND METHODS

Seven patients with locally advanced cervical cancer (LACC) were analyzed through a total of 175 daily cone beam computed tomography (CBCT) scans. Four-beams intensity-modulated proton therapy (IMPT) dose plans were generated targeting the internal target volume (ITV) composed of: primary tumor, elective and pathological nodes. The planned dose was 45 Gy [Relative-Biological-Effectiveness-weighted (RBE)] in 25 fractions and simultaneously integrated boosts of pathologic lymph nodes were 55-57.5 Gy (RBE). In total, 475 modified CTs were generated to evaluate the effect of: 1/gas cavities, 2/outline variations and 3/the two combined. The anatomy of each fraction was simulated by propagating gas cavities contours and body outlines from each daily CBCT to the pCT. Hounsfield units corresponding to gas and fat were assigned to the propagated contours. D98 (least dose received by the hottest 98% of the volume) and D99.9 for targets and V43Gy(RBE) (volume receiving ≥43 Gy(RBE)) for organs at risk (OARs) were recalculated on each modified CT, and total dose was evaluated through dose volume histogram (DVH) addition across all fractions.

RESULTS

Weight changes during radiotherapy were between -3.1% and 1.2%. Gas cavities and outline variations induced a median [range] dose degradation for ITV45 of 1.0% [0.5-3.5%] for D98 and 2.1% [0.8-6.4%] for D99.9. Outline variations had larger dosimetric impact than gas cavities. Worst nodal dose degradation was 2.0% for D98 and 2.3% for D99.9. The impact on bladder, bowel and rectum was limited with V43Gy(RBE) variations ≤3.5 cm³.

CONCLUSION

Bowel gas cavities and outline variations had minor impact on accumulated dose in targets and OAR of four-field IMPT in a LACC population of moderate weight changes.

The role of image-guided intensity modulated proton therapy in glioma

Gliomas represent a broad spectrum of disease with life-expectancy outcomes ranging from months to decades. As our understanding of the molecular profiles of gliomas expands rapidly, practitioners are now better able to identify patients with favorable versus nonfavorable prognoses. Radiation therapy plays a key role in glioma treatment, improving disease control and oftentimes survival. However, for survivors, either long-term or short-term, radiation-induced cognitive impairments may negatively impact their quality of life. For patients with both favorable and unfavorable prognoses, intensity modulated proton therapy (IMPT) may offer significant, yet unproven benefits. IMPT is the newest and most advanced proton delivery technique, one with substantial benefits compared with historical proton techniques. IMPT allows practitioners to maximize the physical benefits of protons, increasing normal tissue sparing and reducing the potential for adverse effects. For more aggressive tumors, the dose conformality and normal tissue sparing afforded by IMPT may also allow for dose escalation to target volumes. However, in order to truly maximize the clinical potential of IMPT, the field of radiation oncology must not only implement the most advanced technologies, but also understand and capitalize on the unique biologic aspects of proton therapy.

Exploratory study of the association of volumetric modulated arc therapy (VMAT) plan robustness with local failure in head and neck cancer

This work is to show which is more relevant to cause local failures (LFs) due to patient setup uncertainty between the planning target volume (PTV) underdosage and the potential target underdosage subject to patient setup uncertainties in head and neck (H&N) cancer treated with volumetric‐modulated arc therapy (VMAT). Thirteen LFs in 10 H&N patients treated by VMAT were analyzed. Measures have been taken to minimize the chances of insufficient target delineation for these patients and the patients were clinically determined to have LF based on the PET/CT scan results by an experienced radiologist and then reviewed by a second experienced radiation oncologist. Two methods were used to identify the possible locations of LF due to underdosage: (a) examining the standard VMAT plan, in which the underdosed volume in the nominal dose distribution (UVN) was generated by subtracting the volumes receiving the prescription doses from PTVs, and (b) plan robustness analysis, in which in addition to the nominal dose distribution, six perturbed dose distributions were created by translating the CT iso‐center in three cardinal directions by the PTV margin. The coldest dose distribution was represented by the minimum of the seven doses in each voxel. The underdosed volume in the coldest dose distribution (UVC) was generated by subtracting the volumes receiving the prescription doses in the coldest dose distribution from the volumes receiving the prescription doses in the nominal dose distribution. UVN and UVC were subsequently examined for spatial association with the locations of LF. The association was tested using the binominal distribution and the Fisher's exact test of independence. We found that of 13 LFs, 11 were associated with UVCs (P = 0.011), while three were associated with UVNs (P = 0.99). We concluded that the possible target underdosage due to patient setup uncertainties appeared to be a more relevant factor associated with LF in VMAT for H&N cancer than the compromised PTV coverage at least for the patients included in this study.

Proton therapy - Present and future

In principle, proton therapy offers a substantial clinical advantage over conventional photon therapy. This is because of the unique depth-dose characteristics of protons, which can be exploited to achieve significant reductions in normal tissue doses proximal and distal to the target volume. These may, in turn, allow escalation of tumor doses and greater sparing of normal tissues, thus potentially improving local control and survival while at the same time reducing toxicity and improving quality of life. Protons, accelerated to therapeutic energies ranging from 70 to 250MeV, typically with a cyclotron or a synchrotron, are transported to the treatment room where they enter the treatment head mounted on a rotating gantry. The initial thin beams of protons are spread laterally and longitudinally and shaped appropriately to deliver treatments. Spreading and shaping can be achieved by electro-mechanical means to treat the patients with "passively-scattered proton therapy" (PSPT) or using magnetic scanning of thin "beamlets" of protons of a sequence of initial energies. The latter technique can be used to treat patients with optimized intensity modulated proton therapy (IMPT), the most powerful proton modality. Despite the high potential of proton therapy, the clinical evidence supporting the broad use of protons is mixed. It is generally acknowledged that proton therapy is safe, effective and recommended for many types of pediatric cancers, ocular melanomas, chordomas and chondrosarcomas. Although promising results have been and continue to be reported for many other types of cancers, they are based on small studies. Considering the high cost of establishing and operating proton therapy centers, questions have been raised about their cost effectiveness. General consensus is that there is a need to conduct randomized trials and/or collect outcomes data in multi-institutional registries to unequivocally demonstrate the advantage of protons. Treatment planning and plan evaluation of PSPT and IMPT require special considerations compared to the processes used for photon treatment planning. The differences in techniques arise from the unique physical properties of protons but are also necessary because of the greater vulnerability of protons to uncertainties, especially from inter- and intra-fractional variations in anatomy. These factors must be considered in designing as well as evaluating treatment plans. In addition to anatomy variations, other sources of uncertainty in dose delivered to the patient include the approximations and assumptions of models used for computing dose distributions for planning of treatments. Furthermore, the relative biological effectiveness (RBE) of protons is simplistically assumed to have a constant value of 1.1. In reality, the RBE is variable and a complex function of the energy of protons, dose per fraction, tissue and cell type, end point, etc. These uncertainties, approximations and current technological limitations of proton therapy may limit the achievement of its true potential. Ongoing research is aimed at better understanding the consequences of the various uncertainties on proton therapy and reducing the uncertainties through image-guidance, adaptive radiotherapy, further study of biological properties of protons and the development of novel dose computation and optimization methods. However, residual uncertainties will remain in spite of the best efforts. To increase the resilience of dose distributions in the face of uncertainties and improve our confidence in dose distributions seen on treatment plans, robust optimization techniques are being developed and implemented. We assert that, with such research, proton therapy will be a commonly applied radiotherapy modality for most types of solid cancers in the near future.