A potential revolution in cancer treatment: A topical review of FLASH radiotherapy

A microscopic oxygen transport model for ultra-high dose rate radiotherapy in vivo: The impact of physiological conditions on FLASH effect

BACKGROUND

Ultra-high dose rate irradiation (≥40 Gy/s, FLASH) has been shown to reduce normal tissue toxicity, while maintaining tumor control compared to conventional dose-rate radiotherapy. The radiolytic oxygen (O2) depletion (ROD) resulting from FLASH has been proposed to explain the normal tissue protection effect; however, in vivo experiments have not confirmed that FLASH induced global tissue hypoxia. Nonetheless, the experiments reported are based on volume-averaged measurement, which have inherent limitations in detecting microscopic phenomena, including the potential preservation of stem cells niches due to local FLASH-induced O2 depletion. Computational modeling offers a complementary approach to understand the ROD caused by FLASH at the microscopic level.

PURPOSE

We developed a comprehensive model to describe the spatial and temporal dynamics of O2 consumption and transport in response to irradiation in vivo. The change of oxygen enhancement ratio (OER) was used to quantify and investigate the FLASH effect as a function of physiological and radiation parameters at microscopic scale.

METHODS

We considered time-dependent O2 supply and consumption in a 3D cylindrical geometry, incorporating blood flow linking the O2 concentration ([O2]) in the capillary to that within the tissue through the Hill equation, radial and axial diffusion of O2, metabolic and zero-order radiolytic O2 consumption, and a pulsed radiation structure. Time-evolved distributions of [O2] were obtained by numerically solving perfusion-diffusion equations. The model enables the computation of dynamic O2 distribution and the relative change of OER (δROD) under various physiological and radiation conditions in vivo.

RESULTS

Initial [O2] level and the subsequent changes during irradiation determined δROD distribution, which strongly depends on physiological parameters, i.e., intercapillary spacing, ultimately determining the tissue area with enhanced radioresistance. We observed that the δROD/FLASH effect is affected by and sensitive to the interplay effect among physiological and radiation parameters. It renders that the FLASH effect can be tissue environment dependent. The saturation of FLASH normal tissue protection upon dose and dose rate was shown. Beyond ∼60 Gy/s, no significant decrease in radiosensitivity within tissue region was observed. In turn, for a given dose rate, the change of radiosensitivity became saturated after a certain dose level. Pulse structures with the same dose and instantaneous dose rate but with different delivery times were shown to have distinguishable δROD thus tissue sparing, suggesting the average dose rate could be a metric assessing the FLASH effect and demonstrating the capability of our model to support experimental findings.

CONCLUSION

On a macroscopic scale, the modeling results align with the experimental findings in terms of dose and dose rate thresholds, and it also indicates that pulse structure can vary the FLASH effect. At the microscopic level, this model enables us to examine the spatially resolved FLASH effect based on physiological and irradiation parameters. Our model thus provides a complementary approach to experimental methods for understanding the underlying mechanism of FLASH radiotherapy. Our results show that physiological conditions can potentially determine the FLASH efficacy in tissue protection. The FLASH effect may be observed under optimal combination of physiological parameters, not limited to radiation conditions alone.

Technical note: Dosimetry and FLASH potential of UHDR proton PBS for small lung tumors: Bragg-peak-based delivery versus transmission beam and IMPT

BACKGROUND

High-energy transmission beams (TBs) are currently the main delivery method for proton pencil beam scanning ultrahigh dose-rate (UHDR) FLASH radiotherapy. TBs place the Bragg-peaks behind the target, outside the patient, making delivery practical and achievement of high dose-rates more likely. However, they lead to higher integral dose compared to conventional intensity-modulated proton therapy (IMPT), in which Bragg-peaks are placed within the tumor. It is hypothesized that, when energy changes are not required and high beam currents are possible, Bragg-peak-based beams can not only achieve more conformal dose distributions than TBs, but also have more FLASH-potential.

PURPOSE

This works aims to verify this hypothesis by taking three different Bragg-peak-based delivery techniques and comparing them with TB and IMPT-plans in terms of dosimetry and FLASH-potential for single-fraction lung stereotactic body radiotherapy (SBRT).

METHODS

For a peripherally located lung target of various sizes, five different proton plans were made using "matRad" and inhouse-developed algorithms for spot/energy-layer/beam reduction and minimum monitor unit maximization: (1) IMPT-plan, reference for dosimetry, (2) TB-plan, reference for FLASH-amount, (3) pristine Bragg-peak plan (non-depth-modulated Bragg-peaks), (4) Bragg-peak plan using generic ridge filter, and (5) Bragg-peak plan using 3D range-modulated ridge filter.

RESULTS

Bragg-peak-based plans are able to achieve sufficient plan quality and high dose-rates. IMPT-plans resulted in lowest OAR-dose and integral dose (also after a FLASH sparing-effect of 30%) compared to both TB-plans and Bragg-peak-based plans. Bragg-peak-based plans vary only slightly between themselves and generally achieve lower integral dose than TB-plans. However, TB-plans nearly always resulted in lower mean lung dose than Bragg-peak-based plans and due to a higher amount of FLASH-dose for TB-plans, this difference increased after including a FLASH sparing-effect.

CONCLUSION

This work indicates that there is no benefit in using Bragg-peak-based beams instead of TBs for peripherally located, UHDR stereotactic lung radiotherapy, if lung dose is the priority.

