Effectiveness of robust optimization in intensity-modulated proton therapy planning for head and neck cancers

Interval Analysis-Based Optimization: A Robust Model for Intensity-Modulated Radiotherapy (IMRT)

Background: Cancer remains one of the leading causes of mortality worldwide, with radiotherapy playing a crucial role in its treatment. Intensity-modulated radiotherapy (IMRT) enables precise dose delivery to tumors while sparing healthy tissues. However, geometric uncertainties such as patient positioning errors and anatomical deformations can compromise treatment accuracy. Traditional methods use safety margins, which may lead to excessive irradiation of healthy organs or insufficient tumor coverage. Robust optimization techniques, such as minimax approaches, attempt to address these uncertainties but can result in overly conservative treatment plans. This study introduces an interval analysis-based optimization model for IMRT, offering a more flexible approach to uncertainty management. Methods: The proposed model represents geometric uncertainties using interval dose influence matrices and incorporates Bertoluzza's metric to balance tumor coverage and organ-at-risk (OAR) protection. The θ parameter allows controlled robustness modulation. The model was implemented in matRad, an open-source treatment planning system, and evaluated on five prostate cancer cases. Results were compared against traditional Planning Target Volume (PTV) and minimax robust optimization approaches. Results: The interval-based model improved tumor coverage by 5.8% while reducing bladder dose by 4.2% compared to PTV. In contrast, minimax robust optimization improved tumor coverage by 25.8% but increased bladder dose by 23.2%. The interval-based approach provided a better balance between tumor coverage and OAR protection, demonstrating its potential to enhance treatment effectiveness without excessive conservatism. Conclusions: This study presents a novel framework for IMRT planning that improves uncertainty management through interval analysis. By allowing adjustable robustness modulation, the proposed model enables more personalized and clinically adaptable treatment plans. These findings highlight the potential of interval analysis as a powerful tool for optimizing radiotherapy outcomes, balancing treatment efficacy and patient safety.

Quantification of beam size impact on intensity-modulated proton therapy with robust optimization in head and neck cancer-comparison with intensity-modulated radiation therapy

We assessed the effect of beam size on plan robustness for intensity-modulated proton therapy (IMPT) of head and neck cancer (HNC) and compared the plan quality including robustness with that of intensity-modulated radiation therapy (IMRT). IMPT plans were generated for six HNC patients using six beam sizes (air-sigma 3-17 mm at isocenter for a 70-230 MeV) and two optimization methods for planning target volume-based non-robust optimization (NRO) and clinical target volume (CTV)-based robust optimization (RO). Worst-case dosimetric parameters and plan robustness for CTV and organs-at-risk (OARs) were assessed under different scenarios, assuming a ± 1-5 mm setup error and a ± 3% range error. Statistical comparisons of NRO-IMPT, RO-IMPT and IMRT plans were performed. In regard to CTV-D99%, RO-IMPT with smaller beam size was more robust than RO-IMPT with larger beam sizes, whereas NRO-IMPT showed the opposite (P < 0.05). There was no significant difference in the robustness of the CTV-D99% and CTV-D95% between RO-IMPT and IMRT. The worst-case CTV coverage of IMRT (±5 mm/3%) for all patients was 96.0% ± 1.4% (D99%) and 97.9% ± 0.3% (D95%). For four out of six patients, the worst-case CTV-D95% for RO-IMPT (±1-5 mm/3%) were higher than those for IMRT. Compared with IMRT, RO-IMPT with smaller beam sizes achieved lower worst-case doses to OARs. In HNC treatment, utilizing smaller beam sizes in RO-IMPT improves plan robustness compared to larger beam sizes, achieving comparable target robustness and lower worst-case OARs doses compared to IMRT.

An AI dose-influence matrix engine for robust pencil beam scanning protons therapy

BACKGROUND

Rapid planning is of tremendous value in proton pencil beam scanning (PBS) therapy in overcoming range uncertainty. However, the dose calculation of the dose influence matrix (Dij) in robust PBS plan optimization is time-consuming and requires substantial acceleration to enhance efficiency.

PURPOSE

To accelerate the Dij calculations in PBS therapy, we developed an AI-Dij engine integrated into our in-house treatment planning system (TPS).

METHODS

The AI-Dij engine calculates spot dose using a transformer-based spot dose calculation model (SDM), which takes CT volumes (CT-bars, 256  × $ \times $  16  × $ \times $  16 voxels, 3 mm resolution) and energy (a float value) as inputs and outputs the spot dose distribution (256  × $ \times $  16  × $ \times $  16). The SDM was trained on over 200 000 CT-bars and Monte Carlo (MC) spot dose (spanning energy levels from 70 to 225 MeV). Clinical-implemented treatment plans for the head, lung, and liver, initially created on Raystation, were replanned using our AI-Dij engine under identical gantry angles and uncertainties settings. After optimizing the spot weight, each in-house plan was recalculated using MCsquare for MC dose evaluation. The dose-volume histogram (DVH) metrics from the in-house TPS and Raystation were compared, evaluating both the optimized and MC doses.

RESULTS

In optimization, the differences of DVH metrics (%, Valuein-house-ValueRaystation) across all uncertainty scenarios between the in-house and Raystation plans were 0.93 ± 2.04% for clinical target volume (CTV) and -5.94 ± 12.19% for organ at risks (OARs). For the MC doses, the differences were 2.48 ± 2.78% for CTV and -5.47 ± 14.16% for OARs. The time cost of a robust AI-Dij calculation can be within 2s on an RTX3090 GPU.

CONCLUSION

We conducted a feasibility study on AI-Dij engine-based robust PBS plan optimization, demonstrating both high planning speed and quality.

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.

Exploring the dosimetric impact of systematic and random setup uncertainties in robust optimization of head and neck IMPT plans

PURPOSE

This study aims to compare the dosimetric impact of incorporating systematic and random setup uncertainties in the robust optimization of head and neck cancer (HNC) Intensity Modulated Proton Therapy (IMPT) plans.

METHODS

Bilateral HNC patients (n = 10) previously treated with conventional photon therapy at our institution were included. Both systematic and random setup uncertainties were incorporated into the robust optimization process of IMPT planning. Dosimetric comparisons were made between plans optimized with systematic (IMPT-S) versus random (IMPT-R) setup uncertainties, assessing both the clinical target volume (CTVs) and organs at risk (OARs) across various dosimetric metrics. Both plans applied a fixed range uncertainty of ± 3 % and a maximum setup uncertainty of ± 3 mm.

RESULTS

Both IMPT-S and IMPT-R plans achieved similar target coverage, meeting robustness criteria for CTVs. On average, the D95% voxel-wise min to the high-risk CTV (CTV_HR) was slightly higher in IMPT-S plans by 1.78 ± 0.72 % compared to IMPT-R plans. However, IMPT-R plans provided better OAR sparing, which was evident in both nominal and voxel-wise maximum values. While random setup errors in robust optimization improved OAR sparing, the clinical impact may be minimal where OAR doses are already well below tolerance levels.

CONCLUSION

Both IMPT-S and IMPT-R techniques met the robustness criteria for CTVs in HNC IMPT planning. Incorporating random setup uncertainties in robust optimization improves OAR sparing compared to systematic setup uncertainties. Further research is needed to explore the broader applicability of random setup errors and to integrate random uncertainties in robustness evaluations for a more comprehensive assessment of treatment plans.

Technical note: Optimizing spot-scanning proton arc therapy with a novel spot sparsity approach

BACKGROUND

One of the main challenges of utilizing spot-scanning proton arc therapy (SPArc) in routine clinics is treatment delivery efficiency. Spot reduction, which relies on spot sparsity optimization (SSO), is crucial for achieving high delivery efficiency in SPArc.

