Assessment of the accuracy of an MCNPX-based Monte Carlo simulation model for predicting three-dimensional absorbed dose distributions

Targeted Inhibition of DNA-PKcs, ATM, ATR, PARP, and Rad51 Modulate Response to X Rays and Protons

Small molecule inhibitors are currently in preclinical and clinical development for the treatment of selected cancers, particularly those with existing genetic alterations in DNA repair and DNA damage response (DDR) pathways. Keen interest has also been expressed in combining such agents with other targeted antitumor strategies such as radiotherapy. Radiotherapy exerts its cytotoxic effects primarily through DNA damage-induced cell death; therefore, inhibiting DNA repair and the DDR should lead to additive and/or synergistic radiosensitizing effects. In this study we screened the response to X-ray or proton radiation in cell lines treated with DDR inhibitors (DDRis) targeting ATM, ATR, DNA-PKcs, Rad51, and PARP, with survival metrics established using clonogenic assays. We observed that DDRis generate significant radiosensitization in cancer and primary cells derived from normal tissue. Existing genetic defects in cancer cells appear to be an important consideration when determining the optimal inhibitor to use for synergistic combination with radiation. We also show that while greater radiosensitization can be achieved with protons (9.9 keV/µm) combined with DDRis, the relative biological effectiveness is unchanged or in some cases reduced. Our results indicate that while targeting the DDR can significantly radiosensitize cancer cells to such combinations, normal cells may also be equally or more severely affected, depending on the DDRi used. These data highlight the importance of identifying genetic defects as predictive biomarkers of response for combination treatment.

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.

Proton Therapy for Prostate Cancer: Challenges and Opportunities

Simple Summary

Reported clinical outcomes of proton therapy (PT) for localized prostate cancer are similar to photon-based external beam radiotherapy. Apparently, the dosimetric advantages of PT have yet to be translated to clinical benefits. The suboptimal clinical outcomes of PT might be attributable to inadequate dose prescription, as indicated by the ASCENDE-RT trial. Moreover, uncertainties involved in the treatment planning and delivery processes, as well as technological limitations in PT treatment systems, may lead to discrepancies between planned doses and actual doses delivered to patients. In this article, we reviewed the current status of PT for prostate cancer and discussed different clinical implementations that could potentially improve the clinical outcome of PT for prostate cancer. Various technological advancements under which uncertainties in dose calculations can be minimized, including MRI-guided PT, dual-energy photon-counting CT and high-resolution Monte Carlo-based treatment planning systems, are highlighted.

Abstract

The dosimetric advantages of proton therapy (PT) treatment plans are demonstrably superior to photon-based external beam radiotherapy (EBRT) for localized prostate cancer, but the reported clinical outcomes are similar. This may be due to inadequate dose prescription, especially in high-risk disease, as indicated by the ASCENDE-RT trial. Alternatively, the lack of clinical benefits with PT may be attributable to improper dose delivery, mainly due to geometric and dosimetric uncertainties during treatment planning, as well as delivery procedures that compromise the dose conformity of treatments. Advanced high-precision PT technologies, and treatment planning and beam delivery techniques are being developed to address these uncertainties. For instance, external magnetic resonance imaging (MRI)-guided patient setup rooms are being developed to improve the accuracy of patient positioning for treatment. In-room MRI-guided patient positioning systems are also being investigated to improve the geometric accuracy of PT. Soon, high-dose rate beam delivery systems will shorten beam delivery time to within one breath hold, minimizing the effects of organ motion and patient movements. Dual-energy photon-counting computed tomography and high-resolution Monte Carlo-based treatment planning systems are available to minimize uncertainties in dose planning calculations. Advanced in-room treatment verification tools such as prompt gamma detector systems will be used to verify the depth of PT. Clinical implementation of these new technologies is expected to improve the accuracy and dose conformity of PT in the treatment of localized prostate cancers, and lead to better clinical outcomes. Improvement in dose conformity may also facilitate dose escalation, improving local control and implementation of hypofractionation treatment schemes to improve patient throughput and make PT more cost effective.

Dose characteristics of Au-198 eye brachytherapy applicator: A Monte Carlo study

PURPOSE

The use of ocular plaques is a promising treatment option for eye melanoma brachytherapy. Although several studies have been done on various ocular plaques, little is known about the dose characterization of ¹⁹⁸Au plaque.

MATERIALS AND METHOD

The full mathematical model of the eye phantom, tumor, ¹⁰⁶Ru/¹⁰⁶Rh CCA, and ¹⁹⁸Au plaque were simulated using the Monte Carlo MCNPX code. The dose distribution was measured in the plaque's central axis direction, and a dose profile was also measured at a distance of 2.5 mm from the plaque surface.

RESULTS

The findings showed that ¹⁹⁸Au plaque has superior dosimetric characteristics than CCA plaque for tumors with a thickness of greater than 3.5 mm, while CCA plaque is better for tumors with a thickness of less than 3.5 mm. The dose to the sclera and choroid is higher in the case of CCA plaque, while the dose to the organs at risk (lens and optic nerve) is greater in the case of ¹⁹⁸Au applicator. In the case of ¹⁹⁸Au plaque, however, the dose to sensitive organs was within their permissible dose range.

CONCLUSION

In the treatment of medium and large tumors, ¹⁹⁸Au plaque is more successful than CCA plaque. It can produce a much more homogeneous lateral dose profile in the target. In the treatment of dome-shaped tumors, ¹⁹⁸Au plaque may be more successful than CCA plaque. As a result, the tumor's shape influences the plaque type selection.

