Feasibility of proton pencil beam scanning treatment of free-breathing lung cancer patients

Predicting the Effect of Proton Beam Therapy Technology on Pulmonary Toxicities for Patients With Locally Advanced Lung Cancer Enrolled in the Proton Collaborative Group Prospective Clinical Trial

PURPOSE

This study aimed to predict the probability of grade ≥2 pneumonitis or dyspnea within 12 months of receiving conventionally fractionated or mildly hypofractionated proton beam therapy for locally advanced lung cancer using machine learning.

METHODS AND MATERIALS

Demographic and treatment characteristics were analyzed for 965 consecutive patients treated for lung cancer with conventionally fractionated or mildly hypofractionated (2.2-3 Gy/fraction) proton beam therapy across 12 institutions. Three machine learning models (gradient boosting, additive tree, and logistic regression with lasso regularization) were implemented to predict Common Terminology Criteria for Adverse Events version 4 grade ≥2 pulmonary toxicities using double 10-fold cross-validation for parameter hyper-tuning without leak of information. Balanced accuracy and area under the curve were calculated, and 95% confidence intervals were obtained using bootstrap sampling.

RESULTS

The median age of the patients was 70 years (range, 20-97), and they had predominantly stage IIIA or IIIB disease. They received a median dose of 60 Gy in 2 Gy/fraction, and 46.4% received concurrent chemotherapy. In total, 250 (25.9%) had grade ≥2 pulmonary toxicity. The probability of pulmonary toxicity was 0.08 for patients treated with pencil beam scanning and 0.34 for those treated with other techniques (P = 8.97e-13). Use of abdominal compression and breath hold were highly significant predictors of less toxicity (P = 2.88e-08). Higher total radiation delivered dose (P = .0182) and higher average dose to the ipsilateral lung (P = .0035) increased the likelihood of pulmonary toxicities. The gradient boosting model performed the best of the models tested, and when demographic and dosimetric features were combined, the area under the curve and balanced accuracy were 0.75 ± 0.02 and 0.67 ± 0.02, respectively. After analyzing performance versus the number of data points used for training, we observed that accuracy was limited by the number of observations.

CONCLUSIONS

In the largest analysis of prospectively enrolled patients with lung cancer assessing pulmonary toxicities from proton therapy to date, advanced machine learning methods revealed that pencil beam scanning, abdominal compression, and lower normal lung doses can lead to significantly lower probability of developing grade ≥2 pneumonitis or dyspnea.

Experimental validation of a 4D dynamic dose calculation model for proton pencil beam scanning without spot time stamp considering free-breathing motion

PURPOSE

We developed a 4-dimensional dynamic dose (4DDD) calculation model for proton pencil beam scanning (PBS). This model incorporates the spill start time for all energies and uses the remaining irradiated spot time model instead of irradiated spot time logs. This study aimed to validate the calculation accuracy of a log file-based 4DDD model by comparing it with dose measurements performed under free-breathing conditions, thereby serving as an alternative approach to the conventional log file-based system.

METHODS

Three cubic verification plans were created using a heterogeneous block phantom; these plans were created using 10 phase 4D-CT datasets of the phantom. The CIRS dynamic platform was used to simulate motion with amplitudes of 2.5, 3.75, and 5.0 mm. These plans consisted of eight- and two-layered rescanning techniques. The lateral profiles were measured using a 2D ionization chamber array (2D-array) and EBT3 Gafchromic films at four starting phases, including three sinusoidal curves (periods of 3, 4, and 6 s) and a representative patient curve during actual treatment. 4DDDs were calculated using in-house scripting that assigned a time stamp to each spot and performed dose accumulation using deformable image registration. Furthermore, to evaluate the impact of parameter selection on our 4DDD model calculations, simulations were performed assuming a ±10% change in irradiation time stamp (0.8 ± 0.08 s) and spot scan speed. We evaluated the 2D gamma index and the absolute point doses between the calculated values and the measurements.

