Computational Modeling and Imaging of the Intracellular Oxygen Gradient

Computational modeling can provide a mechanistic and quantitative framework for describing intracellular spatial heterogeneity of solutes such as oxygen partial pressure (pO₂). This study develops and evaluates a finite-element model of oxygen-consuming mitochondrial bioenergetics using the COMSOL Multiphysics program. The model derives steady-state oxygen (O₂) distributions from Fickian diffusion and Michaelis-Menten consumption kinetics in the mitochondria and cytoplasm. Intrinsic model parameters such as diffusivity and maximum consumption rate were estimated from previously published values for isolated and intact mitochondria. The model was compared with experimental data collected for the intracellular and mitochondrial pO₂ levels in human cervical cancer cells (HeLa) in different respiratory states and under different levels of imposed pO₂. Experimental pO₂ gradients were measured using lifetime imaging of a Förster resonance energy transfer (FRET)-based O₂ sensor, Myoglobin-mCherry, which offers in situ real-time and noninvasive measurements of subcellular pO₂ in living cells. On the basis of these results, the model qualitatively predicted (1) the integrated experimental data from mitochondria under diverse experimental conditions, and (2) the impact of changes in one or more mitochondrial processes on overall bioenergetics.