Outcomes of and treatment planning considerations for a hybrid technique delivering proton pencil-beam scanning radiation to women with metal-containing tissue expanders undergoing post-mastectomy radiation

BACKGROUND

Following mastectomy, immediate breast reconstruction often involves the use of temporary tissue expanders (TEs). TEs contain metallic ports (MPs), which complicate proton pencil-beam scanning (PBS) planning. A technique was implemented for delivering PBS post-mastectomy radiation (PMRT) to patients with TEs and MPs.

METHODS

A protocol utilizing a hybrid single- and multi-field optimization (SFO, MFO) technique was developed. Plans were robustly optimized using a Monte Carlo algorithm. A CTV_eval structure including chest wall (CW) and regional nodal (RNI) targets and excluding the TE was evaluated. Organ at risk (OAR) dosimetry and acute toxicities were analyzed.

RESULTS

Twenty-nine women were treated with this technique. A 2-field SFO technique was used superior and inferior to the MP, with a 3 or 4-field MFO technique used at the level of the MP. Virtual blocks were utilized so that beams did not travel through the MP. A port-to-CW distance of 1 cm was required. Patients underwent daily image-guidance to ensure the port remained within a 0.5 cm internal planning volume (ITV). Median RT dose to CTV_eval was 50.4 Gy (45.0-50.4). Median 95% CTV_eval coverage was 99.5% (95-100). Optically stimulated luminescent dosimeter (OSLD) readings were available for 8 patients and correlated to the dose measurements in the treatment planning system (TPS); median OSLD ratio was 0.99 (range, 0.93-1.02).

CONCLUSIONS

Delivering PMRT with PBS for women with metal-containing TEs using a hybrid SFO/MFO technique is feasible, reproducible, and achieves excellent dose distributions. Specialized planning and image-guidance techniques are required to safely utilize this treatment in the clinic.

A framework for voxel-based assessment of biological effect after proton radiotherapy in pediatric brain cancer patients using multi-modal imaging

INTRODUCTION

The exact dependence of biological effect on dose and linear energy transfer (LET) in human tissue when delivering proton therapy is unknown. In this study, we propose a framework for measuring this dependency using multi-modal image-based assays with deformable registrations within imaging sessions and across time.

MATERIALS AND METHODS

3T MRI scans were prospectively collected from 6 pediatric brain cancer patients before they underwent proton therapy treatment, and every 3 months for a year after treatment. Scans included T1-weighted with contrast enhancement (T1), T2-FLAIR (T2) and fractional anisotropy (FA) images. In addition, the planning CT, dose distributions and Monte Carlo-calculated LET distributions were collected. A multi-modal deformable image registration framework was used to create a dataset of dose, LET and imaging intensities at baseline and follow-up on a voxel-by-voxel basis. We modelled the biological effect of dose and LET from proton therapy using imaging changes over time as a surrogate for biological effect. We investigated various models to show the feasibility of the framework to model imaging changes. To account for interpatient and intrapatient variations, we used a nested generalized linear mixed regression model. The models were applied to predict imaging changes over time as a function of dose and LET for each modality.

RESULTS

Using the nested models to predict imaging changes, we saw a decrease in the FA signal as a function of dose; however, the signal increased with increasing LET. Similarly, we saw an increase in T2 signal as a function of dose, but a decrease in signal with LET. We saw no changes in T1 voxel values as a function of either dose or LET.

CONCLUSIONS

The imaging changes could successfully model biological effect as a function of dose and LET using our proposed framework. Due to the low number of patients, the imaging changes observed for FA and T2 scans were not marked enough to draw any firm conclusions.

A DNA damage multiscale model for NTCP in proton and hadron therapy

PURPOSE

To develop a first principle and multiscale model for normal tissue complication probability (NTCP) as a function of dose and LET for proton and in general for particle therapy with a goal of incorporating nanoscale radio-chemical to macroscale cell biological pathways, spanning from initial DNA damage to tissue late effects.

METHODS

The method is a combination of analytical and multiscale computational steps including (a) derivation of functional dependencies of NTCP on DNA-driven cell lethality in nanometer and mapping to dose and LET in millimeter, and (b) three-dimensional-surface fitting to Monte Carlo data set generated based on postradiation image change and gathered for a cohort of 14 pediatric patients treated by scanning beam of protons for ependymoma. We categorize voxel-based dose and LET associated with development of necrosis in NTCP.

RESULT

Our model fits well the clinical data, generated for postradiation tissue toxicity and necrosis. The fitting procedure results in extraction of in vivo radio-biological α-β indices and their numerical values.

DISCUSSION AND CONCLUSION

The NTCP model, explored in this work, allows to correlate the tissue toxicities to DNA initial damage, cell lethality and the properties and qualities of radiation, dose, and LET.

Automatic phase space generation for Monte Carlo calculations of intensity modulated particle therapy

