Development and benchmarking of a dose rate engine for raster-scanned FLASH helium ions.

BACKGROUND: Radiotherapy with charged particles at high dose and ultra-high dose rate (uHDR) is a promising technique to further increase the therapeutic index of patient treatments. Dose rate is a key quantity to predict the so-called FLASH effect at uHDR settings. However, recent works introduced varying calculation models to report dose rate, which is susceptible to the delivery method, scanning path (in active beam delivery) and beam intensity. PURPOSE: This work introduces an analytical dose rate calculation engine for raster scanned charged particle beams that is able to predict dose rate from the irradiation plan and recorded beam intensity. The importance of standardized dose rate calculation methods is explored here. METHODS: Dose is obtained with an analytical pencil beam algorithm, using pre-calculated databases for integrated depth dose distributions and lateral penumbra. Dose rate is then calculated by combining dose information with the respective particle fluence (i.e., time information) using three dose-rate-calculation models (mean, instantaneous, and threshold-based). Dose rate predictions for all three models are compared to uHDR helium ion beam (145.7 MeV/u, range in water of approximatively 14.6 cm) measurements performed at the Heidelberg Ion Beam Therapy Center (HIT) with a diamond-detector prototype. Three scanning patterns (scanned or snake-like) and four field sizes are used to investigate the dose rate differences. RESULTS: Dose rate measurements were in good agreement with in-silico generated distributions using the here introduced engine. Relative differences in dose rate were below 10% for varying depths in water, from 2.3 to 14.8 cm, as well as laterally in a near Bragg peak area. In the entrance channel of the helium ion beam, dose rates were predicted within 7% on average for varying irradiated field sizes and scanning patterns. Large differences in absolute dose rate values were observed for varying calculation methods. For raster-scanned irradiations, the deviation between mean and threshold-based dose rate at the investigated point was found to increase with the field size up to 63% for a 10 mm × 10 mm field, while no significant differences were observed for snake-like scanning paths. CONCLUSIONS: This work introduces the first dose rate calculation engine benchmarked to instantaneous dose rate, enabling dose rate predictions for physical and biophysical experiments. Dose rate is greatly affected by varying particle fluence, scanning path, and calculation method, highlighting the need for a consensus among the FLASH community on how to calculate and report dose rate in the future. The here introduced engine could help provide the necessary details for the analysis of the sparing effect and uHDR conditions.

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.

Prevalence and determinants of decision regret in long-term prostate cancer survivors following radical prostatectomy.

BACKGROUND: Patients with localized prostate cancer (PC) are faced with a wide spectrum of therapeutic options at initial diagnosis. Following radical prostatectomy (RP), PC patients may experience regret regarding their initial choice of treatment, especially when oncological and functional outcomes are poor. Impacts of psychosocial factors on decision regret, especially after long-term follow-up, are not well understood. This study aimed to investigate the prevalence and determinants of decision regret in long-term PC survivors following RP. METHODS: 3408 PC survivors (mean age 78.8 years, SD = 6.5) from the multicenter German Familial PC Database returned questionnaires after an average of 16.5 (SD = 3.8) years following RP. The outcome of decision regret concerning the initial choice of RP was assessed with one item from the Decision Regret Scale. Health-related quality of life (HRQoL), PC-anxiety, PSA-anxiety, as well as anxiety and depressive symptoms were considered for independent association with decision regret via multivariable logistic regression. RESULTS: 10.9% (373/3408) of PC survivors reported decision regret. Organ-confined disease at RP (OR 1.39, 95%CI 1.02-1.91), biochemical recurrence (OR 1.34, 1.00-1.80), low HRQoL (OR 1.69,1.28-2.24), depressive symptoms (OR 2.32, 1.52-3.53), and prevalent PSA anxiety (OR 1.88,1.17-3.01) were significantly associated with increased risk of decision regret. Shared decision-making reduced the odds of decision regret by 40% (OR 0.59, 0.41-0.86). CONCLUSIONS: PC survivors may experience decision regret even after 16 years following RP. Promoting shared decision-making in light of both established and novel, potentially less invasive treatments at initial diagnosis may help mitigate long-term regret. Awareness regarding patients showing depressive symptoms or PSA anxiety should be encouraged to identify patients at risk of decision regret in need of additional psychological support.

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.