RESULTS

The 2D-array measurements revealed that the gamma scores for the static plans (no motion) and 4DDD plans exceeded 97.5% and 93.9% at 3%/3 mm, respectively. The average gamma score of the 4DDD plans was at least 96.1%. When using EBT3 films, the gamma scores of the 4DDD model exceeded 92.4% and 98.7% at 2%/2 mm and 3%/3 mm, respectively. Regarding the 4DDD point dose differences, more than 95% of the dose regions exhibited discrepancies within ±5.0% for 97.7% of the total points across all plans. The spot time assignment accuracy of our 4DDD model was acceptable even with ±10% sensitivity. However, the accuracy of the scan speed, when varied within ±10% sensitivity, was not acceptable (minimum gamma scores of 82.6% and maximum dose difference of 12.9%).

CONCLUSIONS

Our 4DDD calculations under free-breathing conditions using amplitudes of less than 5.0 mm were in good agreement with the measurements regardless of the starting phases, breathing curve patterns (between 3 and 6 s periods), and varying numbers of layered rescanning. The proposed system allows us to evaluate actual irradiated doses in various breathing periods, amplitudes, and starting phases, even on PBS machines without the ability to record spot logs.

Report of AAPM Task Group 290: Respiratory motion management for particle therapy

Dose uncertainty induced by respiratory motion remains a major concern for treating thoracic and abdominal lesions using particle beams. This Task Group report reviews the impact of tumor motion and dosimetric considerations in particle radiotherapy, current motion‐management techniques, and limitations for different particle‐beam delivery modes (i.e., passive scattering, uniform scanning, and pencil‐beam scanning). Furthermore, the report provides guidance and risk analysis for quality assurance of the motion‐management procedures to ensure consistency and accuracy, and discusses future development and emerging motion‐management strategies. This report supplements previously published AAPM report TG76, and considers aspects of motion management that are crucial to the accurate and safe delivery of particle‐beam therapy. To that end, this report produces general recommendations for commissioning and facility‐specific dosimetric characterization, motion assessment, treatment planning, active and passive motion‐management techniques, image guidance and related decision‐making, monitoring throughout therapy, and recommendations for vendors. Key among these recommendations are that: (1) facilities should perform thorough planning studies (using retrospective data) and develop standard operating procedures that address all aspects of therapy for any treatment site involving respiratory motion; (2) a risk‐based methodology should be adopted for quality management and ongoing process improvement.

How should we model and evaluate breathing interplay effects in IMPT?

Breathing interplay effects in Intensity Modulated Proton Therapy (IMPT) arise from the interaction between target motion and the scanning beam. Assessing the detrimental effect of interplay and the clinical robustness of several mitigation techniques requires statistical evaluation procedures that take into account the variability of breathing during dose delivery. In this study, we present such a statistical method to model intra-fraction respiratory motion based on breathing signals and assess clinical relevant aspects related to the practical evaluation of interplay in IMPT such as how to model irregular breathing, how small breathing changes affect the final dose distribution, and what is the statistical power (number of different scenarios) required for trustworthy quantification of interplay effects. First, two data-driven methodologies to generate artificial patient-specific breathing signals are compared: a simple sinusoidal model, and a precise probabilistic deep learning model generating very realistic samples of patient breathing. Second, we investigate the highly fluctuating relationship between interplay doses and breathing parameters, showing that small changes in breathing period result in large local variations in the dose. Our results indicate that using a limited number of samples to calculate interplay statistics introduces a bigger error than using simple sinusoidal models based on patient parameters or disregarding breathing hysteresis during the evaluation. We illustrate the power of the presented statistical method by analyzing interplay robustness of 4DCT and Internal Target Volume (ITV) treatment plans for a 8 lung cancer patients, showing that, unlike 4DCT plans, even 33 fraction ITV plans systematically fail to fulfill robustness requirements.

