Online adaptive dose restoration in intensity modulated proton therapy of lung cancer to account for inter-fractional density changes

Robustness evaluation of pencil beam scanning proton therapy treatment planning: A systematic review

Compared to conventional radiotherapy using X-rays, proton therapy, in principle, allows better conformity of the dose distribution to target volumes, at the cost of greater sensitivity to physical, anatomical, and positioning uncertainties. Robust planning, both in terms of plan optimization and evaluation, has gained high visibility in publications on the subject and is part of clinical practice in many centers. However, there is currently no consensus on the methods and parameters to be used for robust optimization or robustness evaluation. We propose to overcome this deficiency by following the modified Delphi consensus method. This method first requires a systematic review of the literature. We performed this review using the PubMed and Web Of Science databases, via two different experts. Potential conflicts were resolved by a third expert. We then explored the different methods before focusing on clinical studies that evaluate robustness on a significant number of patients. Many robustness assessment methods are proposed in the literature. Some are more successful than others and their implementation varies between centers. Moreover, they are not all statistically or mathematically equivalent. The most sophisticated and rigorous methods have seen more limited application due to the difficulty of their implementation and their lack of widespread availability.

Time is NTCP: Should we maximize patient throughput or perform online adaptation on proton therapy systems?

BACKGROUND

Compared to conventional radiotherapy (XT), proton therapy (PT) may improve normal tissue complication probabilities (NTCP). However, PT typically requires higher adaptation rates due to an increased sensitivity to anatomical changes. Systematic online adaptation may address this issue, but it requires additional replanning time, decreasing patient throughput. Therefore, less patients would benefit in such case from PT for a given machine capacity, with results in worse NTCP.

AIM

To investigate the trade-off between PT patient throughput and NTCP gain as a function of the time needed for adaptation.

METHODS

A retrospective database of 14 lung patients with two repeated 4DCTs was used to compare NTCP values between XT and PT for NTCP2ym (2-year mortality), NTCPdysphagia and NTCPpneumonitis. Four scenarios were considered for PT: no adaptation using clinical robustness parameters (4D robust optimization, 3 % range error and PTV-equivalent setup errors); systematic online adaptation with clinical robustness parameters; setup errors reduced to 4 mm and to 2 mm. Dose was accumulated on the planning CT. The number of patients treated with PT depended on the extra time needed for adaptation, assuming an 8-hours capacity (assuming 14 patients a day; thus minimum 34.2 min per treatment session if there is no or instantaneous adaptation).

RESULTS

Baseline NTCP gains (PT against XT without adaptation) equaled 6.9 %, 6.1 %, and 7.7 % for NTCP2ym, NTCPdysphagia and NTCPpneumonitis, respectively. Using instantaneous online adaptation and setup errors of 2 mm, the overall gains were then 10.7 %, 13.6 % and 12.4 %. Taking into account loss of capacity, 13.7 min was the maximum extra-time allowed to complete adaptation and maintain an advantage on all three metrics for the 2-mm setup error scenario.

CONCLUSION

This study highlights the critical importance of keeping short online adaptation times when using systems with limited capacity like PT.

A review of the clinical introduction of 4D particle therapy research concepts

Background and purpose

Many 4D particle therapy research concepts have been recently translated into clinics, however, remaining substantial differences depend on the indication and institute-related aspects. This work aims to summarise current state-of-the-art 4D particle therapy technology and outline a roadmap for future research and developments.

Material and methods

This review focused on the clinical implementation of 4D approaches for imaging, treatment planning, delivery and evaluation based on the 2021 and 2022 4D Treatment Workshops for Particle Therapy as well as a review of the most recent surveys, guidelines and scientific papers dedicated to this topic.

Results

Available technological capabilities for motion surveillance and compensation determined the course of each 4D particle treatment. 4D motion management, delivery techniques and strategies including imaging were diverse and depended on many factors. These included aspects of motion amplitude, tumour location, as well as accelerator technology driving the necessity of centre-specific dosimetric validation. Novel methodologies for X-ray based image processing and MRI for real-time tumour tracking and motion management were shown to have a large potential for online and offline adaptation schemes compensating for potential anatomical changes over the treatment course. The latest research developments were dominated by particle imaging, artificial intelligence methods and FLASH adding another level of complexity but also opportunities in the context of 4D treatments.

Conclusion

This review showed that the rapid technological advances in radiation oncology together with the available intrafractional motion management and adaptive strategies paved the way towards clinical implementation.

