The first probe of a FLASH proton beam by PET

The recently observed FLASH effect related to high doses delivered with high rates has the potential to revolutionize radiation cancer therapy if promising results are confirmed and an underlying mechanism understood. Comprehensive measurements are essential to elucidate the phenomenon. We report the first-ever demonstration of measurements of successive in-spill and post-spill emissions of gammas arising from irradiations by a FLASH proton beam. A small positron emission tomography (PET) system was exposed in an ocular beam of the Proton Therapy Center at MD Anderson Cancer Center to view phantoms irradiated by 3.5 × 1010protons with a kinetic energy of 75.8 MeV delivered in 101.5 ms-long spills yielding a dose rate of 164 Gy s-1. Most in-spill events were due to prompt gammas. Reconstructed post-spill tomographic events, recorded for up to 20 min, yielded quantitative imaging and dosimetric information. These findings open a new and novel modality for imaging and monitoring of FLASH proton therapy exploiting in-spill prompt gamma imaging followed by post-spill PET imaging.

The first PET glimpse of a proton FLASH beam

We demonstrate the first ever recorded PET imaging and dosimetry of a FLASH proton beam at the Proton Center of the MD~Anderson Cancer Center.
Two scintillating LYSO crystal arrays, read out by silicon photomultipliers, were configured with a partial field of view of a~cylindrical PMMA phantom irradiated by a FLASH proton beam. The proton beam had a kinetic energy of 75.8\,MeV and an intensity of about $3.5\times10^{10}$ protons that were extracted over 101.5\,ms-long spills. The radiation environment was characterized by CZT and plastic scintillator counters. Preliminary results indicate that the PET technology used in our tests can efficiently record FLASH beam events. The instrument yielded informative and quantitative imaging and dosimetry of beam-activated isotopes in a PMMA phantoms, as supported by Monte Carlo simulations. These studies open a new PET modality that can lead to improved imaging and monitoring of FLASH proton therapy.

Effect of boron compounds on the biological effectiveness of proton therapy

PURPOSE

We assessed whether adding sodium borocaptate (BSH) or 4-borono-l-phenylalanine (BPA) to cells irradiated with proton beams influenced the biological effectiveness of those beams against prostate cancer cells to investigate if the alpha particles generated through proton-boron nuclear reactions would be sufficient to enhance the biological effectiveness of the proton beams.

METHODS

We measured clonogenic survival in DU145 cells treated with 80.4-ppm BSH or 86.9-ppm BPA, or their respective vehicles, after irradiation with 6-MV X-rays, 1.2-keV/μm (low linear energy transfer [LET]) protons, or 9.9-keV/μm (high-LET) protons. We also measured γH2AX and 53BP1 foci in treated cells at 1 and 24 h after irradiation with the same conditions.

RESULTS

We found that BSH radiosensitized DU145 cells across all radiation types. However, no difference was found in relative radiosensitization, characterized by the sensitization enhancement ratio or the relative biological effectiveness, for vehicle- versus BSH-treated cells. No differences were found in numbers of γH2AX or 53BP1 foci or γH2AX/53BP1 colocalized foci for vehicle- versus BSH-treated cells across radiation types. BPA did not radiosensitize DU145 cells nor induced any significant differences when comparing vehicle- versus BPA-treated cells for clonogenic cell survival or γH2AX and 53BP1 foci or γH2AX/53BP1 colocalized foci.

CONCLUSIONS

Treatment with 11 B, at concentrations of 80.4 ppm from BSH or 86.9 ppm from BPA, had no effect on the biological effectiveness of proton beams in DU145 prostate cancer cells. Our results agree with published theoretical calculations indicating that the contribution of alpha particles from such reactions to the total absorbed dose and biological effectiveness is negligible. We also found that BSH radiosensitized DU145 cells to X-rays, low-LET protons, and high-LET protons but that the radiosensitization was not related to DNA damage.

Adaptation and dosimetric commissioning of a synchrotron-based proton beamline for FLASH experiments