Validation of proton dose calculation on scatter corrected 4D cone beam computed tomography using a porcine lung phantom

Proton therapy treatment for lungs remains challenging as images enabling the detection of inter- and intra-fractional motion, which could be used for proton dose adaptation, are not readily available. 4D computed tomography (4DCT) provides high image quality but is rarely available in-room, while in-room 4D cone beam computed tomography (4DCBCT) suffers from image quality limitations stemming mostly from scatter detection. This study investigated the feasibility of using virtual 4D computed tomography (4DvCT) as a prior for a phase-per-phase scatter correction algorithm yielding a 4D scatter corrected cone beam computed tomography image (4DCBCT_cor), which can be used for proton dose calculation. 4DCT and 4DCBCT scans of a porcine lung phantom, which generated reproducible ventilation, were acquired with matching breathing patterns. Diffeomorphic Morphons, a deformable image registration algorithm, was used to register the mid-position 4DCT to the mid-position 4DCBCT and yield a 4DvCT. The 4DCBCT was reconstructed using motion-aware reconstruction based on spatial and temporal regularization (MA-ROOSTER). Successively for each phase, digitally reconstructed radiographs of the 4DvCT, simulated without scatter, were exploited to correct scatter in the corresponding CBCT projections. The 4DCBCT_corwas then reconstructed with MA-ROOSTER using the corrected CBCT projections and the same settings and deformation vector fields as those already used for reconstructing the 4DCBCT. The 4DCBCT_corand the 4DvCT were evaluated phase-by-phase, performing proton dose calculations and comparison to those of a ground truth 4DCT by means of dose-volume-histograms (DVH) and gamma pass-rates (PR). For accumulated doses, DVH parameters deviated by at most 1.7% in the 4DvCT and 2.0% in the 4DCBCT_corcase. The gamma PR for a (2%, 2 mm) criterion with 10% threshold were at least 93.2% (4DvCT) and 94.2% (4DCBCT_cor), respectively. The 4DCBCT_cortechnique enabled accurate proton dose calculation, which indicates the potential for applicability to clinical 4DCBCT scans.

Treatment planning and 4D robust evaluation strategy for proton therapy of lung tumors with large motion amplitude

Purpose

Intensity‐modulated proton therapy (IMPT) for lung tumors with a large tumor movement is challenging due to loss of robustness in the target coverage. Often an upper cut‐off at 5‐mm tumor movement is used for proton patient selection. In this study, we propose (1) a robust and easily implementable treatment planning strategy for lung tumors with a movement larger than 5 mm, and (2) a four‐dimensional computed tomography (4DCT) robust evaluation strategy for evaluating the dose distribution on the breathing phases.

Materials and methods

We created a treatment planning strategy based on the internal target volume (ITV) concept (aim 1). The ITV was created as a union of the clinical target volumes (CTVs) on the eight 4DCT phases. The ITV expanded by 2 mm was the target during robust optimization on the average CT (avgCT). The clinical plan acceptability was judged based on a robust evaluation, computing the voxel‐wise min and max (VWmin/max) doses over 28 error scenarios (range and setup errors) on the avgCT. The plans were created in RayStation (RaySearch Laboratories, Stockholm, Sweden) using a Monte Carlo dose engine, commissioned for our Mevion S250i Hyperscan system (Mevion Medical Systems, Littleton, MA, USA). We developed a new 4D robust evaluation approach (4DRobAvg; aim 2). The 28 scenario doses were computed on each individual 4DCT phase. For each scenario, the dose distributions on the individual phases were deformed to the reference phase and combined to a weighted sum, resulting in 28 weighted sum scenario dose distributions. From these 28 scenario doses, VWmin/max doses were computed. This new 4D robust evaluation was compared to two simpler 4D evaluation strategies: re‐computing the nominal plan on each individual 4DCT phase (4DNom) and computing the robust VWmin/max doses on each individual phase (4DRobInd). The treatment planning and dose evaluation strategies were evaluated for 16 lung cancer patients with tumor movement of 4–26 mm.