PURPOSE

This study aims to develop a novel SSO approach based on the alternating directions method of multipliers (ADMM) for SPArc to achieve high treatment delivery efficiency and maintain optimal dosimetric plan quality.

METHODS

In this study, SSO for SPArc is based on the least-square dose fidelity term with L0-norm regularization. The novel optimization approach is based on the ADMM framework, in which the minimum monitor unit constraint was considered to improve the plan quality. A state-of-the-art SSO method, the primal-dual active set with continuation (PDASC) algorithm published previously, was utilized as a benchmark. Two SPArc plan groups with the same beam assignment and clinical constraint were generated, in which the former group was SPArc plan with SSO utilizing ADMM, denoted asSPArc ADMM $\text{SPArc}_{\text{ADMM}}$ , and the later group was SPArc with SSO utilizing PDASC, denoted asSPArc PDASC $\text{SPArc}_{\text{PDASC}}$ . Nine clinical cases included five different cancer sites (brain, lung, liver, prostate, and head&neck cancer) were used. The SSO method's performance was evaluated in terms of spot sparsity level (the number of zero-valued elements divided by the total number of elements), beam delivery time, dosimetric plan quality, and plan robustness.

RESULTS

Compared to theSPArc PDASC $\text{SPArc}_{\text{PDASC}}$ plan, theSPArc ADMM $\text{SPArc}_{\text{ADMM}}$ plan exhibits superior sparsity and higher delivery efficiency while maintaining good plan quality.

CONCLUSIONS

This study introduces a novel spot sparsity optimization approach using the ADMM framework to improve the delivery efficiency of SPArc. Compared to the existing state-of-the-art SSO method, such an approach could further enhance delivery efficiency while maintaining good plan quality, which could promote the implementation of SPArc in the clinic's.

Dosimetric Evaluation of Beam-specific PTV and Worst-case Optimization Methods for Liver Proton Therapy

BACKGROUND/AIM

In spot-scanning proton therapy, intra-fractional anatomical changes by organ movement can lead to deterioration in dose distribution due to beam range variation. To explore a more robust treatment planning method, this study evaluated the dosimetric characteristics and robustness of two proton therapy planning methods for liver cancer.

PATIENTS AND METHODS

Two- or three-field treatment plans were created for 11 patients with hepatocellular carcinoma or metastatic liver cancer using a single-field uniform dose (SFUD) technique. The plans were optimized using either beam-specific planning target volume (BSPTV) or worst-case optimization (WCO). The target coverage for the gross tumor volume (GTV), planning target volume (PTV), and organs at risk (OAR) parameters related to toxicity were calculated from the perturbed dose distributions, considering setup and range uncertainties. Statistical analyses of the BSPTV and WCO plans were performed using the Wilcoxon signed-rank sum test (p<0.05). The calculation times for a single optimization process were also recorded and compared.

RESULTS

The robustness of the WCO plans in the worst-case scenario was significantly higher than that of the BSPTV plan in terms of GTV target coverage, prevention of maximum dose increase to the gastrointestinal tract, and the dose received by normal liver regions. However, there were no significant differences in PTV, and the calculation time required to create the WCO plan was considerably longer.

CONCLUSION

In SFUD proton therapy for liver cancer, the WCO plans required a longer optimization time but exhibited superior robustness in GTV coverage and sparing of OARs.

Definition of a framework for volumetric modulated arc therapy plan quality assessment with integration of dose-, complexity-, and robustness metrics

Background and purpose

Conventionally, the quality of radiotherapy treatment plans is assessed through visual inspection of dose distributions and dose-volume histograms. This study developed a framework to evaluate plan quality using dose, complexity, and robustness metrics. Additionally, a method for predicting plan robustness metrics using dose and complexity metrics was introduced for cases where plan robustness evaluation is unavailable or impractical.

Materials and methods

The framework and prediction models were developed and validated using 103-bronchial Volumetric Modulated Arc Therapy (VMAT)-plans. The application of the framework was demonstrated using 25-VMAT-plans. To identify significant metrics for plan evaluation, 122-metrics were analysed and narrowed down using multivariate Spearman correlation. Metric limits were set with Statistical process control (SPC). Robustness metrics were predicted using multivariable or single linear regression models based on dose-and complexity-metrics.

Results

Twenty-five-metrics were selected based on the amount and strength of correlations. R95(dose coverage) and HI95/5(homogeneity index) stood out among the dose-metrics, while the complexity-metrics showed similar correlations. Average scenarios dose at 95 % Clinical Target Volume D95mean(CTV) and Errorbar-based Volume-Histograms (EVH) were notable for robustness metrics. Approximately 99 % of evaluated metrics fell within established SPC limits. The prediction model for D95mean(CTV) showed good performance (adjusted R2 = 0.88, mean squared error (MSE) = 3.84 × 10-6), while the model for EVH demonstrated moderate reliability (adjusted R2 = 0.52, MSE = 0.2). No statistically significant differences were found between the predicted (using dose-and complexity-metrics) and calculated robustness metrics (EVH (p-value = 0.9) and D95mean(CTV) (p-value = 1)).

Conclusions

The developed framework enables early detection of sub-optimal, complex and non-robust treatment plans. The predictive model can be used when robustness evaluations are impractical.

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.

First direct machine-specific parameters incorporated in Spot-scanning Proton Arc (SPArc) optimization algorithm

BACKGROUND

Spot-scanning Proton Arc (SPArc) has been of significant interest in recent years because of its superior plan quality. Currently, the primary focus of research and development is on deliverability and treatment efficiency.

PURPOSE

To address the challenges in generating a deliverable and efficient SPArc plan for a proton therapy system with a massive gantry, we developed a novel SPArc optimization algorithm (SPArcDMPO) by directly incorporating the machine-specific parameters such as gantry mechanical constraints and proton delivery sequence.

METHODS

SPArc delivery sequence model (DSMarc) was built based on the machine-specific parameters of the prototype arc delivery system, IBA ProteusONE®, including mechanical constraint (maximum gantry speed, acceleration, and deceleration) and proton delivery sequence (energy and spot delivery sequence, and irradiation time). SPArcDMPO resamples and adjusts each control point's delivery speed based on the DSMarc calculation through the iterative approach. In SPArcDMPO, users could set a reasonable arc delivery time during the plan optimization, which aims to minimize the gantry momentum changes and improve the delivery efficiency. Ten cases were selected to test SPArcDMPO. Two kinds of SPArc plans were generated using the same planning objective functions: (1) original SPArc plan (SPArcoriginal); (2) SPArcDMPO plan with a user-pre-defined delivery time. Additionally, arc delivery sequence was simulated based on the DSMarc and was compared. Treatment delivery time was compared between SPArcoriginal and SPArcDMPO. Dynamic arc delivery time, the static irradiation time, and its corresponding time differential (time differential = dynamic arc delivery time-static irradiation time) were analyzed, respectively. The total gantry velocity change was accumulated throughout the treatment delivery.

RESULTS

With a similar plan quality, objective value, number of energy layers, and spots, both SPArcoriginal and SPArcDMPO plans could be delivered continuously within the ± 1 degree tolerance window. However, compared to the SPArcoriginal, the strategy of SPArcDMPO is able to reduce the time differential from 30.55 ± 11.42%(90 ± 32 s) to 14.67 ± 6.97%(42 ± 20 s), p < 0.01. Furthermore, the corresponding total variations of gantry velocity during dynamic arc delivery are mitigated (SPArcoriginal vs. SPArcDMPO) from 14.73 ± 9.14 degree/s to 4.28 ± 2.42 degree/s, p < 0.01. Consequently, the SPArcDMPO plans could minimize the gantry momentum change based on the clinical user's input compared to the SPArcoriginal plans, which could help relieve the mechanical challenge of accelerating or decelerating the massive proton gantry.