In vivo measurements of change in tissue oxygen level during irradiation reveal novel dose rate dependence

BACKGROUND AND PURPOSE

This study aimed to investigate the radiochemical oxygen depletion (ROD) in vivo by directly measuring oxygen levels in various mouse tissues during ultra-high dose rate (UHDR) irradiation at clinically relevant doses and dose rates.

MATERIALS AND METHODS

Mice bearing subcutaneous human glioblastoma (U-87 MG) tumors were used for tumor and normal tissue (skin, muscle, brain) measurements. An oxygen-sensitive phosphorescent probe (Oxyphor PtG4) was injected into the tissues, and oxygen levels were monitored using a fiberoptic phosphorometer during UHDR irradiation with a 6 MeV electron linear accelerator (LINAC). Dose escalation experiments (10-40 Gy) were performed at a dose rate of 1300 Gy/s, and dose rate escalation experiments were conducted at a fixed dose of 40 Gy with dose rates ranging from 2 to 101 Gy/s.

RESULTS

Radiation-induced change in tissue oxygenation (ΔpO2) increased linearly with dose and correlated with baseline tissue oxygenation levels in the range of 0 - 30 mmHg. At higher baseline tissue oxygenation levels, such as those observed in muscle and brain, there was no corresponding increase in ΔpO2. When we modulated dose rate, ΔpO2 increased steeply up to ∼ 20 Gy/s and plateaued thereafter. The relationship between ΔpO2 and dose rate showcases the interplay between ROD and reoxygenation.

CONCLUSION

While UHDR irradiation induces measurable oxygen depletion in tissues, the observed changes in oxygenation levels do not support the hypothesis that ROD-induced radioresistance is responsible for the FLASH tissue-sparing effect at clinically relevant doses and dose rates. These findings highlight the need for further investigation into alternative mechanisms underlying the FLASH effect.

Radiation safety of ultra-high dose rate electron accelerators for FLASH radiotherapy

BACKGROUND

An ultra-high dose rate (UHDR) electron accelerator for FLASH radiotherapy (RT) produces very intense bremsstrahlung by the interaction of the electron beam with objects both inside and outside of the accelerator. The bremsstrahlung dose per pulse is typically 1-2 orders of magnitude larger than that of conventional RT x-ray treatment of the same energy, and for electron energies above 10 MeV, the bremsstrahlung produces substantially more induced radioactivity outside the accelerator than for conventional RT. Therefore, a thorough radiation safety assessment is mandatory prior to the operation of a UHDR electron accelerator.

PURPOSE

To evaluate the radiation safety of a prototype FLASH-enabled Varian TrueBeam accelerator and to develop a general framework for assessment of all key radiation safety properties of a UHDR electron accelerator for FLASH RT.

METHODS

Production of bremsstrahlung and induced radioactivity by a UHDR electron accelerator is modeled by various analytical methods. The analytical modeling is compared with National Institute of Standards and Technology (NIST) bremsstrahlung yield data as well as measurements of primary bremsstrahlung outside the bunker and induced radioactivity of irradiated thick targets for a FLASH-enabled 16 MeV Varian TrueBeam electron accelerator. In addition, the analytical modeling is complemented by measurements of secondary bremsstrahlung inside/outside the bunker and neutrons at the maze entrance.

RESULTS

Calculated bremsstrahlung yields deviate maximum 8.5% from NIST data, and all measurements of primary bremsstrahlung and induced radioactivity agree with calculations, validating the analytical tools. In addition, it is found that scattering foil bremsstrahlung dominates primary bremsstrahlung and the main source of secondary bremsstrahlung is the irradiated object outside the accelerator. It follows that primary and secondary bremsstrahlung outside the bunker can be calculated using the same simple formalism as that used for conventional RT. Measured primary bremsstrahlung tenth-value layers for concrete of the simple formalism are in good agreement with NCRP and IAEA data, while measured secondary bremsstrahlung tenth-value layers for concrete are considerably lower than NCRP and IAEA data. All calculations and measurements form a general framework for assessment of all key radiation safety properties of a UHDR electron accelerator.

CONCLUSIONS

The FLASH-enabled Varian TrueBeam accelerator is safe for normal operation (max. 99 pulses per irradiation) in a bunker designed for at least 15 MV conventional x-ray treatment unless the UHDR workload is much larger than the x-ray workload. A similar finding applies to other UHDR electron accelerators. However, during beam tuning, radiation survey, or other tests with extended irradiation time, the UHDR workload may become very large, necessitating the implementation of additional safety measures.

Prediction of the treatment effect of FLASH radiotherapy with synchrotron radiation from the Circular Electron-Positron Collider (CEPC)

The Circular Electron-Positron Collider (CEPC) in China can also work as an excellent powerful synchrotron light source, which can generate high-quality synchrotron radiation. This synchrotron radiation has potential advantages in the medical field as it has a broad spectrum, with energies ranging from visible light to X-rays used in conventional radiotherapy, up to several megaelectronvolts. FLASH radiotherapy is one of the most advanced radiotherapy modalities. It is a radiotherapy method that uses ultra-high dose rate irradiation to achieve the treatment dose in an instant; the ultra-high dose rate used is generally greater than 40 Gy s-1, and this type of radiotherapy can protect normal tissues well. In this paper, the treatment effect of CEPC synchrotron radiation for FLASH radiotherapy was evaluated by simulation. First, a Geant4 simulation was used to build a synchrotron radiation radiotherapy beamline station, and then the dose rate that the CEPC can produce was calculated. A physicochemical model of radiotherapy response kinetics was then established, and a large number of radiotherapy experimental data were comprehensively used to fit and determine the functional relationship between the treatment effect, dose rate and dose. Finally, the macroscopic treatment effect of FLASH radiotherapy was predicted using CEPC synchrotron radiation through the dose rate and the above-mentioned functional relationship. The results show that the synchrotron radiation beam from the CEPC is one of the best beams for FLASH radiotherapy.