Spot-scanning hadron arc (SHArc) therapy: A proof of concept using single- and multi-ion strategies with helium, carbon, oxygen, and neon ions.

PURPOSE: To present particle arc therapy treatments using single and multi-ion therapy optimization strategies with helium (4 He), carbon (12 C), oxygen (16 O), and neon (20 Ne) ion beams. METHODS AND MATERIALS: An optimization procedure and workflow were devised for spot-scanning hadron arc therapy (SHArc) treatment planning in the PRECISE (PaRticle thErapy using single and Combined Ion optimization StratEgies) treatment planning system (TPS). Physical and biological beam models were developed for helium, carbon, oxygen, and neon ions via FLUKA MC simulation. SHArc treatments were optimized using both single-ion (12 C, 16 O, or 20 Ne) and multi-ion therapy (16 O+4 He or 20 Ne+4 He) applying variable relative biological effectiveness (RBE) modeling using a modified microdosimetric kinetic model (mMKM) with (α/ß)x values of 2, 5, and 3.1 Gy, respectively, for glioblastoma, pancreatic adenocarcinoma, and prostate adenocarcinoma patient cases. Dose, effective dose, linear energy transfer (LET), and RBE were computed with the GPU-accelerated dose engine FRoG and dosimetric/biophysical attributes were evaluated in the context of conventional particle and photon-based therapies (e.g., volumetric modulated arc therapy [VMAT]). RESULTS: All SHArc plans met the target optimization goals (3GyRBE) and demonstrated increased target conformity and substantially lower low-dose bath to surrounding normal tissues than VMAT. SHArc plans using a singleion species (12 C, 16 O, or 20 Ne) exhibited favorable LET distributions with the highest-LET components centralized in the target volume, with values ranging from ∼80-170 keV/µm, ∼130-220 keV/µm, and ∼180-350 keV/µm for 12 C, 16 O, or 20 Ne, respectively, exceeding mean target LET of conventional particle therapy (12 C:∼55, 16 O:∼75 20 Ne:∼95 keV/µm). Multi-ion therapy with SHArc delivery (SHArcMIT ) provided a similar level of target LET enhancement as SHArc compared to conventional planning, however, with additional benefits of homogenous physical dose and RBE distributions. CONCLUSION: Here, we demonstrate that arc delivery of light and heavy ion beams, using either a single-ion species (12 C, 16 O, or 20 Ne) or combining two ions in a single fraction (16 O+4 He or 20 Ne+4 He) affords enhanced physical and biological distributions (e.g., LET) compared with conventional delivery with photons or particle beams. SHArc marks the first single- and multi-ion arc therapy treatment optimization approach using light and heavy ions.

Biological Dose Optimization for Particle Arc Therapy Using Helium and Carbon Ions.