CONCLUSIONS

For the first time, clinical users not only could generate a SPArc plan meeting the mechanical constraint of their proton system but also directly control the arc treatment speed and momentum changes of the gantry during the plan optimization process. This work paved the way for the routine clinical implementation of proton arc therapy in the treatment planning system.

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.

Biological optimization for hybrid proton-photon radiotherapy

Objective.Hybrid proton-photon radiotherapy (RT) is a cancer treatment option to broaden access to proton RT. Additionally, with a refined treatment planning method, hybrid RT has the potential to offer superior plan quality compared to proton-only or photon-only RT, particularly in terms of target coverage and sparing organs-at-risk (OARs), when considering robustness to setup and range uncertainties. However, there is a concern regarding the underestimation of the biological effect of protons on OARs, especially those in close proximity to targets. This study seeks to develop a hybrid treatment planning method with biological dose optimization, suitable for clinical implementation on existing proton and photon machines, with each photon or proton treatment fraction delivering a uniform target dose.Approach.The proposed hybrid biological dose optimization method optimized proton and photon plan variables, along with the number of fractions for each modality, minimizing biological dose to the OARs and surrounding normal tissues. To mitigate underestimation of hot biological dose spots, proton biological dose was minimized within a ring structure surrounding the target. Hybrid plans were designed to be deliverable separately and robustly on existing proton and photon machines, with enforced uniform target dose constraints for the proton and photon fraction doses. A probabilistic formulation was utilized for robust optimization of setup and range uncertainties for protons and photons. The nonconvex optimization problem, arising from minimum monitor unit constraint and dose-volume histogram constraints, was solved using an iterative convex relaxation method.Main results.Hybrid planning with biological dose optimization effectively eliminated hot spots of biological dose, particularly in normal tissues surrounding the target, outperforming proton-only planning. It also provided superior overall plan quality and OAR sparing compared to proton-only or photon-only planning strategies.Significance.This study presents a novel hybrid biological treatment planning method capable of generating plans with reduced biological hot spots, superior plan quality to proton-only or photon-only plans, and clinical deliverability on existing proton and photon machines, separately and robustly.

A retrospective study on the investigation of potential dosimetric benefits of online adaptive proton therapy for head and neck cancer

PURPOSE

Proton therapy is sensitive to anatomical changes, often occurring in head-and-neck (HN) cancer patients. Although multiple studies have proposed online adaptive proton therapy (APT), there is still a concern in the radiotherapy community about the necessity of online APT. We have performed a retrospective study to investigate the potential dosimetric benefits of online APT for HN patients relative to the current offline APT.

METHODS

Our retrospective study has a patient cohort of 10 cases. To mimic online APT, we re-evaluated the dose of the in-use treatment plan on patients' actual treatment anatomy captured by cone-beam CT (CBCT) for each fraction and performed a templated-based automatic replanning if needed, assuming that these were performed online before treatment delivery. Cumulative dose of the simulated online APT course was calculated and compared with that of the actual offline APT course and the designed plan dose of the initial treatment plan (referred to as nominal plan). The ProKnow scoring system was employed and adapted for our study to quantify the actual quality of both courses against our planning goals.

RESULTS

The average score of the nominal plans over the 10 cases is 41.0, while those of the actual offline APT course and our simulated online course is 25.8 and 37.5, respectively. Compared to the offline APT course, our online course improved dose quality for all cases, with the score improvement ranging from 0.4 to 26.9 and an average improvement of 11.7.

CONCLUSION

The results of our retrospective study have demonstrated that online APT can better address anatomical changes for HN cancer patients than the current offline replanning practice. The advanced artificial intelligence based automatic replanning technology presents a promising avenue for extending potential benefits of online APT.

A fast and robust constraint-based online re-optimization approach for automated online adaptive intensity modulated proton therapy in head and neck cancer

Objective. In head-and-neck cancer intensity modulated proton therapy, adaptive radiotherapy is currently restricted to offline re-planning, mitigating the effect of slow changes in patient anatomies. Daily online adaptations can potentially improve dosimetry. Here, a new, fully automated online re-optimization strategy is presented. In a retrospective study, this online re-optimization approach was compared to our trigger-based offline re-planning (offlineTBre-planning) schedule, including extensive robustness analyses.Approach. The online re-optimization method employs automated multi-criterial re-optimization, using robust optimization with 1 mm setup-robustness settings (in contrast to 3 mm for offlineTBre-planning). Hard planning constraints and spot addition are used to enforce adequate target coverage, avoid prohibitively large maximum doses and minimize organ-at-risk doses. For 67 repeat-CTs from 15 patients, fraction doses of the two strategies were compared for the CTVs and organs-at-risk. Per repeat-CT, 10.000 fractions with different setup and range robustness settings were simulated using polynomial chaos expansion for fast and accurate dose calculations.Main results. For 14/67 repeat-CTs, offlineTBre-planning resulted in <50% probability ofD98%≥ 95% of the prescribed dose (Dpres) in one or both CTVs, which never happened with online re-optimization. With offlineTBre-planning, eight repeat-CTs had zero probability of obtainingD98%≥ 95%Dpresfor CTV7000, while the minimum probability with online re-optimization was 81%. Risks of xerostomia and dysphagia grade ≥ II were reduced by 3.5 ± 1.7 and 3.9 ± 2.8 percentage point [mean ± SD] (p< 10-5for both). In online re-optimization, adjustment of spot configuration followed by spot-intensity re-optimization took 3.4 min on average.Significance. The fast online re-optimization strategy always prevented substantial losses of target coverage caused by day-to-day anatomical variations, as opposed to the clinical trigger-based offline re-planning schedule. On top of this, online re-optimization could be performed with smaller setup robustness settings, contributing to improved organs-at-risk sparing.

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.

Cochlear-optimized treatment planning in photon and proton radiosurgery for vestibular schwannoma patients

Objective

To investigate the potential to reduce the cochlear dose with robotic photon radiosurgery or intensity-modulated proton therapy planning for vestibular schwannomas.

Materials and Methods

Clinically delivered photon radiosurgery treatment plans were compared to five cochlear-optimized plans: one photon and four proton plans (total of 120). A 1x12 Gy dose was prescribed. Photon plans were generated with Precision (Cyberknife, Accuray) with no PTV margin for set-up errors. Proton plans were generated using an in-house automated multi-criterial planning system with three or nine-beam arrangements, and applying 0 or 3 mm robustness for set-up errors during plan optimization and evaluation (and 3 % range robustness). The sample size was calculated based on a reduction of cochlear Dmean > 1.5 Gy(RBE) from the clinical plans, and resulted in 24 patients.

Results

Compared to the clinical photon plans, a reduction of cochlear Dmean > 1.5 Gy(RBE) could be achieved in 11/24 cochlear-optimized photon plans, 4/24 and 6/24 cochlear-optimized proton plans without set-up robustness for three and nine-beam arrangement, respectively, and in 0/24 proton plans with set-up robustness. The cochlea could best be spared in cases with a distance between tumor and cochlea. Using nine proton beams resulted in a reduced dose to most organs at risk.

Conclusion

Cochlear dose reduction is possible in vestibular schwannoma radiosurgery while maintaining tumor coverage, especially when the tumor is not adjacent to the cochlea. With current set-up robustness, proton therapy is capable of providing lower dose to organs at risk located distant to the tumor, but not for organs adjacent to it. Consequently, photon plans provided better cochlear sparing than proton plans.