Dosimetry: Was and Is an Absolute Requirement for Quality Radiation Research

This review aims to trace the evolution of dosimetry, highlight its significance in the advancement of radiation research, and identify the current trends and methodologies in the field. Key historical milestones, starting with the first publications in the journal in 1954, will be synthesized before addressing contemporary practices in radiation medicine and radiobiological investigation. Finally, possibilities for future opportunities in dosimetry will be offered. The overarching goal is to emphasize the indispensability of accurate and reproducible dosimetry in enhancing the quality of radiation research and practical applications of ionizing radiation.

Secondary neutron dosimetry for conformal FLASH proton therapy

BACKGROUND

Cyclotron-based proton therapy systems utilize the highest proton energies to achieve an ultra-high dose rate (UHDR) for FLASH radiotherapy. The deep-penetrating range associated with this high energy can be modulated by inserting a uniform plate of proton-stopping material, known as a range shifter, in the beam path at the nozzle to bring the Bragg peak within the target while ensuring high proton transport efficiency for UHDR. Aluminum has been recently proposed as a range shifter material mainly due to its high compactness and its mechanical properties. A possible drawback lies in the fact that aluminum has a larger cross-section of producing secondary neutrons compared to conventional plastic range shifters. Accordingly, an increase in secondary neutron contamination was expected during the delivery of range-modulated FLASH proton therapy, potentially heightening neutron-induced carcinogenic risks to the patient.

PURPOSE

We conducted neutron dosimetry using simulations and measurements to evaluate excess dose due to neutron exposure during UHDR proton irradiation with aluminum range shifters compared to plastic range shifters.

METHODS

Monte Carlo simulations in TOPAS were performed to investigate the secondary neutron production characteristics with aluminum range shifter during 225 MeV single-spot proton irradiation. The computational results were validated against measurements with a pair of ionization chambers in an out-of-field region ( ≤ $\le$ 30 cm) and with a Proton Recoil Scintillator-Los Alamos rem meter in a far-out-of-field region (0.5-2.5 m). The assessments were repeated with solid water slabs as a surrogate for the conventional range shifter material to evaluate the impact of aluminum on neutron yield. The results were compared with the International Electrotechnical Commission (IEC) standards to evaluate the clinical acceptance of the secondary neutron yield.

RESULTS

For a range modulation up to 26 cm in water, the maximum simulated and measured values of out-of-field secondary neutron dose equivalent per therapeutic dose with aluminum range shifter were found to be( 0.57 ± 0.02 ) mSv/Gy $(0.57\pm 0.02)\ \text{mSv/Gy}$ and( 0.46 ± 0.04 ) mSv/Gy $(0.46\pm 0.04)\ \text{mSv/Gy}$ , respectively, overall higher than the solid water cases (simulation:( 0.332 ± 0.003 ) mSv/Gy $(0.332\pm 0.003)\ \text{mSv/Gy}$ ; measurement:( 0.33 ± 0.03 ) mSv/Gy $(0.33\pm 0.03)\ \text{mSv/Gy}$ ). The maximum far out-of-field secondary neutron dose equivalent was found to be (8.8 ± 0.5 $8.8 \pm 0.5$ ) μ Sv / Gy $\umu {\rm Sv/Gy}$ and (1.62 ± 0.02 $1.62 \pm 0.02$ ) μ Sv / Gy $\umu {\rm Sv/Gy}$ for the simulations and rem meter measurements, respectively, also higher than the solid water counterparts (simulation: (3.3 ± 0.3 $3.3 \pm 0.3$ ) μ Sv / Gy $\umu {\rm Sv/Gy}$ ; measurement: (0.63 ± 0.03 $0.63 \pm 0.03$ ) μ Sv / Gy $\umu {\rm Sv/Gy}$ ).

CONCLUSIONS

We conducted simulations and measurements of secondary neutron production under proton irradiation at FLASH energy with range shifters. We found that the secondary neutron yield increased when using aluminum range shifters compared to conventional materials while remaining well below the non-primary radiation limit constrained by the IEC regulations.

Infrared microspectroscopy to elucidate the underlying biomolecular mechanisms of FLASH radiotherapy

BACKGROUND

FLASH-radiotherapy (FLASH-RT) is an emerging modality that uses ultra-high dose rates of radiation to enable curative doses to the tumor while preserving normal tissue. The biological studies showed the potential of FLASH-RT to revolutionize radiotherapy cancer treatments. However, the complex biological basis of FLASH-RT is not fully known yet.

AIM

Within this context, our aim is to get deeper insights into the biomolecular mechanisms underlying FLASH-RT through Fourier Transform Infrared Microspectroscopy (FTIRM).

METHODS

