Target volume dose considerations in proton beam treatment planning for lung tumors

Comprehensive beam delivery latency evaluation for gated proton therapy system using customized multi-channel signal acquisition platform

PURPOSE

Beam delivery latency in respiratory-gated particle therapy systems is a crucial issue to dose delivery accuracy. The aim of this study is to develop a multi-channel signal acquisition platform for investigating gating latencies occurring within RPM respiratory gating system (Varian, USA) and ProBeam proton treatment system (Varian, USA) individually.

METHODS

The multi-channel signal acquisition platform consisted of several electronic components, including a string position sensor for target motion detection, a photodiode for proton beam sensing, an interfacing board for accessing the trigger signal between the respiratory gating system and the proton treatment system, a signal acquisition device for sampling and synchronizing signals from the aforementioned components, and a laptop for controlling the signal acquisition device and data storage. RPM system latencies were determined by comparing the expected gating phases extracted from the motion signal with the trigger signal's state turning points. ProBeam system latencies were assessed by comparing the state turning points of the trigger signal with the beam signal. The total beam delivery latencies were calculated as the sum of delays in the respiratory gating system and the cyclotron proton treatment system. During latency measurements, simulated sinusoidal motion were applied at different amplitudes and periods for complete beam delivery latency evaluation under different breathing patterns. Each breathing pattern was repeated 30 times for statistical analysis.

RESULTS

The measured gating ON/OFF latencies in the RPM system were found to be 104.20 ± 13.64 ms and 113.60 ± 14.98 ms, respectively. The measured gating ON/OFF delays in the ProBeam system were 108.29 ± 0.85 ms and 1.20 ± 0.04 ms, respectively. The total beam ON/OFF latencies were determined to be 212.50 ± 13.64 ms and 114.80 ± 14.98 ms.

CONCLUSION

With the developed multi-channel signal acquisition platform, it was able to investigate the gating lags happened in both the respiratory gating system and the proton treatment system. The resolution of the platform is enough to distinguish the delays at the millisecond time level. Both the respiratory gating system and the proton treatment system made contributions to gating latency. Both systems contributed nearly equally to the total beam ON latency, with approximately 100 ms. In contrast, the respiratory gating system was the dominant contributor to the total beam OFF latency.

Inter- and intrafractional 4D dose accumulation for evaluating ΔNTCP robustness in lung cancer

BACKGROUND AND PURPOSE

Model-based selection of proton therapy patients relies on a predefined reduction in normal tissue complication probability (NTCP) with respect to photon therapy. The decision is necessarily made based on the treatment plan, but NTCP can be affected when the delivered treatment deviates from the plan due to delivery inaccuracies. Especially for proton therapy of lung cancer, this can be important because of tissue density changes and, with pencil beam scanning, the interplay effect between the proton beam and breathing motion.

MATERIALS AND METHODS

In this work, we verified whether the expected benefit of proton therapy is retained despite delivery inaccuracies by reconstructing the delivered treatment using log-file based dose reconstruction and inter- and intrafractional accumulation. Additionally, the importance of two uncertain parameters for treatment reconstruction, namely deformable image registration (DIR) algorithm and α/β ratio, was assessed.

RESULTS

The expected benefit or proton therapy was confirmed in 97% of all studied cases, despite regular differences up to 2 percent point (p.p.) NTCP between the delivered and planned treatments. The choice of DIR algorithm affected NTCP up to 1.6 p.p., an order of magnitude higher than the effect of α/β ratio.

CONCLUSION

For the patient population and treatment technique employed, the predicted clinical benefit for patients selected for proton therapy was confirmed for 97.0% percent of all cases, although the NTCP based proton selection was subject to 2 p.p. variations due to delivery inaccuracies.

A proof-of-concept study of an in-situ partial-ring time-of-flight PET scanner for proton beam verification

Development of a PET system capable of in-situ imaging requires a design that can accommodate the proton treatment beam nozzle. Among the several PET instrumentation approaches developed thus far, the dual-panel PET scanner is often used as it is simpler to develop and integrate within the proton therapy gantry. Partial-angle coverage of these systems can however lead to limited-angle artefacts in the reconstructed PET image. We have previously demonstrated via simulations that time-of-flight (TOF) reconstruction reduces the artifacts accompanying limited-angle data, and permits proton range measurement with 1-2 mm accuracy and precision. In this work we show measured results from a small proof-of-concept dual-panel PET system that uses TOF information to reconstruct PET data acquired after proton irradiation. The PET scanner comprises of two detector modules, each comprised of an array of 4×4×30 mm³ lanthanum bromide scintillator. Measurements are performed with an oxygen-rich gel-water, an adipose tissue equivalent material, and in vitro tissue phantoms. For each phantom measurement, 2 Gy dose was deposited using 54 - 100 MeV proton beams. For each phantom, a Monte Carlo simulation generating the expected distribution of PET isotope from the corresponding proton irradiation was also performed. Proton range was calculated by drawing multiple depth-profiles over a central region encompassing the proton dose deposition. For each profile, proton range was calculated using two techniques (a) 50% pick-off from the distal edge of the profile, and (b) comparing the measured and Monte Carlo profile to minimize the absolute sum of differences over the entire profile. A 10 min PET acquisition acquired with minimal delay post proton-irradiation is compared with a 10 min PET scan acquired after a 20 min delay. Measurements show that PET acquisition with minimal delay is necessary to collect ¹⁵O signal, and maximize ¹¹C signal collection with a short PET acquisition. In comparison with the 50% pick-off technique, the shift technique is more robust and offers better precision in measuring the proton range for the different phantoms. Range measurements from PET images acquired with minimal delay, and the shift technique demonstrate the ability to achieve <1.5 mm accuracy and precision in estimating proton range.

Diaphragm-Based Position Verification to Improve Daily Target Dose Coverage in Proton and Photon Radiation Therapy Treatment of Distal Esophageal Cancer

PURPOSE

In modern conformal radiation therapy of distal esophageal cancer, target coverage can be affected by variations in the diaphragm position. We investigated if daily position verification (PV) extended by a diaphragm position correction would optimize target dose coverage for esophageal cancer treatment.

METHODS AND MATERIALS

For 15 esophageal cancer patients, intensity modulated proton therapy (IMPT) and volumetric modulated arc therapy (VMAT) plans were computed. Displacements of the target volume were correlated with diaphragm displacements using repeated 4-dimensional computed tomography images to determine the correction needed to account for diaphragm variations. Afterwards, target coverage was evaluated for 3 PV approaches based on: (1) bony anatomy (PV_B), (2) bony anatomy corrected for the diaphragm position (PV_BD) and (3) target volume (PV_T).

RESULTS

The cranial-caudal mean target displacement was congruent with almost half of the diaphragm displacement (y = 0.459x), which was used for the diaphragm correction in PV_BD. Target dose coverage using PV_B was adequate for most patients with diaphragm displacements up till 10 mm (≥94% of the dose in 98% of the volume [D_98%]). For larger displacements, the target coverage was better maintained by PV_T and PV_BD. Overall, PV_BD accounted best for target displacements, especially in combination with tissue density variations (D_98%: IMPT 94% ± 5%, VMAT 96% ± 5%). Diaphragm displacements of more than 10 mm were observed in 22% of the cases.

CONCLUSIONS

PV_B was sufficient to achieve adequate target dose coverage in case of small deviations in diaphragm position. However, large deviations of the diaphragm were best mitigated by PV_BD. To detect the cases where target dose coverage could be compromised due to diaphragm position variations, we recommend monitoring of the diaphragm position before treatment through online imaging.

Retrospective Planning Study of Patients with Superior Sulcus Tumours Comparing Pencil Beam Scanning Protons to Volumetric-Modulated Arc Therapy

AIMS

Twenty per cent of patients with non-small cell lung cancer present with stage III locally advanced disease. Precision radiotherapy with pencil beam scanning (PBS) protons may improve outcomes. However, stage III is a heterogeneous group and accounting for complex tumour motion is challenging. As yet, it remains unclear as to whom will benefit. In our retrospective planning study, we explored if patients with superior sulcus tumours (SSTs) are a select cohort who might benefit from this treatment.

MATERIALS AND METHODS

