Time-resolved dosimetry for validation of 4D dose calculation in PBS proton therapy

Comparing the Robustness of Intensity-modulated Proton Therapy and Proton-arc Therapy Against Interplay Effects of 4D Robust-optimised Plans for Lung Stereotactic Body Radiotherapy

AIMS

To assess the robustness of 4D-optimised IMPT and PAT plans against interplay effects in non-small cell lung cancer (NSCLC) patients with respiratory motion over 10 mm, and to provide insights into the use of proton-based stereotactic body radiotherapy (SBRT) for lung cancer with significant tumour movement.

MATERIALS AND METHODS

Fourteen patients with early-stage NSCLC and tumour motion >10 mm were selected. Three hypofraction regimens were generated using 4D robust optimisation with the IMPT and PAT techniques. The nominal plan qualities for both techniques were compared, and their robustness against setup and range uncertainties was evaluated. 4D dynamic dose and the 4D static dose were generated to calculate ΔIMR(%) for interplay effects.

RESULTS

PAT plans demonstrated superior target metrics such as D95 and D2, and offered enhanced protection for organs at risk (OARs), particularly in lung metrics, across multiple fractionation schemes (p < 0.05). The robustness of target coverage against setup and range uncertainties was better in PAT plans than IMPT, with average pass rates of 97.8% and 95.4%, respectively (p < 0.01). The interplay effect significantly affected target metrics in single-fraction plans, decreasing with more fractions, while its effect on OAR metrics was minimal. Median values for single-fraction plans were: ΔID98GTV was -3% for IMPT and -0.7% for PAT (p < 0.01); ΔID95GTV was -2.4% for IMPT and -0.6% for PAT (p < 0.01); ΔID2GTV was 3.2% for IMPT and 0.9% for PAT (p < 0.05). The interplay effects resulted in median homogeneity index deviations of 9.1% and 2% for the IMPT and PAT plans, respectively (p < 0.01). Different starting phases affected IMPT more significantly than PAT.

CONCLUSION

PAT demonstrated greater robustness to interplay effects than IMPT for hypofractionated treatments of early-stage NSCLC, particularly in single-fraction schemes. Additionally, PAT showed good resilience to variations in different starting phases.

Quantifying dose uncertainties resulting from cardiorespiratory motion in intensity-modulated proton therapy for cardiac stereotactic body radiotherapy

Background

Cardiac stereotactic body radiotherapy (CSBRT) with photons efficaciously and safely treats cardiovascular arrhythmias. Proton therapy, with its unique physical and radiobiological properties, can offer advantages over traditional photon-based therapies in certain clinical scenarios, particularly pediatric tumors and those in anatomically challenging areas. However, dose uncertainties induced by cardiorespiratory motion are unknown.

Objective

This study investigated the effect of cardiorespiratory motion on intensity-modulated proton therapy (IMPT) and the effectiveness of motion-encompassing methods.

Methods

We retrospectively included 12 patients with refractory arrhythmia who underwent CSBRT with four-dimensional computed tomography (4DCT) and 4D cardiac CT (4DcCT). Proton plans were simulated using an IBA accelerator based on the 4D average CT. The prescription was 25 Gy in a single fraction, with all plans normalized to ensure that 95% of the target volume received the prescribed dose. 4D dose reconstruction was performed to generate 4D accumulated and dynamic doses. Furthermore, dose uncertainties due to the interplay effect of the substrate target and organs at risk (OARs) were assessed. The differences between internal organs at risk volume (IRV) and OARreal (manually contoured on average CT) were compared. In 4D dynamic dose, meeting prescription requirements entails V25 and D95 reaching 95% and 25 Gy, respectively.

Results

The 4D dynamic dose significantly differed from the 3D static dose. The mean V25 and D95 were 89.23% and 24.69 Gy, respectively, in 4DCT and 94.35% and 24.99 Gy, respectively, in 4DcCT. Eleven patients in 4DCT and six in 4DcCT failed to meet the prescription requirements. Critical organs showed varying dose increases. All metrics, except for Dmean and D50, significantly changed in 4DCT; in 4DcCT, only D50 remained unchanged with regards to the target dose uncertainties induced by the interplay effect. The interplay effect was only significant for the Dmax values of several OARs. Generally, respiratory motion caused a more pronounced interplay effect than cardiac pulsation. Neither IRV nor OARreal effectively evaluated the dose discrepancies of the OARs.

