Development and validation of MonteRay, a fast Monte Carlo dose engine for carbon ion beam radiotherapy.

BACKGROUND: Monte Carlo (MC) simulations are considered the gold-standard for accuracy in radiotherapy dose calculation; so far however, no commercial treatment planning system (TPS) provides a fast MC for supporting clinical practice in carbon ion therapy. PURPOSE: To extend and validate the in-house developed fast MC dose engine MonteRay for carbon ion therapy, including physical and biological dose calculation. METHODS: MonteRay is a CPU MC dose calculation engine written in C++ that is capable of simulating therapeutic proton, helium and carbon ion beams. In this work, development steps taken to include carbon ions in MonteRay are presented. Dose distributions computed with MonteRay are evaluated using a comprehensive validation dataset, including various measurements (pristine Bragg peaks, spread out Bragg peaks in water and behind an anthropomorphic phantom) and simulations of a patient plan. The latter includes both physical and biological dose comparisons. Runtimes of MonteRay were evaluated against those of FLUKA MC on a standard benchmark problem. RESULTS: Dosimetric comparisons between MonteRay and measurements demonstrated good agreement. In terms of pristine Bragg peaks, mean errors between simulated and measured integral depth dose distributions were between -2.3% and +2.7%. Comparing SOBPs at 5, 12.5 and 20 cm depth, mean absolute relative dose differences were 0.9%, 0.7% and 1.6% respectively. Comparison against measurements behind an anthropomorphic head phantom revealed mean absolute dose differences of 1.2 % ± 1.1 % $1.2\% \pm 1.1\;\%$ with global 3%/3 mm 3D-γ passing rates of 99.3%, comparable to those previously reached with FLUKA (98.9%). Comparisons against dose predictions computed with the clinical treatment planning tool RayStation 11B for a meningioma patient plan revealed excellent local 1%/1 mm 3D-γ passing rates of 98% for physical and 94% for biological dose. In terms of runtime, MonteRay achieved speedups against reference FLUKA simulations ranging from 14× to 72×, depending on the beam's energy and the step size chosen. CONCLUSIONS: Validations against clinical dosimetric measurements in homogeneous and heterogeneous scenarios and clinical TPS calculations have proven the validity of the physical models implemented in MonteRay. To conclude, MonteRay is viable as a fast secondary MC engine for supporting clinical practice in proton, helium and carbon ion radiotherapy.

Commissioning of Helium Ion Therapy and the First Patient Treatment With Active Beam Delivery.

PURPOSE: Helium ions offer intermediate physical and biological properties to the clinically used protons and carbon ions. This work presents the commissioning of the first clinical treatment planning system (TPS) for helium ion therapy with active beam delivery to prepare the first patients' treatment at the Heidelberg Ion-Beam Therapy Center (HIT). METHODS AND MATERIALS: Through collaboration between RaySearch Laboratories and HIT, absorbed and relative biological effectiveness (RBE)-weighted calculation methods were integrated for helium ion beam therapy with raster-scanned delivery in the TPS RayStation. At HIT, a modified microdosimetric kinetic biological model was chosen as reference biological model. TPS absorbed dose predictions were compared against measurements with several devices, using phantoms of different complexities, from homogeneous to heterogeneous anthropomorphic phantoms. RBE and RBE-weighted dose predictions of the TPS were verified against calculations with an independent RBE-weighted dose engine. The patient-specific quality assurance of the first treatment at HIT using helium ion beam with raster-scanned delivery is presented considering standard patient-specific measurements in a water phantom and 2 independent dose calculations with a Monte Carlo or an analytical-based engine. RESULTS: TPS predictions were consistent with dosimetric measurements and independent dose engines computations for absorbed and RBE-weighted doses. The mean difference between dose measurements to the TPS calculation was 0.2% for spread-out Bragg peaks in water. Verification of the first patient treatment TPS predictions against independent engines for both absorbed and RBE-weighted doses presents differences within 2% in the target and with a maximum deviation of 3.5% in the investigated critical regions of interest. CONCLUSIONS: Helium ion beam therapy has been successfully commissioned and introduced into clinical use. Through comprehensive validation of the absorbed and RBE-weighted dose predictions of the RayStation TPS, the first clinical TPS for helium ion therapy using raster-scanned delivery was employed to plan the first helium patient treatment at HIT.

