Dose delivery reproducibility for PBS proton treatment of breast cancer patients with and without mask immobilization.

BACKGROUND: Setup reproducibility of the tissue in the proton beam path is critical in maintaining the planned clinical target volume (CTV) dose coverage and sparing the organs at risk (OAR). In this study, we retrospectively evaluated radiation therapy dose reproducibility for proton pencil beam scanning (PBS) treatment of breast cancer patients with and without mask immobilization. METHODS: Ninety-four patients treated between January 2019 and September 2022 with at least one verification CT scan (V-CT) in treatment position were included for this study. All patients were set up with arms up using the Orfit AIO patient positioning system, with (69 patients) or without (25 patients) mask immobilization in chin, neck, shoulder, upper arm, and chest areas. Two to three enface or near enface single field uniform dose PBS beams were optimized using a commercial treatment planning system. Prescription doses were 25 to 60 GyRBE in 5 to 45 fractions. Treatment plan doses re-calculated on V-CTs were compared to the corresponding planned doses. Cumulative doses were also calculated for patients with at least 3 V-CTs by deform and weighted sum doses from V-CTs to corresponding P-CTs. CTV D95%, ipsilateral-lung V40%, esophagus D0.01cc, and heart mean dose were evaluated and reported as percentages of prescription doses. Differences were large dose deteriorations (LDD) if: (1) CTV (V-CT/cumulative D95%) - (Planned D95%) < - 5%; or (2) Ipsilateral-lung (V-CT/cumulative V40%) - (Planned V40%) > 5%; or (3) Esophagus (V-CT/cumulative D0.01cc) - (Planned D0.01cc) > 10%; or (4) Heart (V-CT/cumulative mean) - (Planned mean) > 1.5%. RESULTS: On average, V-CT/cumulative and planned CTV/OAR dose parameter differences were less than 2.2%/1.7% and 3.4%/3.7% for masked and maskless patients, respectively. The percentages of patients with at least one CTV or OAR V-CT/cumulative dose LDD were 20.3%/25.0% and 72.0%/54.0% for masked and maskless patients, respectively. CONCLUSIONS: On average, masked/maskless setups achieved delivered and planned CTV/OAR dose parameters agreed within 2.2%/3.7% for PBS treatment of breast cancer patients in this study. Maskless patients had higher rate of CTV/OAR LDDs compared to masked patients. Dosimetric differences large enough to raise clinical concerns in either group were able to be addressed with replannings.

Virtual particle Monte Carlo: A new concept to avoid simulating secondary particles in proton therapy dose calculation.

BACKGROUND: In proton therapy dose calculation, Monte Carlo (MC) simulations are superior in accuracy but more time consuming, compared to analytical calculations. Graphic processing units (GPUs) are effective in accelerating MC simulations but may suffer thread divergence and racing condition in GPU threads that degrades the computing performance due to the generation of secondary particles during nuclear reactions. PURPOSE: A novel concept of virtual particle (VP) MC (VPMC) is proposed to avoid simulating secondary particles in GPU-accelerated proton MC dose calculation and take full advantage of the computing power of GPU. METHODS: Neutrons and gamma rays were ignored as escaping from the human body; doses of electrons, heavy ions, and nuclear fragments were locally deposited; the tracks of deuterons were converted into tracks of protons. These particles, together with primary and secondary protons, are considered to be the realistic particles. Histories of primary and secondary protons were replaced by histories of multiple VPs. Each VP corresponded to one proton (either primary or secondary). A continuous-slowing-down-approximation model, an ionization model, and a large angle scattering event model corresponding to nuclear interactions were developed for VPs by generating probability distribution functions (PDFs) based on simulation results of realistic particles using MCsquare. For efficient calculations, these PDFs were stored in the Compute Unified Device Architecture textures. VPMC was benchmarked with TOPAS and MCsquare in phantoms and with MCsquare in 13 representative patient geometries. Comparisons between the VPMC calculated dose and dose measured in water during patient-specific quality assurance (PSQA) of the selected 13 patients were also carried out. Gamma analysis was used to compare the doses derived from different methods and calculation efficiencies were also compared. RESULTS: Integrated depth dose and lateral dose profiles in both homogeneous and inhomogeneous phantoms all matched well among VPMC, TOPAS, and MCsquare calculations. The 3D-3D gamma passing rates with a criterion of 2%/2 mm and a threshold of 10% was 98.49% between MCsquare and TOPAS and 98.31% between VPMC and TOPAS in homogeneous phantoms, and 99.18% between MCsquare and TOPAS and 98.49% between VPMC and TOPAS in inhomogeneous phantoms, respectively. In patient geometries, the 3D-3D gamma passing rates with 2%/2 mm/10% between dose distributions from VPMC and MCsquare were 98.56 ± 1.09% in patient geometries. The 2D-3D gamma analysis with 3%/2 mm/10% between the VPMC calculated dose distributions and the 2D measured planar dose distributions during PSQA was 98.91 ± 0.88%. VPMC calculation was highly efficient and took 2.84 ± 2.44 s to finish for the selected 13 patients running on four NVIDIA Ampere GPUs in patient geometries. CONCLUSION: VPMC was found to achieve high accuracy and efficiency in proton therapy dose calculation.

Exploratory study of seed spots analysis to characterize dose and linear-energy-transfer effect in adverse event initialization of pencil-beam-scanning proton therapy.

