The Effects of Capillary Transit Time Heterogeneity (CTH) on the Cerebral Uptake of Glucose and Glucose Analogs: Application to FDG and Comparison to Oxygen Uptake

Introduction

Normal brain function depends critically on a constant energy supply, and therefore on moment-to-moment regulation of oxygen and glucose availability in brain tissue. In mammals, under normal physiologic conditions, glucose is the major metabolic fuel in the brain (Berg et al., 2012). Glucose can be metabolized through different metabolic pathways: through glycolysis, each glucose molecule is first metabolized into two molecules of pyruvate, with the concomitant production of two molecules of ATP. Then, in the brain, pyruvate can be converted to lactate in the absence of oxygen (anaerobic glycolysis), or completely oxidized to CO₂ under aerobic conditions, so-called oxidative phosphorylation, generating much more ATP (30 molecules per glucose molecule). In the normal adult brain, studies in many laboratories have established the overall stoichiometry of oxygen and glucose utilization. OGI is approximately equal to 5.5 and is therefore close to the theoretical maximum of 6 which corresponds to the complete oxidation of glucose (Edvinsson and Krause, 2002). This has led to the long held thesis that brain glucose metabolism is mainly oxidative. However, this idea appears to be contradicted by a phenomenon called ‘aerobic glycolysis’ which is the disproportionately higher utilization of glucose than O₂ in the normoxic working brain. This phenomenon suggests that lactate is produced although oxygen level are seemingly sufficient to support oxidative phosphorylation (Edvinsson and Krause, 2002). Some studies indicate that aerobic glycolysis is linked to functional activation. For example, during sensory stimulation or mental tasks in human, subjects have been reported to evoke 30–50% increases in blood flow and CMR_glc with little or no change in CMRO₂ when measured by PET (Fox et al., 1988) or from arterio-venous metabolite differences (Madsen et al., 1995). Traditionally, lactic acid production is believed to be related to a lack of oxygen. Indeed, it takes place in particular in skeletal muscles when energy needs outpace the ability to transport oxygen and in solid cancer tumors, which are known to grow more rapidly than the blood vessels nourish them, and therefore to experience hypoxia. Under these conditions, glycolysis and subsequent lactic acid fermentation becomes the primary source of ATP (Berg et al., 2012). In the brain, it is unclear whether lactate production occurs in conjunction with local hypoxia. Hypotheses have been formulated to provide a deeper understanding of aerobic glycolysis. For example, the astrocyte-neuron lactate shuttle (ANLS) hypothesis (Pellerin and Magistretti, 2003) proposes that an enhancement of aerobic glycolysis occurs in astrocytes in response to neuronal activation. The ANLS hypothesis predicts a reduction in the molar ratio of oxygen to glucose consumption during activation, but the proposed compartmentalization of glucose metabolism among cell types remains controversial (Chih and Roberts, 2003; Hertz, 2004). Moreover, many aspects of aerobic glycolysis remain poorly understood, including its dependency on stimulus type, duration, and magnitude (Edvinsson and Krause, 2002).

Biophysical models of oxygen delivery (Jespersen and Østergaard, 2012; Angleys et al., 2015; Rasmussen et al., 2015) suggest that the redistribution of blood flow across the capillary bed, as indexed by the extent of capillary transit time heterogeneity (CTH), affects the effective permeability surface area of the capillary bed, and hence the extraction efficacy of freely diffusible molecules. Unlike oxygen, however, glucose does not diffuse freely across the capillary membrane, and the extent to which CTH affects glucose delivery is therefore less clear. Glucose extraction is instead mediated by glucose transporters, namely the glucose transporter GLUT-1 at the blood brain barrier (BBB). GLUT-1 transporters have a maximal rate of operation, and facilitated diffusion of glucose can therefore become saturated in capillaries which support a sufficiently high rate of glucose delivery.

To examine whether the biophysical properties of oxygen and glucose extraction, respectively, introduces a need for brain parenchyma to utilize aerobic glycolysis under certain conditions, we set out to develop a model that infers glucose extraction and consumption from cerebral blood flow (CBF) as well as CTH. We then combine our predictions with those yielded by a comprehensive model of oxygen extraction (Angleys et al., 2015). Using data from in vivo rat studies, we predict the extent to which oxygen extraction capacity varies relative to that of glucose, especially among physiological conditions, as characterized by their capillary mean transit time (MTT) and CTH.

To compare our model prediction to in vivo measurements, we partly rely on measurements of glucose analogs uptake rather than native glucose. We therefore extend the model to predict FDG uptake, as well as the lumped constant (LC), which, in steady state, equals the ratio of glucose tracer to native glucose extraction and hence allows to relate the concentration of glucose trapped in the tissue to CMR_glc. Our calculation to determine LC is analog to, but differs from Holden and colleagues' (Holden et al., 1991) in that the effects of CTH are taken into account and glucose transport across the BBB is described differently (see Section Methods). After quantifying the extent to which LC changes between physiological states, we examine whether this re-evaluation of LC can explain the apparent discrepancy between PET and NMRS measurement of CMR_glc in the literature. Finally, we discuss possible clinical applications of our model findings.

Methods

Assumptions to Describe Glucose Transport across the BBB

In this study, we employ reversible, symmetric Michaelis-Menten kinetics to describe the transport of glucose across the capillary membrane. We treat the endothelium as a single membrane and thus neglect the endothelial compartment. In fact, when employing Michaelis-Menten kinetics, it is possible to show that treating the double membrane as a whole is mathematically almost equivalent to considering two identical membranes, each with identical Michaelis-Menten parameter K_T and twice the maximum transport capacity v_{max_t} (exact equivalence is obtained when considering non-reversible Michaelis-Menten). Treating the two membranes as a single barrier is also supported experimentally, for example by the reexamination of non-reversible Michaelis-Menten kinetics across a double-membrane system carried out by Gjedde and Christensen (1984). We discuss the appropriateness of reversible Michaelis-Menten kinetics further below.

Aim of the Model and General Procedure

Figure 1 outlines our procedure for computing glucose extraction fraction (GEF) and CMR_glc. Table 1 summarizes the different scientific questions asked in this study and the conditions under which the model is applied to answer each question. Note that the model is fully specified once the different input parameters in Table 1 are known. Table 2 summarizes the variables and quantities used in this computation.

FIGURE 1

Figure 1. Schematic illustrating the procedure for computing GEF and CMR_glc, given MTT, CTH and a (target) mean glucose concentration in the tissue ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$. In the first step, the true value of the effective maximum rate v_{max_m} is not known. Assuming a value for v_{max_m}, the mean concentration in the tissue ${\overline{\text{C}}}_{\text{t}}$ is compared to ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$. After each iteration, v_{max_m} is adjusted in a direction until convergence of ${\overline{\text{C}}}_{\text{t}}$ to ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$. To achieve this in practice, to reduce the computation time, the mean glucose concentration ${\overline{\text{C}}}_{\text{t}}$ is first computed over a grid of values v_{max_m}. This function is then interpolated to get ${\overline{\text{C}}}_{\text{t}}$ as a function of v_{max_m}. In the second step, the objective function $\left|{\overline{\text{C}}}_{\text{t}}\left({\text{v}}_{\text{max}\_\text{m}}\right)-{\overline{\text{C}}}_{\text{t}}^{\text{t}}\right|$ is minimized until its value is lower than 10⁻⁹. Having determined the specific value of v_{max_m}, we use the previously interpolated C_t(τ, v_{max_m}) and Q(τ, C_t(τ, v_{max_m})) to determine C_t (τ) and Q(τ) for any transit time τ. GEF, glucose extraction fraction; MTT, mean transit time; CTH, capillary transit time heterogeneity.

TABLE 1

Table 1. Scientific questions raised in this study and corresponding conditions under which the model is applied.

TABLE 2

Table 2. Terminology and parameter values for glucose extraction model.