Development and benchmarking of the first fast Monte Carlo engine for helium ion beam dose calculation: MonteRay.

BACKGROUND: Monte Carlo (MC) simulations are considered the gold-standard for accuracy in radiotherapy dose calculation; however, general purpose MC engines are computationally demanding and require long runtimes. For this reason, several groups have recently developed fast MC systems dedicated mainly to photon and proton external beam therapy, affording both speed and accuracy. PURPOSE: To support research and clinical activities at the Heidelberg Ion-beam Therapy Center (HIT) with actively scanned helium ion beams, this work presents MonteRay, the first fast MC dose calculation engine for helium ion therapy. METHODS: MonteRay is a CPU MC dose calculation engine written in C++, capable of simulating therapeutic proton and helium ion beams. In this work, development steps taken to include helium ion beams in MonteRay are presented. A detailed description of the newly implemented physics models for helium ions, for example, for multiple coulomb scattering and inelastic nuclear interactions, is provided. MonteRay dose computations of helium ion beams are evaluated using a comprehensive validation dataset, including measurements of spread-out Bragg peaks (SOBPs) with varying penetration depths/field sizes, measurements with an anthropomorphic phantom and FLUKA simulations of a patient plan. Improvement in computational speed is demonstrated in comparison against reference FLUKA simulations. RESULTS: Dosimetric comparisons between MonteRay and measurements demonstrated good agreement. Comparing SOBPs at 5, 12.5, and 20 cm depth, mean absolute percent dose differences were 0.7%, 0.7%, and 1.4%, respectively. Comparison against measurements behind an anthropomorphic head phantom revealed mean absolute dose differences of about 1.2% (FLUKA: 1.5%) with per voxel errors ranging from -4.5% to 4.1% (FLUKA: -6% to 3%). Computed global 3%/3 mm 3D-gamma passing rates of ∼99% were achieved, exceeding those previously reported for an analytical dose engine. Comparisons against FLUKA simulations for a patient plan revealed local 2%/2 mm 3D-gamma passing rates of 98%. Compared to FLUKA in voxelized geometries, MonteRay saw run-time reductions ranging from 20× to 60×, depending on the beam's energy. CONCLUSIONS: MonteRay, the first fast MC engine dedicated to helium ion therapy, has been successfully developed with a focus on both speed and accuracy. Validations against dosimetric measurements in homogeneous and heterogeneous scenarios and FLUKA MC calculations have proven the validity of the physical models implemented. Timing comparisons have shown significant speedups between 20 and 60 when compared to FLUKA, making MonteRay viable for clinical routine. MonteRay will support research and clinical practice at HIT, for example, TPS development, validation and treatment design for upcoming clinical trials for raster-scanned helium ion therapy.

Deciphering Time-Dependent DNA Damage Complexity, Repair, and Oxygen Tension: A Mechanistic Model for FLASH-Dose-Rate Radiation Therapy.

PURPOSE: Irradiation with ultrahigh dose rates (FLASH) has reemerged as a promising radiation therapy approach to effectively lower potential damage burden on normal tissue without sacrificing tumor control. However, the large number of recent FLASH studies have been conducted under vastly different experimental conditions and circumstances (ie, investigated biological endpoint, radiation quality, and environmental oxygen level), with unverified biological mechanisms of action and unexplored interplay effect of the main dependencies. To facilitate radiobiological investigation of FLASH phenomena and assessment of clinical applicability, we present an extension of the mechanistic radiobiological model "UNified and VERSatile bio response Engine" (UNIVERSE). METHODS AND MATERIALS: The dynamic (time-dependent) extension of UNIVERSE was developed incorporating fundamental temporal mechanisms necessary for dose-rate effect prediction, ie, DNA damage repair kinetics [DDRK], oxygen depletion and reoxygenation during irradiation. Model performance in various experimental conditions is validated based on a large panel of in vitro and in vivo data from the literature. The effect of dose, dose rate, oxygen tension, tissue-type, beam quality and DDRK is analyzed. RESULTS: UNIVERSE adequately reproduces dose-, dose-rate- and oxygen tension-dependent influence on cell killing. For the studied systems, results indicate that the extent of cell/tissue sparing effect, if present at all, strongly depends on DDRK and beam quality used for reference conventional irradiation. A validated mechanistic framework for predicting clinically relevant endpoints comparing conventional and FLASH high-dose-rate effect has been successfully established, relying on time-dependent processing of radiation-induced damage classes taking variable oxygen tension into account. CONCLUSIONS: Highlighted by UNIVERSE itself, the multidimensional nature of this relative sparing effect using high-dose-rate radiation compared with conventional means underlines the importance of robust quantification of biophysical characteristics and consistent, well-documented experimental conditions both in vitro and in vivo before clinical translation. To further elucidate underlying mechanisms and appraise clinical viability, UNIVERSE can provide reliable prediction for biophysical investigations of radiation therapy using ultrahigh dose rate.