Patients with SSTs treated with radical radiotherapy using four-dimensional planning computed tomography between 2010 and 2015 were identified. Tumour motion was assessed and excluded if greater than 5 mm. Photon volumetric-modulated arc therapy (VMAT) and PBS proton single-field optimisation plans, with and without inhomogeneity corrections, were generated retrospectively. Robustness analysis was assessed for VMAT and PBS plans involving: (i) 5 mm geometric uncertainty, with an additional 3.5% range uncertainty for proton plans; (ii) verification plans at maximal inhalation and exhalation. Comparative dosimetric and robustness analyses were carried out.

RESULTS

Ten patients were suitable. The mean clinical target volume D95 was 98.1% ± 0.4 (97.5-98.8) and 98.4% ± 0.2 (98.1-98.9) for PBS and VMAT plans, respectively. All normal tissue tolerances were achieved. The same four PBS and VMAT plans failed robustness assessment. Inhomogeneity corrections minimally impacted proton plan robustness and made it worse in one case. The most important factor affecting target coverage and robustness was the clinical target volume entering the spinal canal. Proton plans significantly reduced the mean lung dose (by 21.9%), lung V5, V10, V20 (by 47.9%, 36.4%, 12.1%, respectively), mean heart dose (by 21.4%) and thoracic vertebra dose (by 29.2%) (P < 0.05).

CONCLUSIONS

In this planning study, robust PBS plans were achievable in carefully selected patients. Considerable dose reductions to the lung, heart and thoracic vertebra were possible without compromising target coverage. Sparing these lymphopenia-related organs may be particularly important in this era of immunotherapy.

Dosimetric impact of geometric distortions in an MRI-only proton therapy workflow for lung, liver and pancreas

In a radiation therapy workflow based on Magnetic Resonance Imaging (MRI), dosimetric errors may arise due to geometric distortions introduced by MRI. The aim of this study was to quantify the dosimetric effect of system-dependent geometric distortions in an MRI-only workflow for proton therapy applied at extra-cranial sites. An approach was developed, in which computed tomography (CT) images were distorted using an MRI displacement map, which represented the MR distortions in a spoiled gradient-echo sequence due to gradient nonlinearities and static magnetic field inhomogeneities. A retrospective study was conducted on 4DCT/MRI digital phantoms and 18 4DCT clinical datasets of the thoraco-abdominal site. The treatment plans were designed and separately optimized for each beam in a beam specific Planning Target Volume on the distorted CT, and the final dose distribution was obtained as the average. The dose was then recalculated in undistorted CT using the same beam geometry and beam weights. The analysis was performed in terms of Dose Volume Histogram (DVH) parameters. No clinically relevant dosimetric impact was observed on organs at risk, whereas in the target structure, geometric distortions caused statistically significant variations in the planned dose DVH parameters and dose homogeneity index (DHI). The dosimetric variations in the target structure were smaller in abdominal cases (ΔD_2%, ΔD_98%, and ΔD_mean all below 0.1% and ΔDHI below 0.003) compared to the lung cases. Indeed, lung patients with tumors isolated inside lung parenchyma exhibited higher dosimetric variations (ΔD_2%≥0.3%, ΔD_98%≥15.9%, ΔD_mean≥3.3% and ΔDHI≥0.102) than lung patients with tumor close to soft tissue (ΔD_2%≤0.4%, ΔD_98%≤5.6%, ΔD_mean≤0.9% and ΔDHI≤0.027) potentially due to higher density variations along the beam path. Results suggest the potential applicability of MRI-only proton therapy, provided that specific analysis is applied for isolated lung tumors.

Organ motion quantification and margins evaluation in carbon ion therapy of abdominal lesions

PURPOSE

In image-guided particle radiotherapy of abdominal lesions, respiratory motion hinders treatment accuracy. In this study, 2D cineMRI data were used to quantify the tumor (GTV) motion and to evaluate the clinical approach based on deriving an internal target volume (ITV) from a planning 4DCT for gating treatments.

METHODS

Seven patients with abdominal lesions were treated with carbon-ion therapy at the National Centre of Oncological Hadron-therapy (Italy). The MR scan was performed on the same day of the 4DCT acquisition. For four patients, an additional MR was acquired approximately after 1 week. The cineMRI combined with deformable image registration algorithm was used to quantify tumor motion. Afterwards, two ITVs were defined considering (1) all phases (ITVFB) and (2) only phases within the gating window (ITVG), and then compared with the clinical (4DCT-derived) ITVs (ITVCG and ITVCFB).

RESULTS

Tumor residual motion estimated by cineMRI data in the two MRI sessions resulted not significantly different from 4DCT, although cineMRI accounted for cycle-to-cycle variations. The ITV normalized for the GTV median values were higher for ITVFB with respect to ITVG, ITVCFB and ITVCG. The Hausdorff distances with respect to the GTV were up to 10.55 mm, 3.13 mm, 5.56 mm and 2.51 mm, for ITVFB, ITVG, ITVCFB and ITVCG, respectively. According to both metrics, ITVCG and ITVG were not found significantly different.

CONCLUSIONS

CineMRI acquisitions allowed to quantify organ motion without delivering additional dose to the patient and to verify treatment margins in gated carbon-ion therapy of abdominal lesions.

Measurement-based study on characterizing symmetric and asymmetric respiratory motion interplay effect on target dose distribution in the proton pencil beam scanning

Pencil beam scanning proton therapy makes possible intensity modulation, resulting in improved target dose conformity and organ‐at‐risk (OAR) dose sparing. This benefit, however, results in increased sensitivity to certain clinical and beam delivery parameters, such as respiratory motion. These effects can cause plan degeneration, which could lead to decreased tumor dose or increased OAR dose. This study evaluated the measurements of proton pencil beam scanning delivery made with a 2D ion chamber array in solid water on a 1D motion platform, where respiratory motion was simulated using sine and cosine⁴ waves representing sinusoidal symmetric and realistic asymmetric breathing motions, respectively. Motion amplitudes were 0.5 cm and 1 cm corresponding to 1 cm and 2 cm of maximum respiratory excursions, respectively, with 5 sec fixed breathing cycle. The treatment plans were created to mimic spherical targets of 3 cm or 10 cm diameter located at 5 cm or 1 cm depth in solid water phantom. A reference RBE dose of 200 cGy per fraction was delivered in 1, 5, 10, and 15 fractions for each dataset. We evaluated dose conformity and uniformity at the center plane of targets by using the Conformation Number and the Homogeneity Index, respectively. Results indicated that dose conformity as well as homogeneity was more affected by motion for smaller targets. Dose conformity was better achieved for symmetric breathing patterns than asymmetric breathing patterns regardless of the number of fractions. The presence of a range shifter with shallow targets reduced the motion effect by improving dose homogeneity. While motion effects are known to be averaged out over the course of multifractional treatments, this might not be true for proton pencil beam scanning under asymmetrical breathing pattern.

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

BACKGROUND AND PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSIONS

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

Validation of a model for physical dose variations in irregularly moving targets treated with carbon ion beams

PURPOSE

In particle therapy, conventional treatment planning systems rely on an imaging representation of the irradiated region to compute the dose. For irregular breathing, when an imaging dataset describing the actual motion is not available, a different approach for dose estimation is needed. To this aim, we validate a method for the estimation of physical dose variations in gated carbon ion treatments, providing also a demonstration of the feasibility of physical dose metrics to assess the method performance. Finally, we describe a sample use case, in which this method is used to assess plan robustness with respect to undetected irregular tumor motion.

METHODS

The method entails the definition of a patient- and beam-specific water equivalent depth (WED) space, the simulation of motion as a translation equal to tumor displacement, and the reconstruction of the altered dose. We validated the approach using four-dimensional computed tomographies (4DCTs) and clinical plans in 12 patients, treated with respiratory gated carbon ion beams at the National Centre for Oncological Hadrontherapy (Pavia, Italy). Using the end-exhale CT and dose distribution as a reference, the physical dose delivered at the end-inhale tumor position was estimated and compared to the ground-truth dose recalculation on the end-inhale CT. Biologically effective and physical dose variations between the plan and the recalculation were compared as well. As a use case, we evaluated dose changes caused by simulated irregular tumor motion, that is, linear and nonlinear baseline shifts and/or amplitude variations with hysteresis.

RESULTS

The ratio between biologically effective and physical equivalent uniform dose (EUD) variations due to end-exhale to end-inhale motion was less than one for 96% of investigated structures. In the validation study, we found a median error corresponding to a 14% EUD overestimation for the tumor and 4% EUD underestimation for a subgroup of organs at risk, together with a high EUD variation due to motion [median 352% EUD variation between end-exhale and end-inhale doses in the planning tumor volume (PTV)]. Considering relevant dose-volume histogram (DVH) metrics, the median difference between estimated and ground truth doses was ≤ 4%. Gamma analysis between estimated and recalculated dose distributions resulted in a pass rate > 80% for 83% of the target volumes. For the two patients selected for the sample use case, a patient-specific assessment of the method performance was performed on the 4DCT and it was possible to relate EUD variations of both tumor and organs at risk to the simulated target motion.