Results

The ratio of the ITV and CTV volumes increased linearly with the tumor amplitude, with an average ratio of 1.4. Despite large ITV volumes, a clinically acceptable plan fulfilling all target and organ at risk (OAR) constraints was feasible for all patients. The 4DNom and 4DRobInd evaluation strategies were found to under‐ or overestimate the dosimetric effect of the tumor movement, respectively. 4DRobInd showed target underdosage for five patients, not observed in the robust evaluation on the avgCT or in 4DRobAvg. The accuracy of dose deformation used in 4DRobAvg was quantified and found acceptable, with differences for the dose‐volume parameters below 1 Gy in most cases.

Conclusion

The proposed ITV‐based planning strategy on the avgCT was found to be a clinically feasible approach with adequate tumor coverage and no OAR overdosage even for large tumor movement. The new proposed 4D robust evaluation, 4DRobAvg, was shown to give an easily interpretable understanding of the effect of respiratory motion dose distribution, and to give an accurate estimate of the dose delivered in the different breathing phases.

Visually guided inspiration breath-hold facilitated with nasal high flow therapy in locally advanced lung cancer

BACKGROUND AND PURPOSE

Reducing breathing motion in radiotherapy (RT) is an attractive strategy to reduce margins and better spare normal tissues. The objective of this prospective study (NCT03729661) was to investigate the feasibility of irradiation of non-small cell lung cancer (NSCLC) with visually guided moderate deep inspiration breath-hold (IBH) using nasal high-flow therapy (NHFT).

MATERIAL AND METHODS

Locally advanced NSCLC patients undergoing photon RT were given NHFT with heated humidified air (flow: 40 L/min with 80% oxygen) through a nasal cannula. IBH was monitored by optical surface tracking (OST) with visual feedback. At a training session, patients had to hold their breath as long as possible, without and with NHFT. For the daily cone beam CT (CBCT) and RT treatment in IBH, patients were instructed to keep their BH as long as it felt comfortable. OST was used to analyze stability and reproducibility of the BH, and CBCT to analyze daily tumor position. Subjective tolerance was measured with a questionnaire at 3 time points.

RESULTS

Of 10 included patients, 9 were treated with RT. Seven (78%) completed the treatment with NHFT as planned. At the training session, the mean BH length without NHFT was 39 s (range 15-86 s), and with NHFT 78 s (range 29-223 s) (p = .005). NHFT prolonged the BH duration by a mean factor of 2.1 (range 1.1-3.9s). The mean overall stability and reproducibility were within 1 mm. Subjective tolerance was very good with the majority of patients having no or minor discomfort caused by the devices. The mean inter-fraction tumor position variability was 1.8 mm (-1.1-8.1 mm;SD 2.4 mm).

CONCLUSION

NHFT for RT treatment of NSCLC in BH is feasible, well tolerated and significantly increases the breath-hold duration. Visually guided BH with OST is stable and reproducible. We therefore consider this an attractive patient-friendly approach to treat lung cancer patients with RT in BH.

Pro-con of proton: Dosimetric advantages of intensity-modulation over passive scatter for thoracic malignancies

•

Intensity Modulated Proton Therapy (IMPT) results in significant reduction of dose to organ at risk.

•

Improving plan robustness mitigates interplay effects.

•

Blanket use of small spots on a group of patients may severely worsen interplay in selected patients.

•

Hypofractionated regimes have fewer interplay effects in both fractional and overall simulations.

•

Randomised control trials are required before any clinical benefit of IMPT can be confirmed.