Plan robustness and RBE influence for proton dose painting by numbers for head and neck cancers

PURPOSE

To investigate the feasibility of dose painting by numbers (DPBN) with respect to robustness for proton therapy for head and neck cancers (HNC), and to study the influence of variable RBE on the TCP and OAR dose burden.

METHODS AND MATERIALS

Data for 19 patients who have been scanned pretreatment with PET-FDG and subsequently treated with photon therapy were used in the study. A dose response model developed for photon therapy was implemented in a TPS, allowing DPBN plans to be created. Conventional homogeneous dose and DPBN plans were created for each patient, optimized with either fixed RBE = 1.1 or a variable RBE model. Robust optimization was used to create clinically acceptable plans. To estimate the maximum potential loss in TCP due to actual SUV variations from the pre-treatment imaging, we applied a test case with randomized SUV distribution.

RESULTS

Regardless of the use of variable RBE for optimization or evaluation, a statistically significant increase (p < 0.001) in TCP was found for DPBN plans as compared to homogeneous dose plans. Randomizing the SUV distribution decreased the TCP for all plans. A correlation between TCP increase and variance of the SUV distribution and target volume was also found.

CONCLUSION

DPBN for protons and HNC is feasible and could lead to a TCP gain. Risks associated with the temporal variation of SUV distributions could be mitigated by imposing minimum doses to targets. The correlation found between TCP increase and SUV variance and target volume may be used for patient selection.

Deep-learning based fast and accurate 3D CT deformable image registration in lung cancer

BACKGROUND

Deformable Image Registration (DIR) is an essential technique required in many applications of radiation oncology. However, conventional DIR approaches typically take several minutes to register one pair of 3D CT images and the resulting deformable vector fields (DVFs) are only specific to the pair of images used, making it less appealing for clinical application.

PURPOSE

A deep-learning-based DIR method using CT images is proposed for lung cancer patients to address the common drawbacks of the conventional DIR approaches and in turn can accelerate the speed of related applications, such as contour propagation, dose deformation, adaptive radiotherapy (ART), etc. METHODS: A deep neural network based on VoxelMorph was developed to generate DVFs using CT images collected from 114 lung cancer patients. Two models were trained with the weighted mean absolute error (wMAE) loss and structural similarity index matrix (SSIM) loss (optional) (i.e., the MAE model and the M+S model). In total, 192 pairs of initial CT (iCT) and verification CT (vCT) were included as a training dataset and the other independent 10 pairs of CTs were included as a testing dataset. The vCTs usually were taken 2 weeks after the iCTs. The synthetic CTs (sCTs) were generated by warping the vCTs according to the DVFs generated by the pre-trained model. The image quality of the synthetic CTs was evaluated by measuring the similarity between the iCTs and the sCTs generated by the proposed methods and the conventional DIR approaches, respectively. Per-voxel absolute CT-number-difference volume histogram (CDVH) and MAE were used as the evaluation metrics. The time to generate the sCTs was also recorded and compared quantitatively. Contours were propagated using the derived DVFs and evaluated with SSIM. Forward dose calculations were done on the sCTs and the corresponding iCTs. Dose volume histograms (DVHs) were generated based on dose distributions on both iCTs and sCTs generated by two models, respectively. The clinically relevant DVH indices were derived for comparison. The resulted dose distributions were also compared using 3D Gamma analysis with thresholds of 3 mm/3%/10% and 2 mm/2%/10%, respectively.

RESULTS

The two models (wMAE and M+S) achieved a speed of 263.7±163 / 265.8±190 ms and a MAE of 13.15±3.8 / 17.52±5.8 HU for the testing dataset, respectively. The average SSIM scores of 0.987±0.006 and 0.988±0.004 were achieved by the two proposed models, respectively. For both models, CDVH of a typical patient showed that less than 5% of the voxels had a per-voxel absolute CT-number-difference larger than 55 HU. The dose distribution calculated based on a typical sCT showed differences of ≤2cGy[RBE] for clinical target volume (CTV) D95 and D5 , within ±0.06% for total lung V5 , ≤1.5cGy[RBE] for heart and esophagus Dmean , and ≤6cGy[RBE] for cord Dmax compared to the dose distribution calculated based on the iCT. The good average 3D Gamma passing rates (> 96% for 3 mm/3%/10% and > 94% for 2 mm/2%/10%, respectively) were also observed.

CONCLUSION

A deep neural network-based DIR approach was proposed and has been shown to be reasonably accurate and efficient to register the initial CTs and verification CTs in lung cancer.

Robustness of hypofractionated breast radiotherapy after breast-conserving surgery with free breathing

Purpose

This study aimed to evaluate the robustness with respect to the positional variations of five planning strategies in free-breathing breast hypofractionated radiotherapy (HFRT) for patients after breast-conserving surgery.

Methods

Twenty patients who received breast HFRT with 42.72 Gy in 16 fractions were retrospectively analyzed. Five treatment planning strategies were utilized for each patient, including 1) intensity-modulated radiation therapy (IMRT) planning (IMRTpure); 2) IMRT planning with skin flash tool extending and filling the fluence outside the skin by 2 cm (IMRTflash); 3) IMRT planning with planning target volume (PTV) extended outside the skin by 2 cm in the computed tomography dataset (IMRTePTV); 4) hybrid planning, i.e., 2 Gy/fraction three-dimensional conformal radiation therapy combined with 0.67 Gy/fraction IMRT (IMRThybrid); and 5) hybrid planning with skin flash (IMRThybrid-flash). All plans were normalized to 95% PTV receiving 100% of the prescription dose. Six additional plans were created with different isocenter shifts for each plan, which were 1 mm, 2 mm, 3 mm, 5 mm, 7 mm, and 10 mm distally in the X (left-right) and Y (anterior-posterior) directions, namely, (X,Y), to assess their robustness, and the corresponding doses were recalculated. Variation of dosimetric parameters with increasing isocenter shift was evaluated.

Results

All plans were clinically acceptable. In terms of robustness to isocenter shifts, the five planning strategies followed the pattern IMRTePTV, IMRThybrid-flash, IMRTflash, IMRThybrid, and IMRTpure in descending order. V 95% of IMRTePTV maintained at 99.6% ± 0.3% with a (5,5) shift, which further reduced to 98.2% ± 2.0% with a (10,10) shift. IMRThybrid-flash yielded the robustness second to IMRTePTV with less risk from dose hotspots, and the corresponding V 95% maintained >95% up until (5,5).

Conclusion

Considering the dosimetric distribution and robustness in breast radiotherapy, IMRTePTV performed best at maintaining high target coverage with increasing isocenter shift, while IMRThybrid-flash would be adequate with positional uncertainty<5 mm.

A spot-specific range uncertainty framework for robust optimization of proton therapy treatments

PURPOSE

An accurate estimation of range uncertainties is essential to exploit the potential of proton therapy. According to Paganetti's study, a value of 2.4% (1.5 standard deviation) is currently recommended for planning robust treatments with Monte Carlo dose engines. This number is based on a dominant contribution from the mean excitation energy of tissues. However, it was recently shown that expressing tissues as a mixture of water and "dry" material in the CT calibration process allowed for a significant reduction of this uncertainty. We thus propose an adapted framework for pencil beam scanning robust optimization. First, we move towards a spot-specific range uncertainty (SSRU) determination. Second, we use the water-based formalism to reduce range uncertainties and, potentially, to spare better the organs at risk.

METHODS