PURPOSE: To present biological dose optimization for particle arc therapy using helium and carbon ions. METHODS AND MATERIALS: Treatment planning and optimization procedures were developed for spot-scanning hadron arc (SHArc) delivery using the RayStation treatment planning system and FRoG dose engine. The SHArc optimization algorithm is applicable for charged particle beams and determines angle dependencies for spot and energy selection with three main initiatives: (i) achieve standard clinical optimization goals and constraints for target and organs at risk (OARs), (ii) target dose robustness, and (iii) increase linear energy transfer (LET) in the target volume. Three patient cases previously treated at the Heidelberg Ion-beam Therapy Center (HIT) were selected for evaluation of conventional versus arc delivery for the two clinical particle beams (helium [4He] and carbon [12C] ions): glioblastoma, prostate adenocarcinoma, and skull-base chordoma. Biological dose and dose-averaged LET (LETd) distributions for SHArc were evaluated against conventional planning techniques (volumetric modulated arc therapy [VMAT] and 2-field intensity modulated particle therapy) applying the modified microdosimetric kinetic model with (α/ß)x = 2 Gy. Clinical viability and deliverability were assessed via evaluation of plan quality, robustness, and irradiation time. RESULTS: For all investigated patient cases, SHArc treatment optimizations met planning goals and constraints for target coverage and OARs, exhibiting acceptable target coverage and reduced normal tissue volumes, with effective dose >10-GyRBE compared with conventional 2F planning. For carbon ions, LETd was increased in the target volume from ∼40-60 to ∼80-140 keV/µm for SHArc compared with conventional treatments. Favorable LETd distributions were possible with the SHArc approach, with maximum LETd in clinical target volume/gross tumor volume and potential reductions of high-LET regions in normal tissues and OARs. Compared with VMAT, SHArc affords substantial reductions in normal tissue dose (40%-70%). CONCLUSIONS: SHArc therapy offers potential treatment benefits such as increased normal tissue sparing from higher doses >10-GyRBE, enhanced target LETd, and potential reduction in high-LET components in OARs. Findings justify further development of robust SHArc treatment planning toward potential clinical translation.

How can we consider variable RBE and LETd prediction during clinical practice? A pediatric case report at the Normandy Proton Therapy Centre using an independent dose engine.

BACKGROUND: To develop an auxiliary GPU-accelerated proton therapy (PT) dose and LETd engine for the IBA Proteus®ONE PT system. A pediatric low-grade glioma case study is reported using FRoG during clinical practice, highlighting potential treatment planning insights using variable RBE dose (DvRBE) and LETd as indicators for clinical decision making in PT. METHODS: The physics engine for FRoG has been modified for compatibility with Proteus®ONE PT centers. Subsequently, FRoG was installed and commissioned at NPTC. Dosimetric validation was performed against measurements and the clinical TPS, RayStation (RS-MC). A head patient cohort previously treated at NPTC was collected and FRoG forward calculations were compared against RS-MC for evaluation of 3D-Γ analysis and dose volume histogram (DVH) results. Currently, treatment design at NPTC is supported with fast variable RBE and LETd calculation and is reported in a representative case for pediatric low-grade glioma. RESULTS: Simple dosimetric tests against measurements of iso-energy layers and spread-out Bragg Peaks in water verified accuracy of FRoG and RS-MC. Among the patient cohort, average 3D-Γ applying 2%/2 mm, 3%/1.5 mm and 5%/1 mm were > 97%. DVH metrics for targets and OARs between FRoG and RayStation were in good agreement, with ∆D50,CTV and ∆D2,OAR both ⪅1%. The pediatric case report demonstrated implications of different beam arrangements on DvRBE and LETd distributions. From initial planning in RayStation sharing identical optimization constraints, FRoG analysis led to plan selection of the most conservative approach, i.e., minimized DvRBE,max and LETd,max in OARs, to avoid optical system toxicity effects (i.e., vision loss). CONCLUSION: An auxiliary dose calculation system was successfully integrated into the clinical workflow at a Proteus®ONE IBA facility, in excellent agreement with measurements and RS-MC. FRoG may lead to further insight on DvRBE and LETd implications to help clinical decision making, better understand unexpected toxicities and establish novel clinical procedures with metrics currently absent from the standard clinical TPS.

Combined DNA Damage Repair Interference and Ion Beam Therapy: Development, Benchmark, and Clinical Implications of a Mechanistic Biological Model.

