Patel, Samir H.

Abstract In radiotherapy, 2D orthogonally projected kV images are used for patient alignment when 3D-on-board imaging(OBI) unavailable. But tumor visibility is constrained due to the projection of patient’s anatomy onto a 2D plane, potentially leading to substantial setup errors. In treatment room with 3D-OBI such as cone beam CT(CBCT), the field of view(FOV) of CBCT is limited with unnecessarily high imaging dose, thus unfavorable for pediatric patients. A solution to this dilemma is to reconstruct 3D CT from kV images obtained at the treatment position. Here, we propose a dual-models framework built with hierarchical ViT blocks. Unlike a proof-of-concept approach, our framework considers kV images as the solo input and can synthesize accurate, full-size 3D CT in real time(within milliseconds). We demonstrate the feasibility of the proposed approach on 10 patients with head and neck(H\&N) cancer using image quality(MAE:<45HU), dosimetric accuracy(Gamma passing rate(2%/2mm/10%):>97%) and patient position uncertainty(shift error:<0.4mm). The proposed framework can generate accurate 3D CT faithfully mirroring real-time patient position, thus significantly improving patient setup accuracy, keeping imaging dose minimum, and maintaining treatment veracity.

Background: Approximately 75% of all head and neck cancer patients are treated with radiotherapy (RT). RT to the oral cavity results in acute and late adverse events which can be severe and detrimental to a patient’s quality of life and function. The purpose of this study was to explore associations between RT dose to a defined oral cavity organ-at-risk (OAR) avoidance structure, provider- and patient-reported outcomes (PROs), opioid use, and hospitalization. Methods: This was a retrospective analysis of prospectively obtained outcomes using multivariable modeling. The study included 196 patients treated with RT involving the oral cavity for a head and neck tumor. A defined oral cavity OAR avoidance structure was used in all patients for RT treatment planning. Validated PROs were collected prospectively. Opioid use and hospitalization were abstracted electronically from medical records. Results: Multivariable modeling revealed the mean dose to the oral cavity OAR was significantly associated with opioid use (p = 0.0082) and hospitalization (p = 0.0356) during and within 30 days of completing RT. Conclusions: The findings of this study may be valuable in RT treatment planning for patients with tumors of the head and neck region to reduce the need for opioid use and hospitalization during treatment.

Abstract Purpose . To enhance an in-house graphic-processing-unit accelerated virtual particle (VP)-based Monte Carlo (MC) proton dose engine (VPMC) to model aperture blocks in both dose calculation and optimization for pencil beam scanning proton therapy (PBSPT)-based stereotactic radiosurgery (SRS). Methods and materials . A module to simulate VPs passing through patient-specific aperture blocks was developed and integrated in VPMC based on simulation results of realistic particles (primary protons and their secondaries). To validate the aperture block module, VPMC was first validated by an opensource MC code, MCsquare, in eight water phantom simulations with 3 cm thick brass apertures: four were with aperture openings of 1, 2, 3, and 4 cm without a range shifter, while the other four were with same aperture opening configurations with a range shifter of 45 mm water equivalent thickness. Then, VPMC was benchmarked with MCsquare and RayStation MC for 10 patients with small targets (average volume 8.4 c.c. with range of 0.4–43.3 c.c.). Finally, 3 typical patients were selected for robust optimization with aperture blocks using VPMC. Results . In the water phantoms, 3D gamma passing rate (2%/2 mm/10%) between VPMC and MCsquare was 99.71 ± 0.23%. In the patient geometries, 3D gamma passing rates (3%/2 mm/10%) between VPMC/MCsquare and RayStation MC were 97.79 ± 2.21%/97.78 ± 1.97%, respectively. Meanwhile, the calculation time was drastically decreased from 112.45 ± 114.08 s (MCsquare) to 8.20 ± 6.42 s (VPMC) with the same statistical uncertainties of ~0.5%. The robustly optimized plans met all the dose–volume-constraints (DVCs) for the targets and OARs per our institutional protocols. The mean calculation time for 13 influence matrices in robust optimization by VPMC was 41.6 s and the subsequent on-the-fly ‘trial-and-error’ optimization procedure took only 71.4 s on average for the selected three patients. Conclusion . VPMC has been successfully enhanced to model aperture blocks in dose calculation and optimization for the PBSPT-based SRS.