The stoichiometric calibration was adapted to provide a molecular decomposition (including water) of each voxel of the CT. The SSRU calculation was implemented in MCsquare, a fast Monte Carlo dose engine dedicated to proton therapy. For each spot, a ray-tracing method was used to propagate molecular I-values uncertainties and obtain the corresponding effective range uncertainty. These were then combined with other sources of range uncertainties, according to Paganetti's study of 2012. The method was then assessed on three head-and-neck patients. Two plans were optimized for each patient: the first one with the classical 2.4% flat range uncertainty (FRU), the second one with the variable range uncertainty. Both plans were then compared in terms of target coverage and OAR mean dose reduction. Robustness evaluations were also performed, using the SSRU for both plans in order to simulate errors as realistically as possible.

RESULTS

For patient 1, it was found that the median SSRU was 1.04% (1.5 standard deviation), yielding, therefore, a very large reduction from the 2.4% FRU. All three SSRU plans were found to have a very good robustness level at a 90% confidence interval while sparing OAR better than the classical plan. For instance, in nominal cases, average reductions in the mean dose of 15.7, 8.4, and 13.2% were observed in the left parotid, right parotid, and pharyngeal constrictor muscle, respectively. As expected, the classical plans showed a higher but unnecessary level of robustness.

CONCLUSIONS

Promising results of the SSRU framework were observed on three head-and-neck cases, and more patients should now be considered. The method could also benefit to other tumor sites and, in the long run, the variable part of the range uncertainty could be generalized to other sources of uncertainty in order to move towards more and more patient-specific treatments.

Validation of robust radiobiological optimization algorithms based on the mixed beam model for intensity-modulated carbon-ion therapy

Currently, treatment planning systems (TPSs) that can compute the intensities of intensity-modulated carbon-ion therapy (IMCT) using scanned carbon-ion beams are limited. In the present study, the computational efficacy of the newly designed IMCT algorithms was analyzed for the first time based on the mixed beam model with respect to the physical and biological doses; moreover, the validity and effectiveness of the robust radiobiological optimization were verified. A dose calculation engine was independently generated to validate a clinical dose determined in the TPS. A biological assay was performed using the HSGc-C5 cell line to validate the calculated surviving fraction (SF). Both spot control (SC) and voxel-wise worst-case scenario (WC) algorithms were employed for robust radiobiological optimization followed by their application in a Radiation Therapy Oncology Group benchmark phantom under homogeneous and heterogeneous conditions and a clinical case for range and position errors. Importantly, for the first time, both SC and WC algorithms were implemented in the integrated TPS platform that can compute the intensities of IMCT using scanned carbon-ion beams for robust radiobiological optimization. For assessing the robustness, the difference between the maximum and minimum values of a dose-volume histogram index in the examined error scenarios was considered as a robustness index. The relative biological effectiveness (RBE) determined by the independent dose calculation engine exhibited a -0.6% difference compared with the RBE defined by the TPS at the isocenter, whereas the measured and the calculated SF were similar. Regardless of the objects, compared with the conventional IMCT, the robust radiobiological optimization enhanced the sensitivity of the examined error scenarios by up to 19% for the robustness index. The computational efficacy of the novel IMCT algorithms was verified according to the mixed beam model with respect to the physical and biological doses. The robust radiobiological optimizations lowered the impact of range and position uncertainties considerably in the examined scenarios. The robustness of the WC algorithm was more enhanced compared with that of the SC algorithm. Nevertheless, the SC algorithm can be used as an alternative to the WC IMCT algorithm with respect to the computational cost.

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.

Online adaptive planning methods for intensity-modulated radiotherapy

Online adaptive radiation therapy aims at adapting a patient's treatment plan to their current anatomy to account for inter-fraction variations before daily treatment delivery. As this process needs to be accomplished while the patient is immobilized on the treatment couch, it requires time-efficient adaptive planning methods to generate a quality daily treatment plan rapidly. The conventional planning methods do not meet the time requirement of online adaptive radiation therapy because they often involve excessive human intervention, significantly prolonging the planning phase. This article reviews the planning strategies employed by current commercial online adaptive radiation therapy systems, research on online adaptive planning, and artificial intelligence's potential application to online adaptive planning.

Deep-Learning-based Fast and Accurate 3D CT Deformable Image Registration in Lung Cancer

PURPOSE

In some proton therapy facilities, patient alignment relies on two 2D orthogonal kV images, taken at fixed, oblique angles, as no 3D on-the-bed imaging is available. The visibility of the tumor in kV images is limited since the patient's 3D anatomy is projected onto a 2D plane, especially when the tumor is behind high-density structures such as bones. This can lead to large patient setup errors. A solution is to reconstruct the 3D CT image from the kV images obtained at the treatment isocenter in the treatment position.

METHODS

An asymmetric autoencoder-like network built with vision-transformer blocks was developed. The data was collected from 1 head and neck patient: 2 orthogonal kV images (1024x1024 voxels), 1 3D CT with padding (512x512x512) acquired from the in-room CT-on-rails before kVs were taken and 2 digitally-reconstructed-radiograph (DRR) images (512x512) based on the CT. We resampled kV images every 8 voxels and DRR and CT every 4 voxels, thus formed a dataset consisting of 262,144 samples, in which the images have a dimension of 128 for each direction. In training, both kV and DRR images were utilized, and the encoder was encouraged to learn the jointed feature map from both kV and DRR images. In testing, only independent kV images were used. The full-size synthetic CT (sCT) was achieved by concatenating the sCTs generated by the model according to their spatial information. The image quality of the synthetic CT (sCT) was evaluated using mean absolute error (MAE) and per-voxel-absolute-CT-number-difference volume histogram (CDVH).

RESULTS

The model achieved a speed of 2.1s and a MAE of <40HU. The CDVH showed that <5% of the voxels had a per-voxel-absolute-CT-number-difference larger than 185 HU.

CONCLUSION

A patient-specific vision-transformer-based network was developed and shown to be accurate and efficient to reconstruct 3D CT images from kV images.

Fraction optimization for hybrid proton-photon treatment planning

BACKGROUND

Hybrid proton-photon radiotherapy (RT) can provide better plan quality than proton or photon only RT, in terms of robustness of target coverage and sparing of organs-at-risk (OAR).

PURPOSE

This work develops a hybrid treatment planning method that can optimize the number of proton and photon fractions as well as proton and photon plan variables, so that the hybrid plans can be clinically delivered day-to-day using either proton or photon machine.

METHODS

In the new hybrid treatment planning method, the total dose distribution (sum of proton dose and photon dose) is optimized for robust target coverage and optimal OAR sparing, by jointly optimizing proton spots and photon fluences, while the target dose uniformity is also enforced individually on both proton dose and photon dose, so that the hybrid plans can be separately and robustly delivered on proton or photon machine. To ensure the target dose uniformity for proton and photon plans, the number of proton and photon fractions is explicitly modeled and optimized, so that the target dose for proton and photon dose components is constrained to be a constant fraction of the total prescription dose while the plan quality based on total dose is optimized. The feasibility of hybrid planning using the proposed method is validated with representative clinical cases including abdomen, lung, head-and-neck (HN), and brain.

RESULTS

For all cases, hybrid plans provided better overall plan quality and OAR sparing than proton-only or photon-only plans, better target dose uniformity and robustness than proton-only plans, quantified by treatment planning objectives and dosimetric parameters. Moreover, for HN and brain cases, hybrid plans also had better target coverage than photon-only plans.

CONCLUSIONS

We have developed a new hybrid treatment planning method that optimizes number of proton and photon fractions as well as proton spots and photon fluences, for generating hybrid plans that can be separately and robustly delivered on proton or photon machines. Preliminary results have demonstrated that hybrid plans generated by the new method have better plan quality than proton-only or photon-only plans.

Automatic detection and classification of treatment deviations in proton therapy from realistically simulated prompt gamma imaging data