Objective. Irradiation with ultra-high dose rates (>40 Gy s−1), also known as FLASH irradiation, has the potential to shift the paradigm of radiation therapy because of its reduced toxicity to normal tissues compared to that of conventional irradiations. The goal of this study was to (1) achieve FLASH irradiation conditions suitable for pre-clinical i n vitro and in vivo biology experiments using our synchrotron-based proton beamline and (2) commission the FLASH irradiation conditions achieved. Approach. To achieve these suitable FLASH conditions, we made a series of adaptations to our proton beamline, including modifying the spill length and size of accelerating cycles, repurposing the reference monitor for dose control, and expanding the field size with a custom double-scattering system. We performed the dosimetric commissioning with measurements using an Advanced Markus chamber and EBT-XD films as well as with Monte Carlo simulations. Main results. Through adaptations, we have successfully achieved FLASH irradiation conditions, with an average dose rate of up to 375 Gy s−1. The Advanced Markus chamber was shown to be appropriate for absolute dose calibration under our FLASH conditions with a recombination factor ranging from 1.002 to 1.006 because of the continuous nature of our synchrotron-based proton delivery within a spill. Additionally, the absolute dose measured using the Advanced Markus chamber and EBT-XD films agreed well, with average and maximum differences of 0.32% and 1.63%, respectively. We also performed a comprehensive temporal analysis for FLASH spills produced by our system, which helped us identify a unique relationship between the average dose rate and the dose in our FLASH irradiation. Significance. We have established a synchrotron-based proton FLASH irradiation platform with accurate and precise dosimetry that is suitable for pre-clinical biology experiments. The unique time structure of the FLASH irradiation produced by our synchrotron-based system may shed new light onto the mechanism behind the FLASH effect.

Roadmap: helium ion therapy

Helium ion beam therapy for the treatment of cancer was one of several developed and studied particle treatments in the 1950s, leading to clinical trials beginning in 1975 at the Lawrence Berkeley National Laboratory. The trial shutdown was followed by decades of research and clinical silence on the topic while proton and carbon ion therapy made debuts at research facilities and academic hospitals worldwide. The lack of progression in understanding the principle facets of helium ion beam therapy in terms of physics, biological and clinical findings persists today, mainly attributable to its highly limited availability. Despite this major setback, there is an increasing focus on evaluating and establishing clinical and research programs using helium ion beams, with both therapy and imaging initiatives to supplement the clinical palette of radiotherapy in the treatment of aggressive disease and sensitive clinical cases. Moreover, due its intermediate physical and radio-biological properties between proton and carbon ion beams, helium ions may provide a streamlined economic steppingstone towards an era of widespread use of different particle species in light and heavy ion therapy. With respect to the clinical proton beams, helium ions exhibit superior physical properties such as reduced lateral scattering and range straggling with higher relative biological effectiveness (RBE) and dose-weighted linear energy transfer (LET_d) ranging from ∼4 keVμm^-1to ∼40 keVμm^-1. In the frame of heavy ion therapy using carbon, oxygen or neon ions, where LET_dincreases beyond 100 keVμm^-1, helium ions exhibit similar physical attributes such as a sharp lateral penumbra, however, with reduced radio-biological uncertainties and without potentially spoiling dose distributions due to excess fragmentation of heavier ion beams, particularly for higher penetration depths. This roadmap presents an overview of the current state-of-the-art and future directions of helium ion therapy: understanding physics and improving modeling, understanding biology and improving modeling, imaging techniques using helium ions and refining and establishing clinical approaches and aims from learned experience with protons. These topics are organized and presented into three main sections, outlining current and future tasks in establishing clinical and research programs using helium ion beams-A. Physics B. Biological and C. Clinical Perspectives.

Optimization of FLASH Proton Beams Using a Track-Repeating Algorithm

BACKGROUND

Radiation with high dose rate (FLASH) has shown to reduce toxicities to normal tissues around the target and maintain tumor control with the same amount of dose compared to conventional radiation. This phenomenon has been widely studied in electron therapy, which is often used for shallow tumor treatment. Proton therapy is considered as a more suitable treatment modality for deep-seated tumors. The feasibility of FLASH proton therapy has recently been demonstrated by a series of pre- and clinical trials. One of the challenges is to efficiently generate wide enough dose distributions in both lateral and longitudinal directions to cover the entire tumor volume. The goal of this paper is to introduce a set of automatic FLASH proton beam optimization algorithms developed recently.

PURPOSE

To develop a fast and efficient optimizer for the design of a passive scattering proton FLASH radiotherapy delivery at The University of Texas M. D. Anderson Proton Therapy Center, based on the Fast Dose Calculator (FDC).

METHODS

