Incorporating variable RBE in IMPT optimization for ependymoma.

PURPOSE: To study the dosimetric impact of incorporating variable relative biological effectiveness (RBE) of protons in optimizing intensity-modulated proton therapy (IMPT) treatment plans and to compare it with conventional constant RBE optimization and linear energy transfer (LET)-based optimization. METHODS: This study included 10 pediatric ependymoma patients with challenging anatomical features for treatment planning. Four plans were generated for each patient according to different optimization strategies: (1) constant RBE optimization (ConstRBEopt) considering standard-of-care dose requirements; (2) LET optimization (LETopt) using a composite cost function simultaneously optimizing dose-averaged LET (LETd ) and dose; (3) variable RBE optimization (VarRBEopt) using a recent phenomenological RBE model developed by McNamara et al.; and (4) hybrid RBE optimization (hRBEopt) assuming constant RBE for the target and variable RBE for organs at risk. By normalizing each plan to obtain the same target coverage in either constant or variable RBE, we compared dose, LETd , LET-weighted dose, and equivalent uniform dose between the different optimization approaches. RESULTS: We found that the LETopt plans consistently achieved increased LET in tumor targets and similar or decreased LET in critical organs compared to other plans. On average, the VarRBEopt plans achieved lower mean and maximum doses with both constant and variable RBE in the brainstem and spinal cord for all 10 patients. To compensate for the underdosing of targets with 1.1 RBE for the VarRBEopt plans, the hRBEopt plans achieved higher physical dose in targets and reduced mean and especially maximum variable RBE doses compared to the ConstRBEopt and LETopt plans. CONCLUSION: We demonstrated the feasibility of directly incorporating variable RBE models in IMPT optimization. A hybrid RBE optimization strategy showed potential for clinical implementation by maintaining all current dose limits and reducing the incidence of high RBE in critical normal tissues in ependymoma patients.

Impact of RBE variations on risk estimates of temporal lobe necrosis in patients treated with intensity-modulated proton therapy for head and neck cancer.

BACKGROUND: Temporal lobe necrosis (TLN) is a potential late effect after radiotherapy for skull base head and neck cancer (HNC). Several photon-derived dose constraints and normal tissue complication probability (NTCP) models have been proposed, however variation in relative biological effectiveness (RBE) may challenge the applicability of these dose constraints and models in proton therapy. The purpose of this study was therefore to investigate the influence of RBE variations on risk estimates of TLN after Intensity-Modulated Proton Therapy for HNC. MATERIAL AND METHODS: Seventy-five temporal lobes from 45 previously treated patients were included in the analysis. Sixteen temporal lobes had radiation associated Magnetic Resonance image changes (TLIC) suspected to be early signs of TLN. Fixed (RWDFix) and variable RBE-weighed doses (RWDVar) were calculated using RBE = 1.1 and two RBE models, respectively. RWDFix and RWDVar for temporal lobes were compared using Friedman's test. Based on RWDFix, six NTCP models were fitted and internally validated through bootstrapping. Estimated probabilities from RWDFix and RWDVar were compared using paired Wilcoxon test. Seven dose constraints were evaluated separately for RWDFix and RWDVar by calculating the observed proportion of TLIC in temporal lobes meeting the specific dose constraints. RESULTS: RWDVar were significantly higher than RWDFix (p < 0.01). NTCP model performance was good (AUC:0.79-0.84). The median difference in estimated probability between RWDFix and RWDVar ranged between 5.3% and 20.0% points (p < 0.01), with V60GyRBE and DMax at the smallest and largest differences, respectively. The proportion of TLIC was higher for RWDFix (4.0%-13.1%) versus RWDVar (1.3%-5.3%). For V65GyRBE ≤ 0.03 cc the proportion of TLIC was less than 5% for both RWDFix and RWDVar. CONCLUSION: NTCP estimates were significantly influenced by RBE variations. Dmax as model predictor resulted in the largest deviations in risk estimates between RWDFix and RWDVar. V65GyRBE ≤ 0.03 cc was the most consistent dose constraint for RWDFix and RWDVar.

Monte Carlo evaluation of target dose coverage in lung stereotactic body radiation therapy with flattening filter-free beams.