CONCLUSIONS

The physical dose distribution was found to be more sensitive to motion with respect to the biologically effective one, suggesting the suitability of the physical dose metrics for the WED-space method validation. We showed that the method can compensate for intra-fractional tumor motion with proper accuracy in the selected patient group, although its use is recommended when limited deformations are expected. In conclusion, the WED-space method can provide simulations of dose alteration due to irregular breathing when imaging data are lacking, and, once integrated with relative biological effectiveness (RBE) modeling, it would be useful in evaluating the robustness of carbon ion treatment plans.

Potential for Improvements in Robustness and Optimality of Intensity-Modulated Proton Therapy for Lung Cancer with 4-Dimensional Robust Optimization

BACKGROUND

Major challenges in the application of intensity-modulated proton therapy (IMPT) for lung cancer patients include the uncertainties associated with breathing motion, its mitigation and its consideration in IMPT optimization. The primary objective of this research was to evaluate the potential of four-dimensional robust optimization (4DRO) methodology to make IMPT dose distributions resilient to respiratory motion as well as to setup and range uncertainties; Methods: The effect of respiratory motion, characterized by different phases of 4D computed tomography (4DCT), was incorporated into an in-house 4DRO system. Dose distributions from multiple setup and range uncertainty scenarios were calculated for each of the ten phases of CT datasets. The 4DRO algorithm optimizes dose distributions to achieve target dose coverage and normal tissue sparing for multiple setup and range uncertainty scenarios as well as for all ten respiratory phases simultaneously. IMPT dose distributions of ten lung cancer patients with different tumor sizes and motion magnitudes were optimized to illustrate our approach and its potential; Results: Compared with treatment plans generated using the conventional planning target volume (PTV)-based optimization and 3D robust optimization (3DRO), plans generated by 4DRO were found to have superior clinical target volume coverage and dose robustness in the face of setup and range uncertainties as well as for respiratory motion. In most of the cases we studied, 4DRO also resulted in more homogeneous target dose distributions. Interestingly, such improvements were found even for cases in which moving diaphragms intruded into the proton beam paths; Conclusion: The incorporation of respiratory motion, along with setup and range uncertainties, into robust optimization, has the potential to improve the resilience of target and normal tissue dose distributions in IMPT plans in the face of the uncertainties considered. Moreover, it improves the optimality of plans compared to PTV-based optimization as well as 3DRO.

A quality assurance for respiratory gated proton irradiation with range modulation wheel

The purpose of this study was to provide periodic quality assurance (QA) methods for respiratory-gated proton beam with a range modulation wheel (RMW) and to clarify the characteristics and long-term stability of the respiratory-gated proton beam. A two-dimensional detector array and a solid water phantom were used to measure absolute dose, spread-out Bragg peak (SOBP) width and proton range for monthly QA. SOBP width and proton range were measured using an oblique incidence beam to the lateral side of a solid water phantom and compared between with and without a gating proton beam. To measure the delay time of beam-on/off for annual QA, we collected the beam-on/off signals and the dose monitor-detected pulse. We analyzed the results of monthly QA over a 15-month period and investigated the delay time by machine signal analysis. The dose deviations at proximal, SOBP center and distal points were -0.083 ± 0.25%, 0.026 ± 0.20%, and -0.083 ± 0.35%, respectively. The maximum dose deviation between with and without respiratory gating was -0.95% at the distal point and other deviations were within ±0.5%. Proximal and SOBP center doses showed the same trend over a 15-month period. Delay times of beam-on/off for 200 MeV/SOBP 16 cm were 140.5 ± 0.8 ms and 22.3 ± 13.0 ms, respectively. Delay times for 160 MeV/SOBP 10 cm were 167.5 ± 15.1 ms and 19.1 ± 9.8 ms. Our beam delivery system with the RMW showed sufficient stability for respiratory-gated proton therapy and the system did not show dependency on the energy and the respiratory wave form. The delay times of beam-on/off were within expectations. The proposed QA methods will be useful for managing the quality of respiratory-gated proton beams and other beam delivery systems.

ITV-Based Robust Optimization for VMAT Planning of Stereotactic Body Radiation Therapy of Lung Cancer

PURPOSE

Using planning target volume (PTV) to account for setup uncertainties in stereotactic body radiation therapy (SBRT) of lung cancer has been questioned because a significant portion of the PTV contains low-density lung tissue. The purpose of this study is to (1) investigate the feasibility of using robust optimization to account for setup uncertainties in volumetric modulated arc therapy plan for lung SBRT and (2) evaluate the potential normal tissue-sparing benefit of a robust optimized plan compared with a conventional PTV-based optimized plan.

METHODS AND MATERIALS

The study was conducted with both phantom and patient cases. For each patient or phantom, 2 SBRT lung volumetric modulated arc therapy plans were generated, including an optimized plan based on the PTV (PTV-based plan) with a 5-mm internal target volume (ITV)-to-PTV margin and a second plan based on robust optimization of ITV (ITV-based plan) with ±5-mm setup uncertainties. The target coverage was evaluated on ITV D99 in 15 scenarios that simulated a 5-mm setup error. Dose-volume information on normal lung tissue, intermediate-to-high dose spillage, and integral dose was evaluated.

RESULTS

Compared with PTV-based plans, ITV-based robust optimized plans resulted in lower normal lung tissue dose, lower intermediate-to-high dose spillage to the body, and lower integral dose, while preserving the dose coverage under setup error scenarios for both phantom and patient cases.

CONCLUSIONS

Using ITV-based robust optimization, we have shown that accounting for setup uncertainty in SBRT planning is feasible. Further clinical studies are warranted to confirm the clinical effectiveness of this novel approach.

Towards a Clinical Decision Support System for External Beam Radiation Oncology Prostate Cancer Patients: Proton vs. Photon Radiotherapy? A Radiobiological Study of Robustness and Stability

We present a methodology which can be utilized to select proton or photon radiotherapy in prostate cancer patients. Four state-of-the-art competing treatment modalities were compared (by way of an in silico trial) for a cohort of 25 prostate cancer patients, with and without correction strategies for prostate displacements. Metrics measured from clinical image guidance systems were used. Three correction strategies were investigated; no-correction, extended-no-action-limit, and online-correction. Clinical efficacy was estimated via radiobiological models incorporating robustness (how probable a given treatment plan was delivered) and stability (the consistency between the probable best and worst delivered treatments at the 95% confidence limit). The results obtained at the cohort level enabled the determination of a threshold for likely clinical benefit at the individual level. Depending on the imaging system and correction strategy; 24%, 32% and 44% of patients were identified as suitable candidates for proton therapy. For the constraints of this study: Intensity-modulated proton therapy with online-correction was on average the most effective modality. Irrespective of the imaging system, each treatment modality is similar in terms of robustness, with and without the correction strategies. Conversely, there is substantial variation in stability between the treatment modalities, which is greatly reduced by correction strategies. This study provides a ‘proof-of-concept’ methodology to enable the prospective identification of individual patients that will most likely (above a certain threshold) benefit from proton therapy.

Density overwrites of internal tumor volumes in intensity modulated proton therapy plans for mobile lung tumors

The purpose of this study was to investigate internal tumor volume density overwrite strategies to minimize intensity modulated proton therapy (IMPT) plan degradation of mobile lung tumors. Four planning paradigms were compared for nine lung cancer patients. Internal gross tumor volume (IGTV) and internal clinical target volume (ICTV) structures were defined encompassing their respective volumes in every 4DCT phase. The paradigms use different planning CT (pCT) created from the average intensity projection (AIP) of the 4DCT, overwriting the density within the IGTV to account for movement. The density overwrites were: (a) constant filling with 100 HU (C100) or (b) 50 HU (C50), (c) maximum intensity projection (MIP) across phases, and (d) water equivalent path length (WEPL) consideration from beam's-eye-view. Plans were created optimizing dose-influence matrices calculated with fast GPU Monte Carlo (MC) simulations in each pCT. Plans were evaluated with MC on the 4DCTs using a model of the beam delivery time structure. Dose accumulation was performed using deformable image registration. Interplay effect was addressed applying 10 times rescanning. Significantly less DVH metrics degradation occurred when using MIP and WEPL approaches. Target coverage ([Formula: see text] Gy(RBE)) was fulfilled in most cases with MIP and WEPL ([Formula: see text] Gy (RBE)), keeping dose heterogeneity low ([Formula: see text] Gy(RBE)). The mean lung dose was kept lowest by the WEPL strategy, as well as the maximum dose to organs at risk (OARs). The impact on dose levels in the heart, spinal cord and esophagus were patient specific. Overall, the WEPL strategy gives the best performance and should be preferred when using a 3D static geometry for lung cancer IMPT treatment planning. Newly available fast MC methods make it possible to handle long simulations based on 4D data sets to perform studies with high accuracy and efficiency, even prior to individual treatment planning.