C57Bl/6J female mice were whole brain irradiated at 10 Gy with the eRT6-Oriatron system. 10 Gy FLASH-RT was delivered in 1 pulse of 1.8μs and conventional irradiations at 0.1 Gy/s. Brains were sampled and prepared for analysis 24 h post-RT. FTIRM was performed at the MIRAS beamline of ALBA Synchrotron. Infrared raster scanning maps of the whole mice brain sections were collected for each sample condition. Hyperspectral imaging and Principal Component Analysis (PCA) were performed in several regions of the brain.

RESULTS

PCA results evidenced a clear separation between conventional and FLASH irradiations in the 1800-950 cm-1 region, with a significant overlap between FLASH and Control groups. An analysis of the loading plots revealed that most of the variance accounting for the separation between groups was associated to modifications in the protein backbone (Amide I). This protein degradation and/or conformational rearrangement was concomitant with nucleic acid fragmentation/condensation. Cluster separation between FLASH and conventional groups was also present in the 3000-2800 cm-1 region, being correlated with changes in the methylene and methyl group concentrations and in the lipid chain length. Specific vibrational features were detected as a function of the brain region.

CONCLUSION

This work provided new insights into the biomolecular effects involved in FLASH-RT through FTIRM. Our results showed that beyond nucleic acid investigations, one should take into account other dose-rate responsive molecules such as proteins, as they might be key to understand FLASH effect.

FLASH Radiotherapy: Mechanisms of Biological Effects and the Therapeutic Potential in Cancer

Radiotherapy is an important treatment for many unresectable advanced malignant tumors, and radiotherapy-associated inflammatory reactions to radiation and other toxic side effects are significant reasons which reduce the quality of life and survival of patients. FLASH-radiotherapy (FLASH-RT), a prominent topic in recent radiation therapy research, is an ultra-high dose rate treatment known for significantly reducing therapy time while effectively targeting tumors. This approach minimizes radiation side effects on at-risk organs and maximally protects surrounding healthy tissues. Despite decades of preclinical exploration and some notable achievements, the mechanisms behind FLASH effects remain debated. Standardization is still required for the type of FLASH-RT rays and dose patterns. This review addresses the current state of FLASH-RT research, summarizing the biological mechanisms behind the FLASH effect. Additionally, it examines the impact of FLASH-RT on immune cells, cytokines, and the tumor immune microenvironment. Lastly, this review will discuss beam characteristics, potential clinical applications, and the relevance and applicability of FLASH-RT in treating advanced cancers.

Pixelwise Gradient Model for Image Fusion (PGMIF): a multi-sequence magnetic resonance imaging (MRI) fusion model for tumor contrast enhancement of nasopharyngeal carcinoma

Background

Different image modalities capture different aspects of a patient. It is desirable to produce images that capture all such features in a single image. This research investigates the potential of multi-modal image fusion method to enhance magnetic resonance imaging (MRI) tumor contrast and its consistency across different patients, which can capture both the anatomical structures and tumor contrast clearly in one image, making MRI-based target delineation more accurate and efficient.

Methods

T1-weighted (T1-w) and T2-weighted (T2-w) magnetic resonance (MR) images from 80 nasopharyngeal carcinoma (NPC) patients were used. A novel image fusion method, Pixelwise Gradient Model for Image Fusion (PGMIF), which is based on the pixelwise gradient to capture the shape and a generative adversarial network (GAN) term to capture the image contrast, was introduced. PGMIF is compared with several popular fusion methods. The performance of fusion methods was quantified using two metrics: the tumor contrast-to-noise ratio (CNR), which aims to measure the contrast of the edges, and a Generalized Sobel Operator Analysis, which aims to measure the sharpness of edge.

Results

The PGMIF method yielded the highest CNR [median (mdn) =1.208, interquartile range (IQR) =1.175-1.381]. It was a statistically significant enhancement compared to both T1-w (mdn =1.044, IQR =0.957-1.042, P<5.60×10-4) and T2-w MR images (mdn =1.111, IQR =1.023-1.182, P<2.40×10-3), and outperformed other fusion models: Gradient Model with Maximum Comparison among Images (GMMCI) (mdn =0.967, IQR =0.795-0.982, P<5.60×10-4), Deep Learning Model with Weighted Loss (DLMWL) (mdn =0.883, IQR =0.832-0.943, P<5.60×10-4), Pixelwise Weighted Average (PWA) (mdn =0.875, IQR =0.806-0.972, P<5.60×10-4) and Maximum of Images (MoI) (mdn =0.863, IQR =0.823-0.991, P<5.60×10-4). In terms of the Generalized Sobel Operator Analysis, a measure based on Sobel operator to measure contrast enhancement, PGMIF again exhibited the highest Generalized Sobel Operator (mdn =0.594, IQR =0.579-0.607; mdn =0.692, IQR =0.651-0.718 for comparison with T1-w and T2-w images), compared to: GMMCI (mdn =0.491, IQR =0.458-0.507, P<5.60×10-4; mdn =0.495, IQR =0.487-0.533, P<5.60×10-4), DLMWL (mdn =0.292, IQR =0.248-0.317, P<5.60×10-4; mdn =0.191, IQR =0.179-0.243, P<5.60×10-4), PWA (mdn =0.423, IQR =0.383-0.455, P<5.60×10-4; mdn =0.448, IQR =0.414-0.463, P<5.60×10-4) and MoI (mdn =0.437, IQR =0.406-0.479, P<5.60×10-4; mdn =0.540, IQR =0.521-0.636, P<5.60×10-4), demonstrating superior contrast enhancement and sharpness compared to other methods.