Aim: Previous studies showed that replacing conventional flattened beams (FF) with flattening filter-free (FFF) beams improves the therapeutic ratio in lung stereotactic body radiation therapy (SBRT), but these findings could have been impacted by dose calculation uncertainties caused by the heterogeneity of the thoracic anatomy and by respiratory motion, which were particularly high for target coverage. In this study, we minimized such uncertainties by calculating doses using high-spatial-resolution Monte Carlo and four-dimensional computed tomography (4DCT) images. We aimed to evaluate more reliably the benefits of using FFF beams for lung SBRT. Materials and methods: For a cohort of 15 patients with early stage lung cancer that we investigated in a previous treatment planning study, we recalculated dose distributions with Monte Carlo using 4DCT images. This included fifteen FF and fifteen FFF treatment plans. Results: Compared to Monte Carlo, the treatment planning system (TPS) over-predicted doses in low-dose regions of the planning target volume. For most patients, replacing FF beams with FFF beams improved target coverage, tumor control, and uncomplicated tumor control probabilities. Conclusions: Monte Carlo tends to reveal deficiencies in target coverage compared to coverage predicted by the TPS. Our data support previously reported benefits of using FFF beams for lung SBRT.

A comparison of two methodologies for radiotherapy treatment plan optimization and QA for clinical trials.

BACKGROUND AND PURPOSE: The efficacy of clinical trials and the outcome of patient treatment are dependent on the quality assurance (QA) of radiation therapy (RT) plans. There are two widely utilized approaches that include plan optimization guidance created based on patient-specific anatomy. This study examined these two techniques for dose-volume histogram predictions, RT plan optimizations, and prospective QA processes, namely the knowledge-based planning (KBP) technique and another first principle (FP) technique. METHODS: This analysis included 60, 44, and 10 RT plans from three Radiation Therapy Oncology Group (RTOG) multi-institutional trials: RTOG 0631 (Spine SRS), RTOG 1308 (NSCLC), and RTOG 0522 (H&N), respectively. Both approaches were compared in terms of dose prediction and plan optimization. The dose predictions were also compared to the original plan submitted to the trials for the QA procedure. RESULTS: For the RTOG 0631 (Spine SRS) and RTOG 0522 (H&N) plans, the dose predictions from both techniques have correlation coefficients of >0.9. The RT plans that were re-optimized based on the predictions from both techniques showed similar quality, with no statistically significant differences in target coverage or organ-at-risk sparing. The predictions of mean lung and heart doses from both methods for RTOG1308 patients, on the other hand, have a discrepancy of up to 14 Gy. CONCLUSIONS: Both methods are valuable tools for optimization guidance of RT plans for Spine SRS and Head and Neck cases, as well as for QA purposes. On the other hand, the findings suggest that KBP may be more feasible in the case of inoperable lung cancer patients who are treated with IMRT plans that have spatially unevenly distributed beam angles.

Geometric and dosimetric accuracy of deformable image registration between average-intensity images for 4DCT-based adaptive radiotherapy for non-small cell lung cancer.

PURPOSE: Re-planning for four-dimensional computed tomography (4DCT)-based lung adaptive radiotherapy commonly requires deformable dose mapping between the planning average-intensity image (AVG) and the newly acquired AVG. However, such AVG-AVG deformable image registration (DIR) lacks accuracy assessment. The current work quantified and compared geometric accuracies of AVG-AVG DIR and corresponding phase-phase DIRs, and subsequently investigated the clinical impact of such AVG-AVG DIR on deformable dose mapping. METHODS AND MATERIALS: Hybrid intensity-based AVG-AVG and phase-phase DIRs were performed between the planning and mid-treatment 4DCTs of 28 non-small cell lung cancer patients. An automated landmark identification algorithm detected vessel bifurcation pairs in both lungs. Target registration error (TRE) of these landmark pairs was calculated for both DIR types. The correlation between TRE and respiratory-induced landmark motion in the planning 4DCT was analyzed. Global and local dose metrics were used to assess the clinical implications of AVG-AVG deformable dose mapping with both DIR types. RESULTS: TRE of AVG-AVG and phase-phase DIRs averaged 3.2 ± 1.0 and 2.6 ± 0.8 mm respectively (p < 0.001). Using AVG-AVG DIR, TREs for landmarks with <10 mm motion averaged 2.9 ± 2.0 mm, compared to 3.1 ± 1.9 mm for the remaining landmarks (p < 0.01). Comparatively, no significant difference was demonstrated for phase-phase DIRs. Dosimetrically, no significant difference in global dose metrics was observed between doses mapped with AVG-AVG DIR and the phase-phase DIR, but a positive linear relationship existed (p = 0.04) between the TRE of AVG-AVG DIR and local dose difference. CONCLUSIONS: When the region of interest experiences <10 mm respiratory-induced motion, AVG-AVG DIR may provide sufficient geometric accuracy; conversely, extra attention is warranted, and phase-phase DIR is recommended. Dosimetrically, the differences in geometric accuracy between AVG-AVG and phase-phase DIRs did not impact global lung-based metrics. However, as more localized dose metrics are needed for toxicity assessment, phase-phase DIR may be required as its lower mean TRE improved voxel-based dosimetry.