Planning strategies for inter-fractional robustness in pancreatic patients treated with scanned carbon therapy

BACKGROUND

Managing inter-fractional anatomy changes is a challenging task in radiotherapy of pancreatic tumors, especially in scanned carbon-ion delivery. This treatment planning study aims to focus on clinically feasible solutions, such as the beam angle selection and margin design to increase the robustness against inter-fractional uncertainties.

METHODS

This study included 10 patients with weekly 3D-CT imaging and physician-approved Clinical Target Volume (CTV). The study was directed to keep the CTV-coverage using six beam angle configurations in combination with different Internal Target Volume (ITV) concepts. These were: geometric-margin (symmetric 3 and 5 mm margin); range-equivalent margins with an isotropic HU replacement; and to evaluate the need of asymmetric margins the water-equivalent range path (WEPL) was determined per patient from the set of CTs. Plan optimization and forward dose calculation in each week-CT were performed with the research treatment planning system TRiP98 and the plan quality evaluated in terms of CTV coverage (V95CTV) and homogeneity dose (HCTV = D5-D95).

RESULTS

The beam geometry had a substantial impact on the target irradiation over the treatment course, with the single posterior or two beams showing the best average coverage of the CTV. The use of geometric margins for the more robust beam geometries showed acceptable results, with a V95CTV of (99.2 ± 1.2)% for the 5 mm-margin. For the non-robust configurations, due to substantial changes in the radiological depth, the use of this margin results in a V95CTV that might be below 80%, only showing improvement when the range changes are included.

CONCLUSIONS

Selection of adequate beam configurations and treatment margins in ion-beam therapy of pancreatic tumors is of great importance. For a single posterior beam or two beam configurations, application of geometrical margins compensate for dose degradation induced by inter-fractional anatomy changes for the majority of the analyzed treatment fractions.

Robust Intensity Modulated Proton Therapy (IMPT) Increases Estimated Clinical Benefit in Head and Neck Cancer Patients

PURPOSE

To compare the clinical benefit of robust optimized Intensity Modulated Proton Therapy (minimax IMPT) with current photon Intensity Modulated Radiation Therapy (IMRT) and PTV-based IMPT for head and neck cancer (HNC) patients. The clinical benefit is quantified in terms of both Normal Tissue Complication Probability (NTCP) and target coverage in the case of setup and range errors.

METHODS AND MATERIALS

For 10 HNC patients, PTV-based IMRT (7 fields), minimax and PTV-based IMPT (2, 3, 4, 5 and 7 fields) plans were tested on robustness. Robust optimized plans differed from PTV-based plans in that they target the CTV and penalize possible error scenarios, instead of using the static isotropic CTV-PTV margin. Perturbed dose distributions of all plans were acquired by simulating in total 8060 setup (±3.5 mm) and range error (±3%) combinations. NTCP models for xerostomia and dysphagia were used to predict the clinical benefit of IMPT versus IMRT.

RESULTS

The robustness criterion was met in the IMRT and minimax IMPT plans in all error scenarios, but this was only the case in 1 of 40 PTV-based IMPT plans. Seven (out of 10) patients had relatively large NTCP reductions in minimax IMPT plans compared to IMRT. For these patients, xerostomia and dysphagia NTCP values were reduced by 17.0% (95% CI; 13.0-21.1) and 8.1% (95% CI; 4.9-11.2) on average with minimax IMPT. Increasing the number of fields did not contribute to plan robustness, but improved organ sparing.

CONCLUSIONS

The estimated clinical benefit in terms of NTCP of robust optimized (minimax) IMPT is greater than that of IMRT and PTV-based IMPT in HNC patients. Furthermore, the target coverage of minimax IMPT plans in the presence of errors was comparable to IMRT plans.

Changes in the radiological depth correlate with dosimetric deterioration in particle therapy for stage I NSCLC patients under high frequency jet ventilation

BACKGROUND

Particle dose distributions are highly sensitive to anatomy changes in the beam path, which may lead to substantial dosimetric deviations. Robust planning and dedicated image guidance together with strategies for online decision making to counteract dosimetric deterioration are thus mandatory. We aimed to develop methods to quantify anatomical discrepancies as depicted by repeated computed tomography (CT) imaging and to test whether they can predict deviations in target coverage.

MATERIAL AND METHODS

Dedicated software tools allowed for voxel-based calculations of changes in the water equivalent path length (WEPL) in beam directions. We prepared proton and carbon ion plans with different coplanar beam angle settings on a series of lung cancer patients, for which planning and localization CT scans under high frequency jet ventilation (HFJV) for tumor fixation were performed. We investigated the reproducibility of target coverage between the optimized and recalculated treatment plans. We then studied how different raster scan and planning settings influence the robustness. Finally, we carried out a systematic analysis of the variations in the WEPL along different coplanar beam angles to find beam directions, which could minimize such variations.

RESULTS

The Spearman's correlations for the GTV ΔV95 and ΔV98 with the ΔWEPL for the proton plans with a 0° and -45° two-field configuration were 0.701 (p = 0.02) and 0.719 (p = 0.08), respectively. For beam configurations 0° and -90°, or 0° and + 45°, with lower ΔWEPL, the correlations were no significant. The same trends were observed for the carbon ion plans. Increased beam spot overlap reduced dosimetric deterioration in case of large ΔWEPL.

CONCLUSION

Software tools for fast online analysis of WEPL changes might help supporting clinical decision making of image guidance. Raster scan and treatment planning settings can help to compensate for anatomical deviations.

Tumour Movement in Proton Therapy: Solutions and Remaining Questions: A Review

Movement of tumours, mostly by respiration, has been a major problem for treating lung cancer, liver tumours and other locations in the abdomen and thorax. Organ motion is indeed one component of geometrical uncertainties that includes delineation and target definition uncertainties, microscopic disease and setup errors. At present, minimising motion seems to be the easiest to implement in clinical practice. If combined with adaptive approaches to correct for gradual anatomical variations, it may be a practical strategy. Other approaches such as repainting and tracking could increase the accuracy of proton therapy delivery, but advanced 4D solutions are needed. Moreover, there is a need to perform clinical studies to investigate which approach is the best in a given clinical situation. The good news is that existing and emerging technology and treatment planning systems as will without doubt lead in the forthcoming future to practical solutions to tackle intra-fraction motion in proton therapy. These developments may also improve motion management in photon therapy as well.

Beam-specific planning volumes for scattered-proton lung radiotherapy

This work describes the clinical implementation of a beam-specific planning treatment volume (bsPTV) calculation for lung cancer proton therapy and its integration into the treatment planning process. Uncertainties incorporated in the calculation of the bsPTV included setup errors, machine delivery variability, breathing effects, inherent proton range uncertainties and combinations of the above. Margins were added for translational and rotational setup errors and breathing motion variability during the course of treatment as well as for their effect on proton range of each treatment field. The effect of breathing motion and deformation on the proton range was calculated from 4D computed tomography data. Range uncertainties were considered taking into account the individual voxel HU uncertainty along each proton beamlet. Beam-specific treatment volumes generated for 12 patients were used: a) as planning targets, b) for routine plan evaluation, c) to aid beam angle selection and d) to create beam-specific margins for organs at risk to insure sparing. The alternative planning technique based on the bsPTVs produced similar target coverage as the conventional proton plans while better sparing the surrounding tissues. Conventional proton plans were evaluated by comparing the dose distributions per beam with the corresponding bsPTV. The bsPTV volume as a function of beam angle revealed some unexpected sources of uncertainty and could help the planner choose more robust beams. Beam-specific planning volume for the spinal cord was used for dose distribution shaping to ensure organ sparing laterally and distally to the beam.

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