BACKGROUND: Both dose and linear energy transfer (LET) could play a substantial role in adverse event (AE) initialization of cancer patients treated with pencil-beam-scanning (PBS) proton therapy. However, not all the voxels within the AE regions are directly induced from the dose and LET effect. It is important to study the synergistic effect of dose and LET in AE initialization by only including a subset of voxels that are dosimetrically important. PURPOSE: To perform exploratory investigation of the dose and LET effects upon AE initialization in PBS using seed spots analysis. METHODS: A total of 113 head-and-neck (H&N) cancer patients receiving curative PBS were included. Among them, 20 patients experienced unanticipated CTCAEv4.0 grade ≥3 AEs (AE group) and 93 patients did not (control group). Within the AE group, 13 AE patients were included in the seed spot analysis to derive the descriptive features of AE initialization and the remaining 7 mandible osteoradionecrosis patients and 93 control patients were used to derive the feature-based volume constraint of mandible osteoradionecrosis. The AE regions were contoured and the corresponding dose-LET volume histograms of AE regions were generated for all patients in the AE group. We selected high LET voxels (the highest 5% of each dose bin) with a range of moderate-to-high dose (≥∼40-Gy relative biological effectiveness) as critical voxels. Critical voxels that were contiguous with each other were grouped into clusters. Each cluster was considered a potential independent seed spot for AE initialization. Seed spots were displayed in a 2D dose-LET plane based on their mean dose and LET to derive the descriptive features of AE initialization. A volume constraint of mandible osteoradionecrosis was then established based on the extracted features using a receiver operating characteristic curve. RESULTS: The product of dose and LET (xBD) was found to be a descriptive feature of seed spots leading to AE initialization in this preliminary study. The derived xBD volume constraint for mandible osteoradionecrosis showed good performance with an area under curve of 0.87 (sensitivity of 0.714 and specificity of 0.807 in the leave-one-out cross-validation) for the very limited patient data included in this study. CONCLUSION: Our exploratory study showed that both dose and LET were observed to be important in AE initializations. The derived xBD volume constraint could predict mandible osteoradionecrosis reasonably well in the very limited H&N cancer patient data treated with PBS included in this study.

Comprehensive Commissioning and Clinical Implementation of GammaTiles STaRT for Intracranial Brain Cancer.

Purpose: To validate the dose calculation accuracy and dose distribution of GammaTiles for brain tumors, and to suggest a surgically targeted radiation therapy (STaRT) workflow for planning, delivery, radiation safety documentation, and posttreatment validation. Methods and Materials: Novel surgically targeted radiation therapy, GammaTiles, uses Cs-131 radiation isotopes embedded in collagen-based tiles that can be resorbed after surgery. GammaTile target delineation and dose calculation were performed on MIM Symphony software. Point-based and complex seed distribution calculations in MIM Symphony were verified with hand calculations and BrachyVision calculations. Vendor-provided 2-dimensional dose distribution calculation accuracy was validated using gafchromic EBT3 film measurements at various depths. A workflow was established for safe and effective GammaTile implants. Results: Good agreement was observed between different calculations. Calculation accuracy of less than 0.5% was achieved for all points except one, which had rounding issues for very low doses and resulted in just below 5% difference. Differences in anisotropy and geometry positioning were noticed in the delineation of Cs-131 IsoRay seeds in the compared systems, resulting in minor discrepancies in the calculated dosimetry distributions. Film measurements showed profiles with relatively good agreement of 0% to 5% in nongradient regions with higher differences between 5% to 10% in the sharp dose fall-off regions. Conclusions: A comprehensive evaluation of GammaTile geometry, dose distribution, and clinical workflow was conducted. Safe intro-operative implantation of GammaTiles requires extensive preplanning and interdisciplinary collaboration. A STaRT workflow was outlined to provide a guideline for an accurate treatment planning and safe implant process at other institutions.

GPU-accelerated Monte Carlo-based online adaptive proton therapy: A feasibility study.

PURPOSE: To develop an online graphic processing unit (GPU)-accelerated Monte Carlo-based adaptive radiation therapy (ART) workflow for pencil beam scanning (PBS) proton therapy to address interfraction anatomical changes in patients treated with PBS. METHODS AND MATERIALS: A four-step workflow was developed using our in-house developed GPU-accelerated Monte Carlo-based treatment planning system to implement online Monte Carlo-based ART for PBS. The first step conducts diffeomorphic demon-based deformable image registration (DIR) to propagate contours on the initial planning CT (pCT) to the verification CT (vCT) to form a new structure set. The second step performs forward dose calculation of the initial plan on the vCT with the propagated contours after manual approval (possible modifications involved). The third step triggers a reoptimization of the plan depending on whether the verification dose meets the clinical requirements or not. A robust evaluation will be done for both the verification plan in the second step and the reopotimized plan in the third step. The fourth step involves a two-stage (before and after delivery) patient-specific quality assurance (PSQA) of the reoptimized plan. The before-delivery PSQA is to compare the plan dose to the dose calculated using an independent fast open-source Monte Carlo code, MCsquare. The after-delivery PSQA is to compare the plan dose to the dose recalculated using the log file (spot MU, spot position, and spot energy) collected during the delivery. Jaccard index (JI), dice similarity coefficients (DSCs), and Hausdorff distance (HD) were used to assess the quality of the propagated contours in the first step. A commercial plan evaluation software, ClearCheck™, was integrated into the workflow to carry out efficient plan evaluation. 3D Gamma analysis was used during the fourth step to ensure the accuracy of the plan dose from reoptimization. Three patients with three different disease sites were chosen to evaluate the feasibility of the online ART workflow for PBS. RESULTS: For all three patients, the propagated contours were found to have good volume conformance [JI (lowest-highest: 0.833-0.983) and DSC (0.909-0.992)] but suboptimal boundary coincidence [HD (2.37-20.76 mm)] for organs-at-risk. The verification dose evaluated by ClearCheck™ showed significant degradation of the target coverage due to the interfractional anatomical changes. Reoptimization on the vCT resulted in great improvement of the plan quality to a clinically acceptable level. 3D Gamma analyses of PSQA confirmed the accuracy of the plan dose before delivery (mean Gamma index = 98.74% with a threshold of 2%/2 mm/10%), and after delivery based on the log files (mean Gamma index = 99.05% with a threshold of 2%/2 mm/10%). The average time cost for the complete execution of the workflow was around 858 s, excluding the time for manual intervention. CONCLUSION: The proposed online ART workflow for PBS was demonstrated to be efficient and effective by generating a reoptimized plan that significantly improved the plan quality.

