May oxygen depletion explain the FLASH effect? A chemical track structure analysis

Monte Carlo in the mechanistic modelling of the FLASH effect: a review

FLASH radiotherapy employs ultra-high dose rates of>40Gy s-1, which may reduce normal tissue complication as compared to conventional dose rate treatments, while still ensuring the same level of tumour control. The potential benefit this can offer to patients has been the cause of great interest within the radiation oncology community, but this has not translated to a direct understanding of the FLASH effect. The oxygen depletion and inter-track interaction hypotheses are currently the leading explanations as to the mechanisms behind FLASH, but these are still not well understood, with many questions remaining about the exact underpinnings of FLASH and the treatment parameters required to optimally induce it. Monte Carlo simulations may hold the key to unlocking the mystery behind FLASH, allowing for analysis of the underpinning mechanisms at a fundamental level, where the interactions between individual radiation particles, DNA strands and chemical species can be studied. Currently, however, there is still a great deal of disagreement in simulation findings and the importance of the different mechanisms they support. This review discusses current studies into the mechanisms of FLASH using the Monte Carlo method. The simulation parameters and results for all major investigations are provided. Discussion primarily revolves around the oxygen depletion and inter-track interactions hypotheses, though other, more novel, theories are also mentioned. A general list of recommendations for future simulations is provided, informed by the articles discussed. This review highlights some of the useful parameters and simulation methodologies that may be required to finally understand the FLASH effect.

Probing Spatiotemporal Effects of Intertrack Recombination with a New Implementation of Simultaneous Multiple Tracks in TRAX-CHEM

Among the most investigated hypotheses for a radiobiological explanation of the mechanism behind the FLASH effect in ultra-high dose rate radiotherapy, intertrack recombination between particle tracks arriving at a close spatiotemporal distance has been suggested. In the present work, we examine these conditions for different beam qualities and energies, defining the limits of both space and time where a non-negligible chemical effect is expected. To this purpose the TRAX-CHEM chemical track structure Monte Carlo code has been extended to handle several particle tracks at the same time, separated by pre-defined spatial and temporal distances. We analyzed the yields of different radicals as compared to the non-interacting track conditions and we evaluated the difference. We find a negligible role of intertrack for spatial distances larger than 1 μm, while for temporal distances up to μs, a non-negligible interaction is observed especially at higher LET. In addition, we emphasize the non-monotonic behavior of some relative yield as a function of the time separation, in particular of H2O2, due to the onset of a different reaction involving solvated electrons besides well-known OH· recombination.

Exploring the Metabolic Impact of FLASH Radiotherapy

FLASH radiotherapy (FLASH RT) is an innovative modality in cancer treatment that delivers ultrahigh dose rates (UHDRs), distinguishing it from conventional radiotherapy (CRT). FLASH RT has demonstrated the potential to enhance the therapeutic window by reducing radiation-induced damage to normal tissues while maintaining tumor control, a phenomenon termed the FLASH effect. Despite promising outcomes, the precise mechanisms underlying the FLASH effect remain elusive and are a focal point of current research. This review explores the metabolic and cellular responses to FLASH RT compared to CRT, with particular focus on the differential impacts on normal and tumor tissues. Key findings suggest that FLASH RT may mitigate damage in healthy tissues via altered reactive oxygen species (ROS) dynamics, which attenuate downstream oxidative damage. Studies indicate the FLASH RT influences iron metabolism and lipid peroxidation pathways differently than CRT. Additionally, various studies indicate that FLASH RT promotes the preservation of mitochondrial integrity and function, which helps maintain apoptotic pathways in normal tissues, attenuating damage. Current knowledge of the metabolic influences following FLASH RT highlights its potential to minimize toxicity in normal tissues, while also emphasizing the need for further studies in biologically relevant, complex systems to better understand its clinical potential. By targeting distinct metabolic pathways, FLASH RT could represent a transformative advance in RT, ultimately improving the therapeutic window for cancer treatment.

FLASH Radiotherapy: Benefits, Mechanisms, and Obstacles to Its Clinical Application

Radiotherapy (RT) has been shown to be a cornerstone of both palliative and curative tumor care. RT has generally been reported to be sharply limited by ionizing radiation (IR)-induced toxicity, thereby constraining the control effect of RT on tumor growth. FLASH-RT is the delivery of ultra-high dose rate (UHDR) several orders of magnitude higher than what is presently used in conventional RT (CONV-RT). The FLASH-RT clinical trials have been designed to examine the UHDR deliverability, the effectiveness of tumor control, the dose tolerance of normal tissue, and the reproducibility of treatment effects across several institutions. Although it is still in its infancy, FLASH-RT has been shown to have potential to rival current RT in terms of safety. Several studies have suggested that the adoption of FLASH-RT is very limited, and the incorporation of this new technique into routine clinical RT will require the use of accurate dosimetry methods and reproducible equipment that enable the reliable and robust measurements of doses and dose rates. The purpose of this review is to highlight the advantages of this technology, the potential mechanisms underpinning the FLASH-RT effect, and the major challenges that need to be tackled in the clinical transfer of FLASH-RT.

The oxygen puzzle in FLASH radiotherapy: A comprehensive review and experimental outlook

FLASH radiotherapy is attracting increasing interest because it maintains tumor control while inflicting less damage to normal tissues compared to conventional radiotherapy. This sparing effect, the so-called FLASH effect, is achieved when radiation is delivered at ultra-high dose rates (≥40 Gy/s). Although the FLASH effect has already been demonstrated in several preclinical models, a complete mechanistic description explaining why tumors and normal tissues respond differently is still missing. None of the current hypotheses fully explains the experimental evidence. A common point between many of these is the role of oxygen, which is described as a major factor, either through transient hypoxia in the form of dissolved molecules, or reactive oxygen species (ROS). Therefore, this review focuses on both forms of this molecule, retracing old and more recent theories, while proposing new mechanisms that could provide a complete description of the FLASH effect based on preclinical and experimental evidence. In addition, this manuscript describes a set of experiments designed to provide the FLASH community with new tools for exploring the post-irradiation fate of ROS and their potential biological implications.

An in-silico study of conventional and FLASH radiotherapy iso-effectiveness: potential impact of radiolytic oxygen depletion on tumor growth curves and tumor control probability

Objective. This work aims to investigate the iso-effectiveness of conventional and FLASH radiotherapy on tumors through in-silico mathematical models. We focused on the role of radiolytic oxygen depletion (ROD), which has been argued as a possible factor to explain the FLASH effect.Approach. We used a spatiotemporal reaction-diffusion model, including ROD, to simulate tumor oxygenation and response. From those oxygen distributions we obtained surviving fractions (SFs) using the linear-quadratic (LQ) model with the oxygen enhancement ratios (OERs). We then employed the calculated SFs to describe the evolution of preclinical tumor volumes through a mathematical model of tumor response, and we also extrapolated those results to calculate tumor control probabilities (TCPs) using the Poisson-LQ approach.Main results. Our study suggests that the ROD effect may cause differences in SF between FLASH and conventional radiotherapy, especially in lowα/βandpoorly oxygenatedcells. However, a statistical analysis showed that these changes in SF generally do not result in significant differences in the evolution of preclinical tumor growth curves when the sample size is small, because such differences in SF may not be noticeable in the heterogeneity of the population of animals. Nonetheless, when extrapolating this effect to TCP curves, we observed important differences between both techniques (TCP is lower in FLASH radiotherapy). When analyzing the response of tumors with heterogeneous oxygenations, differences in TCP are more important forwell oxygenatedtumors. This apparent contradiction with the results obtained for homogeneously oxygenated cells is explained by the complex interplay between the heterogeneity of tumor oxygenation, the OER effect, and the ROD effect.Significance. This study supports the experimentally observed iso-effectiveness of FLASH and conventional radiotherapy when analyzing the volume evolution of preclinical tumors (that are far from control). However, this study also hints that tumor growth curves may be less sensitive to small variations in SF than tumor control probability: ROD may lead to increased SF in FLASH radiotherapy, which while not large enough to cause significant differences in tumor growth curves, could lead to important differences in clinical TCPs. Nonetheless, it cannot be discarded that other effects not modeled in this work, like radiation-induced immune effects, can contribute to tumor control and maintain the iso-effectiveness of FLASH radiotherapy. The study of tumor growth curves may not be the ideal experiment to test the iso-effectiveness of FLASH, and experiments reporting TCP orD50may be preferred.