Comparison of Ethos template-based planning and AI-based dose prediction: General performance, patient optimality, and limitations

PURPOSE

Ethos proposes a template-based automatic dose planning (Etb) for online adaptive radiotherapy. This study evaluates the general performance of Etb for prostate cancer, as well as the ability to generate patient-optimal plans, by comparing it with another state-of-the-art automatic planning method, i.e., deep learning dose prediction followed by dose mimicking (DP + DM).

MATERIALS

General performances and capability to produce patient-optimal plan were investigated through two studies: Study-S1 generated plans for 45 patients using our initial Ethos clinical goals template (EG_init), and compared them to manually generated plans (MG). For study-S2, 10 patients which showed poor performances at study-S1 were selected. S2 compared the quality of plans generated with four different methods: 1) Ethos initial template (EG_init_selected), 2) Ethos updated template-based on S1 results (EG_upd_selected), 3) DP + DM, and 4) MG plans.

RESULTS

EG_init plans showed satisfactory performance for dose level above 50 Gy: reported mean metrics differences (EG_init minus MG) never exceeded 0.6 %. However, lower dose levels showed loosely optimized metrics, mean differences for V30Gy to rectum and V20Gy to anal canal were of 6.6 % and 13.0 %. EG_init_selected showed amplified differences in V30Gy to rectum and V20Gy to anal canal: 8.5 % and 16.9 %, respectively. These dropped to 5.7 % and 11.5 % for EG_upd_selected plans but strongly increased V60Gy to rectum for 2 patients. DP + DM plans achieved differences of 3.4 % and 4.6 % without compromising any V60Gy.

CONCLUSION

General performances of Etb were satisfactory. However, optimizing with template of goals might be limiting for some complex cases. Over our test patients, DP + DM outperformed the Etb approach.

Patient selection for proton therapy using Normal Tissue Complication Probability with deep learning dose prediction for oropharyngeal cancer

BACKGROUND

In cancer care, determining the most beneficial treatment technique is a key decision affecting the patient's survival and quality of life. Patient selection for proton therapy (PT) over conventional radiotherapy (XT) currently entails comparing manually generated treatment plans, which requires time and expertise.

PURPOSE

We developed an automatic and fast tool, AI-PROTIPP (Artificial Intelligence Predictive Radiation Oncology Treatment Indication to Photons/Protons), that assesses quantitatively the benefits of each therapeutic option. Our method uses deep learning (DL) models to directly predict the dose distributions for a given patient for both XT and PT. By using models that estimate the Normal Tissue Complication Probability (NTCP), namely the likelihood of side effects to occur for a specific patient, AI-PROTIPP can propose a treatment selection quickly and automatically.

METHODS

A database of 60 patients presenting oropharyngeal cancer, obtained from the Cliniques Universitaires Saint Luc in Belgium, was used in this study. For every patient, a PT plan and an XT plan were generated. The dose distributions were used to train the two dose DL prediction models (one for each modality). The model is based on U-Net architecture, a type of convolutional neural network currently considered as the state of the art for dose prediction models. A NTCP protocol used in the Dutch model-based approach, including grades II and III xerostomia and grades II and III dysphagia, was later applied in order to perform automatic treatment selection for each patient. The networks were trained using a nested cross-validation approach with 11-folds. We set aside three patients in an outer set and each fold consists of 47 patients in training, five in validation and five for testing. This method allowed us to assess our method on 55 patients (five patients per test times the number of folds).

RESULTS

The treatment selection based on the DL-predicted doses reached an accuracy of 87.4% for the threshold parameters set by the Health Council of the Netherlands. The selected treatment is directly linked with these threshold parameters as they express the minimal gain brought by the PT treatment for a patient to be indicated to PT. To validate the performance of AI-PROTIPP in other conditions, we modulated these thresholds, and the accuracy was above 81% for all the considered cases. The difference in average cumulative NTCP per patient of predicted and clinical dose distributions is very similar (less than 1% difference).

CONCLUSIONS

AI-PROTIPP shows that using DL dose prediction in combination with NTCP models to select PT for patients is feasible and can help to save time by avoiding the generation of treatment plans only used for the comparison. Moreover, DL models are transferable, allowing, in the future, experience to be shared with centers that would not have PT planning expertise.

Scatter correction of 4D cone beam computed tomography to detect dosimetric effects due to anatomical changes in proton therapy for lung cancer

BACKGROUND

The treatment of moving tumor entities is expected to have superior clinical outcomes, using image-guided adaptive intensity-modulated proton therapy (IMPT).

PURPOSE

