A Computational Model for Oxygen Depletion Hypothesis in FLASH Effect

Experiments have reported low normal tissue toxicities during FLASH irradiation, but the mechanism has not been elaborated. Several hypotheses have been proposed to explain the mechanism. One hypothesis is oxygen depletion. We analyze the time-dependent change of oxygen concentration in the tissue to study the oxygen depletion hypothesis using a computational model. The effects of physical, chemical and physiological parameters on oxygen depletion were explored. The kinetic equation of the model is solved numerically using the finite difference method with rational boundary conditions. Results of oxygen distribution is supported by the experiments of oxygen-sensitivity electrodes and experiments on the expression and distribution of the hypoxia-inducible factors. The analysis of parameters shows that the steady-state oxygen distribution before irradiation is determined by the oxygen consumption rate of the tissue and the microvessel density. The change of oxygen concentration after irradiation has been found to follow a negative exponential function, and the time constant is mainly determined by the microvessel density. The change of oxygen during exposure increases with dose rate and tends to be saturated because of oxygen diffusion. When the dose rate is high enough, the same dose results in the same reduction of oxygen concentration regardless of dose rate. The analysis of the FLASH effect in the brain tissue based on this model does not support the explanation of the oxygen depletion hypothesis. The oxygen depletion hypothesis remains controversial because the oxygen in most normal tissues cannot be depleted to radiation resistance level by FLASH irradiation.