Evaluation of Helium Ion Radiotherapy in Combination with Gemcitabine in Pancreatic Cancer In Vitro.

BACKGROUND: Pancreatic cancer is one of the most aggressive and lethal cancers. New treatment strategies are highly warranted. Particle radiotherapy could offer a way to overcome the radioresistant nature of pancreatic cancer because of its biological and physical characteristics. Within particles, helium ions represent an attractive therapy option to achieve the highest possible conformity while at the same time protecting the surrounding normal tissue. The aim of this study was to evaluate the cytotoxic efficacy of helium ion irradiation in pancreatic cancer in vitro. METHODS: Human pancreatic cancer cell lines AsPC-1, BxPC-3 and Panc-1 were irradiated with photons and helium ions at various doses and treated with gemcitabine. Photon irradiation was performed with a biological cabin X-ray irradiator, and helium ion irradiation was performed with a spread-out Bragg peak using the raster scanning technique at the Heidelberg Ion Beam Therapy Center (HIT). The cytotoxic effect on pancreatic cancer cells was measured with clonogenic survival. The survival curves were compared to the predicted curves that were calculated via the modified microdosimetric kinetic model (mMKM). RESULTS: The experimental relative biological effectiveness (RBE) of helium ion irradiation ranged from 1.0 to 1.7. The predicted survival curves obtained via mMKM calculations matched the experimental survival curves. Mainly additive cytotoxic effects were observed for the cell lines AsPC-1, BxPC-3 and Panc-1. CONCLUSION: Our results demonstrate the cytotoxic efficacy of helium ion radiotherapy in pancreatic cancer in vitro as well as the capability of mMKM calculation and its value for biological plan optimization in helium ion therapy for pancreatic cancer. A combined treatment of helium irradiation and chemotherapy with gemcitabine leads to mainly additive cytotoxic effects in pancreatic cancer cell lines. The data generated in this study may serve as the radiobiological basis for future experimental and clinical works using helium ion radiotherapy in pancreatic cancer treatment.

Robustness of carbon-ion radiotherapy against DNA damage repair associated radiosensitivity variation based on a biophysical model.

BACKGROUND: Interpatient variation of tumor radiosensitivity is rarely considered during the treatment planning process despite its known significance for the therapeutic outcome. PURPOSE: To apply our mechanistic biophysical model to investigate the biological robustness of carbon ion radiotherapy (CIRT) against DNA damage repair interference (DDRi) associated patient-to-patient variability in radiosensitivity and its potential clinical advantages against conventional radiotherapy approaches. METHODS AND MATERIALS: The "UNIfied and VERSatile bio response Engine" (UNIVERSE) was extended by carbon ions and its predictions were compared to a panel of in vitro and in vivo data including various endpoints and DDRi settings within clinically relevant dose and linear energy transfer (LET) ranges. The implications of UNIVERSE predictions were then assessed in a clinical patient scenario considering DDRi variance. RESULTS: UNIVERSE tests well against the applied benchmarks. While in vitro survival curves were predicted with an R2 > 0.92, deviations from in vivo RBE data were less than 5.6% The conducted paradigmatic patient plan study implies a markedly reduced significance of DDRi based radiosensitivity variability in CIRT (13% change of D 50 ${{D}_{50}}$ in target) compared to conventional radiotherapy (62%) and that boosting the LET within the target further amplifies this robustness of CIRT (8%). In the case of heightened tumor radiosensitivity, a dose de-escalation strategy for photons allows a reduction of the maximum effective dose within the normal tissue (NT) from a D 2 ${{D}_2}$ of 2.65 to 1.64 Gy, which lies below the level found for CIRT ( D 2 ${{D}_2}$  = 2.41 Gy) for the analyzed plan and parameters. However, even after de-escalation, the integral effective dose in the NT is found to be substantially higher for conventional radiotherapy in comparison to CIRT ( D m e a n ${{D}_{mean}}$ of 0.75, 0.46, and 0.24 Gy for the conventional plan, its de-escalation and CIRT, respectively). CONCLUSIONS: The framework offers adequate predictions of in vitro and in vivo radiation effects of CIRT while allowing the consideration of DRRi based solely on parameters derived from photon data. The results of the patient planning study underline the potential of CIRT to minimize important sources of interpatient divergence in therapy outcome, especially when combined with techniques that allow to maximize the LET within the tumor. Despite the potential of de-escalation strategies for conventional radiotherapy to reduce the maximum effective dose in the NT, CIRT appears to remain a more favorable option due to its ability to reduce the integral effective dose within the NT.

