Investigation of intra-fractionated range guided adaptive proton therapy (RGAPT): II. Range-shift compensated on-line treatment adaptation and verification.

We previously proposed range-guided adaptive proton therapy (RGAPT) that uses mid-range treatment beams as probing beams and intra-fractionated range measurements for online adaptation. In this work, we demonstrated experimental verification and reported the dosimetric accuracy for RGAPT. A STEEV phantom was used for the experiments, and a 3 × 3 × 3 cm3cube inside the phantom was assigned to be the treatment target. We simulated three online range shift scenarios: reference, overshoot, and undershoot, by placing upstream Lucite sheets, 4, 0, and 8 that corresponded to changes of 0, 6.8, and -6.8 mm, respectively, in water-equivalent path length. The reference treatment plan was to deliver single-field uniform target doses in pencil beam scanning mode and generated on the Eclipse treatment planning system. Different numbers of mid-range layers, including single, three, and five layers, were selected as probing beams to evaluate beam range (BR) measurement accuracy in positron emission tomography (PET). Online plans were modified to adapt to BR shifts and compensate for probing beam doses. In contrast, non-adaptive plans were also delivered and compared to adaptive plans by film measurements. The mid-range probing beams of three (5.55MU) and five layers (8.71MU) yielded accurate range shift measurements in 60 s of PET acquisition with uncertainty of 0.5 mm while the single-layer probing (1.65MU) was not sufficient for measurements. The adaptive plans achieved an average gamma (2%/2 mm) passing rate of 95%. In contrast, the non-adaptive plans only had an average passing rate of 69%. RGAPT planning and delivery are feasible and verified by the experiments. The probing beam delivery, range measurements, and adaptive planning and delivery added a small increase in treatment delivery workflow time but resulted in substantial dose improvement. The three-layer mid-range probing was most suitable considering the balance of high range measurement accuracy and the low number of probing beam layers.

Investigation of intra-fractionated range guided adaptive proton therapy: I. On-line PET imaging and range measurement.

Objective.Develop a prototype on-line positron emission tomography (PET) scanner and evaluate its capability of on-line imaging and intra-fractionated proton-induced radioactivity range measurement.Approach.Each detector consists of 32 × 32 array of 2 × 2 × 30 mm3Lutetium-Yttrium Oxyorthosilicate scintillators with single-scintillator-end readout through a 20 × 20 array of 3 × 3 mm2Silicon Photomultipliers. The PET can be configurated with a full-ring of 20 detectors for conventional PET imaging or a partial-ring of 18 detectors for on-line imaging and range measurement. All detector-level readout and processing electronics are attached to the backside of the system gantry and their output signals are transferred to a field-programable-gate-array based system electronics and data acquisition that can be placed 2 m away from the gantry. The PET imaging performance and radioactivity range measurement capability were evaluated by both the offline study that placed a radioactive source with known intensity and distribution within a phantom and the online study that irradiated a phantom with proton beams under different radiation and imaging conditions.Main results.The PET has 32 cm diameter and 6.5 cm axial length field-of-view (FOV), ∼2.3-5.0 mm spatial resolution within FOV, 3% sensitivity at the FOV center, 18%-30% energy resolution, and ∼9 ns coincidence time resolution. The offline study shows the PET can determine the shift of distal falloff edge position of a known radioactivity distribution with the accuracy of 0.3 ± 0.3 mm even without attenuation and scatter corrections, and online study shows the PET can measure the shift of proton-induced positron radioactive range with the accuracy of 0.6 ± 0.3 mm from the data acquired with a short-acquisition (60 s) and low-dose (5 MU) proton radiation to a human head phantom.Significance.This study demonstrated the capability of intra-fractionated PET imaging and radioactivity range measurement and will enable the investigation on the feasibility of intra-fractionated, range-shift compensated adaptive proton therapy.

Intensity modulated proton arc therapy via geometry-based energy selection for ependymoma.