PURPOSE

Intensity-modulated proton therapy (IMPT) is highly sensitive to uncertainties in beam range and patient setup. Conventionally, these uncertainties are dealt using geometrically expanded planning target volume (PTV). In this paper, the authors evaluated a robust optimization method that deals with the uncertainties directly during the spot weight optimization to ensure clinical target volume (CTV) coverage without using PTV. The authors compared the two methods for a population of head and neck (H&N) cancer patients.

METHODS

Two sets of IMPT plans were generated for 14 H&N cases, one being PTV-based conventionally optimized and the other CTV-based robustly optimized. For the PTV-based conventionally optimized plans, the uncertainties are accounted for by expanding CTV to PTV via margins and delivering the prescribed dose to PTV. For the CTV-based robustly optimized plans, spot weight optimization was guided to reduce the discrepancy in doses under extreme setup and range uncertainties directly, while delivering the prescribed dose to CTV rather than PTV. For each of these plans, the authors calculated dose distributions under various uncertainty settings. The root-mean-square dose (RMSD) for each voxel was computed and the area under the RMSD-volume histogram curves (AUC) was used to relatively compare plan robustness. Data derived from the dose volume histogram in the worst-case and nominal doses were used to evaluate the plan optimality. Then the plan evaluation metrics were averaged over the 14 cases and were compared with two-sided paired t tests.

RESULTS

CTV-based robust optimization led to more robust (i.e., smaller AUCs) plans for both targets and organs. Under the worst-case scenario and the nominal scenario, CTV-based robustly optimized plans showed better target coverage (i.e., greater D95%), improved dose homogeneity (i.e., smaller D5% - D95%), and lower or equivalent dose to organs at risk.

CONCLUSIONS

CTV-based robust optimization provided significantly more robust dose distributions to targets and organs than PTV-based conventional optimization in H&N using IMPT. Eliminating the use of PTV and planning directly based on CTV provided better or equivalent normal tissue sparing.

Physics controversies in proton therapy

The physical characteristics of proton beams are appealing for cancer therapy. The rapid increase in operational and planned proton therapy facilities may suggest that this technology is a "plug-and-play" valuable addition to the arsenal of the radiation oncologist and medical physicist. In reality, the technology is still evolving, so planning and delivery of proton therapy in patients face many practical challenges. This review article discusses the current status of proton therapy treatment planning and delivery techniques, indicates current limitations in dealing with range uncertainties, and proposes possible developments for proton therapy and supplementary technology to try to realize the actual potential of proton therapy.

The use of proton-beam therapy in the treatment of non-small-cell lung cancer

Lung cancer is the most common cause of cancer death worldwide. Surgical resection has played a major role in the treatment of non-small-cell lung cancer (NSCLC); however, the disease is often detected in a progressive and inoperable form. Surgical resection may also be impossible for early-stage NSCLC due to medical conditions, such as pulmonary or cardiovascular disease and old age. Radiotherapy plays an important role for these patients. Proton-beam therapy is a particle radiotherapy with an excellent dose localization that permits treatment of lung cancer by administering a high dose to the tumor while minimizing damage to the surrounding normal tissues. Thus, proton beams are increasingly being used for lung cancer. In this context, the authors review the current knowledge on proton-beam therapy for the treatment of NSCLC.

Statistical assessment of proton treatment plans under setup and range uncertainties

PURPOSE

To evaluate a method for quantifying the effect of setup errors and range uncertainties on dose distribution and dose-volume histogram using statistical parameters; and to assess existing planning practice in selected treatment sites under setup and range uncertainties.

METHODS AND MATERIALS

Twenty passively scattered proton lung cancer plans, 10 prostate, and 1 brain cancer scanning-beam proton plan(s) were analyzed. To account for the dose under uncertainties, we performed a comprehensive simulation in which the dose was recalculated 600 times per given plan under the influence of random and systematic setup errors and proton range errors. On the basis of simulation results, we determined the probability of dose variations and calculated the expected values and standard deviations of dose-volume histograms. The uncertainties in dose were spatially visualized on the planning CT as a probability map of failure to target coverage or overdose of critical structures.

RESULTS

The expected value of target coverage under the uncertainties was consistently lower than that of the nominal value determined from the clinical target volume coverage without setup error or range uncertainty, with a mean difference of -1.1% (-0.9% for breath-hold), -0.3%, and -2.2% for lung, prostate, and a brain cases, respectively. The organs with most sensitive dose under uncertainties were esophagus and spinal cord for lung, rectum for prostate, and brain stem for brain cancer.

CONCLUSIONS

A clinically feasible robustness plan analysis tool based on direct dose calculation and statistical simulation has been developed. Both the expectation value and standard deviation are useful to evaluate the impact of uncertainties. The existing proton beam planning method used in this institution seems to be adequate in terms of target coverage. However, structures that are small in volume or located near the target area showed greater sensitivity to uncertainties.

Compensation method for respiratory motion in proton treatment planning for mobile liver cancer

We evaluated the dosimetric effect of a respiration motion, and sought an effective planning strategy to compensate the motion using four-dimensional computed tomography (4D CT) dataset of seven selected liver patients. For each patient, we constructed four different proton plans based on: (1) average (AVG) CT, (2) maximum-intensity projection (MIP) CT, (3) AVG CT with density override of tumor volume (OVR), and (4) AVG CT with field-specific proton margin which was determined by the range difference between AVG and MIP plans (mAVG). The overall effectiveness of each planning strategy was evaluated by calculating the cumulative dose distribution over an entire breathing cycle. We observed clear differences between AVG and MIP CT-based plans, with significant underdosages at expiratory and inspiratory phases, respectively. Only the mAVG planning strategy was fully successful as the field-specific proton margin applied in the planning strategy complemented both the limitations of AVG and MIP CT-based strategies. These results demonstrated that respiration motion induced significant changes in dose distribution of 3D proton plans for mobile liver cancer and the changes can be effectively compensated by applying field-specific proton margin to each proton field.

Proton therapy radiation pneumonitis local dose-response in esophagus cancer patients

PURPOSE

This study quantifies pulmonary radiation toxicity in patients who received proton therapy for esophagus cancer.

MATERIALS/METHODS

We retrospectively studied 100 esophagus cancer patients treated with proton therapy. The linearity of the enhanced FDG uptake vs. proton dose was evaluated using the Akaike Information Criterion (AIC). Pneumonitis symptoms (RP) were assessed using the Common Toxicity Criteria for Adverse Events version 4.0 (CTCAEv4). The interaction of the imaging response with dosimetric parameters and symptoms was evaluated.

RESULTS

The RP scores were: 0 grade 4/5, 7 grade 3, 20 grade 2, 37 grade 1, and 36 grade 0. Each dosimetric parameter was significantly higher for the symptomatic group. The AIC winning models were 30 linear, 52 linear quadratic, and 18 linear logarithmic. There was no significant difference in the linear coefficient between models. The slope of the FDG vs. proton dose response was 0.022 for the symptomatic and 0.012 for the asymptomatic (p=0.014). Combining dosimetric parameters with the slope did not improve the sensitivity or accuracy in identifying symptomatic cases.

CONCLUSIONS

The proton radiation dose response on FDG PET/CT imaging exhibited a predominantly linear dose response on modeling. Symptomatic patients had a higher dose response slope.

Proton SBRT for medically inoperable stage I NSCLC

INTRODUCTION

The physical properties of proton beam radiation may offer advantages for treating patients with non-small-cell lung cancer (NSCLC). However, its utility for the treatment of medically inoperable stage I NSCLC patients with stereotactic body radiation therapy (SBRT) is unknown.

METHODS

Outcomes for patients with medically inoperable stage I NSCLC treated with proton SBRT were retrospectively analyzed. Proton SBRT was selected as the treatment modality based on pulmonary comorbidities (n = 5), prior chest radiation or/and multiple primary tumors (n = 7), or other reasons (n = 3). Treatments were administered using 2 to 3 proton beams. Treatment toxicity was scored according to common toxicity criteria for adverse events version 4 criteria.

RESULTS

Fifteen consecutive patients and 20 tumors were treated with proton SBRT to 42 to 50 Gy(relative biological effectiveness) in 3 to 5 fractions between July 2008 and September 2010. Treatments were well tolerated with only one case of grade 2 fatigue, one case of grade 2 dermatitis, three cases of rib fracture (maximum grade 2), and one case of grade 3 pneumonitis in a patient with severe chronic obstructive pulmonary disease. With a median follow-up of 24.1 months, 2-year overall survival and local control rates were 64% (95% confidence limits, 34%-83%) and 100% (83%-100%), respectively.