Understanding Relative Biological Effectiveness and Clinical Outcome of Prostate Cancer Therapy Using Particle Irradiation: Analysis of Tumor Control Probability With the Modified Microdosimetric Kinetic Model.

PURPOSE: Recent experimental studies and clinical trial results might indicate that-at least for some indications-continued use of the mechanistic model for relative biological effectiveness (RBE) applied at carbon ion therapy facilities in Europe for several decades (LEM-I) may be unwarranted. We present a novel clinical framework for prostate cancer treatment planning and tumor control probability (TCP) prediction based on the modified microdosimetric kinetic model (mMKM) for particle therapy. METHODS AND MATERIALS: Treatment plans of 91 patients with prostate tumors (proton: 46, carbon ions: 45) applying 66 GyRBE [RBE = 1.1 for protons and LEM-I, (α/ß)x = 2.0 Gy, for carbon ions] in 20 fractions were recalculated using mMKM [(α/ß)x = 3.1 Gy]). Based solely on the response data of photon-irradiated patient groups stratified according to risk and usage of androgen deprivation therapy, we derived parameters for an mMKM-based Poisson-TCP model. Subsequently, new carbon and helium ion plans, adhering to prescribed biological dose criteria, were generated. These were systematically compared with the clinical experience of Japanese centers employing an analogous fractionation scheme and existing proton plans. RESULTS: mMKM predictions suggested significant biological dose deviation between the proton and carbon ion arms. Patients irradiated with protons received (3.25 ± 0.08) GyRBEmMKM/Fx, whereas patients treated with carbon ions received(2.51 ± 0.05) GyRBEmMKM/Fx. TCP predictions were (86 ± 3)% for protons and (52 ± 4)% for carbon ions, matching the clinical outcome of 85% and 50%. Newly optimized carbon ion plans, guided by the mMKM/TCP model, effectively replicated clinical data from Japanese centers. Using mMKM, helium ions exhibited similar target coverage as proton and carbon ions and improved rectum and bladder sparing compared with proton. CONCLUSIONS: Our mMKM-based model for prostate cancer treatment planning and TCP prediction was validated against clinical data for proton and carbon ion therapy, and its application was extended to helium ion therapy. Based on the data presented in this work, mMKM seems to be a good candidate for clinical biological calculations in carbon ion therapy for prostate cancer.

Exploring Helium Ions' Potential for Post-Mastectomy Left-Sided Breast Cancer Radiotherapy.

Proton therapy presents a promising modality for treating left-sided breast cancer due to its unique dose distribution. Helium ions provide increased conformality thanks to a reduced lateral scattering. Consequently, the potential clinical benefit of both techniques was explored. An explorative treatment planning study involving ten patients, previously treated with VMAT (Volumetric Modulated Arc Therapy) for 50 Gy in 25 fractions for locally advanced, node-positive breast cancer, was carried out using proton pencil beam therapy with a fixed relative biological effectiveness (RBE) of 1.1 and helium therapy with a variable RBE described by the mMKM (modified microdosimetric kinetic model). Results indicated that target coverage was improved with particle therapy for both the clinical target volume and especially the internal mammary lymph nodes compared to VMAT. Median dose value analysis revealed that proton and helium plans provided lower dose on the left anterior descending artery (LAD), heart, lungs and right breast than VMAT. Notably, helium therapy exhibited improved ipsilateral lung sparing over protons. Employing NTCP models as available in the literature, helium therapy showed a lower probability of grade ≤ 2 radiation pneumonitis (22% for photons, 5% for protons and 2% for helium ions), while both proton and helium ions reduce the probability of major coronary events with respect to VMAT.

Diamond detectors for dose and instantaneous dose-rate measurements for ultra-high dose-rate scanned helium ion beams.