Monte Carlo (MC) is generally considered as the most accurate dose calculation tool for particle therapy. However, a proper description of the beam particle kinematics is a necessary input for a realistic simulation. Such a description can be stored in phase space (PS) files for different beam energies. A PS file contains kinetic information such as energies, positions and travelling directions for particles traversing a plane perpendicular to the beam direction. The accuracy of PS files plays a critical role in the performance of the MC method for dose calculations. A PS file can be generated with a set of parameters describing analytically the beam kinematics. However, determining such parameters can be tedious and time consuming. Thus, we have developed an algorithm to obtain those parameters automatically and efficiently. In this paper, we presented such an algorithm and compared dose calculations using PS automatically generated for the Shanghai Proton and Heavy Ion Center (SPHIC) with measurements. The gamma-index for comparing calculated depth dose distributions (DDD) with measurements are above 96.0% with criterion 0.6%/0.6 mm. For each single energy, the mean difference percentage between calculated lateral spot sizes at 5 different locations along beam direction and measurements are below 3.5%.

Sensitivity analysis of Monte Carlo model of a gantry-mounted passively scattered proton system

Purpose

This study aimed to present guidance on the correlation between treatment nozzle and proton source parameters, and dose distribution of a passive double scattering compact proton therapy unit, known as Mevion S250.

Methods

All 24 beam options were modeled using the MCNPX MC code. The calculated physical dose for pristine peak, profiles, and spread out Bragg peak (SOBP) were benchmarked with the measured data. Track‐averaged LET (LET_t) and dose‐averaged LET (LET_d) distributions were also calculated. For the sensitivity investigations, proton beam line parameters including Average Energy (AE), Energy Spread (ES), Spot Size (SS), Beam Angle (BA), Beam Offset (OA), and Second scatter Offset (SO) from central Axis, and also First Scatter (FS) thickness were simulated in different stages to obtain the uncertainty of the derived results on the physical dose and LET distribution in a water phantom.

Results

For the physical dose distribution, the MCNPX MC model matched measurements data for all the options to within 2 mm and 2% criterion. The Mevion S250 was found to have a LET_t between 0.46 and 8.76 keV.μm^–1 and a corresponding LET_d between 0.84 and 15.91 keV.μm^–1. For all the options, the AE and ES had the greatest effect on the resulting depth of pristine peak and peak‐to‐plateau ratio respectively. BA, OA, and SO significantly decreased the flatness and symmetry of the profiles. The LETs were found to be sensitive to the AE, ES, and SS, especially in the peak region.

Conclusions

This study revealed the importance of considering detailed beam parameters, and identifying those that resulted in large effects on the physical dose distribution and LETs for a compact proton therapy machine.

Systematic microdosimetric data for protons of therapeutic energies calculated with Geant4-DNA

The purpose of this study was to generate physical data needed for microdosimetry-based models of proton RBE. Our focus was on the frequency and dose average lineal energies, y  F and y  D . We report data for proton energies from 0.1 to 100 MeV, for spherical volumes 2-103 nm in diameter. These data were calculated using Geant4-DNA Monte Carlo software. The physics implemented in Geant4-DNA has been extensively tested for this type of calculations but data on y  F and y  D for protons generated with this code have been very limited. An innovative aspect of our study is that we introduced a straightforward procedure for calculation of y  F and y  D for polyenergetic beams and presented the data in a format that simplifies these calculations. We compared our data with previous studies that used different Monte Carlo codes and with experimental data.

Ambient neutron and photon dose equivalent H*(10) around a pencil beam scanning proton therapy facility

OBJECTIVES

To measure leakage ambient dose equivalent H*(10) from stray secondary neutron and photon radiation around proton therapy (PT) facility and evaluate adequacy of shielding design.

METHODS AND MATERIALS

H*(10) measurement were carried out at 149 locations around cyclotron vault (CV), beam transport system (BTS) and first treatment room (GTR3) of a multiroom PT facility using WENDI-II and SmartIon survey meter. Measurement were performed under extreme case scenarios wherein maximum secondary neutrons and photons were produced around CV, BTS and GTR3 by stopping 230MeV proton of 300nA on beam degrader, end of BTS and isocenter of GTR3. Weekly time average dose rate (TADR) were calculated from H*(10) value measured at selective hot spots by irradiating actual treatment plans of mix clinical sites.

RESULTS

The maximum total H*(10) were within 2 µSv/hr around CV, 5 µSv/hr around outer wall of BTS which increases up to 62 µSv/hr at the end of inside BTS corridor. Maximum H*(10) of 20.8 µSv/hr in treatment control console (P125), 23.4 µSv/hr behind the common wall between GTR3 and GTR2 (P132) and 25.7 µSv/hr above isocenter (P99) were observed around GTR3. Reduction of beam current from 6 to 3 nA and 1 nA at nozzle exit lead to decrease in total H*(10) at P125 from 20.8 to 11.35 and 4.62 µSv/hr. In comparison to extreme case scenario, H*(10) value at P125, P132 and P99 from clinically relevant irradiation parameters were reduce by a factor ranging from 8.6 for high range cube to 46.4 for brain clinical plan. The maximum weekly TADR per fraction was highest for large volume, sacral chordoma patient at 8.5 µSv/hr compare to 0.3 µSv/hr for brain patient. The calculated weekly TADR for 30 mix clinical cases and 15 fractions of 1 L cube resulted total weekly TADR of 83-84 µSv/hr at P125, P132 and P99. The maximum annual dose level at these hot spots were estimated at 4.37 mSv/Yr.

CONCLUSION

We have carried out an extensive measurement of H*(10) under different conditions. The shielding thickness of our PT facility is adequate to limit the dose to occupational worker and general public within the permissible stipulated limit. The data reported here can bridge the knowledge gap in ambient dose around PT facility and can also be used as a reference for any new and existing proton facility for intercomparison and validation.