Using FFF Beams to Improve the Therapeutic Ratio of Lung SBRT.

Aim: To investigate the extent to which lung stereotactic body radiotherapy (SBRT) treatment plans can be improved by replacing conventional flattening filter (FF) beams with flattening filter-free (FFF) beams. Material and Methods: We selected 15 patients who had received SBRT with conventional 6-MV photon beams for early-stage lung cancer. We imported the patients' treatment plans into the Eclipse 13.6 treatment planning system, in which we configured the AAA dose calculation model using representative beam data for a TrueBeam accelerator operated in 6-MV FFF mode. We then created new treatment plans by replacing the conventional FF beams in the original plans with FFF beams. Results: The FFF plans had better target coverage than the original FF plans did. For the planning target volume, FFF plans significantly improved the D98, D95, D90, homogeneity index, and uncomplicated tumor control probability. In most cases, the doses to organs at risk were lower in FFF plans. FFF plans significantly reduced the mean lung dose, V10, V20, V30, and normal tissue complication probability for the total lung and improved the dosimetric indices for the ipsilateral lung. For most patients, FFF beams achieved lower maximum doses to the esophagus, heart, and the spinal cord; and a lower chest wall V30. Findings: Compared with FF beams, FFF beams achieved lower doses to organs at risk, especially the lung, without compromising tumor coverage; in fact, FFF beams improved coverage in most cases. Thus, replacing FF beams with FFF beams can achieve a better therapeutic ratio.

Impact of Intra-Fractional Motion on Dose Distributions in Lung IMRT.

AIM: To investigate the impact of intra-fractional motion on dose distribution in patients treated with intensity-modulated radiation therapy (IMRT) for lung cancer. MATERIALS AND METHODS: Twenty patients who had undergone IMRT for non-small cell lung cancer were selected for this retrospective study. For each patient, a four-dimensional computed tomography (CT) image set was acquired and clinical treatment plans were developed using the average CT. Dose distributions were then re-calculated for each of the 10 phases of respiratory cycle and combined using deformable image registration to produce cumulative dose distributions that were compared with the clinical treatment plans. RESULTS: Intra-fractional motion reduced planning target volume (PTV) coverage in all patients. The median reduction of PTV volume covered by the prescription isodose was 3.4%; D98 was reduced by 3.1 Gy. Changes in the mean lung dose were within ±0.7 Gy. V20 for the lung increased in most patients; the median increase was 1.6%. The dose to the spinal cord was unaffected by intra-fractional motion. The dose to the heart was slightly reduced in most patients. The median reduction in the mean heart dose was 0.22 Gy, and V30 was reduced by 2.5%.The maximum dose to the esophagus was also reduced in most patients, by 0.74 Gy, whereas V50 did not change significantly. The median number of points in which dose differences exceeded 3%/3 mm was 6.2%. FINDINGS: Intra-fractional anatomical changes reduce PTV coverage compared to the coverage predicted by clinical treatment planning systems that use the average CT for dose calculation. Doses to organs at risk were mostly over-predicted.

Proton therapy for adults with mediastinal lymphomas: the International Lymphoma Radiation Oncology Group guidelines.

Among adult lymphoma survivors, radiation treatment techniques that increase the excess radiation dose to organs at risk (OARs) put patients at risk for increased side effects, especially late toxicities. Minimizing radiation to OARs in adults patients with Hodgkin and non-Hodgkin lymphomas involving the mediastinum is the deciding factor for the choice of treatment modality. Proton therapy may help to reduce the radiation dose to the OARs and reduce toxicities, especially the risks for cardiac morbidity and second cancers. Because proton therapy may have some disadvantages, identifying the patients and the circumstances that may benefit the most from proton therapy is important. We present modern guidelines to identify adult lymphoma patients who may derive the greatest benefit from proton therapy, along with an analysis of the advantages and disadvantages of proton treatment.

Impact of respiratory motion on variable relative biological effectiveness in 4D-dose distributions of proton therapy.