BACKGROUND: The possible emergence of the FLASH effect-the sparing of normal tissue while maintaining tumor control-after irradiations at dose-rates exceeding several tens of Gy per second, has recently spurred a surge of studies attempting to characterize and rationalize the phenomenon. Investigating and reporting the dose and instantaneous dose-rate of ultra-high dose-rate (UHDR) particle radiotherapy beams is crucial for understanding and assessing the FLASH effect, towards pre-clinical application and quality assurance programs. PURPOSE: The purpose of the present work is to investigate a novel diamond-based detector system for dose and instantaneous dose-rate measurements in UHDR particle beams. METHODS: Two types of diamond detectors, a microDiamond (PTW 60019) and a diamond detector prototype specifically designed for operation in UHDR beams (flashDiamond), and two different readout electronic chains, were investigated for absorbed dose and instantaneous dose-rate measurements. The detectors were irradiated with a helium beam of 145.7 MeV/u under conventional and UHDR delivery. Dose-rate delivery records by the monitoring ionization chamber and diamond detectors were studied for single spot irradiations. Dose linearity at 5 cm depth and in-depth dose response from 2 to 16 cm were investigated for both measurement chains and both detectors in a water tank. Measurements with cylindrical and plane-parallel ionization chambers as well as Monte-Carlo simulations were performed for comparisons. RESULTS: Diamond detectors allowed for recording the temporal structure of the beam, in good agreement with the one obtained by the monitoring ionization chamber. A better time resolution of the order of few µs was observed as compared to the approximately 50 µs of the monitoring ionization chamber. Both diamonds detectors show an excellent linearity response in both delivery modalities. Dose values derived by integrating the measured instantaneous dose-rates are in very good agreement with the ones obtained by the standard electrometer readings. Bragg peak curves confirmed the consistency of the charge measurements by the two systems. CONCLUSIONS: The proposed novel dosimetric system allows for a detailed investigation of the temporal evolution of UHDR beams. As a result, reliable and accurate determinations of dose and instantaneous dose-rate are possible, both required for a comprehensive characterization of UHDR beams and relevant for FLASH effect assessment in clinical treatments.

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.

Characteristics of breathing-adapted gating using surface guidance for use in particle therapy: A phantom-based end-to-end test from CT simulation to dose delivery.

To account for intra-fractional tumor motion during dose delivery in radiotherapy, various treatment strategies are clinically implemented such as breathing-adapted gating and irradiating the tumor during specific breathing phases. In this work, we present a comprehensive phantom-based end-to-end test of breathing-adapted gating utilizing surface guidance for use in particle therapy. A commercial dynamic thorax phantom was used to reproduce regular and irregular breathing patterns recorded by the GateRT respiratory monitoring system. The amplitudes and periods of recorded breathing patterns were analysed and compared to planned patterns (ground-truth). In addition, the mean absolute deviations (MAD) and Pearson correlation coefficients (PCC) between the measurements and ground-truth were assessed. Measurements of gated and non-gated irradiations were also analysed with respect to dosimetry and geometry, and compared to treatment planning system (TPS). Further, the latency time of beam on/off was evaluated. Compared to the ground-truth, measurements performed with GateRT showed amplitude differences between 0.03 ± 0.02 mm and 0.26 ± 0.03 mm for regular and irregular breathing patterns, whilst periods of both breathing patterns ranged with a standard deviation between 10 and 190 ms. Furthermore, the GateRT software precisely acquired breathing patterns with a maximum MAD of 0.30 ± 0.23 mm. The PCC constantly ranged between 0.998 and 1.000. Comparisons between TPS and measured dose profiles indicated absolute mean dose deviations within institutional tolerances of ±5%. Geometrical beam characteristics also varied within our institutional tolerances of 1.5 mm. The overall time delays were <60 ms and thus within both recommended tolerances published by ESTRO and AAPM of 200 and 100 ms, respectively. In this study, a non-invasive optical surface-guided workflow including image acquisition, treatment planning, patient positioning and gated irradiation at an ion-beam gantry was investigated, and shown to be clinically viable. Based on phantom measurements, our results show a clinically-appropriate spatial, temporal, and dosimetric accuracy when using surface guidance in the clinical setting, and the results comply with international and institutional guidelines and tolerances.