Empirical Relative Biological Effectiveness (RBE) for Mandible Osteoradionecrosis (ORN) in Head and Neck Cancer Patients Treated With Pencil-Beam-Scanning Proton Therapy (PBSPT): A Retrospective, Case-Matched Cohort Study.

Purpose: To retrospectively investigate empirical relative biological effectiveness (RBE) for mandible osteoradionecrosis (ORN) in head and neck (H&N) cancer patients treated with pencil-beam-scanning proton therapy (PBSPT). Methods: We included 1,266 H&N cancer patients, of which, 931 patients were treated with volumetric-modulated arc therapy (VMAT) and 335 were treated with PBSPT. Among them, 26 VMAT and 9 PBSPT patients experienced mandible ORN (ORN group), while all others were included in the control group. To minimize the impact of the possible imbalance in clinical factors between VMAT and PBSPT patients in the dosimetric comparison between these two modalities and the resulting RBE quantification, we formed a 1:1 case-matched patient cohort (335 VMAT patients and 335 PBSPT patients including both the ORN and control groups) using the greedy nearest neighbor matching of propensity scores. Mandible dosimetric metrics were extracted from the case-matched patient cohort and statistically tested to evaluate the association with mandibular ORN to derive dose volume constraints (DVCs) for VMAT and PBSPT, respectively. We sought the equivalent constraint doses for VMAT so that the critical volumes of VMAT were equal to those of PBSPT at different physical doses. Empirical RBEs of PBSPT for ORN were obtained by calculating the ratio between the derived equivalent constraint doses and physical doses of PBSPT. Bootstrapping was further used to get the confidence intervals. Results: Clinical variables of age, gender, tumor stage, prescription dose, chemotherapy, hypertension or diabetes, dental extraction, smoking history, or current smoker were not statistically related to the incidence of ORN in the overall patient cohort. Smoking history was found to be significantly associated with the ORN incidence in PBSPT patients only. V40Gy[RBE], V50Gy[RBE], and V60Gy[RBE] were statistically different (p<0.05) between the ORN and control group for VMAT and PBSPT. Empirical RBEs of 1.58(95%CI: 1.34-1.64), 1.34(95%CI: 1.23-1.40), and 1.24(95%: 1.15-1.26) were obtained for proton dose at 40 Gy[RBE=1.1], 50 Gy[RBE=1.1] and 60 Gy[RBE=1.1], respectively. Conclusions: Our study suggested that RBEs were larger than 1.1 at moderate doses (between 40 and 60 Gy[RBE=1.1]) with high LET for mandible ORN. RBEs are underestimated in current clinical practice in PBSPT. The derived DVCs can be used for PBSPT plan evaluation and optimization to minimize the incidence rate of mandible ORN.

Implementation of Photon Treatment Back-up Workflow at a High-Volume Proton Center: Safety, Quality, and Patient Considerations.

PURPOSE: A successful proton beam therapy (PBT) center relies heavily on the proper function and maintenance of a proton beam therapy machine. However, when a PBT machine is non-operational, a proton facility is hindered with delays that can potentially lead to inferior treatment outcome due to treatment interruption. This article reports a viable solution for a photon back-up plan in a proton down event. METHODS AND MATERIALS: The implementation of a workflow for which proton plans are converted to photon plans so that patients can be treated using photons has been a successful strategy to reduce delays and mitigate its effect on patient care. This workflow was established through collaboration of physicians, physicists, dosimetrists, therapists, nurses, and schedulers. RESULTS AND CONCLUSIONS: A tiered system established by disease site, number of fractions, and individualized circumstances is used to prioritize patients. Proton-photon backup planning strategy and physics check details were described. This article provides an overview of workflow of conversion of PBT to photon when the PBT machine is down. Specific needs of patients undergoing proton beam therapy are addressed.

Per-voxel constraints to minimize hot spots in linear energy transfer-guided robust optimization for base of skull head and neck cancer patients in IMPT.