ADVANCES IN KNOWLEDGE

First extensive investigation of neutron and photon H*(10) around PT facility and can bridge the knowledge gap on ambient dose.

Proton beam therapy for cancer in the era of precision medicine

Precision radiotherapy, which accurately delivers the dose on a tumor and confers little or no irradiation to the surrounding normal tissue and organs, results in maximum tumor control and decreases the toxicity to the utmost extent. Proton beam therapy (PBT) provides superior dose distributions and has a dosimetric advantage over photon beam therapy. Initially, the clinical practice and study of proton beam therapy focused on ocular tumor, skull base, paraspinal tumors (chondrosarcoma and chordoma), and unresectable sarcomas, which responded poorly when treated with photon radiotherapy. Then, it is widely regarded as an ideal mode for reirradiation and pediatrics due to reducing unwanted side effects by lessening the dose to normal tissue. During the past decade, the application of PBT has been rapidly increasing worldwide and gradually expanding for the treatment of various malignancies. However, to date, the role of PBT in clinical settings is still controversial, and there are considerable challenges in its application. We systematically review the latest advances of PBT and the challenges for patient treatment in the era of precision medicine.

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.

Evaluating the accuracy of a three-term pencil beam algorithm in heterogeneous media

The goal of this work was to evaluate the accuracy of our in-house analytical dose calculation code against MCNPX data in heterogeneous phantoms. The analytical model utilizes a pencil beam model based on Fermi-Eyges theory to account for multiple Coulomb scattering and a least-squares fit to Monte Carlo data to account for nonelastic nuclear interactions as well as any remaining, uncharacterized scatter (the 'nuclear halo'). The model characterized dose accurately (up to 1% of maximum dose in broad fields (4  ×  4 cm² and 10  ×  10 cm²) and up to 0.01% in a narrow field (0.1  ×  0.1 cm²) fit to MCNPX data). The accuracy of the model was benchmarked in three types of stylized phantoms: (1) homogeneous, (2) laterally infinite slab heterogeneities, and (3) laterally finite slab heterogeneities. Results from homogeneous phantoms and laterally infinite slab heterogeneities showed high levels of accuracy (>98% of points within 2% or 0.1 cm distance-to-agreement (DTA)). However, because range straggling and secondary particle production were not included in our model, central-axis dose differences of 2-4% were observed in laterally infinite slab heterogeneities when compared to Monte Carlo dose. In the presence of laterally finite slab heterogeneities, the analytical model resulted in lower pass rates (>96% of points within 2% or 0.1 cm DTA), which was attributed to the use of the central-axis approximation.

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.

The physics of proton therapy

The physics of proton therapy has advanced considerably since it was proposed in 1946. Today analytical equations and numerical simulation methods are available to predict and characterize many aspects of proton therapy. This article reviews the basic aspects of the physics of proton therapy, including proton interaction mechanisms, proton transport calculations, the determination of dose from therapeutic and stray radiations, and shielding design. The article discusses underlying processes as well as selected practical experimental and theoretical methods. We conclude by briefly speculating on possible future areas of research of relevance to the physics of proton therapy.

Implementation of an analytical model for leakage neutron equivalent dose in a proton radiotherapy planning system

Equivalent dose from neutrons produced during proton radiotherapy increases the predicted risk of radiogenic late effects. However, out-of-field neutron dose is not taken into account by commercial proton radiotherapy treatment planning systems. The purpose of this study was to demonstrate the feasibility of implementing an analytical model to calculate leakage neutron equivalent dose in a treatment planning system. Passive scattering proton treatment plans were created for a water phantom and for a patient. For both the phantom and patient, the neutron equivalent doses were small but non-negligible and extended far beyond the therapeutic field. The time required for neutron equivalent dose calculation was 1.6 times longer than that required for proton dose calculation, with a total calculation time of less than 1 h on one processor for both treatment plans. Our results demonstrate that it is feasible to predict neutron equivalent dose distributions using an analytical dose algorithm for individual patients with irregular surfaces and internal tissue heterogeneities. Eventually, personalized estimates of neutron equivalent dose to organs far from the treatment field may guide clinicians to create treatment plans that reduce the risk of late effects.

Feasibility of proton-activated implantable markers for proton range verification using PET

Proton beam range verification using positron emission tomography (PET) currently relies on proton activation of tissue, the products of which decay with a short half-life and necessitate an on-site PET scanner. Tissue activation is, however, negligible near the distal dose fall-off region of the proton beam range due to their high interaction energy thresholds. Therefore Monte Carlo simulation is often supplemented for comparison with measurement; however, this also may be associated with systematic and statistical uncertainties. Therefore, we sought to test the feasibility of using long-lived proton-activated external materials that are inserted or infused into the target volume for more accurate proton beam range verification that could be performed at an off-site PET scanner. We irradiated samples of ≥98% (18)O-enriched water, natural Cu foils, and >97% (68)Zn-enriched foils as candidate materials, along with samples of tissue-equivalent materials including (16)O water, heptane (C7H16), and polycarbonate (C16H14O3)n, at four depths (ranging from 100% to 3% of center of modulation (COM) dose) along the distal fall-off of a modulated 160 MeV proton beam. Samples were irradiated either directly or after being embedded in Plastic Water® or balsa wood. We then measured the activity of the samples using PET imaging for 20 or 30 min after various delay times. Measured activities of candidate materials were up to 100 times greater than those of the tissue-equivalent materials at the four distal dose fall-off depths. The differences between candidate materials and tissue-equivalent materials became more apparent after longer delays between irradiation and PET imaging, due to the longer half-lives of the candidate materials. Furthermore, the activation of the candidate materials closely mimicked the distal dose fall-off with offsets of 1 to 2 mm. Also, signals from the foils were clearly visible compared to the background from the activated Plastic Water® and balsa wood phantoms. These results indicate that markers made from these candidate materials could be used for in vivo proton range verification using an off-site PET scanner.