Proton Therapy for Advanced Juvenile Nasopharyngeal Angiofibroma.

PURPOSE: To provide the first report on proton radiotherapy (PRT) in the management of advanced nasopharyngeal angiofibroma (JNA) and evaluate potential benefits compared to conformal photon therapy (XRT). METHODS: We retrospectively reviewed 10 consecutive patients undergoing PRT for advanced JNA in a definitive or postoperative setting with a relative biological effectiveness weighted dose of 45 Gy in 25 fractions between 2012 and 2022 at the Heidelberg Ion Beam Therapy Center. Furthermore, dosimetric comparisons and risk estimations for short- and long-term radiation-induced complications between PRT plans and helical XRT plans were conducted. RESULTS: PRT was well tolerated, with only low-grade acute toxicities (CTCAE I-II) being reported. The local control rate was 100% after a median follow-up of 27.0 (interquartile range 13.3-58.0) months. PRT resulted in considerable tumor shrinkage, leading to complete remission in five patients and bearing the potential to provide partial or complete symptom relief. Favorable dosimetric outcomes in critical brain substructures by the use of PRT translated into reduced estimated risks for neurocognitive impairment and radiation-induced CNS malignancies compared to XRT. CONCLUSIONS: PRT is an effective treatment option for advanced JNA with minimal acute morbidity and the potential for reduced radiation-induced long-term complications.

Carbon-ion radiotherapy (CIRT) as treatment of pancreatic cancer at HIT: initial radiation plan analysis of the prospective phase II PACK-study.

PURPOSE: To analyze the dose objectives and constraints applied at the prospective phase II PACK-study at Heidelberg ion therapy center (HIT) for different radiobiological models. METHODS: Treatment plans of 14 patients from the PACK-study were analyzed and recomputed in terms of physical, biological dose and dose-averaged linear energy transfer (LETd). Both LEM-I (local effect model 1) and the adapted NIRS-MKM (microdosimetric kinetic model), were used for relative biological effectiveness (RBE)-weighted dose calculations (DBio|HIT and DBio|NIRS). A new constraint to the gastrointestinal (GI) tract was derived from the National Institute of Radiological Science (NIRS) clinical experience and considered for plan reoptimization (DBio|NIRS-const_48Gy and DBio|NIRS-const_50.4Gy). The Lyman-Kutcher-Burman (LKB) model of Normal Tissue Complication Probability (NTCP) for GI toxicity endpoints was computed. Furthermore, the computed LETd distribution was evaluated and correlated with Local Control (LC). RESULTS: Only two patients showed a LETd98% in the GTV greater than 44 keV/µm. A HIT-dose constraint to the GI of [Formula: see text] was derived from the NIRS experience, in alternative to the standard at HIT Dmax = 45.6 GyRBEHIT. In comparison with the original DBio|HIT,DBio|NIRS-const_48GyandDBio|NIRS-const_50.4Gy resulted in an increase in the ITV's D98% of 8.7% and 11.3%. The NTCP calculation resulted in a probability for gastrointestinal bleeding of 4.5%, 12.3% and 13.0%, for DBio|NIRS, DBio|NIRS-const_48Gy and DBio|NIRS-const_50.4Gy, respectively. CONCLUSION: The results indicate that the current standards applied at HIT for CIRT closely align with the Japanese experience. However, to enhance tumor coverage, a more relaxed constraint on the GI tract may be considered. As the PACK-trial progresses, further analyses of various clinical endpoints are anticipated.

Dual-energy CT-based stopping power prediction for dental materials in particle therapy.

