Brief neural stimulation results in a stereotypical pattern of vascular and metabolic response that is the basis for popular brainimaging methods such as functional magnetic resonance imagine. However, the mechanisms of transient oxygen transport and its coupling to cerebral blood flow (CBF) and oxygen metabolism (CMRO 2 ) are poorly understood. Recent experiments show that brief stimulation produces prompt arterial vasodilation rather than venous vasodilation. This work provides a neurovascular response model for brief stimulation based on transient arterial effects using one-dimensional convection-diffusion transport. Hemoglobin oxygen dissociation is included to enable predictions of absolute oxygen concentrations. Arterial CBF response is modeled using a lumped linear flow model, and CMRO 2 response is modeled using a gamma function. Using six parameters, the model successfully fit 161/166 measured extravascular oxygen time courses obtained for brief visual stimulation in cat cerebral cortex. Results show how CBF and CMRO 2 responses compete to produce the observed features of the hemodynamic response: initial dip, hyperoxic peak, undershoot, and ringing. Predicted CBF and CMRO 2 response amplitudes are consistent with experimental measurements. This model provides a powerful framework to quantitatively interpret oxygen transport in the brain; in particular, its intravascular oxygen concentration predictions provide a new model for fMRI responses. Flow & Metabolism (2013) 33, 1429-1439 doi:10.1038/jcbfm.2013; published online 12 June 2013
Journal of Cerebral BloodKeywords: cerebral blood flow; cerebral hemodynamics; mathematical modeling; neurovascular coupling; neurovascular unit INTRODUCTION Neural activity causes local changes in cerebral oxygen consumption (CMRO 2 ), blood flow (CBF), and deoxyhemoglobin concentration. 1,2 This combination of neurometabolic and neurovascular coupling enables the use of hemodynamic imaging methods to infer brain activity. Of particular interest is the response to brief periods (few seconds) of neural activity, the hemodynamic response function (HRF). The stereotypical HRF permits the use of powerful linear analysis methods.Despite the utility of the HRF for quantifying brain activity, the mechanisms of transient oxygen transport remain poorly understood. There are at least three controversial issues. First, the role of venous volume changes is not clear, particularly in response to brief neural stimulation. [3][4][5] Second, there is disagreement about the relative role of capillaries and arterioles in oxygen transport to tissue. [6][7][8] Third, experiments in oxygen mass balance have sometimes shown that metabolic consumption is insufficient to account for global oxygen loss from the vasculature. It was therefore suggested that tissues at or adjacent to the vascular wall may consume the missing oxygen. 9,10 Thus, questions remain about oxygen transport both in steady state and as a consequence of brief stimulation.Various computational models have attempted to explain t...