The energetics of cerebral activity critically relies on the functional and metabolic interactions between neurons and astrocytes. Important open questions include the relation between neuronal versus astrocytic energy demand, glucose uptake and intercellular lactate transfer, as well as their dependence on the level of activity.
We have developed a large-scale, constraint-based network model of the metabolic partnership between astrocytes and glutamatergic neurons that allows for a quantitative appraisal of the extent to which stoichiometry alone drives the energetics of the system. We find that the velocity of the glutamate-glutamine cycle (Vcyc) explains part of the uncoupling between glucose and oxygen utilization at increasing Vcyc levels. Thus, we are able to characterize different activation states in terms of the tissue oxygen-glucose index (OGI). Calculations show that glucose is taken up and metabolized according to cellular energy requirements, and that partitioning of the sugar between different cell types is not significantly affected by Vcyc. Furthermore, both the direction and magnitude of the lactate shuttle between neurons and astrocytes turn out to depend on the relative cell glucose uptake while being roughly independent of Vcyc.
These findings suggest that, in absence of ad hoc activity-related constraints on neuronal and astrocytic metabolism, the glutamate-glutamine cycle does not control the relative energy demand of neurons and astrocytes, and hence their glucose uptake and lactate exchange.
Sustained cerebral activity is crucially dependent on the functional and metabolic interplay of neurons and glial cells. Spectroscopic and imaging methods have indeed shown that the brain accommodates a wealth of cell-to-cell interactions, which ultimately have contributed to displace the decades-old notion that merely coupled whole brain activity to neuronal glucose oxidation (for a comprehensive review, see ). In particular, carbohydrate metabolism is compartmentalized among neurons and astrocytes, which, together with the interstitial space, represent nearly 90% of the tissue. Although there is evidence for the trafficking of metabolic intermediates between the two cell types, its significance and dependence on the activation state are not fully elucidated. More than 15 years ago it was hypothesized that astrocytes may support the energetics of brain function by the provision of glucose-derived lactate to neurons, in an activity-dependent manner . However, the idea of a metabolically significant astrocyte-to-neuron lactate shuttle (ANLS), as well as the activity-dependent increase in astrocytic glucose uptake, has proven to be rather difficult to confirm in vivo, while indirect and not always reproducible experimental proof was mainly obtained from experiments on cell cultures (see [3, 4] and the excellent reviews [5, 6]).
The difficult interpretation and integration of the experimental findings produced a substantial theoretical effort aimed at characterizing intra- and inter-cellular metabolic fluxes [7–12]. So far, mathematical models of transport and metabolism of glucose in neurons and astrocytes using either kinetic [7, 9–11] or stoichiometry-based [13–15] approaches have provided conflicting results about the relevance of the cell-to-cell lactate shuttle (CCLS) (see  for a recent review). This, in turn, raised some debate, especially concerning the partitioning of glucose between neurons and astrocytes and the potentially resulting intercellular lactate flow . A recent flux-balance-analysis (FBA) study indicates that the direction and magnitude of the CCLS between neurons and astrocytes depends critically on the relative uptake of blood-borne glucose . The sharing of glucose between the two cell types is itself governed by the internal energetic demand of cells, implying that glucose partitioning alone cannot be used to draw any conclusion on the functional variations of the CCLS . A critical reassessment of previous modeling results suggests that the CCLS might remain of minor significance in terms of transferred carbon equivalents .
On the other hand, the known regulations of enzyme-catalyzed reactions implemented in dynamical models have so far proved insufficient to justify a fundamental energetic role for the CCLS . In particular, the differences between metabolic pathways of neuronal and astrocytic networks do not imply the occurrence of lactate exchange between cells, most likely because neurons and astrocytes do possess a relatively high self-sufficiency for both glycolytic and oxidative glucose metabolism (see  and references therein). This means that lactate is oxidatively metabolized in the same compartment where it is produced by glycolytic processing of glucose.
The aim of this study is to examine the activity-dependent metabolic cooperation of glutamatergic neurons and astrocytes from a network-based perspective. Specifically, two issues lie at the core of our work: (i) the correlation between partitioning of glucose and lactate shuttling; and (ii) their functional modulation across varying levels of glutamate-glutamine cycle. We have employed a constraint-based setting where an extensive and controlled sampling of the solution space is possible [21, 22] on a large-scale model of compartmentalized brain energy metabolism. At odds with previous studies employing constraint-based schemes to analyze the neuron-glia system in specific conditions (different from those considered here, see ), our approach does not rely on an objective function (which in our case would be hard to design) to define the relevant states. In addition, it allows to analyze in detail the feasible metabolic states for networks whose sizes are beyond those covered by other approaches like Bayesian Flux Balance Analysis (BFBA)  or Montecarlo sampling of mass-balance equations . Finally, we have not made any special assumption on the regulation of biochemical pathways with respect to the activation level, nor have we imposed specific constraints on transport fluxes (except for the uptake of glucose that we use to fine-tune the oxygen-to-glucose index (OGI), see below). In short, we show that, within our stoichiometric approach:
the OGI is able to distinguish states characterized by different levels of neurotransmission, as flux configurations with larger OGIs typically carry smaller values for the velocity of the glutamate-glutamine cycle;
the partitioning of glucose between neurons and astrocytes is roughly independent of the level of activity;
the magnitude and direction of the CCLS depend strongly on glucose partitioning while being roughly independent of the level of neurotransmission.
In other terms, within a purely stoichiometric model, the system’s energetics is determined to a significant extent by the sharing of glucose. These results support the idea that neurotransmission does not impose significant constraints on glucose partitioning or CCLS. In addition, we show that (d) the overall degree of correlation among metabolic reaction fluxes between and within cells changes drastically in the presence of neurotransmission, pointing to an extended metabolic and possibly functional partnership between neurons and astrocytes.
The reaction network we considered (a reduced sketch of which is given in Figure 1) is composed of four main compartments: neuron (n), astrocyte (a), extracellular space (e) and blood capillary (c). Within the neuronal and astrocytic elements we also distinguished the cytosol (nc and ac, respectively), mitochondria (nm and am) and synaptic vesicles (nv, only in neurons). Transport of nutrients from the blood to the brain parenchyma is provided by the capillary. We assumed that under resting conditions glucose and oxygen irreversibly enter the brain, while lactate is not significantly exchanged [25, 26]. Specifically, glucose can be taken up directly by astrocytes via the basal lamina or can diffuse into the extracellular space . The latter, in turn, is a common compartment for glucose uptake by neurons and astrocytes, as well as for lactate shuttling between the two cell types. Oxygen can freely diffuse from the capillary to cells. We lumped together the endothelium and basal lamina with the capillary compartment, which also means that we assume a negligible metabolism for endothelial cells.
The neuronal and astrocytic compartments are equipped with the enzymatic machinery to carry out the main pathways of carbohydrate metabolism (glycolysis, pentose phosphate shunt, TCA cycle, oxidative phosphorylation) . Both cell types indirectly transport reducing equivalents (i.e. NADH) from cytosol to mitochondria via malate-aspartate shuttle (MAS) (see ). We made the simplifying assumption that only astrocytes are capable of glutathione synthesis because neurons, unlike astrocytes, are unable to efficiently transport cystine  and, importantly, they cannot increase the substrate flow through glutamate-cystine ligase , the rate limiting step in glutathione synthesis. Yet, the antioxidant system is equally present in neurons and astrocytes to detoxify the reactive oxygen species (ROS) produced by oxidative phosphorylation. The stoichiometry of ROS production by oxidative phosphorylation was chosen assuming that 5% to 15% of glucose is processed through the pentose phosphate pathway to regenerate the NADPH required for reducing oxidized glutathione [29, 30]. Anaplerosis of TCA cycle intermediates is performed by pyruvate carboxylation, which is confined to astrocytes , as well as by the activity of the neuronal and astrocytic malic enzyme .
The functional portion of the metabolic network includes glutamatergic neurotransmission, transmitter recycling and ionic movements, that together establish the coupling between activity and metabolism through the action of the Na/K-ATPase . Specifically, the glutamate stored in neuronal synaptic vesicles can be released in the extracellular space, from where it is taken up by astrocytes in co-transport with three Na+ ions and counter-transport of one K+ ion. Glutamate is amidated to glutamine by astrocytic glutamine synthetase (GS) with the concurrent hydrolysis of one molecule of ATP. Glutamine is then exported to neurons where it is eventually converted back to glutamate and loaded into synaptic vesicles again, which costs another ATP. Astrocytic uptake of glutamate and release of glutamine, together with neuronal uptake of glutamine and release of glutamate configure the so-called glutamate-glutamine cycle. In this way, the clearance of neuronally released glutamate from the extracellular space is mostly accomplished by astrocytes , although a fraction of the neurotransmitter can be taken up by neurons, especially in synapses not associated with astrocytic processes . At odds with previous mass-balance modeling works [15, 17, 18], we included the ionic currents related to membrane depolarization, albeit these were not explicitly linked to glutamate release. In particular, neurons possess Na+ and K+ channels that mimic voltage-gated ion channels and astrocytes can also take up potassium from the extracellular space with the Na-K cotransporter. Overall, the fluxes of Na+ and K+ activate Na/K-ATPase, which consumes one ATP to transport three Na+ out of the cell and two K+ inside the cell. Importantly, not all the glutamate which is taken up by astrocytes is channeled via the glutamate-glutamine cycle. Glutamate in astrocytes can be used for energy production by entering the TCA cycle after conversion to alpha-ketoglutarate through transamination by aspartate aminotransferase (AAT) or dehydrogenation by glutamate dehydrogenase (GDH) . We did not include the action of other transaminase, e.g. alanine aminotransferase. This choice precludes testing the exchange of lactate and of alanine between neurons and astrocytes for maintaining ammonia homeostasis during glutamate-glutamine cycle . However, the role of this shuttle was experimentally found to be activity-independent in neuronal-astrocytic cultures . Finally, we conformed to other mass-balance modeling works [15, 17] in excluding the pathways involved in the synthesis and degradation of nucleic and amino acids. This is justified by the different characteristic time-scales of processes underlying energy metabolism and gene expression, and does not rule out the possibility of any change in flux velocity brought about by e.g. protein translocation.
The network altogether consists of 139 reactions processing 108 different chemical species. The full lists of reactions and chemical species is reported in the Additional file 1: Supporting Text.
We assume that the reaction network described above operates at stationarity, i.e. that reaction fluxes in feasible configurations are constant. More precisely, we postulate that the system is kept in a non-equilibrium steady-state (NESS) by the boundary conditions (in our case, by the fluxes of glucose and oxygen into the capillary). Although a steady-state approach for cerebral metabolism will clearly be unable to capture transient or kinetic effects, it can be justified by several considerations. In first place, a typical experiment is performed on tissue volumes containing a large number of cells, and a standard outcome will roughly represent an average over cells in the entire sample. Such averaging can be reasonably approximated with a steady–state assumption, provided the environmental conditions, including stimulation and activation, are stationary. This excludes from the analysis the time intervals associated to the transitions from one state to another (e.g. stimulation onset), which commonly last for a few tenths of a second before a steady–state is attained . Related to this is the fact that in many cases the duration of a stimulus largely exceeds the equilibration time of metabolite concentrations. It has indeed been shown that sustained stimulation induces, after a short transient, a switch to different stationary states for metabolism, neuronal activity, and hemodynamic responses [39–41]. Essentially, the steady-state approach allows for the study of brain metabolism on a time scale lying between the fast adaptation to the change in the activation condition and the slow adjustment of regulatory mechanisms. Finally, within this approach it is possible to treat systems much larger than those accessible to kinetic modeling (see for example ), where only a few nodes of the metabolic networks are usually included.
Constraint-based models provide a standard framework for the analysis of biochemical networks in NESS. In Flux-Balance Analysis (FBA), for instance, one imposes that the vector of concentrations of intracellular metabolites c, which in general would vary in time according to
(where S denotes the M × N stoichiometric matrix, N is the number of reactions, M that of metabolites, ν the vector of reaction fluxes, and b the vector of in- and out-takes that govern the transport of chemical species to and from the system) is constant. In turn, fluxes need to adjust to satisfy simple mass-balance conditions for the individual chemical species, amounting to the set of M equations
Note that the elements of b are non-zero only for metabolites that are exchanged with the environment. For sakes of definiteness, bounds of variability for each flux νi (i = 1,…,N) need to be specified. Usually, such bounds account for reversibility assignment, i.e. they are either of the form -∞ < νi < ∞ (for reversible reactions) or of the form 0 ≤ νi < ∞ (for irreversible reactions), although in some cases physiological considerations may lead to consider more complicated cases, e.g. ν0 ≤ νi < ∞ with ν0 > 0. (In the present study, we shall only consider bounds for the putative reversibility of reactions, as detailed in the network reconstruction reported in the Additional file 1: Supporting text, except for the glucose uptake flux to the capillary which is taken to be fixed. See below for details.)
The system (2) now defines a solution space as a polytope of dimension N - M (or more precisely, N - rank(S); note that, typically, N > M). In the absence of a refinement criterion, like an ad hoc optimization prescription (see e.g.  for an excellent introduction to this modeling perspective), the set of solutions should ideally be sampled uniformly to extract both the individual solutions as well as the statistics of fluxes (averages, correlations, etc.). This is indeed the type of information we are interested in retrieving in the present case. Unluckily, exact sampling algorithms (e.g. Monte Carlo) are still inapplicable to genome-scale flux models when the dimension of the solution space exceeds a few tens because of their high computational costs . (See however  for a promising set of alternative techniques.) In addition, the straightforward application of FBA-type of constraints in the our case is also made difficult by the fact that our reconstruction is largely incomplete. This means that the pathways we do not include may have a considerable cross-talk with the core carbon pathways on which we focus, so that the solutions of (2) might depend strongly on the choice of the boundary fluxes that represent the interaction between pathways included in the reconstruction and the rest of the metabolic network.
We therefore took a step back with respect to FBA and considered a broader type of conditions, inspired by Von Neumann’s model of production networks . In essence we simply replace (2) with
for all intracellular metabolites, while keeping (2) for in-takes (nutrients). Clearly, the main difference with (2) is that (3) allows for flux vectors generating a net production of chemical species (corresponding to the metabolites for which the strict inequality holds in (3)). In a nutshell, the rationale behind this is that a net production of certain compounds might be expected in cells if they need to be employed in macromolecular processes (e.g. proteinogenesis) lying out of the domain of metabolism or, more relevantly to our case, in portions of the network that are not included in the reconstruction. (In other terms, the presence of gaps and their impact on the flux organization of the core network may be smoothened out by softening the constraints.) Therefore, strictly speaking, a NESS where the concentration of certain metabolites is formally increasing in time (corresponding to the positive components of the vector Sν) can be physiologically viable. Once the nutrient availability is fixed through the boundary fluxes, the cell’s metabolic production and nutrient usage profiles can be determined self-consistently from the solutions of (3).
The solutions of (2) and (3) will obviously coincide if all inequalities in (3) become equalities, though in general this does not need to be the case. The main technical advantage of using (3) lies in the existence of an effective and scalable relaxation algorithm that allows to obtain a statistically controlled sampling of its solution space in very modest CPU times. Such a method has been defined in  and employed in e.g.  and  to analyze the metabolic capabilities and the energy balance of the bacterium E. coli. The statistical properties of the solution space sampling thus obtained are discussed in  and further explained in the Additional file 1: Supporting text. In brief, the algorithm makes use of a prior probability distribution of fluxes to initialize the flux variables and generates solutions such that the average Euclidean distance between the solutions and the prior is minimized, the average being carried out over initial flux states sampled from the prior. A sufficiently unbiased prior (e.g. a set of uncorrelated uniform distributions, as we have employed here) then injects minimal a priori information in the solution space and therefore provides a reasonably unbiased, statistically controlled sampling of the feasible flux states of the network. In turn, such an information allows to extract physiological details of individual solutions, as well as statistical properties of the solution space (e.g. probability distributions of fluxes and correlations). This is the calculation scheme we have employed. Further details about the network (i.e. the matrix S), the flux model and the algorithm used in the present work are found in the Additional file 1: Supporting text.
In the following we shall denote the flux of an intracellular reaction or of a transport process respectively by the acronym of the corresponding enzyme or the name of the transported metabolite. We shall also highlight the compartment in which the reaction occurs (so that e.g. νPDH(n) will stand for the pyruvate dehydrogenase–catalyzed reaction taking place in the neuron) or the source/destination compartments that are involved in a transport (for instance, νO2(c → a) will denote the transport of molecular oxygen from the capillary to the astrocyte). Shorthands like νO2(→ c) will instead represent the supply of metabolites (oxygen in this case) to the capillary. Unless otherwise stated, fluxes are expressed in arbitrary units and error bars correspond to one standard error.
Validation of the model: activation states
In the network model we consider, capillaries are supplied with two compounds, namely glucose and oxygen, which can then be transferred to the other compartments. We have characterized the metabolic activity of the brain by fixing only the uptake of glucose to the capillary, νGlc(→ c), while leaving the oxygen influx free. In these conditions, each fixed value of νGlc(→ c) generates a different solution space for (3), where νO2(→ c) fluctuates across solutions (i.e., across feasible flux configurations). This in turn yields, for each selected value of νGlc(→ c), a distribution of values for the oxygen-glucose index (OGI), defined as the ratio between cerebral metabolic rates of oxygen () and glucose (CMRGlc). As both nutrients do not accumulate in the tissue at steady-state, the OGI can be defined as
where nutrient influxes (both the numerator and the denominator) obey, as said above, mass-balance conditions:
Note that the last equation does not involve the neuronal compartment because glucose enters neurons via the extracellular space only (as said before, the endothelium and basal lamina, which mediate the transport of glucose to astrocytes and extracellular space, are included in the capillary compartment). Note also that at steady-state one has
Generically, larger OGIs imply larger fluxes through aerobic pathways, with OGI=6 as the physiologic maximum value for the steady state aerobic oxidation of glucose (corresponding to the fact that 6 oxygen molecules are required to metabolize glucose to water and carbon dioxide). Nevertheless, OGI values larger than 6 are possible in cells whenever the carbon supply for cellular respiration exceeds glucose processing through glycolysis (as happens, for instance, during lactate uptake from the bloodstream). In the flux model (3) it is possible to obtain OGI values slightly above 6, as a consequence of the fact that the oxygen intake to the capillary is a free variable, not bounded (within the model) by the condition OGI ≤ 6. From a modeling viewpoint, this may correspond to a small accumulation of intracellular oxygen due to e.g. the absence of some oxygen-consuming pathways in the network (this condition might not have a physiological counterpart). Experimental in vivo measurements show that the OGI decreases with increasing cerebral activation, from values around 5.5 (almost complete glucose oxidation) under awake resting conditions to values generally ranging from 4 to slightly above 5 during focal brain activity, depending for example on the stimulation paradigm, on the involved brain area, on the experimental procedure (reviewed in [47, 48]). Since the contribution of individual physiologic processes under different conditions is not known, we sought to model the level of activation by using the OGI as a proxy. This assumption stems from the notion that different metabolic states can be characterized in term of their energy expenditure . Furthermore, the transition to a more glycolytic than oxidative metabolism is thought to identify the transition from basal to activated conditions . These arguments suggests that brain metabolism approaches full glucose oxidation as the overall signaling activity decreases.
We have hence solved (3) for different values of νGlc(→ c) recording the resulting OGI distributions. Figure 2 shows four distributions of OGI corresponding to different glucose consumption rates, starting from lower values of νGlc(→ c) corresponding to a larger average OGI. It should be noted that there is not a clear consensus about the quantitative degree of OGI decrease during activity. Indeed, while a value under activation around 5.1 is suggested by several works [47, 48], many others point towards lower values (see e.g. [51, 52] and the recent ). Because of this, and because the OGI distributions we found at a given νGlc(→ c) are rather broad, we preferred to explore a relatively broad range of OGI values.
We have further characterized the model by monitoring the degree of activation in terms of the velocity of the glutamine synthetase (GS) catalyzed reaction. In our simulated network, this univocally represents the rate of the so–called glutamate–glutamine cycle Vcyc, or
It should be noted that, in vivo, there is a residual rate of Gln synthesis unrelated to neurotransmitters cycling , but under physiological conditions νGS(a)≃Vcyc, and thus we conformed to our general choice of neglecting neurotransmission–unrelated amino acids synthesis also in this specific case. The assumption that the residual rate of Gln synthesis is independent of glutamatergic neurotransmission is made also in the original experimental work that reported Vcyc, and the relevant implications are discussed therein.
As happens for the OGI, each choice of νGlc(→ c) leads to a distribution of values of Vcyc. Generically, solutions with Vcyc = 0 (i.e. no neurotransmissiona) will coexist with solutions carrying a non-zero level of activation for any choice of νGlc(→ c). In order to highlight the quantitative changes induced by activation, we performed a correlation analysis between the fluctuations of flow rates of reaction or transport pathways, first in all the sampled solutions (i.e. with Vcyc ≥ 0) and then in solutions with Vcyc > 0. It turns out (see Figure 3) that the transition from uncostrained neurotransmission to presence of neurotransmission has profound consequences on the flux correlations within and between cells. The solutions obtained for Vcyc > 0 underlie a substantial increase of both the global and the regional correlation level of the network. This is especially true for the correlation between neuronal and astrocytic metabolism, confirming that the condition Vcyc > 0 identifies the functional and metabolic partnership between the two cell types. This finding shows that intercellular signaling has a major role for the overall metabolic regulation at tissue level, constraining the catabolism of each cell in a concerted range (see Discussion).
The above results suggest that, in order to capture the quantitative changes that occur in solutions at higher levels of activation (recall that our model does not, per se, constrain the magnitude of neurotransmission), it is useful to analyze the behavior of the conditional average Vcyc versus the conditional average OGI. In Figure 4 we display the results obtained by retaining solutions with Vcyc ≥ 0 (i.e. all solutions, returning the standard average) and Vcyc > 0, respectively. One sees that as the average OGI decreases, the average Vcyc increases to an enhanced activation level. The slope of the curve is less negative for the condition Vcyc > 0, i.e. when the sample is restricted by filtering–out the states with suppressed neurotransmission. This changes the basal state at OGI = 5.5 to a rate of glutamate-glutamine cycle larger by about 30%, from 0.06 to 0.08 (roughly).
In summary, in agreement with the literature, we find that the rise in anaerobic glucose consumption during activation turns out to be out of proportion to oxygen utilization, as evidenced by the decrease in the (average) OGI for increasing values of Vcyc and CMRGlc. This indicates that the glutamate-glutamine cycle by itself suffices to explain part of the uncoupling between glucose and oxygen utilization. It should be emphasized that the reported OGI reduction, although potentially significant if caused by a specific subset of glycolytically-served energetic demands, does not change the overall strategy of brain energy metabolism, which remains largely aerobic because of the higher ATP yield of respiration .
Neuronal oxidative glucose metabolism versus glutamatergic activity
The relationship between glutamate-glutamine cycle and neuronal glucose oxidation was experimentally reported to be close to a 1:1 relation . In particular, the rate of neuronal oxidative metabolism of glucose (i.e. the level of activity of pyruvate dehydrogenase, PDH) increases linearly with the velocity of transmitter cycling, with a slope close to one. Our framework allows to address the dependence of glucose oxidation in neurons on the rate of glutamate-glutamine cycle for various states of activation. Considering that PDH is the primary entry point of glucose-derived pyruvate into the TCA cycle, we define
where the factor 1/2 is required as glycolysis produces two molecules of pyruvate for each glucose molecule. Plotting the average CMRGlc(ox)(n) against the average Vcyc we find approximately two different regimes around the physiologic range corresponding to OGI ≃ 5.5 (see Figure 5). These regimes are characterized by almost linear relations, in agreement with the experimentally reported constant stoichiometry between aerobic Glc oxidation in neurons and glutamate cycling [54, 56]. However, the coupling pattern changes as one explores states with values of Vcyc departing from the basal level to higher activity. The slope of the curve clearly increases if only strictly positive Vcyc are considered, although it stays below one. Interestingly, in the high Vcyc regime (upstream the awake value, for which no experimental data exist) the plot features a slight attenuation of the curve. The fact that the slope is lower than the linear extrapolation at high activity indicates that some source of energy consumption adds to the glutamate–glutamine cycle during the transition from normal to high neurotransmission levels. In this region, the discrepancy between simulated and experimental data is explained by the relatively low increase in ionic fluxes through neuronal voltage-gated Na+ and K+ channels obtained in our simulations (data not shown), suggesting that the glutamate–glutamine cycle alone is insufficient to account for the rise in brain glucose utilization. To this end, it is mandatory for stoichiometric models to incorporate energy use by action and synaptic potentials in addition to glutamatergic neurotransmission.
Modulation of glucose uptake and cell-to-cell lactate shuttling by glutamate-glutamine cycle
The fate of carbons undergoing oxidative phosphorylation in neurons and astrocytes provides a quantitative hint of the relative amount of energy produced by aerobic pathways in the two cell types. We found that anaerobic and aerobic metabolism is similarly increased in both cell types at increasing activity (not shown). In particular, the fraction of cerebral oxidative metabolism in astrocytes is about 35% of the total, which is consistent with a substantial astrocytic contribution to functional brain energy metabolism (reviewed in ).
Analysis of glucose fluxes showed large fluctuations, so that both the neuron and the astrocyte may be the primary sites of glucose consumption at fixed values of Glu/Gln cycling. We indeed observed that both cell compartments can absorb from 10% to 90% of the total glucose uptake respectively. Therefore, at this level of detail, the network does not place significant constraints on the cellular utilization of glucose. The failure of up-regulation of ionic fluxes that we report here might play a role if their contribution is substantially different for neurons and astrocytes, which unfortunately has not been yet experimentally determined.
The direction of the shuttle of lactate depends in a robust way on the sharing of glucose between neurons and astrocytes, resulting in ANLS when the relative astrocytic glucose uptake becomes larger than about 65% (implying that states supporting ANLS can be obtained by an ad hoc adjustment of Glc partitioning) and NALS otherwise (Figure 6). As the latter is also the mean value for the fraction of neuronal versus astrocytic oxygen utilization, it turns out that, on average, the contribution of CCLS to cell metabolism is very low compared to the lactate generated intracellularly by glucose. The strict dependence of CCLS on cellular glucose uptake supports previous results obtained through various modeling approaches [11, 18, 19].
A closer look at the individual solutions reveals the absence of a significant correlation between the glutamate–glutamine cycle and both the uptake of glucose and the shuttle of lactate (Figure 7). It should however be stressed that, as can be seen in Figure 7, the emerging scenario presents large fluctuations, in the sense that, even within the physiological range for the OGI, solutions with ANLS and NALS coexist. Together, these simulation outcomes indicate that no preferential route is undertaken by lactate at any given rate of glutamate-glutamine cycle. Thus, the determinants for lactate accumulation and shuttling, if any, must reside elsewhere, for example in the balance between spiking and synaptic activity .
Conserved moieties, transcellular aspartate shuttling and glutathione-ascorbate cycle
To conclude, we discuss the flux organization of several pathways that were included in the metabolic network reconstruction, but whose involvement is not strictly related to lactate trafficking. The ATP buffering systems of creatine and adenylate kinase in neurons and astrocytes, as well as astrocytic glycogen metabolism were found to have negligible net flux, consistently with the role of adenylates, creatine and glycogen as conserved moieties in the steady-state. We are unable to test hypotheses about the role of brain glycogen [59–61], because our approach does not allow to describe the consumption of previously stored metabolites, such as glycogen in this case. We found that the contribution to NADH shuttling by the malate-aspartate shuttle (MAS) dominates over the glycerol-3-phosphate pathway (relative flux is a few percent). We could not support the recent hypothesis of a significant steady state neuron-to-astrocyte transport of aspartate . However, this finding was expected due to the presence of mitochondrial aspartate-glutamate carrier in astrocytes, which allows for a shunting pathway, alternative to neuron-to-astrocyte aspartate transport that is sufficient to sustain NADH shuttling from cytosol to mitochondria in astrocytes, as previously suggested . In order to examine a possible role of (dehydro)ascorbate transfer between neurons and astrocytes , we included the detoxification of reactive oxygen species (ROS) via glutathione/ascorbate cycle. Unluckily, we were unable to find statistically significant exchange fluxes of the two forms of ascorbic acids between cells. This is possibly correlated with the fact that, as said above, we did not include the entire network section associated with glutathione synthesis. Theoretical analysis of transcellular ascorbate cycling will thus require further refinements of the network reconstruction.
This work is concerned with a model of brain energy metabolism consisting of four compartments, representing neuron, astrocyte, extracellular space and the capillary. We used a steady state approach based on Von Neumann’s theoretical framework for the analysis of production networks , successfully applied to cellular metabolism in previous studies . The steady-state assumption for the metabolic coupling between oxygen and glucose consumption underlying different cortical states is consistent with (i) the fact that experimental parameters have been measured under stationary conditions during suppressed brain activity , and (ii) the establishment of a new metabolic steady-state during enhanced brain activity . At odds with standard schemes based on flux-balance analysis though, the frame employed here doesn’t constrain the net metabolite production to zero, but, rather, aims at recovering the steady state self-consistently from minimal stability requirements. The 'soft’ type of constraints thus arising makes it possible to sample the solution space corresponding to our large-scale network model in a statistically controlled manner, and returns a full range of feasible values for each reaction flux in the network, as well as detailed information on correlations. Based on this, one can elucidate the extent to which stoichiometry alone constrains the operation of brain metabolism, since the emerging picture is obtained without imposing specific functional constraints.
First, we found that the rate of glutamate-glutamine cycle distinguishes different activation states according to the fraction of glucose that is processed via glycolysis versus respiration. In particular, decreasing OGI values predict increasing velocity of the cycle, and viceversa. Notably, ATP-consuming Na + and K + fluxes across voltage- and ligand-gated ion channels are not directly dependent on glutamate-glutamine cycle in our theoretical account. Pathways analysis showed that their reaction rates do not “automatically” up-regulate along with transmitter cycling. This indicates that the absence of causal changes (i.e. constraints) in these different aspects of neuronal signaling strongly underestimates the glucose utilization at high activity levels. Unfortunately, it is presently unknown how to model the exact cause-and-effect relationship between glutamate-glutamine cycle and the ionic currents that generated neurotransmission on one hand, and those that are generated by neurotransmission on the other. The latter will constitute a primary target of future studies.
Second, the results of the present model support previous kinetic analyses indicating that lactate derived from astrocytes may not provide an important source of carbon compounds for neuronal energy metabolism in an activation-dependent manner (Figure 6) . The model also supports the conclusion that the direction and magnitude of CCLS are secondary to glucose partitioning between cell types [11, 19]. We found that glucose uptake, glycolysis and respiration change proportionally in a wide range of glutamate recycling rates. This implies that the glucose taken up by each cell is completely metabolized and not significantly converted to lactate. Thus, lactate does not accumulate in a specific cell type, i.e. the lactate concentration gradient, and hence lactate shuttle, remains small.
Third, simulations showed that glutamate-glutamine cycle is correlated with overall tissue glucose utilization and lactate production, but not with specific patterns of cellular glucose uptake and lactate shuttle (Figure 7). These results agree well with the experimental knowledge [54, 64] and recover features of different (and sometimes conflicting) numerical studies performed previously [17, 19]. Most importantly, they add further arguments to the idea that the trafficking of molecules between neurons and astrocytes underlies a broad functional, rather than strictly energetic partnership (see, e.g. ). This is evident in our modeling perspective, as possibly energy-related changes in lactate fluxes are independent of concomitant function-related variations in glutamate-glutamine cycle (see Figure 7, that essentially shows a lack of correlation between lactate shuttling and transmitter cycling). Alternative functions for lactate include the discrimination between dilation and constriction of cerebral arterioles , or the modulation of GABAergic inhibitory activity of specific neuronal populations . These functional roles of lactate might be still secondary to its accumulation in the tissue. Our simulations suggest that this accumulation likely results from up-regulation of non-oxidative metabolism in both neurons and astrocytes, as previously suggested [11, 60].
We found that intercellular shuttling of lactate increases only when glucose partitioning between neurons and astrocytes is significantly uneven (Figure 6). Therefore, lactate transfer can be interpreted as a local biochemical shunt that allows for the optimal use of carbon supply in correspondence of variable environmental challenges.
Interestingly, if the functional partnership between neurons and astrocytes is suppressed via the zeroing of neurotransmission, the anticorrelation between glycolysis in the two cellular compartments further increases (i.e. correlation becomes more negative), a feature shared with many other fluxes related to energetics, including TCA cycle and oxidative phosphorylation. So, even when there is no functional relationship between neurons and astrocytes, their energy (primarily glycolytic) metabolism is anticorrelated. The reason for this behavior is that the two cell types share the same environment, thereby the same glucose availability. Glucose availability thus represents the primary drive on glucose partitioning. If the two cells are functionally independent, oxidative metabolism follows glycolysis, as also evidenced by the negative correlation between neuronal and astrocytic oxidative metabolism. As soon as neurons and astrocytes become coupled by neurotransmission, the network shifts to correlated patterns of activity, both within the same compartment and between different compartments. Accordingly, glycolysis in one cell type becomes slightly positively correlated with oxidative metabolism in the other. However, glycolysis in either neurons or astrocytes, which essentially reflect partitioning of tissue glucose, remains poorly correlated with the glutamate-glutamine cycle. Overall, glutamate-glutamine cycle positively correlates with cell TCA cycle and oxidative phosphorylation much more than it does with glycolysis, which is especially significant for the astrocytic compartment.
In conclusion, we developed a large-scale model for compartmentalized brain energy metabolism including the core carbon pathways of neurons and astrocytes, as well as further compartments (capillary, extracellular space). A constraint-based scheme was then employed in order to define feasible configurations of reaction fluxes. Our numerical analysis was based on a relaxation algorithm which allows to obtain a statistically controlled sampling of the solution space without any prior assumption on the behaviour of energy producing/consuming pathways. Results have shown that only a large imbalance of cell glucose uptake can explain the occurrence of a significant lactate shuttle between neurons and astrocytes, the latter being roughly independent of the rates of transmitter cycling. Our results therefore do not support a link between glutamate-glutamine cycle and CCLS as a mechanism for activation-dependent transfer of energy compounds within the brain. On the other hand, CCLS can be found by assuming that neurons have limited access to glucose and/or by bounding their glucose uptake flux. The lack of correlation we observe stems from the fact that the distribution of activity-related energy stress between neurons and astrocytes cannot be estimated by the stoichiometry of the metabolic network and should be a primary target of current experimental research. Future developments will focus on introducing minimal ad hoc constraints on neurotransmitter cycling and ionic fluxes in the hope to capture the non-stoichiometric side of the energetics of brain activity.
a Because of finite numerical precision, for our putposes a flux is null when it is smaller than 10-6.
This work is supported by the DREAM Seed Project of the Italian Institute of Technology (IIT) and by the joint IIT/Sapienza Lab “Nanomedicine”. The IIT Platform “Computation” is gratefully acknowledged. FAM acknowledges financial support from European Union Grants PIRG-GA-2010-277166 and PIRG-GA-2010-268342.
Departament d’Enginyeria Quimica, Universitat Rovira i Virgili
Magnetic Resonance for Brain Investigation Lab, Enrico Fermi Center
Dipartimento di Fisica, Sapienza Università di Roma
Fondazione Santa Lucia
Department of Mathematics, King’s College London
Center for Life Nano Science@Sapienza, Istituto Italiano di Tecnologia
CNR–IPCF, Unità di Roma Sapienza
Mangia S, Giove F, Tkác I, Logothetis NK, Olman CA, Maraviglia B, Di Salle F, Ugurbil K: Metabolic and hemodynamic events after changes in neuronal activity: current hypotheses, theoretical predictions and in vivo NMR experimental findings.J Cereb Blood Flow Metab 2009, 29:441–463.PubMedView Article
Pellerin L, Magistretti P: Glutamate uptake into astrocytes stimulates aerobic glycolysis: a mechanism coupling neuronal activity to glucose utilization.Proc Natl Acad Sci USA 1994, 91:10625–10629.PubMedView Article
Dienel GA, Cruz NF: Astrocyte activation in working brain: energy supplied by minor substrates.Neurochem Int 2006, 48:586–595.PubMedView Article
Bak LK, Schousboe A, Sonnewald U, Waagepetersen HS: Glucose is necessary to maintain neurotransmitter homeostasis during synaptic activity in cultured glutamatergic neurons.J Cereb Blood Flow Metab 2006, 26:1285–1297.PubMedView Article
Dienel GA: Astrocytic energetics during excitatory neurotransmission: What are contributions of glutamate oxidation and glycolysis?Neurochem Int 2013, 63:244–258.PubMedView Article
Mergenthaler P, Lindauer U, Dienel GA, Meisel A: Sugar for the brain: the role of glucose in physiological and pathological brain function (2013).Trends Neurosci 2013,. in press
Aubert A, Costalat R: Interaction between astrocytes and neurons studied using a mathematical model of compartmentalized energy metabolism.J Cereb Blood Flow Metab 2005, 25:1476–1490.PubMedView Article
Aubert A, Pellerin L, Magistretti PJ, Costalat R: A coherent neurobiological framework for functional neuroimaging provided by a model integrating compartmentalized energy metabolism.Proc Natl Acad Sci USA 2007, 104:4188–4193.PubMedView Article
Simpson IA, Carruthers A, Vannucci SJ: Supply and demand in cerebral energy metabolism: the role of nutrient transporters.J Cereb Blood Flow Metab 2007, 27:1766–1791.PubMedView Article
Mangia S, Simpson IA, Vannucci SJ, Carruthers A: The in vivo neuron-to-astrocyte lactate shuttle in human brain: evidence from modeling of measured lactate levels during visual stimulation.J Neurochem 2009,109(Suppl 1):55–62.PubMedView Article
DiNuzzo M, Mangia S, Maraviglia B, Giove F: Changes in glucose uptake rather than lactate shuttle take center stage in subserving neuroenergetics: evidence from mathematical modeling.J Cereb Blood Flow Metab 2010, 30:586–602.PubMedView Article
Somersalo E, Cheng Y, Calvetti D: The metabolism of neurons and astrocytes through mathematical models.Ann Biomed Eng 2012, 40:2328–2344.PubMedView Article
Calvetti D, Somersalo E: Dynamic activation model for a glutamatergic neurovascular unit.J Theor Biol 2011, 274:12–29.PubMedView Article
Jolivet R, Magistretti PJ, Weber B: Deciphering neuron-glia compartmentalization in cortical energy metabolism.Front Neuroenergetics 2009, 1:4.PubMedView Article
Occhipinti R, Somersalo E, Calvetti D: Astrocytes as the glucose shunt for glutamatergic neurons at high activity: an in silico study.J Neurophysiol 2009, 101:2528–2538.PubMedView Article
Dienel GA: Brain lactate metabolism: the discoveries and the controversies.J Cereb Blood Flow Metab 2012, 32:1107–1138.PubMedView Article
Jolivet R, Allaman I, Pellerin L, Magistretti PJ, Weber B: Comment on recent modeling studies of astrocyte-neuron metabolic interactions.J Cereb Blood Flow Metab 2010, 30:1982–1986.PubMedView Article
Calvetti D, Somersalo E: Ménage à trois: the role of neurotransmitters in the energy metabolism of astrocytes, glutamatergic, and GABAergic neurons.J Cereb Blood Flow Metab 2012, 32:1472–1783.PubMedView Article
Mangia S, DiNuzzo M, Carruthers A, Simpson IA, Vannucci SJ, GioveF: Response to 'Comment on recent modeling studies of astrocyte-neuron metabolic interactions’: much ado about nothing.J Cereb Blood Flow Metab 2011, 31:1346–1353.PubMedView Article
Martelli C, De Martino A, Marinari E, Marsili M, Castillo IP: Identifying essential genes in escherichia coli from a metabolic optimization principle.Proc Natl Acad Sci USA 2009, 106:2607–2611.PubMedView Article
De Martino D, Figliuzzi M, De Martino A, Marinari E: A scalable algorithm to explore the gibbs energy landscape of genome-scale metabolic networks.PLoS Comput Biol 2012, 8:e1002562.PubMedView Article
Tunahan Ç, Selma A, Hale S, Ata A, Kutlu Ü: Reconstruction and flux analysis of coupling between metabolic pathways of astrocytes and neurons: application to cerebral hypoxia.Theor Biol Med Model 2007, 4:48.View Article
Price ND, Schellenberger J, Palsson BØ: Uniform sampling of steady-state flux spaces: means to design experiments and to interpret enzymopathies.Biophys J 2004, 87:2172–2186.PubMedView Article
Hyder F, Shulman RG, Rothman DL: A model for the regulation of cerebral oxygen delivery.J Appl Physiol 1998, 85:554–564.PubMed
Gjedde A: The relation between brain function and cerebral blood flow and metabolism. In Cerebrovascular Disease. Edited by: Batjer HH. Philadelphia: Lippincott–Raven; 1997:23–40.
Lovatt D, Sonnewald U, Waagepetersen HS, Schousboe A, He W, Lin JH, Han X, Takano T, Wang S, Sim FJ, Goldman SA, Nedergaard M: The transcriptome and metabolic gene signature of protoplasmic astrocytes in the adult murine cortex.J Neurosci 2007, 27:12255–12266.PubMedView Article
Dringen R, Pfeiffer B, Hamprecht B: Synthesis of the antioxidant glutathione in neurons: supply by astrocytes of cysgly as precursor for neuronal glutathione.J Neurosci 1999, 19:562–569.PubMed
Gegg ME, Beltran B, Salas-Pino S, Bolanos JP, Clark JB, Moncada S, Heales SJ: Differential effect of nitric oxide on glutathione metabolism and mitochondrial function in astrocytes and neurones: implications for neuroprotection/neurodegeneration?J Neurochem 2003, 86:228–237.PubMedView Article
Yu AC, Drejer J, Hertz L, Schousboe A: Pyruvate carboxylase activity in primary cultures of astrocytes and neurons.J Neurochem 1983, 41:1484–1487.PubMedView Article
Huang YH, Bergles DE: Glutamate transporters bring competition to the synapse.Curr Opin Neurobiol 2004, 14:346–352.PubMedView Article
Mangia S, Giove F, DiNuzzo M: Metabolic pathways and activity-dependent modulation of glutamate concentration in the human brain.Neurochem Res 2012, 37:2554–2561.PubMedView Article
Waagepetersen HS, Sonnewald U, Larsson OM, Schousboe A: A possible role of alanine for ammonia transfer between astrocytes and glutamatergic neurons.J Neurochem 2000, 75:471–479.PubMedView Article
Bak LK, Sickmann HM, Schousboe A, Waagepetersen HS: Activity of the lactate-alanine shuttle is independent of glutamate-glutamine cycle activity in cerebellar neuronal-astrocytic cultures.J Neurosci Res 2005, 79:88–96.PubMedView Article
Mangia S, Tkác I, Gruetter R, Van de Moortele PF, Maraviglia B, Ugurbil K: Sustained neuronal activation raises oxidative metabolism to a new steady-state level: evidence from (1)H NMR spectroscopy in the human visual cortex.J Cereb Blood Flow Metab 2007, 27:1055–1063.PubMed
Bandettini PA, Kwong KK, Davis TL, Tootell RB, Wong EC, Fox PT, Belliveau JW, Weisskoff RM, Rosen BR: Characterization of cerebral blood oxygenation and flow changes during prolonged brain activation.Hum Brain Mapp 1997, 5:93–109.PubMedView Article
Logothetis NK, Pauls J, Augath M, Trinath T, Oeltermann A: Neurophysiological investigation of the basis of the fMRI signal.Nature 2001, 412:150–157.PubMedView Article
Palsson BØ: Systems Biology: Properties of Reconstructed Networks. Cambridge: Cambridge University Press; 2006.View Article
Schellenberger J, Palsson B: Use of randomized sampling for analysis of metabolic networks.J Biol Chem 2009, 284:5457–5461.PubMedView Article
Braunstein A, Mulet R, Pagnani A: Estimating the size of the solution space of metabolic networks.BMC Bioinformatics 2008, 9:240.PubMedView Article
Gale D: The Theory of Linear Economic Models. Chicago: The University of Chicago Press; 1989.
De Martino A, Martelli C, Monasson R, Perez Castillo I: Von Neumann’s expanding model on random graphs.J Stat Mech 2007, 2007:P05012.View Article
Shulman RG, Hyder F, Rothman DL: Lactate efflux and the neuroenergetic basis of brain function.NMR Biomed 2001, 14:389–396.PubMedView Article
Giove F, Mangia S, Bianciardi M, Garreffa G, Di Salle F, Morrone R, Maraviglia B: The physiology and metabolism of neuronal activation: in vivo studies by NMR and other methods.Magn Reson Imaging 2003, 21:1283–1293.PubMedView Article
Shulman RG: A philosophical analysis of neuroenergetics.Front Neuroenergetics 2011, 3:6.PubMedView Article
Raichle ME, Gusnard DA: Appraising the brain’s energy budget.Proc Natl Acad Sci USA 2002, 99:10237–10239.PubMedView Article
Fox PT, Raichle ME: Focal physiological uncoupling of cerebral blood flow and oxidative metabolism during somatosensory stimulation in human subjects.Proc Natl Acad Sci USA 1986, 83:1140–1144.PubMedView Article
Fox PT, Raichle ME, Mintun MA, Dence C: Nonoxydative glucose consumption during focal physiologic neural activity.Science 1988, 241:462–464.PubMedView Article
Vafaee MS, Vang K, Bergersen LH, Gjedde A: Oxygen consumption and blood flow coupling in human motor cortex during intense finger tapping: implication for a role of lactate.J Cereb Blood Flow Metab 2012, 32:1859–1868.PubMedView Article
Sibson NR, Dhankhar A, Mason GF, Rothman DL, Behar KL, Shulman RG: Stoichiometric coupling of brain glucose metabolism and glutamatergic neuronal activity.Proc Natl Acad Sci USA 1998, 95:316–321.PubMedView Article
Sibson NR, Dhankhar A, Mason GF, Rothman DL, Behar KL, Shulman RG: In vivo 13C NMR measurements of cerebral glutamine synthesis as evidence for glutamate-glutamine cycling.Proc Natl Acad Sci USA 1997, 94:2699–2704.PubMedView Article
Hyder F, Patel AB, Gjedde A, Rothman DL, Behar KL, Shulman RG: Neuronal-glial glucose oxidation and glutamatergic-GABAergic function.J Cereb Blood Flow Metab 2006, 26:865–877.PubMedView Article
Hertz L, Peng L, Dienel GA: Energy metabolism in astrocytes: high rate of oxidative metabolism and spatiotemporal dependence on glycolysis/glycogenolysis.J Cereb Blood Flow Metab 2007, 27:219–249.PubMedView Article
DiNuzzo M, Giove F: Activity-dependent energy budget for neocortical signaling: effect of short-term synaptic plasticity on the energy expended by spiking and synaptic activity.J Neurosci Res 2012, 90:2094–2102.PubMedView Article
Shulman RG, Hyder F, Rothman DL: Cerebral energetics and the glycogen shunt: neurochemical basis of functional imaging.Proc Natl Acad Sci USA 2001, 98:6417–6422.PubMedView Article
DiNuzzo M, Mangia S, Maraviglia B, Giove F: Glycogenolysis in astrocytes supports blood-borne glucose channeling not glycogen-derived lactate shuttling to neurons: evidence from mathematical modeling.J Cereb Blood Flow Metab 2010, 30:1895–1904.PubMedView Article
DiNuzzo M, Maraviglia B, Giove F: Why does the brain (not) have glycogen?Bioessays 2011, 33:319–326.PubMedView Article
Pardo B, Rodrigues TB, Contreras L, Garzon M, Llorente-Folch I, Kobayashi K, Saheki T, Cerdan S, Satrústegui J: Brain glutamine synthesis requires neuronal-born aspartate as amino donor for glial glutamate formation.J Cereb Blood Flow Metab 2011, 31:90–101.PubMedView Article
Castro MA, Beltrán FA, Brauchi S, Concha II: A metabolic switch in brain: glucose and lactate metabolism modulation by ascorbic acid.J Neurochem 2009, 110:423–440.PubMedView Article
Bak L, Walls A, Schousboe A, Ring A, Sonnewald U, Waagepetersen HS: Neuronal glucose but not lactate utilization is positively correlated with NMDA-induced neurotransmission and fluctuations in cytosolic Ca2+ levels.J Neurochem 2009, 109:87–93.PubMedView Article
Bergersen LH, Gjedde A: Is lactate a volume transmitter of metabolic states of the brain?Front Neuroenerg 2012, 4:5.View Article
Gordon GRJ, Choi HB, Rungta RL, Ellis-Davies GCR, MacVicar BA: Brain metabolism dictates the polarity of astrocyte control over arterioles.Nature 2008, 456:745–749.PubMedView Article
Shimizu H, Watanabe E, Hiyama TY, Nagakura A, Fujikawa A, Okado H, Yanagawa Y, Obata K, Noda M: Glial Nax channels control lactate signaling to neurons for brain [Na+] sensing.Neuron 2007, 54:59–72.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.