A track repeating algorithm, FDC, was validated vs. Geant4 simulations and applied to calculate dose distributions in various beamline setups. The design of the components was optimized to deliver homogeneous fields with well-defined diameters between 11.0 mm and 20.5 mm, as well as a spread-out Bragg peak (SOBP) with modulations between 8.5 and 39.0 mm. A ridge filter, a high-Z material scatterer, and a collimator with range compensator were inserted in the beam path, and their shapes and sizes were optimized to spread out the Bragg peak, widen the beam, and reduce the penumbra. The optimizer was developed and tested using two proton energies (87.0 MeV and 159.5 MeV) in a variety of beam line arrangements. Dose rates of the optimized beams were estimated by scaling their doses to those of unmodified beams.

RESULTS

The optimized 87.0 MeV beams, with a distance from the beam pipe window to the phantom surface (window-to-surface distance) of 550 mm, produced an 8.5-mm-wide spread-out Bragg peak (proximal 90% to distal 90% of the maximum dose); 14.5 mm, 12.0 mm, and 11.0 mm lateral widths at the 50%, 80%, and 90% dose location, respectively; and a 2.5 mm penumbra from 80% to 20% in the lateral profile. The 159.5 MeV beam had a SOBP of 39.0 mm and lateral widths of 20.5 mm, 15.0 mm, and 12.5 mm at 50%, 80%, and 90% dose location, respectively, when the window-to-surface distance (WSD) was 550 mm. Wider lateral widths were obtained with increased WSD. The SOBP modulations changed when the ridge filters with different characteristics were inserted. Dose rates on the beam central axis for all optimized beams (other than the 87.0 MeV beam with 2000 mm-WSD) were above that needed for the FLASH effect threshold (40 Gy/s) except at the very end of the depth dose profile scaling with a dose rate of 1400 Gy/s at the Bragg peak in the unmodified beams. The optimizer was able to instantly design the individual beamline components for each of the beam line setups, without the need of time intensive iterative simulations.

CONCLUSION

An efficient system, consisting of an optimizer and a Fast Dose Calculator have been developed and validated in a variety of beam line setups, comprising two proton energies, several WSDs and SOBPs. The set of automatic optimization algorithms produces beam shaping element designs efficiently and with excellent quality. This article is protected by copyright. All rights reserved.

Monte Carlo study of the mechanism of proton-boron fusion therapy

PURPOSE

Proton beam therapy has been found to have enhanced biological effectiveness in targets that contain the boron isotope ¹¹ B, with the alpha particles resulting from the p + ¹¹ B → 3α reaction being hypothesized as the mechanism; in this study, we aimed to elucidate the causes of the enhanced biological effectiveness of proton-boron fusion therapy by performing a detailed Monte Carlo study of the p + ¹¹ B → 3α reaction in a phantom geometry.

METHODS

We utilized the Geant4 toolkit to create Monte Carlo particle physics simulations. These simulations consisted of a proton beam with a range 30 mm, creating a Spread-Out Bragg Peak with a modulation width of 10 mm, directed into a water phantom containing a region of boron material. Energy deposition, particle energy, and particle fluence were scored along the path of the beam and grouped by particle species. The scoring was performed using a series of cylindrical volumes with a radius of 2.5 mm and depth of 0.1 mm, constructed such that the depth was parallel to the proton beam. Root was then used to perform the data analysis.

RESULTS

Our simulations showed that the dose delivered by alpha particles produced by p + ¹¹ B → 3α was several orders of magnitude lower than the dose delivered directly by protons, even when the boron uptake region was comprised entirely of natural boron or pure ¹¹ B.

CONCLUSIONS

Our findings do not support the theory that an alpha particle-based mechanism is responsible for the enhanced biological effectiveness of proton-boron fusion therapy. We conclude that any enhanced biological effect seen in experimental studies was not caused by fusion reactions between protons and ¹¹ B nuclei. However, it is necessary to reproduce the past experiments that indicated significant dose enhancement.

Design and validation of a synchrotron proton beam line for FLASH radiotherapy preclinical research experiments

PURPOSE

The main purpose of this work was to generate and validate the dosimetric accuracy of proton beams of dimensions that are appropriate for in vivo small animal and in vitro ultrahigh dose rate (FLASH) radiotherapy experiments using a synchrotron-based treatment delivery system. This study was performed to enable future investigations of the relevance of a spread-out Bragg peak (SOBP) under FLASH conditions.

METHODS