FRoG: An independent dose and LETd prediction tool for proton therapy at ProBeam® facilities.

PURPOSE: Particle therapy is becoming increasingly available world-wide for precise tumor targeting, its favorable depth dose deposition, and increased biological damage to tumor tissue compared to conventional photon therapy. As demand increases for improved robustness and conformality, next-generation secondary dose calculation engines are needed to verify treatment plans independently and provide estimates for clinical decision-making factors, such as dose-averaged linear energy transfer (LETd ) and relative biological effectiveness (RBE). METHOD: FRoG (Fast dose Recalculation on GPU) has been installed and commissioned at the Danish Centre for Particle Therapy (DCPT). FRoG was developed for synchrotron-based facilities and has previously demonstrated good agreement with gold-standard Monte Carlo simulations and measurements. In this work, additions and modifications to FRoG's pencil beam algorithm to support the ion beam delivery with cyclotron-based technology as used at the DCPT, range shifter (RS) implementation, and robustness analysis methods are presented. FRoG dose predictions are compared to measurements and predictions of the clinical treatment planning system (TPS) Eclipse (Varian Medical Systems, Palo Alto, United States of America, CA, v.13.7.16) in both homogenous and heterogeneous scenarios using a solid-water/water and a half-head anthropomorphic phantom, respectively. Additional capabilities of FRoG are explored by performing a plan robustness analysis, analyzing dose and LETd for ten patients. RESULTS: Mid-target measurements in spread-out Bragg Peaks (SOBP) were on average within -0.19% ± 0.30% and ≤0.5% of FRoG predictions for irradiations without and with RS, respectively. Average 3%/2mm 3D γ-analysis passing rates were 99.1% for ~200 patient plan QA comparisons. Measurement with an anthropomorphic head-phantom yielded a γ-passing rate >98%. Overall, maximum target differences in D02% of <2% between the TPS and FRoG were observed for patient plans. The robustness analysis study accounting for range, delivery, and positioning uncertainties revealed small differences in target dose and a maximum LETd VH02% (LETd received by 2% of the volume having dose larger than 1% of maximum dose) values below 10.1 keV/µm to the brain stem. CONCLUSION: We demonstrate that auxiliary dose calculation systems like FRoG can yield excellent agreement to measurements comparable to clinical beam models. Through this work, application of FRoG as a secondary engine at third party cyclotron-based particle treatment facilities is now established for dose verification as well as providing further insight on LETd and variable RBE distributions for protons, currently absent from the standard clinical TPS.

Assessment of RBE-Weighted Dose Models for Carbon Ion Therapy Toward Modernization of Clinical Practice at HIT: In Vitro, in Vivo, and in Patients.