Radiotherapy with protons or light ions can offer accurate and precise treatment delivery. Accurate knowledge of the stopping power ratio (SPR) distribution of the tissues in the patient is crucial for improving dose prediction in patients during planning. However, materials of uncertain stoichiometric composition such as dental implant and restoration materials can substantially impair particle therapy treatment planning due to related SPR prediction uncertainties. This study investigated the impact of using dual-energy computed tomography (DECT) imaging for characterizing and compensating for commonly used dental implant and restoration materials during particle therapy treatment planning. Radiological material parameters of ten common dental materials were determined using two different DECT techniques: sequential acquisition CT (SACT) and dual-layer spectral CT (DLCT). DECT-based direct SPR predictions of dental materials via spectral image data were compared to conventional single-energy CT (SECT)-based SPR predictions obtained via indirect CT-number-to-SPR conversion. DECT techniques were found overall to reduce uncertainty in SPR predictions in dental implant and restoration materials compared to SECT, although DECT methods showed limitations for materials containing elements of a high atomic number. To assess the influence on treatment planning, an anthropomorphic head phantom with a removable tooth containing lithium disilicate as a dental material was used. The results indicated that both DECT techniques predicted similar ranges for beams unobstructed by dental material in the head phantom. When ion beams passed through the lithium disilicate restoration, DLCT-based SPR predictions using a projection-based method showed better agreement with measured reference SPR values (range deviation: 0.2 mm) compared to SECT-based predictions. DECT-based SPR prediction may improve the management of certain non-tissue dental implant and restoration materials and subsequently increase dose prediction accuracy.

Do We Preserve Tumor Control Probability (TCP) in FLASH Radiotherapy? A Model-Based Analysis.

Reports of concurrent sparing of normal tissue and iso-effective treatment of tumors at ultra-high dose-rates (uHDR) have fueled the growing field of FLASH radiotherapy. However, iso-effectiveness in tumors is often deduced from the absence of a significant difference in their growth kinetics. In a model-based analysis, we investigate the meaningfulness of these indications for the clinical treatment outcome. The predictions of a previously benchmarked model of uHDR sparing in the "UNIfied and VERSatile bio response Engine" (UNIVERSE) are combined with existing models of tumor volume kinetics as well as tumor control probability (TCP) and compared to experimental data. The potential TCP of FLASH radiotherapy is investigated by varying the assumed dose-rate, fractionation schemes and oxygen concentration in the target. The developed framework describes the reported tumor growth kinetics appropriately, indicating that sparing effects could be present in the tumor but might be too small to be detected with the number of animals used. The TCP predictions show the possibility of substantial loss of treatment efficacy for FLASH radiotherapy depending on several variables, including the fractionation scheme, oxygen level, and DNA repair kinetics. The possible loss of TCP should be seriously considered when assessing the clinical viability of FLASH treatments.

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.

Radiotherapy with Helium Ions Has the Potential to Improve Both Endocrine and Neurocognitive Outcome in Pediatric Patients with Ependymoma.

Ependymomas are the third most-frequent pediatric brain tumors. To prevent local recurrence, the resection site should be irradiated. Compared to photon radiation treatment, proton therapy often achieves even better results regarding target coverage and organ-sparing. Due to their physical properties, helium ions could further reduce side effects, providing better protection of healthy tissue despite similar target coverage. In our in silico study, 15 pediatric ependymoma patients were considered. All patients underwent adjuvant radiotherapeutic treatment with active-scanned protons at Heidelberg Ion Beam Therapy Center (HIT). Both helium ion and highly conformal IMRT plans were calculated to evaluate the potential dosimetric advantage of ion beam therapy compared to the current state-of-the-art photon-based treatments. To estimate the potential clinical benefit of helium ions, normal tissue complication probabilities (NTCP) were calculated. Target coverage was comparable in all three modalities. As expected, the integral dose absorbed by healthy brain tissue could be significantly reduced with protons by up to -48% vs. IMRT. Even compared to actively scanned protons, relative dose reductions for critical neuronal structures of up to another -39% were achieved when using helium ions. The dose distribution of helium ions is significantly superior when compared to proton therapy and IMRT due to the improved sparing of OAR. In fact, previous studies could clearly demonstrate that the dosimetric advantage of protons translates into a measurable clinical benefit for pediatric patients with brain tumors. Given the dose-response relationship of critical organs at risk combined with NTCP calculation, the results of our study provide a strong rationale that the use of helium ions has the potential to even further reduce the risk for treatment related sequelae.

Ways to unravel the clinical potential of carbon ions for head and neck cancer reirradiation: dosimetric comparison and local failure pattern analysis as part of the prospective randomized CARE trial.