The spill characteristics of the small field fixed horizontal beam line were modified to deliver accelerated protons in times as short as 2 ms and to control the dose delivered. A Gaussian-like transverse beam profile was transformed into a square uniform one at FLASH dose rates, while avoiding low-dose regions, a crucial requirement to protect normal tissue during FLASH irradiation. Novel beam-shaping devices were designed using Monte Carlo techniques to produce up to about 6 cm³ of uniform dose in SOBPs while maximizing the dose rate. These included a scattering foil, a conical flattening filter to maximize the flux of protons into the region of interest, energy filters, range compensators, and collimators. The shapes, sizes, and positions of the components were varied to provide the required field sizes and SOBPs.

RESULTS

The designed and fabricated devices were used to produce 10-, 15-, and 20-mm diameter, circular field sizes and 10-, 15-, and 20-mm SOBP modulation widths at uniform physical dose rates of up to 375 Gy/s at the center of the SOBP and a minimum dose rate of about 255 Gy/s at the entrance, respectively, in cylindrical volumes. The flatness of lateral dose profiles at the center could be adjusted to within ±1.5% at the center of the SOBP. Assessment of systematic uncertainties, such as impact of misalignments and positioning uncertainties, was performed using simulations, and the results were used to provide appropriate adjustments to ensure high-accuracy FLASH beam delivery for both in vitro and in vivo preclinical experiments.

CONCLUSIONS

It is feasible to use synchrotron-generated proton beams of sufficient dimensions for FLASH radiobiology experiments. We expect to use the system we developed to acquire in vitro and in vivo small animal FLASH radiobiology data as a function of dose, dose rate, oxygen content, and linear energy transfer to help us understand the underlying mechanisms of the FLASH phenomenon.

Mapping the Relative Biological Effectiveness of Proton, Helium and Carbon Ions with High-Throughput Techniques

Large amounts of high quality biophysical data are needed to improve current biological effects models but such data are lacking and difficult to obtain. The present study aimed to more efficiently measure the spatial distribution of relative biological effectiveness (RBE) of charged particle beams using a novel high-accuracy and high-throughput experimental platform. Clonogenic survival was selected as the biological endpoint for two lung cancer cell lines, H460 and H1437, irradiated with protons, carbon, and helium ions. Ion-specific multi-step microplate holders were fabricated such that each column of a 96-well microplate is spatially situated at a different location along a particle beam path. Dose, dose-averaged linear energy transfer (LET_d), and dose-mean lineal energy (y_d) were calculated using an experimentally validated Geant4-based Monte Carlo system. Cells were irradiated at the Heidelberg Ion Beam Therapy Center (HIT). The experimental results showed that the clonogenic survival curves of all tested ions were y_d-dependent. Both helium and carbon ions achieved maximum RBEs within specific y_d ranges before biological efficacy declined, indicating an overkill effect. For protons, no overkill was observed, but RBE increased distal to the Bragg peak. Measured RBE profiles strongly depend on the physical characteristics such as y_d and are ion specific.

Monte Carlo Simulations of Neutron Ambient Dose Equivalent in a Novel Proton Therapy Facility Design

Purpose

The neutron shielding properties of the concrete structures of a proposed proton therapy facility were evaluated with help of the Monte Carlo technique. The planned facility's design omits the typical maze-structured entrances to the treatment rooms to facilitate more efficient access and, instead, proposes the use of massive concrete/steel doors. Furthermore, straight conduits in the treatment room walls were used in the design of the facility, necessitating a detailed investigation of the neutron radiation outside the rooms to determine if the design can be applied without violating existing radiation protection regulations. This study was performed to investigate whether the operation of a proton therapy unit using such a facility design will be in compliance with radiation protection requirements.

Methods

A detailed model of the planned proton therapy expansion project of the University of Texas, M. D. Anderson Cancer Center in Houston, Texas, was produced to simulate secondary neutron production from clinical proton beams using the MCNPX Monte Carlo radiation transport code. Neutron spectral fluences were collected at locations of interest and converted to ambient dose equivalents using an in-house code based on fluence to dose-conversion factors provided by the International Commission on Radiological Protection.

Results and Conclusions

