Modeling of cellular response after FLASH irradiation: a quantitative analysis based on the radiolytic oxygen depletion hypothesis

Possible mechanisms and simulation modeling of FLASH radiotherapy

FLASH radiotherapy (FLASH-RT) has great potential to improve patient outcomes. It delivers radiation doses at an ultra-high dose rate (UHDR: ≥ 40 Gy/s) in a single instant or a few pulses. Much higher irradiation doses can be administered to tumors with FLASH-RT than with conventional dose rate (0.01-0.40 Gy/s) radiotherapy. UHDR irradiation can suppress toxicity in normal tissues while sustaining antitumor efficiency, which is referred to as the FLASH effect. However, the mechanisms underlying the effects of the FLASH remain unclear. To clarify these mechanisms, the development of simulation models that can contribute to treatment planning for FLASH-RT is still underway. Previous studies indicated that transient oxygen depletion or augmented reactions between secondary reactive species produced by irradiation may be involved in this process. To discuss the possible mechanisms of the FLASH effect and its clinical potential, we summarized the physicochemical, chemical, and biological perspectives as well as the development of simulation modeling for FLASH-RT.

The dose-related plateau effect of surviving fraction in normal tissue during the ultra-high-dose-rate radiotherapy

Objective.Ultra-high-dose-rate radiotherapy, referred to as FLASH therapy, has been demonstrated to reduce the damage of normal tissue as well as inhibiting tumor growth compared with conventional dose-rate radiotherapy. The transient hypoxia may be a vital explanation for sparing the normal tissue. The heterogeneity of oxygen distribution for different doses and dose rates in the different radiotherapy schemes are analyzed. With these results, the influence of doses and dose rates on cell survival are evaluated in this work.Approach.The two-dimensional reaction-diffusion equations are used to describe the heterogeneity of the oxygen distribution in capillaries and tissue. A modified linear quadratic model is employed to characterize the surviving fraction at different doses and dose rates.Main results.The reduction of the damage to the normal tissue can be observed if the doses exceeds a minimum dose threshold under the ultra-high-dose-rate radiation. Also, the surviving fraction exhibits the 'plateau effect' under the ultra-high dose rates radiation, which signifies that within a specific range of doses, the surviving fraction either exhibits minimal variation or increases with the dose. For a given dose, the surviving fraction increases with the dose rate until tending to a stable value, which means that the protection in normal tissue reaches saturation.Significance.The emergence of the 'plateau effect' allows delivering the higher doses while minimizing damage to normal tissue. It is necessary to develop appropriate program of doses and dose rates for different irradiated tissue to achieve more efficient protection.

TRAX-CHEMxt: Towards the Homogeneous Chemical Stage of Radiation Damage

The indirect effect of radiation plays an important role in radio-induced biological damages. Monte Carlo codes have been widely used in recent years to study the chemical evolution of particle tracks. However, due to the large computational efforts required, their applicability is typically limited to simulations in pure water targets and to temporal scales up to the µs. In this work, a new extension of TRAX-CHEM is presented, namely TRAX-CHEMxt, able to predict the chemical yields at longer times, with the capability of exploring the homogeneous biochemical stage. Based on the species coordinates produced around one track, the set of reaction-diffusion equations is solved numerically with a computationally light approach based on concentration distributions. In the overlapping time scale (500 ns-1 µs), a very good agreement to standard TRAX-CHEM is found, with deviations below 6% for different beam qualities and oxygenations. Moreover, an improvement in the computational speed by more than three orders of magnitude is achieved. The results of this work are also compared with those from another Monte Carlo-based algorithm and a fully homogeneous code (Kinetiscope). TRAX-CHEMxt will allow for studying the variation in chemical endpoints at longer timescales with the introduction, as the next step, of biomolecules, for more realistic assessments of biological response under different radiation and environmental conditions.

Comet Assay Profiling of FLASH-Induced Damage: Mechanistic Insights into the Effects of FLASH Irradiation