PURPOSE: Due to the employment of quadratic programming using soft constraints to implement dose volume constraints and the "trial-and-error" procedure needed to achieve a clinically acceptable plan, conventional dose volume constraints (upper limit) are not adequately effective in controlling small and isolated hot spots in the dose/linear energy transfer (LET) distribution. Such hot spots can lead to adverse events. In order to mitigate the risk of brain necrosis, one of the most clinically significant adverse events in patients receiving intensity-modulated proton therapy (IMPT) for base of skull (BOS) cancer, we propose per-voxel constraints to minimize hot spots in LET-guided robust optimization. METHODS AND MATERIALS: Ten BOS cancer patients treated with IMPT were carefully selected by meeting one of the following conditions: (1) diagnosis of brain necrosis during follow-up; and (2) considered high risk for brain necrosis by not meeting dose constraints to the brain. An optimizing structure (BrainOPT) and an evaluating structure (BrainROI) that both contained the aforementioned hot dose regions in the brain were generated for optimization and evaluation, respectively. Two plans were generated for every patient: one using conventional dose-only robust optimization, the other using LET-guided robust optimization. The impact of LET was integrated into the optimization via a term of extra biological dose (xBD). A novel optimization tool of per-voxel constraints to control small and isolated hot spots in either the dose, LET, or combined (dose/LET) distribution was developed and used to minimize dose/LET hot spots of the selected structures. Indices from dose-volume histogram (DVH) and xBD dose-volume histogram (xBDVH) were used in the plan evaluation. A newly developed tool of the dose-LET-volume histogram (DLVH) was also adopted to illustrate the underlying mechanism. Wilcoxon signed-rank test was used for statistical comparison of the DVH and xBDVH indices between the conventional dose-only and the LET-guided robustly optimized plans. RESULTS: Per-voxel constraints effectively and efficiently minimized dose hot spots in both dose-only and LET-guided robust optimization and LET hot spots in LET-guided robust optimization. Compared to the conventional dose-only robust optimization, the LET-guided robust optimization could generate plans with statistically lower xBD hot spots in BrainROI (VxBD,50 Gy[RBE], p = 0.009; VxBD,60 Gy[RBE], p = 0.025; xBD1cc, p = 0.017; xBD2cc, p = 0.022) with comparable dose coverage, dose hot spots in the target, and dose hot spots in BrainROI. DLVH analysis indicated that LET-guided robust optimization could either reduce LET at the same dose level or redistribute high LET from high dose regions to low dose regions. CONCLUSION: Per-voxel constraint is a powerful tool to minimize dose/LET hot spots in IMPT. The LET-guided robustly optimized plans outperformed the conventional dose-only robustly optimized plans in terms of xBD hot spots control.

Using Novel Statistical Techniques to Accurately Determine the Predictive Dose Range in a Study of Overall Survival after Definitive Radiotherapy for Stage III Non-Small Cell Lung Cancer in Association with Heart Dose.

PURPOSE: Recent studies of radiotherapy (RT) for stage III non-small-cell lung cancer (NSCLC) have associated high dose to the heart with cardiac toxicity and decreased overall survival (OS). We used advanced statistical techniques to account for correlations between dosimetric variables and more accurately determine the range of heart doses which are associated with reduced OS in patients receiving RT for stage III NSCLC. METHODS: From 2006 to 2013, 119 patients with stage III NSCLC received definitive RT at our institution. OS data was obtained from institutional tumor registry. We used multivariate Cox model to determine patient specific covariates predictive for reduced overall survival. We examined age, prescription dose, mean lung dose, lung V20, RT technique, stage, chemotherapy, tumor laterality, tumor volume, and tumor site as candidate covariates. We subsequently used novel statistical techniques within multivariate Cox model to systematically search the whole heart dose-volume histogram (DVH) for dose parameters associated with OS. RESULTS: Patients were followed until death or 2.5 to 81.2 months (median 30.4 months) in those alive at last follow up. On multivariate analysis of whole heart DVH, the dose of 51 Gy was identified as a threshold dose above which the dose volume relationship becomes predictive for OS. We identified V55Gy (percentage of the whole heart volume receiving at least 55 Gy) as the best single DVH index which can be used to set treatment optimization constraints (Hazard Ratio = 1.044 per 1% increase in heart volume exposed to at least 55 Gy, P = 0.03). Additional characteristics correlated with OS on multivariate analysis were age, stage (IIIA/IIIB), and administration of chemotherapy. CONCLUSION: Doses above 51 Gy, applied to small volumes of the heart, are associated with worse OS in stage III NSCLC patients treated with definitive RT. Higher stage, older age and lack of chemotherapy were also associated with reduced OS.

Impact of Cardiac Dose on Overall Survival in Lung Stereotactic Body Radiotherapy (SBRT) Compared to Conventionally Fractionated Radiotherapy for Locally Advanced Non-Small Cell Lung Cancer (LA-NSCLC).

PURPOSE: To examine possible association between heart irradiation and Overall Survival (OS) in lung SBRT patients and to compare observed associations with cardiac toxicity models previously derived in LA-NSCLC patient studies. MATERIALS AND METHODS: 197 Patients treated with lung SBRT at Mayo Clinic Arizona were selected for this IRB-approved study. Multivariate Cox model with Akaike Information Criterion (AIC) was used to select patient specific covariates associated with OS. Heart dosimetry was represented by VD indices, which is a percentage of volume exposed to dose D or greater. Multivariate Cox models with patient specific covariates and single VD index per model was used to find a range of doses which were predictive for OS. A digital subdivision of the heart was further used to determine the spatial distribution of doses which were predictive for OS. A coarse subdivision divided heart into 4 segments, while the fine subdivision divided heart into 64 segments. Knowledge constrained Fused Lasso operator was used to derive a more complete model which correlated heart dosimetry with OS. Results of statistical analysis were compared to predictions of a model of cardiac toxicity in LA-NSCLC patients. RESULTS: Higher age (p < 0.001), higher stage (p < 0.001) and squamous cell histology (p = 0.001) were associated with reduced OS. Whole heart DVH analysis did not reveal associations between heart irradiation and reduced OS. Coarse subdivision of the heart into four segments revealed that the irradiation of two inferior segments of the heart with low doses was associated with reduced OS, V 2Gy in the right-inferior segment (HR = 1.012/1%, p = 0.02), and V 1Gy in the left-inferior segment (HR = 1.01/1%, p = 0.04). Maximum dose in the right-inferior segment of the heart was also associated with reduced OS (HR = 1.02/Gy, p = 0.02). Fine subdivision of the heart into 64 segments revealed that approximately 25% of heart volume in the inferior part of the heart (15/64 segments), when irradiated to doses in the 1 Gy - 5 Gy range, were predictive for reduced OS (HR = 1.01/1%, p = 0.01). A previously derived model of cardiac toxicity in LA-NSCLC patients did not predict a reduction of OS due to heart irradiation in lung SBRT patients, because of relatively low doses to the heart in most lung SBRT patients. CONCLUSIONS: Doses lower than 5 Gy in the inferior segments of the heart may be associated with reduced overall survival in patients treated for lung lesions with SBRT. Stage and histology of the disease, as well as patients' age, were also associated with overall survival. Comparisons of cardiac toxicity patterns in LA-NSCLC patients and lung SBRT patients suggest different etiology of cardiac toxicity in the two groups.