Our model aims to determine CMR_glc in steady state from MTT, CTH, and the mean glucose concentration in the tissue ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$, where the superscript t stands for “target value.” CMR_glc is derived directly from GEF using the formula CMR_glc = C_A · CBF · GEF (Fick's principle) where C_A is the arterial glucose concentration, and from the central volume theorem which relates CBF to the MTT and the capillary volume V_cap through the relation CBF = V_cap/MTT.

To compute the mean value of any function over the capillary network, we sum the contribution of the function for every capillary, weighted by the transit time distribution h(τ). In particular, GEF corresponds simply to the mean of the single capillary glucose extraction fraction Q:

$\begin{array}{ll}\text{GEF }\left(\text{MTT},\text{CTH}\right)=\overline{\text{Q}}={\displaystyle {\int }_0^{+\infty }\text{d}\tau }·\text{Q}\left(\tau \right)·\text{h}\left(\tau ,\text{MTT},\text{CTH}\right)& (1)\end{array}$

The transit time distribution h(τ) is assumed to be a gamma distribution. This choice has been discussed extensively in Angleys et al. (2015). The gamma distribution is a two parameter distribution and is therefore fully and uniquely specified through the dependence on its mean (MTT) and standard deviation (CTH).

The individual steps to determine Q as a function of τ in order to determine GEF with Equation (1) are explained in details below. To simplify the notation, we will not indicate the dependence of the different parameters on MTT and CTH, which are fixed during the entire procedure.

Derivation of the System Equations

Network Level

At the capillary network level, the mean glucose concentration in the tissue is given and equal to ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$. In this model, we assume that glucose transfer among capillaries is negligible. As a result, at the capillary level, glucose concentrations C_t are not necessarily identical around the capillaries in our tissue compartment, and glucose concentration in tissue may therefore vary at the intercapillary distance scale.

Tissue concentrations in each compartment C_t must fulfill the equation:

$\begin{array}{ll}{\overline{\text{C}}}_{\text{t}}^{\text{t}}-{\displaystyle \underset{0}{\overset{+\infty }{\int }}\text{d}\tau ·\tilde{\text{h}}\left(\tau \right)·{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}=0& (2)\end{array}$

where v_{max_m} is a constant which will be determined later. $\tilde{\text{h}}$ is derived from h, according to the relation:

$\begin{array}{ll}\tilde{\text{h}}\left(\tau \right)=\text{h}\left(\tau \right)·\frac{\tau }{\text{MTT}}& (3)\end{array}$

Here, $\tilde{\text{h}}\left(\tau \right)d\tau$ represents a volume fraction of capillaries, as opposed to a fraction of the flow (implicit for h). Please see the Supplementary Material for more details concerning the derivation of this distribution.

Capillary Scale

Equilibrium concentration in the tissue compartment

At the capillary scale, in steady state, there is no glucose accumulation in the tissue compartment and the net rate of glucose uptake from the plasma equals the rate M at which glucose is phosphorylated by hexokinase:

$\begin{array}{ll}\frac{{\text{C}}_{\text{A}}·\text{Q}\left(\tau ,{\text{C}}_{\text{t}}\right)}{\tau }·{\text{V}}_{\text{cap}}-\text{M}\left(\tau ,{\text{C}}_{\text{t}},{\text{v}}_{\max \_\text{m}}\right)=0& (4)\end{array}$

where C_A · Q(τ, C_t) · V_cap/τ represents the net flux of glucose across a single capillary membrane, with C_A being the arterial glucose concentration in plasma, Q the extraction fraction for a single capillary and the capillary transit time. V_cap is assumed to be constant and equal to 1.4%.

We assume that M is governed by Michaelis-Menten kinetics (Michaelis and Menten, 1913):

$\begin{array}{ll}\text{M }\left(\tau ,{\text{C}}_{\text{t}},{\text{v}}_{\max \_\text{m}}\right)={\text{v}}_{\max \_\text{m}}·\frac{\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}{{\text{K}}_{\text{M}}+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}& (5)\end{array}$

where K_M is a parameter such that V_d · K_M is the concentration at which metabolism equals v_{max_m}/2. v_{max_m} is the effective maximum rate at which hexokinase can metabolize glucose to glucose-6-phosphate.

Hence, Equation (4) can be rewritten:

$\begin{array}{ll}\frac{{\text{C}}_{\text{A}}·\text{Q}\left(\tau ,{\text{C}}_{\text{t}}\right)}{\tau }·{\text{V}}_{\text{cap}}-{\text{v}}_{\max \_\text{m}}·\frac{\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}{{\text{K}}_{\text{M}}+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}=0& (6)\end{array}$

Glucose extraction fraction from the plasma

In this section, we detail how we derive the equations to express the glucose extraction fraction for a single capillary as a function of the transit time and of the concentration gradient. Glucose is considered in two compartments: plasma and extravascular tissue. Glucose is transported across the capillary membrane by facilitated diffusion via the glucose transporter GLUT-1, and the unidirectional flux of glucose across the BBB is assumed to be governed by reversible Michaelis-Menten kinetics (Cunningham, 1986). We assume that within the capillary, axial diffusion can be neglected compared to advective transport. Considering only steady-state conditions and choosing the z-axis to be oriented along the capillary, we let C_p(z) and C_t denote plasma and tissue glucose concentration, respectively. The system is then described by the following differential equation:

$\begin{array}{ll}\frac{\partial {\text{C}}_{\text{p}}\left(\text{x},\tau ,{\text{C}}_{\text{t}}\right)}{\partial \text{x}}=-\frac{\tau }{{\text{V}}_{\text{cap}}}·{\text{v}}_{\max \_\text{t}}.\frac{{\text{C}}_{\text{p}}\left(\text{x},\tau ,{\text{C}}_{\text{t}}\right)-\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}{{\text{K}}_{\text{T}}+{\text{C}}_{\text{p}}\left(\text{x},\tau ,{\text{C}}_{\text{t}}\right)+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}& (7)\end{array}$

Here,

$-\frac{\tau }{{\text{V}}_{\text{cap}}}·{\text{v}}_{\max \_\text{t}}.\frac{{\text{C}}_{\text{p}}\left(\text{x},\tau ,{\text{C}}_{\text{t}}\right)}{{\text{K}}_{\text{T}}+{\text{C}}_{\text{p}}\left(\text{x},\tau ,{\text{C}}_{\text{t}}\right)+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}$

corresponds to glucose efflux from the plasma to the tissue, while:

$\frac{\tau }{{\text{V}}_{\text{cap}}}·{\text{v}}_{\max \_\text{t}}.\frac{\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}{{\text{K}}_{\text{T}}+{\text{C}}_{\text{p}}\left(\text{x},\tau ,{\text{C}}_{\text{t}}\right)+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}$

corresponds to glucose influx, back from the tissue to the plasma.

In equation (7), x ∈ [0;1] is a normalized axial coordinate, i.e., x = z/L, with L being the capillary length; V_d is the physical distribution space of glucose in the brain, equal to 0.77 (Lund-Andersen, 1979). Throughout this study, the two Michaelis-Menten parameters, namely v_{max_t}, which corresponds to the maximum rate at which glucose can be transported unidirectionally across the BBB, and K_T, have been assigned values as explained in the Section Calibration of the Model Parameters. The plasma glucose concentration at the point x = 0 is assumed to be equal to glucose arterial plasma concentration C_A. The glucose extraction fraction for a single capillary is defined by the ratio

$\begin{array}{ll}\text{Q}\left(\tau ,{\text{C}}_{\text{t}}\right)=\frac{{\text{C}}_{\text{p}}\left(0,\tau ,{\text{C}}_{\text{t}}\right)-{\text{C}}_{\text{p}}\left(1,\tau ,{\text{C}}_{\text{t}}\right)}{{\text{C}}_{\text{p}}\left(0,\tau ,{\text{C}}_{\text{t}}\right)}& (8)\end{array}$

and depends on the transit time τ.

Summarizing, we must solve the following system of coupled equations:

$\begin{array}{l}\{\begin{array}{l}\text{(a)}{\overline{\text{C}}}_{\text{t}}^{\text{t}}-{\displaystyle \underset{0}{\overset{+\infty }{\int }}\tilde{\text{h}}\left(\tau \right)·{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}·d\tau =0\\ \text{(b)}\frac{-{\text{v}}_{\max \_\text{m}}·\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}{{\text{K}}_{\text{M}}+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}+\frac{{\text{C}}_{\text{A}}·Q\left(\tau ,{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)\right)}{\tau }·\\ {\text{V}}_{\text{cap}}=0\text{,}\\ \text{with}\text{Q}\left(\tau ,{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)\right)=1-\frac{{\text{C}}_{\text{p}}\left(1,\tau ,{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)\right)}{{\text{C}}_{\text{p}}\left(0,\tau ,{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)\right)}\\ \text{(c)}\frac{\partial {\text{C}}_{\text{p}}\left(\text{x},\tau ,{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)\right)}{\partial x}=-\frac{\tau }{{\text{V}}_{\text{cap}}}·\\ {\text{v}}_{\max \_\text{t}}.\frac{{\text{C}}_{\text{p}}\left(\text{x},\tau ,{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)\right)-\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}{{\text{K}}_{\text{T}}+{\text{C}}_{\text{p}}\left(\text{x},\tau ,{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)\right)+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}\text{,}\\ \text{with}{\text{C}}_{\text{p}}\left(0,\tau ,{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)\right)={\text{C}}_{\text{A}}& (9)\end{array}\end{array}$

for C_p, C_t and v_{max_m}, for any value of the transit time τ and relevant values of ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$. The computation is performed numerically in several steps, as no analytical solution exists for this system.

Solving the Equation System

In this section, we detail the steps to solve system (9). Briefly, in the first step, we will use Equation (9)c to express Q as a function of τ and C_t. In the second step, we will use Equation (9)b to express C_t as a function of τ and v_{max_m}. Finally, we use Equation (9)a to determine the value of v_{max_m} to get explicitly Q as a function of τ.

In the first step, we solve Equation (9)c independently over a grid of values (τ, C_t). We compute the corresponding Q function on the same grid:

$\begin{array}{ll}\text{Q}\left(\tau ,{\text{C}}_{\text{t}}\right)=1-\frac{{\text{C}}_{\text{p}}\left(1,\tau ,{\text{C}}_{\text{t}}\right)}{{\text{C}}_{\text{p}}\left(0,\tau ,{\text{C}}_{\text{t}}\right)}& (10)\end{array}$

This function is then appropriately interpolated to get sufficiently high resolution of Q(τ, C_t) while minimizing the amount of numerical computation.

In the second step, we numerically solve the equation:

$\begin{array}{ll}-{\text{v}}_{\max \_\text{m}}·\frac{\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}{{\text{K}}_{\text{M}}+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}+\frac{{\text{C}}_{\text{A}}·\text{Q}\left(\tau ,{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)\right)}{\tau }·{\text{V}}_{\text{cap}}=0& (11)\end{array}$

for C_t(τ, v_{max_m}) over an array of values v_{max_m}, and for relevant values of τ, to obtain C_t (τ, v_{max_m}) in steady-state on the same grid (τ, v_{max_m}). This function is then appropriately interpolated.

The last step consists in solving numerically the equation:

$\begin{array}{ll}{\overline{\text{C}}}_{\text{t}}^{\text{t}}-{\displaystyle \underset{0}{\overset{+\infty }{\int }}\text{d}\tau ·\tilde{\text{h}}\left(\tau \right)·{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}=0& (12)\end{array}$

in order to determine the v_{max_m} which fulfills this relation.

Having determined v_{max_m}, we use the previously interpolated C_t (τ, v_{max_m}) and Q(τ, C_t (τ, v_{max_m})) to determine C_t (τ) and Q(τ) for any transit time.

Predicted Changes in OGI during Functional Activation Based on Experimental Data from the Literature

We combine CMR_glc values with CMRO₂ values obtained with our previous model of oxygen extraction (Angleys et al., 2015) to predict the extent to which the OGI changes as a function of MTT, CTH and physiological state. Transit time characteristics in baseline and during activation are obtained from in vivo rat studies by Schulte et al. (2003) and Stefanovic et al. (2008). They appear as symbols on CMR_glc and CMRO₂ maps of Figure 4 and are listed in Table 2 in Angleys et al. (2015). We predict and compare CMRO₂ and CMR_glc in the different physiological states, and evaluate their dependence on the transit time distribution. The model of oxygen extraction is used under the assumption that the maximum metabolic rate of oxygen v_max increases by 10% from baseline condition to stimulation. See discussion in Angleys et al. (2015). Note, however, that in this latter model, CMRO₂ itself is predicted to increase by about 20% between baseline and stimulation.

Calibration of the Model Parameters

The model parameters K_M, v_{max_t}, K_T, v_{max_m} as well as the input function ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ must be fixed based on literature values and realistic assumptions. K_M was inferred from literature reports, and we set its value to 5 μmol/100 mL_brain (McIlwain and Bachelard, 1985). We tested our model using two different sets of Michaelis-Menten parameters to describe glucose transport. One parameter set, S_R, was obtained by calibrating our model to two studies conducted in rats. The first (Cunningham et al., 1986) utilized an intravenous infusion technique in rats and the second study (Choi et al., 2001) involves NMRS in rats. While CMR_glc is reported in the first study, we assume that it is equal to 45μmol/100mL/min (Choi and Gruetter, 2012) in the second, where α-chloralose-anesthetized rats were used. We set v_{max_t} to 136μmol/100mL/min, the mean of the value reported in the two studies.

The second parameter set, S_H, was obtained by calibrating our model to two studies involving NMRS in human (Gruetter et al., 1998; Seaquist et al., 2001). We assumed CMR_glc = 30μmol/100mL/min in these two studies and set v_{max_t} to 68μmol/100mL/min.

The values reported for K_T vary considerably between studies, even when identical protocols are used (see for example, Choi et al., 2001; Seaquist et al., 2001, for reported values). Even negative values have been reported (Seaquist et al., 2001), although the physical meaning of such findings remain unclear. As a result, the calibration of K_T is somewhat uncertain. When predicting CMR_glc as a function of plasma concentration (Figure 3), K_T primarily influences the plasma concentration at which the metabolism reaches its asymptotical value, that is, the slope of the curve in Figure 3. Several reports suggest that CMR_glc is relatively independent of plasma concentration for concentrations ranging between mild hypoglycemic and hyperglycemic levels (Bryan et al., 1986; Orzi et al., 1988; Suda et al., 1990; Hasselbalch et al., 2001b). Accordingly, we chose K_T equal to 50 μmol/100 mL and 150 μmol/100 mL when the model was applied to humans and rats, respectively. These values are of the same order of magnitude as those often reported in the literature.

Several reports suggest that the phosphorylation catalyzed by hexokinase is an important control step in the regulation of glucose metabolism in the brain (Clarke et al., 1989). In resting conditions, hexokinase is strongly inhibited, in particular by its product, glucose-6-phosphate, so that its operation rate at rest is only 3–10% of its maximum value (Clarke et al., 1989), suggesting that glucose metabolism can increase to accommodate higher energy demand. Accordingly, v_{max_m} is calibrated in this model to yield the glucose tissue concentration target value ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$, as described in the previous section (Solving the Equation System). Consequently, v_{max_m} is the maximum rate at which hexokinase can metabolize glucose to glucose-6-phosphate when inhibited, and is therefore the effective maximum rate of hexokinase.

In the following, the state (MTT_bl = 1.4 s, CTH_bl = 1.33 s) is taken as a reference for resting state, while the state (MTT_st = 0.81 s, CTH_st = 0.52 s) is taken as a reference for stimulation, in accordance with an experiment involving functional activation in rat (Stefanovic et al., 2008), and two studies modeling cerebral oxygen consumption using these same states as a reference for resting state and stimulation (Jespersen and Østergaard, 2012; Angleys et al., 2015). We will refer to these two states as resting state or baseline condition, and stimulation or activation, respectively.

${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ in Baseline Condition

Part of this study aims to determine the extent to which CMR_glc and LC are predicted to change as a function of the plasma concentration (see Figures 3, 6). For this part of the study, we assumed that ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ in the baseline condition varies with plasma concentration according to measurements reported in two studies, where plasma concentrations ranged from hypoglycemic to hyperglycemic levels. The first study reported direct measurements of glucose concentration in cerebral tissue in conscious rats (Dienel et al., 1991) and the second study was based on NMRS measurements in humans (Gruetter et al., 1998). Arterial and tissue concentrations were thus assumed to obey relations:

$\begin{array}{ll}{\overline{\text{C}}}_{\text{t}}^{\text{t}}=\{\begin{array}{c}0.24·{\text{C}}_{\text{A}}-72\text{if}{\text{C}}_{\text{A}}>400\mu \text{mol/100}\text{mL}\_\text{ plasma}\hfill \\ \text{24}·\frac{\text{exp}\left({\text{C}}_{\text{A}}/102\right)-1}{\text{exp}\left(400/102\right)-1}\text{ otherwise}\hfill \end{array}& (13)\end{array}$

$\begin{array}{ll}{\overline{\text{C}}}_{\text{t}}^{\text{t}}=\{\begin{array}{c}0.30·{\text{C}}_{\text{A}}-18\text{ if}{\text{C}}_{\text{A}}>250\mu \text{mol/100}\text{mL}\_\text{ plasma}\hfill \\ \text{ 57}·\frac{\text{exp}\left({\text{C}}_{\text{A}}/434\right)-1}{\text{exp}\left(250/434\right)-1}\text{ otherwise}\hfill \end{array}& (14)\end{array}$

when the model was applied to rats and humans, respectively.

In Equations (13) and (14), ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ is expressed in μmol/100 mL_brain and C_A in μmol/100 mL_plasma. Another part of the study aims to quantify the extent to which CMR_glc and the OGI change between physiological conditions, under condition of fixed plasma concentration (Figures 4, 5) and when the model is applied to human. We fixed ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ and C_A such that the molar ratio CMRO₂:CMR_glc in the resting state (MTT = 1.4 s, CTH = 1.33 s) is equal to 5.5, as reported in Madsen et al. (1995) and Hasselbalch et al. (1996b). In this state, we assume CMRO₂ to be 3.8μmol/100 mL_brain/min (i.e., 158 μmol/100 mL_brain/min), as in Jespersen and Østergaard (2012) and Angleys et al. (2015). Accordingly, we set ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ to 129 μmol/100 mL_brain and C_A = 500μmol/100 mL_plasma, which yields a corresponding CMR_glc in baseline condition equal to 29 μmol/100 mL/min.

${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ during Stimulation

Based on experimental data, we assume that ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ decreases by 30% from baseline condition to stimulation. For states intermediate between baseline and stimulation, we assume that ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ is a function of MTT only, and we determine its value by linear interpolation. We also assume that the linear relation between MTT and ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ is preserved for MTT outside of the interval [MTT_bl MTT_st]. These assumptions are discussed further below.

Glucose Tracer Kinetics and the Lumped Constant

We use our model to predict the value of the LC, which corresponds to the ratio between glucose analog (tracer) and native glucose extraction in steady state. In the following, models were calibrated to predict the LC for FDG, and “glucose analog” therefore refers to FDG throughout the manuscript unless otherwise specified. Accordingly, we determine the mean extraction fraction of glucose tracer (GEF′) over the capillary network. GEF′ is determined as GEF previously:

$\begin{array}{ll}{\text{GEF}}^{\prime }={\displaystyle {\int }_0^{+\infty }\text{d}\tau }·{\text{Q}}^{\prime }\left(\tau \right)·\text{h}\left(\tau \right)& (15)\end{array}$

where Q′ is the single capillary glucose analog extraction fraction that we must determine.

In steady state, there is no accumulation of unphosphorylated glucose analog in the extravascular compartment and the net rate of glucose analog uptake from the plasma equals the rate at which glucose is phosphorylated.

Hence, at the single capillary level:

$\begin{array}{ll}\frac{{\text{C}}_{\text{A}}\prime ·{\text{Q}}^{\prime }\left(\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)}{\tau }·{\text{V}}_{\text{cap}}-\text{M}\prime \left({\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime ,{\text{v}}_{\max \_\text{m}},{\text{v}}_{\max \_\text{m}}\prime \right)=0& (16)\end{array}$

Where M′ is the rate at which glucose tracer is phosphorylated and ${\text{C}}_{\text{A}}^{\prime }$ is the arterial glucose analog concentration.

Single Capillary Scale

The glucose analog is assumed to follow the same kinetics as native glucose, that is:

$\begin{array}{l}{\text{M}}^{\prime }\left({\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime ,{\text{v}}_{\max \_\text{m}},{\text{v}}_{\max \_\text{m}}\prime \right)\\ ={\text{v}}_{\max \_\text{m}}\prime ·\frac{{\text{C}}_{\text{t}}\prime \left(\tau ,{\text{v}}_{\max \_\text{m}},{\text{v}}_{\max \_\text{m}}\prime \right)/{\text{V}}_{\text{d}}\prime }{{\text{K}}_{\text{M}}\prime ·\left(1+\frac{{\text{C}}_{\text{t}}/{\text{V}}_{\text{d}}}{{\text{K}}_{\text{M}}}\right)+\frac{{\text{C}}_{\text{t}}\prime }{{\text{V}}_{\text{d}}\prime }}& (17)\end{array}$

With v_{max_m}′ being the effective maximum rate at which glucose analog can be phosphorylated, C_t′ the concentration of glucose analog in the tissue, V_d′ the physical distribution space of glucose analog in the brain, and ${{\text{K}}_{\text{M}}}^{\prime }$ such that ${\text{V}}_{\text{d}}·{{\text{K}}_{\text{M}}}^{\prime }$ is the concentration at which ${\text{M}}^{\prime }={{\text{v}}_{max\text{\_}\text{m}}}^{\prime }/2$ under condition no native glucose in the tissue (C_t = 0). In the following, we set ${{\text{v}}_{max\text{\_}\text{m}}}^{\prime }=0.3·{\text{v}}_{max\text{\_}\text{m}}$ (Kuwabara et al., 1990; Kuwabara and Gjedde, 1991). Equation (17) can be simplified to:

${\text{M}}^{\prime }\left({\text{C}}_{\text{t}},{\text{v}}_{\max \_\text{m}}\right)={\text{v}}_{\max \_\text{m}}\prime ·\frac{{\text{C}}_{\text{t}}\prime \left(\tau ,{\text{v}}_{\max \_\text{m}}\right)/{\text{V}}_{\text{d}}\prime }{{\text{K}}_{\text{M}}+{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)/{\text{V}}_{\text{d}}}$

in the limit where C_t is much greater than ${\text{C}}_{\text{t}}^{\prime }$ and setting ${{\text{K}}_{\text{M}}}^{\prime }$ = K_M, assuming identical hexokinase affinities for the two substrates. Note that, because of the concentrations involved, glucose analog transport is influenced by native glucose concentration, whereas the contrary is not true.

Extraction fraction from the plasma

We consider glucose analog in the same two compartments as native glucose, and we make the same assumptions to determine its concentration along the capillary axis:

$\begin{array}{l}\frac{\partial {\text{C}}_{\text{p}}\prime \left(\text{x},\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)}{\partial \text{x}}=-\frac{\tau }{{\text{V}}_{\text{cap}}}·{\text{v}}_{\max \_\text{t}}\prime ·\\ \frac{{\text{C}}_{\text{p}}\prime \left(\text{x},\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)-\frac{{\text{C}}_{\text{t}}\prime \left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}\prime }}{{\text{K}}_{\text{T}}\prime ·\left(1+\frac{{\text{C}}_{\text{p}}\left(\text{x},\tau ,{\text{C}}_{\text{t}}\right)}{{\text{K}}_{\text{T}}}+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}·{\text{K}}_{\text{T}}}\right)+{\text{C}}_{\text{p}}\prime \left(\text{x},\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)+\frac{{\text{C}}_{\text{t}}\prime \left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}\prime }}& (18)\end{array}$

Where ${\text{C}}_{\text{p}}^{\prime }$ denotes the glucose analog concentration in plasma, ${{\text{v}}_{max\text{\_}\text{t}}}^{\prime }$ the maximum rate at which it can be transported across the capillary membrane, and ${\text{K}}_{\text{T}}^{\prime }$ the Michaelis-Menten parameters for glucose tracer transport. The glucose analog concentration at the point x = 0 is equal to glucose tracer arterial concentration ${\text{C}}_{\text{A}}^{\prime }$, which is assumed to be negligible compared to C_A (tracer concentration).

Hence, Equation (18) can be simplified

$\begin{array}{l}\frac{\partial {\text{C}}_{\text{p}}\prime \left(\text{x},\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)}{\partial \text{x}}=\\ -\frac{\tau }{{\text{V}}_{\text{cap}}}·{\text{v}}_{\max \_\text{t}}\prime .\frac{{\text{C}}_{\text{p}}\prime \left(\text{x},\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)-\frac{{\text{C}}_{\text{t}}\prime \left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}\prime }}{{\text{K}}_{\text{T}}\prime ·\left(1+\frac{{\text{C}}_{\text{p}}\left(\tau ,{\text{C}}_{\text{t}}\right)}{{\text{K}}_{\text{T}}}+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}·{\text{K}}_{\text{T}}}\right)}& (19)\end{array}$

Where we set ${{\text{v}}_{max\text{\_}\text{t}}}^{\prime }=1.4·{\text{v}}_{max\text{\_}\text{t}}$ (Crane et al., 1983; Kuwabara et al., 1990; Hasselbalch et al., 1996a), and approximate in the following ${{\text{K}}_{\text{T}}}^{\prime }$ and ${{\text{V}}_{\text{d}}}^{\prime }$ to be equal to K_T and V_d, respectively.

Hence, for glucose analog, we must solve the following system of coupled equation:

$\begin{array}{ll}\{\begin{array}{l}\text{(a)}-{\text{v}}_{\max \_\text{m}}\prime ·\frac{{\text{C}}_{\text{t}}\prime \left(\tau ,{\text{v}}_{\max \_\text{m}}\right)/{\text{V}}_{\text{d}}}{{\text{K}}_{\text{M}}+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}+\frac{{\text{C}}_{\text{A}}\prime ·{\text{Q}}^{\prime }\left(\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)}{\tau }·\\ {\text{V}}_{\text{cap}}=0,\text{ with }{\text{Q}}^{\prime }\left(\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)=1-\frac{{\text{C}}_{\text{p}}\prime \left(1;\tau \right)}{{\text{C}}_{\text{p}}\prime \left(0;\tau \right)}\\ \text{(b)}\frac{\partial {\text{C}}_{\text{p}}\prime \left(\text{x},\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)}{\partial \text{x}}=-\frac{\tau }{{\text{V}}_{\text{cap}}}·{\text{v}}_{\max \_\text{t}}\prime .\frac{{\text{C}}_{\text{p}}\prime \left(\text{x},\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)-\frac{{\text{C}}_{\text{t}}\prime \left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}}{{\text{K}}_{\text{T}}+{\text{C}}_{\text{p}}\left(\tau ,{\text{C}}_{\text{t}}\right)+\frac{{\text{C}}_{\text{t}}\left(\tau ,{\text{v}}_{\max \_\text{m}}\right)}{{\text{V}}_{\text{d}}}},\\ \text{ with }{\text{C}}_{\text{p}}\prime \left(0,\tau ,{\text{C}}_{\text{t}},{\text{C}}_{\text{t}}\prime \right)={\text{C}}_{\text{A}}\prime \end{array}& (20)\end{array}$

for ${{\text{C}}_{\text{p}}}^{\prime }$ and ${{\text{C}}_{\text{t}}}^{\prime }$, for any value of the transit time τ.

In the following, we will not indicate the dependence of ${{\text{C}}_{\text{p}}}^{\prime }$ and ${{\text{C}}_{\text{t}}}^{\prime }$ on v_{max_m} and C_t, as these functions have been determined previously. We follow the same strategy to solve this system as we did for native glucose. Accordingly, we solve Equation (20)b independently over a grid of values $\left(\tau ,{{\text{C}}_{\text{t}}}^{\prime }\right)$. We compute the corresponding Q′ function on the same grid $\left(\tau ,{{\text{C}}_{\text{t}}}^{\prime }\right)$:

$\begin{array}{ll}{\text{Q}}^{\prime }\left(\tau ,{\text{C}}_{{\text{t}}^{\prime }}\right)=1-\frac{{\text{C}}_{{\text{p}}^{\prime }}\left(1;\tau ,{\text{C}}_{{\text{t}}^{\prime }}\right)}{{\text{C}}_{{\text{p}}^{\prime }}\left(0;\tau ,{\text{C}}_{{\text{t}}^{\prime }}\right)}& (21)\end{array}$

In the second step, we numerically solve the equation:

$\begin{array}{ll}-{\text{v}}_{\max \_\text{m}}\prime ·\frac{\frac{{\text{C}}_{\text{t}}\prime \left(\tau \right)}{{\text{V}}_{\text{d}}}}{{\text{K}}_{\text{M}}+\frac{{\text{C}}_{\text{t}}\prime \left(\tau \right)}{{\text{V}}_{\text{d}}}}+\frac{{\text{C}}_{\text{A}}\prime ·{\text{Q}}^{\prime }\left(\tau ,{\text{C}}_{\text{t}}\prime \right)}{\tau }·{\text{V}}_{\text{cap}}=0& (22)\end{array}$

for relevant values of τ, to obtain ${\text{C}}_{\text{t}}^{\prime }\left(\tau \right)$ in steady-state. Having determined the specific value of ${\text{C}}_{\text{t}}^{\prime }\left(\tau ,{\text{v}}_{max\text{\_}\text{m}}\right)$, we use the previously interpolated Q′ to determine ${\text{C}}^{\prime }\left(\tau ,{{\text{C}}_{\text{t}}}^{\prime }\right)$ for any transit time τ.

From Equation (15), we can determine the mean glucose tracer extraction fraction (GEF′). The LC is then simply derived as the ratio $\text{LC}=\frac{{\text{GEF}}^{\prime }}{\text{GEF}}$.

Lactate Production

We used our model to predict lactate production in different physiological states. To do so, we decompose CMR_glc in two parts. CMR_glc = M_a + M_s, where M_a corresponds to glucose fully oxidized to CO₂, and M_s corresponds to glucose non-oxidatively metabolized to lactate. As each molecule of glucose metabolized non-oxidatively produces two molecules of lactate, the rate at which lactate is produced, P_l, fulfills:

$\begin{array}{ll}{\text{P}}_{\text{l}}=2·{\text{M}}_{\text{s}}& (23)\end{array}$

Let us introduce OGI, the ratio between CMRO₂ and CMR_glc, with OGI < 6.

$\begin{array}{ll}\left({\text{M}}_{\text{a}}+{\text{M}}_{\text{s}}\right)·\text{OGI}={\text{CMRO}}_2& (24)\end{array}$

As M_a corresponds to glucose metabolized oxidatively, CMRO₂/M_a = 6.

We can rewrite Equation (24):

$\begin{array}{ll}\left(\frac{{\text{CMRO}}_2}{6}+{\text{M}}_{\text{s}}\right)·\text{OGI}={\text{CMRO}}_2& (25)\end{array}$

That is,

$\begin{array}{ll}{\text{M}}_{\text{s}}={\text{CMRO}}_2·\frac{6-\text{OGI}}{6·\text{OGI}}& (26)\end{array}$

Finally,

$\begin{array}{ll}{\text{P}}_{\text{l}}={\text{CMRO}}_2·\frac{6-\text{OGI}}{3·\text{OGI}}={\text{CMR}}_{\text{glc}}·\frac{6-\text{OGI}}{3}& (27)\end{array}$

with P_l being the rate at which lactate is produced, expressed in the same units as CMR_glc.

Results

In this section, we first present results regarding the dependency of baseline CMR_glc on plasma concentration. We then show how CMR_glc is predicted to vary between physiological conditions and compare these changes to those reported in the literature. Experimental evidence suggests that CBF and CMR_glc change in proportion (see Paulson et al., 2010, for references and an in-depth discussion). To examine whether our model produces predictions that are consistent with these observations, we determined CMR_glc and the phosphorylation rate of FDG (which we refer to as CMR_FDG) as a function of CBF for (MTT, CTH) states between baseline and activation. Figure 2 shows that an almost linear relation between CBF and CMR_glc, as well as between CBF and CMR_FDG is obtained, which is consistent with the experimental findings, although increases in CBF are not a prerequisite in the support of increase glucose metabolism (Leithner and Royl, 2014). We also compare changes in CMR_glc to that of CMRO₂, to determine the OGI. Finally, we assess the effects of a change in plasma concentration on the LC and quantify the extent to which it is predicted to change between physiological states.

FIGURE 2

Figure 2. Relation between CBF and CMR_glc (red curve) and between CBF and the phosphorylation rate of FDG, that is, CMR_FDG (blue curve), for physiological states between baseline (MTT = 1.4 s, CTH = 1.33 s) and stimulation (MTT = 0.81 s, CTH = 0.52 s), where CTH is assumed to vary linearly with MTT, using parameter set S_H: v_{max_t} = 68 μmol/100 mL/min and K_T = 50 μmol/100 mL_brain. CBF is related to MTT through the relation CBF = V_cap/MTT. CMR_glc in baseline condition is equal to 29μmol/100 mL/min. bl, baseline; CBF, cerebral blood flow; V_cap, cerebral blood volume.

CMR_glc in Baseline Condition

Figure 3 shows the relation between glucose concentration in plasma and CMR_glc in the baseline condition. The CMR_glc predicted by our model using parameter set S_R (see Methods: Calibration of the Model Parameters) was 72 μmol/100 mL/min, slightly lower than the value reported in conscious rats, see for example (Choi et al., 2002) for review. This underestimation could be caused by the fact that we derived our transport parameters from anesthetized animal data, see discussion in Section Calibration of the Model Parameters. We used a plasma glucose concentration of 1000 μmol/100 mL_plasma in our calculations, as reported for rats under euglycemic conditions. The CMR_glc predicted by our model using parameter set S_H was 28 μmol/100 mL/min, compared to 30μmol/100 mL/min assumed in human when we derived transport parameters. In this case, we used a plasma glucose concentration of 500μmol/100 mL_plasma, as reported for humans under euglycemic conditions.

FIGURE 3

Figure 3. Relation between glucose concentration in plasma and CMR_glc in baseline condition. The purple and green solid lines correspond to the mean glucose concentration in the tissue (input) as a function of arterial glucose concentration according to relations (13) and (14), based on experimental studies by Dienel et al. (1991) and Gruetter et al. (1998), respectively. The blue dotted-dashed line corresponds to predictions obtained using parameter set S_R: v_{max_t} = 136 μmol/100 mL/min and K_T = 150 μmol/100 mL_brain. The red dotted-dashed line corresponds to parameter set S_H: v_{max_t} = 68 μmol/100 mL/min and K_T = 50 μmol/100 mL_brain. The gray area corresponds to conditions of normoglycemia in humans (plasma glucose concentration ranging from 400 to 600 μmol/100 mL).

For both set of Michaelis-Menten parameters, our model shows that CMR_glc in baseline conditions does not depend on arterial glucose concentration until it reaches values as low as 200–300 μmol/100 mL_plasma, concentration at which the mean glucose concentration in the tissue approaches zero. In particular, CMR_glc under condition of mild hypoglycemia (C_A = 300–400μmol/100 mL_plasma) and hyperglycemia is found to be almost equal to CMR_glc under condition of normoglycemia (relative difference lower than 10%).

Relative Changes in CMR_glc between Physiological States

Figure 4A shows a contour plot of CMR_glc, using parameter set S_H, that is, when the model is applied to human. Figures 4B,C show CMRO₂ contour plots obtained with the model from Angleys et al. (2015) with a maximum metabolic rate for oxygen metabolism (v_max) equal to 4.75 mL/100 mL/min and 5.23 mL/100 mL/min (that is, 10% higher), respectively. Note that there is no straightforward way to illustrate CMR_glc during the transition from the baseline condition (bl) characterized by MTT_bl, CTH_bl, and v_max, to an activated state (act) characterized by MTT_act, CTH_act, and v_max+10% in a single contour plot without specific knowledge on the relation between MTT, CTH, and v_max. In Figure 5, v_max is assumed to increase in proportion to the stimulus intensity, by 10% between baseline and stimulation, which corresponds to the value that yielded the most realistic results in Angleys et al. (2015).

FIGURE 4

Figure 4. (A) CMR_glc contour plot, assuming a glucose arterial plasma concentration equal to 500μmol/100 mL, a mean glucose concentration in the tissue under baseline condition equal to 129μmol/100 mL, and Michaelis-Menten parameters for transport of glucose v_{max_t} and K_T equal to 68μmol/100 mL/min and 50μmol/100 mL, respectively. (B) CMRO₂ contour plot (Angleys et al., 2015) assuming oxygen metabolism to be governed by Michaelis-Menten kinetics, with parameters K_M = 2.71mmHg (3.5 μmol/L) and v_max = 4.75 mL/100 mL/min. (C) CMRO₂ contour plot assuming v_max = 5.23 mL/100 mL/min [i.e., 10% higher than in (B)] and making otherwise the same assumptions as in (B). The parameter v_max in (B,C) is assumed to be constant. The yellow line separates states where a blood flow increase (decrease in MTT) given a fixed CTH will lead to an increased (right side of the line) or decreased (left side of the line) glucose (A) or oxygen (B,C) consumption, respectively. The roman numeral accompanying each symbol corresponds to physiological data. The numeral “0” stands for resting state, whereas other numerals refer to state of altered basal physiology. Note that the CMR_glc and CMRO₂ iso-contours do not show the same slope for the experimental data used in this figure, indicating that a change in MTT (resp. CTH) will have a strong (resp. moderate) influence on CMR_glc, and inversely for CMRO₂. Symbols: +, functional activation (Stefanovic et al., 2008); ×, cortical electrical stimulation (Schulte et al., 2003).

FIGURE 5

Figure 5. OGI (A) and lactate production (B) in different physiological conditions indicated by symbols on Figure 4. We assume a glucose arterial plasma concentration equal to 500 μmol/100 mL, a mean glucose concentration in the tissue under baseline condition equal to 129 μmol/100 mL, and Michaelis-Menten parameters for transport of glucose v_{max_t} and K_T equal to 68μmol/100 mL/min and 150μmol/100 mL (Cunningham et al., 1986), respectively. We furthermore assume that oxygen metabolism is governed by Michaelis-Menten kinetics, with parameters K_M = 2.71 mmHg (3.5μmol/L) and v_max = 4.75 mL/100 mL/min, and that v_max increases by 10% from baseline condition to state I (+) or state V (×). In (B) lactate production is derived from CMR_glc from the formula P_I = CMR_glc · (6−OGI)/3, with P_I being lactate production, and OGI being computed in (A). OGI: oxygen-glucose index: CMRO₂/CMR_glc.

One important observation in Figure 4A is that variations in CMR_glc as function of MTT are expected to be higher than variations in CMRO₂. The OGI ratio therefore depends on the physiological state considered and is expected to decrease during activation, as illustrated in Figure 5A. When we consider changes in MTT and CTH derived from Stefanovic et al. (2008) and Schulte et al. (2003) (Figure 5A), our model predicts that CMR_glc increases by 31% (43%) from baseline condition to stimulated state, while CMRO₂ is expected to increase by only 19% (18%) between baseline and stimulation, resulting in an OGI decreasing by 10–20%, from 5.5 (5.8) to 5.0 (4.7). As a result, lactate production is expected to increase from 3 to 5 μmol/100 mL/min to about 15 to 20μmol/100 mL/min (Figure 5B).

Additional information that can be inferred from Figure 4 is that, for the considered physiological states, CMR_glc is influenced primarily by changes in MTT, while a change in CTH has little influence on this variable. This is in contrast to CMRO₂, which is expected to be influenced primarily by CTH and more moderately by MTT.

Models from Jespersen and Østergaard (2012) and Angleys et al. (2015) predict that for large CTH values and under a condition of fixed CTH, a blood flow increase leads to a decrease in oxygen delivery. This phenomenon which has been referred to as malignant CTH is observed in this model with glucose as well. Accordingly, for states on the left hand side of the yellow line in Figure 4A, CMR_glc decreases if flow increases under condition of constant CTH.

LC Changes in Response to a Change in Plasma Concentration and between Physiological Conditions

Figure 6 shows the relation between plasma glucose concentration and the expected value for the LC in baseline condition and during stimulation, assuming the relations in Equations (13) and (14) between ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ and C_A. Accordingly, the LC was found to be equal to 0.61 and 0.76 with parameter sets S_H and S_R, respectively under baseline, normoglycemic conditions (C_A = 500μmol/100 mL and 1000μmol/100 mL, for set S_H and S_R, respectively). Our model predicts that LC increases under conditions of severe hypoglycemia (C_A = 210 μmol/100 mL) to 0.68 (1.11), and decreases under conditions of hyperglycemia to 0.58 (0.69) with parameter set S_H (S_R).

FIGURE 6

Figure 6. Relation between glucose concentration in plasma and the lumped constant (LC) for FDG in baseline condition (dashed-dotted line) and during stimulation (dotted line). The black and the red lines correspond to predictions obtained using parameter set S_R: v_{max_t} = 136μmol/100mL/min and K_T = 150μmol/100 mL_brain. The yellow and the blue lines correspond to predictions obtained using parameter set S_H: v_{max_t} = 68 μmol/100 mL/min and K_T = 50 μmol/100 mL_brain. The purple and green solid lines correspond to the mean glucose concentration in the tissue (input) as a function of arterial glucose concentration according to relations (13) and (14), based on experimental studies by Dienel et al. (1991) and Gruetter et al. (1998), respectively. The gray area corresponds to conditions of normoglycemia in humans (plasma glucose concentration ranging from 400 to 600μmol/100 mL). bl, baseline condition; st, stimulation.

For arterial glucose concentrations ranging from 400 to 3000 μmol/100 mL, and assuming that changes in MTT and CTH between baseline condition and stimulation are accompanied by a 30% decrease in tissue glucose concentration, our model predicts that LC increases by about 15 and 20% when using Michaelis-Menten parameter sets S_R and S_H, respectively. Figure 7 shows the extent to which these changes in LC would impact the computation of relative increase in CMR_glc in human when neglected, provided that the uptake of glucose analog from the blood is measured accurately. In the following, CMR_{glc, app} denotes the apparent metabolic rate of glucose, when neglecting changes in LC value between physiological conditions. Relative change in CMR_glc, CMR_{glc, app}, and LC are related according to the equation:

$\begin{array}{ll}\frac{\Delta {\text{CMR}}_{\text{glc},\text{app}}\left(\text{d}\right)}{{\text{CMR}}_{\text{glc},\text{0},\text{app}}}=\left(1+\frac{\Delta \text{LC}\left(\text{d}\right)}{{\text{LC}}_0}\right)·\frac{\Delta {\text{CMR}}_{\text{glc}}\left(\text{d}\right)}{{\text{CMR}}_{\text{glc},\text{0}}}+\frac{\Delta \text{LC}\left(\text{d}\right)}{{\text{LC}}_0}& (28)\end{array}$

It should be noted that Equation (28) predicts CMR_{glc, app} and CMR_glc to be equal only when LC remains constant during a change in physiological state (ΔLC = 0). In Equation (28), LC₀ and CMR_{glc, 0} is the LC and CMR_glc value, respectively, in the baseline condition, and ΔLC and ΔCMR_glc the subsequent increase in LC and CMR_glc, respectively. Here, d denotes the relative decrease in tissue glucose concentration between conditions. Note that ΔCMR_glc, ΔLC, and therefore ΔCMR_{glc, app} depend on d. The blue curve in Figure 7 shows the relation between CMR_glc and CMR_{glc, app}, using parameter set S_H, with d varying from 0 to 40%. This range allows ΔCMR_{glc, app} to vary in the range from 7 to 83%, therefore covering the relative increase in glucose metabolism reported in the literature (see Table 3).

FIGURE 7

Figure 7. Relation between the relative (true) CMR_glc increase, and the apparent relative CMR_{glc, app} increase computed when assuming that the lumped constant does not vary between physiologic states (blue line). The relationship between CMR_{glc, app}, CMR_glc, and LC is given in Equation (28). Note that the increase in glucose metabolism is overestimated by more than 50% for any “true” relative CMR_glc increases of 12% or more if one neglects the state-dependency of LC. In this plot, parameter set S_H is used to describe glucose transport, and we assume that the uptake of glucose tracer from the blood is measured accurately. The green line shows the resulting relative overestimation. The black line is the line of equation y = x.

TABLE 3

Table 3. Comparison between relative increases in glucose metabolism reported in experiments involving NMRS and PET in human.

The function that associates CMR_{glc, app} to CMR_glc shows a derivative which is increasing with CMR_glc. The relative overestimation (green curve in Figure 7) therefore increases as a function of CMR_glc and takes particularly high values for high CMR_glc relative increases. Because of the non-linearity between relative increases in CMR_glc and LC, Equation (28) leads to surprisingly large overestimations of the increase in CMR_glc if even small changes in LC are neglected. For example, considering conservative values for LC and CMR_glc increases of 10% (see Dienel et al., 1991, for experimental evidence) and 20%, respectively, would lead to an apparent 32% increase in CMR_glc (Equation 28), thus overestimating the increase in CMR_glc by 60%. Also note that Equation (28) is model-independent and therefore does not rely on simplifying assumptions used in our model. Please refer to the Supplementary Material for details about its derivation.

Discussion

In this study, we developed a model that takes the effects of tissue glucose concentration and CTH into account when describing glucose extraction in the brain. We employed reversible Michaelis-Menten kinetics, which has previously been shown to support cerebral glucose utilization across a range of arterial glucose/tissue concentration (Choi et al., 2001; van de Ven et al., 2012).

The first main finding in our study is that, similar to oxygen extraction, glucose extraction is not only a function of the CBF and concentration in plasma and tissue, but also depends on capillary transit time heterogeneity. However, as glucose and oxygen transport involve different mechanisms, the changes in CBF and CTH between physiological states do not affect the transport of these two substrates to the same extent. In particular, Figure 4A shows that for physiological CTH/MTT ratio (typically smaller than 2), CMR_glc iso-contours are almost vertical, meaning that CTH does not influence glucose uptake to a very large extent. In contrast, for the same CTH/MTT ratio, CMRO₂ iso-contour slopes are lower (Figures 4B,C), although they would appear to be slightly higher if v_max was not kept constant in the figures (see Section Relative Changes in CMRglc between Physiological States for details about this point). Blood flow increases are therefore expected to be less efficient as a means of increasing oxygen than glucose consumption for a given CTH value. As a consequence, glucose delivery can increase more than that of oxygen during enhanced energy demand, favoring non-oxidative glucose consumption. To better illustrate that CTH has a larger influence on oxygen than glucose delivery, we made a simulation in which we assumed that no change in CTH occurred between physiological states, keeping other parameters unchanged. Instead of increasing by 19% (18%) as reported when CTH decreases in parallel with MTT, our model predicts that CMRO₂ would increase by only 3.6% (7.4%) when using MTT and CTH derived from Stefanovic et al. (2008) and Schulte et al. (2003). In contrast, CMR_glc would be less affected, increasing by 23% (38%), compared to 31% (43%) reported earlier in this manuscript. As a consequence, decreases in OGI (16% (22%)) would be larger than predicted when CTH varies (9% (18%)). In other words, during activation, blood flow homogenization is expected to increase oxygen extraction capacity to a greater extent than that of glucose, thus limiting OGI reduction.

The second main finding is that the ratio of glucose tracer to native glucose extraction at steady state (LC) depends on the physiological state we consider, and varies accordingly between baseline and stimulated conditions. This finding is in contrast to previous studies, where this ratio is considered to be constant. We show that neglecting variations in this ratio could lead to overestimations of the relative increase in CMR_glc (ΔCMR_glc/CMR_{glc, 0}) between different physiological states of 50% or more. In the following, we discuss how such variations in the LC could reconcile what previously seemed to be incompatible measurements obtained with PET and NMRS, respectively.

Changes in Plasma Concentration

Data from the literature reports that glucose metabolism is largely insensitive to changes in plasma concentration. It was therefore important to show that our model is consistent with this data. Accordingly, we tested our model over a range of plasma concentrations to simulate conditions of hypoglycemia and hyperglycemia, keeping the model parameters which describe glucose transport unchanged. We employed two sets of parameters to describe glucose transport for rats and humans, respectively. Our assumption that glucose transport is described by reversible Michaelis-Menten kinetics requires also to assume that tissue glucose concentration varies linearly with plasma concentration in the range from hypoglycemia to hyperglycemia levels (see Equations (13) and (14)), in which CMR_glc has been reported to remain almost constant. This linear relation is supported by several experimental studies (Gruetter et al., 1998; Choi et al., 2001; Seaquist et al., 2001; van de Ven et al., 2012). When glucose plasma concentration increases above normoglycemic levels, our model predicts that CMR_glc is essentially unchanged. The same observation applies under conditions of hypoglycemia, where the model predicts that CMR_glc only decreases slightly until plasma concentrations reaches concentration as low as 200–300 μmol/100 mL_plasma; levels at which glucose concentration in tissue approaches zero. When plasma concentration decreases further, CMR_glc is expected to decrease sharply.

Our results are in good agreement with several studies showing that CMR_glc only decreases slightly under conditions of moderate hypoglycemia in rats (Bryan et al., 1986). Moreover, several studies report a sharp increase in CBF when glucose plasma concentration decreases to concentration lower than 200μmol/100 mL (Choi et al., 2001), which corresponds to concentration in the tissue close to zero. It suggests that glucose metabolism is more severely impaired when plasma concentration is lower than 200μmol/100 mL_plasma than during mild hypoglycemia, and that CMR_glc cannot be maintained at levels observed during normoglycemia. Finally, although CMR_glc predicted by our model when applied to rats is slightly lower than reported values (see Section CMR_glc in Baseline Condition for explanation of this underestimation), predictions for human were in good agreement with the value we assumed in our calibration.

In summary, our model captures crucial characteristics of glucose delivery to the brain, in particular the remarkable ability of the brain to maintain sufficient glucose metabolism across a wide range of plasma and tissue glucose concentrations. This is in contrast to oxygen metabolism, which is much more sensitive to changes in plasma concentration and CBF. As a result, the brain is much more vulnerable to change in oxygen than in glucose levels (Leithner and Royl, 2014).

Changes in CMR_glc between Physiological States

We employed our model to predict CMR_glc in different physiological states. These predictions are based on the assumption that the mean glucose concentration in the tissue (i) decreases by 30% between the two physiological states that we took as reference for baseline and stimulation (See sections Calibration of the Models Parameters in Methods and in Discussion for further discussion about this choice), and (ii) is linearly related to MTT. While the first assumption is based on experimental data, a precise relationship between tissue glucose concentration and MTT has yet to be determined, and future work should therefore test the validity of the second assumption. With recent technical developments, tissue glucose concentrations may indeed be measured (Lugo-Morales et al., 2013), and this aspect of our models could therefore be tested experimentally.

Our model predicts that CMR_glc increases relatively more than CMRO₂ between baseline and stimulation. Accordingly, the OGI is predicted to decrease by about 10–20% and lactate production to increase, which is in good agreement with literature reports (Prichard et al., 1991; Madsen et al., 1995, 1998; Frahm et al., 1996). The reasonable number of parameters employed in our model may provide insights into the understanding of the regulation of glucose and oxygen metabolism. Indeed, the strong inhibition of hexokinase in baseline condition allows glucose metabolism to increase more than that of oxygen metabolism, which has been reported to be approximately 80–85% of its maximum value, already at rest (Gjedde et al., 2005). It is also worth noting that the transport capacity of glucose across the BBB varies with the concentration gradient in a non-linear fashion because of glucose transporters properties. In particular, if an increase in concentration gradient is caused by a decrease in tissue glucose levels, as observed during functional activation, glucose transport capacity will increase relatively more than does the gradient. In contrast if the increase in concentration gradient is caused by an increase in arterial glucose concentrations, as observed during normal physiological variations or during hyperglycemia, glucose transport capacity will increase much less than does the gradient. For example, with reversible Michaelis-Menten kinetics (e.g., parameter set S_H) and the parameters used in our study, a 30% reduction in tissue glucose levels from 132 to 92 μmol/100 mL_brain will increase the concentration gradient by 11% and the net transport rate by 25%. If the same increase in concentration gradient was introduced through an increase in plasma concentration, then the transport rate would increase by only 6.2%. The transport capacity of glucose can thus amplify or attenuate the effect of changes in concentration gradient, and thereby accommodate both high metabolic demands and changes in plasma concentration without the risk of a tissue energy crisis. These uptake properties are reflected by our model response to a change in concentration gradient introduced through a change in plasma concentration and a decrease in tissue concentration, respectively. While a relatively small decrease in tissue concentration is indeed accompanied by a large increase in the net transport rate and hence in CMR_glc (Figure 4), a large change in plasma concentration is compensated by a comparatively reduced change in tissue concentration (Figure 3), leading to a stable transport rate and CMR_glc. These properties are in contrast to free diffusion, which is characterized by a transport rate across the BBB that is proportional to the concentration gradient. These observations further support that non-oxidative glucose consumption may therefore be an inherent consequence of glucose and oxygen extraction, especially during functional activation. Needless to say, the way in which the brains cell types have specialized to utilize the available glucose and lactate as substrates in their metabolism is beyond the scope of this paper.

Lumped Constant

With our model, we can compute the ratio of glucose analog to native glucose extraction (LC). We tested our model over a large range of glucose concentrations, to quantify the extent to which this ratio is predicted to change as a function of the plasma concentration, and to assess whether this is in good agreement with literature reports.

We observe that LC increases from 0.76 with parameter set S_R to more than one under severely hypoglycemic conditions, and decreases to 0.68 under hyperglycemic conditions. Variations in LC with parameter set S_H are predicted to show the same pattern but with smaller amplitude. These LC values under baseline and normoglycemic conditions are in good agreement with some experiments obtained for FDG in rats (Tokugawa et al., 2007) and in humans (Hasselbalch et al., 2001a), but about 50% higher than others obtained in rats (Huang et al., 1980). See Hasselbalch et al. (2001a) for in-depth discussions of this variability. The variation of the LC under hypo- and hyperglycemic conditions for rats is in line with experimental reports (Suda et al., 1990; Dienel et al., 1991) and modeling (Pardridge et al., 1982; Holden et al., 1991). As pointed out by Pardridge et al. (1982), Crane et al. (1983) and Holden et al. (1991), these variations in LC can be explained by differences between the rate at which native glucose and glucose analog are metabolized and transported across the BBB. Indeed, under hyperglycemia/normoglycemia, glucose metabolism is not limited by transport across the capillary membrane. Under severe hypoglycemia, however, transport capacity limits glucose delivery. As a result, the LC ranges between two extreme values, namely (i) the ratio between the maximal operation rate of hexokinase for glucose analog and native glucose, ${{\text{v}}_{max\text{\_}\text{m}}}^{\prime }/{\text{v}}_{max\text{\_}\text{m}}$, equal to 0.3 in our study, which is approached under hyperglycemic conditions, and (ii) the ratio between the maximal rate at which glucose analog and native glucose are transported across the BBB, that is ${{\text{v}}_{max\text{\_}\text{t}}}^{\prime }/{\text{v}}_{max\text{\_}\text{t}}$, equal to 1.4, which is approached under hypoglycemic conditions.

In summary, both CMR_glc and LC values yielded by our model, as well as their relative changes between two glycemic conditions are consistent with the literature.

We also quantified the change in the value of LC between two physiological states (baseline condition and activation). Our model predicts that, under conditions of normoglycemia and constant blood glucose concentration, the ratio of glucose tracer to native glucose extraction (so-called LC) increases by 15–20% from baseline condition to stimulation, which is in line with earlier predictions by Dienel et al. (1991) based on typical decrease in tissue glucose concentration during activation. As explained in the Supplementary Material and established experimentally (Dienel et al., 1991), we show that this variation is partly due to changes in glucose concentration in the tissue, which does vary substantially from baseline condition to stimulated state (e.g., see Merboldt et al., 1992; Chen et al., 1993; Adachi et al., 1995; Frahm et al., 1996; Mangia et al., 2006; Lin et al., 2012). Accordingly, we believe it is crucial to adjust LC when quantifying the relative increase in CMR_glc between two physiological states. Changes in LC between physiological states would have to be determined experimentally (e.g., Dienel et al., 1991), or from model predictions.

Comparison of Relative Increases in CMR_glc Obtained with NMRS and PET

In Table 3, we list several experiments where NMRS and FDG PET were used to quantify relative increases in CMR_glc between baseline and an activated state. All experiments involved visual stimulation and all but one (Phelps et al., 1981) used the same kind of stimulus, that is, a checkboard or red diodes flashing at 8 or 10 Hz (see Table 3). The table reveals tendency for PET measurements to yield higher CMR_glc increases than those obtained with NMRS (mean relative increase = 43% with PET, 22% with NMRS). When taking into account the change in LC between baseline and stimulation by using the relation between apparent glucose metabolism increase and true glucose increase illustrated in Figure 7, the mean relative increase in CMR_glc is reduced to 24% when assessed with PET. The relative standard deviation (standard deviation to mean ratio) is furthermore decreased by 40% when applying the correction, which reflects the better agreement between measures obtained with NMRS and PET after correction than before. Experiments involving NMRS use bigger regions of interest (ROIs) and therefore may contain more white matter than experiments with PET, which might lead to a slight underestimation of the signal during activation. However, in these experiments, ROIs are of comparable sizes and locations, so that signal underestimation with NMRS, if any, is likely to be small. Although more data would be needed to calibrate our model more precisely and hence to achieve a more accurate correction, we believe that such corrections are crucial when inferring glucose uptake from FDG uptake in different physiological states, without which CMR_glc estimate may be misleading.

Applying Our Model to Disease Conditions

We applied our model to tumor cells, as an example of possible clinical application. More specifically, we assessed the consequences on the glucose uptake and on the LC of an altered glucose transport or metabolism, such as it has been reported in the literature on tumors: Tumor cells express an isozyme of hexokinase, hexokinase II, which is less susceptible to feedback inhibition by its product, glucose-6-phosphate, than hexokinase in healthy tissue (Bustamante and Pedersen, 1977; Bustamante et al., 1981). When applying our model to tumor cells, to take this reduced inhibition into account, we have increased the maximum rate at which glucose and FDG can be metabolized (v_{max_m} and ${{\text{v}}_{max\text{\_}\text{m}}}^{\prime }$). We find that the LC increases by 36 and 66%, from 0.61 in baseline condition to 0.83 and 1.01, respectively, under condition of reduced inhibition, while the mean glucose concentration in the tissue is reduced to one half and one fourth of healthy tissue values, respectively. This increase in the LC is accompanied by a 45 and 75% increase in CMR_glc compared to healthy tissue. These results are in line with the literature, reporting that most tumor tissues are known to be highly metabolic and depend on aerobic glycolysis (the Warburg effect Warburg, 1956). Spence and colleagues furthermore showed that the LC in tumors was higher than in healthy tissue and was generally found to exceed unity (Spence et al., 1990, 1998).

Our model therefore suggests that smaller inhibition of hexokinase could be a possible factor leading to both increased glucose metabolism and LC value. As discussed in previous sections, underestimation of the LC leads to overestimation of CMR_glc as assessed by FDG PET. When neglecting LC increases, quantitative CMR_glc estimates in tumors therefore contains little extra information compared to qualitative measures, and could even turn out to be misleading. Overestimating CMR_glc in tumors would in turn lead to OGI underestimation and hence to an overestimation of non-oxidative glucose consumption.

Our model also allows us to understand the biophysical mechanisms that lead to a higher LC value. Because hexokinase is less inhibited in tumor than in healthy tissues, glucose phosphorylation is quicker and the equilibrium concentrations in tissues are lower. Consequently, once a glucose molecule crosses the blood-tumor-barrier, its life time in the tissue compartment is reduced and glucose metabolism is limited to a greater extent by glucose transport across the BBB than in healthy tissue. This suggests that glucose transport capacity in tumor cells is reduced relatively to the phosphorylation rate. The LC value therefore becomes weighted to a greater extent by the ratio between native glucose and glucose analog transport rates across the BBB (${{\text{v}}_{max\text{\_}\text{t}}}^{\prime }/{\text{v}}_{max\text{\_}\text{t}}$) than by their phosphorylation rates (${{\text{v}}_{max\text{\_}\text{m}}}^{\prime }/{\text{v}}_{max\text{\_}\text{m}}$), and hence increases essentially for the same reasons as under conditions of hypoglycemia.

In future work, we could employ our model to better understand disease states related to aerobic glycolysis. For example, in Alzheimer disease (AD), several studies show that areas of the normal human brain with elevated aerobic glycolysis are nearly identical to those that accumulate amyloid and exhibit atrophy and disrupted metabolism in AD (Buckner et al., 2005; Vlassenko et al., 2010). It has therefore been suggested that there might be a link between dependence on aerobic glycolysis and AD. It would be interesting to employ our model to assess the extent to which microvascular dysfunction such as it has been reported in AD and parameters such as OGI are related.

Calibration of the Model Parameters

Our model involves several parameters to describe glucose transport. We calibrated most of the parameters with values from the literature. However, the parameters derived in different studies are not easy to compare, because they involve several kinetic models, such as irreversible and reversible Michaelis-Menten kinetics. Choi et al. (2001) made a comparison between parameter values derived when assuming reversible and irreversible Michaelis-Menten, respectively, in human and in rats, while Cunningham et al. (1986) made such as comparison in rats only. Although the parameters of the reversible Michaelis-Menten model take lower values than the non-reversible version, there is no general relation between them. Some studies have included a non-saturable component, often called K_d, making the comparison between parameters even more difficult. Based on four studies, we applied our model to rats and to humans, using reversible Michaelis-Menten kinetics, in particular to assess how our results depend on a particular choice of parameters. Anesthetics affect the cerebral metabolic rate of glucose (see for example Choi and Gruetter, 2012) and therefore the relationship between glucose concentration in plasma and tissue. While calibrating S_R, we had to rely on previous studies using reversible Michaelis-Menten kinetics in anesthetized rats (Cunningham et al., 1986; Choi et al., 2001). Choi et al. (2001), however, determined v_{max_t}/CMR_glc ratio rather than v_{max_t} alone. Provided that anesthesia impacts metabolism rather than glucose transport across the BBB itself, the parameter v_{max_t} is in principle unaffected by the effects of anesthesia. This is in contrast to the second study (Cunningham et al., 1986), where v_{max_t} is determined independently of CMR_glc. As discussed by Cunningham and colleagues', the value determined for v_{max_t} would have been higher if the experiment had been performed on unanesthetized animals. Our study may therefore have underestimated v_{max_t} in parameter set S_R, and thus the CMR_glc predicted for rats. Note that parameter set S_H was based on data from unanesthetized volunteers and the corresponding CMR_glc predictions thus apply to awake humans.

Reported values for the Michaelis-Menten constant K_T range from −98 to 330 μmol/100 mL in the four studies we used to calibrate our parameters. Moreover, one study (Cunningham et al., 1986) involves rats anesthetized with pentobarbital, which is known to inhibit glucose transport by binding to the glucose transporter itself (el-Barbary et al., 1996; Haspel et al., 1999). The value determined for this parameter (194 μmol/100 mL) is therefore likely to be lower in conscious rats. Accordingly, we chose a slightly lower value for K_T, 150 μmol/100 mL.

While some of our conclusions are found to be largely insensitive to the choice of parameters, others depend more on the values assigned to describe glucose transport. For example, employing different sets of parameters leads to the same conclusion when one computes CMR_glc as a function of arterial plasma concentration, namely that CMR_glc is essentially independent to the plasma concentration over a large range of concentration (300–3000 μmol/100 mL) and decreases only when the plasma concentration is so low that glucose concentration in the tissue approaches zero (Figure 3). However, the predicted CMR_glc in baseline condition is more dependent on Michaelis-Menten parameters. This is reflected by the results shown in Figure 3, where CMR_glc for rats is predicted to be more than two-fold higher than for humans.

In this study, we use our model to assess the extent to which CMR_glc and OGI vary between physiological states, and we assume that ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ decreases by 30% from baseline condition to stimulation. The extent to which tissue concentration decreases during, e.g., photic stimulation has been widely debated and reported values range between 0 and 50% (Collins et al., 1987; Merboldt et al., 1992; Chen et al., 1993; Adachi et al., 1995; Frahm et al., 1996; Mangia et al., 2006; Lin et al., 2012). These differences may be attributed to differences in stimuli involved, experimental protocols and measurement techniques. We estimated our model's sensitivity to this assumption by assuming different ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ reduction, ranging between 10 and 40%. CMR_glc is predicted to increases by 41% (59%), 31% (43%), 22% (29%), and 14% (16%) when ${\overline{\text{C}}}_{\text{t}}^{\text{t}}$ is assumed to decrease by 40, 30, 20, and 10%, and using changes in MTT and CTH derived from Stefanovic et al. (2008) and Schulte et al. (2003). While this large range of predicted CMR_glc increases reflects the variability observed in literature reports, it would be valuable to have additional simultaneous recordings of local CMRO₂, CMR_glc, CBF, and CTH to sharpen our predictions. We also assume in our model that glucose concentration in plasma and in the tissue are related according to Equations (13) and (14), when the model is applied to rats and humans, respectively. For a given plasma concentration, a large range of tissue concentration have been reported, see for example (Dienel et al., 1991; Madsen et al., 1999) for measures in rats, and therefore other relations could equally have been used. Although our overall conclusions do not depend on these particular choices, more accurate estimates would improve our predictions.

Glucose metabolism is commonly assessed by two different glucose analogs: FDG and 2-DG. FDG is mainly used in humans with PET, while 2-DG is used in rodents with autoradiographic methods. In this study, we focused on the LC in human and chose parameters to describe FDG transport and metabolism rather than 2-DG. While the use of literature values for 2-DG would result in different rate constants, we expect that 2-DG and FDG transport are sufficiently similar for the overall conclusions of our study to hold for both tracers. As discussed above (please see Section Lumped constant in the discussion), the LC can be seen as a weighted average between transport capacity and phosphorylation rate ratios of glucose analogs to native glucose, which in our case are equal to 0.3 and 1.4, respectively. The LC therefore varies between these two values, depending on the physiological condition, which affects the weight given to the transport and to the phosphorylation rates. Consequently, increasing the ratio ${{\text{v}}_{max\text{\_}\text{t}}}^{\prime }/{\text{v}}_{max\text{\_}\text{t}}$ to 1.5 as it is sometimes reported for FDG would make LC vary by 22% between baseline condition and stimulation, and by 17% if this ratio is equal to 1.1, compared to 21% in our study, using otherwise the same parameters as in Figure 6, with C_A = 500μmol/100 mL and parameter set S_H.

Reversible Michaelis-Menten Kinetics

Although modeling of glucose transport across the BBB has been dominated by non-reversible Michaelis-Menten kinetics, we chose to employ reversible Michaelis-Menten kinetics in our study.

Glucose transport across the BBB involves transporters (essentially GLUT-1), which have been studied extensively, and several models have been derived to describe their kinetics (Carruthers and Helgerson, 1991; Baldwin, 1993; Mueckler, 1994; Cloherty et al., 1996). These models are characterized by a high degree of complexity, and all of them allow glucose to bind back to the transporter just after its release, which is supported experimentally (Carruthers and Helgerson, 1991; Cloherty et al., 1996). Although reversible Michaelis-Menten is an oversimplification of these models, it includes this latter possibility, in contrast to non-reversible Michaelis-Menten kinetics, and seems therefore more suitable to describe glucose transport.

Although other studies employ reversible Michaelis-Menten kinetics, only two to our knowledge assessed the relevance of employing reversible instead of non-reversible Michaelis-Menten kinetics (Cunningham et al., 1986; Gruetter et al., 1998). In these studies, it has been shown that employing reversible Michaelis-Menten generally leads to better agreement with experimental data (Cunningham et al., 1986; Gruetter et al., 1998; Seaquist et al., 2001). For example, under constant CMR_glc, non-reversible Michaelis-Menten-kinetics predicts a non-linear relation between plasma and tissue concentration, the latter being limited by a saturation value. In contrast, reversible Michaelis-Menten predicts a linear relation between tissue and plasma concentration (Gruetter et al., 1998). Experimentally, it has been shown in several studies that concentrations in the tissue and plasma are linearly related to one another over a large range of glucose plasma concentration (250–3000μmol/100 mL) (Gruetter et al., 1998; Choi et al., 2001; Seaquist et al., 2001; van de Ven et al., 2012), suggesting that reversible Michaelis-Menten kinetics is more suitable than non-reversible Michaelis-Menten to describe glucose transport. Finally, reversible Michaelis-Menten kinetics does not involve the use of non-saturable component of unidirectional influx as is the case in several studies employing non-reversible Michaelis-Menten (Cremer and Cunningham, 1979; Pardridge et al., 1982), and therefore better reflects the underlying process involved in facilitated diffusion.

Neglecting Phosphatase Activity

In our model, we assume that glucose tracer is trapped in the tissue once it has been phosphorylated, and we therefore neglect the back reaction which is catalyzed by glucose-6-phosphatase, allowing glucose tracer to be dephosphorylated. In studies employing first order rate constants, the constant for this latter reaction is often designated as ${{\text{k}}_4}^{*}$ and its value is reported to be equal to approximatively one tenth of ${{\text{k}}_3}^{*}$, the first order rate constant for the phosphorylation of glucose by hexokinase. Several studies addressed the effects of phosphatase activity on the estimation of LC and CMR_glc (Nelson et al., 1986; Dienel et al., 1988; Schmidt et al., 1992; Gotoh et al., 2000; Hasselbalch et al., 2001b) and suggested that they are negligible for experimental periods shorter than 45 min, and begin to appear with increasing time. Including glucose-phosphatase activity in our model would require the use of additional parameters and to make more assumptions, including assumptions regarding the setup of the experiment (tracer concentration, time of the infusion). While we think that it would be valuable to get quantitatively more precise results, the overall conclusions of our study are not believed to be sensitive to this particular choice.

Neglecting Glucose Transfer between Extravascular Compartments

In our model, we neglect any effects of glucose transfer between extravascular compartments. In reality, glucose might diffuse from nearby microvessels to counteract any effects of CTH, meaning that our model might overestimate tissue concentration-, and therefore concentration gradient heterogeneity. On the basis of CMR_glc, ${\overline{\text{C}}}_{\text{t}}$, and the diffusion coefficient for glucose in water, the diffusion time of glucose before glycolysis is estimated to be 5 min corresponding to a diffusion distance of 400 μm, which is approximately 10 times that of oxygen, and significantly exceeds the intercapillary distance. This suggests that glucose diffuses not only to the tissue immediately surrounding the nearest capillary, as assumed in our model, but also to the tissue further away. To estimate how this assumption impacts our conclusions, we tested an alternative version of the model, where we assumed no tissue concentration heterogeneity by considering only one extravascular compartment where glucose is well stirred. The results obtained with these two different assumptions lead to almost same results (not shown), suggesting that glucose mixing between extravascular compartments would not affect our conclusion. We also quantified the extent to which CTH introduces heterogeneity at the tissue concentration level, by considering the distribution of extravascular glucose concentrations in the state (MTT_bl, CTH_bl) and with the parameter set S_H. We observed that the highest concentration in the decile (resp. centile) showing the lowest concentration differs by only 22% (resp. 59%) from the lowest concentration in the decile (resp. centile) showing the highest concentration. As a comparison, the transit time of the blood flowing through these compartments differs by 7 (resp. 43)-fold. We speculate therefore that tissue glucose concentration heterogeneity induced by CTH is too small to impact glucose delivery.

Author Contributions

HA developed the biophysical model used in this study, contributed to interpretation of the results and wrote the first draft of the manuscript. SJ developed the study concept, contributed to the development of the biophysical model used in this study and to the interpretation of the results. LØ developed the study concept and contributed to the interpretation of the results.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This study was supported by the Danish National Research Foundation (CFIN), the Danish Ministry of Science, Innovation, and Education (MINDLab), and the VELUX Foundation (ARCADIA). The authors wish to thank Olaf Paulson for helpful discussions and suggestions during the preparation of the manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fncom.2016.00103

The model source code can be found online at: https://github.com/HugoAngleys/glucose_extraction

Abbreviations

ANLS, astrocyte-neuron lactate shuttle; BBB, blood brain barrier; CBF, cerebral blood flow; CMR_glc, cerebral metabolic rate of glucose; CMRO₂, cerebral metabolic rate of oxygen; CTH, capillary transit time heterogeneity; GEF, glucose extraction fraction; LC, lumped constant; MTT, mean transit time; NMRS, nuclear magnetic resonance spectroscopy; OEF, oxygen extraction fraction; OGI, oxygen-glucose index; PET, positron emission tomography.

References

Adachi, K., Cruz, N. F., Sokoloff, L., and Dienel, G. A. (1995). Labeling of metabolic pools by [6-14C] glucose during K+-induced stimulation of glucose utilization in rat brain. J. Cereb. Blood Flow Metab. 15, 97–110. doi: 10.1038/jcbfm.1995.11

CrossRef Full Text | Google Scholar

Angleys, H., Østergaard, L., and Jespersen, S. N. (2015). The effects of capillary transit time heterogeneity (CTH) on brain oxygenation. J. Cereb. Blood Flow Metab. 35, 806–817. doi: 10.1038/jcbfm.2014.254

PubMed Abstract | CrossRef Full Text | Google Scholar

Baldwin, S. A. (1993). Mammalian passive glucose transporters: members of an ubiquitous family of active and passive transport proteins. Biochim. Biophys. Acta BBA Rev. Biomembr. 1154, 17–49. doi: 10.1016/0304-4157(93)90015-G

PubMed Abstract | CrossRef Full Text

Berg, J. M., Tymoczko, J. L., and Stryer, L. (2012). Biochemistry: [This Edition is for Use Outside the USA and Canada], 7th Edn. International [Nachdr.]. New York, NY: Freeman, Palgrave Macmillan.

Bryan, R. M., Keefer, K. A., and MacNeill, C. (1986). Regional cerebral glucose utilization during insulin-induced hypoglycemia in unanesthetized rats. J. Neurochem. 46, 1904–1911. doi: 10.1111/j.1471-4159.1986.tb08512.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckner, R. L., Snyder, A. Z., Shannon, B. J., LaRossa, G., Sachs, R., Fotenos, A. F., et al. (2005). Molecular, structural, and functional characterization of Alzheimer's disease: evidence for a relationship between default activity, amyloid, and memory. J. Neurosci. 25, 7709–7717. doi: 10.1523/JNEUROSCI.2177-05.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Bustamante, E., Morris, H. P., and Pedersen, P. L. (1981). Energy metabolism of tumor cells. Requirement for a form of hexokinase with a propensity for mitochondrial binding. J. Biol. Chem. 256, 8699–8704.

PubMed Abstract | Google Scholar

Bustamante, E., and Pedersen, P. L. (1977). High aerobic glycolysis of rat hepatoma cells in culture: role of mitochondrial hexokinase. Proc. Natl. Acad. Sci. U.S.A. 74, 3735–3739. doi: 10.1073/pnas.74.9.3735

PubMed Abstract | CrossRef Full Text | Google Scholar

Carruthers, A., and Helgerson, A. L. (1991). Inhibitions of sugar transport produced by ligands binding at opposite sides of the membrane. Evidence for simultaneous occupation of the carrier by maltose and cytochalasin B. Biochemistry (Mosc.) 30, 3907–3915. doi: 10.1021/bi00230a015

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Novotny, E. J., Zhu, X.-H., Rothman, D. L., and Shulman, R. G. (1993). Localized 1H NMR measurement of glucose consumption in the human brain during visual stimulation. Proc. Natl. Acad. Sci. U.S.A. 90, 9896–9900. doi: 10.1073/pnas.90.21.9896

PubMed Abstract | CrossRef Full Text | Google Scholar

Chih, C.-P., and Roberts, E. L. (2003). Energy substrates for neurons during neural activity: a critical review of the astrocyte-neuron lactate shuttle hypothesis. J. Cereb. Blood Flow Metab. 23, 1263–1281. doi: 10.1097/01.WCB.0000081369.51727.6F

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, I.-Y., and Gruetter, R. (2012). Neural Metabolism In vivo, Advances in Neurobiology. New York, NY: Springer Science & Business Media.

Choi, I.-Y., Lee, S.-P., Kim, S.-G., and Gruetter, R. (2001). In vivo measurements of brain glucose transport using the reversible michaelis–menten model and simultaneous measurements of cerebral blood flow changes during hypoglycemia. J. Cereb. Blood Flow Metab. 21, 653–663. doi: 10.1097/00004647-200106000-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, I.-Y., Lei, H., and Gruetter, R. (2002). Effect of deep pentobarbital anesthesia on neurotransmitter metabolism in vivo: on the correlation of total glucose consumption with glutamatergic action. J. Cereb. Blood Flow Metab. 22, 1343–1351. doi: 10.1097/01.WCB.0000040945.89393.46

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, D. D., Lajtha, A. L., and Maker, H. S. (1989). “Intermediary metabolism,” in Basic Neurochemistry (New York, NY: Raven Press), 542–550.

Cloherty, E. K., Heard, K. S., and Carruthers, A. (1996). Human erythrocyte sugar transport is incompatible with available carrier models. Biochemistry (Mosc.) 35, 10411–10421. doi: 10.1021/bi953077m

PubMed Abstract | CrossRef Full Text | Google Scholar

Collins, R. C., McCandless, D. W., and Wagman, I. L. (1987). Cerebral glucose utilization: comparison of [14C] deoxyglucose and [6-14C] glucose quantitative autoradiography. J. Neurochem. 49, 1564–1570. doi: 10.1111/j.1471-4159.1987.tb01028.x

CrossRef Full Text | Google Scholar

Crane, P. D., Pardridge, W. M., Braun, L. D., and Oldendorf, W. H. (1983). Kinetics of transport and phosphorylation of 2-fluoro-2-deoxy-d-glucose in rat brain. J. Neurochem. 40, 160–167. doi: 10.1111/j.1471-4159.1983.tb12666.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cremer, J. E., and Cunningham, V. J. (1979). Effects of Some Chlorinated Sugar Derivatives on the Hexose Transport System of the Blood/Brain Barrier. Available online at: http://www.biochemj.org/bj/180/bj1800677.htm [Accessed May 22, 2015].

Cunningham, V. J. (1986). The influence of transport and metabolism on brain glucose content. Ann. N.Y. Acad. Sci. 481, 161–173. doi: 10.1111/j.1749-6632.1986.tb27148.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cunningham, V. J., Hargreaves, R. J., Pelling, D., and Moorhouse, S. R. (1986). Regional blood-brain glucose transfer in the rat: a novel double-membrane kinetic analysis. J. Cereb. Blood Flow Metab. Off. J. Int. Soc. Cereb. Blood Flow Metab. 6, 305–314. doi: 10.1038/jcbfm.1986.53

PubMed Abstract | CrossRef Full Text | Google Scholar

Dienel, G. A., Cruz, N. F., Mori, K., Holden, J. E., and Sokoloff, L. (1991). Direct measurement of the λ of the lumped constant of the deoxyglucose method in rat brain: determination of λ and lumped constant from tissue glucose concentration or equilibrium brain/plasma distribution ratio for methylglucose. J. Cereb. Blood Flow Metab. 11, 25–34. doi: 10.1038/jcbfm.1991.3

PubMed Abstract | CrossRef Full Text | Google Scholar

Dienel, G. A., Nelson, T., Cruz, N. F., Jay, T., Crane, A. M., and Sokoloff, L. (1988). Over-estimation of glucose-6-phosphatase activity in brain in vivo. Apparent difference in rates of [2-3H] glucose and [U-14C] glucose utilization is due to contamination of precursor pool with 14C-labeled products and incomplete recovery of 14C-labeled metabolites. J. Biol. Chem. 263, 19697–19708.

Google Scholar

Edvinsson, L., and Krause, D. N. (eds.). (2002). Cerebral Blood Flow and Metabolism, 2nd Edn. Philadelphia: Lippincott Williams & Wilkins.

el-Barbary, A., Fenstermacher, J. D., and Haspel, H. C. (1996). Barbiturate inhibition of GLUT-1 mediated hexose transport in human erythrocytes exhibits substrate dependence for equilibrium exchange but not unidirectional sugar flux. Biochemistry (Mosc.) 35, 15222–15227. doi: 10.1021/bi962050f

PubMed Abstract | CrossRef Full Text | Google Scholar

Fox, P. T., Raichle, M. E., Mintun, M. A., and Dence, C. (1988). Nonoxidative glucose consumption during focal physiologic neural activity. Sci. Wash. 241, 462–462. doi: 10.1126/science.3260686

PubMed Abstract | CrossRef Full Text

Frahm, J., Krüger, G., Merboldt, K. D., and Kleinschmidt, A. (1996). Dynamic uncoupling and recoupling of perfusion and oxidative metabolism during focal brain activation in man. Magn. Reson. Med. Off. J. Soc. Magn. Reson. Med. Soc. Magn. Reson. Med. 35, 143–148. doi: 10.1002/mrm.1910350202

PubMed Abstract | CrossRef Full Text | Google Scholar

Gjedde, A., and Christensen, O. (1984). Estimates of michaelis-menten constants for the two membranes of the brain endothelium. J. Cereb. Blood Flow Metab. 4, 241–249. doi: 10.1038/jcbfm.1984.33

PubMed Abstract | CrossRef Full Text | Google Scholar

Gjedde, A., Johannsen, P., Cold, G. E., and Østergaard, L. (2005). Cerebral metabolic response to low blood flow: possible role of cytochrome oxidase inhibition. J. Cereb. Blood Flow Metab. 25, 1183–1196. doi: 10.1038/sj.jcbfm.9600113

PubMed Abstract | CrossRef Full Text | Google Scholar

Gotoh, J., Itoh, Y., Kuang, T.-Y., Cook, M., Law, M. J., and Sokoloff, L. (2000). Negligible Glucose-6-Phosphatase Activity in Cultured Astroglia. J. Neurochem. 74, 1400–1408. doi: 10.1046/j.1471-4159.2000.0741400.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gruetter, R., Ugurbil, K., and Seaquist, E. R. (1998). Steady-state cerebral glucose concentrations and transport in the human brain. J. Neurochem. 70, 397–408. doi: 10.1046/j.1471-4159.1998.70010397.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Haspel, H. C., Stephenson, K. N., Davies-Hill, T., El-Barbary, A., Lobo, J. F., Croxen, R. L., et al. (1999). Effects of barbiturates on facilitative glucose transporters are pharmacologically specific and isoform selective. J. Membr. Biol. 169, 45–53. doi: 10.1007/PL00005900

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasselbalch, S. G., Holm, S., Pedersen, H. S., Svarer, C., Knudsen, G. M., Madsen, P. L., et al. (2001a). The 18F-fluorodeoxyglucose lumped constant determined in human brain from extraction fractions of 18f-fluorodeoxyglucose and glucose. J. Cereb. Blood Flow Metab. 21, 995–1002. doi: 10.1097/00004647-200108000-00012

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasselbalch, S. G., Knudsen, G. M., Capaldo, B., Postiglione, A., and Paulson, O. B. (2001b). Blood-brain barrier transport and brain metabolism of glucose during acute hyperglycemia in humans. J. Clin. Endocrinol. Metab. 86, 1986–1990. doi: 10.1210/jc.86.5.1986

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasselbalch, S. G., Knudsen, G. M., Holm, S., Pinborg Hageman, L., Capaldo, B., and Paulson, O. B. (1996a). Transport of D-glucose and 2-fluorodeoxyglucose across the blood-brain barrier in humans. J. Cereb. Blood Flow Metab. 16, 659–666. doi: 10.1097/00004647-199607000-00017

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasselbalch, S. G., Madsen, P. L., Hageman, L. P., Olsen, K. S., Justesen, N., Holm, S., et al. (1996b). Changes in cerebral blood flow and carbohydrate metabolism during acute hyperketonemia. Am. J. Physiol. Endocrinol. Metab. 270, E746–E751.

PubMed Abstract | Google Scholar

Hertz, L. (2004). The astrocyte-neuron lactate shuttle: a challenge of a challenge. J. Cereb. Blood Flow Metab. 24, 1241–1248. doi: 10.1097/00004647-200411000-00008

PubMed Abstract | CrossRef Full Text | Google Scholar

Holden, J. E., Mori, K., Dienel, G. A., Cruz, N. F., Nelson, T., and Sokoloff, L. (1991). Modeling the dependence of hexose distribution volumes in brain on plasma glucose concentration: implications for estimation of the local 2-deoxyglucose lumped constant. J. Cereb. Blood Flow Metab. 11, 171–182. doi: 10.1038/jcbfm.1991.50

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, S.-C., Phelps, M. E., Hoffman, E. J., Sideris, K., Selin, C. J., and Kuhl, D. E. (1980). Noninvasive determination of local cerebral metabolic rate of glucose in man. Am. J. Physiol. Endocrinol. Metab. 238, E69–E82.

PubMed Abstract | Google Scholar

Jespersen, S. N., and Østergaard, L. (2012). The roles of cerebral blood flow, capillary transit time heterogeneity, and oxygen tension in brain oxygenation and metabolism. J. Cereb. Blood Flow Metab. Off. J. Int. Soc. Cereb. Blood Flow Metab. 32, 264–277. doi: 10.1038/jcbfm.2011.153

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuwabara, H., Evans, A. C., and Gjedde, A. (1990). Michaelis-menten constraints improved cerebral glucose metabolism and regional lumped constant measurementswith [18F]fluorodeoxyglucose. J. Cereb. Blood Flow Metab. 10, 180–189. doi: 10.1038/jcbfm.1990.33

CrossRef Full Text | Google Scholar

Kuwabara, H., and Gjedde, A. (1991). Measurements of glucose phosphorylation with FDG and PET are not reduced by dephosphorylation of FDG-6-phosphate. J. Nucl. Med. 32, 692–698.

Google Scholar

Leithner, C., and Royl, G. (2014). The oxygen paradox of neurovascular coupling. J. Cereb. Blood Flow Metab. 34, 19–29. doi: 10.1038/jcbfm.2013.181

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Y., Stephenson, M. C., Xin, L., Napolitano, A., and Morris, P. G. (2012). Investigating the metabolic changes due to visual stimulation using functional proton magnetic resonance spectroscopy at 7 T. J. Cereb. Blood Flow Metab. 32, 1484–1495. doi: 10.1038/jcbfm.2012.33

PubMed Abstract | CrossRef Full Text | Google Scholar

Lugo-Morales, L. Z., Loziuk, P. L., Corder, A. K., Toups, J. V., Roberts, J. G., McCaffrey, K. A., et al. (2013). Enzyme-modified carbon-fiber microelectrode for the quantification of dynamic fluctuations of nonelectroactive analytes using fast-scan cyclic voltammetry. Anal. Chem. 85, 8780–8786. doi: 10.1021/ac4017852

PubMed Abstract | CrossRef Full Text | Google Scholar

Lund-Andersen, H. (1979). Transport of glucose from blood to brain. Physiol. Rev. 59, 305–352.

PubMed Abstract | Google Scholar

Madsen, P. L., Cruz, N. F., Sokoloff, L., and Dienel, G. A. (1999). Cerebral Oxygen/glucose ratio is low during sensory stimulation and rises above normal during recovery: excess glucose consumption during stimulation is not accounted for by lactate efflux from or accumulation in brain tissue. J. Cereb. Blood Flow Metab. 19, 393–400. doi: 10.1097/00004647-199904000-00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Madsen, P. L., Hasselbalch, S. G., Hagemann, L. P., Olsen, K. S., Bülow, J., Holm, S., et al. (1995). Persistent resetting of the cerebral oxygen/glucose uptake ratio by brain activation: evidence obtained with the Kety-Schmidt technique. J. Cereb. Blood Flow Metab. Off. J. Int. Soc. Cereb. Blood Flow Metab. 15, 485–491. doi: 10.1038/jcbfm.1995.60

PubMed Abstract | CrossRef Full Text | Google Scholar

Madsen, P. L., Linde, R., Hasselbalch, S. G., Paulson, O. B., and Lassen, N. A. (1998). Activation-induced resetting of cerebral oxygen and glucose uptake in the rat. J. Cereb. Blood Flow Metab. 18, 742–748. doi: 10.1097/00004647-199807000-00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Mangia, S., Tkáč, I., Gruetter, R., Van de Moortele, P.-F., Maraviglia, B., and Uğurbil, K. (2006). Sustained neuronal activation raises oxidative metabolism to a new steady-state level: evidence from 1H NMR spectroscopy in the human visual cortex. J. Cereb. Blood Flow Metab. 27, 1055–1063. doi: 10.1038/sj.jcbfm.9600401

PubMed Abstract | CrossRef Full Text | Google Scholar

McIlwain, H., and Bachelard, H. S. (1985). Biochemistry and the Central Nervous System, 5th Edn. Edinburgh; New York, NY: Churchill Livingstone.

Merboldt, K.-D., Bruhn, H., Hanicke, W., Michaelis, T., and Frahm, J. (1992). Decrease of glucose in the human visual cortex during photic stimulation. Magn. Reson. Med. 25, 187–194. doi: 10.1002/mrm.1910250119

PubMed Abstract | CrossRef Full Text | Google Scholar

Michaelis, L., and Menten, M. L. (1913). Die kinetik der invertinwirkung. Biochem. Z 49, 352.

Google Scholar

Mueckler, M. (1994). Facilitative glucose transporters. Eur. J. Biochem. 219, 713–725. doi: 10.1111/j.1432-1033.1994.tb18550.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelson, T., Lucignani, G., Goochee, J., Crane, A. M., and Sokoloff, L. (1986). Invalidity of criticisms of the deoxyglucose method based on alleged glucose-6-phosphatase activity in brain. J. Neurochem. 46, 905–919. doi: 10.1111/j.1471-4159.1986.tb13057.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Orzi, F., Lucignani, G., Dow-Edwards, D., Namba, H., Nehlig, A., Patlak, C. S., et al. (1988). Local cerebral glucose utilization in controlled graded levels of hyperglycemia in the conscious rat. J. Cereb. Blood Flow Metab. 8, 346–356. doi: 10.1038/jcbfm.1988.70

PubMed Abstract | CrossRef Full Text | Google Scholar

Pardridge, W. M., Crane, P. D., Mietus, L. J., and Oldendorf, W. H. (1982). Kinetics of regional blood–brain barrier transport and brain phosphorylation of glucose and 2-deoxyglucose in the barbiturate–anesthetized rat. J. Neurochem. 38, 560–568. doi: 10.1111/j.1471-4159.1982.tb08663.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Paulson, O. B., Hasselbalch, S. G., Rostrup, E., Knudsen, G. M., and Pelligrino, D. (2010). Cerebral blood flow response to functional activation. J. Cereb. Blood Flow Metab. 30, 2–14. doi: 10.1038/jcbfm.2009.188

PubMed Abstract | CrossRef Full Text | Google Scholar

Pellerin, L., and Magistretti, P. J. (2003). Food for thought: challenging the dogmas. J. Cereb. Blood Flow Metab. 23, 1282–1286. doi: 10.1097/01.WCB.0000096064.12129.3D

PubMed Abstract | CrossRef Full Text | Google Scholar

Phelps, M. E., Kuhl, D. E., and Mazziota, J. C. (1981). Metabolic mapping of the brain's response to visual stimulation: studies in humans. Science 211, 1445–1448. doi: 10.1126/science.6970412

PubMed Abstract | CrossRef Full Text | Google Scholar

Prichard, J., Rothman, D., Novotny, E., Petroff, O., Kuwabara, T., Avison, M., et al. (1991). Lactate rise detected by 1H NMR in human visual cortex during physiologic stimulation. Proc. Natl. Acad. Sci. U.S.A. 88, 5829–5831. doi: 10.1073/pnas.88.13.5829

PubMed Abstract | CrossRef Full Text | Google Scholar

Rasmussen, P. M., Jespersen, S. N., and Østergaard, L. (2015). The effects of transit time heterogeneity on brain oxygenation during rest and functional activation. J. Cereb. Blood Flow Metab. 35, 432–442. doi: 10.1038/jcbfm.2014.213

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, K., Lucignani, G., Moresco, R. M., Rizzo, G., Gilardi, M. C., Messa, C., et al. (1992). Errors introduced by tissue heterogeneity in estimation of local cerebral glucose utilization with current kinetic models of the [18F]fluorodeoxyglucose method. J. Cereb. Blood Flow Metab. 12, 823–834. doi: 10.1038/jcbfm.1992.114

PubMed Abstract | CrossRef Full Text | Google Scholar

Schulte, M. L., Wood, J. D., and Hudetz, A. G. (2003). Cortical electrical stimulation alters erythrocyte perfusion pattern in the cerebral capillary network of the rat. Brain Res. 963, 81–92. doi: 10.1016/S0006-8993(02)03848-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Seaquist, E. R., Damberg, G. S., Tkac, I., and Gruetter, R. (2001). The effect of insulin on in vivo cerebral glucose concentrations and rates of glucose transport/metabolism in humans. Diabetes 50, 2203–2209. doi: 10.2337/diabetes.50.10.2203

PubMed Abstract | CrossRef Full Text | Google Scholar

Spence, A. M., Graham, M. M., Muzi, M., Abbott, G. L., Krohn, K. A., Kapoor, R., et al. (1990). Deoxyglucose lumped constant estimated in a transplanted rat astrocytic glioma by the hexoseutilization index. J. Cereb. Blood Flow Metab. 10, 190–198. doi: 10.1038/jcbfm.1990.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Spence, A. M., Muzi, M., Graham, M. M., O'Sullivan, F., Krohn, K. A., Link, J. M., et al. (1998). Glucose metabolism in human malignant gliomas measured quantitatively with PET, 1-[C-11]glucose and FDG: analysis of the FDG lumped constant. J. Nucl. Med. Off. Publ. Soc. Nucl. Med. 39, 440–448.

PubMed Abstract | Google Scholar

Stefanovic, B., Hutchinson, E., Yakovleva, V., Schram, V., Russell, J. T., Belluscio, L., et al. (2008). Functional reactivity of cerebral capillaries. J. Cereb. Blood Flow Metab. Off. J. Int. Soc. Cereb. Blood Flow Metab. 28, 961–972. doi: 10.1038/sj.jcbfm.9600590

PubMed Abstract | CrossRef Full Text | Google Scholar

Suda, S., Shinohara, M., Miyaoka, M., Lucignani, G., Kennedy, C., and Sokoloff, L. (1990). The lumped constant of the deoxyglucose method in hypoglycemia: effects of moderate hypoglycemia on local cerebral glucose utilization in the rat. J. Cereb. Blood Flow Metab. 10, 499–509. doi: 10.1038/jcbfm.1990.92

PubMed Abstract | CrossRef Full Text | Google Scholar

Tokugawa, J., Ravasi, L., Nakayama, T., Schmidt, K. C., and Sokoloff, L. (2007). Operational lumped constant for FDG in normal adult male rats. J. Nucl. Med. 48, 94–99.

PubMed Abstract | Google Scholar

van de Ven, K. C. C., van der Graaf, M., Tack, C. J., Heerschap, A., and de Galan, B. E. (2012). Steady-state brain glucose concentrations during hypoglycemia in healthy humans and patients with type 1 diabetes. Diabetes 61, 1974–1977. doi: 10.2337/db11-1778

PubMed Abstract | CrossRef Full Text | Google Scholar

Villien, M., Wey, H.-Y., Mandeville, J. B., Catana, C., Polimeni, J. R., Sander, C. Y., et al. (2014). Dynamic functional imaging of brain glucose utilization using fPET-FDG. Neuroimage 100, 192–199. doi: 10.1016/j.neuroimage.2014.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Vlassenko, A. G., Rundle, M. M., and Mintun, M. A. (2006). Human brain glucose metabolism may evolve during activation: findings from a modified FDG PET paradigm. Neuroimage 33, 1036–1041. doi: 10.1016/j.neuroimage.2006.06.065

PubMed Abstract | CrossRef Full Text | Google Scholar

Vlassenko, A. G., Vaishnavi, S. N., Couture, L., Sacco, D., Shannon, B. J., Mach, R. H., et al. (2010). Spatial correlation between brain aerobic glycolysis and amyloid- (A) deposition. Proc. Natl. Acad. Sci. U.S.A. 107, 17763–17767. doi: 10.1073/pnas.1010461107

PubMed Abstract | CrossRef Full Text | Google Scholar

Warburg, O. (1956). On the origin of cancer cells. Science 123, 309–314. doi: 10.1126/science.123.3191.309

PubMed Abstract | CrossRef Full Text | Google Scholar