Conclusions

Complex cardiorespiratory motion can introduce dose uncertainties during IMPT. Motion-encompassing techniques may mitigate but cannot entirely compensate for the dose discrepancies. Individualized 4D dose assessments are recommended to verify the effectiveness and safety of CSBRT.

Motion-induced dose perturbations in photon radiotherapy and proton therapy measured by deformable liver-shaped 3D dosimeters in an anthropomorphic phantom

Background and purpose

The impact of intrafractional motion and deformations on clinical radiotherapy delivery has so far only been investigated by simulations as well as point and planar dose measurements. The aim of this study was to combine anthropomorphic 3D dosimetry with a deformable abdominal phantom to measure the influence of intra-fractional motion and gating in photon radiotherapy and evaluate the applicability in proton therapy.

Material and methods

An abdominal phantom was modified to hold a deformable anthropomorphic 3D dosimeter shaped as a human liver. A liver-specific photon radiotherapy and a proton pencil beam scanning therapy plan were delivered to the phantom without motion as well as with 12 mm sinusoidal motion while using either no respiratory gating or respiratory gating.

Results

Using the stationary irradiation as reference the local 3 %/2 mm 3D gamma index pass rate of the motion experiments in the planning target volume (PTV) was above 97 % (photon) and 78 % (proton) with gating whereas it was below 74 % (photon) and 45 % (proton) without gating.

Conclusions

For the first time a high-resolution deformable anthropomorphic 3D dosimeter embedded in a deformable abdominal phantom was applied for experimental validation of both photon and proton treatments of targets exhibiting respiratory motion. It was experimentally shown that gating improves dose coverage and the geometrical accuracy for both photon radiotherapy and proton therapy.

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.

A phantom to simulate organ motion and its effect on dose distribution in carbon ion therapy for pancreatic cancer

Objective. Carbon ion radiotherapy is a promising radiation technique for malignancies like pancreatic cancer. However, organs' motion imposes challenges for achieving homogeneous dose delivery. In this study, an anthropomorphicPancreasPhantom forIon-beamTherapy (PPIeT) was developed to simulate breathing and gastrointestinal motion during radiotherapy.Approach. The developed phantom contains a pancreas, two kidneys, a duodenum, a spine and a spinal cord. The shell of the organs was 3D printed and filled with agarose-based mixtures. Hounsfield Units (HU) of PPIeTs' organs were measured by CT. The pancreas motion amplitude in cranial-caudal (CC) direction was evaluated from patients' 4D CT data. Motions within the obtained range were simulated and analyzed in PPIeT using MRI. Additionally, GI motion was mimicked by changing the volume of the duodenum and quantified by MRI. A patient-like treatment plan was calculated for carbon ions, and the phantom was irradiated in a static and moving condition. Dose measurements in the organs were performed using an ionization chamber and dosimetric films.Main results. PPIeT presented tissue equivalent HU and reproducible breathing-induced CC displacements of the pancreas between (3.98 ± 0.36) mm and a maximum of (18.19 ± 0.44) mm. The observed maximum change in distance of (14.28 ± 0.12) mm between pancreas and duodenum was consistent with findings in patients. Carbon ion irradiation revealed homogenous coverage of the virtual tumor at the pancreas in static condition with a 1% deviation from the treatment plan. Instead, the dose delivery during motion with the maximum amplitude yielded an underdosage of 21% at the target and an increased uncertainty by two orders of magnitude.Significance. A dedicated phantom was designed and developed for breathing motion assessment of dose deposition during carbon ion radiotherapy. PPIeT is a unique tool for dose verification in the pancreas and its organs at risk during end-to-end tests.

Parameter based 4D dose calculations for proton therapy

