Energy metabolism and glutamate-glutamine cycle in the brain: a stoichiometric modeling perspective
- Francesco A Massucci1,
- Mauro DiNuzzo2, 3,
- Federico Giove2, 3,
- Bruno Maraviglia3, 4,
- Isaac Perez Castillo5,
- Enzo Marinari†3, 6 and
- Andrea De Martino†3, 6, 7Email author
© Massucci et al.; licensee BioMed Central Ltd. 2013
Received: 16 May 2013
Accepted: 30 September 2013
Published: 10 October 2013
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.
KeywordsBrain energetics Lactate shuttle Metabolic modeling Glutamate-glutamine cycle Glucose partitioning OGI
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 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 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.
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.
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
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.
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.
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
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.
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.
- 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. 10.1038/jcbfm.2008.134.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1073/pnas.91.22.10625.PubMedPubMed CentralView ArticleGoogle Scholar
- Dienel GA, Cruz NF: Astrocyte activation in working brain: energy supplied by minor substrates. Neurochem Int. 2006, 48: 586-595. 10.1016/j.neuint.2006.01.004.PubMedView ArticleGoogle Scholar
- 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. 10.1038/sj.jcbfm.9600281.PubMedView ArticleGoogle Scholar
- Dienel GA: Astrocytic energetics during excitatory neurotransmission: What are contributions of glutamate oxidation and glycolysis?. Neurochem Int. 2013, 63: 244-258. 10.1016/j.neuint.2013.06.015.PubMedPubMed CentralView ArticleGoogle Scholar
- 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 pressGoogle Scholar
- 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. 10.1038/sj.jcbfm.9600144.PubMedView ArticleGoogle Scholar
- 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. 10.1073/pnas.0605864104.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1038/sj.jcbfm.9600521.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1038/jcbfm.2009.232.PubMedPubMed CentralView ArticleGoogle Scholar
- Somersalo E, Cheng Y, Calvetti D: The metabolism of neurons and astrocytes through mathematical models. Ann Biomed Eng. 2012, 40: 2328-2344. 10.1007/s10439-012-0643-z.PubMedView ArticleGoogle Scholar
- Calvetti D, Somersalo E: Dynamic activation model for a glutamatergic neurovascular unit. J Theor Biol. 2011, 274: 12-29. 10.1016/j.jtbi.2010.12.007.PubMedView ArticleGoogle Scholar
- Jolivet R, Magistretti PJ, Weber B: Deciphering neuron-glia compartmentalization in cortical energy metabolism. Front Neuroenergetics. 2009, 1: 4-PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Dienel GA: Brain lactate metabolism: the discoveries and the controversies. J Cereb Blood Flow Metab. 2012, 32: 1107-1138. 10.1038/jcbfm.2011.175.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1038/jcbfm.2010.132.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1038/jcbfm.2012.31.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1038/jcbfm.2011.29.PubMedPubMed CentralView ArticleGoogle Scholar
- Dienel GA: Fueling and imaging brain activation. ASN Neuro ASN Neuro. 2012, 4: e00093-PubMedGoogle Scholar
- 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. 10.1073/pnas.0813229106.PubMedPubMed CentralView ArticleGoogle Scholar
- 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-10.1371/journal.pcbi.1002562.PubMedPubMed CentralView ArticleGoogle Scholar
- 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-10.1186/1742-4682-4-48.View ArticleGoogle Scholar
- 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. 10.1529/biophysj.104.043000.PubMedPubMed CentralView ArticleGoogle Scholar
- Hyder F, Shulman RG, Rothman DL: A model for the regulation of cerebral oxygen delivery. J Appl Physiol. 1998, 85: 554-564.PubMedGoogle Scholar
- Gjedde A: The relation between brain function and cerebral blood flow and metabolism. Cerebrovascular Disease. Edited by: Batjer HH. 1997, Philadelphia: Lippincott–Raven, 23-40.Google Scholar
- 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. 10.1523/JNEUROSCI.3404-07.2007.PubMedView ArticleGoogle Scholar
- Hertz L: Brain glutamine synthesis requires neuronal aspartate: a commentary. J Cereb Blood Flow Metab. 2011, 31: 384-387. 10.1038/jcbfm.2010.199.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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 ArticleGoogle Scholar
- Yu AC, Drejer J, Hertz L, Schousboe A: Pyruvate carboxylase activity in primary cultures of astrocytes and neurons. J Neurochem. 1983, 41: 1484-1487. 10.1111/j.1471-4159.1983.tb00849.x.PubMedView ArticleGoogle Scholar
- Hassel B, Bråthe A: Neuronal pyruvate carboxylation supports formation of transmitter glutamate. J Neurosci. 2000, 20: 1342-1347.PubMedGoogle Scholar
- Ames A: CNS energy metabolism as related to function. Brain Res Brain Res Rev. 2000, 34: 42-68. 10.1016/S0165-0173(00)00038-2.PubMedView ArticleGoogle Scholar
- Danbolt NC: Glutamate uptake. Prog Neurobiol. 2001, 65: 1-105. 10.1016/S0301-0082(00)00067-8.PubMedView ArticleGoogle Scholar
- Huang YH, Bergles DE: Glutamate transporters bring competition to the synapse. Curr Opin Neurobiol. 2004, 14: 346-352. 10.1016/j.conb.2004.05.007.PubMedView ArticleGoogle Scholar
- 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. 10.1007/s11064-012-0848-4.PubMedPubMed CentralView ArticleGoogle Scholar
- 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 ArticleGoogle Scholar
- 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. 10.1002/jnr.20319.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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. 10.1002/(SICI)1097-0193(1997)5:2<93::AID-HBM3>3.0.CO;2-H.PubMedView ArticleGoogle Scholar
- Logothetis NK, Pauls J, Augath M, Trinath T, Oeltermann A: Neurophysiological investigation of the basis of the fMRI signal. Nature. 2001, 412: 150-157. 10.1038/35084005.PubMedView ArticleGoogle Scholar
- Palsson BØ: Systems Biology: Properties of Reconstructed Networks. 2006, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Schellenberger J, Palsson B: Use of randomized sampling for analysis of metabolic networks. J Biol Chem. 2009, 284: 5457-5461.PubMedView ArticleGoogle Scholar
- Braunstein A, Mulet R, Pagnani A: Estimating the size of the solution space of metabolic networks. BMC Bioinformatics. 2008, 9: 240-10.1186/1471-2105-9-240.PubMedPubMed CentralView ArticleGoogle Scholar
- Gale D: The Theory of Linear Economic Models. 1989, Chicago: The University of Chicago PressGoogle Scholar
- De Martino A, Martelli C, Monasson R, Perez Castillo I: Von Neumann’s expanding model on random graphs. J Stat Mech. 2007, 2007: P05012-10.1088/1742-5468/2007/05/P05012.Google Scholar
- Shulman RG, Hyder F, Rothman DL: Lactate efflux and the neuroenergetic basis of brain function. NMR Biomed. 2001, 14: 389-396. 10.1002/nbm.741.PubMedView ArticleGoogle Scholar
- 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. 10.1016/j.mri.2003.08.028.PubMedView ArticleGoogle Scholar
- Shulman RG: A philosophical analysis of neuroenergetics. Front Neuroenergetics. 2011, 3: 6-PubMedPubMed CentralView ArticleGoogle Scholar
- Raichle ME, Gusnard DA: Appraising the brain’s energy budget. Proc Natl Acad Sci USA. 2002, 99: 10237-10239. 10.1073/pnas.172399499.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1073/pnas.83.4.1140.PubMedPubMed CentralView ArticleGoogle Scholar
- Fox PT, Raichle ME, Mintun MA, Dence C: Nonoxydative glucose consumption during focal physiologic neural activity. Science. 1988, 241: 462-464. 10.1126/science.3260686.PubMedView ArticleGoogle Scholar
- 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. 10.1038/jcbfm.2012.89.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1073/pnas.95.1.316.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1073/pnas.94.6.2699.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1038/sj.jcbfm.9600263.PubMedView ArticleGoogle Scholar
- 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. 10.1038/sj.jcbfm.9600343.PubMedView ArticleGoogle Scholar
- 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. 10.1002/jnr.23098.PubMedView ArticleGoogle Scholar
- 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. 10.1073/pnas.101129298.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1038/jcbfm.2010.151.PubMedPubMed CentralView ArticleGoogle Scholar
- DiNuzzo M, Maraviglia B, Giove F: Why does the brain (not) have glycogen?. Bioessays. 2011, 33: 319-326. 10.1002/bies.201000151.PubMedView ArticleGoogle Scholar
- 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. 10.1038/jcbfm.2010.146.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1111/j.1471-4159.2009.06151.x.PubMedView ArticleGoogle Scholar
- 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 ArticleGoogle Scholar
- Bergersen LH, Gjedde A: Is lactate a volume transmitter of metabolic states of the brain?. Front Neuroenerg. 2012, 4: 5-View ArticleGoogle Scholar
- 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. 10.1038/nature07525.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1016/j.neuron.2007.03.014.PubMedView ArticleGoogle Scholar