Technical Note: Multiple energy extraction techniques for synchrotron-based proton delivery systems may exacerbate motion interplay effects in lung cancer treatments.

PURPOSE: The multiple energy extraction (MEE) delivery technique for synchrotron-based proton delivery systems reduces beam delivery time by decelerating the beam multiple times during one accelerator spill, but this might cause additional plan quality degradation due to intrafractional motion. We seek to determine whether MEE causes significantly different plan quality degradation compared to single energy extraction (SEE) for lung cancer treatments due to the interplay effect. METHODS: Ten lung cancer patients treated with IMPT at our institution were nonrandomly sampled based on a representative range of tumor motion amplitudes, tumor volumes, and respiratory periods. Dose-volume histogram (DVH) indices from single-fraction SEE and MEE four-dimensional (4D) dynamic dose distributions were compared using the Wilcoxon signed-rank test. Distributions of monitor units (MU) to breathing phases were investigated for features associated with plan quality degradation. SEE and MEE DVH indices were compared in fractionated deliveries of the worst-case patient treatment scenario to evaluate the impact of fractionation. RESULTS: There were no clinically significant differences in target mean dose, target dose conformity, or dose to organs-at-risk between SEE and MEE in single-fraction delivery. Three patients had significantly worse dose homogeneity with MEE compared to SEE (single-fraction mean D5% -D95% increased by up to 9.6% of prescription dose), and plots of MU distribution to breathing phases showed synchronization patterns with MEE but not SEE. However, after 30 fractions the patient in the worst-case scenario had clinically acceptable target dose homogeneity and coverage with MEE (mean D5% -D95% increased by 1% compared to SEE). CONCLUSIONS: For some patients with breathing periods close to the mean spill duration, MEE resulted in significantly worse single-fraction target dose homogeneity compared to SEE due to the interplay effect. However, this was mitigated by fractionation, and target dose homogeneity and coverage were clinically acceptable after 30 fractions with MEE.

Technical Note: 4D robust optimization in small spot intensity-modulated proton therapy (IMPT) for distal esophageal carcinoma.

PURPOSE: To compare the dosimetric performances of small-spot three-dimensional (3D) and four-dimensional (4D) robustly optimized intensity-modulated proton (IMPT) plans in the presence of uncertainties and interplay effect simultaneously for distal esophageal carcinoma. METHOD AND MATERIALS: Thirteen (13) patients were selected and re-planned with small-spot ( σ  ~ 2-6 mm) 3D and 4D robust optimization in IMPT, respectively. The internal clinical target volumes (CTVhigh3d , CTVlow3d ) were used in 3D robust optimization. Different CTVs (CTVhigh4d , CTVlow4d ) were generated by subtracting an inner margin of the motion amplitudes in three cardinal directions from the internal CTVs and used in 4D robust optimization. All patients were prescribed the same dose to CTVs (50 Gy[RBE] for CTVhigh3d /CTVhigh4d and 45 Gy[RBE] for CTVlow3d /CTVlow4d ). Dose-volume histogram (DVH) indices were calculated to assess plan quality. Comprehensive plan robustness evaluations that consisted of 300 perturbed scenarios (10 different motion patterns to consider irregular motion (sampled from a Gaussian distribution) and 30 different uncertainties scenarios (sampled from a 4D uniform distribution) combined), were performed to quantify robustness to uncertainties and interplay effect simultaneously. Wilcoxon signed-rank test was used for statistical analysis. RESULTS: Compared to 3D robustly optimized plans, 4D robustly optimized plans had statistically improved target coverage and better sparing of lungs and heart (heart Dmean , P = 0.001; heart V30Gy[RBE] , P = 0.001) in the nominal scenario. 4D robustly optimized plans had better robustness in target dose coverage (CTVhigh4d V100% , P = 0.002) and the protection of lungs and heart (heart Dmean , P = 0.001; heart V30Gy[RBE] , P = 0.001) when uncertainties and interplay effect were considered simultaneously. CONCLUSIONS: Even with small spots in IMPT, 4D robust optimization outperformed 3D robust optimization in terms of normal tissue protection and robustness to uncertainties and interplay effect simultaneously. Our findings support the use of 4D robust optimization to treat distal esophageal carcinoma with small spots in IMPT.