The use of passively scattered proton therapy (PSPT) or intensity modulated proton therapy (IMPT) opens the potential for dose escalation or critical structure sparing in thoracic malignancies. While the latter offers greater dose conformality, dose distributions are subjected to greater uncertainties, especially due to interplay effects. Exploration in this area is warranted to determine if there is any dosimetric advantages in using IMPT for thoracic malignancies. This review aims to both compare organs-at-risk sparing and plan robustness between PSPT and IMPT and examine the mitigation strategies for the reduction of interplay effects currently available. Early evidence suggests that IMPT is dosimetrically superior to PSPT in thoracic malignancies. Randomised control trials are required before any clinical benefit of IMPT can be confirmed.

Investigation of inter-fraction target motion variations in the context of pencil beam scanned proton therapy in non-small cell lung cancer patients

Purpose

For locally advanced‐stage non‐small cell lung cancer (NSCLC), inter‐fraction target motion variations during the whole time span of a fractionated treatment course are assessed in a large and representative patient cohort. The primary objective is to develop a suitable motion monitoring strategy for pencil beam scanning proton therapy (PBS‐PT) treatments of NSCLC patients during free breathing.

Methods

Weekly 4D computed tomography (4DCT; 41 patients) and daily 4D cone beam computed tomography (4DCBCT; 10 of 41 patients) scans were analyzed for a fully fractionated treatment course. Gross tumor volumes (GTVs) were contoured and the 3D displacement vectors of the centroid positions were compared for all scans. Furthermore, motion amplitude variations in different lung segments were statistically analyzed. The dosimetric impact of target motion variations and target motion assessment was investigated in exemplary patient cases.

Results

The median observed centroid motion was 3.4 mm (range: 0.2–12.4 mm) with an average variation of 2.2 mm (range: 0.1–8.8 mm). Ten of 32 patients (31.3%) with an initial motion <5 mm increased beyond a 5‐mm motion amplitude during the treatment course. Motion observed in the 4DCBCT scans deviated on average 1.5 mm (range: 0.0–6.0 mm) from the motion observed in the 4DCTs. Larger motion variations for one example patient compromised treatment plan robustness while no dosimetric influence was seen due to motion assessment biases in another example case.

Conclusions

Target motion variations were investigated during the course of radiotherapy for NSCLC patients. Patients with initial GTV motion amplitudes of < 2 mm can be assumed to be stable in motion during the treatment course. For treatments of NSCLC patients who exhibit motion amplitudes of > 2 mm, 4DCBCT should be considered for motion monitoring due to substantial motion variations observed.

Photons or protons for reirradiation in (non-)small cell lung cancer: Results of the multicentric ROCOCO in silico study

OBJECTIVE

Locally recurrent disease is of increasing concern in (non-)small cell lung cancer [(N)SCLC] patients. Local reirradiation with photons or particles may be of benefit to these patients. In this multicentre in silico trial performed within the Radiation Oncology Collaborative Comparison (ROCOCO) consortium, the doses to the target volumes and organs at risk (OARs) were compared when using several photon and proton techniques in patients with recurrent localised lung cancer scheduled to undergo reirradiation.

METHODS

24 consecutive patients with a second primary (N)SCLC or recurrent disease after curative-intent, standard fractionated radio(chemo)therapy were included in this study. The target volumes and OARs were centrally contoured and distributed to the participating ROCOCO sites. Remaining doses to the OARs were calculated on an individual patient's basis. Treatment planning was performed by the participating site using the clinical treatment planning system and associated beam characteristics.

RESULTS

Treatment plans for all modalities (five photon and two proton plans per patient) were available for 22 patients (N = 154 plans). 3D-conformal photon therapy and double-scattered proton therapy delivered significantly lower doses to the target volumes. The highly conformal techniques, i.e., intensity modulated radiation therapy (IMRT), volumetric modulated arc therapy (VMAT), CyberKnife, TomoTherapy and intensity-modulated proton therapy (IMPT), reached the highest doses in the target volumes. Of these, IMPT was able to statistically significantly decrease the radiation doses to the OARs.

CONCLUSION

Highly conformal photon and proton beam techniques enable high-dose reirradiation of the target volume. They, however, significantly differ in the dose deposited in the OARs. The therapeutic options, i.e., reirradiation or systemic therapy, need to be carefully weighed and discussed with the patients.