Conclusions

Based on the tumor CNR and Generalized Sobel Operator Analysis, the proposed PGMIF method demonstrated its capability of enhancing MRI tumor contrast while keeping the anatomical structures of the input images. It holds promises for NPC tumor delineation in radiotherapy.

MRI-only based material mass density and relative stopping power estimation via deep learning for proton therapy: a preliminary study

Magnetic Resonance Imaging (MRI) is increasingly being used in treatment planning due to its superior soft tissue contrast, which is useful for tumor and soft tissue delineation compared to computed tomography (CT). However, MRI cannot directly provide mass density or relative stopping power (RSP) maps, which are required for calculating proton radiotherapy doses. Therefore, the integration of artificial intelligence (AI) into MRI-based treatment planning to estimate mass density and RSP directly from MRI has generated significant interest. A deep learning (DL) based framework was developed to establish a voxel-wise correlation between MR images and mass density as well as RSP. To facilitate the study, five tissue substitute phantoms were created, representing different tissues such as skin, muscle, adipose tissue, 45% hydroxyapatite (HA), and spongiosa bone. The composition of these phantoms was based on information from ICRP reports. Additionally, two animal tissue phantoms, simulating pig brain and liver, were prepared for DL training purposes. The phantom study involved the development of two DL models. The first model utilized clinical T1 and T2 MRI scans as input, while the second model incorporated zero echo time (ZTE) MRI scans. In the patient application study, two more DL models were trained: one using T1 and T2 MRI scans as input, and another model incorporating synthetic dual-energy computed tomography (sDECT) images to provide accurate bone tissue information. The DECT empirical model was used as a reference to evaluate the proposed models in both phantom and patient application studies. The DECT empirical model was selected as the reference for evaluating the proposed models in both phantom and patient application studies. In the phantom study, the DL model based on T1, and T2 MRI scans demonstrated higher accuracy in estimating mass density and RSP for skin, muscle, adipose tissue, brain, and liver. The mean absolute percentage errors (MAPE) were 0.42%, 0.14%, 0.19%, 0.78%, and 0.26% for mass density, and 0.30%, 0.11%, 0.16%, 0.61%, and 0.23% for RSP, respectively. The DL model incorporating ZTE MRI further improved the accuracy of mass density and RSP estimation for 45% HA and spongiosa bone, with MAPE values of 0.23% and 0.09% for mass density, and 0.19% and 0.07% for RSP, respectively. These results demonstrate the feasibility of using an MRI-only approach combined with DL methods for mass density and RSP estimation in proton therapy treatment planning. By employing this approach, it is possible to obtain the necessary information for proton radiotherapy directly from MRI scans, eliminating the need for additional imaging modalities.

Mechanisms of Action in FLASH Radiotherapy: A Comprehensive Review of Physicochemical and Biological Processes on Cancerous and Normal Cells

The advent of FLASH radiotherapy (FLASH-RT) has brought forth a paradigm shift in cancer treatment, showcasing remarkable normal cell sparing effects with ultra-high dose rates (>40 Gy/s). This review delves into the multifaceted mechanisms underpinning the efficacy of FLASH effect, examining both physicochemical and biological hypotheses in cell biophysics. The physicochemical process encompasses oxygen depletion, reactive oxygen species, and free radical recombination. In parallel, the biological process explores the FLASH effect on the immune system and on blood vessels in treatment sites such as the brain, lung, gastrointestinal tract, skin, and subcutaneous tissue. This review investigated the selective targeting of cancer cells and the modulation of the tumor microenvironment through FLASH-RT. Examining these mechanisms, we explore the implications and challenges of integrating FLASH-RT into cancer treatment. The potential to spare normal cells, boost the immune response, and modify the tumor vasculature offers new therapeutic strategies. Despite progress in understanding FLASH-RT, this review highlights knowledge gaps, emphasizing the need for further research to optimize its clinical applications. The synthesis of physicochemical and biological insights serves as a comprehensive resource for cell biology, molecular biology, and biophysics researchers and clinicians navigating the evolution of FLASH-RT in cancer therapy.

FLASH Radiotherapy: What Can FLASH's Ultra High Dose Rate Offer to the Treatment of Patients With Sarcoma?

FLASH is an emerging treatment paradigm in radiotherapy (RT) that utilizes ultra-high dose rates (UHDR; >40 Gy)/s) of radiation delivery. Developing advances in technology support the delivery of UHDR using electron and proton systems, as well as some ion beam units (eg, carbon ions), while methods to achieve UHDR with photons are under investigation. The major advantage of FLASH RT is its ability to increase the therapeutic index for RT by shifting the dose response curve for normal tissue toxicity to higher doses. Numerous preclinical studies have been conducted to date on FLASH RT for murine sarcomas, alongside the investigation of its effects on relevant normal tissues of skin, muscle, and bone. The tumor control achieved by FLASH RT of sarcoma models is indistinguishable from that attained by treatment with standard RT to the same total dose. FLASH's high dose rates are able to mitigate the severity or incidence of RT side effects on normal tissues as evaluated by endpoints ranging from functional sparing to histological damage. Large animal studies and clinical trials of canine patients show evidence of skin sparing by FLASH vs. standard RT, but also caution against delivery of high single doses with FLASH that exceed those safely applied with standard RT. Also, a human clinical trial has shown that FLASH RT can be delivered safely to bone metastasis. Thus, data to date support continued investigations of clinical translation of FLASH RT for the treatment of patients with sarcoma. Toward this purpose, hypofractionated irradiation schemes are being investigated for FLASH effects on sarcoma and relevant normal tissues.