PURPOSE: Present-day treatment planning in carbon ion therapy is conducted with assumptions for a limited number of tissue types and models for effective dose. Here, we comprehensively assess relative biological effectiveness (RBE) in carbon ion therapy and associated models toward the modernization of current clinical practice in effective dose calculation. METHODS: Using 2 human (A549, H460) and 2 mouse (B16, Renca) tumor cell lines, clonogenic cell survival assay was performed for examination of changes in RBE along the full range of clinical-like spread-out Bragg peak (SOBP) fields. Prediction power of the local effect model (LEM1 and LEM4) and the modified microdosimetric kinetic model (mMKM) was assessed. Experimentation and analysis were carried out in the frame of a multidimensional end point study for clinically relevant ranges of physical dose (D), dose-averaged linear energy transfer (LETd), and base-line photon radio-sensitivity (α/ß)x. Additionally, predictions were compared against previously reported RBE measurements in vivo and surveyed in patient cases. RESULTS: RBE model prediction performance varied among the investigated perspectives, with mMKM prediction exhibiting superior agreement with measurements both in vitro and in vivo across the 3 investigated end points. LEM1 and LEM4 performed their best in the highest LET conditions but yielded overestimations and underestimations in low/midrange LET conditions, respectively, as demonstrated by comparison with measurements. Additionally, the analysis of patient treatment plans revealed substantial variability across the investigated models (±20%-30% uncertainty), largely dependent on the selected model and absolute values for input tissue parameters αx and ßx. CONCLUSION: RBE dependencies in vitro, in vivo, and in silico were investigated with respect to various clinically relevant end points in the context of tumor-specific tissue radio-sensitivity assignment and accurate RBE modeling. Discovered model trends and performances advocate upgrading current treatment planning schemes in carbon ion therapy and call for verification via clinical outcome analysis with large patient cohorts.

Dosimetric validation of Monte Carlo and analytical dose engines with raster-scanning 1H, 4He, 12C, and 16O ion-beams using an anthropomorphic phantom.

With high-precision radiotherapy on the rise towards mainstream healthcare, comprehensive validation procedures are essential, especially as more sophisticated technologies emerge. In preparation for the upcoming translation of novel ions, case-/disease-specific ion-beam selection and advanced multi-particle treatment modalities at the Heidelberg Ion-beam Therapy Center (HIT), we quantify the accuracy limits in particle therapy treatment planning under complex heterogeneous conditions for the four ions (1H, 4He, 12C, 16O) using a Monte Carlo Treatment Planning platform (MCTP), an independent GPU-accelerated analytical dose engine developed in-house (FRoG) and the clinical treatment planning system (Syngo RT Planning). Attaching an anthropomorphic half-head Alderson RANDO phantom to entrance window of a dosimetric verification water tank, a cubic target spread-out Bragg peak (SOBP) was optimized using the MCTP to best resolve effects of anatomic heterogeneities on dose homogeneity. Subsequent forward calculations were executed in FRoG and Syngo. Absolute and relative dosimetry was performed in the experimental beam room using 1D and 2D array ionization chamber detectors. Mean absolute percent deviation in dose (|%Δ|) between predictions and PinPoint ionization chamber measurements were within ∼2% for all investigated ions for both MCTP and FRoG. For protons and carbon ions, |%Δ| values were ∼4% for Syngo. For the four ions, 3D-γ analysis (3%/3mm criteria) of FLUKA and FRoG presented mean passing rates of 97.0(±2.4)% and 93.6(±4.2)%. FRoG demonstrated satisfactory agreement with gold standard Monte Carlo simulation and measurement, superior to the commercial system. Our pre-clinical trial landmarks the first measurements taken in anthropomorphic settings for helium, carbon and oxygen ion-beam therapy.

Biophysical modeling and experimental validation of relative biological effectiveness (RBE) for 4He ion beam therapy.

BACKGROUND: Helium (4He) ion beam therapy provides favorable biophysical characteristics compared to currently administered particle therapies, i.e., reduced lateral scattering and enhanced biological damage to deep-seated tumors like heavier ions, while simultaneously lessened particle fragmentation in distal healthy tissues as observed with lighter protons. Despite these biophysical advantages, raster-scanning 4He ion therapy remains poorly explored e.g., clinical translational is hampered by the lack of reliable and robust estimation of physical and radiobiological uncertainties. Therefore, prior to the upcoming 4He ion therapy program at the Heidelberg Ion-beam Therapy Center (HIT), we aimed to characterize the biophysical phenomena of 4He ion beams and various aspects of the associated models for clinical integration. METHODS: Characterization of biological effect for 4He ion beams was performed in both homogenous and patient-like treatment scenarios using innovative models for estimation of relative biological effectiveness (RBE) in silico and their experimental validation using clonogenic cell survival as the gold-standard surrogate. Towards translation of RBE models in patients, the first GPU-based treatment planning system (non-commercial) for raster-scanning 4He ion beams was devised in-house (FRoG). RESULTS: Our data indicate clinically relevant uncertainty of ±5-10% across different model simulations, highlighting their distinct biological and computational methodologies. The in vitro surrogate for highly radio-resistant tissues presented large RBE variability and uncertainty within the clinical dose range. CONCLUSIONS: Existing phenomenological and mechanistic/biophysical models were successfully integrated and validated in both Monte Carlo and GPU-accelerated analytical platforms against in vitro experiments, and tested using pristine peaks and clinical fields in highly radio-resistant tissues where models exhibit the greatest RBE uncertainty. Together, these efforts mark an important step towards clinical translation of raster-scanning 4He ion beam therapy to the clinic.

