Wan Chan Tseung, Hok Seum

Background: Post-treatment prognosis and monitoring are critical for determining the timing of salvage treatment in glioblastoma patients but has been challenging due to difficulties differentiating progression from treatment effects in conventional images. This exploratory study aimed to establish the correlation of radiomics image features from time series of amino acid tracer 18F-DOPA PET images, with outcomes, using machine learning and dimension reduction analysis. Methods: 18F-DOPA PET images were collected for a patient cohort with wild-type IDH and unmethylated MGMT who underwent dose-escalated radiation therapy. Quantitative features were derived from the high uptake region (T/N > 2.0) in pre- and post-radiation therapy follow-up 18F-DOPA PET images. A customized workflow was utilized for pre-selecting predictive features, followed by manifold learning. Machine learning algorithms were employed to establish associations between imaging features and remaining survival (RS), defined as the time between a follow-up scan and date of death. Results: The ML models exhibited 81–83% ROC_AUC in predicting RS evaluated on an independent test dataset. A RS map is proposed for monitoring tumor alterations through serial 18F-DOPA PET scans, demonstrating superior sensitivity and better correlation with survival compared to the RANO criteria. Conclusions: Our study demonstrates that ML models utilizing FU 18F-DOPA PET images have the potential to effectively predict future survival outcomes in patients with glioblastoma treated with dose-escalated radiation therapy. The capability to assess changes in tumor over time through imaging can potentially assist in patient stratification and the selection of salvage treatments, while also aiding in distinguishing treatment effects from genuine tumor progression.

AbstractPurposeDental fixtures are commonplace in an aging, radiation treatment population. The current, local standard of practice in particle therapy is to employ treatment geometries to avoid delivery through implanted dental fixtures. The present study aims to observe the physical effect of delivering therapeutic proton beams through common dental fixture materials as prelude to an eventual goal of assessing the feasibility of using treatment geometries not specified for avoidance of oral implants. A sampling of common dental materials was selected based on prosthodontic consult and was evaluated in terms of relative stopping power and three‐dimensional (3D) dose perturbation.MethodsAmalgams, porcelain‐fused‐to‐metal (PFM) crowns consisting of zirconia and non‐noble base metals, and lithium disilicate implants were chosen for analysis. Theoretical stopping power (S) and mass stopping power (S/ρ) were calculated using the Stopping and Range of Ions in Matter (SRIM) application, basing stoichiometric compositions of each fixture on published materials data. S and S/ρ were calculated for a range of historically available compositions of amalgams from 1900 until the current era. The perturbance of S and S/ρ as a function of clinically relevant ranges of amalgam compositions for the modern era was analyzed. Water equivalent thickness (WET) and relative stopping power (Srel) of each material was measured for a clinical spot‐scanning proton beam with monoenergies of 159.9 and 228.8 MeV with a multi‐layer ionization chamber (MLIC). Subsequently, 3D dose perturbation was assessed by delivering proton beams through a custom phantom designed to simulate both en‐face and on‐edge treatment geometries through the selected materials. A treatment plan mimicking the experimental delivery was constructed in the institutional treatment planning system and calculated using TOPAS‐based Monte Carlo simulation (MCS). Experimental results were used to validate the MCS. Finally, treatment planning system (TPS) outputs were compared to MCS to determine the accuracy of the dose calculation model.ResultsHistorical compositions of amalgams ranged in S from 44.8 to 42.9 MeV/cm, with the greatest deviation being observed for the 1900–1959 era. Deviation as a function of amalgam composition from the modern era was most sensitive to proportion of Hg, accounting for deviations up to ‐4.2% at the greatest clinically relevant concentration. S/ρ was not found to vary greatly between each porcelain and metal alloy material for PFM type crowns. Relative stopping powers ranged between 1.3 and 5.4 for all studied materials, suggesting substantial changes in proton range with respect to water. Film measurements of pristine spots confirm dose perturbance and shortening of proton range, with an upstream shift of each Bragg peak being observed directly behind the installed fixture. At high energies, cold spots were found in all cases directly behind each material feature with a medial fill‐in of dose occurring distally. Qualitative agreement of spot perturbance was confirmed between film measurements and MCS. Finally, when comparing integrated depth doses (IDD) by summing over all axial directions, good agreement is observed between TPS and MCS.ConclusionsAll dental materials studied substantially perturbed the dosimetry of pristine proton spots both in terms of WET/Srel as well as the spatial distribution of dose. Proton range was quantifiably shortened, and each dental material affected a cold spot directly behind the object with medial dose back‐filling was observed distally. MCS and Eclipse dose calculations exhibited good agreement with measurements, suggesting that treatment planning without employing avoidance strategies may be possible with further investigation.