Benchmark measurements and simulations of dose perturbations due to metallic spheres in proton beams

Monte Carlo simulations are increasingly used for dose calculations in proton therapy due to its inherent accuracy. However, dosimetric deviations have been found using Monte Carlo code when high density materials are present in the proton beam line. The purpose of this work was to quantify the magnitude of dose perturbation caused by metal objects. We did this by comparing measurements and Monte Carlo predictions of dose perturbations caused by the presence of small metal spheres in several clinical proton therapy beams as functions of proton beam range, spread-out Bragg peak width and drift space. Monte Carlo codes MCNPX, GEANT4 and Fast Dose Calculator (FDC) were used. Generally good agreement was found between measurements and Monte Carlo predictions, with the average difference within 5% and maximum difference within 17%. The modification of multiple Coulomb scattering model in MCNPX code yielded improvement in accuracy and provided the best overall agreement with measurements. Our results confirmed that Monte Carlo codes are well suited for predicting multiple Coulomb scattering in proton therapy beams when short drift spaces are involved.

The predicted relative risk of premature ovarian failure for three radiotherapy modalities in a girl receiving craniospinal irradiation

In girls and young women, irradiation of the ovaries can reduce the number of viable ovarian primordial follicles, which may lead to premature ovarian failure (POF) and subsequently to sterility. One strategy to minimize this late effect is to reduce the radiation dose to the ovaries. A primary means of reducing dose is to choose a radiotherapy technique that avoids irradiating nearby normal tissue; however, the relative risk of POF (RRPOF) due to the various therapeutic options has not been assessed. This study compared the predicted RRPOF after craniospinal proton radiotherapy, conventional photon radiotherapy (CRT) and intensity-modulated photon radiotherapy (IMRT). We calculated the equivalent dose delivered to the ovaries of an 11-year-old girl from therapeutic and stray radiation. We then predicted the percentage of ovarian primordial follicles killed by radiation and used this as a measure of the RRPOF; we also calculated the ratio of the relative risk of POF (RRRPOF) among the three radiotherapies. Proton radiotherapy had a lower RRPOF than either of the other two types. We also tested the sensitivity of the RRRPOF between photon and proton therapies to the anatomic position of the ovaries, i.e., proximity to the treatment field (2 ≤ RRRPOF ≤ 10). We found that CRT and IMRT have higher risks of POF than passive-scattering proton radiotherapy (PRT) does, regardless of uncertainties in the ovarian location. Overall, PRT represents a lower RRPOF over the two other modalities.

Effect of inhomogeneity in a patient's body on the accuracy of the pencil beam algorithm in comparison to Monte Carlo

The pencil beam algorithm (PBA) is reasonably accurate and fast. It is, therefore, the primary method used in routine clinical treatment planning for proton radiotherapy; still, it needs to be validated for use in highly inhomogeneous regions. In our investigation of the effect of patient inhomogeneity, PBA was compared with Monte Carlo (MC). A software framework was developed for the MC simulation of radiotherapy based on Geant4. Anatomical sites selected for the comparison were the head/neck, liver, lung and pelvis region. The dose distributions calculated by the two methods in selected examples were compared, as well as a dose volume histogram (DVH) derived from the dose distributions. The comparison of the off-center ratio (OCR) at the iso-center showed good agreement between the PBA and MC, while discrepancies were seen around the distal fall-off regions. While MC showed a fine structure on the OCR in the distal fall-off region, the PBA showed smoother distribution. The fine structures in MC calculation appeared downstream of very low-density regions. Comparison of DVHs showed that most of the target volumes were similarly covered, while some OARs located around the distal region received a higher dose when calculated by MC than the PBA.

Comparison of MCNPX and Geant4 proton energy deposition predictions for clinical use

Several different Monte Carlo codes are currently being used at proton therapy centers to improve upon dose predictions over standard methods using analytical or semi-empirical dose algorithms. There is a need to better ascertain the differences between proton dose predictions from different available Monte Carlo codes. In this investigation Geant4 and MCNPX, the two most-utilized Monte Carlo codes for proton therapy applications, were used to predict energy deposition distributions in a variety of geometries, comprising simple water phantoms, water phantoms with complex inserts and in a voxelized geometry based on clinical CT data. The Gamma analysis was used to evaluate the differences of the predictions between the codes. The results show that in all the cases the agreement was better than clinical acceptance criteria.

Range uncertainties in proton therapy and the role of Monte Carlo simulations