Pencil Beam Scanning Proton Bragg Peak Conformal FLASH in Prostate Cancer Stereotactic Body Radiotherapy

Bragg peak FLASH radiotherapy (RT) uses a distal tracking method to eliminate exit doses and can achieve superior OAR sparing. This study explores the application of this novel method in stereotactic body radiotherapy prostate FLASH-RT. An in-house platform was developed to enable intensity-modulated proton therapy (IMPT) planning using a single-energy Bragg peak distal tracking method. The patients involved in the study were previously treated with proton stereotactic body radiotherapy (SBRT) using the pencil beam scanning (PBS) technique to 40 Gy in five fractions. FLASH plans were optimized using a four-beam arrangement to generate a dose distribution similar to the conventional opposing beams. All of the beams had a small angle of two degrees from the lateral direction to increase the dosimetry quality. Dose metrics were compared between the conventional PBS and the Bragg peak FLASH plans. The dose rate histogram (DRVH) and FLASH metrics of 40 Gy/s coverage (V40Gy/s) were investigated for the Bragg peak plans. There was no significant difference between the clinical and Bragg peak plans in rectum, bladder, femur heads, large bowel, and penile bulb dose metrics, except for Dmax. For the CTV, the FLASH plans resulted in a higher Dmax than the clinical plans (116.9% vs. 103.3%). For the rectum, the V40Gy/s reached 94% and 93% for 1 Gy dose thresholds in composite and single-field evaluations, respectively. Additionally, the FLASH ratio reached close to 100% after the application of the 5 Gy threshold in composite dose rate assessment. In conclusion, the Bragg peak distal tracking method can yield comparable plan quality in most OARs while preserving sufficient FLASH dose rate coverage, demonstrating that the ultra-high dose technique can be applied in prostate FLASH SBRT.

Initial experience with an electron FLASH research extension (FLEX) for the Clinac system

PURPOSE

Radiotherapy delivered at ultra-high-dose-rates (≥40 Gy/s), that is, FLASH, has the potential to effectively widen the therapeutic window and considerably improve the care of cancer patients. The underlying mechanism of the FLASH effect is not well understood, and commercial systems capable of delivering such dose rates are scarce. The purpose of this study was to perform the initial acceptance and commissioning tests of an electron FLASH research product for preclinical studies.

METHODS

A linear accelerator (Clinac 23EX) was modified to include a non-clinical FLASH research extension (the Clinac-FLEX system) by Varian, a Siemens Healthineers company (Palo Alto, CA) capable of delivering a 16 MeV electron beam with FLASH and conventional dose rates. The acceptance, commissioning, and dosimetric characterization of the FLEX system was performed using radiochromic film, optically stimulated luminescent dosimeters, and a plane-parallel ionization chamber. A radiation survey was conducted for which the shielding of the pre-existing vault was deemed sufficient.

RESULTS

The Clinac-FLEX system is capable of delivering a 16 MeV electron FLASH beam of approximately 1 Gy/pulse at isocenter and reached a maximum dose rate >3.8 Gy/pulse near the upper accessory mount on the linac gantry. The percent depth dose curves of the 16 MeV FLASH and conventional modes for the 10 × 10 cm2 applicator agreed within 0.5 mm at a range of 50% of the maximum dose. Their respective profiles agreed well in terms of flatness but deviated for field sizes >10 × 10 cm2 . The output stability of the FLASH system exhibited a dose deviation of <1%. Preliminary cell studies showed that the FLASH dose rate (180 Gy/s) had much less impact on the cell morphology of 76N breast normal cells compared to the non-FLASH dose rate (18 Gy/s), which induced large-size cells.

CONCLUSION

Our studies characterized the non-clinical Clinac-FLEX system as a viable solution to conduct FLASH research that could substantially increase access to ultra-high-dose-rate capabilities for scientists.

Pulse parameter optimizer: an efficient tool for achieving prescribed dose and dose rate with electron FLASH platforms

Purpose. Commercial electron FLASH platforms deliver ultra-high dose rate doses at discrete combinations of pulse parameters including pulse width (PW), pulse repetition frequency (PRF) and number of pulses (N), which dictate unique combinations of dose and dose rates. Additionally, collimation, source to surface distance, and airgaps also vary the dose per pulse (DPP). Currently, obtaining pulse parameters for the desired dose and dose rate is a cumbersome manual process involving creating, updating, and looking up values in large spreadsheets for every treatment configuration. This work presents a pulse parameter optimizer application to match intended dose and dose rate precisely and efficiently.Methods. Dose and dose rate calculation methods have been described for a commercial electron FLASH platform. A constrained optimization for the dose and dose rate cost function was modelled as a mixed integer problem in MATLAB (The MathWorks Inc., Version9.13.0 R2022b, Natick, Massachusetts). The beam and machine data required for the application were acquired using GafChromic film and alternating current current transformers (ACCTs). Variables for optimization included DPP for every collimator, PW and PRF measured using ACCT and airgap factors.Results. Using PW, PRF,Nand airgap factors as parameters, a software was created to optimize dose and dose rate, reaching the closest match if exact dose and dose rates are not achievable. Optimization took 20 s or less to converge to results. This software was validated for accuracy of dose calculation and precision in matching prescribed dose and dose rate.Conclusion. A pulse parameter optimization application was built for a commercial electron FLASH platform to increase efficiency in dose, dose rate, and pulse parameter prescription process. Automating this process reduces safety concerns associated with manual look up and calculation of these parameters, especially when many subjects at different doses and dose rates are to be safely managed.