Proton and helium ion radiotherapy for meningioma tumors: a Monte Carlo-based treatment planning comparison.

BACKGROUND: Due to their favorable physical and biological properties, helium ion beams are increasingly considered a promising alternative to proton beams for radiation therapy. Hence, this work aims at comparing in-silico the treatment of brain and ocular meningiomas with protons and helium ions, using for the first time a dedicated Monte Carlo (MC) based treatment planning engine (MCTP) thoroughly validated both in terms of physical and biological models. METHODS: Starting from clinical treatment plans of four patients undergoing proton therapy with a fixed relative biological effectiveness (RBE) of 1.1 and a fraction dose of 1.8 Gy(RBE), new treatment plans were optimized with MCTP for both protons (with variable and fixed RBE) and helium ions (with variable RBE) under the same constraints derived from the initial clinical plans. The resulting dose distributions were dosimetrically compared in terms of dose volume histograms (DVH) parameters for the planning target volume (PTV) and the organs at risk (OARs), as well as dose difference maps. RESULTS: In most of the cases helium ion plans provided a similar PTV coverage as protons with a consistent trend of superior OAR sparing. The latter finding was attributed to the ability of helium ions to offer sharper distal and lateral dose fall-offs, as well as a more favorable differential RBE variation in target and normal tissue. CONCLUSIONS: Although more studies are needed to investigate the clinical potential of helium ions for different tumour entities, the results of this work based on an experimentally validated MC engine support the promise of this modality with state-of-the-art pencil beam scanning delivery, especially in case of tumours growing in close proximity of multiple OARs such as meningiomas.

Next generation multi-scale biophysical characterization of high precision cancer particle radiotherapy using clinical proton, helium-, carbon- and oxygen ion beams.

The growing number of particle therapy facilities worldwide landmarks a novel era of precision oncology. Implementation of robust biophysical readouts is urgently needed to assess the efficacy of different radiation qualities. This is the first report on biophysical evaluation of Monte Carlo simulated predictive models of prescribed dose for four particle qualities i.e., proton, helium-, carbon- or oxygen ions using raster-scanning technology and clinical therapy settings at HIT. A high level of agreement was found between the in silico simulations, the physical dosimetry and the clonogenic tumor cell survival. The cell fluorescence ion track hybrid detector (Cell-Fit-HD) technology was employed to detect particle traverse per cell nucleus. Across a panel of radiobiological surrogates studied such as late ROS accumulation and apoptosis (caspase 3/7 activation), the relative biological effectiveness (RBE) chiefly correlated with the radiation species-specific spatio-temporal pattern of DNA double strand break (DSB) formation and repair kinetic. The size and the number of residual nuclear γ-H2AX foci increased as a function of linear energy transfer (LET) and RBE, reminiscent of enhanced DNA-damage complexity and accumulation of non-repairable DSB. These data confirm the high relevance of complex DSB formation as a central determinant of cell fate and reliable biological surrogates for cell survival/ RBE. The multi-scale simulation, physical and radiobiological characterization of novel clinical quality beams presented here constitutes a first step towards development of high precision biologically individualized radiotherapy.

Clinical implementation and range evaluation of in vivo PET dosimetry for particle irradiation in patients with primary glioma.