The main advantages of proton therapy are the reduced total energy deposited in the patient as compared to photon techniques and the finite range of the proton beam. The latter adds an additional degree of freedom to treatment planning. The range in tissue is associated with considerable uncertainties caused by imaging, patient setup, beam delivery and dose calculation. Reducing the uncertainties would allow a reduction of the treatment volume and thus allow a better utilization of the advantages of protons. This paper summarizes the role of Monte Carlo simulations when aiming at a reduction of range uncertainties in proton therapy. Differences in dose calculation when comparing Monte Carlo with analytical algorithms are analyzed as well as range uncertainties due to material constants and CT conversion. Range uncertainties due to biological effects and the role of Monte Carlo for in vivo range verification are discussed. Furthermore, the current range uncertainty recipes used at several proton therapy facilities are revisited. We conclude that a significant impact of Monte Carlo dose calculation can be expected in complex geometries where local range uncertainties due to multiple Coulomb scattering will reduce the accuracy of analytical algorithms. In these cases Monte Carlo techniques might reduce the range uncertainty by several mm.

Range uncertainties in proton therapy and the role of Monte Carlo simulations

The main advantages of proton therapy are the reduced total energy deposited in the patient as compared to photon techniques and the finite range of the proton beam. The latter adds an additional degree of freedom to treatment planning. The range in tissue is associated with considerable uncertainties caused by imaging, patient setup, beam delivery and dose calculation. Reducing the uncertainties would allow a reduction of the treatment volume and thus allow a better utilization of the advantages of protons. This paper summarizes the role of Monte Carlo simulations when aiming at a reduction of range uncertainties in proton therapy. Differences in dose calculation when comparing Monte Carlo with analytical algorithms are analyzed as well as range uncertainties due to material constants and CT conversion. Range uncertainties due to biological effects and the role of Monte Carlo for in vivo range verification are discussed. Furthermore, the current range uncertainty recipes used at several proton therapy facilities are revisited. We conclude that a significant impact of Monte Carlo dose calculation can be expected in complex geometries where local range uncertainties due to multiple Coulomb scattering will reduce the accuracy of analytical algorithms. In these cases Monte Carlo techniques might reduce the range uncertainty by several mm.

Clinical implementation of a GPU-based simplified Monte Carlo method for a treatment planning system of proton beam therapy

We implemented the simplified Monte Carlo (SMC) method on graphics processing unit (GPU) architecture under the computer-unified device architecture platform developed by NVIDIA. The GPU-based SMC was clinically applied for four patients with head and neck, lung, or prostate cancer. The results were compared to those obtained by a traditional CPU-based SMC with respect to the computation time and discrepancy. In the CPU- and GPU-based SMC calculations, the estimated mean statistical errors of the calculated doses in the planning target volume region were within 0.5% rms. The dose distributions calculated by the GPU- and CPU-based SMCs were similar, within statistical errors. The GPU-based SMC showed 12.30-16.00 times faster performance than the CPU-based SMC. The computation time per beam arrangement using the GPU-based SMC for the clinical cases ranged 9-67 s. The results demonstrate the successful application of the GPU-based SMC to a clinical proton treatment planning.

GRID-ENABLED TREATMENT PLANNING FOR PROTON THERAPY USING MONTE CARLO SIMULATIONS

Grid computing is an emerging technology that enables computational tasks to be accomplished in a collaborative approach by using a distributed network of computers. The grid approach is especially important for computationally intensive problems that are not tractable with a single computer or even with a small cluster of computers, e.g., radiation transport calculations for cancer therapy. The objective of this work was to extend a Monte Carlo (MC) transport code used for proton radiotherapy to utilize grid computing techniques and demonstrate its promise in reducing runtime from days to minutes. As proof of concept we created the Medical Grid between Texas Tech University and Rice University. Preliminary computational experiments were carried out in the GEANT4 simulation environment for transport of 25 ×106 200 MeV protons in a prostate cancer treatment plan. The simulation speedup was approximately linear; deviations were attributed to the spectrum of parallel runtimes and communication overhead due to Medical Grid computing. The results indicate that ~3 × 105 to 5 × 105 proton events with processor core would result in 65 to 83% efficiency. Extrapolation of our results indicates that about 103 processor cores of the class used here would reduce the MC simulation runtime from 18.3 days to ~1 h.

Application of a fast proton dose calculation algorithm to a thorax geometry

Treatment planning in proton therapy requires the calculation of absorbed dose distributions on beam shaping components and the patient anatomy. Analytical pencil-beam dose algorithms commonly used are not always accurate enough. The Monte Carlo approach is more accurate but extremely computationally intensive. The Fast Dose Calculator, a track-repeating algorithm, has been proposed as an alternative fast and accurate dose calculation. In this work FDC is applied to a proton therapy patient thoracic anatomy.

A GPU implementation of a track-repeating algorithm for proton radiotherapy dose calculations

An essential component in proton radiotherapy is the algorithm to calculate the radiation dose to be delivered to the patient. The most common dose algorithms are fast but they are approximate analytical approaches. However their level of accuracy is not always satisfactory, especially for heterogeneous anatomical areas, like the thorax. Monte Carlo techniques provide superior accuracy; however, they often require large computation resources, which render them impractical for routine clinical use. Track-repeating algorithms, for example the fast dose calculator, have shown promise for achieving the accuracy of Monte Carlo simulations for proton radiotherapy dose calculations in a fraction of the computation time. We report on the implementation of the fast dose calculator for proton radiotherapy on a card equipped with graphics processor units (GPUs) rather than on a central processing unit architecture. This implementation reproduces the full Monte Carlo and CPU-based track-repeating dose calculations within 2%, while achieving a statistical uncertainty of 2% in less than 1 min utilizing one single GPU card, which should allow real-time accurate dose calculations.