BACKGROUND: Organ motion during radiation therapy with scanned protons leads to deviations between the planned and the delivered physical dose. Using a constant relative biological effectiveness (RBE) of 1.1 linearly maps these deviations into RBE-weighted dose. However, a constant value cannot account for potential nonlinear variations in RBE suggested by variable RBE models. Here, we study the impact of motion on recalculations of RBE-weighted dose distributions using a phenomenological variable RBE model. MATERIAL AND METHODS: 4D-dose calculation including variable RBE was implemented in the open source treatment planning toolkit matRad. Four scenarios were compared for one field and two field proton treatments for a liver cancer patient assuming (α∕ß)x = 2 Gy and (α∕ß)x = 10 Gy: (A) the optimized static dose distribution with constant RBE, (B) a static recalculation with variable RBE, (C) a 4D-dose recalculation with constant RBE and (D) a 4D-dose recalculation with variable RBE. For (B) and (D), the variable RBE was calculated by the model proposed by McNamara. For (C), the physical dose was accumulated with direct dose mapping; for (D), dose-weighted radio-sensitivity parameters of the linear quadratic model were accumulated to model synergistic irradiation effects on RBE. RESULTS: Dose recalculation with variable RBE led to an elevated biological dose at the end of the proton field, while 4D-dose recalculation exhibited random deviations everywhere in the radiation field depending on the interplay of beam delivery and organ motion. For a single beam treatment assuming (α∕ß)x = 2 Gy, D95% was 1.98 Gy (RBE) (A), 2.15 Gy (RBE) (B), 1.81 Gy (RBE) (C) and 1.98 Gy (RBE) (D). The homogeneity index was 1.04 (A), 1.08 (B), 1.23 (C) and 1.25 (D). CONCLUSION: For the studied liver case, intrafractional motion did not reduce the modulation of the RBE-weighted dose postulated by variable RBE models for proton treatments.

Radiobiological issues in proton therapy.

BACKGROUND: The relative biological effectiveness (RBE) for particle therapy is a complex function of particle type, radiation dose, linear energy transfer (LET), cell type, endpoint, etc. In the clinical practice of proton therapy, the RBE is assumed to have a fixed value of 1.1. This assumption, along with the effects of physical uncertainties, may mean that the biologically effective dose distributions received by the patient may be significantly different from what is seen on treatment plans. This may contribute to unforeseen toxicities and/or failure to control the disease. Variability of Proton RBE: It has been shown experimentally that proton RBE varies significantly along the beam path, especially near the end of the particle range. While there is now an increasing acceptance that proton RBE is variable, there is an ongoing debate about whether to change the current clinical practice. Clinical Evidence: A rationale against the change is the uncertainty in the models of variable RBE. Secondly, so far there is no clear clinical evidence of the harm of assuming proton RBE to be 1.1. It is conceivable, however, that the evidence is masked partially by physical uncertainties. It is, therefore, plausible that reduction in uncertainties and their incorporation in the estimation of dose actually delivered may isolate and reveal the variability of RBE in clinical practice. Nevertheless, clinical evidence of RBE variability is slowly emerging as more patients are treated with protons and their response data are analyzed. Modelling and Incorporation of RBE in the Optimization of Proton Therapy: The improvement in the knowledge of RBE could lead to better understanding of outcomes of proton therapy and in the improvement of models to predict RBE. Prospectively, the incorporation of such models in the optimization of intensity-modulated proton therapy could lead to improvements in the therapeutic ratio of proton therapy.

Comparison of linear and nonlinear programming approaches for "worst case dose" and "minmax" robust optimization of intensity-modulated proton therapy dose distributions.

Robust optimization of intensity-modulated proton therapy (IMPT) takes uncertainties into account during spot weight optimization and leads to dose distributions that are resilient to uncertainties. Previous studies demonstrated benefits of linear programming (LP) for IMPT in terms of delivery efficiency by considerably reducing the number of spots required for the same quality of plans. However, a reduction in the number of spots may lead to loss of robustness. The purpose of this study was to evaluate and compare the performance in terms of plan quality and robustness of two robust optimization approaches using LP and nonlinear programming (NLP) models. The so-called "worst case dose" and "minmax" robust optimization approaches and conventional planning target volume (PTV)-based optimization approach were applied to designing IMPT plans for five patients: two with prostate cancer, one with skull-based cancer, and two with head and neck cancer. For each approach, both LP and NLP models were used. Thus, for each case, six sets of IMPT plans were generated and assessed: LP-PTV-based, NLP-PTV-based, LP-worst case dose, NLP-worst case dose, LP-minmax, and NLP-minmax. The four robust optimization methods behaved differently from patient to patient, and no method emerged as superior to the others in terms of nominal plan quality and robustness against uncertainties. The plans generated using LP-based robust optimization were more robust regarding patient setup and range uncertainties than were those generated using NLP-based robust optimization for the prostate cancer patients. However, the robustness of plans generated using NLP-based methods was superior for the skull-based and head and neck cancer patients. Overall, LP-based methods were suitable for the less challenging cancer cases in which all uncertainty scenarios were able to satisfy tight dose constraints, while NLP performed better in more difficult cases in which most uncertainty scenarios were hard to meet tight dose limits. For robust optimization, the worst case dose approach was less sensitive to uncertainties than was the minmax approach for the prostate and skull-based cancer patients, whereas the minmax approach was superior for the head and neck cancer patients. The robustness of the IMPT plans was remarkably better after robust optimization than after PTV-based optimization, and the NLP-PTV-based optimization outperformed the LP-PTV-based optimization regarding robustness of clinical target volume coverage. In addition, plans generated using LP-based methods had notably fewer scanning spots than did those generated using NLP-based methods.