FLASH Radiotherapy and the Use of Radiation Dosimeters

Radiotherapy (RT) using ultra-high dose rate (UHDR) radiation, known as FLASH RT, has shown promising results in reducing normal tissue toxicity while maintaining tumor control. However, implementing FLASH RT in clinical settings presents technical challenges, including limited depth penetration and complex treatment planning. Monte Carlo (MC) simulation is a valuable tool for dose calculation in RT and has been investigated for optimizing FLASH RT. Various MC codes, such as EGSnrc, DOSXYZnrc, and Geant4, have been used to simulate dose distributions and optimize treatment plans. Accurate dosimetry is essential for FLASH RT, and radiation detectors play a crucial role in measuring dose delivery. Solid-state detectors, including diamond detectors such as microDiamond, have demonstrated linear responses and good agreement with reference detectors in UHDR and ultra-high dose per pulse (UHDPP) ranges. Ionization chambers are commonly used for dose measurement, and advancements have been made to address their response nonlinearities at UHDPP. Studies have proposed new calculation methods and empirical models for ion recombination in ionization chambers to improve their accuracy in FLASH RT. Additionally, strip-segmented ionization chamber arrays have shown potential for the experimental measurement of dose rate distribution in proton pencil beam scanning. Radiochromic films, such as GafchromicTM EBT3, have been used for absolute dose measurement and to validate MC simulation results in high-energy X-rays, triggering the FLASH effect. These films have been utilized to characterize ionization chambers and measure off-axis and depth dose distributions in FLASH RT. In conclusion, MC simulation provides accurate dose calculation and optimization for FLASH RT, while radiation detectors, including diamond detectors, ionization chambers, and radiochromic films, offer valuable tools for dosimetry in UHDR environments. Further research is needed to refine treatment planning techniques and improve detector performance to facilitate the widespread implementation of FLASH RT, potentially revolutionizing cancer treatment.

Deep learning-based Fast Volumetric Image Generation for Image-guided Proton FLASH Radiotherapy

Objective

FLASH radiotherapy leverages ultra-high dose-rate radiation to enhance the sparing of organs at risk without compromising tumor control probability. This may allow dose escalation, toxicity mitigation, or both. To prepare for the ultra-high dose-rate delivery, we aim to develop a deep learning (DL)-based image-guide framework to enable fast volumetric image reconstruction for accurate target localization for proton FLASH beam delivery.

Approach

The proposed framework comprises four modules, including orthogonal kV x-ray projection acquisition, DL-based volumetric image generation, image quality analyses, and water equivalent thickness (WET) evaluation. We investigated volumetric image reconstruction using kV projection pairs with four different source angles. Thirty patients with lung targets were identified from an institutional database, each patient having a four-dimensional computed tomography (CT) dataset with ten respiratory phases. Leave-phase-out cross-validation was performed to investigate the DL model's robustness for each patient.

Main results

The proposed framework reconstructed patients' volumetric anatomy, including tumors and organs at risk from orthogonal x-ray projections. Considering all evaluation metrics, the kV projections with source angles of 135° and 225° yielded the optimal volumetric images. The patient-averaged mean absolute error, peak signal-to-noise ratio, structural similarity index measure, and WET error were 75±22 HU, 19±3.7 dB, 0.938±0.044, and -1.3%±4.1%.

Significance

The proposed framework has been demonstrated to reconstruct volumetric images with a high degree of accuracy using two orthogonal x-ray projections. The embedded WET module can be used to detect potential proton beam-specific patient anatomy variations. This framework can rapidly deliver volumetric images to potentially guide proton FLASH therapy treatment delivery systems.

Twenty Years of Advancements in a Radiotherapy Facility: Clinical Protocols, Technology, and Management

BACKGROUND

Hypo-fractionation can be an effective strategy to lower costs and save time, increasing patient access to advanced radiation therapy. To demonstrate this potential in practice within the context of temporal evolution, a twenty-year analysis of a representative radiation therapy facility from 2003 to 2022 was conducted. This analysis utilized comprehensive data to quantitatively evaluate the connections between advanced clinical protocols and technological improvements. The findings provide valuable insights to the management team, helping them ensure the delivery of high-quality treatments in a sustainable manner.

METHODS

Several parameters related to treatment technique, patient positioning, dose prescription, fractionation, equipment technology content, machine workload and throughput, therapy times and patients access counts were extracted from departmental database and analyzed on a yearly basis by means of linear regression.

RESULTS