PURPOSE: We developed and tested a novel method of creating intensity modulated proton arc therapy (IMPAT) plans that uses computing resources similar to those for regular intensity-modulated proton therapy (IMPT) plans and may offer a dosimetric benefit for patients with ependymoma or similar tumor geometries. METHODS: Our IMPAT planning method consists of a geometry-based energy selection step with major scanning spot contributions as inputs computed using ray-tracing and single-Gaussian approximation of lateral spot profiles. Based on the geometric relation of scanning spots and dose voxels, our energy selection module selects a minimum set of energy layers at each gantry angle such that each target voxel is covered by sufficient scanning spots as specified by the planner, with dose contributions above the specified threshold. Finally, IMPAT plans are generated by robustly optimizing scanning spots of the selected energy layers using a commercial proton treatment planning system (TPS). The IMPAT plan quality was assessed for four ependymoma patients. Reference three-field IMPT plans were created with similar planning objective functions and compared with the IMPAT plans. RESULTS: In all plans, the prescribed dose covered 95% of the clinical target volume (CTV) while maintaining similar maximum doses for the brainstem. While IMPAT and IMPT achieved comparable plan robustness, the IMPAT plans achieved better homogeneity and conformity than the IMPT plans. The IMPAT plans also exhibited higher relative biological effectiveness (RBE) enhancement than did the corresponding reference IMPT plans for the CTV in all four patients and brainstem in three of them. CONCLUSIONS: The proposed method demonstrated potential as an efficient technique for IMPAT planning and may offer a dosimetric benefit for patients with ependymoma or tumors in close proximity to critical organs. IMPAT plans created using this method had elevated RBE enhancement associated with increased linear energy transfer (LET) in both targets and abutting critical organs.

Our Experience Leading a Large Medical Physics Practice During the COVID-19 Pandemic.

PURPOSE: To provide a series of suggestions for other Medical Physics practices to follow in order to provide effective radiation therapy treatments during the COVID-19 pandemic. METHODS AND MATERIALS: We reviewed our entire Radiation Oncology infrastructure to identify a series of workflows and policy changes that we implemented during the pandemic that yielded more effective practices during this time. RESULTS: We identified a structured list of several suggestions that can help other Medical Physics practices overcome the challenges involved in delivering high quality radiotherapy services during this pandemic. CONCLUSIONS: Our facility encompasses 4 smaller Houston Area Locations (HALs), a main campus with 8 distinct services based on treatment site (ie. Thoracic, Head and Neck, Breast, Gastrointestinal, Gynecology, Genitourinary, Hematologic Malignancies, Melanoma and Sarcoma and Central Nervous System/Pediatrics), a Proton Center facility, an MR-Linac, a Gamma Knife clinic and an array of brachytherapy services. Due to the scope of our services, we have gained experience in dealing with the rapidly changing pandemic effects on our clinical practice. Our paper provides a resource to other Medical Physics practices in search of workflows that have been resilient during these challenging times.

Patterns of protein expression in human head and neck cancer cell lines differ after proton vs photon radiotherapy.

BACKGROUND: Proton radiotherapy (PRT) may be a less toxic alternative to photon radiotherapy (XRT) for patients with head and neck squamous cell carcinoma (HNSCC). However, the molecular responses of HNSCC cells to PRT vs XRT are unclear. METHODS: Proteomics analyses of protein expression profiles by reverse-phase protein arrays were done for two human papillomavirus [HPV]-negative and two HPV+ cell lines. Expression patterns of 175 proteins involved in several signaling pathways were tested. RESULTS: Compared with PRT, XRT tended to induce lower expression of DNA damage repair-and cell cycle arrest-related proteins and higher expression of cell survival- and proliferation-related proteins. CONCLUSIONS: Under these experimental conditions, PRT and XRT induced different protein expression and activation profiles. Further preclinical verification is needed, as are studies of tumor pathway mutations as biomarkers for choice of treatment or as radiosensitization targets to improve the response of HNSCC to PRT or XRT.

Quantifying the accuracy of deformable image registration for cone-beam computed tomography with a physical phantom.