Background and purpose Retrospective log file-based analysis provides the actual dose delivered based on the patient's breathing and the daily beam-delivery dynamics. To predict the motion sensitivity of the treatment plan on a patient-specific basis before treatment start a prospective tool is required. Such a parameter-based tool has been investigated with the aim to be used in clinical routine. Materials and Methods 4D dose calculations (4DDC) were performed for seven cancer patients with small breathing motion treated with scanned pulsed proton beams. Validation of the parameter-based 4DDC (p-4DDC) method was performed with an anthropomorphic phantom and patient data employing measurements and a log file-based 4DDC tool. The dose volume histogram parameters (Dx%) were investigated for the target and the organs at risk, compared to static and the file-based approach. Results The difference between the measured and the p-4DDC dose was within the deviation of the measurements. The maximum deviation was 0.4Gy. For the planning target volume D98% varied up to 15% compared to the static scenario, while the results from the log file and p-4DDC agreed within 2%. For the liver patients, D33%liver deviated up to 35% compared to static and 10% comparing the two 4DDC tools, while for the pancreas patients the D1%stomach varied up to 45% and 11%, respectively. Conclusion The results showed that p-4DDC could be used prospectively. The next step will be the clinical implementation of the p-4DDC tool, which can support a decision to either adapt the treatment plan or apply motion mitigation strategies.

Patient Breathing Motion and Delivery Specifics Influencing the Robustness of a Proton Pancreas Irradiation

Motion compensation strategies in particle therapy depend on the anatomy, motion amplitude and underlying beam delivery technology. This retrospective study on pancreas patients with small moving tumours analysed existing treatment concepts and serves as a basis for future treatment strategies for patients with larger motion amplitudes as well as the transition towards carbon ion treatments. The dose distributions of 17 hypofractionated proton treatment plans were analysed using 4D dose tracking (4DDT). The recalculation of clinical treatment plans employing robust optimisation for mitigating different organ fillings was performed on phased-based 4D computed tomography (4DCT) data considering the accelerator (pulsed scanned pencil beams delivered by a synchrotron) and the breathing-time structure. The analysis confirmed the robustness of the included treatment plans concerning the interplay of beam and organ motion. The median deterioration of D_50% (ΔD_50%) for the clinical target volume (CTV) and the planning target volume (PTV) was below 2%, while the only outlier was observed for ΔD_98% with -35.1%. The average gamma pass rate over all treatment plans (2%/ 2 mm) was 88.8% ± 8.3, while treatment plans for motion amplitudes larger than 1 mm performed worse. For organs at risk (OARs), the median ΔD_2% was below 3%, but for single patients, essential changes, e.g., up to 160% for the stomach were observed. The hypofractionated proton treatment for pancreas patients based on robust treatment plan optimisation and 2 to 4 horizontal and vertical beams showed to be robust against intra-fractional movements up to 3.7 mm. It could be demonstrated that the patient's orientation did not influence the motion sensitivity. The identified outliers showed the need for continuous 4DDT calculations in clinical practice to identify patient cases with more significant deviations.

Limitations of phase-sorting based pencil beam scanned 4D proton dose calculations under irregular motion

Objective.4D dose calculation (4DDC) for pencil beam scanned (PBS) proton therapy is typically based on phase-sorting of individual pencil beams onto phases of a single breathing cycle 4DCT. Understanding the dosimetric limitations and uncertainties of this approach is essential, especially for the realistic treatment scenario with irregular free breathing motion.Approach.For three liver and three lung cancer patient CTs, the deformable multi-cycle motion from 4DMRIs was used to generate six synthetic 4DCT(MRI)s, providing irregular motion (11/15 cycles for liver/lung; tumor amplitudes ∼4-18 mm). 4DDCs for two-field plans were performed, with the temporal resolution of the pencil beam delivery (4-200 ms) or with 8 phases per breathing cycle (500-1000 ms). For the phase-sorting approach, the tumor center motion was used to determine the phase assignment of each spot. The dose was calculated either using the full free breathing motion or individually repeating each single cycle. Additionally, the use of an irregular surrogate signal prior to 4DDC on a repeated cycle was simulated. The CTV volume with absolute dose differences >5% (V_dosediff>5%) and differences in CTVV_95%andD_5%-D_95%compared to the free breathing scenario were evaluated.Main results.Compared to 4DDC considering the full free breathing motion with finer spot-wise temporal resolution, 4DDC based on a repeated single 4DCT resulted inV_dosediff>5%of on average 34%, which resulted in an overestimation ofV_95%up to 24%. However, surrogate based phase-sorting prior to 4DDC on a single cycle 4DCT, reduced the averageV_dosediff>5%to 16% (overestimationV_95%up to 19%). The 4DDC results were greatly influenced by the choice of reference cycle (V_dosediff>5%up to 55%) and differences due to temporal resolution were much smaller (V_dosediff>5%up to 10%).Significance.It is important to properly consider motion irregularity in 4D dosimetric evaluations of PBS proton treatments, as 4DDC based on a single 4DCT can lead to an underestimation of motion effects.