PURPOSE: The physical and biological properties of ion-beams offer various advantages in comparison to conventional radiotherapy, though uncertainties concerning quality assurance are still left. Due to the inverted depth dose profile, range accuracy is of paramount importance. We investigated the range deviations between planning simulation and post-fractional PET/CT measurement from particle therapy in primary glioblastoma. METHODS AND MATERIALS: 20 patients with glioblastoma undergoing particle therapy at our institution were selected. 10 received a proton-boost, 10 a carbon-ion-boost in addition to standard treatment. After two fractions, we performed a PET/CT-scan of the brain. We compared the resulting range deviation based on the Most-likely-shift method between the two measurements, and the measurements with corresponding expectations, calculated with the Monte-Carlo code FLUKA. RESULTS: A patient's two measurements deviated by 0.7mm (±0.7mm). Overall comparison between measurements and simulation resulted in a mean range deviation of 3.3mm (±2.2mm) with significant lower deviations in the (12)C-arm. CONCLUSION: The used planning concepts display the actual dose distributions adequately. The carbon ion group's results are below the used PTV safety margins (3mm). Further adjustments to the simulation are required for proton irradiations. Some anatomical situations require particular attention to ensure highest accuracy and safety.

Implementation and initial clinical experience of offline PET/CT-based verification of scanned carbon ion treatment.

BACKGROUND AND PURPOSE: We report on the implementation of offline PET/CT-based treatment verification at the Heidelberg Ion Beam Therapy Centre (HIT) and present first clinical cases for post-activation measurements after scanned carbon ion irradiation. Key ingredient of this in-vivo treatment verification is the comparison of irradiation-induced patient activation measured by a PET scanner with a prediction simulated by means of Monte Carlo techniques. MATERIAL AND METHODS: At HIT, a commercial full-ring PET/CT scanner has been installed in close vicinity to the treatment rooms. After selected irradiation fractions, the patient either walks to the scanner for acquisition of the activation data or is transported using a shuttle system. The expected activity distribution is obtained from the production of ß(+)-active isotopes simulated by the FLUKA code on the basis of the patient-specific treatment plan, post-processed considering the time course of the respective treatment fraction, the estimated biological washout of the induced activity and a simplified model of the imaging process. RESULTS: We present four patients with different indications of head, head/neck, liver and pelvic tumours. A clear correlation between the measured PET signal and the simulated activity pattern is observed for all patients, thus supporting a proper treatment delivery. In the case of a pelvic tumour patient it was possible to detect minor treatment delivery inaccuracies. CONCLUSIONS: The initial clinical experience proves the feasibility of the implemented strategy for offline confirmation of scanned carbon ion irradiation and therefore constitutes a first step towards a comprehensive PET/CT-based treatment verification in the clinical routine at HIT.

Assessment of early toxicity and response in patients treated with proton and carbon ion therapy at the Heidelberg ion therapy center using the raster scanning technique.

UNLABELLED: PUROPOSE: To asses early toxicity and response in 118 patients treated with scanned ion beams to validate the safety of intensity-controlled raster scanning at the Heidelberg Ion Therapy Center. PATIENTS AND METHODS: Between November 2009 and June 2010, we treated 118 patients with proton and carbon ion radiotherapy (RT) using active beam delivery. The main indications included skull base chordomas and chondrosarcomas, salivary gland tumors, and gliomas. We evaluated early toxicity within 6 weeks after RT and the initial clinical and radiologic response for quality assurance in our new facility. RESULTS: In all 118 patients, few side effects were observed, in particular, no high numbers of severe acute toxicity were found. In general, the patients treated with particle therapy alone showed only a few single side effects, mainly Radiation Therapy Oncology Group/Common Terminology Criteria grade 1. The most frequent side effects and cumulative incidence of single side effects were observed in the head-and-neck patients treated with particle therapy as a boost and photon intensity-modulated RT. The toxicities included common radiation-attributed reactions known from photon RT, including mucositis, dysphagia, and skin erythema. The most predominant imaging responses were observed in patients with high-grade gliomas and those with salivary gland tumors. For skull base tumors, imaging showed a stable tumor outline in most patients. Thirteen patients showed improvement of pre-existing clinical symptoms. CONCLUSIONS: Side effects related to particle treatment were rare, and the overall tolerability of the treatment was shown. The initial response was promising. The data have confirmed the safe delivery of carbon ions and protons at the newly opened Heidelberg facility.