Abstract Introduction Historical reservations regarding stereotactic radiosurgery (SRS) for small-cell lung cancer (SCLC) brain metastases include concerns for short-interval and diffuse central nervous system (CNS) progression, poor prognoses, and increased neurological mortality specific to SCLC histology. We compared SRS outcomes for SCLC and non-small cell lung cancer (NSCLC) where SRS is well established. Methods Multicenter first-line SRS outcomes for SCLC and NSCLC from 2000 to 2022 were retrospectively collected (n = 892 SCLC, n = 4785 NSCLC). Data from the prospective Japanese Leksell Gamma Knife Society (JLGK0901) clinical trial of first-line SRS were analyzed as a comparison cohort (n = 98 SCLC, n = 814 NSCLC). Overall survival (OS) and CNS progression were analyzed using Cox proportional hazard and Fine-Gray models, respectively, with multivariable adjustment for cofactors including age, sex, performance status, year, extracranial disease status, and brain metastasis number and volume. Mutation-stratified analyses were performed in propensity score–matched retrospective cohorts of epidermal growth factor receptor (EGFR) and anaplastic lymphoma kinase (ALK) positive NSCLC, mutation-negative NSCLC, and SCLC. Results OS was superior for patients with NSCLC compared to SCLC in the retrospective dataset (median OS = 10.5 vs 8.6 months; P &lt; .001) and in the JLGK0901 dataset. Hazard estimates for first CNS progression favoring NSCLC were similar in both datasets but reached statistical significance in the retrospective dataset only (multivariable hazard ratio = 0.82, 95% confidence interval = 0.73 to 0.92, P = .001). In the propensity score–matched cohorts, there were continued OS advantages for NSCLC patients (median OS = 23.7 [EGFR and ALK positive NSCLC] vs 13.6 [mutation-negative NSCLC] vs 10.4 months [SCLC], pairwise P values &lt; 0.001), but no statistically significant differences in CNS progression were observed in the matched cohorts. Neurological mortality and number of lesions at CNS progression were similar for NSCLC and SCLC patients. Leptomeningeal progression was increased in patients with NSCLC compared to SCLC in the retrospective dataset only (multivariable hazard ratio = 1.61, 95% confidence interval = 1.14 to 2.26, P = .007). Conclusions After SRS, SCLC histology was associated with shorter OS compared to NSCLC. CNS progression occurred earlier in SCLC patients overall but was similar in patients matched on baseline factors. SCLC was not associated with increased neurological mortality, number of lesions at CNS progression, or leptomeningeal progression compared to NSCLC. These findings may better inform clinical expectations and individualized decision making regarding SRS for SCLC patients.

Abstract OBJECTIVE Immune checkpoint inhibitors (ICI) and stereotactic radiosurgery (SRS) are commonly utilized in the management of brain metastases. Treatment-related imaging changes (TRIC) are a frequently observed clinical manifestation and are commonly classified as radiographic radiation necrosis. However, these findings are not well characterized and may predict for response to SRS and ICI. METHODS The diagnosis-specific graded prognostic assessment was used to guide variable selection. TRIC were determined based upon MRI, PET/CT, or MR spectroscopy and a consensus by local clinical providers was required. RESULTS The analysis included 697 patients with 4,536 brain metastases across 11 institutions in 4 countries. The median follow-up after SRS was 13.6 months. The median age was 66 years, 54.1% of patients were male, and 57.3%, 36.3%, and 6.4% were non-small cell lung cancer, melanoma, and renal cell carcinoma (RCC) histology, respectively. TRIC were observed in 9.8%. On univariable analysis, Karnofsky Performance Status (KPS) (hazard ratio [HR]: 0.98; p &lt; 0.001), presence of TRIC (HR: 0.67; p = 0.03), female sex (HR: 0.67; p &lt; 0.001), and prior resection (HR: 0.60; p = 0.03) were associated with improved OS. On multivariable analysis, KPS (HR: 0.98; p &lt; 0.001) and the presence of TRIC (HR: 0.66; p = 0.03) were associated with improved OS. A V12 Gy ≥ 10 cm3 (Odds Ratio [OR]: 2.78; p &lt; 0.001), prior whole brain radiation therapy (OR: 3.46; p = 0.006), and RCC histology (OR: 3.10; p = 0.01) were associated with an increased probability of developing TRIC. The median OS in patients with and without TRIC was 29.0 and 23.1 months, respectively (log-rank p = 0.03). CONCLUSION TRIC following ICI and SRS are associated with a median OS benefit of approximately 6 months. Further prospective study is warranted to further elucidate the role and etiology of this common clinical scenario.