For 21 lung cancer patients, IMPT dose calculations were performed on scatter-corrected 4D cone beam CTs (4DCBCTcor ) to evaluate their potential for triggering treatment adaptation. Additional dose calculations were performed on corresponding planning 4DCTs and day-of-treatment 4D virtual CTs (4DvCTs).

METHODS

A 4DCBCT correction workflow, previously validated on a phantom, generates 4DvCT (CT-to-CBCT deformable registration) and 4DCBCTcor images (projection-based correction using 4DvCT as a prior) with 10 phase bins, using day-of-treatment free-breathing CBCT projections and planning 4DCT images as input. Using a research planning system, robust IMPT plans administering eight fractions of 7.5 Gy were created on a free-breathing planning CT (pCT) contoured by a physician. The internal target volume (ITV) was overridden with muscle tissue. Robustness settings for range and setup uncertainties were 3% and 6 mm, and a Monte Carlo dose engine was used. On every phase of planning 4DCT, day-of-treatment 4DvCT, and 4DCBCTcor , the dose was recalculated. For evaluation, image analysis as well as dose analysis were performed using mean error (ME) and mean absolute error (MAE) analysis, dose-volume histogram (DVH) parameters, and 2%/2-mm gamma pass rate analysis. Action levels (1.6% ITV D98 and 90% gamma pass rate) based on our previous phantom validation study were set to determine which patients had a loss of dosimetric coverage.

RESULTS

Quality enhancements of 4DvCT and 4DCBCTcor over 4DCBCT were observed. ITV D98% and bronchi D2% had its largest agreement for 4DCBCTcor -4DvCT, and the largest gamma pass rates (>94%, median 98%) were found for 4DCBCTcor -4DvCT. Deviations were larger and gamma pass rates were smaller for 4DvCT-4DCT and 4DCBCTcor -4DCT. For five patients, deviations were larger than the action levels, suggesting substantial anatomical changes between pCT and CBCT projections acquisition.

CONCLUSIONS

This retrospective study shows the feasibility of daily proton dose calculation on 4DCBCTcor for lung tumor patients. The applied method is of clinical interest as it generates up-to-date in-room images, accounting for breathing motion and anatomical changes. This information could be used to trigger replanning.

Comparison of 3D and 4D robustly optimized proton treatment plans for non-small cell lung cancer patients with tumour motion amplitudes larger than 5 mm

Background and purpose

There is no consensus about an ideal robust optimization (RO) strategy for proton therapy of targets with large intrafractional motion. We investigated the plan robustness of 3D and different 4D RO strategies.

Materials and methods

For eight non-small cell lung cancer patients with clinical target volume (CTV) motion >5 mm, different RO approaches were investigated: 3DRO considering the average CT (AvgCT) with a target density override, 4DRO considering three/all 4DCT phases, and 4DRO considering the AvgCT and three/all 4DCT phases. Robustness against setup/range errors, interplay effects based on breathing and machine log file data for deliveries with/without rescanning, and interfractional anatomical changes were analyzed for target coverage and OAR sparing.

Results

All nominal plans fulfilled the clinical requirements with individual CTV coverage differences <2pp; 4DRO without AvgCT generated the most conformal dose distributions. Robustness against setup/range errors was best for 4DRO with AvgCT (18% more passed error scenarios than 3DRO). Interplay effects caused fraction-wise median CTV coverage loss of 3pp and missed maximum dose constraints for heart and esophagus in 18% of scenarios. CTV coverage and OAR sparing fulfilled requirements in all cases when accumulating four interplay scenarios. Interfractional changes caused less target misses for RO with AvgCT compared to 4DRO without AvgCT (≤42%/33% vs. ≥56%/44% failed single/accumulated scenarios).

Conclusions

All RO strategies provided acceptable plans with equally low robustness against interplay effects demanding other mitigation than rescanning to ensure fraction-wise target coverage. 4DRO considering three phases and the AvgCT provided best compromise on planning effort and robustness.

Image guidance and interfractional anatomical variation in paediatric abdominal radiotherapy

OBJECTIVES

To identify variables predicting interfractional anatomical variations measured with cone-beam CT (CBCT) throughout abdominal paediatric radiotherapy, and to assess the potential of surface-guided radiotherapy (SGRT) to monitor these changes.

METHODS