Current views on mechanisms of the FLASH effect in cancer radiotherapy

FLASH radiotherapy (FLASH-RT) is a new modality of radiotherapy that delivers doses with ultra-high dose rates. The FLASH effect was defined as the ability of FLASH-RT to suppress tumor growth while sparing normal tissues. Although the FLASH effect has been proven to be valid in various models by different modalities of irradiation and clinical trials of FLASH-RT have achieved promising initial success, the exact underlying mechanism is still unclear. This article summarizes mainstream hypotheses of the FLASH effect at physicochemical and biological levels, including oxygen depletion and free radical reactions, nuclear and mitochondria damage, as well as immune response. These hypotheses contribute reasonable explanations to the FLASH effect and are interconnected according to the chronological order of the organism's response to ionizing radiation. By collating the existing consensus, evidence and hypotheses, this article provides a comprehensive overview of potential mechanisms of the FLASH effect and practical guidance for future investigation in the field of FLASH-RT.

Modelling radiobiology

Radiotherapy has played an essential role in cancer treatment for over a century, and remains one of the best-studied methods of cancer treatment. Because of its close links with the physical sciences, it has been the subject of extensive quantitative mathematical modelling, but a complete understanding of the mechanisms of radiotherapy has remained elusive. In part this is because of the complexity and range of scales involved in radiotherapy-from physical radiation interactions occurring over nanometres to evolution of patient responses over months and years. This review presents the current status and ongoing research in modelling radiotherapy responses across these scales, including basic physical mechanisms of DNA damage, the immediate biological responses this triggers, and genetic- and patient-level determinants of response. Finally, some of the major challenges in this field and potential avenues for future improvements are also discussed.

Proton FLASH: Impact of Dose Rate and Split Dose on Acute Skin Toxicity in a Murine Model

PURPOSE

Preclinical studies have shown a preferential normal tissue sparing effect of FLASH radiation therapy with ultra-high dose rates. The aim of the present study was to use a murine model of acute skin toxicity to investigate the biologic effect of varying dose rates, time structure, and introducing pauses in the dose delivery.

METHODS AND MATERIALS

The right hind limbs of nonanaesthetized mice were irradiated in the entrance plateau of a pencil beam scanning proton beam with 39.3 Gy. Experiment 1 was with varying field dose rates (0.7-80 Gy/s) without repainting, experiment 2 was with varying field dose rates (0.37-80 Gy/s) with repainting, and in experiment 3, the dose was split into 2, 3, 4, or 6 identical deliveries with 2-minute pauses. In total, 320 mice were included, with 6 to 25 mice per group. The endpoints were skin toxicity of different levels up to 25 days after irradiation.

RESULTS

The dose rate50, which is the dose rate to induce a response in 50% of the animals, depended on the level of skin toxicity, with the higher toxicity levels displaying a FLASH effect at 0.7-2 Gy/s. Repainting resulted in higher toxicity for the same field dose rate. Splitting the dose into 2 deliveries reduced the FLASH effect, and for 3 or more deliveries, the FLASH effect was almost abolished for lower grades of toxicity.

CONCLUSIONS

The dose rate that induced a FLASH effect varied for different skin toxicity levels, which are characterized by a differing degree of sensitivity to radiation dosage. Conclusions on a threshold for the dose rate needed to obtain a FLASH effect can therefore be influenced by the dose sensitivity of the used endpoint. Splitting the total dose into more deliveries compromised the FLASH effect. This can have an impact for fractionation as well as for regions where 2 or more FLASH fields overlap within the same treatment session.

Particle Beam Radiobiology Status and Challenges: A PTCOG Radiobiology Subcommittee Report

Particle therapy (PT) represents a significant advancement in cancer treatment, precisely targeting tumor cells while sparing surrounding healthy tissues thanks to the unique depth-dose profiles of the charged particles. Furthermore, their linear energy transfer and relative biological effectiveness enhance their capability to treat radioresistant tumors, including hypoxic ones. Over the years, extensive research has paved the way for PT's clinical application, and current efforts aim to refine its efficacy and precision, minimizing the toxicities. In this regard, radiobiology research is evolving toward integrating biotechnology to advance drug discovery and radiation therapy optimization. This shift from basic radiobiology to understanding the molecular mechanisms of PT aims to expand the therapeutic window through innovative dose delivery regimens and combined therapy approaches. This review, written by over 30 contributors from various countries, provides a comprehensive look at key research areas and new developments in PT radiobiology, emphasizing the innovations and techniques transforming the field, ranging from the radiobiology of new irradiation modalities to multimodal radiation therapy and modeling efforts. We highlight both advancements and knowledge gaps, with the aim of improving the understanding and application of PT in oncology.

Oxygen Enhancement Ratio-Weighted Dose Quantitatively Describes Acute Skin Toxicity Variations in Mice After Pencil Beam Scanning Proton FLASH Irradiation With Changing Doses and Time Structures

PURPOSE

The aim of this work was to investigate the ability of a biological oxygen enhancement ratio-weighted dose, DOER, to describe acute skin toxicity variations observed in mice after proton pencil beam scanning irradiations with changing doses and beam time structures.

METHODS AND MATERIALS

In five independent experiments, the right hind leg of a total of 621 CDF1 mice was irradiated previously in the entrance plateau of a pencil beam scanning proton beam. The incidence of acute skin toxicity (of level 1.5-2.0-2.5-3.0-3.5) was scored for 47 different mouse groups that mapped toxicity as function of dose for conventional and FLASH dose rate, toxicity as function of field dose rate with and without repainting, and toxicity when splitting the treatment into 1 to 6 identical deliveries separated by 2 minutes. DOER was calculated for all mouse groups using a simple oxygen kinetics model to describe oxygen depletion. The three independent model parameters (oxygen-depletion rate, oxygen-recovery rate, oxygen level without irradiation) were fitted to the experimental data. The ability of DOER to describe the toxicity variations across all experiments was investigated by comparing DOER-response curves across the five independent experiments.

RESULTS

After conversion from the independent variable tested in each experiment to DOER, all five experiments had similar MDDOER50 (DOER giving 50% toxicity incidence) with standard deviations of 0.45 - 1.6 Gy for the five toxicity levels. DOER could thus describe the observed toxicity variations across all experiments.

CONCLUSIONS

DOER described the varying FLASH-sparing effect observed for a wide range of conditions. Calculation of DOER for other irradiation conditions can quantitatively estimate the FLASH-sparing effect for arbitrary irradiations for the investigated murine model. With appropriate fitting parameters DOER also may be able to describe FLASH effect variations with dose and dose rate for other assays and endpoints.

How quickly does FLASH need to be delivered? A theoretical study of radiolytic oxygen depletion kinetics in tissues