Perturbation of water-equivalent thickness as a surrogate for respiratory motion in proton therapy.

Respiratory motion is traditionally assessed using tumor motion magnitude. In proton therapy, respiratory motion causes density variations along the beam path that result in uncertainties of proton range. This work has investigated the use of water-equivalent thickness (WET) to quantitatively assess the effects of respiratory motion on calculated dose in passively scattered proton therapy (PSPT). A cohort of 29 locally advanced non-small cell lung cancer patients treated with 87 PSPT treatment fields were selected for analysis. The variation in WET (ΔWET) along each field was calculated between exhale and inhale phases of the simulation four-dimensional computed tomography. The change in calculated dose (ΔDose) between full-inhale and full-exhale phase was quantified for each field using dose differences, 3D gamma analysis, and differential area under the curve (ΔAUC) analysis. Pearson correlation coefficients were calculated between ΔDose and ΔWET. Three PSPT plans were redesigned using field angles to minimize variations in ΔWET and compared to the original plans. The median ΔWET over 87 treatment fields ranged from 1-9 mm, while the ΔWET 95th percentile value ranged up to 42 mm. The ΔWET was significantly correlated (p < 0.001) to the ΔDose for all metrics analyzed. The patient plans that were redesigned using ΔWET analysis to select field angles were more robust to the effects of respiratory motion, as ΔAUC values were reduced by more than 60% in all three cases. The tumor motion magnitude alone does not capture the potential dosimetric error due to respiratory motion because the proton range is sensitive to the motion of all patient anatomy. The use of ΔWET has been demonstrated to identify situations where respiratory motion can impact the calculated dose. Angular analysis of ΔWET may be capable of designing radiotherapy plans that are more robust to the effects of respiratory motion.

Is there an impact of heart exposure on the incidence of radiation pneumonitis? Analysis of data from a large clinical cohort.

BACKGROUND: The goal of the present study was to determine, in a large clinical cohort, whether incidental radiation exposure to the heart during definitive radiotherapy of inoperable non-small cell lung cancer (NSCLC) detectably increased the risk of radiation pneumonitis (RP) beyond that resulting from radiation exposure to lung. MATERIAL AND METHODS: Data were analyzed from all patients who received definitive three-dimensional (3D) concurrent radiotherapy or intensity-modulated radiotherapy for the treatment of NSCLC over a 10-year period at our institution, except those who had previous lung cancer or for whom radiation treatment plans were unavailable for calculation of heart and lung dose-volume histograms (DVHs). Parameters computed from heart and lung DVHs included mean lung dose (MLD), effective lung dose computed using volume parameter n = 0.5 (Deff), mean heart dose (MHD), percentage of heart receiving > 65 Gy (V65), and minimum dose to the hottest 10% of heart (D10). Univariate and multivariate normal-tissue complication probability (NTCP) models were used to analyze incidence of Grade ≥ 2 or Grade ≥ 3 RP as a function of these and other parameters. RESULTS: The study cohort included 629 patients, with crude rates of Grade ≥ 2 RP and Grade ≥ 3 RP of N = 263 (42%) and N = 124 (20%), respectively. Univariate NTCP models based on dosimetric lung parameters (MLD and Deff) fit the data better than models based on univariate heart parameters (heart D10, heart V65 or MHD). In multivariate modeling, incorporation of heart parameters did not significantly improve the fit of RP risk models based on lung parameters alone (p > 0.38 in each case). CONCLUSIONS: In this large clinical cohort, there was no evidence that incidental heart exposure during radiotherapy of NSCLC had a detectable impact on the occurrence of moderate or severe RP.