CONCLUSIONS

We conclude that proton SBRT is effective and well tolerated in this unfavorable group of patients. Prospective clinical trials testing the utility of proton SBRT in stage I NSCLC are warranted.

Fast range-corrected proton dose approximation method using prior dose distribution

For robust plan optimization and evaluation purposes, one needs a computationally efficient way to calculate dose distributions and dose-volume histograms (DVHs) under various changes in the variables associated with beam delivery and images. In this study, we report an approximate method for rapid calculation of dose when setup errors and anatomical changes occur during proton therapy. This fast dose approximation method calculates new dose distributions under various circumstances based on the prior knowledge of dose distribution from a reference setting. In order to validate the method, we calculated and compared the dose distributions from our approximation method to the dose distributions calculated from a clinically commissioned treatment planning system which was used as the ground truth. The overall accuracy of the proposed method was tested against varying degrees of setup error and anatomical deformation for selected patient cases. The setup error was simulated by rigid shifts of the patient; while the anatomical deformation was introduced using weekly acquired repeat CT data sets. We evaluated the agreement between the dose approximation method and full dose recalculation using a 3D gamma index and the root-mean-square (RMS) and maximum deviation of the cumulative dose volume histograms (cDVHs). The average passing rate of 3D gamma analysis under 3% dose and 3 mm distance-to-agreement criteria were 96% and 89% for setup errors and severe anatomy changes, respectively. The average of RMS and maximum deviation of the cDVHs under the setup error was 0.5% and 1.5%, respectively for all structures considered. Similarly, the average of RMS and maximum deviations under the weekly anatomical change were 0.6% and 2.7%, respectively. Our results show that the fast dose approximation method was able to account for the density variation of the patient due to the setup and anatomical changes with acceptable accuracy while significantly improving the computation time.

Results of a multicentric in silico clinical trial (ROCOCO): comparing radiotherapy with photons and protons for non-small cell lung cancer

INTRODUCTION

This multicentric in silico trial compares photon and proton radiotherapy for non-small cell lung cancer patients. The hypothesis is that proton radiotherapy decreases the dose and the volume of irradiated normal tissues even when escalating to the maximum tolerable dose of one or more of the organs at risk (OAR).

METHODS

Twenty-five patients, stage IA-IIIB, were prospectively included. On 4D F18-labeled fluorodeoxyglucose-positron emission tomography-computed tomography scans, the gross tumor, clinical and planning target volumes, and OAR were delineated. Three-dimensional conformal radiotherapy (3DCRT) and intensity-modulated radiotherapy (IMRT) photon and passive scattered conformal proton therapy (PSPT) plans were created to give 70 Gy to the tumor in 35 fractions. Dose (de-)escalation was performed by rescaling to the maximum tolerable dose.

RESULTS

Protons resulted in the lowest dose to the OAR, while keeping the dose to the target at 70 Gy. The integral dose (ID) was higher for 3DCRT (59%) and IMRT (43%) than for PSPT. The mean lung dose reduced from 18.9 Gy for 3DCRT and 16.4 Gy for IMRT to 13.5 Gy for PSPT. For 10 patients, escalation to 87 Gy was possible for all 3 modalities. The mean lung dose and ID were 40 and 65% higher for photons than for protons, respectively.

CONCLUSIONS

The treatment planning results of the Radiation Oncology Collaborative Comparison trial show a reduction of ID and the dose to the OAR when treating with protons instead of photons, even with dose escalation. This shows that PSPT is able to give a high tumor dose, while keeping the OAR dose lower than with the photon modalities.

Status of particle therapy for lung cancer

The current standard treatment for lung cancer, the most common type of cancer worldwide, depends on disease stage. Surgery is the treatment of choice for early-stage tumors, but radiotherapy is a good option for those with early-stage tumors who cannot undergo surgery, and radiotherapy in conjunction with chemotherapy is the standard of care for locally advanced tumors. Although advances in photon (x-ray)-based radiotherapy involving three-dimensional conformal radiotherapy and intensity-modulated radiotherapy allow radiation doses to be escalated beyond the traditional limit of 60 Gy, this dose is not considered to be sufficient for tumor eradication. Moreover, the improvements in local control and survival conferred by concurrent chemotherapy come at the cost of considerable toxicity owing to inadvertent irradiation of surrounding normal tissues, and this toxicity often limits the radiation dose that can be delivered. Unfortunately for patients with locally advanced lung cancer, local control and survival remain poor. Attempts to improve clinical outcomes for patients with lung cancer have led to the use of charged particle therapy in an effort to exploit the physical properties of such particles to escalate the dose to the tumor while simultaneously limiting the dose to nearby structures, thereby enhancing the therapeutic ratio and potentially improving cancer cure rates. This review summarizes the rationale for and challenges associated with the use of charged particles for lung cancer therapy and reviews the clinical experience to date with using protons and carbon ions for early-stage and locally advanced stage non-small cell lung cancer.

Estimated radiation pneumonitis risk after photon versus proton therapy alone or combined with chemotherapy for lung cancer

BACKGROUND

Traditionally, radiation therapy plans are optimized without consideration of chemotherapy. Here, we model the risk of radiation pneumonitis (RP) in the presence of a possible interaction between chemotherapy and radiation dose distribution.

MATERIAL AND METHODS

Three alternative treatment plans are compared in 18 non-small cell lung cancer patients previously treated with helical tomotherapy; the tomotherapy plan, an intensity modulated proton therapy plan (IMPT) and a three dimensional conformal radiotherapy (3D-CRT) plan. All plans are optimized without consideration of the chemotherapy effect. The effect of chemotherapy is modeled as an independent cell killing process using a uniform chemotherapy equivalent radiation dose (CERD) added to the entire organ at risk. We estimate the risk of grade 3 or higher RP (G3RP) using the critical volume model.

RESULTS

The mean risk of clinical G3RP at zero CERD is 5% for tomotherapy (range: 1-18 %) and 14% for 3D-CRT (range 2-49%). When the CERD exceeds 9 Gy, however, the risk of RP with the tomotherapy plans become higher than the 3D-CRT plans. The IMPT plans are less toxic both at zero CERD (mean 2%, range 1-5%) and at CERD = 10 Gy (mean 7%, range 1-28%). Tomotherapy yields a lower risk of RP than 3D-CRT for 17/18 patients at zero CERD, but only for 7/18 patients at CERD = 10 Gy. IMPT gives the lowest risk of all plans for 17/18 patients at zero CERD and for all patients with CERD = 10 Gy.

CONCLUSIONS

The low dose bath from highly conformal photon techniques may become relevant for lung toxicity when radiation is combined with cytotoxic chemotherapy as shown here. Proton therapy allows highly conformal delivery while minimizing the low dose bath potentially interacting with chemotherapy. Thus, intensive drug-radiation combinations could be an interesting indication for selecting patients for proton therapy. It is likely that the IMRT plans would perform better if the CERD was accounted for during optimization, but more clinical data is required to facilitate evidence-based plan optimization in the multi-modality setting.

A beam-specific planning target volume (PTV) design for proton therapy to account for setup and range uncertainties

PURPOSE

To report a method for explicitly designing a planning target volume (PTV) for treatment planning and evaluation in heterogeneous media for passively scattered proton therapy and scanning beam proton therapy using single-field optimization (SFO).

METHODS AND MATERIALS

A beam-specific PTV (bsPTV) for proton beams was derived by ray-tracing and shifting ray lines to account for tissue misalignment in the presence of setup error or organ motion. Range uncertainties resulting from inaccuracies in computed tomography-based range estimation were calculated for proximal and distal surfaces of the target in the beam direction. The bsPTV was then constructed based on local heterogeneity. The bsPTV thus can be used directly as a planning target as if it were in photon therapy. To test the robustness of the bsPTV, we generated a single-field proton plan in a virtual phantom. Intentional setup and range errors were introduced. Dose coverage to the clinical target volume (CTV) under various simulation conditions was compared between plans designed based on the bsPTV and a conventional PTV.

RESULTS

The simulated treatment using the bsPTV design performed significantly better than the plan using the conventional PTV in maintaining dose coverage to the CTV. With conventional PTV plans, the minimum coverage to the CTV dropped from 99% to 67% in the presence of setup error, internal motion, and range uncertainty. However, plans using the bsPTV showed minimal drop of target coverage from 99% to 94%.