AbstractBackgroundWe seek to inform radiotherapy (RT) delivery for adenoid cystic carcinoma of the head and neck (ACC) by evaluating RT techniques and recurrence patterns.MethodsWe identified patients with ACC treated with curative‐intent RT from 2005 to 2021. Imaging was reviewed to determine local recurrence (LR).ResultsNinety‐one patients were included. The 5‐year LR risk was 12.2% (6.6–22.7). One patient each experienced a marginal and out‐of‐field recurrence. Patients receiving >60 Gy postoperatively had a 5‐year LR risk of 0% compared to 10.7% (4.2–27.2) with ≤60 Gy. Those receiving 70 and <70 Gy definitively had a 5‐year LR risk of 15.2% (2.5–91.6) and 33.3% (6.7–100.0), respectively. No patients had regional nodal failure.ConclusionsModern, conformal RT for ACC results in low rates of LR. Doses >60 and 70 Gy may improve control in the postoperative and definitive settings, respectively. Elective nodal treatment can be omitted in well‐selected patients.

AbstractBackgroundIn 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.PurposeA 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.MethodsNeutrons 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.ResultsIntegrated 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.ConclusionVPMC was found to achieve high accuracy and efficiency in proton therapy dose calculation.

AbstractBackgroundBoth 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.PurposeTo perform exploratory investigation of the dose and LET effects upon AE initialization in PBS using seed spots analysis.MethodsA 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.ResultsThe 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.ConclusionOur 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.

AbstractPurposeTo 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 MaterialsA 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.ResultsFor 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.ConclusionThe proposed online ART workflow for PBS was demonstrated to be efficient and effective by generating a reoptimized plan that significantly improved the plan quality.

Abstract Purpose To evaluate the accuracy of the RayStation Monte Carlo dose engine (RayStation MC) in modeling small‐field block apertures in proton pencil beam scanning. Furthermore, we evaluate the suitability of MCsquare as a second check for RayStation MC. Methods We have enhanced MCsquare to model block apertures. To test the accuracy of both RayStation MC and the newly enhanced MCsquare, we compare the dose predictions of each to in‐water dose measurements obtained using diode detectors and radiochromic film. Nine brass apertures with openings of 1, 2, 3, 4, and 5 cm and either 2 cm or 4 cm thickness were used in the irradiation of a water phantom. Two measurement setups were used, one with a range shifter and 119.7 MeV proton beam energy and the other with no range shifter and 147 MeV proton beam energy. To further test the validity of RayStation MC and MCsquare in modeling block apertures and to evaluate MCsquare as a second check tool, 10 small‐field (average target volume 8.3 cm 3 ) patient treatment plans were calculated by each dose engine followed by a statistical comparison. Results Comparing to the absolute dose measurements in water, RayStation MC differed by 1.2% ± 1.0% while MCsquare differed by −1.8% ± 3.7% in the plateau region of a pristine Bragg peak. Compared to the in‐water film measurements, RayStation MC and MCsquare both performed well with an average 2D‐3D gamma passing rate of 99.4% and 99.7% (3%/3 mm), respectively. A t ‐test comparing the agreement with the film measurements between RayStation MC and MCsquare suggested that the relative spatial dose distributions calculated by MCsquare and RayStation MC were statistically indistinguishable. Directly comparing the dose calculations between MCsquare and RayStation MC over 10 patients resulted in an average 3D‐3D gamma passing rates of 98.5% (3%/3 mm) and 94.1% (2%/2 mm), respectively. Conclusion The validity of RayStation MC algorithm for use with patient‐specific apertures has been expanded to include small apertures. MCsquare has been enhanced to model apertures and was found to be an adequate second check of RayStation MC in this scenario.