ADVANCES IN KNOWLEDGE

Highly conformal photon and proton beam techniques enable high-dose reirradiation of the target volume. In light of the abilities of the various highly conformal techniques to spare specific OARs, the therapeutic options need to be carefully weighed and patients included in the decision-making process.

An analytical model for the upper bound estimation of respiratory motion-induced dose uncertainty in spot-scanning proton beam therapy

PURPOSE

We developed an analytical model of a spot-scanning beam delivery system to estimate the upper bound of respiratory motion-induced dose uncertainty for a given treatment plan.

METHODS

The effective delivery time for each spot position in the treatment plan was calculated on the basis of the parameters of the delivery system. The upper bound of the dose uncertainty was then calculated as a function of the effective delivery time. Two-dimensional (2D) measurements with a detector array on a one-dimensional moving platform were obtained to validate the model.

RESULTS

We performed 351 two-dimensional measurements on a moving platform for different delivery sequences of a single-layer uniform pattern and patient treatment field. The measured dose uncertainty was a strong function of the effective delivery time: The shortest effective delivery time resulted in a maximum absolute dose error of >90%, while the longest ones resulted in a maximum absolute dose error of 4.9% for a single layer and 9.7% for a patient field with heterogeneity. The relationship of the effective delivery time and the measured dose uncertainty followed the analytical formula.

CONCLUSIONS

With our analytical model, the upper bound of the dose uncertainty due to motion can be estimated in spot-scanning proton therapy without four-dimensional dynamic dose calculation.

Energy layer optimization strategies for intensity-modulated proton therapy of lung cancer patients

PURPOSE

When treating lung cancer patients with intensity-modulated proton therapy (IMPT), target coverage can only be guaranteed when utilizing motion mitigation. The three motion mitigation techniques, gating, breath-hold, and dose repainting, all benefit from a more rapid application of the treatment plan. A lower limit for the ungated treatment time is defined by the number of energy layers in the IMPT plan. By limiting this number during treatment planning, IMPT could become more viable for lung cancer patients. We investigate to what extend the number of layers can be reduced in single-field optimization (SFO) and multifield optimization (MFO) plans and which implications it has on the plan quality and robustness.

METHODS

We have implemented three distinct layer-reducing strategies in the treatment planning system Hyperion; constant energy steps, exponential energy steps, and an adaptive strategy, where the spot weights are exposed to a group sparsity penalty in combination with layer exclusion during optimization. Four levels of increasing layer removal are planned for each strategy. SFO and MFO plans with three treatment fields are created for eleven locally advanced NSCLC patients on the midventilation 4DCT phase to simulate a breath-hold. A minimum dose to the target is ensured for each degree of layer reduction, reflecting the plan quality in the homogeneity index (HI). Plan quality was also assessed by a robustness evaluation, where the patient setup was shifted 2 mm or 4 mm in six directions.

RESULTS

The three strategies result in very similar target coverages and robustness levels as a function of removed layers. The HI increases unacceptably for all the SFO plans after 50% layer removal as compared to the reference plan, while all the MFO plans are clinically acceptable with up to a highest removed percentage of 75%. The robustness level is constant as a function of removed layers. The SFO plans are significantly more robust than the MFO plans with all P-values below 0.001 (Wilcoxon signed-rank). The overall mean D98% CTV dose difference is at 2-mm setup error amplitude: 0.7 Gy (SFO) and 1.9 Gy (MFO), and at 4 mm: 3.2 Gy (SFO) and 5.4 Gy (MFO), respectively.

CONCLUSIONS

The number of layers in MFO plans can be reduced substantially more than in SFO plans without compromising plan quality. Furthermore, as the robustness is independent of the number of layers, it follows that if the level of robustness is acceptable or enforced via robust optimization, MFO plans could be candidates for treatment time reductions via energy layer reductions.