At all investigated locations of interest, the ambient dose equivalent values were below the occupational dose limits and the dose limits for individual members of the public. The impact of straight conduits was negligible because their location and orientation were such that no line of sight to the neutron sources (ie, the isocenter locations) was established. Finally, the treatment room doors were specially designed to provide spatial efficiency and, compared with traditional maze designs, showed that while it would be possible to achieve a lower neutron ambient dose equivalent with a maze, the increased spatial (and financial) requirements may offset this advantage.

Comparing 2 Monte Carlo Systems in Use for Proton Therapy Research

Purpose

Several Monte Carlo transport codes are available for medical physics users. To ensure confidence in the accuracy of the codes, they must be continually cross-validated. This study provides comparisons between MC² and Tool for Particle Simulation (TOPAS) simulations, that is, between medical physics applications for Monte Carlo N-Particle Transport Code (MCNPX) and Geant4.

Materials and Methods

Monte Carlo simulations were repeated with 2 wrapper codes: TOPAS (based on Geant4) and MC² (based on MCNPX). Simulations increased in geometrical complexity from a monoenergetic beam incident on a water phantom, to a monoenergetic beam incident on a water phantom with a bone or tissue slab at various depths, to a spread-out Bragg peak incident on a voxelized computed tomography (CT) geometry. The CT geometry cases consisted of head and neck tissue and lung tissue. The results of the simulations were compared with one another through dose or energy deposition profiles, r ₉₀ calculations, and γ-analyses.

Results

Both codes gave very similar results with monoenergetic beams incident on a water phantom. Systematic differences were observed between MC² and TOPAS simulations when using a lung or bone slab in a water phantom, particularly in the r ₉₀ values, where TOPAS consistently calculated r ₉₀ to be deeper by about 0.4%. When comparing the performance of the 2 codes in a CT geometry, the results were still very similar, exemplified by a 3-dimensional γ-analysis pass rate > 95% at the 2%-2-mm criterion for tissues from both head and neck and lung.

Conclusion

Differences between TOPAS and MC² were minor and were not considered clinically relevant.

Nonhomologous End Joining Is More Important Than Proton Linear Energy Transfer in Dictating Cell Death

PURPOSE

This study seeks to identify biological factors that may yield a therapeutic advantage of proton therapy versus photon therapy. Specifically, we address the role of nonhomologous end-joining (NHEJ) and homologous recombination (HR) in the survival of cells in response to clinical photon and proton beams.

METHODS AND MATERIALS

We irradiated HT1080, M059K (DNA-PKcs^+/+), and HCC1937 human cancer cell lines and their isogenic counterparts HT1080-shDNA-PKcs, HT1080-shRAD51^IND, M059J (DNA-PKcs^-/-), and HCC1937-BRCA1 (BRCA1 complemented) to assess cell clonogenic survival and γ-H2AX radiation-induced foci. Cells were irradiated with either clinically relevant photons or 1 of 3 proton linear energy transfer (LET) values.

RESULTS

Our results indicate that NHEJ deficiency is more important in dictating cell survival than proton LET. Cells with disrupted HR through BRCA1 mutation showed increased radiosensitivity only for high-LET protons whereas RAD51 depletion showed increased radiosensitivity for both photons and protons. DNA double strand breaks, assessed by γ-H2AX radiation-induced foci, showed greater numbers after 24 hours in cells exposed to higher LET protons. We also observed that NHEJ-deficient cells were unable to repair the vast majority of double strand breaks after 24 hours.

CONCLUSIONS

BRCA1 mutation significantly sensitizes cells to protons, but not photons. Loss of NHEJ renders cells hypersensitive to radiation, whereas the relative importance of HR increases with LET across several cell lines. This may be attributable to the more clustered damage induced by higher LET protons, which are harder to repair through NHEJ. This highlights the importance of tumor biology in dictating treatment modality and suggests BRCA1 as a potential biomarker for proton therapy response. Our data also support the use of pharmacologic inhibitors of DNA repair to enhance the sensitivity to different radiation types, although this raises issues for normal tissue toxicity.

Fixed- versus Variable-RBE Computations for Intensity Modulated Proton Therapy

Purpose

To evaluate how using models of proton therapy that incorporate variable relative biological effectiveness (RBE) versus the current practice of using a fixed RBE of 1.1 affects dosimetric indices on treatment plans for large cohorts of patients treated with intensity modulated proton therapy (IMPT).

Methods and Materials