Metrics of variation in gastrointestinal (GI) gas volume and separation of the body contour and abdominal wall were calculated from 21 planning CTs and 77 weekly CBCTs for 21 abdominal neuroblastoma patients (median 4 years, range: 2 - 19 years). Age, sex, feeding tubes, and general anaesthesia (GA) were explored as predictive variables for anatomical variation. Furthermore, GI gas variation was correlated with changes in body and abdominal wall separation, as well as simulated SGRT metrics of translational and rotational corrections between CT/CBCT.

RESULTS

GI gas volumes varied 74 ± 54 ml across all scans, while body and abdominal wall separation varied 2.0 ± 0.7 mm and 4.1 ± 1.5 mm from planning, respectively. Patients < 3.5 years (p = 0.04) and treated under GA (p < 0.01) experienced greater GI gas variation; GA was the strongest predictor in multivariate analysis (p < 0.01). Absence of feeding tubes was linked to greater body contour variation (p = 0.03). GI gas variation correlated with body (R = 0.53) and abdominal wall (R = 0.63) changes. The strongest correlations with SGRT metrics were found for anterior-posterior translation (R = 0.65) and rotation of the left-right axis (R = -0.36).

CONCLUSIONS

Young age, GA, and absence of feeding tubes were linked to stronger interfractional anatomical variation and are likely indicative of patients benefiting from adaptive/robust planning pathways. Our data suggest a role for SGRT to inform the need for CBCT at each treatment fraction in this patient group.

ADVANCES IN KNOWLEDGE

This is the first study to suggest the potential role of SGRT for the management of internal interfractional anatomical variation in paediatric abdominal radiotherapy.

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.

A survey of practice patterns for adaptive particle therapy for interfractional changes

Background and purpose

Anatomical changes may compromise the planned target coverage and organs-at-risk dose in particle therapy. This study reports on the practice patterns for adaptive particle therapy (APT) to evaluate current clinical practice and wishes and barriers to further implementation.

Materials and methods

An institutional questionnaire was distributed to PT centres worldwide (7/2020-6/2021) asking which type of APT was used, details of the workflow, and what the wishes and barriers to implementation were. Seventy centres from 17 countries participated. A three-round Delphi consensus analysis (10/2022) among the authors followed to define recommendations on required actions and future vision.

Results

Out of the 68 clinically operational centres, 84% were users of APT for at least one treatment site with head and neck being most common. APT was mostly performed offline with only two online APT users (plan-library). No centre used online daily re-planning. Daily 3D imaging was used for APT by 19% of users. Sixty-eight percent of users had plans to increase their use or change their technique for APT. The main barrier was "lack of integrated and efficient workflows". Automation and speed, reliable dose deformation for dose accumulation and higher quality of in-room volumetric imaging were identified as the most urgent task for clinical implementation of online daily APT.

Conclusion

Offline APT was implemented by the majority of PT centres. Joint efforts between industry research and clinics are needed to translate innovations into efficient and clinically feasible workflows for broad-scale implementation of online APT.

Deformed dose restoration to account for tumor deformation and position changes for adaptive proton therapy

BACKGROUND

Online adaptation during intensity-modulated proton therapy (IMPT) can minimize the effect of inter-fractional anatomical changes, but remains challenging because of the complex workflow. One approach for fast and automated online IMPT adaptation is dose restoration, which restores the initial dose distribution on the updated anatomy. However, this method may fail in cases where tumor deformation or position changes occur.

PURPOSE

To develop a fast and robust IMPT online adaptation method named "deformed dose restoration (DDR)" that can adjust for inter-fractional tumor deformation and position changes.

METHODS

The DDR method comprises two steps: (1) calculation of the deformed dose distribution, and (2) restoration of the deformed dose distribution. First, the deformable image registration (DIR) between the initial clinical target volume (CTV) and the new CTV were performed to calculate the vector field. To ensure robustness for setup and range uncertainty and the ability to restore the deformed dose distribution, an expanded CTV-based registration to maintain the dose gradient outside the CTV was developed. The deformed dose distribution was obtained by applying the vector field to the initial dose distribution. Then, the voxel-by-voxel dose difference optimization was performed to calculate beam parameters that restore the deformed dose distribution on the updated anatomy. The optimization function was the sum of total dose differences and dose differences of each field to restore the initial dose overlap of each field. This method only requires target contouring, which eliminates the need for organs at risk (OARs) contouring. Six clinical cases wherein the tumor deformation and/or position changed on repeated CTs were selected. DDR feasibility was evaluated by comparing the results with those from three other strategies, namely, not adapted (continuing the initial plan), adapted by previous dose restoration, and fully optimized.

RESULTS

