Monte Carlo simulations of microdosimetry and radiolytic species production at long time post proton irradiation using GATE and Geant4-DNA.

BACKGROUND: Radiobiological effectiveness of radiation in cancer treatment can be studied at different scales (molecular till organ scale) and different time post irradiation. The production of free radicals and reactive oxygen species during water radiolysis is particularly relevant to understand the fundamental mechanisms playing a role in observed biological outcomes. The development and validation of Monte Carlo tools integrating the simulation of physical, physico-chemical and chemical stages after radiation is very important to maintain with experiments. PURPOSE: Therefore, in this study, we propose to validate a new Geant4-DNA chemistry module through the simulation of water radiolysis and Fricke dosimetry experiments on a proton preclinical beam line. MATERIAL AND METHODS: In this study, we used the GATE Monte Carlo simulation platform (version 9.3) to simulate a 67.5 MeV proton beam produced with the ARRONAX isochronous cyclotron (IBA Cyclone 70XP) at conventional dose rate (0.2 Gy/s) to simulate the irradiation of ultra-pure liquid water samples and Fricke dosimeter. We compared the depth dose profile with measurements performed with a plane parallel Advanced PTW 34045 Markus ionization chamber. Then, a new Geant4-DNA chemistry application proposed from Geant4 version 11.2 has been used to assess the evolution of HO • ${\mathrm{HO}}^ \bullet $ , e aq - ${\mathrm{e}}_{{\mathrm{aq}}}^ - $ , H 3 O + ${{\mathrm{H}}}_3{{\mathrm{O}}}^ + $ , H 2 O 2 ${{\mathrm{H}}}_2{{\mathrm{O}}}_2$ , H 2 ${{\mathrm{H}}}_2$ , HO 2 • ${\mathrm{HO}}_2^ \bullet $ , HO 2 - , O 2 • - ${\mathrm{HO}}_2^ - ,{\mathrm{\ O}}_2^{ \bullet - }$ and HO - ${\mathrm{HO}}^ - $ reactive species along time until 1-h post-irradiation. In particular, the effect of oxygen and pH has been investigated through comparisons with experimental measurements of radiolytic yields for H 2 O 2 ${{\mathrm{H}}}_2{{\mathrm{O}}}_2$ and Fe3+. RESULTS: GATE simulations reproduced, within 4%, the depth dose profile in liquid water. With Geant4-DNA, we were able to reproduce experimental H 2 O 2 ${{\mathrm{H}}}_2{{\mathrm{O}}}_2$ radiolytic yields 1-h post-irradiation in aerated and deaerated conditions, showing the impact of small changes in oxygen concentrations on species evolution along time. For the Fricke dosimeter, simulated G(Fe3+) is 15.97 ± 0.2 molecules/100 eV which is 11% higher than the measured value (14.4 ± 04 molecules/100 eV). CONCLUSIONS: These results aim to be consolidated by new comparisons involving other radiolytic species, such as e aq - ${\mathrm{e}}_{{\mathrm{aq}}}^ - $ or , O 2 • - $,{\mathrm{\ O}}_2^{ \bullet - }$ to further study the mechanisms underlying the FLASH effect observed at ultra-high dose rates (UHDR).

Methodology for small animals targeted irradiations at conventional and ultra-high dose rates 65 MeV proton beam.

As part of translational research projects, mice may be irradiated on radiobiology platforms such as the one at the ARRONAX cyclotron. Generally, these platforms do not feature an integrated imaging system. Moreover, in the context of ultra-high dose-rate radiotherapy (FLASH-RT), treatment planning should consider potential changes in the beam characteristics and internal movements in the animal. A patient-like set-up and methodology has been implemented to ensure target coverage during conformal irradiations of the brain, lungs and intestines. In addition, respiratory cycle amplitudes were quantified by fluoroscopic acquisitions on a mouse, to ensure organ coverage and to assess the impact of respiration during FLASH-RT using the 4D digital phantom MOBY. Furthermore, beam incidence direction was studied from mice µCBCT and Monte Carlo simulations. Finally,in vivodosimetry with dose-rate independent radiochromic films (OC-1) and their LET dependency were investigated. The immobilization system ensures that the animal is held in a safe and suitable position. The geometrical evaluation of organ coverage, after the addition of the margins around the organs, was satisfactory. Moreover, no measured differences were found between CONV and FLASH beams enabling a single model of the beamline for all planning studies. Finally, the LET-dependency of the OC-1 film was determined and experimentally verified with phantoms, as well as the feasibility of using these filmsin vivoto validate the targeting. The methodology developed ensures accurate and reproducible preclinical irradiations in CONV and FLASH-RT without in-room image guidance in terms of positioning, dose calculation andin vivodosimetry.

First evidence of in vivo effect of FLASH radiotherapy with helium ions in zebrafish embryos.

The ability to reduce toxicity of ultra-high dose rate (UHDR) helium ion irradiation has not been reported in vivo. Here, we tested UHDR helium ion irradiation in an embryonic zebrafish model. Our results show that UHDR helium ions spare body development and reduce spine curvature, compared to conventional dose rate.

Ultrahigh-Dose-Rate Proton Irradiation Elicits Reduced Toxicity in Zebrafish Embryos.