PURPOSE: Our purpose was to develop a mechanistic model that describes and predicts radiation response after combined DNA damage repair interference (DDRi) and particle radiation therapy. METHODS AND MATERIALS: The heterogeneous dose distributions of protons and 4He ions were implemented into the "UNIfied and VERSatile bio-response Engine" (UNIVERSE). Predictions for monoenergetic and mixed fields over clinically relevant dose and linear energy transfer range were compared with experimental in vitro survival data measured in this work as well as data available in the literature, including different cell lines and DDR interferences. Ultimately, UNIVERSE predictions were investigated in a patient plan. RESULTS: UNIVERSE accurately predicts survival of cell lines with and without DDRi in clinical settings of ion beam therapy based only on 3 parameters derived from photon data. With increasing dose or linear energy transfer, the radiosensitizing effect of DDRi decreases, resulting in diminished relative biological effect of ion beam radiation for cells subjected to DDRi in comparison to cells that are not. Similar trends were observed in patient plan recalculations; however, this analysis also suggests that DDRi + particle radiation therapy may better preserve the therapeutic window in comparison to DDRi + photon radiation therapy. CONCLUSIONS: The presented framework represents the first mechanistic model of combined DDRi and particle radiation therapy comprehensively benchmarked in clinically relevant scenarios and a step toward more personalized treatment. It reveals potential differences between DDRi + photon radiation therapy versus DDRi + particle radiation therapy, which have not been described so far. UNIVERSE could aid in appraising the clinical viability of combined administration of radiosensitizing drugs and charged particle therapy, as well as the identification of patients with known DDR deficiencies in the tumor who might benefit from therapy with light ions, freeing limited space at heavy ion therapy centers.

FRoG dose computation meets Monte Carlo accuracy for proton therapy dose calculation in lung.

PURPOSE: To benchmark and evaluate the clinical viability of novel analytical GPU-accelerated and CPU-based Monte Carlo (MC) dose-engines for spot-scanning intensity-modulated-proton-therapy (IMPT) towards the improvement of lung cancer treatment. METHODS: Nine patient cases were collected from the CNAO clinical experience and The Cancer Imaging Archive-4D-Lung-Database for in-silico study. All plans were optimized with 2 orthogonal beams in RayStation (RS) v.8. Forward calculations were performed with FRoG, an independent dose calculation system using a fast robust approach to the pencil beam algorithm (PBA), RS-MC (CPU for v.8) and general-purpose MC (gp-MC). Dosimetric benchmarks were acquired via irradiation of a lung-like phantom and ionization chambers for both a single-field-uniform-dose (SFUD) and IMPT plans. Dose-volume-histograms, dose-difference and γ-analyses were conducted. RESULTS: With respect to reference gp-MC, the average dose to the GTV was 1.8% and 2.3% larger for FRoG and the RS-MC treatment planning system (TPS). FRoG and RS-MC showed a local γ-passing rate of ~96% and ~93%. Phantom measurements confirmed FRoG's high accuracywith a deviation < 0.1%. CONCLUSIONS: Dose calculation performance using the GPU-accelerated analytical PBA, MC-TPS and gp-MC code were well within clinical tolerances. FRoG predictions were in good agreement with both the full gp-MC and experimental data for proton beams optimized for thoracic dose calculations. GPU-accelerated dose-engines like FRoG may alleviate current issues related to deficiencies in current commercial analytical proton beam models. The novel approach to the PBA implemented in FRoG is suitable for either clinical TPS or as an auxiliary dose-engine to support clinical activity for lung patients.

Spot-Scanning Hadron Arc (SHArc) Therapy: A Study With Light and Heavy Ions.

PURPOSE: To evaluate the clinical potential of spot-scanning hadron arc (SHArc) therapy with a heavy-ion gantry. METHODS AND MATERIALS: A series of in silico studies was conducted via treatment plan optimization in FRoG and the RayStation TPS to compare SHArc therapy against reference plans using conventional techniques with single, parallel-opposed, and 3-field configurations for 3 clinical particle beams (protons [p], helium [4He], and carbon [12C] ions). Tests were performed on water-equivalent cylindrical phantoms for simple targets and clinical-like scenarios with an organ-at-risk in proximity of the target. Effective dose and dose-averaged linear energy transfer (LETD) distributions for SHArc were evaluated against conventional planning techniques applying the modified microdosimetric kinetic model for considering bio-effect with (α/ß)x = 2 Gy. A model for hypoxia-induced tumor radio-resistance was developed for particle therapy with dependence on oxygen concentration and particle species/energy (Zeff/ß)2 to investigate the impact on effective dose. RESULTS: SHArc plans exhibited similar target coverage with unique treatment attributes and distributions compared with conventional planning, with carbon ions demonstrating the greatest potential for tumor control and normal tissue sparing among the arc techniques. All SHArc plans exhibited a low-dose bath outside the target volume with a reduced maximum dose in normal tissues compared with single, parallel-opposed, and 3-field configuration plans. Moreover, favorable LETD distributions were made possible using the SHArc approach, with maximum LETD in the r = 5 mm tumor core (~8 keVµm-1, ~30 keVµm-1, and ~150 keVµm-1 for p, 4He, and 12C ions, respectively) and reductions of high-LET regions in normal tissues and organs-at-risk compared with static treatment beam delivery. CONCLUSION: SHArc therapy offers potential treatment benefits such as increased normal tissue sparing. Without explicit consideration of oxygen concentration during treatment planning and optimization, SHArc-C may mitigate tumor hypoxia-induced loss of efficacy. Findings justify further development of robust SHArc treatment planning toward potential clinical translation.

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.