PURPOSE: Kilo-voltage cone-beam computed tomography (CBCT) is widely used for patient alignment, contour propagation, and adaptive treatment planning in radiation therapy. In this study, we evaluated the accuracy of deformable image registration (DIR) for CBCT under various imaging protocols with different noise and patient dose levels. METHODS: A physical phantom previously developed to facilitate end-to-end testing of the DIR accuracy was used with Varian Velocity v4.0 software to evaluate the performance of image registration from CT to CT, CBCT to CT, and CBCT to CBCT. The phantom is acrylic and includes several inserts that simulate different tissue shapes and properties. Deformations and anatomic changes were simulated by changing the rotations of both the phantom and the inserts. CT images (from a head and neck protocol) and CBCT images (from pelvis, head and "Image Gently" protocols) were obtained with different image noise and dose levels. Large inserts were filled with Mobil DTE oil to simulate soft tissue, and small inserts were filled with bone materials. All inserts were contoured before the DIR process to provide a ground truth contour size and shape for comparison. After the DIR process, all deformed contours were compared with the originals using Dice similarity coefficient (DSC) and mean distance to agreement (MDA). Both large and small volume of interests (VOIs) for DIR volume selection were tested by simulating a DIR process that included whole patient image volume and clinical target volumes (CTV) only (for CTVs propagation). RESULTS: For cross-modality DIR registration (CT to CBCT), the DSC were >0.8 and the MDA were <3 mm for CBCT pelvis, and CBCT head protocols. For CBCT to CBCT and CT to CT, the DIR accuracy was improved relative to the cross-modality tests. For smaller VOIs, the DSC were >0.8 and MDA <2 mm for all modalities. CONCLUSIONS: The accuracy of DIR depends on the quality of the CBCT image at different dose and noise levels.

3D treatment planning system-Varian Eclipse for proton therapy planning.

The capabilities of the Eclipse treatment planning system (TPS) (Varian Medical Systems, Palo Alto, CA) for proton therapy treatment planning are described. Various steps involved in the planning process to produce a 3-dimensional (3D) dose distribution both for the passive scattering and pencil beam scanning proton beam therapy are outlined. Mitigation of range and setup uncertainties through robust optimization is discussed. Use of verification plans for patient treatment field dosimetry quality assurance (QA) is presented. Limitations of the Eclipse TPS and future developments are discussed.

Technical Note: Dosimetric characteristics of the ocular beam line and commissioning data for an ocular proton therapy planning system at the Proton Therapy Center Houston.

PURPOSE: To systematically analyze and present the properties of a small-field, double-scattering proton beam line intended to be used for the treatment of ocular cancer, and to provide configuration data for commission of the Eclipse Ocular Proton Planning System. METHODS: Measurements were made using ionization chambers, diodes, and films to determine dose profiles and output factors of the proton beams of the beam line at the Proton Therapy Center Houston. In parallel, Monte Carlo simulations were performed to validate the measured data and to provide additional insight into detailed beam parameters that are hard to measure, such as field size factors and a comparison of output factors as a function of circular and rectangular fields. RESULTS: The presented data comprise depth dose profiles, including distal and proximal profiles used to configure the Eclipse Ocular Proton Planning system, distal fall-off widths, lateral profiles and penumbrae sizes, as well as output factors as a function of field size, SOBP width, range shifter thickness, snout position, and source-to-surface distance. CONCLUSIONS: We have completed a comprehensive characterization of the beam line. The data will be useful to characterize proton beams in clinical and experimental small-field applications.

Consensus Guidelines for Implementing Pencil-Beam Scanning Proton Therapy for Thoracic Malignancies on Behalf of the PTCOG Thoracic and Lymphoma Subcommittee.

Pencil-beam scanning (PBS) proton therapy (PT), particularly intensity modulated PT, represents the latest advanced PT technology for treating cancers, including thoracic malignancies. On the basis of virtual clinical studies, PBS-PT appears to have great potential in its ability to tightly tailor the dose to the target while sparing critical structures, thereby reducing treatment-related toxicities, particularly for tumors in areas with complicated anatomy. However, implementing PBS-PT for moving targets has several additional technical challenges compared with intensity modulated photon radiation therapy or passive scattering PT. Four-dimensional computed tomography-based motion management and robust optimization and evaluation are crucial for minimizing uncertainties associated with beam range and organ motion. Rigorous quality assurance is required to validate dose delivery both before and during the course of treatment. Active motion management (eg, breath hold), beam gating, rescanning, tracking, or adaptive planning may be needed for cases involving significant motion or changes in motion or anatomy over the course of treatment.

