Abstract:Recent progresses in intravital imaging have enabled highly-resolved measurements of periarteriolar oxygen gradients (POGs) within the brain parenchyma. POGs are increasingly used as proxies to estimate the local baseline oxygen consumption, which is a hallmark of cell activity. However, the oxygen profile around a given arteriole arises from an interplay between oxygen consumption and delivery, not only by this arteriole but also by distant capillaries. Integrating such interactions across scales while accounting for the complex architecture of the microvascular network remains a challenge from a modelling perspective. This limits our ability to interpret the experimental oxygen maps and constitutes a key bottleneck toward the inverse determination of metabolic rates of oxygen. We revisit the problem of parenchymal oxygen transport and metabolism and introduce a simple, conservative, accurate and scalable direct numerical method going beyond canonical Krogh-type models and their associated geometrical simplifications. We focus on a two-dimensional formulation, and introduce the concepts needed to combine an operator-splitting and a Green's function approach. Oxygen concentration is decomposed into a slowly-varying contribution, discretized by Finite Volumes over a coarse cartesian grid, and a rapidly-varying contribution, approximated analytically in grid-cells surrounding each vessel. Starting with simple test cases, we thoroughly analyze the resulting errors by comparison with highly-resolved simulations of the original transport problem, showing considerable improvement of the computational-cost/accuracy balance compared to previous work. We then demonstrate the model ability to flexibly generate synthetic data reproducing the spatial dynamics of oxygen in the brain parenchyma, with sub-grid resolution. Based on these synthetic data, we show that capillaries distant from the arteriole cannot be overlooked when interpreting POGs, thus reconciling recent measurements of POGs across cortical layers with the fundamental idea that variations of vascular density within the depth of the cortex may reveal underlying differences in neuronal organization and metabolic load. The cerebral microvascular network is the logistics system that provides energy to brain cells at the right time and place. Blood flow and oxygen can now be observed dynamically in living rodents, which transformed our knowledge of the system and its role in ageing and disease. However, oxygen concentration at a given location is the result of a subtle balance between local cellular consumption, supply by neighboring vessels and their interconnections to distant ones. Thus, measurements are difficult to interpret without integrating this multi-scale component, which requires advanced computational models. This hinders our ability to bridge the gap between experiments in rodents and clinical applications in humans. In this work, we focus on oxygen transport between vessels, leveraging recent advances in multi-scale modelling and their mathematical foundations. By this way, we formulate for the first time a simple, conservative, accurate and scalable computational model for cerebral oxygen across scales, that is able to integrate the spatially heterogenous distribution of vessels. We illustrate how this model, combined to imaging, will pave the way towards better estimates of oxygen consumption, a hallmark of neural activity that cannot be directly measured.

A spatially distributed model of brain metabolism highlights the role of diffusion in brain energy metabolism

Energy Demand and Metabolite Partitioning in Spatially Lumped and Distributed Models of Neuron-Astrocyte Complex

Energy metabolism and glutamate-glutamine cycle in the brain: a stoichiometric modeling perspective

An electrodiffusive network model with multicompartmental neurons and synaptic connections

Application of a multicomponent model of convectional reaction-diffusion to description of glucose gradients in a neurovascular unit

An electrodiffusive neuron-extracellular-glia model for exploring the genesis of slow potentials in the brain

Brain Metabolism in Health and Neurodegeneration: The Interplay Among Neurons and Astrocytes

Metabolic Energy Consumption and Information Transmission of a Two-Compartment Neuron Model and Its Cortical Network

A Computational Model of Interactions Between Neuronal and Astrocytic Networks: The Role of Astrocytes in the Stability of the Neuronal Firing Rate

Analytic Models of Oxygen and Nutrient Diffusion, Metabolism Dynamics, and Architecture Optimization in Three-Dimensional Tissue Constructs with Applications and Insights in Cerebral Organoids

Analytic and Numerical Models of Oxygen and Nutrient Diffusion, Metabolism Dynamics, and Architecture Optimization in Three-Dimensional Tissue Constructs with Applications and Insights in Cerebral Organoids

Electrodiffusive model for astrocytic and neuronal ion concentration dynamics

Modeling oxygen transport in the brain: An efficient coarse-grid approach to capture perivascular gradients in the parenchyma

Using NEURON for Reaction-Diffusion Modeling of Extracellular Dynamics

Bridging Glucose Metabolism and Intrinsic Functional Organization of the Human Cortex

Spatiotemporal Neuronal Signaling by Nitric Oxide: Diffusion-Reaction Modeling and Analysis

Effect of Ionic Diffusion on Extracellular Potentials in Neural Tissue

Unveiling the role of local metabolic constraints on the structure and activity of spiking neural networks

A 3D Atlas of Functional Human Brain Energetic Connectome Based on Neuropil Distribution

Machine Learning-based Predictions of Spatial Metabolic Profiles Demonstrate the Impact of Morphology on Astrocytic Energy Metabolism

Accurate numerical simulation of electrodiffusion and water movement in brain tissue