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 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 MD Anderson Proton Therapy Center, based on the fast dose calculator (FDC). METHODS: A track-repeating algorithm, FDC, was validated versus 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 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 and 159.5 MeV) in a variety of beamline 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 [WSD]) of 550 mm, produced an 8.5-mm-wide SOBP (proximal 90% to distal 90% of the maximum dose); 14.5, 12.0, 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 an SOBP of 39.0 mm and lateral widths of 20.5, 15.0, and 12.5 mm at 50%, 80%, and 90% dose location, respectively, when the 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 beamline setups, without the need of time intensive iterative simulations. CONCLUSION: An efficient system, consisting of an optimizer and an FDC have been developed and validated in a variety of beamline 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.

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-clinicalin vitroandin vivobiology 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.

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.

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 (LETd) ranging from ∼4 keVµm-1to ∼40 keVµm-1. In the frame of heavy ion therapy using carbon, oxygen or neon ions, where LETdincreases 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.

Technical note: 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 11 B, with the alpha particles resulting from the p + 11 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 + 11 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 + 11 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 11 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 11 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 cm3 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 (LETd), and dose-mean lineal energy (yd) 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 yd-dependent. Both helium and carbon ions achieved maximum RBEs within specific yd 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 yd 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 MC2 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 MC2 (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 90 calculations, and γ-analyses. RESULTS: Both codes gave very similar results with monoenergetic beams incident on a water phantom. Systematic differences were observed between MC2 and TOPAS simulations when using a lung or bone slab in a water phantom, particularly in the r 90 values, where TOPAS consistently calculated r 90 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 MC2 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-shRAD51IND, 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 137Cs 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.

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.

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 (χ2 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 (LETd ). These studies were extended to helium and carbon ion beams. RESULTS: Sensitivity analysis of the MC parameterization revealed substantial dependence of the dose and LETd values on both the choice of physics list and the production threshold values. While the dose and LETd 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 LETd 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 LETd 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 LETd . The optimized MC parameterization results in a different LETd 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.

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.

Clinical evidence of variable proton biological effectiveness in pediatric patients treated for ependymoma.

BACKGROUND AND PURPOSE: A constant relative biological effectiveness (RBE) is used for clinical proton therapy; however, experimental evidence indicates that RBE can vary. We analyzed pediatric ependymoma patients who received proton therapy to determine if areas of normal tissue damage indicated by post-treatment image changes were associated with increased biological dose effectiveness. MATERIAL AND METHODS: Fourteen of 34 children showed T2-FLAIR hyperintensity on post-treatment magnetic resonance (MR) images. We delineated regions of treatment-related change and calculated dose and linear energy transfer (LET) distributions with Monte Carlo. Voxel-level image change data were fit to a generalized linear model incorporating dose and LET. Cross-validation was used to determine model parameters and for receiver operating characteristic curve analysis. Tolerance dose (TD50; dose at which 50% of patients would experience toxicity) was interpolated from the model. RESULTS: Image changes showed dependence on increasing LET and dose. TD50 decreased with increasing LET, indicating an increase in biological dose effectiveness. The cross-validated area under the curve for the model was 0.91 (95% confidence interval 0.88-0.94). CONCLUSIONS: Our correlation of changes on MR images after proton therapy with increased LET constitutes the first clinical evidence of variable proton biological effectiveness.