Purpose: Recently, ultrahigh-dose-rate radiation therapy (UHDR-RT) has emerged as a promising strategy to increase the benefit/risk ratio of external RT. Extensive work is on the way to characterize the physical and biological parameters that control the so-called "Flash" effect. However, this healthy/tumor differential effect is observable in in vivo models, which thereby drastically limits the amount of work that is achievable in a timely manner. Methods and Materials: In this study, zebrafish embryos were used to compare the effect of UHDR irradiation (8-9 kGy/s) to conventional RT dose rate (0.2 Gy/s) with a 68 MeV proton beam. Viability, body length, spine curvature, and pericardial edema were measured 4 days postirradiation. Results: We show that body length is significantly greater after UHDR-RT compared with conventional RT by 180 µm at 30 Gy and 90 µm at 40 Gy, while pericardial edema is only reduced at 30 Gy. No differences were obtained in terms of survival or spine curvature. Conclusions: Zebrafish embryo length appears as a robust endpoint, and we anticipate that this model will substantially fasten the study of UHDR proton-beam parameters necessary for "Flash."

Proton Irradiations at Ultra-High Dose Rate vs. Conventional Dose Rate: Strong Impact on Hydrogen Peroxide Yield.

During ultra-high dose rate (UHDR) external radiation therapy, healthy tissues appear to be spared while tumor control remains the same compared to conventional dose rate. However, the understanding of radiochemical and biological mechanisms involved are still to be discussed. This study shows how the hydrogen peroxide (H2O2) production, one of the reactive oxygen species (ROS), could be controlled by early heterogenous radiolysis processes in water during UHDR proton-beam irradiations. Pure water was irradiated in the plateau region (track-segment) with 68 MeV protons under conventional (0.2 Gy/s) and several UHDR conditions (40 Gy/s to 60 kGy/s) at the ARRONAX cyclotron. Production of H2O2 was then monitored using the Ghormley triiodide method. New values of GTS(H2O2) were added in conventional dose rate. A substantial decrease in H2O2 production was observed from 0.2 to 1.5 kGy/s with a more dramatic decrease below 100 Gy/ s. At higher dose rate, up to 60 kGy/s, the H2O2 production stayed stable with a mean decrease of 38% ± 4%. This finding, associated to the decrease in the production of hydroxyl radical (•OH) already observed in other studies in similar conditions can be explained by the well-known spur theory in radiation chemistry. Thus, a two-step FLASH-RT mechanism can be envisioned: an early step at the microsecond scale mainly controlled by heterogenous radiolysis, and a second, slower, dominated by O2 depletion and biochemical processes. To validate this hypothesis, more measurements of radiolytic species will soon be performed, including radicals and associated lifetimes.

Technical note: Proton beam dosimetry at ultra-high dose rates (FLASH): Evaluation of GAFchromic™ (EBT3, EBT-XD) and OrthoChromic (OC-1) film performances.

PURPOSE: The ARRONAX cyclotron facility offers the possibility to deliver proton beams from low to ultra-high dose rates (UHDR). As a good control of the dosimetry is a prerequisite of UHDR experimentations, we evaluated in different conditions the usability and the dose rate dependency of several radiochromic films commonly used for dosimetry in radiotherapy. METHODS: We compared the dose rate dependency of three types of radiochromic films: GAFchromic™ EBT3 and GAFchromic™ EBT-XD (Ashland Inc., Wayne, NJ, USA), and OrthoChromic OC-1 (OrthoChrome Inc., Hillsborough, NJ, USA), after proton irradiations at various mean dose rates (0.25, 40, 1500, and 7500 Gy/s) and for 10 doses (2-130 Gy). We also evaluated the dose rate dependency of each film considering beam structures, from single pulse to multiple pulses with various frequencies. RESULTS: EBT3 and EBT-XD films showed differences of response between conventional (0.25 Gy/s) and UHDR (7500 Gy/s) conditions, above 10 Gy. On the contrary, OC-1 films did not present overall difference of response for doses except below 3 Gy. We observed an increase of the netOD with the mean dose rate for EBT3 and EBT-XD films. OC-1 films did not show any impact of the mean dose rate up to 7500 Gy/s, above 3 Gy. No difference was found based on the beam structure, for all three types of films. CONCLUSIONS: EBT3 and EBT-XD radiochromic films should be used with caution for the dosimetry of UHDR proton beams over 10 Gy. Their overresponse, which increases with mean dose rate and dose, could lead to non-negligible overestimations of the absolute dose. OC-1 films are dose rate independent up to 7500 Gy/s in proton beams. Films response is not impacted by the beam structure. A broader investigation of the usability of OC-1 films in UHDR conditions should be conducted at intermediate and higher mean dose rates and other beam energies.

A Monte Carlo Determination of Dose and Range Uncertainties for Preclinical Studies with a Proton Beam.

Proton therapy (PRT) is an irradiation technique that aims at limiting normal tissue damage while maintaining the tumor response. To study its specificities, the ARRONAX cyclotron is currently developing a preclinical structure compatible with biological experiments. A prerequisite is to identify and control uncertainties on the ARRONAX beamline, which can lead to significant biases in the observed biological results and dose-response relationships, as for any facility. This paper summarizes and quantifies the impact of uncertainty on proton range, absorbed dose, and dose homogeneity in a preclinical context of cell or small animal irradiation on the Bragg curve, using Monte Carlo simulations. All possible sources of uncertainty were investigated and discussed independently. Those with a significant impact were identified, and protocols were established to reduce their consequences. Overall, the uncertainties evaluated were similar to those from clinical practice and are considered compatible with the performance of radiobiological experiments, as well as the study of dose-response relationships on this proton beam. Another conclusion of this study is that Monte Carlo simulations can be used to help build preclinical lines in other setups.