Exploratory Investigation of Dose-Linear Energy Transfer (LET) Volume Histogram (DLVH) for Adverse Events Study in Intensity Modulated Proton Therapy (IMPT).

PURPOSE: We proposed a novel tool-a dose linear energy transfer (LET)-volume histogram (DLVH)-and performed an exploratory study to investigate rectal bleeding in prostate cancer treated with intensity modulated proton therapy. METHODS AND MATERIALS: The DLVH was constructed with dose and LET as 2 axes, and the normalized volume of the structure was contoured in the dose-LET plane as isovolume lines. We defined the DLVH index, DLv%(d,l) (ie, v% of the structure) to have a dose of ≥d Gy and an LET of ≥l keV/µm, similar to the dose-volume histogram index Dv%. Nine patients with prostate cancer with rectal bleeding (Common Terminology Criteria for Adverse Events grade ≥2) were included as the adverse event group, and 48 patients with no complications were considered the control group. A P value map was constructed by comparison of the DLVH indices of all patients between the 2 groups using the Mann-Whitney U test. Dose-LET volume constraints (DLVCs) were derived based on the P value map with a manual selection procedure facilitated by Spearman's correlation tests. The obtained DLVCs were further cross-validated using a multivariate support vector machine (SVM)-based normal tissue complication probability (NTCP) model with an independent testing data set composed of 8 adverse event and 13 control patients. RESULTS: We extracted 2 DLVC constraints. One DLVC was obtained, Vdose/LETboundary:2.5keVµmat 75 Gy to 3.2keVµmat8.65Gy <1.27% (DLVC1), revealing a high LET volume effect. The second DLVC, V(72.2Gy,0keVµm) < 2.23% (DVLC2), revealed a high dose volume effect. The SVM-based NTCP model with 2 DLVCs provided slightly superior performance than using dose only, with an area under the curve of 0.798 versus 0.779 for the testing data set. CONCLUSIONS: Our results demonstrated the importance of rectal "hot spots" in both high LET (DLVC1) and high dose (DLVC2) in inducing rectal bleeding. The SVM-based NTCP model confirmed the derived DLVCs as good predictors for rectal bleeding when intensity modulated proton therapy is used to treat prostate cancer.

Beam angle comparison for distal esophageal carcinoma patients treated with intensity-modulated proton therapy.

PURPOSE: To compare the dosimetric performances of intensity-modulated proton therapy (IMPT) plans generated with two different beam angle configurations (the Right-Left oblique posterior beams and the Superior-Inferior oblique posterior beams) for the treatment of distal esophageal carcinoma in the presence of uncertainties and interplay effect. METHODS AND MATERIALS: Twenty patients' IMPT plans were retrospectively selected, with 10 patients treated with the R-L oblique posterior beams (Group R-L) and the other 10 patients treated with the S-I oblique posterior beams (Group S-I). Patients in both groups were matched by their clinical target volumes (CTVs-high and low dose levels) and respiratory motion amplitudes. Dose-volume-histogram (DVH) indices were used to assess plan quality. DVH bandwidth was calculated to evaluate plan robustness. Interplay effect was quantified using four-dimensional (4D) dynamic dose calculation with random respiratory starting phase of each fraction. Normal tissue complication probability (NTCP) for heart, liver, and lung was calculated, respectively, to estimate the clinical outcomes. Wilcoxon signed-rank test was used for statistical comparison between the two groups. RESULTS: Compared with plans in Group R-L, plans in Group S-I resulted in significantly lower liver Dmean and lung V30Gy[RBE] with slightly higher but clinically acceptable spinal cord Dmax . Similar plan robustness was observed between the two groups. When interplay effect was considered, plans in Group S-I performed statistically better for heart Dmean and V30Gy[RBE] , lung Dmean and V5Gy[RBE] , and liver Dmean , with slightly increased but clinically acceptable spinal cord Dmax . NTCP for liver was significantly better in Group S-I. CONCLUSIONS: IMPT plans in Group S-I have better sparing of liver, heart, and lungs at the slight cost of spinal cord maximum dose protection, and are more interplay-effect resilient compared to IMPT plans in Group R-L. Our study supports the routine use of the S-I oblique posterior beams for the treatments of distal esophageal carcinoma.

Intensity-modulated proton therapy (IMPT) interplay effect evaluation of asymmetric breathing with simultaneous uncertainty considerations in patients with non-small cell lung cancer.