Personalized Composite Dosimetric Score-Based Machine Learning Model of Severe Radiation-Induced Lymphopenia Among Patients With Esophageal Cancer.

PURPOSE: Radiation-induced lymphopenia (RIL) is common among patients undergoing radiation therapy (RT)' Severe RIL has been linked to adverse outcomes. The severity and risk of RIL can be predicted from baseline clinical characteristics and dosimetric parameters. However, dosimetric parameters, e.g. dose-volume (DV) indices, are highly correlated with one another and are only weakly associated with RIL. Here we introduce the novel concept of "composite dosimetric score" (CDS) as the index that condenses the dose distribution in immune tissues of interest to study the dosimetric dependence of RIL. We derived an improved multivariate classification scheme for risk of grade 4 RIL (G4RIL), based on this novel RT dosimetric feature, for patients receiving chemo RT for esophageal cancer. METHODS AND MATERIALS: DV indices were extracted for 734 patients who received chemo RT for biopsy-proven esophageal cancer. Nonnegative matrix factorization was used to project the DV indices of lung, heart, and spleen into a single CDS; XGBoost was employed to explore significant interactions among predictors; and logistic regression was applied to combine the resultant CDS with baseline clinical factors and interaction terms to facilitate individualized prediction of immunotoxicity. Five-fold cross-validation was applied to evaluate the model performance. RESULTS: The CDS for selected immune organs at risk (ie, heart, lung, and spleen) (OR 1.791; 95 CI [1.350, 2.377]) was a statistically significant risk determinant for G4RIL. Pearson correlation coefficients for CDS versus G4RIL risk for individual immune organs at risk were greater than any single DV indicx. Personalized prediction of G4RIL based on CDS and 4 clinical risk factors yielded an area under the curve value of 0.78. Interaction between age and CDS revealed that G4RIL risk increased more sharply with increasing CDS for patients aged ≥65 years. CONCLUSIONS: Risk of immunotoxicity for patients undergoing chemo RT for esophageal cancer can be predicted by CDS. The CDS concept can be extended to immunotoxicity in other cancer types and in dose-response models currently based on DV indices. Personalized treatment planning should leverage composite dosimetric scoring methods rather than using individual or subsets of DV indices.

Severe Lymphopenia During Chemoradiation Therapy for Esophageal Cancer: Comprehensive Analysis of Randomized Phase 2B Trial of Proton Beam Therapy Versus Intensity Modulated Radiation Therapy.

PURPOSE: Lymphocytes play an important role in antitumor immunity; however, they are also especially vulnerable to depletion during chemoradiation therapy (CRT). The purpose of this study was to compare the incidence of grade 4 lymphopenia (G4L) between proton beam therapy (PBT) and intensity modulated photon radiation therapy (IMRT) in patients with esophageal cancer treated with CRT in a completed randomized trial and to ascertain patient heterogeneity to G4L risk based on treatment and established prognostic factors. METHODS AND MATERIALS: Between April 2012 and March 2019, a single-institution, open-label, nonblinded, phase 2 randomized trial (NCT01512589) was conducted at the University of Texas MD Anderson Cancer Center. Patients were randomly assigned to IMRT or PBT, either definitively or preoperatively. This secondary analysis of the randomized trial was G4L during concurrent CRT according to Common Terminology Criteria for Adverse Events version 5.0. RESULTS: Among 105 patients evaluable for analysis, 44 patients (42%) experienced G4L at a median of 28 days after the start date of concurrent CRT. Induction chemotherapy (P = .003), baseline absolute lymphocyte count (P < .001), radiation therapy modality (P = .002), and planning treatment volume (P = .033) were found to be significantly associated with G4L. Multivariate classification analysis partitioned patients into 5 subgroups for whom the incidence of G4L was observed in 0%, 14%, 35%, 70%, and 100% of patients. The benefit of PBT over IMRT was most pronounced in patients with an intermediate baseline absolute lymphocyte count and large planning treatment volume (P = .011). CONCLUSIONS: This is the first prospective evidence that limiting dose scatter by PBT significantly reduced the incidence of G4L, especially in the intermediate-risk patients. The implication of this immune-sparing effect of PBT, especially in the context of standard adjuvant immunotherapy, needs further examination in the current phase 3 randomized trials.

Comparing 90-Day Postoperative Mortality After Neoadjuvant Proton-Based Versus Photon-Based Chemoradiotherapy for Esophageal Cancer.

