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.

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.

Spectroscopic study of prompt-gamma emission for range verification in proton therapy.

We present the results of an investigation of the prompt-gamma emission from an interaction of a proton beam with phantom materials. Measurements were conducted with a novel setup allowing the precise selection of the investigated depth in the phantom, featuring three different materials composed of carbon, oxygen and hydrogen. We studied two beam energies of 70.54 and 130.87MeV and two detection angles: 90° and 120°. The results are presented in form of profiles of the prompt-gamma yield as a function of depth. In the analysis we focused on the transitions with the largest cross sections: 12C4.44→g.s. and 16O6.13→g.s.. We compare the profiles obtained under various irradiation conditions, with emphasis on the shape of the distal fall-off. The results are also compared to calculations including different cross-section models. They are in agreement with the model exploiting published cross-section data, but the comparison with the Talys model shows discrepancies.