BACKGROUND

A clinical study regarding the potential of range verification in proton therapy (PT) by prompt gamma imaging (PGI) is carried out at our institution. Manual interpretation of the detected spot-wise range shift information is time-consuming, highly complex, and therefore not feasible in a broad routine application.

PURPOSE

Here, we present an approach to automatically detect and classify treatment deviations in realistically simulated PGI data for head-and-neck cancer (HNC) treatments using convolutional neural networks (CNNs) and conventional machine learning (ML) approaches.

METHODS

For 12 HNC patients and 1 anthropomorphic head phantom (n = 13), pencil beam scanning (PBS) treatment plans were generated, and 1 field per plan was assumed to be monitored with a PGI slit camera system. In total, 386 scenarios resembling different relevant or non-relevant treatment deviations were simulated on planning and control CTs and manually classified into 7 classes: non-relevant changes (NR) and relevant changes (RE) triggering treatment intervention due to range prediction errors (±RP), setup errors in beam direction (±SE), anatomical changes (AC), or a combination of such errors (CB). PBS spots with reliable PGI information were considered with their nominal Bragg peak position for the generation of two 3D spatial maps of 16 × 16 × 16 voxels containing PGI-determined range shift and proton number information. Three complexity levels of simulated PGI data were investigated: (I) optimal PGI data, (II) realistic PGI data with simulated Poisson noise based on the locally delivered proton number, and (III) realistic PGI data with an additional positioning uncertainty of the slit camera following an experimentally determined distribution. For each complexity level, 3D-CNNs were trained on a data subset (n = 9) using patient-wise leave-one-out cross-validation and tested on an independent test cohort (n = 4). Both the binary task of detecting RE and the multi-class task of classifying the underlying error source were investigated. Similarly, four different conventional ML classifiers (logistic regression, multilayer perceptron, random forest, and support vector machine) were trained using five previously established handcrafted features extracted from the PGI data and used for performance comparison.

RESULTS

On the test data, the CNN ensemble achieved a binary accuracy of 0.95, 0.96, and 0.93 and a multi-class accuracy of 0.83, 0.81, and 0.76 for the complexity levels (I), (II), and (III), respectively. In the case of binary classification, the CNN ensemble detected treatment deviations in the most realistic scenario with a sensitivity of 0.95 and a specificity of 0.88. The best performing ML classifiers showed a similar test performance.

CONCLUSIONS

This study demonstrates that CNNs can reliably detect relevant changes in realistically simulated PGI data and classify most of the underlying sources of treatment deviations. The CNNs extracted meaningful features from the PGI data with a performance comparable to ML classifiers trained on previously established handcrafted features. These results highlight the potential of a reliable, automatic interpretation of PGI data for treatment verification, which is highly desired for a broad clinical application and a prerequisite for the inclusion of PGI in an automated feedback loop for online adaptive PT.

Potential benefits of using radioactive ion beams for range margin reduction in carbon ion therapy

Sharp dose gradients and high biological effectiveness make ions such as ¹²C an ideal tool to treat deep-seated tumors, however, at the same time, sensitive to errors in the range prediction. Tumor safety margins mitigate these uncertainties, but during the irradiation they lead to unavoidable damage to the surrounding healthy tissue. To fully exploit the Bragg peak benefits, a large effort is put into establishing precise range verification methods. Despite positron emission tomography being widely in use for this purpose in ¹²C therapy, the low count rates, biological washout, and broad activity distribution still limit its precision. Instead, radioactive beams used directly for treatment would yield an improved signal and a closer match with the dose fall-off, potentially enabling precise in vivo beam range monitoring. We have performed a treatment planning study to estimate the possible impact of the reduced range uncertainties, enabled by radioactive ¹¹C ions treatments, on sparing critical organs in tumor proximity. Compared to ¹²C treatments, (i) annihilation maps for ¹¹C ions can reflect sub- millimeter shifts in dose distributions in the patient, (ii) outcomes of treatment planning with ¹¹C significantly improve and (iii) less severe toxicities for serial and parallel critical organs can be expected.

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.

An online adaptive plan library approach for intensity modulated proton therapy for head and neck cancer

BACKGROUND AND PURPOSE

In intensity modulated proton therapy (IMPT), the impact of setup errors and anatomical changes is commonly mitigated by robust optimization with population-based setup robustness (SR) settings and offline replanning. In this study we propose and evaluate an alternative approach based on daily plan selection from patient-specific pre-treatment established plan libraries (PLs). Clinical implementation of the PL strategy would be rather straightforward compared to daily online re-planning.

MATERIALS AND METHODS

For 15 head-and-neck cancer patients, the planning CT was used to generate a PL with 5 plans, robustly optimized for increasing SR: 0, 1, 2, 3, 5 mm, and 3% range robustness. Repeat CTs (rCTs) and realistic setup and range uncertainty distributions were used for simulation of treatment courses for the PL approach, treatments with fixed SR (fSR₃) and a trigger-based offline adaptive schedule for 3 mm SR (fSR₃OfA). Daily plan selection in the PL approach was based only on recomputed dose to the CTV on the rCT.

RESULTS

Compared to using fSR₃ and fSR₃OfA, the risk of xerostomia grade ≥ II & III and dysphagia ≥ grade III were significantly reduced with the PL. For 6/15 patients the risk of xerostomia and/or dysphagia ≥ grade II could be reduced by > 2% by using PL. For the other patients, adherence to target coverage constraints was often improved. fSR₃OfA resulted in significantly improved coverage compared to PL for selected patients.

CONCLUSION

The proposed PL approach resulted in overall reduced NTCPs compared to fSR₃ and fSR₃OfA at limited cost in target coverage.

Virtual particle Monte Carlo: A new concept to avoid simulating secondary particles in proton therapy dose calculation

BACKGROUND

In proton therapy dose calculation, Monte Carlo (MC) simulations are superior in accuracy but more time consuming, compared to analytical calculations. Graphic processing units (GPUs) are effective in accelerating MC simulations but may suffer thread divergence and racing condition in GPU threads that degrades the computing performance due to the generation of secondary particles during nuclear reactions.

PURPOSE

A novel concept of virtual particle (VP) MC (VPMC) is proposed to avoid simulating secondary particles in GPU-accelerated proton MC dose calculation and take full advantage of the computing power of GPU.

METHODS

Neutrons and gamma rays were ignored as escaping from the human body; doses of electrons, heavy ions, and nuclear fragments were locally deposited; the tracks of deuterons were converted into tracks of protons. These particles, together with primary and secondary protons, are considered to be the realistic particles. Histories of primary and secondary protons were replaced by histories of multiple VPs. Each VP corresponded to one proton (either primary or secondary). A continuous-slowing-down-approximation model, an ionization model, and a large angle scattering event model corresponding to nuclear interactions were developed for VPs by generating probability distribution functions (PDFs) based on simulation results of realistic particles using MCsquare. For efficient calculations, these PDFs were stored in the Compute Unified Device Architecture textures. VPMC was benchmarked with TOPAS and MCsquare in phantoms and with MCsquare in 13 representative patient geometries. Comparisons between the VPMC calculated dose and dose measured in water during patient-specific quality assurance (PSQA) of the selected 13 patients were also carried out. Gamma analysis was used to compare the doses derived from different methods and calculation efficiencies were also compared.

RESULTS

