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.

Depth-Dependent Contributions of Various Vascular Zones to Cerebral Autoregulation and Functional; Hyperemia: An In-Silico Analysis

Regional and Depth-Dependence of Cortical Blood-Flow Assessed with High-Resolution Arterial Spin Labeling (ASL)

Brain Microvasculature Has a Common Topology with Local Differences in Geometry That Match Metabolic Load.

A new computational model for quantifying blood flow dynamics across myogenically-active cerebral arterial networks

Influence of hydrodynamic and functional nonlinearities of blood flow in the cerebral vasculature on cerebral perfusion and autoregulation pressure reserve

The post-arteriole transitional zone: a specialized capillary region that regulates blood flow within the CNS microvasculature

The conducted vascular response as a mediator of hypercapnic cerebrovascular reactivity: A modelling study

Arterial vasodilation drives convective fluid flow in the brain: a poroelastic model

Signal processing methods for the analysis of cerebral blood flow and metabolism

A Theoretical Model of Cerebral Hemodynamics: Application to the Study of Arteriovenous Malformations

The contribution of the vascular architecture and cerebrovascular reactivity to the BOLD signal formation across cortical depth

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

Cerebral hemodynamics in acute stroke: pathophysiology and clinical implications

Development of a globally optimised model of the cerebral arteries

Modeling the mechanisms of non-neurogenic dynamic cerebral autoregulation

Vascular synthesis based on hemodynamic efficiency principle recapitulates measured cerebral circulation properties in the human brain

Integration of Clinical, Biological, and Computational Perspectives to Support Cerebral Autoregulatory Informed Clinical Decision Making Decomposing Cerebral Autoregulation using Mechanistic Timescales to Support Clinical Decision-Making

Precapillary sphincters and pericytes at first-order capillaries as key regulators for brain capillary perfusion

A multi-domain model for microcirculation in optic nerve: blood flow and oxygen transport

A fully synthetic three-dimensional human cerebrovascular model based on histological characteristics to investigate the hemodynamic fingerprint of the layer BOLD fMRI signal formation

Modelling the Role of Nitric Oxide in Cerebral Autoregulation