Adjustment of the lateral and longitudinal size of scanned proton beam spots using a pre-absorber to optimize penumbrae and delivery efficiency

In scanned-beam proton therapy, the beam spot properties, such as the lateral and longitudinal size and the minimum achievable range, are influenced by beam optics, scattering media and drift spaces in the treatment unit. Currently available spot scanning systems offer few options for adjusting these properties. We investigated a method for adjusting the lateral and longitudinal spot size that utilizes downstream plastic pre-absorbers located near a water phantom. The spot size adjustment was characterized using Monte Carlo simulations of a modified commercial scanned-beam treatment head. Our results revealed that the pre-absorbers can be used to reduce the lateral full width at half maximum (FWHM) of dose spots in water by up to 14 mm, and to increase the longitudinal extent from about 1 mm to 5 mm at residual ranges of 4 cm and less. A large factor in manipulating the lateral spot sizes is the drift space between the pre-absorber and the water phantom. Increasing the drift space from 0 cm to 15 cm leads to an increase in the lateral FWHM from 2.15 cm to 2.87 cm, at a water-equivalent depth of 1 cm. These findings suggest that this spot adjustment method may improve the quality of spot-scanned proton treatments.

An MCNPX Monte Carlo model of a discrete spot scanning proton beam therapy nozzle

PURPOSE

The purposes of this study were to validate a discrete spot scanning proton beam nozzle using the Monte Carlo (MC) code MCNPX and use the MC validated model to investigate the effects of a low-dose envelope, which surrounds the beam's central axis, on measurements of integral depth dose (IDD) profiles.

METHODS

An accurate model of the discrete spot scanning beam nozzle from The University of Texas M. D. Anderson Cancer Center (Houston, Texas) was developed on the basis of blueprints provided by the manufacturer of the nozzle. The authors performed simulations of single proton pencil beams of various energies using the standard multiple Coulomb scattering (MCS) algorithm within the MCNPX source code and a new MCS algorithm, which was implemented in the MCNPX source code. The MC models were validated by comparing calculated in-air and in-water lateral profiles and percentage depth dose profiles for single pencil beams with their corresponding measured values. The models were then further tested by comparing the calculated and measured three-dimensional (3-D) dose distributions. Finally, an IDD profile was calculated with different scoring radii to determine the limitations on the use of commercially available plane-parallel ionization chambers to measure IDD.

RESULTS

The distance to agreement, defined as the distance between the nearest positions of two equivalent distributions with the same value of dose, between measured and simulated ranges was within 0.13 cm for both MCS algorithms. For low and intermediate pencil beam energies, the MC simulations using the standard MCS algorithm were in better agreement with measurements. Conversely, the new MCS algorithm produced better results for high-energy single pencil beams. The IDD profile calculated with cylindrical tallies with an area equivalent to the area of the largest commercially available ionization chamber showed up to 7.8% underestimation of the integral dose in certain depths of the IDD profile.

CONCLUSIONS

The authors conclude that a combination of MCS algorithms is required to accurately reproduce experimental data of single pencil beams and 3-D dose distributions for the scanning beam nozzle. In addition, the MC simulations showed that because of the low-dose envelope, ionization chambers with radii as large as 4.08 cm are insufficient to accurately measure IDD profiles for a 221.8 MeV pencil beam in the scanning beam nozzle.

Monte Carlo simulation of the neutron spectral fluence and dose equivalent for use in shielding a proton therapy vault

Neutron production is of principal concern when designing proton therapy vault shielding. Conventionally, neutron calculations are based on analytical methods, which do not accurately consider beam shaping components and nozzle shielding. The goal of this study was to calculate, using Monte Carlo modeling, the neutron spectral fluence and neutron dose equivalent generated by a realistic proton therapy nozzle and evaluate how these data could be used in shielding calculations. We modeled a contemporary passive scattering proton therapy nozzle in detail with the MCNPX simulation code. The neutron spectral fluence and dose equivalent at various locations in the treatment room were calculated and compared to those obtained from a thick iron target bombarded by parallel proton beams, the simplified geometry on which analytical methods are based. The neutron spectral fluence distributions were similar for both methods, with deeply penetrating high-energy neutrons (E > 10 MeV) being most prevalent along the beam central axis, and low-energy neutrons predominating the neutron spectral fluence in the lateral region. However, unlike the inverse square falloff used in conventional analytical methods, this study shows that the neutron dose equivalent per therapeutic dose in the treatment room decreased with distance approximately following a power law, with an exponent of about -1.63 in the lateral region and -1.73 in the downstream region. Based on the simulated data according to the detailed nozzle modeling, we developed an empirical equation to estimate the neutron dose equivalent at any location and distance in the treatment vault, e.g. for cases in which detailed Monte Carlo modeling is not feasible. We applied the simulated neutron spectral fluence and dose equivalent to a shielding calculation as an example.

Monte Carlo modelling of the treatment line of the Proton Therapy Center in Orsay