Patients increased by 121 ± 6 new per year (NPY). Since 2010, the incidence of hypo-fractionation protocols grew thanks to increasing Linac technology. In seven years, both the average number of fractions and daily machine workload decreased by -0.84 ± 0.12 fractions/year and -1.61 ± 0.35 patients/year, respectively. The implementation of advanced dose delivery techniques, image guidance and high dose rate beams for high fraction doses, currently systematically used, has increased the complexity and reduced daily treatment throughput since 2010 from 40 to 32 patients per 8 h work shift (WS8). Thanks to hypo-fractionation, such an efficiency drop did not affect NPY, estimating 693 ± 28 NPY/WS8, regardless of the evaluation time. Each newly installed machine was shown to add 540 NPY, while absorbing 0.78 ± 0.04 WS8. The COVID-19 pandemic brought an overall reduction of 3.7% of patients and a reduction of 0.8 fractions/patient, to mitigate patient crowding in the department.

CONCLUSIONS

The evolution of therapy protocols towards hypo-fractionation was supported by the use of proper technology. The characteristics of this process were quantified considering time progression and organizational aspects. This strategy optimized resources while enabling broader access to advanced radiation therapy. To truly value the benefit of hypo-fractionation, a reimbursement policy should focus on the patient rather than individual treatment fractionation.

Feasibility study of hybrid inverse planning with transmission beams and single-energy spread-out Bragg peaks for proton FLASH radiotherapy

BACKGROUND

Ultra-high dose rate (FLASH) proton planning with only transmission beams (TBs) has limitations in normal tissue sparing. The single-energy spread-out Bragg peaks (SESOBPs) of the FLASH dose rate have been demonstrated feasible for proton FLASH planning.

PURPOSE

To investigate the feasibility of combining TBs and SESOBPs for proton FLASH treatment.

METHODS

A hybrid inverse optimization method was developed to combine the TBs and SESOBPs (TB-SESOBP) for FLASH planning. The SESOBPs were generated field-by-field from spreading out the BPs by pre-designed general bar ridge filters (RFs) and placed at the central target by range shifters (RSs) to obtain a uniform dose within the target. The SESOBPs and TBs were fully placed field-by-field allowing automatic spot selection and weighting in the optimization process. A spot reduction strategy was conducted in the optimization process to push up the minimum MU/spot assuring the plan deliverability at beam current of 165 nA. The TB-SESOBP plans were validated in comparison with the TB only (TB-only) plans and the plans with the combination of TBs and BPs (TB-BP plans) regarding 3D dose and dose rate (dose-averaged dose rate) distributions for five lung cases. The FLASH dose rate coverage (V40Gy/s ) was evaluated in the structure volume receiving > 10% of the prescription dose.

RESULTS

Compared to the TB-only plans, the mean spinal cord D1.2cc drastically reduced by 41% (P < 0.05), the mean lung V7Gy and V7.4 Gy moderately reduced by up to 17% (P < 0.05), and the target dose homogeneity slightly increased in the TB-SESOBP plans. Comparable dose homogeneity was achieved in both TB-SESOBP and TB-BP plans. Besides, prominent improvements were achieved in lung sparing for the cases of relatively large targets by the TB-SESOBP plans compared to the TB-BP plans. The targets and the skin were fully covered with the FLASH dose rate in all three plans. For the OARs, V40Gy/s  = 100% was achieved by the TB-only plans while V40Gy/s  > 85% was obtained by the other two plans.

CONCLUSION

We have demonstrated that the hybrid TB-SESOBP planning was feasible to achieve FLASH dose rate for proton therapy. With pre-designed general bar RFs, the hybrid TB-SESOBP planning could be implemented for proton adaptive FLASH radiotherapy. As an alternative FLASH planning approach to TB-only planning, the hybrid TB-SESOBP planning has great potential in dosimetrically improving OAR sparing while maintaining high target dose homogeneity.

Characterization of 250 MeV Protons from the Varian ProBeam PBS System for FLASH Radiation Therapy

Shoot-through proton FLASH radiation therapy has been proposed where the highest energy is extracted from a cyclotron to maximize the dose rate (DR). Although our proton pencil beam scanning system can deliver 250 MeV (the highest energy), this energy is not used clinically, and as such, 250 MeV has yet to be characterized during clinical commissioning. We aim to characterize the 250-MeV proton beam from the Varian ProBeam system for FLASH and assess the usability of the clinical monitoring ionization chamber (MIC) for FLASH use. We measured the following data for beam commissioning: integral depth dose curve, spot sigma, and absolute dose. To evaluate the MIC, we measured output as a function of beam current. To characterize a 250 MeV FLASH beam, we measured (1) the central axis DR as a function of current and spot spacing and arrangement, (2) for a fixed spot spacing, the maximum field size that achieves FLASH DR (ie, > 40 Gy/s), and (3) DR reproducibility. All FLASH DR measurements were performed using an ion chamber for the absolute dose, and irradiation times were obtained from log files. We verified dose measurements using EBT-XD films and irradiation times using a fast, pixelated spectral detector. R90 and R80 from integral depth dose were 37.58 and 37.69 cm, and spot sigma at the isocenter were σx = 3.336 and σy = 3.332 mm, respectively. The absolute dose output was measured as 0.343 Gy*mm2/MU for the commissioning conditions. Output was stable for beam currents up to 15 nA and gradually increased to 12-fold for 115 nA. Dose and DR depended on beam current, spot spacing, and arrangement and could be reproduced with 6.4% and 4.2% variations, respectively. Although FLASH was achieved and the largest field size that delivers FLASH DR was determined as 35 × 35 mm2, the current MIC has DR dependence, and users should measure dose and DR independently each time for their FLASH applications.

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.