BACKGROUND: Carbon ion radiotherapy (CIRT) yields biophysical advantages compared to photons but randomized studies for the reirradiation setting are pending. The aim of the current project was to evaluate potential clinical benefits and drawbacks of CIRT compared to volumetric modulated arc therapy (VMAT) in recurrent head and neck cancer. METHODS: Dose-volume parameters and local failure patterns of CIRT compared to VMAT were evaluate in 16 patients from the randomized CARE trial on head and neck cancer reirradiation. RESULTS: Despite an increased target dose, CIRT resulted in significantly reduced organ at risk (OAR) dose across all patients (- 8.7% Dmean). The dose-volume benefits were most pronounced in the brainstem (- 20.7% Dmax) and the optic chiasma (- 13.0% Dmax). The most frequent local failure was type E (extraneous; 50%), followed type B (peripheral; 33%) and type A (central; 17%). In one patient with type A biological and/or dosimetric failure after CIRT, mMKM dose recalculation revealed reduced target coverage. CONCLUSIONS: CIRT resulted in highly improved critical OAR sparing compared to VMAT across all head and neck cancer reirradiation scenarios despite an increased prescription dose. Local failure pattern analysis revealed further potential CIRT specific clinical benefits and potential pitfalls with regard to image-guidance and biological dose-optimization.

Impact of DNA Repair Kinetics and Dose Rate on RBE Predictions in the UNIVERSE.

Accurate knowledge of the relative biological effectiveness (RBE) and its dependencies is crucial to support modern ion beam therapy and its further development. However, the influence of different dose rates of the reference radiation and ion beam are rarely considered. The ion beam RBE-model within our "UNIfied and VERSatile bio response Engine" (UNIVERSE) is extended by including DNA damage repair kinetics to investigate the impact of dose-rate effects on the predicted RBE. It was found that dose-rate effects increase with dose and biological effects saturate at high dose-rates, which is consistent with data- and model-based studies in the literature. In a comparison with RBE measurements from a high dose in-vivo study, the predictions of the presented modification were found to be improved in comparison to the previous version of UNIVERSE and existing clinical approaches that disregard dose-rate effects. Consequently, DNA repair kinetics and the different dose rates applied by the reference and ion beams might need to be considered in biophysical models to accurately predict the RBE. Additionally, this study marks an important step in the further development of UNIVERSE, extending its capabilities in giving theoretical guidance to support progress in ion beam 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.

Potential of a Second-Generation Dual-Layer Spectral CT for Dose Calculation in Particle Therapy Treatment Planning.

In particle therapy treatment planning, dose calculation is conducted using patient-specific maps of tissue ion stopping power ratio (SPR) to predict beam ranges. Improving patient-specific SPR prediction is therefore essential for accurate dose calculation. In this study, we investigated the use of the Spectral CT 7500, a second-generation dual-layer spectral computed tomography (DLCT) system, as an alternative to conventional single-energy CT (SECT) for patient-specific SPR prediction. This dual-energy CT (DECT)-based method allows for the direct prediction of SPR from quantitative measurements of relative electron density and effective atomic number using the Bethe equation, whereas the conventional SECT-based method consists of indirect image data-based prediction through the conversion of calibrated CT numbers to SPR. The performance of the Spectral CT 7500 in particle therapy treatment planning was characterized by conducting a thorough analysis of its SPR prediction accuracy for both tissue-equivalent materials and common non-tissue implant materials. In both instances, DLCT was found to reduce uncertainty in SPR predictions compared to SECT. Mean deviations of 0.7% and 1.6% from measured SPR values were found for DLCT- and SECT-based predictions, respectively, in tissue-equivalent materials. Furthermore, end-to-end analyses of DLCT-based treatment planning were performed for proton, helium, and carbon ion therapies with anthropomorphic head and pelvic phantoms. 3D gamma analysis was performed with ionization chamber array measurements as the reference. DLCT-predicted dose distributions revealed higher passing rates compared to SECT-predicted dose distributions. In the DLCT-based treatment plans, measured distal-edge evaluation layers were within 1 mm of their predicted positions, demonstrating the accuracy of DLCT-based particle range prediction. This study demonstrated that the use of the Spectral CT 7500 in particle therapy treatment planning may lead to better agreement between planned and delivered dose compared to current clinical SECT systems.

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.