Intensity-Modulated Proton Therapy Adaptive Planning for Patients with Oropharyngeal Cancer.

PURPOSE: The authors aimed to illustrate the potential dose differences to clinical target volumes (CTVs) and organs-at-risk (OARs) volumes after proton adaptive treatment planning was used. PATIENTS AND METHODS: The records of 10 patients with oropharyngeal cancer were retrospectively reviewed. Each patient's treatment plan was generated by using the Eclipse treatment planning system. Verification computed tomography (CT) scan was performed during the fourth week of treatment. Deformable image registrations were performed between the 2 CT image sets, and the CTVs and major OARs were transferred to the verification CT images to generate the adaptive plan. We compared the accumulated doses to CTVs and OARs between the original and adaptive plans, as well as between the adaptive and verification plans to simulate doses that would have been delivered if the adaptive plans were not used. RESULTS: Body contours were different on planning and week-4 verification CTs. Mean volumes of all CTVs were reduced by 4% to 8% (P ≤ .04), and the volumes of left and right parotid glands also decreased (by 11% to 12%, P ≤ .004). Brainstem and oral cavity volumes did not significantly differ (all P ≥ .14). All mean doses to the CTV were decreased for up to 7% (P ≤ .04), whereas mean doses to the right parotid and oral cavity increased from a range of 5% to 8% (P ≤ .03), respectively. CONCLUSION: Verification and adaptive planning should be recommended during the course of proton therapy for patients with head and neck cancer to ensure adequate dose deliveries to the planned CTVs, while safe doses to OARs can be respected.

Consensus Statement on Proton Therapy in Early-Stage and Locally Advanced Non-Small Cell Lung Cancer.

Radiation dose escalation has been shown to improve local control and survival in patients with non-small cell lung cancer in some studies, but randomized data have not supported this premise, possibly owing to adverse effects. Because of the physical characteristics of the Bragg peak, proton therapy (PT) delivers minimal exit dose distal to the target volume, resulting in better sparing of normal tissues in comparison to photon-based radiation therapy. This is particularly important for lung cancer given the proximity of the lung, heart, esophagus, major airways, large blood vessels, and spinal cord. However, PT is associated with more uncertainty because of the finite range of the proton beam and motion for thoracic cancers. PT is more costly than traditional photon therapy but may reduce side effects and toxicity-related hospitalization, which has its own associated cost. The cost of PT is decreasing over time because of reduced prices for the building, machine, maintenance, and overhead, as well as newer, shorter treatment programs. PT is improving rapidly as more research is performed particularly with the implementation of 4-dimensional computed tomography-based motion management and intensity modulated PT. Given these controversies, there is much debate in the oncology community about which patients with lung cancer benefit significantly from PT. The Particle Therapy Co-operative Group (PTCOG) Thoracic Subcommittee task group intends to address the issues of PT indications, advantages and limitations, cost-effectiveness, technology improvement, clinical trials, and future research directions. This consensus report can be used to guide clinical practice and indications for PT, insurance approval, and clinical or translational research directions.

Is there a clinical benefit with a smooth compensator design compared with a plunged compensator design for passive scattered protons?

In proton therapy, passive scattered proton plans use compensators to conform the dose to the distal surface of the planning volume. These devices are custom made from acrylic or wax for each treatment field using either a plunge-drilled or smooth-milled compensator design. The purpose of this study was to investigate if there is a clinical benefit of generating passive scattered proton radiation treatment plans with the smooth compensator design. We generated 4 plans with different techniques using the smooth compensators. We chose 5 sites and 5 patients for each site for the range of dosimetric effects to show adequate sample. The plans were compared and evaluated using multicriteria (MCA) plan quality metrics for plan assessment and comparison using the Quality Reports [EMR] technology by Canis Lupus LLC. The average absolute difference for dosimetric metrics from the plunged-depth plan ranged from -4.7 to +3.0 and the average absolute performance results ranged from -6.6% to +3%. The manually edited smooth compensator plan yielded the best dosimetric metric, +3.0, and performance, + 3.0% compared to the plunged-depth plan. It was also superior to the other smooth compensator plans. Our results indicate that there are multiple approaches to achieve plans with smooth compensators similar to the plunged-depth plans. The smooth compensators with manual compensator edits yielded equal or better target coverage and normal tissue (NT) doses compared with the other smooth compensator techniques. Further studies are under investigation to evaluate the robustness of the smooth compensator design.

