Automated VMAT treatment planning using sequential convex programming: algorithm development and clinical implementation.

Objective.To develop and clinically implement a fully automated treatment planning system (TPS) for volumetric modulated arc therapy (VMAT).Approach.We solve two constrained optimization problems sequentially. The tumor coverage is maximized at the first step while respecting all maximum/mean dose clinical criteria. The second step further reduces the dose at the surrounding organs-at-risk as much as possible. Our algorithm optimizes the machine parameters (leaf positions and monitor units) directly and the resulting mathematical non-convexity is handled using thesequential convex programmingby solving a series of convex approximation problems. We directly integrate two novel convex surrogate metrics to improve plan delivery efficiency and reduce plan complexity by promoting aperture shape regularity and neighboring aperture similarity. The entire workflow is automated using the Eclipse TPS application program interface scripting and provided to users as a plug-in, requiring the users to solely provide the contours and their preferred arcs. Our program provides the optimal machine parameters and does not utilize the Eclipse optimization engine, however, it utilizes the Eclipse final dose calculation engine. We have tested our program on 60 patients of different disease sites and prescriptions for stereotactic body radiotherapy (paraspinal (24 Gy × 1, 9 Gy × 3), oligometastis (9 Gy × 3), lung (18 Gy × 3, 12 Gy × 4)) and retrospectively compared the automated plans with the manual plans used for treatment. The program is currently deployed in our clinic and being used in our daily clinical routine to treat patients.Main results.The automated plans found dosimetrically comparable or superior to the manual plans. For paraspinal (24 Gy × 1), the automated plans especially improved tumor coverage (the average PTV (Planning Target Volume) 95% from 96% to 98% and CTV100% from 95% to 97%) and homogeneity (the average PTV maximum dose from 120% to 116%). For other sites/prescriptions, the automated plans especially improved the duty cycle (23%-39.4%).Significance.This work proposes a fully automated approach to the mathematically challenging VMAT problem. It also shows how the capabilities of the existing (Food and Drug Administration)FDA-approved commercial TPS can be enhanced using an in-house developed optimization algorithm that completely replaces the TPS optimization engine. The code and pertained models along with a sample dataset will be released on our ECHO-VMAT GitHub (https://github.com/PortPy-Project/ECHO-VMAT).

Automated and Clinically Optimal Treatment Planning for Cancer Radiotherapy.

Each year, approximately 18 million new cancer cases are diagnosed worldwide, and about half must be treated with radiotherapy. A successful treatment requires treatment planning with the customization of penetrating radiation beams to sterilize cancerous cells without harming nearby normal organs and tissues. This process currently involves extensive manual tuning of parameters by an expert planner, making it a time-consuming and labor-intensive process, with quality and immediacy of critical care dependent on the planner's expertise. To improve the speed, quality, and availability of this highly specialized care, Memorial Sloan Kettering Cancer Center developed and applied advanced optimization tools to this problem (e.g., using hierarchical constrained optimization, convex approximations, and Lagrangian methods). This resulted in both a greatly improved radiotherapy treatment planning process and the generation of reliable and consistent high-quality plans that reflect clinical priorities. These improved techniques have been the foundation of high-quality treatments and have positively impacted over 4,000 patients to date, including numerous patients in severe pain and in urgent need of treatment who might have otherwise required longer hospital stays or undergone unnecessary surgery to control the progression of their disease. We expect that the wide distribution of the system we developed will ultimately impact patient care more broadly, including in resource-constrained countries.

Dose to Highly Functional Ventilation Zones Improves Prediction of Radiation Pneumonitis for Proton and Photon Lung Cancer Radiation Therapy.

PURPOSE: We hypothesized that the radiation dose in high-ventilation portions of the lung better predicts radiation pneumonitis (RP) outcome for patients treated with proton radiation therapy (PR) and photon radiation therapy (PH). METHODS AND MATERIALS: Seventy-four patients (38 protons, 36 photons) with locally advanced non-small cell lung cancer treated with concurrent chemoradiation therapy were identified, of whom 24 exhibited RP (graded using Common Terminology Criteria for Adverse Events v4.0) after PR or PH, and 50 were negative controls. The inhale and exhale simulation computed tomography scans were deformed using Advanced Normalization Tools. The 3-dimensional lung ventilation maps were derived from the deformation matrix and partitioned into low- and high-ventilation zones for dosimetric analysis. Receiver operating curve analysis was used to study the power of relationship between RP and ventilation zones to determine an optimal ventilation cutoff. Univariate logistic regression was used to correlate dose in high- and low-ventilation zones with risk of RP. A nonparametric random forest process was used for multivariate importance assessment. RESULTS: The optimal high-ventilation zone definition was determined to be the higher 45% to 60% of the ventilation values. The parameter vV20Gy_high (high ventilation volume receiving ≥20 Gy) was found to be a significant indicator for RP (PH: P = .002, PR: P = .035) with improved areas under the curve compared with the traditional V20Gy for both photon and proton cohorts. The relationship of RP with dose to the low-ventilation zone of the lung was insignificant (PH: P = .123, PR: P = .661). Similar trends were observed for ventilation mean lung dose and ventilation V5Gy. Multivariate importance assessment determined that vV20Gy_high, vV5_high, and mean lung dose were the most significant parameters for the proton cohort with a combined area under the curve of 0.78. CONCLUSION: Dose to the high-ventilated regions of the lung can improve predictions of RP for both PH and PR.

A hybrid proton and hyperpolarized gas tagging MRI technique for lung respiratory motion imaging: a feasibility study.

The aim of this work was to develop a novel hybrid 3D hyperpolarized (HP) gas tagging MRI (t-MRI) technique and to evaluate it for lung respiratory motion measurement with comparison to deformable image registrations (DIR) methods. Three healthy subjects underwent a hybrid MRI which combines 3D HP gas t-MRI with a low resolution (Low-R, 4.5 mm isotropic voxels) 3D proton MRI (p-MRI), plus a high resolution (High-R, 2.5 mm isotropic voxels) 3D p-MRI, during breath-holds at the end-of-inhalation (EOI) and the end-of-exhalation (EOE). Displacement vector field (DVF) of the lung motion was determined from the t-MRI images by tracking tagging grids and from the High-R p-MRI using three DIR methods (B-spline based method implemented by Velocity, Free Form Deformation by MIM, and B-spline by an open source software Elastix: denoted as A, B, and C, respectively), labeled as tDVF and dDVF, respectively. The tDVF from the HP gas t-MRI was used as ground-truth reference to evaluate performance of the three DIR methods. Differences in both magnitude and angle between the tDVF and dDVFs were analyzed. The mean lung motion of the three subjects was 37.3 mm, 8.9 mm and 12.9 mm, respectively. Relatively large discrepancies were observed between the tDVF and the dDVFs as compared to previously reported DIR errors. The mean ± standard deviation (SD) DVF magnitude difference was 8.3 ± 5.6 mm, 9.2 ± 4.5 mm, and 9.3 ± 6.1 mm, and the mean ± SD DVF angular difference was 29.1 ± 12.1°, 50.1 ± 28.6°, and 39.0 ± 6.3°, for the DIR Methods A, B, and C, respectively. These preliminary results showed that the hybrid HP gas t-MRI technique revealed different lung motion patterns as compared to the DIR methods. It may provide unique perspectives in developing and evaluating DIR of the lungs. Novelty and Significance We designed a MRI protocol that includes a novel hybrid MRI technique (3D HP gas t-MRI with a low resolution 3D p-MRI) plus a high resolution 3D p-MRI. We tested the novel hybrid MRI technique on three healthy subjects for measuring regional lung respiratory motion with comparison to deformable image registrations (DIR) methods, and observed relatively large discrepancies in lung motion between HP gas t-MRI and DIR methods.

An initial investigation of hyperpolarized gas tagging magnetic resonance imaging in evaluating deformable image registration-based lung ventilation.

BACKGROUND: Deformable image registration (DIR)-based lung ventilation mapping is attractive due to its simplicity, and also challenging due to its susceptibility to errors and uncertainties. In this study, we explored the use of 3D Hyperpolarized (HP) gas tagging MRI to evaluate DIR-based lung ventilation. METHOD AND MATERIAL: Three healthy volunteers included in this study underwent both 3D HP gas tagging MRI (t-MRI) and 3D proton MRI (p-MRI) using balanced steady-state free precession pulse sequence at end of inhalation and end of exhalation. We first obtained the reference displacement vector fields (DVFs) from the t-MRIs by tracking the motion of each tagging grid between the exhalation and the inhalation phases. Then, we determined DIR-based DVFs from the p-MRIs by registering the images at the two phases with two commercial DIR algorithms. Lung ventilations were calculated from both the reference DVFs and the DIR-based DVFs using the Jacobian method and then compared using cross correlation and mutual information. RESULTS: The DIR-based lung ventilations calculated using p-MRI varied considerably from the reference lung ventilations based on t-MRI among all three subjects. The lung ventilations generated using Velocity AI were preferable for the better spatial homogeneity and accuracy compared to the ones using MIM, with higher average cross correlation (0.328 vs 0.262) and larger average mutual information (0.528 vs 0.323). CONCLUSION: We demonstrated that different DIR algorithms resulted in different lung ventilation maps due to underlining differences in the DVFs. HP gas tagging MRI provides a unique platform for evaluating DIR-based lung ventilation.

Dosimetric feasibility of 4DCT-ventilation imaging guided proton therapy for locally advanced non-small-cell lung cancer.

BACKGROUND: The principle aim of this study is to incorporate 4DCT ventilation imaging into functional treatment planning that preserves high-functioning lung with both double scattering and scanning beam techniques in proton therapy. METHODS: Eight patients with locally advanced non-small-cell lung cancer were included in this study. Deformable image registration was performed for each patient on their planning 4DCTs and the resultant displacement vector field with Jacobian analysis was used to identify the high-, medium- and low-functional lung regions. Five plans were designed for each patient: a regular photon IMRT vs. anatomic proton plans without consideration of functional ventilation information using double scattering proton therapy (DSPT) and intensity modulated proton therapy (IMPT) vs. functional proton plans with avoidance of high-functional lung using both DSPT and IMPT. Dosimetric parameters were compared in terms of tumor coverage, plan heterogeneity, and avoidance of normal tissues. RESULTS: Our results showed that both DSPT and IMPT plans gave superior dose advantage to photon IMRTs in sparing low dose regions of the total lung in terms of V5 (volume receiving 5Gy). The functional DSPT only showed marginal benefit in sparing high-functioning lung in terms of V5 or V20 (volume receiving 20Gy) compared to anatomical plans. Yet, the functional planning in IMPT delivery, can further reduce the low dose in high-functioning lung without degrading the PTV dosimetric coverages, compared to anatomical proton planning. Although the doses to some critical organs might increase during functional planning, the necessary constraints were all met. CONCLUSIONS: Incorporating 4DCT ventilation imaging into functional proton therapy is feasible. The functional proton plans, in intensity modulated proton delivery, are effective to further preserve high-functioning lung regions without degrading the PTV coverage.