Integrated depth dose and lateral dose profiles in both homogeneous and inhomogeneous phantoms all matched well among VPMC, TOPAS, and MCsquare calculations. The 3D-3D gamma passing rates with a criterion of 2%/2 mm and a threshold of 10% was 98.49% between MCsquare and TOPAS and 98.31% between VPMC and TOPAS in homogeneous phantoms, and 99.18% between MCsquare and TOPAS and 98.49% between VPMC and TOPAS in inhomogeneous phantoms, respectively. In patient geometries, the 3D-3D gamma passing rates with 2%/2 mm/10% between dose distributions from VPMC and MCsquare were 98.56 ± 1.09% in patient geometries. The 2D-3D gamma analysis with 3%/2 mm/10% between the VPMC calculated dose distributions and the 2D measured planar dose distributions during PSQA was 98.91 ± 0.88%. VPMC calculation was highly efficient and took 2.84 ± 2.44 s to finish for the selected 13 patients running on four NVIDIA Ampere GPUs in patient geometries.

CONCLUSION

VPMC was found to achieve high accuracy and efficiency in proton therapy dose calculation.

Assessment of IMPT versus VMAT plans using different uncertainty scenarios for prostate cancer

Background

To assess the impact of systematic setup and range uncertainties for robustly optimized (RO) intensity modulated proton therapy (IMPT) and volumetric modulated arc therapy (VMAT) plans in patients with localized prostate cancer.

Methods

Twenty-six localized prostate patients previously treated with VMAT (CTV to PTV expansion of 3-5 mm) were re-planned with RO-IMPT with 3 mm and 5 mm geometrical uncertainties coupled with 3% range uncertainties. Robust evaluations (RE) accounting for the geometrical uncertainties of 3 and 5 mm were evaluated for the IMPT and VMAT plans. Clinical target volume (CTV), anorectum, and bladder dose metrics were analyzed between the nominal plans and their uncertainty perturbations.

Results

With geometric uncertainties of 5 mm and accounting for potential inter-fractional perturbations, RO-IMPT provided statistically significant (p < 0.05) sparing at intermediate doses (V_4000cGy) to the anorectum and bladder and high dose sparring (V_8000cGy) to the bladder compared to VMAT. Decreasing the RO and RE parameters to 3 mm improved IMPT sparing over VMAT at all OAR dose levels investigated while maintaining equivalent coverage to the CTV.

Conclusions

For localized prostate treatments, if geometric uncertainties can be maintained at or below 3 mm, RO-IMPT provides clear dosimetric advantages in anorectum and bladder sparing compared to VMAT. This advantage remains even under uncertainty scenarios. As geometric uncertainties increase to 5 mm, RO-IMPT still provides dosimetric advantages, but to a smaller magnitude.

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.

Range uncertainty reductions in proton therapy may lead to the feasibility of novel beam arrangements which improve organ-at-risk sparing

PURPOSE

In proton therapy, dose distributions are currently often conformed to organs at risk (OARs) using the less sharp dose fall-off at the lateral beam edge to reduce the effects of uncertainties in the in vivo proton range. However, range uncertainty reductions may make greater use of the sharper dose fall-off at the distal beam edge feasible, potentially improving OAR sparing. We quantified the benefits of such novel beam arrangements.

METHODS

For each of 10 brain or skull base cases, five treatment plans robust to 2 mm setup and 0%-4% range uncertainty were created for the traditional clinical beam arrangement and a novel beam arrangement making greater use of the distal beam edge to conform the dose distribution to the brainstem. Metrics including the brainstem normal tissue complication probability (NTCP) with the endpoint of necrosis were determined for all plans and all setup and range uncertainty scenarios.

RESULTS

For the traditional beam arrangement, reducing the range uncertainty from the current level of approximately 4% to a potentially achievable level of 1% reduced the brainstem NTCP by up to 0.9 percentage points in the nominal and up to 1.5 percentage points in the worst-case scenario. Switching to the novel beam arrangement at 1% range uncertainty improved these values by a factor of 2, that is, to 1.8 percentage points and 3.2 percentage points, respectively. The novel beam arrangement achieved a lower brainstem NTCP in all cases starting at a range uncertainty of 2%.

CONCLUSION

The benefits of novel beam arrangements may be of the same magnitude or even exceed the direct benefits of range uncertainty reductions. Indirect effects may therefore contribute markedly to the benefits of reducing proton range uncertainties.

Redefine the Role of Proton Beam Therapy for the Locally-Advanced Non-Small Cell Lung Cancer Assisting the Reduction of Acute Hematologic Toxicity

Purpose

To investigate the potential clinical benefit of utilizing intensity-modulated proton therapy (IMPT) to reduce acute hematologic toxicity for locally advanced non-small cell lung cancer (LA-NSCLC) patients and explore the feasibility of a model-based patient selection approach via the normal tissue complication probability (NTCP).

Methods

Twenty patients with LA-NSCLC were retrospectively selected. Volumetric modulated arc photon therapy (VMAT) and IMPT plans were generated with a prescription dose of 60 Gy in 30 fractions. A wide range of cases with varied tumor size, location, stations of metastatic lymph nodes were selected to represent the general cancer group. Contouring and treatment planning followed RTOG-1308 protocol. Doses to thoracic vertebral bodies (TVB) and other organ at risks were compared. Risk of grade ≥ 3 acute hematologic toxicity (HT3+) were calculated based on the NTCP model, and patients with a reduction on NTCP of HT3+ from VMAT to IMPT (△NTCP_HT3+) ≥ 10% were considered to ‘significantly benefit from proton therapy.’

Results

Compared to VMAT, IMPT significantly reduced the dose to the TVB, the lung, the heart, the esophagus and the spinal cord. Tumor distance to TVB was significantly associated with △NTCP _HT3+ ≥ 10%. For the patients with tumor distance ≤ 0.7 cm to TVB, the absolute reduction of dose (mean, V30 and V40) to TVB was significantly lower than that in patients with tumor distance &gt; 0.7 cm.

Conclusion

IMPT decreased the probability of HT3+ compared to VMAT by reducing the dose to the TVB in LA-NSCLC patients. Patients with tumor distance to TVB less than 0.7 cm are likely to benefit most from proton over photon therapy.

Proton Therapy for Squamous Cell Carcinoma of the Head and Neck: Early Clinical Experience and Current Challenges

Simple Summary

Proton therapy is a promising type of radiation therapy used to destroy tumor cells. It has the potential to further improve the outcomes for patients with head and neck cancer since it allows to minimize the radiation dose to vital structures around the tumor, leading to less toxicity. This paper describes the current experience worldwide with proton therapy in head and neck cancer.

Abstract

Proton therapy (PT) is a promising development in radiation oncology, with the potential to further improve outcomes for patients with squamous cell carcinoma of the head and neck (HNSCC). By utilizing the finite range of protons, healthy tissue can be spared from beam exit doses that would otherwise be irradiated with photon-based treatments. Current evidence on PT for HNSCC is limited to comparative dosimetric analyses and retrospective single-institution series. As a consequence, the recognized indications for the reimbursement of PT remain scarce in most countries. Nevertheless, approximately 100 PT centers are in operation worldwide, and initial experiences for HNSCC are being reported. This review aims to summarize the results of the early clinical experience with PT for HNSCC and the challenges that are currently faced.

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.

Evaluation of the clinical value of automatic online dose restoration for adaptive proton therapy of head and neck cancer

INTRODUCTION

Intensity modulated proton therapy (IMPT) is highly sensitive to anatomical variations which can cause inadequate target coverage during treatment. This study compares not-adapted (NA) robust plans to two adaptive IMPT methods - a fully-offline adaptive (FOA) and a simplified automatic online adaptive strategy (dose restoration (DR)) to determine the benefit of DR, in head and neck cancer (HNC).

MATERIAL/METHODS