In-beam PET imaging for on-line adaptive proton therapy: an initial phantom study.

We developed and investigated a positron emission tomography (PET) system for use with on-line (both in-beam and intra-fraction) image-guided adaptive proton therapy applications. The PET has dual rotating depth-of-interaction measurable detector panels by using solid-state photomultiplier (SSPM) arrays and LYSO scintillators. It has a 44 mm diameter trans-axial and 30 mm axial field-of-view (FOV). A 38 mm diameter polymethyl methacrylate phantom was placed inside the FOV. Both PET and phantom axes were aligned with a collimated 179.2 MeV beam. Each beam delivered ∼50 spills (0.5 s spill and 1.5 s inter-spill time, 3.8 Gy at Bragg peak). Data from each beam were acquired with detectors at a given angle. Nine datasets for nine beams with detectors at nine different angles over 180° were acquired for full-tomographic imaging. Each dataset included data both during and 5 min after irradiations. The positron activity-range was measured from the PET image reconstructed from all nine datasets and compared to the results from simulated images. A (22)Na disc-source was also imaged after each beam to monitor the PET system's performance. PET performed well except for slight shifts of energy photo-peak positions (<1%) after each beam, due mainly to the neutron exposure of SSPM that increased the dark-count noise. This minor effect was corrected offline with a shifting 350-650 keV energy window for each dataset. The results show a fast converging of activity-ranges measured by the prototype PET with high sensitivity and uniform resolution. Sub-mm activity-ranges were achieved with minimal 6 s acquisition time and three spill irradiations. These results indicate the feasibility of PET for intra-fraction beam-range verification. Further studies are needed to develop and apply a novel clinical PET system for on-line image-guided adaptive proton therapy.

Parameterization of multiple Bragg curves for scanning proton beams using simultaneous fitting of multiple curves.

Although Bortfeld's analytical formula is useful for describing Bragg curves, measured data can deviate from the values predicted by the model. Thus, we sought to determine the parameters of a closed analytical expression of multiple Bragg curves for scanning proton pencil beams using a simultaneous optimization algorithm and to determine the minimum number of energies that need to be measured in treatment planning so that complete Bragg curves required by the treatment planning system (TPS) can be accurately predicted. We modified Bortfeld's original analytical expression of Bragg curves to accurately describe the dose deposition resulting from secondary particles. The parameters of the modified analytical expression were expressed as the parabolic cylinder function of the ranges of the proton pencil beams in water. Thirty-nine discrete Bragg curves were measured in our center using a PTW Bragg Peak chamber during acceptance and commission of the scanning beam proton delivery system. The coefficients of parabolic function were fitted by applying a simultaneous optimization algorithm to seven measured curves. The required Bragg curves for 45 energies in the TPS were calculated using our parameterized analytical expression. Finally, the 10 cm width of spread-out Bragg peaks (SOBPs) of beams with maximum energies of 221.8 and 121.2 MeV were then calculated in the TPS and compared with measured data. Compared with Bortfeld's original formula, our modified formula improved fitting of the measured depth dose curves at depths around three-quarters of the maximum range and in the beam entrance region. The parabolic function described the relationship between the parameters of the analytic expression of different energies. The predicted Bragg curves based on the parameters fitted using the seven measured curves accurately described the Bragg curves of proton pencil beams of 45 energies configured in our TPS. When we used the calculated Bragg curves as the input to TPS, the standard deviations of the measured and calculated data points along the 10 cm SOBPs created with proton pencil beams with maximum energies of 221.8 and 121.2 MeV were 1.19% and 1.18%, respectively, using curves predicted by the algorithm generated from the seven measured curves. Our method would be a valuable tool to analyze measured Bragg curves without the need for time-consuming measurements and correctly describe multiple Bragg curves using a closed analytical expression.