Treatment plans for 4 groups of patients who received IMPT for brain, head-and-neck, thoracic, or prostate cancer were selected. Dose distributions were recalculated in 4 ways: 1 with a fast-dose Monte Carlo calculator with fixed RBE and 3 with RBE calculated to 3 different models—McNamara, Wedenberg, and repair-misrepair-fixation. Differences among dosimetric indices (D02, D50, D98, and mean dose) for target volumes and organs at risk (OARs) on each plan were compared between the fixed-RBE and variable-RBE calculations.

Results

In analyses of all target volumes, for which the main concern is underprediction or RBE less than 1.1, none of the models predicted an RBE less than 1.05 for any of the cohorts. For OARs, the 2 models based on linear energy transfer, McNamara and Wedenberg, systematically predicted RBE >1.1 for most structures. For the mean dose of 25% of the plans for 2 OARs, they predict RBE equal to or larger than 1.4, 1.3, 1.3, and 1.2 for brain, head-and-neck, thorax, and prostate, respectively. Systematically lower increases in RBE are predicted by repair-misrepair-fixation, with a few cases (eg, femur) in which the RBE is less than 1.1 for all plans.

Conclusions

The variable-RBE models predict increased doses to various OARs, suggesting that strategies to reduce high-dose linear energy transfer in critical structures should be developed to minimize possible toxicity associated with IMPT.

RBE Model-Based Biological Dose Optimization for Proton Radiobiology Studies

Purpose

The purpose of the current study was (1) to develop a straightforward and rapid method to incorporate a dose-averaged linear energy transfer (LET _d )-based biological effect model into a dose optimization algorithm for scanned protons; and (2) to apply a novel beam delivery strategy with increased LET _d within the target, thereby enhancing the biological effect predicted using the selected relative biological effectiveness (RBE) model.

Materials and Methods

We first generated pristine dose Bragg curves in water and their corresponding LET _d distributions for 94 groups of proton beams, using experimentally validated Geant4 Monte Carlo simulations. Next, we developed 1-dimensional dose optimization algorithms by using the Python programming language. To calculate the RBE of protons for biological dose optimization, we invoked the McNamara RBE model and applied the radiobiological parameters of the lung cancer H460 cell line with ¹³⁷Cs reference photons.

Results

High-accuracy optimization results were obtained. The relative difference between the delivered dose and the prescribed dose was approximately within ±1.0% in the target. In addition, we obtained the RBE enhancement within the target by applying the LET-painting technique. For example, considering a simple case in which 2 opposed downslope dose fields were superimposed to form a uniform dose in the 5- to 10-cm target region, the center RBE was 1.23 ± 0.01, which was greater than the center RBE of 1.16 ± 0.01 found when using the traditional method of delivering 2 opposed flat dose fields.

Conclusion

We have successfully developed an easy-to-implement method to perform the biological dose optimization procedure by invoking the McNamara RBE model in the iteration process using the Python programming language. According to the RBE model predictions, we conclude that the increased target LET _d enhances the RBE. The accuracy of the RBE model predictions needs to be validated in radiobiological experiments.

Comparison of Monte Carlo and analytical dose computations for intensity modulated proton therapy

To evaluate the effect of approximations in clinical analytical calculations performed by a treatment planning system (TPS) on dosimetric indices in intensity modulated proton therapy. TPS calculated dose distributions were compared with dose distributions as estimated by Monte Carlo (MC) simulations, calculated with the fast dose calculator (FDC) a system previously benchmarked to full MC. This study analyzed a total of 525 patients for four treatment sites (brain, head-and-neck, thorax and prostate). Dosimetric indices (D02, D05, D20, D50, D95, D98, EUD and Mean Dose) and a gamma-index analysis were utilized to evaluate the differences. The gamma-index passing rates for a 3%/3 mm criterion for voxels with a dose larger than 10% of the maximum dose had a median larger than 98% for all sites. The median difference for all dosimetric indices for target volumes was less than 2% for all cases. However, differences for target volumes as large as 10% were found for 2% of the thoracic patients. For organs at risk (OARs), the median absolute dose difference was smaller than 2 Gy for all indices and cohorts. However, absolute dose differences as large as 10 Gy were found for some small volume organs in brain and head-and-neck patients. This analysis concludes that for a fraction of the patients studied, TPS may overestimate the dose in the target by as much as 10%, while for some OARs the dose could be underestimated by as much as 10 Gy. Monte Carlo dose calculations may be needed to ensure more accurate dose computations to improve target coverage and sparing of OARs in proton therapy.

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.