PURPOSE: Intensity-modulated proton therapy (IMPT) is sensitive to uncertainties from patient setup and proton beam range, as well as interplay effect. In addition, respiratory motion may vary from cycle to cycle, and also from day to day. These uncertainties can severely degrade the original plan quality and potentially affect patient's outcome. In this work, we developed a new tool to comprehensively consider the impact of all these uncertainties and provide plan robustness evaluation under them. METHODS: We developed a comprehensive plan robustness evaluation tool that considered both uncertainties from patient setup and proton beam range, as well as respiratory motion simultaneously. To mimic patients' respiratory motion, the time spent in each phase was randomly sampled based on patient-specific breathing pattern parameters as acquired during the four-dimensional (4D)-computed tomography (CT) simulation. Spots were then assigned to one specific phase according to the temporal relationship between spot delivery sequence and patients' respiratory motion. Dose in each phase was calculated by summing contributions from all the spots delivered in that phase. The final 4D dynamic dose was obtained by deforming all doses in each phase to the maximum exhalation phase. Three hundred (300) scenarios (10 different breathing patterns with 30 different setup and range uncertainty scenario combinations) were calculated for each plan. The dose-volume histograms (DVHs) band method was used to assess plan robustness. Benchmarking the tool as an application's example, we compared plan robustness under both three-dimensional (3D) and 4D robustly optimized IMPT plans for 10 nonrandomly selected patients with non-small cell lung cancer. RESULTS: The developed comprehensive plan robustness tool had been successfully applied to compare the plan robustness between 3D and 4D robustly optimized IMPT plans for 10 lung cancer patients. In the presence of interplay effect with uncertainties considered simultaneously, 4D robustly optimized plans provided significantly better CTV coverage (D95% , P = 0.002), CTV homogeneity (D5% -D95% , P = 0.002) with less target hot spots (D5% , P = 0.002), and target coverage robustness (CTV D95% bandwidth, P = 0.004) compared to 3D robustly optimized plans. Superior dose sparing of normal lung (lung Dmean , P = 0.020) favoring 4D plans and comparable normal tissue sparing including esophagus, heart, and spinal cord for both 3D and 4D plans were observed. The calculation time for all patients included in this study was 11.4 ± 2.6 min. CONCLUSION: A comprehensive plan robustness evaluation tool was successfully developed and benchmarked for plan robustness evaluation in the presence of interplay effect, setup and range uncertainties. The very high efficiency of this tool marks its clinical adaptation, highly practical and versatile nature, including possible real-time intra-fractional interplay effect evaluation as a potential application for future use.

Three-Dimensionally Printed On-Skin Radiation Shields Using High-Density Filament.

PURPOSE: Custom-fabricated lead shields are often used for superficial radiation treatments to reduce radiation doses to adjacent healthy tissue. However, the process for fabricating these lead shields is time consuming, labor intensive, and uncomfortable for patients. Alternatively, patient-specific shields can be 3-dimensionally (3D) printed from a high-density bronze-based filament to address these concerns. This study was performed to assess the shielding characteristics of 3D-printed bronze (3DPB) shields, demonstrate their clinical viability, and report the first ever published case of a patient treated with a 3DPB shield. METHODS AND MATERIALS: The transmission of 6 and 9 MeV electron beams through varying thicknesses of 3DPB was first measured. Percent depth doses and beam profiles were measured with flat 3DPB shields and equivalent lead shields to determine surface dose enhancement, output factors, and field widths. Two 3DPB shields were designed and fabricated for an anthropomorphic phantom, and phantom measurements were performed using optically stimulated luminescence dosimeters and film. Finally, 3DPB shields were used during the treatment of 7 patients' skin lesions. RESULTS: Ten and 15 mm of 3DPB were sufficient to shield 6 and 9 MeV electrons, respectively, by 95%. The 3DPB and lead shields had nearly identical beam widths (within 1%). Output factors were on average within 0.8% for bronze shields and 1.2% for lead shields relative to an unshielded field. The skin enhancement for bronze was higher than for lead by an average of 6.3%. Phantom measurements using 3DPB shields generally showed less than 3% transmission of the primary beam under the 3DPB shield. The patients' shields fit as designed and were all deemed clinically acceptable by their physicians. CONCLUSIONS: The 3DPB shields fit better than lead shields, are easier to design and manufacture, and have similar dosimetric properties. 3DPB shields are a viable clinical option for patient-specific superficial shielding.

Technical Note: Integrating an open source Monte Carlo code "MCsquare" for clinical use in intensity-modulated proton therapy.

PURPOSE: To commission an open source Monte Carlo (MC) dose engine, "MCsquare" for a synchrotron-based proton machine, integrate it into our in-house C++-based I/O user interface and our web-based software platform, expand its functionalities, and improve calculation efficiency for intensity-modulated proton therapy (IMPT). METHODS: We commissioned MCsquare using a double Gaussian beam model based on in-air lateral profiles, integrated depth dose of 97 beam energies, and measurements of various spread-out Bragg peaks (SOBPs). Then we integrated MCsquare into our C++-based dose calculation code and web-based second check platform "DOSeCHECK." We validated the commissioned MCsquare based on 12 different patient geometries and compared the dose calculation with a well-benchmarked GPU-accelerated MC (gMC) dose engine. We further improved the MCsquare efficiency by employing the computed tomography (CT) resampling approach. We also expanded its functionality by adding a linear energy transfer (LET)-related model-dependent biological dose calculation. RESULTS: Differences between MCsquare calculations and SOBP measurements were <2.5% (<1.5% for ~85% of measurements) in water. The dose distributions calculated using MCsquare agreed well with the results calculated using gMC in patient geometries. The average 3D gamma analysis (2%/2 mm) passing rates comparing MCsquare and gMC calculations in the 12 patient geometries were 98.0 ± 1.0%. The computation time to calculate one IMPT plan in patients' geometries using an inexpensive CPU workstation (Intel Xeon E5-2680 2.50 GHz) was 2.3 ± 1.8 min after the variable resolution technique was adopted. All calculations except for one craniospinal patient were finished within 3.5 min. CONCLUSIONS: MCsquare was successfully commissioned for a synchrotron-based proton beam therapy delivery system and integrated into our web-based second check platform. After adopting CT resampling and implementing LET model-dependent biological dose calculation capabilities, MCsquare will be sufficiently efficient and powerful to achieve Monte Carlo-based and LET-guided robust optimization in IMPT, which will be done in the future studies.

Detecting spatial susceptibility to cardiac toxicity of radiation therapy for lung cancer.

