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

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 ofD 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 aD 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.

Optimizing the traversal time for gantry trajectories for proton arc therapy treatment plans

Background.Proton arc therapy (PAT) is an emerging radiation therapy technique where either the gantry or the patient continuously rotates during the irradiation treatment. One of the perceived advantages of PAT is the reduced treatment time, but it is still unclear exactly how long these treatment times will be, given that no machine capable of its delivery is available on the market at the time of writing.Objective.We introduce the algorithm arc trajectory optimization method (ATOM), which aims to determine an efficient velocity profile for the gantry for rapid delivery of a given proton arc treatment plan. This algorithm could be used to minimize the delivery time of a proton arc plan without changing the plan or updating the machine.Approach.ATOM computes the trajectory with the shortest delivery time while ensuring there is enough time to deliver all spots in each energy layer and switch energy between layers. The feasibility of the dynamic gantry movement was assured by enforcing maximum and minimum limits for velocity, acceleration, and jerk. This was achieved by discretizing the gantry velocity and combining theA* algorithm with the open-source motion generation library Ruckig. The algorithm was tested on a synthetic data set as well as a liver case, a prostate case and a head and neck case.Main results.Arc trajectories for plans with 360 energy layers were calculated in under a second using 256 discrete velocities. The delivery time of the liver case, the prostate case and the head and neck case were 284 s, 288 s and 309 s respectively, for 180 energy layers.Significance.ATOM is an open-source C++ library with a Python interface that rapidly generates velocity profiles, making it a highly efficient tool for determining proton arc delivery times, which could be integrated into the treatment planning process.

High-LET charged particles: radiobiology and application for new approaches in radiotherapy

The number of patients treated with charged-particle radiotherapy as well as the number of treatment centers is increasing worldwide, particularly regarding protons. However, high-linear energy transfer (LET) particles, mainly carbon ions, are of special interest for application in radiotherapy, as their special physical features result in high precision and hence lower toxicity, and at the same time in increased efficiency in cell inactivation in the target region, i.e., the tumor. The radiobiology of high-LET particles differs with respect to DNA damage repair, cytogenetic damage, and cell death type, and their increased LET can tackle cells' resistance to hypoxia. Recent developments and perspectives, e.g., the return of high-LET particle therapy to the US with a center planned at Mayo clinics, the application of carbon ion radiotherapy using cost-reducing cyclotrons and the application of helium is foreseen to increase the interest in this type of radiotherapy. However, further preclinical research is needed to better understand the differential radiobiological mechanisms as opposed to photon radiotherapy, which will help to guide future clinical studies for optimal exploitation of high-LET particle therapy, in particular related to new concepts and innovative approaches. Herein, we summarize the basics and recent progress in high-LET particle radiobiology with a focus on carbon ions and discuss the implications of current knowledge for charged-particle radiotherapy. We emphasize the potential of high-LET particles with respect to immunogenicity and especially their combination with immunotherapy.

Partitioning of discrete proton arcs into interlaced subplans can bring proton arc advances to existing proton facilities

BACKGROUND

Proton arcs have shown potential to reduce the dose to organs at risks (OARs) by delivering the protons from many different directions. While most previous studies have been focused on dynamic arcs (delivery during rotation), an alternative approach is discrete arcs, where step-and-shoot delivery is used over a large number of beam directions. The major advantage of discrete arcs is that they can be delivered at existing proton facilities. However, this advantage comes at the expense of longer treatment times.

PURPOSE

To exploit the dosimetric advantages of proton arcs, while achieving reasonable delivery times, we propose a partitioning approach where discrete arc plans are split into subplans to be delivered over different fractions in the treatment course.

METHODS

For three oropharyngeal cancer patients, four different arc plans have been created and compared to the corresponding clinical IMPT plan. The treatment plans are all planned to be delivered in 35 fractions, but with different delivery approaches over the fractions. The first arc plan (1×30) has 30 directions to be delivered every fraction, while the others are partitioned into subplans with 10 and 6 beam directions, each to be delivered every third (3×10), fifth fraction (5×6), or seventh fraction (7×10). All plans are assessed with respect to delivery time, target robustness over the treatment course, doses to OARs and NTCP for dysphagia and xerostomia.

RESULTS