Experimental comparison of clinically used ion beams for imaging applications using a range telescope.

In particle therapy, the x-ray based treatment planning converting photon attenuation values to relative stopping power ratio (RSP) introduces clinically relevant range uncertainties. Recently, novel imaging technologies using transmission ion beams have been investigated to directly assess the water equivalent thickness (WET) of tissue, showing improved accuracy in RSP reconstruction, while potentially reducing the imaging dose. Due to their greater availability, protons have been mostly used for ion imaging. To this end, in this work, the influence of three ion species (protons, helium and carbon ions) on the image quality of radiographic WET retrieval has been explored with a dedicated experimental setup and compared to Monte Carlo (MC) simulations. Three phantom setups with different tissue interfaces and features have been irradiated with clinically validated proton, helium and carbon ion pencil beams under comparable imaging dose and beam settings at the Heidelberg Ion-Beam Therapy Center. Ion radiographies (iRADs) were acquired with an integration mode detector, that functions as a range telescope with 61 parallel plate ionization chambers. For comparison, experiments were reproduced in-silico with FLUKA MC simulations. Carbon ions provide iRADs with highest image quality in terms of normalized root mean square error, followed by helium ions and protons. All ions show similar capabilities of resolving WET for the considered phantoms, as shown by the similar average relative error < 3%. Besides for the slab phantom, MC simulations yielded better results than the experiment, indicating potential improvement of the experimental setup. Our results showed that the ability to resolve the WET is similar for all particles, intrinsically limited by the granularity of the detector system. While carbon ions are best suited for acquiring iRADs with the investigated integration mode detector, helium ions are put forward as a less technical challenging alternative.

Development and Validation of Single Field Multi-Ion Particle Therapy Treatments.

PURPOSE: To develop and validate combined ion-beam with constant relative biological effectiveness (RBE) (CICR) particle therapy in single field arrangements for improved treatment efficacy, robustness, and normal tissue sparing. METHODS AND MATERIALS: The PRECISE (PaRticle thErapy using single and Combined Ion optimization StratEgies) treatment planning system was developed to investigate clinical viability of CICR treatments. Single-field uniform dose (SFUD) with a single ion (proton [p], helium [He], or carbon [C]) and CICR (C-p and C-He) treatments were generated for 3 patient cases with a clinically prescribed dose of 3 Gy (RBE) per fraction. Spread-out Bragg peak plans were irradiated in homogenous and clinical-like settings using an anthropomorphic head phantom. A dosimetric and biological verification of CICRC-p treatments using a murine glioma cell line (GL261) was performed. RESULTS: CICR treatment plans for the 3 patients presented highly uniform physical dose while reducing high dose-averaged linear energy transfer gradients compared with carbon ions alone. When considering uncertainty in tissue parameter (α/ß)x assignment and RBE modeling, the CICRC-p treatment exhibited enhanced biophysical stability within the target volume, similar to protons alone. CICR treatments reduced dose to normal tissue surrounding the target, exhibiting similar or improved dosimetric features compared with SFUDHe. For both CICRC-p and SFUD treatments, measurements verified the planned dose in the target within ∼3%. Planned versus measured target RBE values were 1.38 ± 0.02 and 1.39 ± 0.07 (<1% deviation), respectively, for the CICRC-p treatment in heterogenous settings. CONCLUSIONS: Here, we demonstrate that by combining 2 (or more) ions in a single field arrangement, more robust biological and more conformal dose distributions can be delivered compared with conventional particle therapy treatment planning. This work constitutes the first dosimetric and biological verification of multi-ion particle therapy in homogeneous as well as heterogenous settings.

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.