This paper presents the main results of a Monte Carlo simulation describing the Orsay Proton Therapy Center (CPO) beam line. The project aimed to obtain a prediction of the dose distribution in a water phantom within 2% accuracy in the dose value and a 2 mm of range. The simulation tool used was MCNPX, version 2.5.0, and included all the elements of the CPO beam line. A new algorithm of multiple Coulomb scattering has been incorporated in MCNPX, resulting in a better prediction of the spatial dose distribution and absolute values of the deposited energy. The simulations of 3D dose profiles in water show a very good agreement with measured data to within 2%. We first performed a comparative analysis of the dosimetry in heterogeneous phantoms between the pencil beam algorithm and MCNPX. The simulations give a better agreement with experimental data compared to the pencil beam approach. In a second phase, we simulated the patient-dependent fields along with the spatial dose distributions in a water phantom. The simulated response of a Pixel chamber located 2 m upstream of the water phantom revealed a good agreement with the measured data to within 1%. The results presented herein support the applicability of Monte Carlo models for absolute dosimetry and for design purposes regarding existing and new beam lines at CPO. This work completes a series of publications reporting the progress in the development of a Monte Carlo simulation tool for the CPO beam line dedicated for the treatment of head and neck tumours.

The risk of developing a second cancer after receiving craniospinal proton irradiation

The purpose of this work was to compare the risk of developing a second cancer after craniospinal irradiation using photon versus proton radiotherapy by means of simulation studies designed to account for the effects of neutron exposures. Craniospinal irradiation of a male phantom was calculated for passively-scattered and scanned-beam proton treatment units. Organ doses were estimated from treatment plans; for the proton treatments, the amount of stray radiation was calculated separately using the Monte Carlo method. The organ doses were converted to risk of cancer incidence using a standard formalism developed for radiation protection purposes. The total lifetime risk of second cancer due exclusively to stray radiation was 1.5% for the passively scattered treatment versus 0.8% for the scanned proton beam treatment. Taking into account the therapeutic and stray radiation fields, the risk of second cancer from intensity-modulated radiation therapy and conventional radiotherapy photon treatments were 7 and 12 times higher than the risk associated with scanned-beam proton therapy, respectively, and 6 and 11 times higher than with passively scattered proton therapy, respectively. Simulations revealed that both passively scattered and scanned-beam proton therapies confer significantly lower risks of second cancers than 6 MV conventional and intensity-modulated photon therapies.

Contemporary Proton Therapy Systems Adequately Protect Patients from Exposure to Stray Radiation

Proton beam therapy has provided safe and effective treatments for a variety of adult cancers. In recent years, there has been increasing interest in utilizing proton therapy for pediatric cancers because it allows better sparing of healthy tissues. Minimizing exposures of normal tissues is especially important in children because they are highly susceptible to consequential late effects, including the development of a radiogenic second cancer, which may occur years or even decades after treatment of the first cancer. While the dosimetric advantage of therapeutic proton beams is well understood, relatively little attention has been paid to the whole-body exposure to stray neutron radiation that is inherent in proton therapy. In this report, we review the physical processes that lead to neutron exposures, discuss the potential for mitigating these exposures using advanced proton beam delivery systems, and present a comparative analysis of predicted second cancer incidence following various external beam therapies. In addition, we discuss uncertainties in the relative biological effectiveness of neutrons for carcinogenesis and the impact that these uncertainties have on second-cancer risk predictions for survivors of adult and childhood cancer who receive proton therapy.

Effective Dose from Stray Radiation for a Patient Receiving Proton Therapy for Liver Cancer

Because of its advantageous depth-dose relationship, proton radiotherapy is an emerging treatment modality for patients with liver cancer. Although the proton dose distribution conforms to the target, healthy tissues throughout the body receive low doses of stray radiation, particularly neutrons that originate in the treatment unit or in the patient. The aim of this study was to calculate the effective dose from stray radiation and estimate the corresponding risk of second cancer fatality for a patient receiving proton beam therapy for liver cancer. Effective dose from stray radiation was calculated using detailed Monte Carlo simulations of a double-scattering proton therapy treatment unit and a voxelized human phantom. The treatment plan and phantom were based on CT images of an actual adult patient diagnosed with primary hepatocellular carcinoma. For a prescribed dose of 60 Gy to the clinical target volume, the effective dose from stray radiation was 370 mSv; 61% of this dose was from neutrons originating outside of the patient while the remaining 39% was from neutrons originating within the patient. The excess lifetime risk of fatal second cancer corresponding to the total effective dose from stray radiation was 1.2%. The results of this study establish a baseline estimate of the stray radiation dose and corresponding risk for an adult patient undergoing proton radiotherapy for liver cancer and provide new evidence to corroborate the suitability of proton beam therapy for the treatment of liver tumors.

Prompt gamma-ray emission from biological tissues during proton irradiation: a preliminary study

In this paper, we present the results of a preliminary study of secondary 'prompt' gamma-ray emission produced by proton-nuclear interactions within tissue during proton radiotherapy. Monte Carlo simulations were performed for mono-energetic proton beams, ranging from 2.5 MeV to 250 MeV, irradiating elemental and tissue targets. Calculations of the emission spectra from different biological tissues and their elemental components were made. Also, prompt gamma rays emitted during delivery of a clinical proton spread-out Bragg peak (SOBP) in a homogeneous water phantom and a water phantom containing heterogeneous tissue inserts were calculated to study the correlation between prompt gamma-ray production and proton dose delivery. The results show that the prompt gamma-ray spectra differ significantly for each type of tissue studied. The relative intensity of the characteristic gamma rays emitted from a given tissue was shown to be proportional to the concentration of each element in that tissue. A strong correlation was found between the delivered SOBP dose distribution and the characteristic prompt gamma-ray production. Based on these results, we discuss the potential use of prompt gamma-ray emission as a method to verify the accuracy and efficacy of doses delivered with proton radiotherapy.