In all cases, continuing the initial plan was largely distorted on the repeated CTs and the dose-volume histogram (DVH) metrics for the target were reduced due to the tumor deformation or position changes. On the other hand, DDR improved DVH metrics for the target to the same level as the initial dose distribution. Dose increase was seen for some OARs because tumor growth had reduced the relative distance between CTVs and OARs. Robustness evaluation for setup and range uncertainty (3 mm/3.5%) showed that deviation in DVH-bandwidth for CTV D_95% from the initial plan was 0.4% ± 0.5% (Mean ± S.D.) for DDR. The calculation time was 8.1 ± 6.4 min.

CONCLUSIONS

An online adaptation algorithm was developed that improved the treatment quality for inter-fractional anatomical changes and retained robustness for intra-fractional setup and range uncertainty. The main advantage of this method is that it only requires target contouring alone and saves the time for OARs contouring. The fast and robust adaptation method for tumor deformation and position changes described here can reduce the need for offline adaptation and improve treatment efficiency.

Low-dose CT allows for accurate proton therapy dose calculation and plan optimization

Objective. Protons offer a more conformal dose delivery compared to photons, yet they are sensitive to anatomical changes over the course of treatment. To minimize range uncertainties due to anatomical variations, a new CT acquisition at every treatment session would be paramount to enable daily dose calculation and subsequent plan adaptation. However, the series of CT scans results in an additional accumulated patient dose. Reducing CT radiation dose and thereby decreasing the potential risk of radiation exposure to patients is desirable, however, lowering the CT dose results in a lower signal-to-noise ratio and therefore in a reduced quality image. We hypothesized that the signal-to-noise ratio provided by conventional CT protocols is higher than needed for proton dose distribution estimation. In this study, we aim to investigate the effect of CT imaging dose reduction on proton therapy dose calculations and plan optimization. Approach. To verify our hypothesis, a CT dose reduction simulation tool has been developed and validated to simulate lower-dose CT scans from an existing standard-dose scan. The simulated lower-dose CTs were then used for proton dose calculation and plan optimization and the results were compared with those of the standard-dose scan. The same strategy was adopted to investigate the effect of CT dose reduction on water equivalent thickness (WET) calculation to quantify CT noise accumulation during integration along the beam. Main results. The similarity between the dose distributions acquired from the low-dose and standard-dose CTs was evaluated by the dose-volume histogram and the 3D Gamma analysis. The results on an anthropomorphic head phantom and three patient cases indicate that CT imaging dose reduction up to 90% does not have a significant effect on proton dose calculation and plan optimization. The relative error was employed to evaluate the similarity between WET maps and was found to be less than 1% after reducing the CT imaging dose by 90%. Significance. The results suggest the possibility of using low-dose CT for proton therapy dose estimation, since the dose distributions acquired from the standard-dose and low-dose CTs are clinically equivalent.

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.

Evaluation of alternative parameter settings for dose restoration and full plan adaptation in IMPT for prostate cancer

BACKGROUND/PURPOSE

Intensity-modulated proton therapy is highly sensitive to anatomical variations. A dose restoration method and a full plan adaptation method have been developed earlier, both requiring several parameter settings. This study evaluates the validity of the previously selected settings by systematically comparing them to alternatives.

MATERIALS/METHODS

The dose restoration method takes a prior plan and uses an energy-adaptation followed by a spot-intensity re-optimization to restore the plan to its initial state. The full adaptation method uses an energy-adaptation followed by the addition of new spots and a spot-intensity optimization to fit the new anatomy. We varied: 1) The margins and robustness settings of the prior plan, 2) the spot-addition sample size, i.e. the number of added spots, 3) the spot-addition stopping criterion, and 4) the spot-intensity optimization approach. The last three were evaluated only for the full plan adaptation. Evaluations were done on 88 CT scans of 11 prostate cancer patients. Dose was prescribed as 55 Gy(RBE) to the lymph nodes and seminal vesicles with a boost to 74 Gy(RBE) to the prostate.

RESULTS

For the dose restoration method, changing the applied CTV-to-PTV margins and plan robustness in the prior plans yielded insufficient target coverage or increased OAR doses. For the full plan adaptation, more spot-addition iterations and using a different optimization approach resulted in lower OAR doses compared to the default settings while maintaining target coverage. However, the calculation times increased by up to 20 times, making these variations infeasible for online-adaptation.

CONCLUSION

We recommend maintaining the default setting for the dose restoration approach. For the full plan adaptation we recommend to focus on fine-tuning the optimization-parameters, and apart from this using the default settings.