Differences in Normal Tissue Response in the Esophagus Between Proton and Photon Radiation Therapy for Non-Small Cell Lung Cancer Using In Vivo Imaging Biomarkers

PURPOSE

To determine whether there exists any significant difference in normal tissue toxicity between intensity modulated radiation therapy (IMRT) or proton therapy for the treatment of non-small cell lung cancer.

METHODS AND MATERIALS

A total of 134 study patients (n=49 treated with proton therapy, n=85 with IMRT) treated in a randomized trial had a previously validated esophageal toxicity imaging biomarker, esophageal expansion, quantified during radiation therapy, as well as esophagitis grade (Common Terminology Criteria for Adverse Events version 3.0), on a weekly basis during treatment. Differences between the 2 modalities were statically analyzed using the imaging biomarker metric value (Kruskal-Wallis analysis of variance), as well as the incidence and severity of esophagitis grade (χ² and Fisher exact tests, respectively). The dose-response of the imaging biomarker was also compared between modalities using esophageal equivalent uniform dose, as well as delivered dose to an isotropic esophageal subvolume.

RESULTS

No statistically significant difference in the distribution of esophagitis grade, the incidence of grade ≥3 esophagitis (15 and 11 patients treated with IMRT and proton therapy, respectively), or the esophageal expansion imaging biomarker between cohorts (P>.05) was found. The distribution of imaging biomarker metric values had similar distributions between treatment arms, despite a slightly higher dose volume in the proton arm (P>.05). Imaging biomarker dose-response was similar between modalities for dose quantified as esophageal equivalent uniform dose and delivered esophageal subvolume dose. Regardless of treatment modality, there was high variability in imaging biomarker response, as well as esophagitis grade, for similar esophageal doses between patients.

CONCLUSIONS

There was no significant difference in esophageal toxicity from either proton- or photon-based radiation therapy as quantified by esophagitis grade or the esophageal expansion imaging biomarker.

Optimization of Monte Carlo particle transport parameters and validation of a novel high throughput experimental setup to measure the biological effects of particle beams

PURPOSE

Accurate modeling of the relative biological effectiveness (RBE) of particle beams requires increased systematic in vitro studies with human cell lines with care towards minimizing uncertainties in biologic assays as well as physical parameters. In this study, we describe a novel high-throughput experimental setup and an optimized parameterization of the Monte Carlo (MC) simulation technique that is universally applicable for accurate determination of RBE of clinical ion beams. Clonogenic cell-survival measurements on a human lung cancer cell line (H460) are presented using proton irradiation.

METHODS

Experiments were performed at the Heidelberg Ion Therapy Center (HIT) with support from the Deutsches Krebsforschungszentrum (DKFZ) in Heidelberg, Germany using a mono-energetic horizontal proton beam. A custom-made variable range selector was designed for the horizontal beam line using the Geant4 MC toolkit. This unique setup enabled a high-throughput clonogenic assay investigation of multiple, well defined dose and linear energy transfer (LETs) per irradiation for human lung cancer cells (H460) cultured in a 96-well plate. Sensitivity studies based on application of different physics lists in conjunction with different electromagnetic constructors and production threshold values to the MC simulations were undertaken for accurate assessment of the calculated dose and the dose-averaged LET (LET_d ). These studies were extended to helium and carbon ion beams.

RESULTS

Sensitivity analysis of the MC parameterization revealed substantial dependence of the dose and LET_d values on both the choice of physics list and the production threshold values. While the dose and LET_d calculations using FTFP_BERT_LIV, FTFP_BERT_EMZ, FTFP_BERT_PEN and QGSP_BIC_EMY physics lists agree well with each other for all three ions, they show large differences when compared to the FTFP_BERT physics list with the default electromagnetic constructor. For carbon ions, the dose corresponding to the largest LET_d value is observed to differ by as much as 78% between FTFP_BERT and FTFP_BERT_LIV. Furthermore, between the production threshold of 700 μm and 5 μm, proton dose varies by as much as 19% corresponding to the largest LET_d value sampled in the current investigation. Based on the sensitivity studies, the FTFP_BERT physics list with the low energy Livermore electromagnetic constructor and a production threshold of 5 μm was employed for determining accurate dose and LET_d . The optimized MC parameterization results in a different LET_d dependence of the RBE curve for 10% SF of the H460 cell line irradiated with proton beam when compared with the results from a previous study using the same cell line. When the MC parameters are kept consistent between the studies, the proton RBE results agree well with each other within the experimental uncertainties.