Numerous studies have demonstrated the normal tissue-sparing effects of ultra-high dose rate 'FLASH' irradiation in vivo, with an associated reduction in damage burden being reported in vitro. Towards this, two key radiochemical mechanisms have been proposed: radical-radical recombination (RRR) and transient oxygen depletion (TOD), with both being proposed to lead to reduced levels of induced damage. Previously, we reported that FLASH induces lower levels of DNA strand break damage in whole-blood peripheral blood lymphocytes (WB-PBL) ex vivo, but our study failed to distinguish the mechanism(s) involved. A potential outcome of RRR is the formation of crosslink damage (particularly, if any organic radicals recombine), whilst a possible outcome of TOD is a more anoxic profile of induced damage resulting from FLASH. Therefore, the aim of the current study was to profile FLASH-induced damage via the Comet assay, assessing any DNA crosslink formation as a putative marker of RRR and/or anoxic DNA damage formation as an indicative marker of TOD, to determine the extent to which either mechanism contributes to the "FLASH effect". Following FLASH irradiation, we see no evidence of any crosslink formation; however, FLASH irradiation induces a more anoxic profile of induced damage, supporting the TOD mechanism. Furthermore, treatment of WB-PBLs pre-irradiation with BSO abrogates the reduced strand break damage burden mediated by FLASH exposures. In summary, we do not see any experimental evidence to support the RRR mechanism contributing to the reduced damage burden induced by FLASH. However, the observation of a greater anoxic profile of damage following FLASH irradiation, together with the BSO abrogation of the reduced strand break damage burden mediated by FLASH, lends further support to TOD being a driver of the reduced damage burden plus a change in the damage profile mediated by FLASH.

Proton FLASH effects on mouse skin at different oxygen tensions

Objective. Irradiation at FLASH dose rates (>40 Gy s-1) has received great attention due to its reported normal tissue sparing effect. The FLASH effect was originally observed in electron irradiations but has since been shown to also occur with both photon and proton beams. Several mechanisms have been proposed to explain the tissue sparing at high dose rates, including effects involving oxygen, such as depletion of oxygen within the irradiated cells. In this study, we investigated the protective role of FLASH proton irradiation on the skin when varying the oxygen concentration.Approach. Our double scattering proton system provided a 1.2 × 1.6 cm2elliptical field at a dose rate of ∼130 Gy s-1. The conventional dose rate was ∼0.4 Gy s-1. The legs of the FVB/N mice were marked with two tattooed dots and fixed in a holder for exposure. To alter the skin oxygen concentration, the mice were breathing pure oxygen or had their legs tied to restrict blood flow. The distance between the two dots was measured to analyze skin contraction over time.Main results. FLASH irradiation mitigated skin contraction by 15% compared to conventional dose rate irradiation. The epidermis thickness and collagen deposition at 75 d following 25 to 30 Gy exposure suggested a long-term protective function in the skin from FLASH irradiation. Providing the mice with oxygen or reducing the skin oxygen concentration removed the dose-rate-dependent difference in response.Significance. FLASH proton irradiation decreased skin contraction, epidermis thickness and collagen deposition compared to standard dose rate irradiations. The observed oxygen-dependence of the FLASH effect is consistent with, but not conclusive of, fast oxygen depletion during the exposure.

CPU-GPU coupling independent reaction times method in NASIC and application in water radiolysis by FLASH irradiation

The mechanism of the FLASH effect remains unclear and could be revealed by studying chemical reactions during irradiation. Monte Carlo simulation of the radiolytic species is an effective tool to analyze chemical reactions, but the simulation is limited by computing costs of the step-by-step simulation of radiolytic species, especially when considering beam with complex time structure. The complexity of the time structure of beams from accelerators in FLASH radiotherapy requires a high-performance Monte Carlo code. In this work, we develop a CPU-GPU coupling accelerating code with the independent reaction times (IRT) method to extend the chemical module of our nanodosimetry Monte Carlo code NASIC. Every chemical molecule in the microenvironment contains time information to consider the reactions from different tracks and simulate beams with complex time structures. Performance test shows that our code significantly improved the computing efficiency of the chemical module by four orders of magnitude. Then the code is used to study the oxygen depletion hypothesis in FLASH radiotherapy for different conditions by setting different parameters. The transient oxygen consumption rate values in the water are calculated when the pulses width ranges from 2 ps to 2 μs, the total dose ranges from 0.5 Gy to 100 Gy and the initial oxygen concentration ranges from 0.1% to 21%. The time evolution curves are simulated to study the effect of the time structure of an electron linear accelerator. Results show that the total dose in several microseconds is a better indicator reflecting the radiolytic oxygen consumption rate than the dose rate. The initial oxygen greatly affects the oxygen consumption rate because of the reaction competition. The diffusion of oxygen determined by the physiological parameters is the key factor affecting oxygen depletion during the radiation using electron linear accelerators. Our code provides an efficient tool for simulating water radiolysis in different conditions