PurposeThe presence of respiratory motion during radiation treatment leads to degradation of the expected dose distribution, both for target coverage and healthy tissue sparing, particularly for techniques like pencil beam scanning proton therapy which have dynamic delivery systems. While tools exist to estimate this degraded four‐dimensional (4D) dose, they typically have one or more deficiencies such as not including the particular effects from a dynamic delivery, using analytical dose calculations, and/or using nonphysical dose‐accumulation methods. This work presents a clinically useful 4D‐dose calculator that addresses each of these shortcomings.MethodsTo quickly compute the 4D dose, the three main tasks of the calculator were run on graphics processing units (GPUs). These tasks were (a) simulating the delivery of the plan using measured delivery parameters to distribute the plan amongst 4DCT phases characterizing the patient breathing, (b) using an in‐house Monte Carlo simulation (MC) dose calculator to determine the dose delivered to each breathing phase, and (c) accumulating the doses from the various breathing phases onto a single phase for evaluation. The accumulation was performed by individually transferring the energy and mass of dose‐grid subvoxels, a technique that models the transfer of dose in a more physically realistic manner. The calculator was run on three test cases, with lung, esophagus, and liver targets, respectively, to assess the various uncertainties in the beam delivery simulation as well as to characterize the dose‐accumulation technique.ResultsFour‐dimensional doses were successfully computed for the three test cases with computation times ranging from 4–6 min on a server with eight NVIDIA Titan X graphics cards; the most time‐consuming component was the MC dose engine. The subvoxel‐based dose‐accumulation technique produced stable 4D‐dose distributions at subvoxel scales of 0.5–1.0 mm without impairing the total computation time. The uncertainties in the beam delivery simulation led to moderate variations of the dose–volume histograms for these cases; the variations were reduced by implementing repainting or phase‐gating motion mitigation techniques in the calculator.ConclusionsA MC‐based and GPU‐accelerated 4D‐dose calculator was developed to estimate the effects of respiratory motion on pencil beam scanning proton therapy treatments. After future validation, the calculator could be used to assess treatment plans and its quick runtime would make it easily usable in a future 4D‐robust optimization system.

PurposeAccuracy of dose calculation models and robustness under various uncertainties are key factors influencing the quality of intensity modulated proton therapy (IMPT) plans. To mitigate the effects of uncertainties and to improve the dose calculation accuracy, an all‐scenario robust IMPT optimization based on accurate Monte Carlo (MC) dose calculation was developed.MethodsIn the all‐scenario robust IMPT optimization, dose volume histograms (DVHs) were computed for the nominal case and for each uncertainty scenario. All scenarios were weighted by DVH values dynamically throughout optimization iterations. In contrast, probabilistic approach weighted scenarios with fixed scenario weights and worst case optimizations picked one single scenario — the worst scenario for each iteration. Uncertainties in patient setup and proton range were considered in all clinical cases studied. Graphics processing unit (GPU) computation was employed to reduce the computational time in both the MC dose calculation and optimization stages. A previously published adaptive quasi‐Newton method for proton optimization was extended to include robustness. To validate the all‐scenario algorithm extension, it was compared with the single scenario optimization target volume (OTV) based approach in clinical cases of three different disease sites. Additional comparisons with worst case optimization methods were conducted to evaluate the performance of the all‐scenario robust optimization against other robust optimizations.ResultsThe all‐scenario robust IMPT optimization spared organs at risk (OARs) better than the OTV‐based method while maintaining target coverage and improving the robustness of targets and OARs. Compared with composite and voxel‐wise worst case optimization, the all‐scenario robust optimization converged faster, and arrived at solutions of tighter DVH robustness spread, better target coverage and lower OAR dose.ConclusionAn all‐scenario robust IMPT treatment planning system was developed using an adaptive quasi‐Newton optimization method. The optimization system was GPU‐accelerated and based on MC dose calculation. Improved performance was observed in clinical cases when compared with worst case optimization methods.