Radiation therapy (RT) is a commonly used approach for treating lung cancer. Because the lungs are close to the heart, radiation dose may inevitably spill to the heart, causing heart damage and diminishing treatment outcomes. There is an urgent need to better understand how treatment outcomes are affected by radiation dose spilled to the heart in order to optimize RT planning. However, despite the fact that dose distribution on the heart is 3-D, most existing research collapses the 3-D dose map into a 1-D histogram to be linked with outcomes. This ignores the spatial information. We propose a novel method that automatically searches for subregions of the heart that are susceptible to radiation toxicity, called Toxicity-Susceptible Subregions (TSSs), based on the 3-D dose distribution. We apply the proposed method to a real-world dataset and find TSSs that harbor the sinoatrial node of the electronic conduction system of the heart. Damage of the sinoatrial node by radiation toxicity disrupts the crucial function of the heart, leading to shortening of the overall survival. Our finding suggests that protective strategies may be developed to spare the TSSs, and thus helping RT planning achieve optimal results in treating lung cancer patients.

Impact of planned dose reporting methods on Gamma pass rates for IROC lung and liver motion phantoms treated with pencil beam scanning protons.

PURPOSE: The purpose of this study is to evaluate the impact of two methods of reporting planned dose distributions on the Gamma analysis pass rates for comparison with measured 2D film dose and simulated delivered 3D dose for proton pencil beam scanning treatment of the Imaging and Radiation Oncology Core (IROC) proton lung and liver mobile phantoms. METHODS AND MATERIALS: Four-dimensional (4D) computed-tomography (CT) image sets were acquired for IROC proton lung and liver mobile phantoms, which include dosimetry inserts that contains targets, thermoluminescent dosimeters and EBT2 films for plan dose verification. 4DCT measured fixed motion magnitudes were 1.3 and 1.0 cm for the lung and liver phantoms, respectively. To study the effects of motion magnitude on the Gamma analysis pass rate, three motion magnitudes for each phantom were simulated by creating virtual 4DCT image sets with motion magnitudes scaled from the scanned phantom motion by 50, 100, and 200%. The internal target volumes were contoured on the maximum intensity projection CTs of the 4DCTs for the lung phantom and on the minimum intensity projection CTs of the 4DCTs for the liver phantom. Treatment plans were optimized on the average intensity projection (AVE) CTs of the 4DCTs using the RayStation treatment planning system. Plan doses were calculated on the AVE CTs, which was defined as the planned AVE dose (method one). Plan doses were also calculated on all 10 phase CTs of the 4DCTs and were registered using target alignment to and equal-weight-summed on the 50% phase (T50) CT, which was defined as the planned 4D dose (method two). The planned AVE doses and 4D doses for phantom treatment were reported to IROC, and the 2D-2D Gamma analysis pass rates for measured film dose relative to the planned AVE and 4D doses were compared. To evaluate motion interplay effects, simulated delivered doses were calculated for each plan by sorting spots into corresponding respiratory phases using spot delivery time recorded in the log files by the beam delivery system to calculate each phase dose and accumulate dose to the T50 CTs. Ten random beam starting phases were used for each beam to obtain the range of the simulated delivered dose distributions. 3D-3D Gamma analyses were performed to compare the planned 4D/AVE doses with simulated delivered doses. RESULTS: The planned 4D dose matched better with the measured 2D film dose and simulated delivered 3D dose than the planned AVE dose. Using planned 4D dose as institution reported planned dose to IROC improved IROC film dose 2D-2D Gamma analysis pass rate from 92 to 96% on average for three films for the lung phantom (7% 5 mm), and from 92 to 94% in the sagittal plane for the liver phantom (7% 4 mm), respectively, compared with using the planned AVE dose. The 3D-3D Gamma analysis (3% 3 mm) pass rate showed that the simulated delivered doses for lung and liver phantoms using 10 random beam starting phases for each delivered beam matched the planned 4D dose significantly better than the planned AVE dose for phantom motions larger than 1 cm (p ≤ 0.04). CONCLUSIONS: It is recommended to use the planned 4D dose as the institution reported planned dose to IROC to compare with the measured film dose for proton mobile phantoms to improve film Gamma analysis pass rate in the IROC credentialing process.

Automation of routine elements for spot-scanning proton patient-specific quality assurance.

PURPOSE: At our institution, all proton patient plans undergo patient-specific quality assurance (PSQA) prior to treatment delivery. For intensity-modulated proton beam therapy, quality assurance is complex and time consuming, and it may involve multiple measurements per field. We reviewed our PSQA workflow and identified the steps that could be automated and developed solutions to improve efficiency. METHODS: We used the treatment planning system's (TPS) capability to support C# scripts to develop an Eclipse scripting application programming interface (ESAPI) script and automate the preparation of the verification phantom plan for measurements. A local area network (LAN) connection between our measurement equipment and shared database was established to facilitate equipment control, measurement data transfer, and storage. To improve the analysis of the measurement data, a Python script was developed to automatically perform a 2D-3D γ-index analysis comparing measurements in the plane of a two-dimensional detector array with TPS predictions in a water phantom for each acquired measurement. RESULTS: Device connection via LAN granted immediate access to the plan and measurement information for downstream analysis using an online software suite. Automated scripts applied to verification plans reduced time from preparation steps by at least 50%; time reduction from automating γ-index analysis was even more pronounced, dropping by a factor of 10. On average, we observed an overall time savings of 55% in completion of the PSQA per patient plan. CONCLUSIONS: The automation of the routine tasks in the PSQA workflow significantly reduced the time required per patient, reduced user fatigue, and frees up system users from routine and repetitive workflow steps allowing increased focus on evaluating key quality metrics.