CONCLUSIONS

The conventional geometry-based PTV concept used in photon therapy does not work well for proton therapy. We investigated and validated a beam-specific PTV method for designing and evaluating proton plans.

Uncertainties in planned dose due to the limited voxel size of the planning CT when treating lung tumors with proton therapy

Dose calculation for lung tumors can be challenging due to the low density and the fine structure of the geometry. The latter is not fully considered in the CT image resolution used in treatment planning causing the prediction of a more homogeneous tissue distribution. In proton therapy, this could result in predicting an unrealistically sharp distal dose falloff, i.e. an underestimation of the distal dose falloff degradation. The goal of this work was the quantification of such effects. Two computational phantoms resembling a two-dimensional heterogeneous random lung geometry and a swine lung were considered applying a variety of voxel sizes for dose calculation. Monte Carlo simulations were used to compare the dose distributions predicted with the voxel size typically used for the treatment planning procedure with those expected to be delivered using the finest resolution. The results show, for example, distal falloff position differences of up to 4 mm between planned and expected dose at the 90% level for the heterogeneous random lung (assuming treatment plan on a 2 × 2 × 2.5 mm(3) grid). For the swine lung, differences of up to 38 mm were seen when airways are present in the beam path when the treatment plan was done on a 0.8 × 0.8 × 2.4 mm(3) grid. The two-dimensional heterogeneous random lung phantom apparently does not describe the impact of the geometry adequately because of the lack of heterogeneities in the axial direction. The differences observed in the swine lung between planned and expected dose are presumably due to the poor axial resolution of the CT images used in clinical routine. In conclusion, when assigning margins for treatment planning for lung cancer, proton range uncertainties due to the heterogeneous lung geometry and CT image resolution need to be considered.

Investigation of an implantable dosimeter for single-point water equivalent path length verification in proton therapy

PURPOSE

In vivo range verification in proton therapy is highly desirable. A recent study suggested that it was feasible to use point dose measurement for in vivo beam range verification in proton therapy, provided that the spread-out Bragg peak dose distribution is delivered in a different and rather unconventional manner. In this work, the authors investigate the possibility of using a commercial implantable dosimeter with wireless reading for this particular application.

METHODS

The traditional proton treatment technique delivers all the Bragg peaks required for a SOBP field in a single sequence, producing a constant dose plateau across the target volume. As a result, a point dose measurement anywhere in the target volume will produce the same value, thus providing no information regarding the water equivalent path length to the point of measurement. However, the same constant dose distribution can be achieved by splitting the field into a complementary pair of subfields, producing two oppositely "sloped" depth-dose distributions, respectively. The ratio between the two distributions can be a sensitive function of depth and measuring this ratio at a point inside the target volume can provide the water equivalent path length to the dosimeter location. Two types of field splits were used in the experiment, one achieved by the technique of beam current modulation and the other by manipulating the location and width of the beam pulse relative to the range modulator track. Eight MOSFET-based implantable dosimeters at four different depths in a water tank were used to measure the dose ratios for these field pairs. A method was developed to correct the effect of the well-known LET dependence of the MOSFET detectors on the depth-dose distributions using the columnar recombination model. The LET-corrected dose ratios were used to derive the water equivalent path lengths to the dosimeter locations to be compared to physical measurements.

RESULTS

The implantable dosimeters measured the dose ratios with a reasonable relative uncertainty of 1%-3% at all depths, except when the ratio itself becomes very small. In total, 55% of the individual measurements reproduced the water equivalent path lengths to the dosimeters within 1 mm. For three dosimeters, the difference was consistently less than 1 mm. Half of the standard deviations over the repeated measurements were equal or less than 1 mm.

CONCLUSIONS

With a single fitting parameter, the LET-correction method worked remarkably well for the MOSFET detectors. The overall results were very encouraging for a potential method of in vivo beam range verification with millimeter accuracy. This is sufficient accuracy to expand range of clinical applications in which the authors could use the distal fall off of the proton depth dose for tight margins.

The precision of respiratory-gated delivery of synchrotron-based pulsed beam proton therapy

A synchrotron-based proton therapy system operates in a low repetition rate pulsed beam delivery mode. Unlike cyclotron-based beam delivery, there is no guarantee that a synchrotron beam can be delivered effectively or precisely under the respiratory-gated mode. To evaluate the performance of gated synchrotron treatment, we simulated proton beam delivery in the synchrotron-based respiratory-gated mode using realistic patient breathing signals. Parameters used in the simulation were respiratory motion traces (70 traces from 24 patients), respiratory gate levels (10%, 20% and 30% duty cycles at the exhalation phase) and synchrotron magnet excitation cycles (T(cyc)) (fixed T(cyc) mode: 2.7, 3.0-6.0 s and each patient breathing cycle, and variable T(cyc) mode). The simulations were computed according to the breathing trace in which the proton beams were delivered. In the shorter fixed T(cyc) (<4 s), most of the proton beams were delivered uniformly to the target during the entire expiration phase of the respiratory cycle. In the longer fixed T(cyc) (>4 s) and the variable T(cyc) mode, the proton beams were not consistently delivered during the end-expiration phase of the respiratory cycle. However we found that the longer and variable T(cyc) operation modes delivered proton beams more precisely during irregular breathing.

Evaluation of patient residual deviation and its impact on dose distribution for proton radiotherapy

The residual deviations after final patient repositioning based on bony anatomy and the impact of such deviations on the proton dose distributions was investigated. Digitally reconstructed radiographs (DRRs) and kilovoltage (kV) "portal verification" images from 10 patients treated with passively scattered proton radiotherapy was used to estimate the residual deviation. These changes were then applied to the location of isocenter points that, in effect, moved the isocenter relative to the apertures and compensators. A composite verification plan was obtained and compared with the original clinical treatment plan to evaluate any changes in dose distributions. The residual deviations were fitted to a Gaussian distribution with μ = -0.9 ± 0.1 mm and σ = 2.55 ± 0.07 mm. The dose distribution showed under- and overcovered dose spots with complex dose distributions both in the target volumes and in the organs at risk. In some cases, this amounts to 63.5% above the intended clinical plan. Although patient positioning is carefully verified before treatment delivery and setup uncertainties are accounted for by using compensator smearing and aperture margins, a residual shift in a patient's position can considerably affect the dose distribution.

Proton radiotherapy: the biological effect of treating alternating subsets of fields for different treatment fractions

PURPOSE

Common practice in proton radiotherapy is to deliver a subset of all fields in the treatment plan on any given treatment day. We investigate using biological modeling if the resulting variation in daily dose to normal tissues has a relevant detrimental biological effect.

METHODS AND MATERIALS

For four patient groups, the cumulative normalized total dose (NTD) was determined for normal tissues (OARs) of each patient using the clinically delivered fractionation schedule (FS(clin)), and for hypothetical fractionation schedules delivering all fields every day (FS(all)) or only a single field each day (FS(single)). Cumulative three-dimensional NTD distributions were summarized using the generalized equivalent uniform dose (gEUD) model.

RESULTS

For the skull base/cervical spine chordoma group, the largest effect is a 4-Gy increase in gEUD of the chiasm when treating only a subset of fields on any day. For lung cancer and pancreatic cancer patients, the variation in the gEUD of normal tissues is <0.2 Gy. For the prostate group, FS(clin) increases the gEUD of the femoral heads by 9 Gy compared with FS(all). Use of FS(single) resulted in the highest NTD to normal tissues for any patient. FS(all) resulted in an integral NTD to the patient that is on average 5% lower than FS(clin) and 10% lower than FS(single).

CONCLUSION

The effects of field set of the day treatment delivery depend on the tumor site and number of fields treated each day. Modeling these effects may be important for accurate risk assessment.

High-dose proton therapy and carbon-ion therapy for stage I nonsmall cell lung cancer

BACKGROUND

A study was undertaken to evaluate the clinical outcome of particle therapy for stage I nonsmall cell lung cancer (NSCLC).

METHODS

From April 2003 to April 2007, 80 patients with stage I NSCLC were treated with proton therapy or carbon-ion therapy (57 with proton therapy and 23 with carbon-ion therapy) using 3 treatment protocols. In the first protocol, 80 gray equivalents (GyE) of proton therapy was given in 20 fractions, and the second proton therapy protocol used 60 GyE in 10 fractions. For carbon-ion therapy, 52.8 GyE was given in 4 fractions. After achieving promising preliminary results for the first protocol, the authors started to use the second proton therapy protocol to shorten the overall treatment time. Carbon-ion therapy was started in 2005, and thereafter, both proton and carbon-ion therapy plans were made for each patient, and the 1 that appeared superior was adopted. Patient age ranged from 48 to 89 years (median, 76 years). Thirty-seven patients were medically inoperable, and 43 refused surgery. Forty-two patients had T1 tumors, and 38 had T2 tumors.