The delivery time (including an additional delay of 30 s between the discrete directions to simulate manual interaction with the treatment control system) is reduced from on average 25.2 min for the 1×30 plan to 9.2 min for the 3×10 and 7×10 plans and 5.7 min for the 5×6 plans. The delivery time for the IMPT plan is 7.9 min. When accounting for the combination of delivery time, target robustness, OAR sparing, and NTCP reduction, the plans with 10 directions in each fraction are the preferred choice. Both the 3×10 and 7×10 plans show improved target robustness compared to the 1×30 plans, while keeping OAR doses and NTCP values at almost as low levels as for the 1×30 plans. For all patients the NTCP values for dysphagia are lower for the partitioned plans with 10 directions compared to the IMPT plans. NTCP reduction for xerostomia compared to IMPT is seen in two of the three patients. The best results are seen for the first patient, where the NTCP reductions for the 7×10 plan are 1.6 p.p. (grade 2 xerostomia) and 1.5 p.p. (grade 2 dysphagia). The corresponding NTCP reductions for the 1×30 plan are 2.7 p.p. (xerostomia, grade 2) and 2.0 p.p. (dysphagia, grade 2).

CONCLUSIONS

Discrete proton arcs can be implemented at any proton facility with reasonable treatment times using a partitioning approach. The technique also makes the proton arc treatments more robust to changes in the patient anatomy.

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.

Roadmap: helium ion therapy

Helium ion beam therapy for the treatment of cancer was one of several developed and studied particle treatments in the 1950s, leading to clinical trials beginning in 1975 at the Lawrence Berkeley National Laboratory. The trial shutdown was followed by decades of research and clinical silence on the topic while proton and carbon ion therapy made debuts at research facilities and academic hospitals worldwide. The lack of progression in understanding the principle facets of helium ion beam therapy in terms of physics, biological and clinical findings persists today, mainly attributable to its highly limited availability. Despite this major setback, there is an increasing focus on evaluating and establishing clinical and research programs using helium ion beams, with both therapy and imaging initiatives to supplement the clinical palette of radiotherapy in the treatment of aggressive disease and sensitive clinical cases. Moreover, due its intermediate physical and radio-biological properties between proton and carbon ion beams, helium ions may provide a streamlined economic steppingstone towards an era of widespread use of different particle species in light and heavy ion therapy. With respect to the clinical proton beams, helium ions exhibit superior physical properties such as reduced lateral scattering and range straggling with higher relative biological effectiveness (RBE) and dose-weighted linear energy transfer (LET_d) ranging from ∼4 keVμm^-1to ∼40 keVμm^-1. In the frame of heavy ion therapy using carbon, oxygen or neon ions, where LET_dincreases beyond 100 keVμm^-1, helium ions exhibit similar physical attributes such as a sharp lateral penumbra, however, with reduced radio-biological uncertainties and without potentially spoiling dose distributions due to excess fragmentation of heavier ion beams, particularly for higher penetration depths. This roadmap presents an overview of the current state-of-the-art and future directions of helium ion therapy: understanding physics and improving modeling, understanding biology and improving modeling, imaging techniques using helium ions and refining and establishing clinical approaches and aims from learned experience with protons. These topics are organized and presented into three main sections, outlining current and future tasks in establishing clinical and research programs using helium ion beams-A. Physics B. Biological and C. Clinical Perspectives.

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

Treatment plan planning and optimization procedures were developed for spot-scanning hadron arc (SHArc) delivery using the RayStation TPS and a GPU-accelerated dose engine (^†TPS-XXX). The SHArc optimization algorithm is applicable for charged particle beams and determines angle-dependencies for spot/energy selection with three main initiatives: i) achieve standard clinical optimization goals and constraints for target and OARs, ii) target dose robustness and iii) increasing LET in the target volume. Three patient cases previously treated at the ^†INSTITUTION-XXX were selected for evaluation of conventional versus arc delivery for the two clinical particle beams (helium [⁴He] and carbon [¹²C] ions): glioblastoma, prostate-adenocarcinoma and skull-base chordoma. Biological dose and dose-averaged linear energy transfer (LET_d) distributions for SHArc were evaluated against conventional planning techniques (VMAT and IMPT_2F) applying the modified microdosimetric kinetic model (mMKM) for considering bio-effect with (α/β)_x=2Gy. 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 >10GyRBE compared to conventional 2F planning. For carbon ions, LET_d was increased in the target volume from ∼40-60keV/µm to ∼80-140keV/µm for SHArc compared to conventional treatments. Favorable LET_d distributions were possible with the SHArc approach, with maximum LET_d in CTV/GTV and potential reductions of high-LET regions in normal tissues and OARs. Compared to VMAT, SHArc affords substantial reductions in normal tissue dose (40-70%).

CONCLUSION

SHArc therapy offers potential treatment benefits such as increased normal tissue sparing from higher doses >10GyRBE, enhanced target LET_d, and potential reduction in high-LET components in OARs. Findings justify further development of robust SHArc treatment planning towards potential clinical translation.