Purpose. Radiation delivered over ultra-short timescales ('FLASH' radiotherapy) leads to a reduction in normal tissue toxicities for a range of tissues in the preclinical setting. Experiments have shown this reduction occurs for total delivery times less than a 'critical' time that varies by two orders of magnitude between brain (∼0.3 s) and skin (⪆10 s), and three orders of magnitude across different bowel experiments, from ∼0.01 to ⪆(1-10) s. Understanding the factors responsible for this broad variation may be important for translation of FLASH into the clinic and understanding the mechanisms behind FLASH.Methods.Assuming radiolytic oxygen depletion (ROD) to be the primary driver of FLASH effects, oxygen diffusion, consumption, and ROD were evaluated numerically for simulated tissues with pseudorandom vasculatures for a range of radiation delivery times, capillary densities, and oxygen consumption rates (OCR's). The resulting time-dependent oxygen partial pressure distribution histograms were used to estimate cell survival in these tissues using the linear quadratic model, modified to incorporate oxygen-enhancement ratio effects.Results. Independent of the capillary density, there was a substantial increase in predicted cell survival when the total delivery time was less than the capillary oxygen tension (mmHg) divided by the OCR (expressed in units of mmHg/s), setting the critical delivery time for FLASH in simulated tissues. Using literature OCR values for different normal tissues, the predicted range of critical delivery times agreed well with experimental values for skin and brain and, modifying our model to allow for fluctuating perfusion, bowel.Conclusions. The broad three-orders-of-magnitude variation in critical irradiation delivery times observed inin vivopreclinical experiments can be accounted for by the ROD hypothesis and differences in the OCR amongst simulated normal tissues. Characterization of these may help guide future experiments and open the door to optimized tissue-specific clinical protocols.

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.

Mini-GRID radiotherapy on the CLEAR very-high-energy electron beamline: collimator optimization, film dosimetry, and Monte Carlo simulations

Objective.Spatially-fractionated radiotherapy (SFRT) delivered with a very-high-energy electron (VHEE) beam and a mini-GRID collimator was investigated to achieve synergistic normal tissue-sparing through spatial fractionation and the FLASH effect.Approach.A tungsten mini-GRID collimator for delivering VHEE SFRT was optimized using Monte Carlo (MC) simulations. Peak-to-valley dose ratios (PVDRs), depths of convergence (DoCs, PVDR ≤ 1.1), and peak and valley doses in a water phantom from a simulated 150 MeV VHEE source were evaluated. Collimator thickness, hole width, and septal width were varied to determine an optimal value for each parameter that maximized PVDR and DoC. The optimized collimator (20 mm thick rectangular prism with a 15 mm × 15 mm face with a 7 × 7 array of 0.5 mm holes separated by 1.1 mm septa) was 3D-printed and used for VHEE irradiations with the CERN linear electron accelerator for research beam. Open beam and mini-GRID irradiations were performed at 140, 175, and 200 MeV and dose was recorded with radiochromic films in a water tank. PVDR, central-axis (CAX) and valley dose rates and DoCs were evaluated.Main results.Films demonstrated peak and valley dose rates on the order of 100 s of MGy/s, which could promote FLASH-sparing effects. Across the three energies, PVDRs of 2-4 at 13 mm depth and DoCs between 39 and 47 mm were achieved. Open beam and mini-GRID MC simulations were run to replicate the film results at 200 MeV. For the mini-GRID irradiations, the film CAX dose was on average 15% higher, the film valley dose was 28% higher, and the film PVDR was 15% lower than calculated by MC.Significance.Ultimately, the PVDRs and DoCs were determined to be too low for a significant potential for SFRT tissue-sparing effects to be present, particularly at depth. Further beam delivery optimization and investigations of new means of spatial fractionation are warranted.

Proton and Electron Ultrahigh-Dose-Rate Isodose Irradiations Produce Differences in Reactive Oxygen Species Yields

Purpose:

Investigations into ultra-high dose rate (UHDR) radiotherapy have dramatically risen because of the observed normal tissue sparing FLASH effect without sacrificing tumor control. The purpose of this study was to provide a direct beamline comparison of protons and electrons to determine where UHDR to conventional dose rates (CDR) changes affect the resultant radiochemistry.

Methods and Materials:

We used well characterized assays of reactive oxygen species (ROS) and oxygen consumption to assess the radiolysis in protein solutions. Three optical reporters related to ROS (CellROX Deep Red, reflects highly reactive radicals; Amplex Red reflects H2O2; and Oxyphor reflects partial pressure loss (ΔpO2)). A Varian ProBeam proton cyclotron and a converted Varian Trilogy electron linac were used for irradiation at both their CDR and UHDR capable level, to assess the assay change per unit dose.

Results:

For both protons and electrons an expected reduction in ROS was noted going from CDR to UHDR, and results interpreted as a reduction in the ratio of UHDR/CDR yield. The CellROX assay showed no difference between beamlines, each showing ~80% reduction in ROS from CDR to UHDR. The Amplex assay showed the largest inter-beamline difference, with ~5% loss using protons vs ~69% loss with electrons, in protein solution. The Oxyphor assay of ΔpO2 showed a small difference in CDR to UHDR with a 23% loss with protons and 43% loss with electrons.

Conclusion:

Interpretation of ROS assays and oxygen consumption is notoriously challenging. These assays might be interpreted by their most activating species’ lifetime. The assay for highly reactive OH●, appeared independent of beamline, whereas the assays for the longer lived H2O2 species and ΔpO2 assay showed differences between beamlines via the UHDR/CDR ratio. This work can be used for FLASH hypothesis testing by comparing these assays to isodose biological FLASH effects in vivo.

Oxygen supplementation in anesthesia can block FLASH effect and anti-tumor immunity in conventional proton therapy

BACKGROUND

Radiation-induced neurocognitive dysfunction is a major adverse effect of brain radiation therapy and has specific relevance in pediatric oncology, where serious cognitive deficits have been reported in survivors of pediatric brain tumors. Moreover, many pediatric patients receive proton therapy under general anesthesia or sedation to guarantee precise ballistics with a high oxygen content for safety. The present study addresses the relevant question of the potential effect of supplemental oxygen administered during anesthesia on normal tissue toxicity and investigates the anti-tumor immune response generated following conventional and FLASH proton therapy.

METHODS

Rats (Fischer 344) were cranially irradiated with a single high dose of proton therapy (15 Gy or 25 Gy) using FLASH dose rate proton irradiation (257 ± 2 Gy/s) or conventional dose rate proton irradiation (4 ± 0.02 Gy/s), and the toxicities in the normal tissue were examined by histological, cytometric and behavioral analysis. Glioblastoma-bearing rats were irradiated in the same manner and tumor-infiltrating leukocytes were quantified by flow cytometry.

RESULTS

Our findings indicate that supplemental oxygen has an adverse impact on both functional and anatomical evaluations of normal brain following conventional and FLASH proton therapy. In addition, oxygen supplementation in anesthesia is particularly detrimental for anti-tumor immune response by preventing a strong immune cell infiltration into tumoral tissues following conventional proton therapy.

CONCLUSIONS

These results demonstrate the need to further optimize anesthesia protocols used in radiotherapy with the goal of preserving normal tissues and achieving tumor control, specifically in combination with immunotherapy agents.

The DRESDEN PLATFORM is a research hub for ultra-high dose rate radiobiology

The recently observed FLASH effect describes the observation of normal tissue protection by ultra-high dose rates (UHDR), or dose delivery in a fraction of a second, at similar tumor-killing efficacy of conventional dose delivery and promises great benefits for radiotherapy patients. Dedicated studies are now necessary to define a robust set of dose application parameters for FLASH radiotherapy and to identify underlying mechanisms. These studies require particle accelerators with variable temporal dose application characteristics for numerous radiation qualities, equipped for preclinical radiobiological research. Here we present the DRESDEN PLATFORM, a research hub for ultra-high dose rate radiobiology. By uniting clinical and research accelerators with radiobiology infrastructure and know-how, the DRESDEN PLATFORM offers a unique environment for studying the FLASH effect. We introduce its experimental capabilities and demonstrate the platform's suitability for systematic investigation of FLASH by presenting results from a concerted in vivo radiobiology study with zebrafish embryos. The comparative pre-clinical study was conducted across one electron and two proton accelerator facilities, including an advanced laser-driven proton source applied for FLASH-relevant in vivo irradiations for the first time. The data show a protective effect of UHDR irradiation up to [Formula: see text] and suggests consistency of the protective effect even at escalated dose rates of [Formula: see text]. With the first clinical FLASH studies underway, research facilities like the DRESDEN PLATFORM, addressing the open questions surrounding FLASH, are essential to accelerate FLASH's translation into clinical practice.

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.

A feasibility study of ultra-high dose rate mini-GRID therapy using very-high-energy electron beams for a simulated pediatric brain case

Ultra-high dose rate (UHDR, >40 Gy/s), spatially-fractionated minibeam GRID (mini-GRID) therapy using very-high-energy electrons (VHEE) was investigated using Monte Carlo simulations. Multi-directional VHEE treatments with and without mini-GRID-fractionation were compared to a clinical 6 MV volumetric modulated arc therapy (VMAT) plan for a pediatric glioblastoma patient using dose-volume histograms, volume-averaged dose rates in critical patient structures, and planning target volume D98s. Peak-to-valley dose ratios (PVDRs) and dose rates in organs at risk (OARs) were evaluated due to their relevance for normal-tissue sparing in FLASH and spatially-fractionated techniques. Depths of convergence, defined where the PVDR is first ≤1.1, and depths at which dose rates fall below the UHDR threshold were also evaluated. In a water phantom, the VHEE mini-GRID treatments presented a surface (5 mm depth) PVDR of (51±2) and a depth of convergence of 42 mm at 150 MeV and a surface PVDR of (33±1) with a depth of convergence of 57 mm at 250 MeV. For a pediatric GBM case, VHEE treatments without mini-GRID-fractionation produced 25% and 22% lower volume-averaged doses to OARs compared to the 6 MV VMAT plan and 8/9 and 9/9 of the patient structures were exposed to volume-averaged dose rates >40 Gy/s for the 150 MeV and 250 MeV plans, respectively. The 150 MeV and 250 MeV mini-GRID treatments produced 17% and 38% higher volume-averaged doses to OARs and 3/9 patient structures had volume-averaged dose rates above 40 Gy/s. VHEE mini-GRID plans produced many comparable dose metrics to the clinical VMAT plan, encouraging further optimization.

GEANT4-DNA simulation of temperature-dependent and pH-dependent yields of chemical radiolytic species

Objective.GEANT4-DNA can simulate radiation chemical yield (G-value) for radiolytic species such as the hydrated electron (eaq-) with the independent reaction times (IRT) method, however, only at room temperature and neutral pH. This work aims to modify the GEANT4-DNA source code to enable the calculation ofG-values for radiolytic species at different temperatures and pH values.Approach.In the GEANT4-DNA source code, values of chemical parameters such as reaction rate constant, diffusion coefficient, Onsager radius, and water density were replaced by corresponding temperature-dependent polynomials. The initial concentration of hydrogen ion (H+)/hydronium ion (H3O+) was scaled for a desired pH using the relationship pH = -log10[H+]. To validate our modifications, two sets of simulations were performed. (A) A water cube with 1.0 km sides and a pH of 7 was irradiated with an isotropic electron source of 1 MeV. The end time was 1μs. The temperatures varied from 25 °C to 150 °C. (B) The same setup as (A) was used, however, the temperature was set to 25 °C while the pH varied from 5 to 9. The results were compared with published experimental and simulated work.Main results.The IRT method in GEANT4-DNA was successfully modified to simulateG-values for radiolytic species at different temperatures and pH values. Our temperature-dependent results agreed with experimental data within 0.64%-9.79%, and with simulated data within 3.52%-12.47%. The pH-dependent results agreed well with experimental data within 0.52% to 3.19% except at a pH of 5 (15.99%) and with simulated data within 4.40%-5.53%. The uncertainties were below ±0.20%. Overall our results agreed better with experimental than simulation data.Significance.Modifications in the GEANT4-DNA code enabled the calculation ofG-values for radiolytic species at different temperatures and pH values.

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.

The general-purpose Geant4 Monte Carlo toolkit and its Geant4-DNA extension to investigate mechanisms underlying the FLASH effect in radiotherapy: Current status and challenges

FLASH radiotherapy is a promising approach to cancer treatment that offers several advantages over conventional radiotherapy. With this novel technique, high doses of radiation are delivered in a short period of time, inducing the so-called FLASH effect - a phenomenon characterized by healthy tissue sparing without alteration of tumor control. The mechanisms behind the FLASH effect remain unknown. One way to approach this problem is to gain insight into the initial parameters that can distinguish FLASH from conventional irradiation by simulating particle transport in aqueous media using the general-purpose Geant4 Monte Carlo toolkit and its Geant4-DNA extension. This review article discusses the current status of Geant4 and Geant4-DNA simulations to investigate mechanisms underlying the FLASH effect, as well as the challenges faced in this research field. One of the primary challenges is to accurately simulate the experimental irradiation parameters. Another challenge is the temporal extension of the simulations. This review also focuses on two hypotheses to explain the FLASH effect - namely the oxygen depletion hypothesis and the inter-track interactions hypothesis - and discusses how the Geant4 toolkit can be used to investigate them. The aim of this review is to provide an overview of Geant4 and Geant4-DNA simulations for FLASH radiotherapy and to highlight the challenges that need to be overcome in order to better study the FLASH effect.

The current status of FLASH particle therapy: a systematic review

Particle therapies are becoming increasingly available clinically due to their beneficial energy deposition profile, sparing healthy tissues. This may be further promoted with ultra-high dose rates, termed FLASH. This review comprehensively summarises current knowledge based on studies relevant to proton- and carbon-FLASH therapy. As electron-FLASH literature presents important radiobiological findings that form the basis of proton and carbon-based FLASH studies, a summary of key electron-FLASH papers is also included. Preclinical data suggest three key mechanisms by which proton and carbon-FLASH are able to reduce normal tissue toxicities compared to conventional dose rates, with equipotent, or enhanced, tumour kill efficacy. However, a degree of caution is needed in clinically translating these findings as: most studies use transmission and do not conform the Bragg peak to tumour volume; mechanistic understanding is still in its infancy; stringent verification of dosimetry is rarely provided; biological assays are prone to limitations which need greater acknowledgement.

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.

Reinventing Radiobiology in the Light of FLASH Radiotherapy

Ultrahigh-dose rate FLASH radiotherapy (FLASH-RT) is a potentially paradigm-shifting treatment modality that holds the promise of expanding the therapeutic index for nearly any cancer. At the heart of this exciting technology comes the capability to ameliorate major normal tissue complications without compromising the efficacy of tumor killing. This combination of benefits has now been termed the FLASH effect and relies on an in vivo validation to rigorously demonstrate the absence of normal tissue toxicity. The FLASH effect occurs when the overall irradiation time is extremely short (<500 ms), and in this review we attempt to understand how FLASH-RT can kill tumors but spare normal tissues-likely the single most pressing question confronting the field today.

Multicomponent mathematical model for tumor volume calculation with setup error using single-isocenter stereotactic radiotherapy for multiple brain metastases

We evaluated the tumor residual volumes considering six degrees-of-freedom (6DoF) patient setup errors in stereotactic radiotherapy (SRT) with multicomponent mathematical model using single-isocenter irradiation for brain metastases. Simulated spherical gross tumor volumes (GTVs) with 1.0 (GTV 1), 2.0 (GTV 2), and 3.0 (GTV 3)-cm diameters were used. The distance between the GTV center and isocenter (d) was set at 0-10 cm. The GTV was simultaneously translated within 0-1.0 mm (T) and rotated within 0°-1.0° (R) in the three axis directions using affine transformation. We optimized the tumor growth model parameters using measurements of non-small cell lung cancer cell lines' (A549 and NCI-H460) growth. We calculated the GTV residual volume at the irradiation's end using the physical dose to the GTV when the GTV size, d, and 6DoF setup error varied. The d-values that satisfy tolerance values (10%, 35%, and 50%) of the GTV residual volume rate based on the pre-irradiation GTV volume were determined. The larger the tolerance value set for both cell lines, the longer the distance to satisfy the tolerance value. In GTV residual volume evaluations based on the multicomponent mathematical model on SRT with single-isocenter irradiation, the smaller the GTV size and the larger the distance and 6DoF setup error, the shorter the distance that satisfies the tolerance value might need to be.

Modeling of scavenging systems in water radiolysis with Geant4-DNA

PURPOSE

This paper presents the capabilities of the Geant4-DNA Monte Carlo toolkit to simulate water radiolysis with scavengers using the step-by-step (SBS) or the independent reaction times (IRT) methods. It features two examples of application areas: (1) computing the escape yield of H₂O₂ following a ⁶⁰Co γ-irradiation and (2) computing the oxygen depletion in water irradiated with 1 MeV electrons.

METHODS

To ease the implementation of the chemical stage in Geant4-DNA, we developed a user interface that helps define the chemical reactions and set the concentration of scavengers. The first application area example required two computational steps to perform water radiolysis using NO₂^- and NO₃^- as scavengers and a ⁶⁰Co irradiation. The oxygen depletion computation technique for the second application area example consisted of simulating track segments of 1 MeV electrons and determining the radio-induced loss and gain of oxygen molecules.

RESULTS

The production of H₂O₂ under variable scavenging levels is consistent with the literature; the mean relative difference between the SBS and IRT methods is 7.2 % ± 0.5 %. For the oxygen depletion 1 µs post-irradiation, the mean relative difference between both methods is equal to 9.8 % ± 0.3 %. The results in the microsecond scale depend on the initial partial pressure of oxygen in water. In addition, the computed oxygen depletions agree well with the literature.

CONCLUSIONS

The Geant4-DNA toolkit makes it possible to simulate water radiolysis in the presence of scavengers. This feature offers perspectives in radiobiology, with the possibility of simulating cell-relevant scavenging mechanisms.

Flash radiotherapy-gateway to promised land or another mirage

Radiation therapy damages cancer cells with ionizing radiation, leading to their death. However, radiation‑induced toxicity limits the dose delivered to the tumor, thereby constraining the control effect of radiotherapy n tumor growth. In addition, the delayed toxicity caused by radiotherapy significantly harms the physical and mental health of patients. FLASH‑RT, an emerging class of radiotherapy, causes a phenomenon known as the 'FLASH effect', which delivers radiotherapy at an ultra‑high dose rate with lower toxicity to normal tissue than conventional radiotherapy to achieve local tumor control.

Investigating the potential contribution of inter-track interactions within ultra-high dose-rate proton therapy

Objective. Laser-accelerated protons offer an alternative delivery mechanism for proton therapy. This technique delivers dose-rates of ≥10⁹Gy s^-1, many orders of magnitude greater than used clinically. Such ultra-high dose-rates reduce delivery time to nanoseconds, equivalent to the lifetime of reactive chemical species within a biological medium. This leads to the possibility of inter-track interactions between successive protons within a pulse, potentially altering the yields of damaging radicals if they are in sufficient spatial proximity. This work investigates the temporal evolution of chemical species for a range of proton energies and doses to quantify the circumstances required for inter-track interactions, and determine any relevance within ultra-high dose-rate proton therapy.Approach. The TOPAS-nBio Monte Carlo toolkit was used to investigate possible inter-track interactions. Firstly, protons between 0.5 and 100 MeV were simulated to record the radial track dimensions throughout the chemical stage from 1 ps to 1μs. Using the track areas, the geometric probability of track overlap was calculated for various exposures and timescales. A sample of irradiations were then simulated in detail to compare any change in chemical yields for independently and instantaneously delivered tracks, and validate the analytic model.Main results. Track overlap for a clinical 2 Gy dose was negligible for biologically relevant timepoints for all energies. Overlap probability increased with time after irradiation, proton energy and dose, with a minimum 23 Gy dose required before significant track overlap occurred. Simulating chemical interactions confirmed these results with no change in radical yields seen up to 8 Gy for independently and instantaneously delivered tracks.Significance. These observations suggest that the spatial separation between incident protons is too large for physico-chemical inter-track interactions, regardless of the delivery time, indicating such interactions would not play a role in any potential changes in biological response between laser-accelerated and conventional proton therapy.

Induction of DNA strand breaks and oxidative base damages in plasmid DNA by ultra-high dose rate proton irradiation

PURPOSE

Radiation cancer therapy with ultra-high dose rate (UHDR) exposure, so-called FLASH radiotherapy, appears to reduce normal tissue damage without compromising tumor response to therapy. The aim of this study was to clarify whether a 59.5 MeV proton beam at an UHDR of 48.6 Gy/s could effectively reduce the DNA damage of pBR322 plasmid DNA in solution compared to the conventional dose rate (CONV) of 0.057 Gy/s.

MATERIALS AND METHODS

A simple system, consisting of pBR322 plasmid DNA in 1× Tris-EDTA buffer, was initially employed for proton beam exposure. We then used formamidopyrimidine-DNA glycosylase (Fpg) enzymes. which convert oxidative base damages of oxidized purines to DNA strand breaks, to quantify DNA single strand breaks (SSBs) and double strand breaks (DSBs) by agarose gel electrophoresis.

RESULTS

Our findings showed that the SSB induction rate (SSB per plasmid DNA/Gy) at UHDR and the induction of Fpg enzyme sensitive sites (ESS) were significantly reduced in UHDR compared to CONV. However, there was no significant difference in DSB induction and non-DSB cluster damages.

CONCLUSIONS

UHDR of a 59.5 MeV proton beam could reduce non-clustered, non-DSB damages, such as SSB and sparsely distributed ESS. However, this effect may not be significant in reducing lethal DNA damage that becomes apparent only in acute radiation effects of mammalian cells and in vivo studies.

Validation of a High-Throughput Dicentric Chromosome Assay Using Complex Radiation Exposures

Validation of biodosimetry assays is routinely performed using primarily orthovoltage irradiators at a conventional dose rate of approximately 1 Gy/min. However, incidental/ accidental exposures caused by nuclear weapons can be more complex. The aim of this work was to simulate the DNA damage effects mimicking those caused by the detonation of a several kilotons improvised nuclear device (IND). For this, we modeled complex exposures to: 1. a mixed (photons + IND-neutrons) field and 2. different dose rates that may come from the blast, nuclear fallout, or ground deposition of radionuclides (ground shine). Additionally, we assessed whether myeloid cytokines affect the precision of radiation dose estimation by modulating the frequency of dicentric chromosomes. To mimic different exposure scenarios, several irradiation systems were used. In a mixed field study, human blood samples were exposed to a photon field enriched with neutrons (ranging from 10% to 37%) from a source that mimics Hiroshima's A-bomb's energy spectrum (0.2-9 MeV). Using statistical analysis, we assessed whether photons and neutrons act in an additive or synergistic way to form dicentrics. For the dose rates study, human blood was exposed to photons or electrons at dose rates ranging from low (where the dose was spread over 32 h) to extremely high (where the dose was delivered in a fraction of a microsecond). Potential effects of cytokine treatment on biodosimetry dose predictions were analyzed in irradiated blood subjected to Neupogen or Neulasta for 24 or 48 h at the concentration recommended to forestall manifestation of an acute radiation syndrome in bomb survivors. All measurements were performed using a robotic station, the Rapid Automated Biodosimetry Tool II, programmed to culture lymphocytes and score dicentrics in multiwell plates (the RABiT-II DCA). In agreement with classical concepts of radiation biology, the RABiT-II DCA calibration curves suggested that the frequency of dicentrics depends on the type of radiation and is modulated by changes in the dose rate. The resulting dose-response curves suggested an intermediate dicentric yields and additive effects of photons and IND-neutrons in the mixed field. At ultra-high dose rate (600 Gy/s), affected lymphocytes exhibited significantly fewer dicentrics (P < 0.004, t test). In contrast, we did not find the dose-response modification effects of radiomitigators on the yields of dicentrics (Bonferroni corrected P > 0.006, ANOVA test). This result suggests no bias in the dose predictions should be expected after emergency cytokine treatment initiated up to 48 h prior to blood collection for dicentric analysis.

FLASH radiotherapy: A promising new method for radiotherapy

Among the treatments for malignant tumors, radiotherapy is of great significance both as a main treatment and as an adjuvant treatment. Radiation therapy damages cancer cells with ionizing radiation, leading to their death. However, radiation-induced toxicity limits the dose delivered to the tumor, thereby constraining the control effect of radiotherapy on tumor growth. In addition, the delayed toxicity caused by radiotherapy significantly harms the physical and mental health of patients. FLASH-RT, an emerging class of radiotherapy, causes a phenomenon known as the 'FLASH effect', which delivers radiotherapy at an ultra-high dose rate with lower toxicity to normal tissue than conventional radiotherapy to achieve local tumor control. Although its mechanism remains to be fully elucidated, this modality constitutes a potential new approach to treating malignant tumors. In the present review, the current research progress of FLASH-RT and its various particular effects are described, including the status of research on FLASH-RT and its influencing factors. The hypothetic mechanism of action of FLASH-RT is also summarized, providing insight into future tumor treatments.

A phenomenological model of proton FLASH oxygen depletion effects depending on tissue vasculature and oxygen supply

Introduction

Radiation-induced oxygen depletion in tissue is assumed as a contributor to the FLASH sparing effects. In this study, we simulated the heterogeneous oxygen depletion in the tissue surrounding the vessels and calculated the proton FLASH effective-dose-modifying factor (FEDMF), which could be used for biology-based treatment planning.

Methods

The dose and dose-weighted linear energy transfer (LET) of a small animal proton irradiator was simulated with Monte Carlo simulation. We deployed a parabolic partial differential equation to account for the generalized radiation oxygen depletion, tissue oxygen diffusion, and metabolic processes to investigate oxygen distribution in 1D, 2D, and 3D solution space. Dose and dose rates, particle LET, vasculature spacing, and blood oxygen supplies were considered. Using a similar framework for the hypoxic reduction factor (HRF) developed previously, the FEDMF was derived as the ratio of the cumulative normoxic-equivalent dose (CNED) between CONV and UHDR deliveries.

Results

Dynamic equilibrium between oxygen diffusion and tissue metabolism can result in tissue hypoxia. The hypoxic region displayed enhanced radio-resistance and resulted in lower CNED under UHDR deliveries. In 1D solution, comparing 15 Gy proton dose delivered at CONV 0.5 and UHDR 125 Gy/s, 61.5% of the tissue exhibited ≥20% FEDMF at 175 μm vasculature spacing and 18.9 μM boundary condition. This percentage reduced to 34.5% and 0% for 8 and 2 Gy deliveries, respectively. Similar trends were observed in the 3D solution space. The FLASH versus CONV differential effect remained at larger vasculature spacings. A higher FLASH dose rate showed an increased region with ≥20% FEDMF. A higher LET near the proton Bragg peak region did not appear to alter the FLASH effect.

Conclusion

We developed 1D, 2D, and 3D oxygen depletion simulation process to obtain the dynamic HRF and derive the proton FEDMF related to the dose delivery parameters and the local tissue vasculature information. The phenomenological model can be used to simulate or predict FLASH effects based on tissue vasculature and oxygen concentration data obtained from other experiments.

Radical Production with Pulsed Beams: Understanding the Transition to FLASH

Ultra-high dose rate (UHDR) irradiation regimes have the potential to spare normal tissue while keeping equivalent tumoricidal capacity than conventional dose rate radiotherapy (CONV-RT). This has been called the FLASH effect. In this work, we present a new simulation framework aiming to study the production of radical species in water and biological media under different irradiation patterns. The chemical stage (heterogeneous phase) is based on a nonlinear reaction-diffusion model, implemented in GPU. After the first 1 μs, no further radical diffusion is assumed, and radical evolution may be simulated over long periods of hundreds of seconds. Our approach was first validated against previous results in the literature and then employed to assess the influence of different temporal microstructures of dose deposition in the expected biological damage. The variation of the Normal Tissue Complication Probability (NTCP), assuming the model of Labarbe et al., where the integral of the peroxyl radical concentration over time (AUC-ROO) is taken as surrogate for biological damage, is presented for different intra-pulse dose rate and pulse frequency configurations, relevant in the clinical scenario. These simulations yield that overall, mean dose rate and the dose per pulse are the best predictors of biological effects at UHDR.

Potential Molecular Mechanisms behind the Ultra-High Dose Rate "FLASH" Effect

FLASH radiotherapy, or the delivery of a dose at an ultra-high dose rate (>40 Gy/s), has recently emerged as a promising tool to enhance the therapeutic index in cancer treatment. The remarkable sparing of normal tissues and equivalent tumor control by FLASH irradiation compared to conventional dose rate irradiation—the FLASH effect—has already been demonstrated in several preclinical models and even in a first patient with T-cell cutaneous lymphoma. However, the biological mechanisms responsible for the differential effect produced by FLASH irradiation in normal and cancer cells remain to be elucidated. This is of great importance because a good understanding of the underlying radiobiological mechanisms and characterization of the specific beam parameters is required for a successful clinical translation of FLASH radiotherapy. In this review, we summarize the FLASH investigations performed so far and critically evaluate the current hypotheses explaining the FLASH effect, including oxygen depletion, the production of reactive oxygen species, and an altered immune response. We also propose a new theory that assumes an important role of mitochondria in mediating the normal tissue and tumor response to FLASH dose rates.

Monte Carlo optimization of a GRID collimator for preclinical megavoltage ultra-high dose rate spatially-fractionated radiation therapy

Objective. A 2-dimensional pre-clinical SFRT (GRID) collimator was designed for use on the ultra-high dose rate (UHDR) 10 MV ARIEL beamline at TRIUMF. TOPAS Monte Carlo simulations were used to determine optimal collimator geometry with respect to various dosimetric quantities. Approach. The GRID-averaged peak-to-valley dose ratio (PVDR) and mean dose rate of the peaks were investigated with the intent of maximizing both values in a given design. The effects of collimator thickness, focus position, septal width, and hole width on these metrics were found by testing a range of values for each parameter on a cylindrical GRID collimator. For each tested collimator geometry, photon beams with energies of 10, 5, and 1 MV were transported through the collimator and dose rates were calculated at various depths in a water phantom located 1.0 cm from the collimator exit. Main results. In our optimization, hole width proved to be the only collimator parameter which increased both PVDR and peak dose rates. From the optimization results, it was determined that our optimized design would be one which achieves the maximum dose rate for a PVDR ≥ 5 at 10 MV. Ultimately, this was achieved using a collimator with a thickness of 75 mm, 0.8 mm septal and hole widths, and a focus position matched to the beam divergence. This optimized collimator maintained the PVDR of 5 in the phantom between water depths of 0–10 cm at 10 MV and had a mean peak dose rate of 3.06 ± 0.02 Gy s − 1 at 0–1 cm depth. Significance. We have investigated the impact of various GRID-collimator design parameters on the dose rate and spatial fractionation of 10, 5, and 1 MV photon beams. The optimized collimator design for the 10 MV ultra-high dose rate photon beam could become a useful tool for radiobiology studies synergizing the effects of ultra-high dose rate (FLASH) delivery and spatial fractionation.

oChanges in Radical Levels as a Cause for the FLASH effect: Impact of beam structure parameters at ultra-high dose rates on oxygen depletion in water

The influence of different average and bunch dose rates in electron beams on the FLASH effect was investigated. The present study measures O2 content in water at different beam pulse patterns and finds strong correlation with biological data, strengthening the hypothesis of radical-related mechanisms as a reason for the FLASH effect.

Combining FLASH and spatially fractionated radiation therapy: The best of both worlds

FLASH radiotherapy (FLASH-RT) and spatially fractionated radiation therapy (SFRT) are two new therapeutical strategies that use non-standard dose delivery methods to reduce normal tissue toxicity and increase the therapeutic index. Although likely based on different mechanisms, both FLASH-RT and SFRT have shown to elicit radiobiological effects that significantly differ from those induced by conventional radiotherapy. With the therapeutic potential having been established separately for each technique, the combination of FLASH-RT and SFRT could therefore represent a winning alliance. In this review, we discuss the state of the art, advantages and current limitations, potential synergies, and where a combination of these two techniques could be implemented today or in the near future.

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.

Using oxygen dose histograms to quantify voxelised ultra-high dose rate (FLASH) effects in multiple radiation modalities

Purpose.To introduce a methodology to predict tissue sparing effects in pulsed ultra-high dose rate radiation exposures which could be included in a dose-effect prediction system or treatment planning system and to illustrate it by using three published experiments.Methods and materials.The proposed system formalises the variability of oxygen levels as an oxygen dose histogram (ODH), which provides an instantaneous oxygen level at a delivered dose. The histogram concept alleviates the need for a mechanistic approach. At each given oxygen level the oxygen fixation concept is used to calculate the change in DNA-damage induction compared to the fully hypoxic case. Using the ODH concept it is possible to estimate the effect even in the case of multiple pulses, partial oxygen depletion, and spatial oxygen depletion. The system is illustrated by applying it to the seminal results by Town (Nat. 1967) on cell cultures and the pre-clinical experiment on cognitive effects by Montay-Gruelet al(2017Radiother. Oncol.124365-9).Results.The proposed system predicts that a possible FLASH-effect depends on the initial oxygenation level in tissue, the total dose delivered, pulse length and pulse repetition rate. The magnitude of the FLASH-effect is the result of a redundant system, in that it will have the same specific value for a different combination of these dependencies. The cell culture data are well represented, while a correlation between the pre-clinical experiments and the calculated values is highly significant (p < 0.01).Conclusions. A system based only on oxygen related effects is able to quantify most of the effects currently observed in FLASH-radiation.

Modeling the impact of spatial oxygen heterogeneity on radiolytic oxygen depletion during FLASH radiotherapy

Purpose.It has been postulated that the delivery of radiotherapy at ultra-high dose rates ('FLASH') reduces normal tissue toxicities by depleting them of oxygen. The fraction of normal tissue and cancer cells surviving radiotherapy depends on dose and oxygen levels in an exponential manner and even a very small fraction of tissue at low oxygen levels can determine radiotherapy response. To quantify the differential impact of FLASH radiotherapy on normal and tumour tissues, the spatial heterogeneity of oxygenation in tissue should thus be accounted for.Methods.The effect of FLASH on radiation-induced normal and tumour tissue cell killing was studied by simulating oxygen diffusion, metabolism, and radiolytic oxygen depletion (ROD) over domains with simulated capillary architectures. To study the impact of heterogeneity, two architectural models were used: (1) randomly distributed capillaries and (2) capillaries forming a regular square lattice array. The resulting oxygen partial pressure distribution histograms were used to simulate normal and tumour tissue cell survival using the linear quadratic model of cell survival, modified to incorporate oxygen-enhancement ratio effects. The ratio ('dose modifying factors') of conventional low-dose-rate dose and FLASH dose at iso-cell survival was computed and compared with empirical iso-toxicity dose ratios.Results.Tumour cell survival was found to be increased by FLASH as compared to conventional radiotherapy, with a 0-1 order of magnitude increase for expected levels of tumour hypoxia, depending on the relative magnitudes of ROD and tissue oxygen metabolism. Interestingly, for the random capillary model, the impact of FLASH on well-oxygenated (normal) tissues was found to be much greater, with an estimated increase in cell survival by up to 10 orders of magnitude, even though reductions in mean tissue partial pressure were modest, less than ∼7 mmHg for the parameter values studied. The dose modifying factor for normal tissues was found to lie in the range 1.2-1.7 for a representative value of normal tissue oxygen metabolic rate, consistent with preclinical iso-toxicity results.Conclusions.The presence of very small nearly hypoxic regions in otherwise well-perfused normal tissues with high mean oxygen levels resulted in a greater proportional sparing of normal tissue than tumour cells during FLASH irradiation, possibly explaining empirical normal tissue sparing and iso-tumour control results.

Image guidance for FLASH radiotherapy

FLASH radiotherapy (FLASH-RT) is an emerging ultra-high dose (>40 Gy/s) delivery that promises to improve the therapeutic potential by limiting toxicities compared to conventional RT while maintaining similar tumor eradication efficacy. Image guidance is an essential component of modern RT that should be harnessed to meet the special emerging needs of FLASH-RT and its associated high risks in planning and delivering of such ultra-high doses in short period of times. Hence, this contribution will elaborate on the imaging requirements and possible solutions in the entire chain of FLASH-RT treatment, from the planning, through the setup and delivery with online in vivo imaging and dosimetry, up to the assessment of biological mechanisms and treatment response. In patient setup and delivery, higher temporal sampling than in conventional RT should ensure that the short treatment is delivered precisely to the targeted region. Additionally, conventional imaging tools such as cone-beam computed tomography will continue to play an important role in improving patient setup prior to delivery, while techniques based on magnetic resonance imaging or positron emission tomography may be extremely valuable for either linear accelerator (Linac) or particle FLASH therapy, to monitor and track anatomical changes during delivery. In either planning or assessing outcomes, quantitative functional imaging could supplement conventional imaging for more accurate utilization of the biological window of the FLASH effect, selecting for or verifying things such as tissue oxygen and existing or transient hypoxia on the relevant timescales of FLASH-RT delivery. Perhaps most importantly at this time, these tools might help improve the understanding of the biological mechanisms of FLASH-RT response in tumor and normal tissues. The high dose deposition of FLASH provides an opportunity to utilize pulse-to-pulse imaging tools such as Cherenkov or radiation acoustic emission imaging. These could provide individual pulse mapping or assessing the 3D dose delivery superficially or at tissue depth, respectively. In summary, the most promising components of modern RT should be used for safer application of FLASH-RT, and new promising developments could be advanced to cope with its novel demands but also exploit new opportunities in connection with the unique nature of pulsed delivery at unprecedented dose rates, opening a new era of biological image guidance and ultrafast, pulse-based in vivo dosimetry.

Beam pulse structure and dose rate as determinants for the flash effect observed in zebrafish embryo

BACKGROUND AND PURPOSE

Continuing recent experiments at the research electron accelerator ELBE at the Helmholtz-Zentrum Dresden-Rossendorf the influence of beam pulse structure on the Flash effect was investigated.

MATERIALS AND METHODS

The proton beam pulse structure of an isochronous cyclotron (UHDR_iso) and a synchrocyclotron (UHDR_synchro) was mimicked at ELBE by quasi-continuous electron bunches at 13 MHz delivering mean dose rates of 287 Gy/s and 177 Gy/s and bunch dose rates of 10⁶Gy/s and 10⁹ Gy/s, respectively. For UHDR_synchro, 40 ms macro pulses at a frequency of 25 Hz superimposed the bunch delivery. For comparison, a maximum beam intensity (2.5 x 10⁵ Gy/s mean and ∼10⁹ Gy/s bunch dose rate) and a reference irradiation (of ∼8 Gy/min mean dose rate) were applied. Radiation induced changes were assessed in zebrafish embryos over four days post irradiation.

RESULTS

Relative to the reference a significant protecting Flash effect was observed for all electron beam pulse regimes with less severe damage the higher the mean dose rate of the electron beam. Accordingly, the macro pulsing induced prolongation of treatment time at UHDR_synchro regime reduces the protecting effect compared to the maximum regime delivered at same bunch but higher mean dose rate. The Flash effect of the UHDR_iso regime was confirmed at a clinical isochronous cyclotron comparing the damage induced by proton beams delivered at 300 Gy/s and ∼9 Gy/min.

CONCLUSION

The recent findings indicate that the mean dose rate or treatment time are decisive for the normal tissue protecting Flash effect in zebrafish embryo.

FLASH irradiation induces lower levels of DNA damage ex vivo, an effect modulated by oxygen tension, dose, and dose rate

OBJECTIVE

FLASH irradiation reportedly produces less normal tissue toxicity, while maintaining tumour response. To investigate oxygen's role in the 'FLASH effect', we assessed DNA damage levels following irradiation at different oxygen tensions, doses and dose rates.

METHODS

Samples of whole blood were irradiated (20 Gy) at various oxygen tensions (0.25-21%) with 6 MeV electrons at dose rates of either 2 kGy/s (FLASH) or 0.1 Gy/s (CONV), and subsequently with various doses (0-40 Gy) and intermediate dose rates (0.3-1000 Gy/s). DNA damage of peripheral blood lymphocytes (PBL) were assessed by the alkaline comet assay.

RESULTS

Following 20 Gy irradiation, lower levels of DNA damage were induced for FLASH, the difference being significant at 0.25% (p < 0.05) and 0.5% O2 (p < 0.01). The differential in DNA damage at 0.5% O2 was found to increase with total dose and dose rate, becoming significant for doses ≥20 Gy and dose rates ≥30 Gy/s.

CONCLUSION

This study shows, using the alkaline comet assay, that lower levels of DNA damage are induced following FLASH irradiation, an effect that is modulated by the oxygen tension, and increases with the total dose and dose rate of irradiation, indicating that an oxygen related mechanism, e.g. transient radiation-induced oxygen depletion, may contribute to the tissue sparing effect of FLASH irradiation.

ADVANCES IN KNOWLEDGE

This paper is first to directly show that FLASH-induced DNA damage is modulated by oxygen tension, total dose and dose rate, with FLASH inducing significantly lower levels of DNA damage for doses ≥20 Gy and dose rates ≥30 Gy/s, at 0.5% O2.

Development of a portable hypoxia chamber for ultra-high dose rate laser-driven proton radiobiology applications

Background

There is currently significant interest in assessing the role of oxygen in the radiobiological effects at ultra-high dose rates. Oxygen modulation is postulated to play a role in the enhanced sparing effect observed in FLASH radiotherapy, where particles are delivered at 40–1000 Gy/s. Furthermore, the development of laser-driven accelerators now enables radiobiology experiments in extreme regimes where dose rates can exceed 10⁹ Gy/s, and predicted oxygen depletion effects on cellular response can be tested. Access to appropriate experimental enviroments, allowing measurements under controlled oxygenation conditions, is a key requirement for these studies. We report on the development and application of a bespoke portable hypoxia chamber specifically designed for experiments employing laser-driven sources, but also suitable for comparator studies under FLASH and conventional irradiation conditions.

Materials and methods

We used oxygen concentration measurements to test the induction of hypoxia and the maintenance capacity of the chambers. Cellular hypoxia induction was verified using hypoxia inducible factor-1α immunostaining. Calibrated radiochromic films and GEANT-4 simulations verified the dosimetry variations inside and outside the chambers. We irradiated hypoxic human skin fibroblasts (AG01522B) cells with laser-driven protons, conventional protons and reference 225 kVp X-rays to quantify DNA DSB damage and repair under hypoxia. We further measured the oxygen enhancement ratio for cell survival after X-ray exposure in normal fibroblast and radioresistant patient- derived GBM stem cells.

Results

Oxygen measurements showed that our chambers maintained a radiobiological hypoxic environment for at least 45 min and pathological hypoxia for up to 24 h after disconnecting the chambers from the gas supply. We observed a significant reduction in the 53BP1 foci induced by laser-driven protons, conventional protons and X-rays in the hypoxic cells compared to normoxic cells at 30 min post-irradiation. Under hypoxic irradiations, the Laser-driven protons induced significant residual DNA DSB damage in hypoxic AG01522B cells compared to the conventional dose rate protons suggesting an important impact of these extremely high dose-rate exposures. We obtained an oxygen enhancement ratio (OER) of 2.1 ± 0.1 and 2.5 ± 0.1 respectively for the AG01522B and patient-derived GBM stem cells for X-ray irradiation using our hypoxia chambers.

Conclusion

We demonstrated the design and application of portable hypoxia chambers for studying cellular radiobiological endpoints after exposure to laser-driven protons at ultra-high dose, conventional protons and X-rays. Suitable levels of reduced oxygen concentration could be maintained in the absence of external gassing to quantify hypoxic effects. The data obtained provided indication of an enhanced residual DNA DSB damage under hypoxic conditions at ultra-high dose rate compared to the conventional protons or X-rays.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-022-02024-3.

Radiobiology of the FLASH effect

Radiation exposures at ultrahigh dose rates (UHDRs) at several orders of magnitude greater than in current clinical radiotherapy (RT) have been shown to manifest differential radiobiological responses compared to conventional (CONV) dose rates. This has led to studies investigating the application of UHDR for therapeutic advantage (FLASH-RT) that have gained significant interest since the initial discovery in 2014 that demonstrated reduced lung toxicity with equivalent levels of tumor control compared with conventional dose-rate RT. Many subsequent studies have demonstrated the potential protective role of FLASH-RT in normal tissues, yet the underlying molecular and cellular mechanisms of the FLASH effect remain to be fully elucidated. Here, we summarize the current evidence of the FLASH effect and review FLASH-RT studies performed in preclinical models of normal tissue response. To critically examine the underlying biological mechanisms of responses to UHDR radiation exposures, we evaluate in vitro studies performed with normal and tumor cells. Differential responses to UHDR versus CONV irradiation recurrently involve reduced inflammatory processes and differential expression of pro- and anti-inflammatory genes. In addition, frequently reduced levels of DNA damage or misrepair products are seen after UHDR irradiation. So far, it is not clear what signal elicits these differential responses, but there are indications for involvement of reactive species. Different susceptibility to FLASH effects observed between normal and tumor cells may result from altered metabolic and detoxification pathways and/or repair pathways used by tumor cells. We summarize the current theories that may explain the FLASH effect and highlight important research questions that are key to a better mechanistic understanding and, thus, the future implementation of FLASH-RT in the clinic.