Dosimetric robustness of lung tumor photon radiotherapy evaluated from multiple event CT imaging

PURPOSE

Intrafractional respiratory motion is a concern for lung tumor radiotherapy but full evaluation of its impact is hampered by the lack of images representing the true motion. This study presents a novel evaluation using free-breathing images acquired over realistic treatment times to study the dosimetric impact of respiratory motion in photon radiotherapy.

METHODS

Cine-CT images of 14 patients with lung cancer acquired during eight minutes of free-breathing at three occasions were used to simulate dose tracking of four different planning methods. These methods aimed to deliver 54 Gy in three fractions to D_50% of the target and were denoted as robust 4D (RB4), homogeneous fluence to the ITV (FLU) and an isodose prescription to the ITV with a high central dose (ISD), concurrently renormalized (IRN). Differences in dose coverage probability and homogeneity between the methods were quantified. Correlations between underdosage and attributes regarding the tumor and its motion were investigated.

RESULTS

Despite tumor motion amplitudes being larger than in the 4DCT all but FLU achieved the intended CTV D_50% for the cohort average. For all methods but IRN at least 93% of the patients would have received 95% of the intended dose. No differences in D_50% were found between RB4 and ISD nor IRN. However, RB4 led to better homogeneity.

CONCLUSIONS

Tumor motion in free-breathing not covered by the 4DCT had a small impact on dose. The RB4 is recommended for planning of free-breathing treatments. No factor was found that consistently correlated dose degradation with patient or motion attributes.

The Influence of Motion on the Delivery Accuracy When Comparing Actively Scanned Carbon Ions versus Protons at a Synchrotron-Based Radiotherapy Facility

Simple Summary

The interplay of breathing and beam motion reduces the efficacy of particle irradiation in moving tumours. The effect of motion on protons and carbon ion treatments was investigated dosimetrically and the results were benchmarked against each other by employing an anthropomorphic thorax phantom that was able to simulate tumour, rib, and lung motion. The critical question was whether target coverage and organ-at-risk sparing could be maintained when the application of simple motion mitigation was addressed. Special focus was put on unique synchrotron characteristics, such as pulsed beam delivery and beam intensity variations. It could be demonstrated that the effect of motion was greater for carbon ions than for protons. These findings demonstrated the need for applying motion mitigation techniques depending on the motion amplitude, particle type, and treatment prescription considering complex time correlations.

Abstract

Motion amplitudes, in need of mitigation for moving targets irradiated with pulsed carbon ions and protons, were identified to guide the decision on treatment and motion mitigation strategy. Measurements with PinPoint ionisation chambers positioned in an anthropomorphic breathing phantom were acquired to investigate different tumour motion scenarios, including rib and lung movements. The effect of beam delivery dynamics and spot characteristics was considered. The dose in the tumour centre was deteriorated up to 10% for carbon ions but only up to 5% for protons. Dose deviations in the penumbra increased by a factor of two when comparing carbon ions to protons, ranging from 2 to 30% for an increasing motion amplitude that was strongly dependent on the beam intensity. Layer rescanning was able to diminish the dose distortion caused by tumour motion, but an increase in spot size could reduce it even further to 5% within the target and 10% at the penumbra. An increased need for motion mitigation of carbon ions compared to protons was identified to assure target coverage and sparing of adjacent organs at risk in the penumbra region and outside the target. For the clinical implementation of moving target treatments at a synchrotron-based particle facility complex, time dependencies needed to be considered.