CONCLUSIONS

A custom high-throughput, high-accuracy experimental design for accurate in vitro cell survival measurements was employed at a horizontal beam line. High sensitivity of the physics-based optimization establishes the importance of accurate MC parameterization and hence the conditioning of the MC system on a case-by-case basis. The proton RBE results from current investigations are observed to agree with a previous measurement made under different experimental conditions. This establishes the consistency of our experimental findings across different experiments and institutions.

Analysis of the track- and dose-averaged LET and LET spectra in proton therapy using the geant4 Monte Carlo code

PURPOSE

The motivation of this study was to find and eliminate the cause of errors in dose-averaged linear energy transfer (LET) calculations from therapeutic protons in small targets, such as biological cell layers, calculated using the geant 4 Monte Carlo code. Furthermore, the purpose was also to provide a recommendation to select an appropriate LET quantity from geant 4 simulations to correlate with biological effectiveness of therapeutic protons.

METHODS

The authors developed a particle tracking step based strategy to calculate the average LET quantities (track-averaged LET, LETt and dose-averaged LET, LETd) using geant 4 for different tracking step size limits. A step size limit refers to the maximally allowable tracking step length. The authors investigated how the tracking step size limit influenced the calculated LETt and LETd of protons with six different step limits ranging from 1 to 500 μm in a water phantom irradiated by a 79.7-MeV clinical proton beam. In addition, the authors analyzed the detailed stochastic energy deposition information including fluence spectra and dose spectra of the energy-deposition-per-step of protons. As a reference, the authors also calculated the averaged LET and analyzed the LET spectra combining the Monte Carlo method and the deterministic method. Relative biological effectiveness (RBE) calculations were performed to illustrate the impact of different LET calculation methods on the RBE-weighted dose.

RESULTS

Simulation results showed that the step limit effect was small for LETt but significant for LETd. This resulted from differences in the energy-deposition-per-step between the fluence spectra and dose spectra at different depths in the phantom. Using the Monte Carlo particle tracking method in geant 4 can result in incorrect LETd calculation results in the dose plateau region for small step limits. The erroneous LETd results can be attributed to the algorithm to determine fluctuations in energy deposition along the tracking step in geant 4. The incorrect LETd values lead to substantial differences in the calculated RBE.

CONCLUSIONS

When the geant 4 particle tracking method is used to calculate the average LET values within targets with a small step limit, such as smaller than 500 μm, the authors recommend the use of LETt in the dose plateau region and LETd around the Bragg peak. For a large step limit, i.e., 500 μm, LETd is recommended along the whole Bragg curve. The transition point depends on beam parameters and can be found by determining the location where the gradient of the ratio of LETd and LETt becomes positive.

Degradation of proton depth dose distributions attributable to microstructures in lung-equivalent material

PURPOSE

The purpose of the work reported here was to investigate the influence of sub-millimeter size heterogeneities on the degradation of the distal edges of proton beams and to validate Monte Carlo (MC) methods' ability to correctly predict such degradation.

METHODS

A custom-designed high-resolution plastic phantom approximating highly heterogeneous, lung-like structures was employed in measurements and in Monte Carlo simulations to evaluate the degradation of proton Bragg curves penetrating heterogeneous media.

RESULTS

Significant differences in distal falloff widths and in peak dose values were observed in the measured and the Monte Carlo simulated curves compared to pristine proton Bragg curves. Furthermore, differences between simulations of beams penetrating CT images of the phantom did not agree well with the corresponding experimental differences. The distal falloff widths in CT image-based geometries were underestimated by up to 0.2 cm in water (corresponding to 0.8-1.4 cm in lung tissue), and the peak dose values of pristine proton beams were overestimated by as much as ˜35% compared to measured curves or depth-dose curves simulated on the basis of true geometry. The authors demonstrate that these discrepancies were caused by the limited spatial resolution of CT images that served as a basis for dose calculations and lead to underestimation of the impact of the fine structure of tissue heterogeneities. A convolution model was successfully applied to mitigate the underestimation.

CONCLUSIONS

The results of this study justify further development of models to better represent heterogeneity effects in soft-tissue geometries, such as lung, and to correct systematic underestimation of the degradation of the distal edge of proton doses.