Purpose: Evidence suggests that proton-beam therapy (PBT) results in less toxicity and postoperative complications compared to photon-based radiotherapy in patients who receive chemoradiotherapy followed by esophagectomy for cancer. Ninety-day mortality (90DM) is an important measure of the postoperative (nononcologic) outcome as proxy of quality-of-care. We hypothesize that PBT could reduce 90DM compared to photon-based radiotherapy. Materials and Methods: From a single-center retrospective database patients treated with chemoradiotherapy before esophagectomy for cancer were selected (1998-2022). Univariable logistic regression was used to study the association of radiotherapy modality with 90DM. Three separate methods were applied to adjust for confounding bias, including multivariable logistic regression, propensity score matching, and inverse probability of treatment weighting. Stratified analysis for the age threshold that maximized the difference in 90DM (ie, ≥67 vs <67 years) was performed. Results: A total of 894 eligible patients were included and 90DM was 5/202 (2.5%) in the PBT versus 29/692 (4.2%) in the photon-based radiotherapy group (P = .262). After adjustment for age and tumor location, PBT versus photon-based radiotherapy was not significantly associated with 90DM (P = .491). The 90DM was not significantly different for PBT versus photon-based radiotherapy in the propensity score matching (P = .379) and inverse probability of treatment weighting cohort (P = .426). The stratified analysis revealed that in patients aged ≥67 years, PBT was associated with decreased 90DM (1.3% vs 8.8%; P = .026). Higher age significantly increased 90DM risk within the photon-based radiotherapy (8.8% vs 2.7%; P = .001), but not within the PBT group (1.3% vs 3.2%; P = .651). Conclusion: No statistically significant difference was observed in postoperative 90DM after esophagectomy for cancer between PBT and photon-based neoadjuvant chemoradiotherapy. However, among older patients a signal was observed that PBT may reduce 90DM risk.

Quantifying the Effect of 4-Dimensional Computed Tomography-Based Deformable Dose Accumulation on Representing Radiation Damage for Patients with Locally Advanced Non-Small Cell Lung Cancer Treated with Standard-Fractionated Intensity-Modulated Radiation Therapy.

PURPOSE: The aim of this study was to investigate the dosimetric and clinical effects of 4-dimensional computed tomography (4DCT)-based longitudinal dose accumulation in patients with locally advanced non-small cell lung cancer treated with standard-fractionated intensity-modulated radiation therapy (IMRT). METHODS AND MATERIALS: Sixty-seven patients were retrospectively selected from a randomized clinical trial. Their original IMRT plan, planning and verification 4DCTs, and ∼4-month posttreatment follow-up CTs were imported into a commercial treatment planning system. Two deformable image registration algorithms were implemented for dose accumulation, and their accuracies were assessed. The planned and accumulated doses computed using average-intensity images or phase images were compared. At the organ level, mean lung dose and normal-tissue complication probability (NTCP) for grade ≥2 radiation pneumonitis were compared. At the region level, mean dose in lung subsections and the volumetric overlap between isodose intervals were compared. At the voxel level, the accuracy in estimating the delivered dose was compared by evaluating the fit of a dose versus radiographic image density change (IDC) model. The dose-IDC model fit was also compared for subcohorts based on the magnitude of NTCP difference (|ΔNTCP|) between planned and accumulated doses. RESULTS: Deformable image registration accuracy was quantified, and the uncertainty was considered for the voxel-level analysis. Compared with planned doses, accumulated doses on average resulted in <1-Gy lung dose increase and <2% NTCP increase (up to 8.2 Gy and 18.8% for a patient, respectively). Volumetric overlap of isodose intervals between the planned and accumulated dose distributions ranged from 0.01 to 0.93. Voxel-level dose-IDC models demonstrated a fit improvement from planned dose to accumulated dose (pseudo-R2 increased 0.0023) and a further improvement for patients with ≥2% |ΔNTCP| versus for patients with <2% |ΔNTCP|. CONCLUSIONS: With a relatively large cohort, robust image registrations, multilevel metric comparisons, and radiographic image-based evidence, we demonstrated that dose accumulation more accurately represents the delivered dose and can be especially beneficial for patients with greater longitudinal response.

Intensity-modulated radiation therapy - you can have your cake and eat it too!