RESULTS

The median follow-up period for living patients was 35.5 months. For all 80 patients, the 3-year overall survival, cause-specific survival, and local control rates were 75% (IA: 74%; IB: 76%), 86% (IA: 84%; IB: 88%), and 82% (IA: 87%; IB: 77%), respectively. There were no significant differences in treatment results among the 3 protocols. Grade 3 pulmonary toxicity was observed in only 1 patient.

CONCLUSIONS

Proton therapy and carbon-ion therapy are safe and effective for stage I NSCLC. Further investigation of particle therapy for stage I NSCLC is warranted.

Accuracy of an automatic patient-positioning system based on the correlation of two edge images in radiotherapy

We have clinically evaluated the accuracy of an automatic patient-positioning system based on the image correlation of two edge images in radiotherapy. Ninety-six head & neck images from eight patients undergoing proton therapy were compared with a digitally reconstructed radiograph (DRR) of planning CT. Two edge images, a reference image and a test image, were extracted by applying a Canny edge detector algorithm to a DRR and a 2D X-ray image, respectively, of each patient before positioning. In a simulation using a humanoid phantom, performed to verify the effectiveness of the proposed method, no registration errors were observed for given ranges of rotation, pitch, and translation in the x, y, and z directions. For real patients, however, there were discrepancies between the automatic positioning method and manual positioning by physicians or technicians. Using edged head coronal- and sagittal-view images, the average differences in registration between these two methods for the x, y, and z directions were 0.11 cm, 0.09 cm and 0.11 cm, respectively, whereas the maximum discrepancies were 0.34 cm, 0.38 cm, and 0.50 cm, respectively. For rotation and pitch, the average registration errors were 0.95° and 1.00°, respectively, and the maximum errors were 3.6° and 2.3°, respectively. The proposed automatic patient-positioning system based on edge image comparison was relatively accurate for head and neck patients. However, image deformation during treatment may render the automatic method less accurate, since the test image many differ significantly from the reference image.

Breathing interplay effects during proton beam scanning: simulation and statistical analysis

Treatment delivery with active beam scanning in proton radiation therapy introduces the problem of interplay effects when pencil beam motion occurs on a similar time scale as intra-fractional tumor motion. In situations where fractionation may not provide enough repetition to blur the effects of interplay, repeated delivery or 'repainting' of each field several times within a fraction has been suggested. The purpose of this work was to investigate the effectiveness of different repainting strategies in proton beam scanning. To assess the dosimetric impact of interplay effects, we performed a series of simulations considering the following parameters: tumor motion amplitude, breathing period, asymmetry in the motion trajectory for the target and time required to change the beam energy for the delivery system. Several repainting strategies were compared in terms of potential vulnerability to a dose delivery error. Breathing motion perpendicular to the beam direction (representing superior-inferior type tumor motion in patients) was considered and modeled as an asymmetric sine function with a peak-to-peak amplitude of between 10 and 30 mm. The results show that motion effects cause a narrowing of the high-dose profile and widening of the penumbra. The 90% isodose area was reduced significantly when considering a large motion amplitude of 3 cm. The broadening of the penumbra appears to depend only on the amplitude of tumor motion (assuming harmonic motion). The delivered dose exhibits a shift of 10-15% of the tumor amplitude (or 1-5 mm) in the caudal direction due to breathing asymmetry observed for both sin(4)(x) and sin(6)(x) motion. Of the five repainting techniques studied, so-called 'breath sampling' turned out to be most effective in reducing dose errors with a minimal increase in treatment time. In this method, each energy level is repainted at several evenly spaced times within one breathing period. To keep dose delivery errors below 5% while minimizing treatment time, it is recommended that breath sampling repainting be employed using 5-10 paintings per field for an assumed tumor volume of 8.5 x 8.5 x 10 cm(3). For smaller tumor volumes more repaintings will be required, while for larger volumes five repaintings should be sufficient to achieve the required dose accuracy.

On the clinical spatial resolution achievable with protons and heavier charged particle radiotherapy beams

The 'sub-millimetre precision' often claimed to be achievable in protons and light ion beam therapy is analysed using the Monte Carlo code SHIELD-HIT for a broad range of energies. Based on the range of possible values and uncertainties of the mean excitation energy of water and human tissues, as well as of the composition of organs and tissues, it is concluded that precision statements deserve careful reconsideration for treatment planning purposes. It is found that the range of I-values of water stated in ICRU reports 37, 49 and 73 (1984, 1993 and 2005) for the collision stopping power formulae, namely 67 eV, 75 eV and 80 eV, yields a spread of the depth of the Bragg peak of protons and heavier charged particles (carbon ions) of up to 5 or 6 mm, which is also found to be energy dependent due to other energy loss competing interaction mechanisms. The spread is similar in protons and in carbon ions having analogous practical range. Although accurate depth-dose distribution measurements in water can be used at the time of developing empirical dose calculation models, the energy dependence of the spread causes a substantial constraint. In the case of in vivo human tissues, where distribution measurements are not feasible, the problem poses a major limitation. In addition to the spread due to the currently accepted uncertainties of their I-values, a spread of the depth of the Bragg peak due to the varying compositions of soft tissues is also demonstrated, even for cases which could be considered practically identical in clinical practice. For these, the spreads found were similar to those of water or even larger, providing support to international recommendations advising that body-tissue compositions should not be given the standing of physical constants. The results show that it would be necessary to increase the margins of a clinical target volume, even in the case of a water phantom, due to an 'intrinsic basic physics uncertainty', adding to those margins usually considered in normal clinical practice due to anatomical or therapeutic strategy reasons. Individualized patient determination of tissue composition along the complete beam path, rather than CT Hounsfield numbers alone, would also probably be required even to reach 'sub-centimetre precision'.

A point dose method for in vivo range verification in proton therapy

Range uncertainty in proton therapy is a recognized concern. For certain treatment sites, less optimal beam directions are used to avoid the potential risk, but also with reduced benefit. In vivo dosimetry, with implanted or intra-cavity dosimeters, has been widely used for treatment verification in photon/electron therapy. The method cannot, however, verify the beam range for proton treatment, unless we deliver the treatment in a different manner. Specifically, we split the spread-out Bragg peaks in a proton field into two separate fields, each delivering a 'sloped' depth-dose distribution, rather than the usual plateau in a typical proton field. The two fields are 'sloped' in opposite directions so that the total depth-dose distribution retains the constant dose plateau covering the target volume. By measuring the doses received from both fields and calculating the ratio, the water-equivalent path length to the location of the implanted dosimeter can be verified, thus limiting range uncertainty to only the remaining part of the beam path. Production of such subfields has been experimented with a passive scattering beam delivery system. Phantom measurements have been performed to illustrate the application for in vivo beam range verification.

Characteristics of movement-induced dose reduction in target volume: a comparison between photon and proton beam treatment

We compared the main characteristics of movement-induced dose reduction during photon and proton beam treatment, based on an analysis of dose-volume histograms. To simulate target movement, a target contour was delineated in a scanned phantom and displaced by 3 to 20 mm. Although the dose reductions to the target in the 2 treatment systems were similar for transverse (perpendicular to beam direction) target motion, they were completely different for longitudinal (parallel to beam direction) target motion. While both modalities showed a relationship between the degree of target shift and the reduction in dose coverage, dose reduction showed a strong directional dependence in proton beam treatment. Clinical simulation of target movement for a prostate cancer patient showed that, although coverage and conformity indices for a 6-mm lateral movement of the prostate were reduced by 9% and 16%, respectively, for proton beam treatment, they were reduced by only 1% and 7%, respectively, for photon treatment. This difference was greater for a 15-mm target movement in the lateral direction, which lowered the coverage and conformity indices by 34% and 54%, respectively, for proton beam treatment, but changed little during photon treatment. In addition, we found that the equivalent uniform dose (EUD) and homogeneity index show similar characteristics during target movement. These results suggest that movement-induced dose reduction differs significantly between photon and proton beam treatment. Attention should be paid to the target margin in proton beam treatment due to the distinct characteristics of heavy ion beams.