FRoG-A New Calculation Engine for Clinical Investigations with Proton and Carbon Ion Beams at CNAO.

A fast and accurate dose calculation engine for hadrontherapy is critical for both routine clinical and advanced research applications. FRoG is a graphics processing unit (GPU)-based forward calculation tool developed at CNAO (Centro Nazionale di Adroterapia Oncologica) and at HIT (Heidelberg Ion Beam Therapy Center) for fast and accurate calculation of both physical and biological dose. FRoG calculation engine adopts a triple Gaussian parameterization for the description of the lateral dose distribution. FRoG provides dose, dose-averaged linear energy transfer, and biological dose-maps, -profiles, and -volume-histograms. For the benchmark of the FRoG calculation engine, using the clinical settings available at CNAO, spread-out Bragg peaks (SOBPs) and patient cases for both proton and carbon ion beams have been calculated and compared against FLUKA Monte Carlo (MC) predictions. In addition, FRoG patient-specific quality assurance (QA) has been performed for twenty-five proton and carbon ion fields. As a result, for protons, biological dose values, using a relative biological effectiveness (RBE) of 1.1, agree on average with MC within ~1% for both SOBPs and patient plans. For carbon ions, RBE-weighted dose (DRBE) agreement against FLUKA is within ~2.5% for the studied SOBPs and patient plans. Both MKM (Microdosimetric Kinetic Model) and LEM (Local Effect Model) DRBE are implemented and tested in FRoG to support the NIRS (National Institute of Radiological Sciences)-based to LEM-based biological dose conversion. FRoG matched the measured QA dosimetric data within ~2.0% for both particle species. The typical calculation times for patients ranged from roughly 1 to 4 min for proton beams and 3 to 6 min for carbon ions on a NVIDIA® GeForce® GTX 1080 Ti. This works demonstrates FRoG's potential to bolster clinical activity with proton and carbon ion beams at CNAO.

Fast robust dose calculation on GPU for high-precision 1H, 4He, 12C and 16O ion therapy: the FRoG platform.

Radiotherapy with protons and heavier ions landmarks a novel era in the field of high-precision cancer therapy. To identify patients most benefiting from this technologically demanding therapy, fast assessment of comparative treatment plans utilizing different ion species is urgently needed. Moreover, to overcome uncertainties of actual in-vivo physical dose distribution and biological effects elicited by different radiation qualities, development of a reliable high-throughput algorithm is required. To this end, we engineered a unique graphics processing unit (GPU) based software architecture allowing rapid and robust dose calculation. FRoG, Fast Recalculation on GPU, currently operates with four particle beams available at Heidelberg Ion Beam Therapy center, i.e., raster-scanning proton (1H), helium (4He), carbon (12C) and oxygen ions (16O). FRoG enables comparative analysis of different models for estimation of physical and biological effective dose in 3D within minutes and in excellent agreement with the gold standard Monte Carlo (MC) simulation. This is a crucial step towards development of next-generation patient specific radiotherapy.

Comparative Monte Carlo study on the performance of integration- and list-mode detector configurations for carbon ion computed tomography.