Robustly optimized clinical IMPT doses in planning-CTs (pCTs) were available for a cohort of 10 HNC patients. During robust re-optimization, DR used isodose contours, generated from the clinical dose on pCTs, and patient specific objectives to reproduce the clinical dose in every repeated-CT(rCT). For each rCT(n=50), NA, DR and FOA plans were robustly evaluated.

RESULTS

An improvement in DVH-metrics and robustness was seen for DR and FOA plans compared to NA plans. For NA plans, 74%(37/50) of rCTs did not fulfill the CTV coverage criteria (D98%>95%Dprescription). DR improved target coverage, target homogeneity and variability on critical risk organs such as the spinal cord. After DR, 52%(26/50) of rCTs met all clinical goals. Because of large anatomical changes and/or inaccurate patient repositioning, 48%(24/50) of rCTs still needed full offline adaptation to ensure an optimal treatment since dose restoration was not able to re-establish the initial plan quality.

CONCLUSION

Robust optimization together with fully-automatized DR avoided offline adaptation in 52% of the cases. Implementation of dose restoration in clinical routine could ensure treatment plan optimality while saving valuable human and material resources to radiotherapy departments.

Validation of automated complex head and neck treatment planning with pencil beam scanning proton therapy

Background

Pencil beam scanning (PBS) proton therapy offers dosimetric advantages for several treatment sites, including head and neck (H&N). However, to achieve the optimal target coverage and robustness, these plans can be complex and time consuming to develop and optimize. Automating the treatment planning process can ensure a high‐quality and standardized plan, reduce burden to the planner, and decrease time‐to‐treatment. We utilized in‐house scripting to automate a four‐field multi‐field optimization (MFO) H&N planning technique.

Methods and materials

Ten bilateral H&N patients were planned in RayStation v6 with a four‐field modified‐X beam configuration using MFO planning. Automation included creation of avoidance structures to control spot placement and development of standardized beams, PBS spot settings, robust optimization objectives, and patient‐specific predicted planning constraints. Each patient was planned both with and without automation to evaluate differences in planning time, perceived effort and plan quality, plan robustness, and OAR sparing.

Results

On average, scripted plans required 3.2 h, compared to 4.3 h without the script. There was no difference in target coverage or plan robustness with or without automation. Automation significantly reduced mean dose to the oral cavity, parotids, esophagus, trachea, and larynx. Perceived effort was scaled from 1 (minimum effort) to 100 (maximum effort), and automation reduced perceived effort by 42% (p < 0.05). Two non‐scripted plans required re‐planning due to errors.

Conclusions

Automation of this multi‐beam, the MFO proton planning process reduced planning time and improved OAR sparing compared to the same planning process without scripting. Scripting generation of complex structures and planning objectives reduced burden on the planner. With most current treatment planning software, this automation is simple to implement and can standardize quality of care across all treatment planners.

A comparison of proton stopping power measured with proton CT and x-ray CT in fresh postmortem porcine structures

PURPOSE

Currently, calculations of proton range in proton therapy patients are based on a conversion of CT Hounsfield units of patient tissues into proton relative stopping power. Uncertainties in this conversion necessitate larger proximal and distal planned target volume margins. Proton CT can potentially reduce these uncertainties by directly measuring proton stopping power. We aim to demonstrate proton CT imaging with complex porcine samples, to analyze in detail three-dimensional regions of interest, and to compare proton stopping powers directly measured by proton CT to those determined from x-ray CT scans.

METHODS

We have used a prototype proton imaging system with single proton tracking to acquire proton radiography and proton CT images of a sample of porcine pectoral girdle and ribs, and a pig's head. We also acquired close in time x-ray CT scans of the same samples and compared proton stopping power measurements from the two modalities. In the case of the pig's head, we obtained x-ray CT scans from two different scanners and compared results from high-dose and low-dose settings.

RESULTS

Comparing our reconstructed proton CT images with images derived from x-ray CT scans, we find agreement within 1% to 2% for soft tissues and discrepancies of up to 6% for compact bone. We also observed large discrepancies, up to 40%, for cavitated regions with mixed content of air, soft tissue, and bone, such as sinus cavities or tympanic bullae.

CONCLUSIONS

Our images and findings from a clinically realistic proton CT scanner demonstrate the potential for proton CT to be used for low-dose treatment planning with reduced margins.

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.

Probabilistic target definition and planning in patients with prostate cancer

Intro.Current radiation therapy (RT) planning guidelines handle uncertainties in RT using geometric margins. This approach is simple to use but oversimplifies complex underlying processes and is cumbersome for non-homogeneous dose prescriptions. In this work, we characterize the performance of a novel probabilistic target definition and planning (PTP) approach, which uses voxel-level tumor likelihood information in treatment plan optimization.Methods.We expanded a treatment planning system with probabilistic therapy planning functionality that utilizes non-binary target maps (TM) as voxel-level input to dose plan optimization. Different dose plans were calculated and compared for twelve prostate cancer patients with multiparametric magnetic resonance imaging derived TMs. Dose plans were created using both classical and PTP approaches for uniform and integrated dose boost prescriptions. Dose performance between the different approaches was compared using dose benchmarks on target and organ-at-risk (OAR) volumes.Results.Over all dose metrics, PTP was shown to be comparable to classical planning. For plans of uniform dose prescription, the PTP approach created plans within 1 Gy of the classical planning approach across all dose metrics, with no significant differences (p > 0.2). For plans with the integrated dose boost, PTP plans exhibited higher dose heterogeneity, but still showed target doses comparable to the classical approach, without increasing doses to OAR.Conclusion.In this work we introduce direct incorporation of probabilistic target definition into treatment planning. This treatment planning approach can produce both uniform dose plans and plans with integrated dose boosts that are comparable to ones created using classical dose planning. PTP is a flexible way to optimize external beam radiotherapy, as it is not limited by the use of margins. PTP can produce dose plans equivalent to classical planning, while also allows for greater versatility in dose prescription and direct incorporation of patient target definition uncertainty into treatment planning.

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.

Proton range uncertainty reduction benefits for skull base tumors in terms of normal tissue complication probability (NTCP) and healthy tissue doses

PURPOSE

Proton therapy allows for more conformal dose distributions and lower organ at risk and healthy tissue doses than conventional photon-based radiotherapy, but uncertainties in the proton range currently prevent proton therapy from making full use of these advantages. Numerous developments therefore aim to reduce such range uncertainties. In this work, we quantify the benefits of reductions in range uncertainty for treatments of skull base tumors.

METHODS

The study encompassed 10 skull base patients with clival tumors. For every patient, six treatment plans robust to setup errors of 2 mm and range errors from 0% to 5% were created. The determined metrics included the brainstem and optic chiasm normal tissue complication probability (NTCP) with the endpoints of necrosis and blindness, respectively, as well as the healthy tissue volume receiving at least 70% of the prescription dose.

RESULTS

A range uncertainty reduction from the current level of 4% to a potentially achievable level of 1% reduced the probability of brainstem necrosis by up to 1.3 percentage points in the nominal scenario in which neither setup nor range errors occur and by up to 2.9 percentage points in the worst-case scenario. Such a range uncertainty reduction also reduced the optic chiasm NTCP with the endpoint of blindness by up to 0.9 percentage points in the nominal scenario and by up to 2.2 percentage points in the worst-case scenario. The decrease in the healthy tissue volume receiving at least 70% of the prescription dose ranged from -7.8 to 24.1 cc in the nominal scenario and from -3.4 to 38.4 cc in the worst-case scenario.

CONCLUSION

The benefits quantified as part of this study serve as a guideline of the OAR and healthy tissue dose benefits that range monitoring techniques may be able to achieve. Benefits were observed between all levels of range uncertainty. Even smaller range uncertainty reductions may therefore be beneficial.