3D computational model of oxygen depletion kinetics in brain vasculature during FLASH RT and its implications for in vivo oximetry experiments

Abstract

Purpose: Ultra-high-dose-rate irradiation, also known as FLASH, has been shown to improve the therapeutic ratio of radiation therapy (RT). The mechanism behind this effect has been partially explained by the radiochemical oxygen depletion (ROD) hypothesis, which attributes the protection of the normal tissue to the induction of transient hypoxia by ROD. To better understand the contribution of oxygen to the FLASH effect, it is necessary to measure oxygen (O₂ ) in vivo during FLASH irradiation. This study's goal is to determine the temporal resolution required to accurately measure the rapidly changing oxygen concentration immediately after FLASH irradiation.

Methods: We conducted a computational simulation of oxygen dynamics using a real vascular model that was constructed from a public fluorescence microscopy dataset. The dynamic distribution of oxygen tension (p_o2 ) during and after FLASH RT was modeled by a partial differential equation (PDE) considering oxygen diffusion, metabolism, and ROD. The underestimation of ROD due to oxygen recovery was evaluated assuming either complete or partial depletion, and a range of possible values for parameters such as oxygen diffusion, consumption, vascular p_o2 and vessel density.

Result: The O₂ concentration recovers rapidly after FLASH RT. Assuming a temporal resolution of 0.5 s, the estimated ROD is only 50.7% and 36.7% of its actual value in cases of partial and complete depletion, respectively. Additionally, the underestimation of ROD is highly dependent on the vascular density. To estimate ROD rate with 90% accuracy, temporal resolution on the order of milliseconds is required considering the uncertainty in parameters involved, especially, the diverse vascular density of the tissue.

Conclusion: The rapid recovery of O₂ poses a great challenge for in vivo ROD measurements during FLASH RT. Temporal resolution on the order of milliseconds is recommended for ROD measurements in the normal tissue. Further work is warranted to investigate whether the same requirements apply to tumors, given their irregular vasculature.

Keywords: FLASH radiotherapy; computational model; in vivo oximetry; partial differential equations; radiochemical oxygen depletion.