Ion beam therapy offers the possibility of a highly conformal tumor-dose distribution; however, this technique is extremely sensitive to inaccuracies in the treatment procedures. Ambiguities in the conversion of Hounsfield units of the treatment planning x-ray CT to relative stopping power (RSP) can cause uncertainties in the estimated ion range of up to several millimeters. Ion CT (iCT) represents a favorable solution allowing to directly assess the RSP. In this simulation study we investigate the performance of the integration-mode configuration for carbon iCT, in comparison with a single-particle approach under the same set-up. The experimental detector consists of a stack of 61 air-filled parallel-plate ionization chambers, interleaved with 3 mm thick PMMA absorbers. By means of Monte Carlo simulations, this design was applied to acquire iCTs of phantoms of tissue-equivalent materials. An optimization of the acquisition parameters was performed to reduce the dose exposure, and the implications of a reduced absorber thickness were assessed. In order to overcome limitations of integration-mode detection in the presence of lateral tissue heterogeneities a dedicated post-processing method using a linear decomposition of the detector signal was developed and its performance was compared to the list-mode acquisition. For the current set-up, the phantom dose could be reduced to below 30 mGy with only minor image quality degradation. By using the decomposition method a correct identification of the components and a RSP accuracy improvement of around 2.0% was obtained. The comparison of integration- and list-mode indicated a slightly better image quality of the latter, with an average median RSP error below 1.8% and 1.0%, respectively. With a decreased absorber thickness a reduced RSP error was observed. Overall, these findings support the potential of iCT for low dose RSP estimation, showing that integration-mode detectors with dedicated post-processing strategies can provide a RSP accuracy comparable to list-mode configurations.

Analysis of aldehydes via headspace SPME with on-fiber derivatization to their O-(2,3,4,5,6-pentafluorobenzyl)oxime derivatives and comprehensive 2D-GC-MS.

A method for the analysis of the homologous series of alkanals, (E)-2-alkenals, and (E,E)-2,4-alkadienals is described utilizing a headspace solid-phase microextraction (HS-SPME) step and on-fiber derivatization with O-(2,3,4,5,6-pentafluorobenzyl)hydroxylamine (PFBHA) hydrochloride. Oxime derivatives formed on the fiber are desorbed in the gas chromatographic injector and analyzed by comprehensive 2-D GC coupled to quadrupole MS (GC x GC-qMS). Selecting specific fragment ions within the electron impact mass spectra of the oxime derivatives provides a suitable method for the target analysis of these aldehyde classes, which furthermore benefits from the increased separation efficiency by GC x GC. The analysis of higher molecular weight aldehydes is described in wine and grape seed oil as examples. Quantification of the aldehydes utilizes a stable isotope dilution analysis (SIDA) assay with octan-d(16)-al as isotopomeric internal standard. Besides the selectivity and sensitivity of aldehyde analysis using PFBHA derivatives, critical aspects on background level contamination and repeatability of the sample preparation method are discussed. Optimization of GC x GC-qMS parameters allowed a considerable saving of the cryogenic medium, involving additional (unmodulated) conditioning runs, rendering the method more amenable to routine analysis.

Analysis of carbonyl compounds via headspace solid-phase microextraction with on-fiber derivatization and gas chromatographic-ion trap tandem mass spectrometric determination of their O-(2,3,4,5,6-pentafluorobenzyl)oxime derivatives.

An improved method for the analysis of carbonyls is described utilizing a headspace solid-phase microextraction (HS-SPME) step and on-fiber derivatization with O-(2,3,4,5,6-pentafluorobenzyl)hydroxylamine (PFBHA) hydrochloride. Thermal desorption of the oxime derivatives formed on the fiber is followed by gas chromatographic separation coupled to an ion trap tandem mass spectrometer (GC-ITMS). Selecting specific fragment ions within the electron ionization (EI(+)) mass spectra of these oxime derivatives as precursor ions for MS-MS fragmentation provides a suitable method for the target analysis of individual carbonyl classes, such as alkanals, (E)-2-alkenals, (E,E)-2,4-alkadienals, and others. Retention indices on polar as well as on apolar stationary phases along with EI(+) mass spectra patterns are presented for a large set of oxime derivatives, giving valuable information needed for unambiguous assignment of substances in complex sample matrices. The fast sample preparation and derivatization step via HS-SPME can be automated and is applicable to a variety of biological samples and foodstuffs, allowing rapid and sensitive screening analyses of important aldehydic biomarkers and aroma active compounds.