PURPOSE: Arguably, intensity-modulated radiotherapy (IMRT) is one of the most important contributions by physicists in medicine to the treatment of cancer. It enabled the achievement of dose distributions that allowed treatment of targets to high doses and yet spared normal tissues to degrees previously considered virtually impossible. The concept underlying IMRT for photons, originally called "inverse planning," was introduced by Anders Brahme in the early 1980s. Since then, many physicists have been conducting research to advance the state of the art, to overcome obstacles discovered and to develop tools and techniques to translate it clinically. IMRT was first implemented clinically the early 1990's. Soon thereafter, it was also incorporated into intensity modulated particle therapy (IMPT), initially for protons and eventually for ions heavier than protons. METHODS: In IMRT, intensities of small segments or beamlets of multiple incident beams are optimized to produce the best approximation of desired dose distributions that deliver the requisite tumor dose and maximally spare normal tissues. Such dose distributions are delivered using dynamic multi-leaf collimators whose leaves move continuously while the radiation is on or in a step-and-shoot fashion. For IMPT, there is an extra degree of freedom, that of energies of particle beamlets, which are part of the optimization process. IMPT dose distributions allow substantial additional sparing of normal tissues. They are delivered using scanning beams. RESULTS: IMRT and IMPT represents a true paradigm shift in the way cancer patients are treated. Their evolution, especially that of IMPT, is still ongoing, with new advancements being added continually. With the improving understanding of the immunomodulatory effects of radiotherapy, it is quite plausible that a combination of IMPT and immunotherapy will mature into one of the most effective ways of treating not just localized cancers but also systemic disease.

GPU-accelerated calculation of proton microdosimetric spectra as a function of target size, proton energy, and bounding volume size.

Objective.Shortcomings of dose-averaged linear energy transfer (LETD), the quantity which is most commonly used to quantify proton relative biological effectiveness, have long been recognized. Microdosimetric spectra may overcome the limitations of LETDbut are extremely computationally demanding to calculate. A systematic library of lineal energy spectra for monoenergetic protons could enable rapid determination of microdosimetric spectra in a clinical environment. The objective of this work was to calculate and validate such a library of lineal energy spectra.Approach. SuperTrack, a GPU-accelerated CUDA/C++ based application, was developed to superimpose tracks calculated using Geant4 onto targets of interest and to compute microdosimetric spectra. Lineal energy spectra of protons with energies from 0.1 to 100 MeV were determined in spherical targets of diameters from 1 nm to 10µm and in bounding voxels with side lengths of 5µm and 3 mm.Main results.Compared to an analogous Geant4-based application, SuperTrack is up to 3500 times more computationally efficient if each track is resampled 1000 times. Dose spectra of lineal energy and dose-mean lineal energy calculated with SuperTrack were consistent with values published in the literature and with comparison to a Geant4 simulation. Using SuperTrack, we developed the largest known library of proton microdosimetric spectra as a function of primary proton energy, target size, and bounding volume size.Significance. SuperTrack greatly increases the computational efficiency of the calculation of microdosimetric spectra. The elevated lineal energy observed in a 3 mm side length bounding volume suggests that lineal energy spectra determined experimentally or computed in small bounding volumes may not be representative of the lineal energy spectra in voxels of a dose calculation grid. The library of lineal energy spectra calculated in this work could be integrated with a treatment planning system for rapid determination of lineal energy spectra in patient geometries.

Declarations of Independence: How Embedded Multicollinearity Errors Affect Dosimetric and Other Complex Analyses in Radiation Oncology.

The statistical technique of multiple regression, commonly referred to as "multivariable regression," is often used in clinical research to quantify the relationships between multiple predictor variables and a single outcome variable of interest. The foundational theory underpinning multivariable regression assumes that all predictor variables are independent of one another. In other words, the effect of each independent variable is measured by its contribution to the regression equation while all other variables remain unchanged. In the presence of correlations between two or more variables, however, it is impossible to change one variable without a consequent change in the variable(s) it is linked to. This condition, known as "multicollinearity," can introduce errors into multivariable regression models by affecting estimates of the regression coefficients that quantify the relationship between individual predictor variables and the outcome variable. Errors that arise due to violations of the multicollinearity assumption are of special interest to radiation oncology researchers. Because of high levels of correlation among variables derived from points along individual organ dose-volume histogram (DVH) curves, as well as strong intercorrelations among dose-volume parameters in neighboring organs, dosimetric analyses are particularly subject to multicollinearity errors. For example, dose-volume parameters for the heart are strongly correlated not only with other points along the heart DVH curve but are likely also correlated with dose-volume parameters in neighboring organs such as the lung. In this paper, we describe the problem of multicollinearity in accessible terms and discuss examples of violations of the multicollinearity assumption within the radiation oncology literature. Finally, we provide recommendations regarding best practices for identifying and managing multicollinearity in complex data sets.