Computational singular perturbation analysis of brain lactate metabolism
Authors:
Dimitris G. Patsatzis aff001; Efstathios-Al. Tingas aff001; Dimitris A. Goussis aff004; S. Mani Sarathy aff001
Authors place of work:
King Abdullah University of Science and Technology (KAUST), Clean Combustion Research Center (CCRC), Thuwal, Saudi Arabia
aff001; Department of Mechanics, School of Applied Mathematics and Physical Sciences, National Technical University of Athens (NTUA), Athens, Greece
aff002; Perth College, University of the Highlands and Islands, Crieff Rd, Perth PH1 2NX, United Kingdom
aff003; Department of Mechanical Engineering, Khalifa University of Science, Technology and Research (KUSTAR), Abu Dhabi, United Arab Emirates
aff004
Published in the journal:
PLoS ONE 14(12)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0226094
Summary
Lactate in the brain is considered an important fuel and signalling molecule for neuronal activity, especially during neuronal activation. Whether lactate is shuttled from astrocytes to neurons or from neurons to astrocytes leads to the contradictory Astrocyte to Neuron Lactate Shuttle (ANLS) or Neuron to Astrocyte Lactate Shuttle (NALS) hypotheses, both of which are supported by extensive, but indirect, experimental evidence. This work explores the conditions favouring development of ANLS or NALS phenomenon on the basis of a model that can simulate both by employing the two parameter sets proposed by Simpson et al. (J Cereb. Blood Flow Metab., 27:1766, 2007) and Mangia et al. (J of Neurochemistry, 109:55, 2009). As most mathematical models governing brain metabolism processes, this model is multi-scale in character due to the wide range of time scales characterizing its dynamics. Therefore, we utilize the Computational Singular Perturbation (CSP) algorithm, which has been used extensively in multi-scale systems of reactive flows and biological systems, to identify components of the system that (i) generate the characteristic time scale and the fast/slow dynamics, (ii) participate to the expressions that approximate the surfaces of equilibria that develop in phase space and (iii) control the evolution of the process within the established surfaces of equilibria. It is shown that a decisive factor on whether the ANLS or NALS configuration will develop during neuronal activation is whether the lactate transport between astrocytes and interstitium contributes to the fast dynamics or not. When it does, lactate is mainly generated in astrocytes and the ANLS hypothesis is realised, while when it doesn’t, lactate is mainly generated in neurons and the NALS hypothesis is realised. This scenario was tested in exercise conditions.
Keywords:
Glucose metabolism – Algorithms – Evolutionary rate – Neurons – Glucose – Reaction dynamics – Reactants – Astrocytes
Introduction
The human brain uses glucose as its main source of energy. Even though the adult brain accounts for ∼2% of the body weight, it consumes ∼20% of glucose-derived energy. Brain neurons have the highest energy demand [1] and the brain becomes dependent on glucose for proper functioning.
Glucose is metabolized to pyruvate in the presence oxygen, in both astrocytes and neurons. Pyruvate can follow two different pathways: oxidative or aerobic (producing energy in the mitochondria, in the presence of oxygen) and non-oxidative or anaerobic (producing energy and lactate in the cytosol, before the involvement of oxygen). Even though brain metabolism has been the subject of interest and research for decades, there is still a lot of discussion on to what extent energy derived from glucose is consumed by the astrocytes and the neurons themselves, and how much is shuttled between astrocytes and neurons in the form of lactate [2]. Two different hypotheses were introduced for the shuttling of lactate; one advocating the flow of lactate from astrocytes to neurons and another one advocating the flow of lactate from neurons to astrocytes. These two hypotheses are supported by experimental data, which are indirect “because current technology does not have adequate spatiotemporal resolution to quantify metabolic activity in single cells in vivo” [3].
The astrocyte-to-neuron lactate shuttle (ANLS) was proposed on the basis of glutamate-evoked increases in glucose utilization and lactate release by cultured astrocytes [4, 5]. According to ANLS, the neuronal activity triggers a glucose uptake in astrocytes, which leads to large increases in production of lactate. Lactate is then released to the extracellular space and transports to nearby neurons, where it is used as substrate for energy production (for detailed review see Ref. [5]). The validity of the ANLS hypothesis has been supported by experimental evidence over the years; many of them can be found in Refs. [6, 7].
On the other hand, the neuron-to-astrocyte lactate shuttle (NALS) states that lactate is generated by neurons and taken up by astrocytes during activation. Due to the high energy demands in neuronal activation, neurons take up glucose and transfer lactate to astrocytes, since the neuronal glucose transporter, GLUT3, appears to have higher transport rate compared with the astrocytic glucose transporter, GLUT1 [8–12]. Similarly to the ANLS hypothesis, there is experimental evidence that supports the NALS hypothesis [13–15].
In the context of the ANLS hypothesis, neurons are unable to increase their glycolytic activity in response to neuronal activation [16–18]. Therefore, enhanced glycolysis due to neuronal activation can take place only in astrocytes. On the other hand, in the context of the NALS hypothesis, neuronal activation enhances glycolysis in neurons [19–21].
A large portion of the scientific community is currently split between the validity of either the ANLS or NALS hypothesis, while another portion advocates that both hypotheses can be valid depending on the operating conditions. There is a large body of literature advocating the validity of ANLS [6, 22–32], as this is supported by indirect experimental findings [4, 18, 32–44]. The main criticism of the ANLS hypothesis stems from the NALS supporters [8, 9, 45–54], and the related indirect experimental evidence [10, 13–15, 55–64]. On the other hand, recent findings suggest that perhaps brain responds differently to different environmental challenges, depending on the availability of the fuel source. For example, it was shown that under inadequate supply of glucose, alternative substrates contribute to brain energetics [7, 11], such as lactate and/or ketones during exercise, hypoxia, hypoglycemia or intense brain activity [65, 66]. Therefore, under different brain environmental conditions either the ANLS or NALS hypothesis can be valid. The identification of the conditions under which ANLS or NALS manifests is one of the objectives of this manuscript.
The very beginning of the brain chemistry can be traced back to 1884 [67] but early evidence on metabolism in brain tissue can be found in Himwich’s early work [68] dated in 1951. Since then, the brain energy metabolism network has been investigated extensively mainly through in vitro and in vivo experiments and measurements on rats (e.g., Refs. [34, 69–76]) and humans (e.g., Refs. [55, 70, 77–79]).
On the other hand, the development of mathematical models that aim to simulate the metabolic interactions between metabolites into and among the various cerebral compartments has been flourishing very recently. Aubert et al. were among the first to develop mathematical brain metabolism models in order to investigate the connections between electrical activity, energy metabolism and hemodynamics [80, 81]. These models served as baselines for the development of more detailed ones with various investigation objectives like: the ANLS hypothesis [82], the initial rapid decrement in the extracellular lactate concentration [83], the blood oxygenation level-dependent (BOLD) signal in functional MRI (fMRI) [84], the inclusion of glycogen dynamics in astrocytes [85]. The model proposed by Simpson et al. [9] was developed with the purpose of addressing the role of transport in cerebral metabolism and the limitations of the previously developed models related to the concentrations and kinetic properties of the glucose transporters (GLUTs) and monocarboxylate transporters (MCTs) in the brain. With this work, Simpson et al. provided the setup for both the ANLS and the NALS hypotheses. Although, this model was criticized for its simplisity, its increased accuracy compared to in vivo rodent data lead to its utilization (with small parameters’ adjustments) for studying in vivo data from humans at high field functional magnetic resonance spectroscopy (fMRS) by Mangia et al. [10].
Combining the models proposed in [9, 10, 81, 82], Di Nuzzo et al. [8] introduced a new model to investigate the effects of cellular transport, metabolic capacity, specialization for energy use and metabolite trafficking during normal and enhanced cerebral activity. Aiming on the steady state conditions of the metabolic pathways and their relation to neuronal activation, Occhipinti et al. initially proposed a model [86] and used statistical tools to examine the reaction fluxes and the transport rates. Later this model was revised in [87] in order to investigate the preferred energy substrate (i.e., glucose or lactate) in neurons at high activity. They further investigated the reasons of production of γ-aminobutyric acid (GABA), by revising accordingly their latest model in [88]. By modifying and extending the models of Occhipinti et al. [87, 88] for glutamatergic neurons and astrocyte-GABAergic neuron cellular complex, respectively, Calvetti at al. [89] proposed a very detailed steady state model comprising of three cells, each one subdivided into cytosol and mitochondria, where astrocytes, glutamatergic, and GABAergic neurons interact through a common extracellular space (ECS). Finally, the models proposed in [82, 84] were used recently as baseline for the development of a more sophisticated model by Jolivet et al. [22], in order to investigate the ANLS hypothesis and its connection to cerebral blood flow.
Sensitivity analysis and flux analysis are the standard and widely used techniques for the investigation of such complex models. The former technique quantifies the effect of a perturbation in the model’s parameters to the evolution of the system. It is often utilized to obtain insights into the metabolic properties of the system by examining the derivative of fluxes with respect to each model parameters [8, 85], but also to quantify the relative contribution of input parameters on the features of an output attribute, such as BOLD signal [90, 91]. On the other hand, flux balance analysis targets on finding the physiological feasible configuration of reaction fluxes and transport rates capable on maintaining a steady state. It is often utilized to verify the hypotheses on the role of different pathways related to a brain function that have been proposed and to predict the activation levels of pathways at different inputs according to the initial hypotheses [86, 92]. Flux balance analysis can be applied in a mathematical model even when only the stoichiometry is known, and therefore it has been widely used and extended in both deterministic [93, 94] and stochastic frameworks [87, 95, 96]. Flux balance analysis is indeed useful in the interrogation of biochemical systems, but its value is limited since it can lead to conclusions that are only global in character. In addition, sensitivity coefficients (determined through conventional local sensitivity analysis) can be misleading, since a parameter can be related to insignificant sensitivity coefficient and yet be important [97]. Moreover, even in the early 80s works it was acknowledged that there are functional connections between the sensitivity coefficients that they cannot be identified by their mere inspection, thus, additional treatment of the local sensitivity coefficients was proposed [98]. As such, none of these methods can provide details on the temporal evolution of the set of biochemical reactions and metabolites that are most influential on the dynamics and progress of the system [94]. Moreover, they cannot identify the chemical reactions that control the evolution of a metabolite or a variable in general. Such identifications are of paramount importance since they can facilitate the associated physical understanding and eventually can lead to the development of efficient drugs.
Here, the Computational Singular Perturbation (CSP) algorithm will be employed for the analysis of a brain lactate metabolism model. CSP is an algorithmic method for asymptotic analysis developed in the late 1980s by Lam and Goussis [99–101]. It exploits the fact that systems describing the time evolution of biological kinetics models, are driven by processes characterized by a wide range of time scales. In such systems, the processes characterized by the fastest time scales become quickly equilibrated (exhausted). The system is thus confined to evolve within the related surfaces of equilibria that emerge in phase space and is driven by the processes that are characterized by the slow time scales. CSP can provide both mathematical expressions that approximate the generated surfaces of equilibria and the simplified system that approximates the slow evolution along these surfaces. CSP was initially developed to treat large and complex chemical kinetic mechanisms in the context of reacting flows [100, 102], but was later employed for the analysis of other physical problems, like computational mechanics [103], atmospheric science [104], etc., including biological systems [105–108] and pharmacokinetics [109, 110]. CSP exhibits many advantages compared to other similar approaches. Firstly, it is fully algorithmic, requiring no input from the investigator, apart from the detailed model and the accuracy that the simplified model is required to provide. In addition, the CSP algorithm can identify the physical processes in the model that contribute to the emergence of the surfaces of equilibria, those that control the slow (long-term) evolution of the system and those that are responsible for the development of the fast and slow time scales. Finally, CSP can identify the variables whose evolution is characterised by the fast or slow time scales. This understanding allows the investigator to identify the processes that can control not only the slow evolution of the whole system, but also the evolution of each individual metabolite; e.g., see the applications in pharmacokinetics [109, 110].
The CSP algorithm refers to the Geometric Singular Perturbation (GSP) theory, which introduced slow invariant manifolds (low dimensional surfaces on which trajectories evolve according to the slow time scales) and the fast fibers (along which trajectories approach the slow manifold according to the fast time scales) [111–113]. At each point on the manifold, CSP provides basis vectors that span the fast and slow subdomains of the tangent space. These vectors are provided by two iterative procedures; one that improves the accuracy of the basis vectors that are tangent to the manifold and one that improves the accuracy of the basis vectors aligned with the fast fibers. Considering the small parameter ϵ, equal to the ratio of the characteristic fast and slow time scales, each of these iterations increases the accuracy of the CSP basis vectors by O ( ϵ ) [114–117]. With this set of basis vectors being available, the slow manifold is approximated by the relations that result from the negligible projection of the vector field along the fast basis vectors and the system that governs the flow on the manifold is approximated by the projection of the vector field along the slow basis vectors [100]. Due to the multi-scale character and the complexity of the models that are currently of interest, approaches based on the GSP theory are becoming popular in the field of biology [118–122]. The major advantage of CSP is that is algorithmic and does not require the mathematical model to be cast in a non-dimensional form [105–110].
In the current study, CSP and its algorithmic tools are employed in order to analyse brain lactate metabolism, using the model proposed by Simpson et al. [9], as modified by Mangia et al. [10]. This model was selected for the following reasons. Firstly, it is a fairly simple model, thus, an ideal candidate for the demonstration purpose of the proposed method (i.e., CSP and its tools). Secondly, the model of Mangia et al. [10] provides a basis for both the ANLS and the NALS hypotheses. Therefore, both hypotheses can be examined, in view of the same mechanism, which can lead to a fair comparison of the results. Thirdly, it is among the very few models that was fitted against human in vivo data [10, 28, 84], in contrast to the majority of the models that are calibrated against in vitro human data [8, 80–84, 88, 89] or in vivo rat data [9, 28, 85].
In the following, results obtained from the CSP analysis of a simple brain lactate metabolism model will be reported. By investigating the fast and slow dynamics of the model, the couplings among the glucose and lactate pathways are revealed in the context of both the ANLS and the NALS hypotheses. In addition, the rate limiting steps in each of the two pathways are identified. These findings are crucial in order to understand the mechanism for lactate shuttling in the brain, and for the identification of the most conducive point in the pathways for the control of the process. These findings will be demonstrated in normal and exercise conditions of neuronal activation.
Methods
The computational model
The brain lactate metabolic network is analyzed by utilizing the model introduced in Simpson et al. [9], as modified by Mangia et al. [10], which simulates the kinetic behavior of glucose and lactate among the cerebral compartments. As shown in Fig 1 the model consists of 5 compartments: the endothelium (e), the basal lamina (bl), the astrocyte (a), the interstitium (int) and the neuron (n). Glucose (Glc) and lactate (Lac), which are the only chemical species that are accounted for in the model, are oxidized in astrocytes and neurons in order to produce energy. The transport or diffusion of these two substances from one compartment to another is denoted as a reaction. As shown in Fig 1, there are 12 reversible and 6 irreversible reactions in the model; see S1 Text for a detailed description of the model. At this point it suffices to state that the interactions of Glc and Lac through the various compartments is governed by a 10-dimensional system of Ordinary Differential Equations (ODEs):
Computational singular perturbation (CSP) method and algorithmic tools
The Computational Singular Perturbation (CSP) method relies on the multi-scale character of the system under investigation and leads to an algorithm that delivers everything the traditional singular perturbation technique does. As the system evolves, CSP identifies the subspaces in phase space, in which the fast and slow dynamics act. The components of the model that tend to drive the system along the fast directions equilibrate, generating thus constraints in which the system is bound to evolve. This evolution is driven by components that act along the slow directions. With the fast and slow subspaces provided by CSP at each point in time, the components of the model that are mainly responsible for the generated constraints and for driving the system are easily identified.
According to the CSP approach, the vector field g(y) is resolved in N modes, so that Eq (1) is cast in the form [100, 125]:
The M constraints in Eq (3) fr = (βr ⋅ S1)R1 + … + (βr ⋅ SK)RK ≈ 0 (r = 1, …, M) are the result of significant cancellations among some of the additive terms (βr ⋅ Sk)Rk (k = 1, …, K). The reactions that contribute significantly to the formation of each of the M constraints are identified by the Amplitude Participation Index (API):
The formation of the M constraints and the dynamics of the slow system in Eq (3) are characterized by the M fastest time scales and by the fastest of the N − M slow ones, respectively. These time scales are approximated by the inverse of the eigenvalues of the Jacobian J, τn = |λn|−1 (n = 1, …, N). The reactions that contribute significantly to the generation of these time scales are identified by the Time scale Participation Index (TPI):
Each chemical species associates differently to each exhausted CSP mode, e.g., a chemical species can relate mostly to the m-th CSP mode (m = 1, …, M) and much less to the rest. The relation of the m-th CSP mode to the various chemical species is assessed by the CSP Pointer (Po):
The reactions that contribute the most to the evolution of the system within the M constraints, according to the N-dim. system of ODEs in Eq (3), are identified by the slow Importance Index (II):
Results
The dynamics of the system
The evolution of glucose and lactate concentrations in the various cerebral compartments is simulated during neuronal activation by the model described briefly in the previous section and in detail in S1 Text. The variations of the concentrations of metabolites during neuronal activation (electrical and chemical excitability of neurons) have been recorded experimentally by fMRS [10]. It was shown that neuronal activation begins (ends) with a rapid increase (decrease) in lactate and a decrease (increase) in glucose in the human cortex. In order for the model employed to reproduce this response of the concentrations, the parameters in the model were properly adjusted in the periods of neural activation. In particular, the system was allowed to rest on its baseline steady-state conditions for the first 1000 s and then the parameters associated with lactate (glucose) utilization were adjusted in order to generate increased (decreased) levels of lactate (glucose) during the first period of neuronal activation P1 (1020 s < t < 1240 s). These adjustments were removed at the end of P1 and the system was let to evolve towards the next period, P2 (1320 s < t < 1540 s). As in P1, adjustments in the parameters were introduced throughout P2 in order to produce decreased (increased) lactate (glucose) levels. After the completion of P2, the adjustments were removed and the system was allowed to reach steady state conditions. For a detailed explanation of the parameter adjustments during P1 and P2 see Ref. [10].
Multi-scale character and modal decomposition
The multi-scale character of the model is established by wide gaps among the fast and slow time scales τn and by the ensuing negligible amplitude of the fast modes fn.
There are N = 10 variables in the model employed, so that according to Eq (2) there are 10 amplitudes, which are related to 10 time scales. Fig 2 displays the evolution of the N time scales τn that characterise the dynamics of the metabolic network under consideration and the amplitudes fn of the N CSP modes, for both ANLS (top) and NALS (bottom) cases. The period considered includes the periods P1 and P2 of neuronal activation, as indicated in the figure. The rapid “step” deviations either in time scale or in amplitude evolution relate to the parameter adjustments during P1 and P2. These adjustments have a negligible influence in the time scales, but a significant one in the amplitudes.
In particular, in the left part of Fig 2 it is shown that in both the ANLS and NALS cases the time scale spectrum extends from O(1) to O(103) s. These time scales are all dissipative (the related eigenvalues are all real and negative), so they are generated by processes that tend to drive the system towards a fixed point. There is a general agreement in the evolution of the time scales in the ANLS and NALS cases. The largest deviations are recorded in the 5th and 8th time scales: τ5 is somewhat slower in the NALS case, while τ8 is much slower. The most notable difference in the two cases is that the largest time scale gap is encountered between the 8th and the 9th scales in the ANLS case and between the 7th and the 8th scales in the NALS case. This finding suggests that there is one more constraint established in the ANLS case.
The corresponding amplitudes fi of the CSP modes are displayed in the right part of Fig 2, where the 5 first modes are denoted with dotted lines and the last 5 with solid ones. Before the neuronal activation (t < 1020 s), where the system has already attained its steady state conditions, the two largest amplitudes are O(10−23) and O(10−22) and relate to the slowest 9th and 10th CSP modes respectively; indicating that all CSP modes have been exhausted (i.e., all 10 modes are inactive). As the neuronal activation manifests itself at t = 1020 s, all amplitudes increase. The largest amplitudes during the neuronal activation periods P1 and P2 in both ANLS and NALS cases are those of the two slowest 9th and 10th CSP modes, reaching O(10−16), followed by the amplitude of the 8th mode. In the ANLS case, the amplitude of the 8th CSP mode relaxes very quickly after its rapid increase at the start of neuronal activation. In contrast, in the NALS case the amplitude of the 8th CSP mode attains a relatively large value throughout the period of neuronal activation. This is in agreement to the finding that the largest fast/slow time scale gap is encountered between the 8th and the 9th time scales in the ANLS case and between the 7th and the 8th time scales in the ANLS case; i.e., Fig 2 shows that during neuronal activation only the amplitudes of the slow modes are large. After the end of P2, all amplitudes decay according to the related time scale as the system returns to its steady state; the slowest decay exhibited by the amplitude of the slowest 10th mode.
Our study will focus on the periods P1 and P2 of neuronal activation, which are displayed with dashed red arrows in both the left and right panels of Fig 2. Period P1 relates mostly to the glycolytic uptake in astrocytes, while period P2 relates to the Lac oxidation in both neurons and astrocytes [10]. Both periods display quite similar dynamical behavior. In particular, in both the ANLS and NALS cases the amplitudes of all modes grow at the start of the two activation periods, but the majority of them relax very quickly. In the ANLS case the number of exhausted time scales is M = 8 for the largest part of P1 and P2. Therefore, the 9th and 10th CSP modes are considered to be the slow ones (i.e., the ones that drive the system in P1 and P2), while the 1st to the 8th are the exhausted fast ones (i.e., their negligible amplitudes define the confines in which the system evolves). In contrast, in the NALS case, the number of exhausted time scales is M = 7 for the same periods. Now, the 8th, 9th and 10th CSP modes are considered slow, while the 1st to the 7th modes are fast.
In the following, the CSP diagnostic tools discussed previously are utilized in order to acquire physical understanding and derive conclusions about the evolution of the system during neuronal activation. Interested in the long-term evolution of the system, the analysis will focus on the slow dynamics. However, in order to understand the influence of the surface of equilibria, which is generated by the M constraints, on the behavior of the system, a brief presentation of the established constraints is presented first. These constraints involve the most intense components of the system, the equilibration of which allows the less intense ones to drive the system.
Fast dynamics and established constraints
The fast dynamics of the system is responsible for the constraints that are established when the fast CSP modes become exhausted (fr ≈ 0, r = 1, M). As previously stated, during neuronal activation, the number of exhausted modes is identified to be M = 8 and M = 7 in the ANLS and NALS cases, respectively. These fastest CSP modes were analysed with the CSP tools and the resulting diagnostics are displayed in Table 1 for the ANLS case and in Table 2 for the NALS case (the diagnostics of the non-exhausted 8th mode in the NALS case are included in Table 2 for comparison with the exhausted 8th mode in the ANLS case). The results displayed in the two Tables were computed at t = 1100 s, which belongs to P1. These diagnostics are indicative of those throughout P1 and P2; only 1-2% differences in the CSP diagnostics are recorded. Along with the Pointer Po, the Amplitude Participation Index API and the Time scale Participation Index TPI, the product bm⋅Sk is displayed in Tables 1 and 2. When b m · S k = O ( 1 ), the stoichiometric vector of the k-th reaction has a significant component along am. Since am is a fast CSP basis vector, b m · S k = O ( 1 ) indicates that the k-th reaction is a fast one and can potentially contribute to the fast dynamics.
As stated previously, the m-th exhausted CSP mode is distinguished by (i) the reactions that participate in the constraint expressed by the relation fm ≈ 0 (identified by the API index, Eq (4)), (ii) the reactions that generate the time scale τm that characterizes the formation of this constraint (identified by the TPI index, Eq (5)) and (iii) the variables (concentrations) that relate the most to this mode (identified by the Po index, Eq (6)). These CSP diagnostics for the exhausted modes are displayed in Tables 1 and 2.
The main observation from Tables 1 and 2 is that among the eight fastest modes 4 relate to the Glucose path and 4 relate to the Lactate path; i.e., the 1st, 4th, 5th and 8th modes relate to variables and reactions along the Glucose path (Glc-path modes) and the 2nd, 3rd, 6th and 7th ones relate to variables along the Lactate path (Lac-path modes). In particular,
-
The pointed by Po variables of the Glc-path modes (Lac-path modes), relate only to Glc (Lac).
-
The reactants of the reactions identified by TPI to contribute to the development of the four time scales τ1, τ4, τ5, τ8 that relate to the Glc-path modes (τ2, τ3, τ6, τ7 that relate to the Lac-path modes), involve only Glc (Lac).
-
The reactants of the reactions identified by API to contribute to the constraints fm ≈ 0, associated to the four Glc-path modes, involve only Glc. To a large extent, a similar situation arises in the four constraints that relate to the Lac-path modes; i.e., the reactants of almost all contributing reactions involve Lac. The only exception in this case is that in the constraint related to the 6th Lac-path mode the Glc-consuming reaction 7 provides a small contribution in the ANLS case, while the Glc-consuming reaction 9 provides a small contribution in the NALS case.
Therefore, considering the couplings of the two paths in the expressions that approximate the manifold M:
-
the Lactate path couples to the Glucose path through the 6th mode; this coupling takes place inside the astrocyte compartment in the ANLS case (reaction 7) and in the neuron compartment in the NALS case (reaction 9) and
-
the Glucose path is decoupled from the Lactate path in both the ANLS and NALS cases.
The conclusion is that a perturbation in the reactions on the Lactate path will not alter the Glucose-related constraints. On the other hand, a perturbation in the reactions on the Glucose path will alter the Lactate constraints; as stated previously, this influence is exercised inside the astrocyte compartment in the ANLS case and in the neuron compartment in the NALS case.
Let us now examine the four constraints that develop along the Glucose path, on the basis of the CSP diagnostics displayed in Tables 1 and 2; i.e., those related to the 1st, 4th, 5th and 8th modes in the ANLS case and those related to the 1st, 4th and 5th modes in the NALS case. The 1st (fastest) mode denotes in both ANLS and NALS cases an equilibration between the 2nd and the 4th reactions:
A schematic representation of these constraints that develop along the Glucose path in the ANLS and NALS cases is shown in Fig 3. Clearly, the Glc-related constraints involve all five compartments along the Glucose path and all means of transport and diffusion among them, with the contribution of the 3rd reaction being negligible, since its rate is sufficiently small; see S1 Text.
The four constraints that develop along the Lactate path relate to the 2nd, 3rd, 6th and 7th modes, for both the ANLS and NALS cases. On the basis of the CSP diagnostics displayed in Tables 1 and 2, the following conclusions can be reached. The 2nd mode denotes mainly the equilibration between the forward and backward directions of the 14th reaction, corrected by the net contribution of the 16th and 15th reactions:
A schematic representation of these constraints that develop along the Lactate path is shown in Fig 3. As with the Glc-related constraints, the four Lac-related constraints involve all five compartments along the Lactate path and all means of transport and diffusion among them, with the contribution of the 13th reaction being negligible, since its influence is overshadowed by that of the 14th and 16th reactions; see S1 Text.
The system evolves on a surface of equilibria, approximated by these constraints, on which the slow dynamics become the characteristic ones. Having established, from the data displayed in Tables 1 and 2, that the species related the most to the slow dynamics are Glca, Glcn, Laca and Lacn (since they relate the least to the fast dynamics), we proceed forward by analyzing the metabolic profiles of neurons and astrocytes.
Slow dynamics
Interested in examining the dynamics of neurons and astrocytes interaction, the slow evolution of the system during neuronal activation is analyzed with the CSP tools. The capability to algorithmically determine the underlying physics of the slowly evolving system, within the constraints imposed by the exhausted modes, is demonstrated in this section.
Characteristic dynamics and driving reactions during the slow evolution
As stated in the previous section, during neuronal activation (periods P1 and P2), only the 9th and 10th modes are active in the ANLS case and only the 8th, 9th and 10th modes are active in the NALS case. Given that τ9 < τ10 in the ANLS case, the slow evolution of the system is characterized by τ9. However, since τ8 ≈ τ9 < τ10 in the NALS case, the slow evolution of the system is characterized by both τ8 and τ9. Tables 3 and 4 display CSP diagnostics (Po, TPI) of the slow modes in the ANLS and NALS cases, respectively, computed at t = 1100 s.
It is shown in Table 3 that in the ANLS case the pointed variables and the reactions contributing to the related τi of the two slow modes refer to either the Lactate path (9th mode) or to the Glucose path (10th mode). As it will be shown next, this feature is due to the fact that these two modes are practically decoupled. Naturaly, the pointed variables in the two modes are the ones related the least with the fast dynamics, as it is evident by the results displayed in Table 1.
Regarding the 9th mode, the Lacn-consuming reaction 18 and the Laca-consuming reaction 17, which represent the Lac oxidation in neurons and astrocytes, respectively, account for more than 50% of the total contribution to the generation of τ9. On the basis of their negative TPI, it is concluded that these reactions tend to drive the system to a fixed point. On the other hand, the contributions to τ9 of the two directions of the 16th reaction, which represent the transport of Lac from interstitium to neuron and vice versa, are of opposing nature; the Lacint-consuming forward direction promoting the movement towards a fixed point and the backward opposing it. Their contributions cancel each other, so that the net influence of the 16th reaction is negligible. The dominance of the 17th and 18th reactions in the generation of τ9 is reflected in the large Po value of their reactants Laca and Lacint. This dominance implies that the duration of the neuronal activation is mainly controlled by these two reactions. For example, a decrease of the Lac oxidation rate via reaction 18 will lead to the increase of τ9; i.e., it will lead to a slower evolution of the process. This influence of the 18th reaction will be demonstrated later by properly perturbing R18.
The results displayed in Table 3 also show that in the context of the 10th mode in the ANLS case, significant contributions to the generation of τ10 are provided by the Glcint-consuming reaction 4b, the Glcn-consuming reaction 6b and the Glca-consuming reaction 5f. Other reactions provide additional contributions, mainly via the already established constraints. Counting the net contribution of each reaction (by adding the contributions of the forward and backward direction), the most significant influence originates from the 1st reaction, followed by the 2nd and then by the 4th.
Table 4 shows that, as in the ANLS case, the pointed variables and the reactions contributing to the related τi of the three slow modes in the NALS case refer to either the Lactate path (9th mode) or to the Glucose path (8th and 10th mode). The dynamics of the two slowest modes (9th and 10th) are similar to the equivalent modes in the ANLS case with small variations. However, in contrast to the ANLS case, in the NALS case the 8th mode is accounted in the slow component of the vector field and its dynamics characterise the evolution of the system, along with the dynamics of the 9th mode (since τ8 ≈ τ9). It is shown in Table 4 that the largest net contribution in generating τ8 originates from the Glca-consuming 5th reaction, while the net contributions of the 6th and 4th reactions are much smaller. In comparison to the ANLS case, the dynamics of the slowest 10th mode in the NALS case mainly differs in that Glca is more pointed and Glcn is less, while the differences in the 9th mode are minor.
Comparing the slow dynamics of the ANLS and NALS cases, Tables 3 and 4 show that the slow dynamics along the Glc-path shifts from neurons towards astrocytes, as it is evident by the sum of the CSP Pointer of Glcn and Glca due to the slow modes (D G l c n 9 + 10 = 0.526 and D G l c a 9 + 10 = 0 . 280 in the ANLS case and D G l c n 8 + 9 + 10 = 0.970 and D G l c a 8 + 9 + 10 = 0 . 772 in the NALS case; i.e., a 84.5% and a 175.7% increase of the sum of slow Pointers of Glcn and Glca, respectively). This shift of the slow dynamics along the Glc-path is accompanied by an opposite shift of the fast dynamics; i.e. from astrocytes towards neurons. This feature is clearly evident in Fig 3, by the significantly reduced contribution of the Glca-related 5th reaction in the established constraints in the NALS case and the ensuing reduced influence of Glca in the fast dynamics. In turn, as shown in Tables 1 and 2, this last feature results in a reduced involvement in the fast dynamics of the Glca-consuming 7th reaction, that couples the Glc-path to the Lac-path.
The decoupling of the Po and TPI diagnostics along the Glc- and Lac-paths in the context of the slow modes in the ANLS and NALS cases, as indicated by Tables 3 and 4, is reasonable for the Glc-path, since the concentration of Glc in the various compartments is not dependent on that of Lac, as shown in Fig 1. However, this is not that obvious for the Lac-path, since it is influenced by the Glc-path via the 7th and 9th reactions. In order to investigate this issue further, the reactions contributing the most to the evolution of Glca, Glcn, Laca and Lacn (thereafter, Glca,n and Laca,n) will be identified, with the use of the Importance Index II. The results displayed in Table 5 identify (i) the reactions that influence most the evolution of Glca,n and Laca,n and (ii) the corresponding values of II. The positive (negative) sign of the contribution of a reaction displayed in Table 5, denotes its effect in producing (consuming) the related species.
Table 5 shows that in the ANLS case the evolution of Glca,n is mainly determined by the R1f (Glc transport from serum to endothelium), R1b (Glc transport from endothelium to serum) and R2f (Glc transport from endothelium to basal lamina). Smaller influences of the rates R7 and R2b towards depleting Glca,n are also reported, while the even smaller influences of the two rates R4f and R4b cancel each other. Regarding Laca,n, Table 5 shows that its evolution is mainly determined by the rates R7 (glycolytic rate in astrocytes) and R18 and R17 (Lac oxidative rate in neurons and astrocytes, respectively); the first tending to increase Laca,n, while the latter two to decrease it. Smaller influences are reported by R6f and R6b, which cancel each other. Table 5 shows that R9 has a small influence in driving both Glca,n and Laca,n.
Regarding the reactions that drive Glca,n and Laca,n in the NALS case, the differences in the results displayed in Table 5 relative to the ANLS case are mainly due to the decreased influence of fast dynamics on Glca. As it is shown in Table 5, the evolution of Glca is now mainly determined by reactions that either consume it or produce it (R5f, R5b, R7), while the recorded net contributions by the 1st, 2nd and 6th reactions are smaller. The reactions that mainly determine the evolution of Glcn in the NALS case are similar to those identified in the ANLS case. The main difference is the decreased influence of the 7th reaction and the increased one of the 9th. This is also the main difference in the reactions that drive Laca,n; apart from R18 and R17 that play a major role in the ANLS case, now R9 exhibits the strongest influence, while that of R7 is significantly diminished.
From the ANLS to the NALS configuration
The differences between the ANLS and NALS configurations lies in the different values of the parameters of reactions 3, 5, 7, 8, 9, 10, 17 and 18. In order to investigate the influence of these different values, each of the changes in the value of the parameters, leading from the ANLS to the NALS configuration, was implemented separately. It was shown that only the changes in the parameters of the reactions 5, 7, 8, 9 and 10 produce significant deviations from the solution obtained with the ANLS model. The influence of these changes is demonstrated in Fig 4, where the evolution of the net rates R15 = R15,f − R15,b and R16 = R16,f − R16,b is displayed. These two net rates relate to the Lac flow from astrocytes to interstitium and from there to neurons, in the ANLS (A) and NALS (N) cases. As shown in the figure, in the ANLS case R15 and R16 are positive, while in the NALS case are negative. In addition, Fig 4 displays the evolution of these two rates in the case where each of the changes in the parameters of the reactions that lead from the ANLS to the NALS configuration is implemented separately; e.g., “R15 A + 7, 8” indicates the profile of R15 when the ANLS configuration is considered with the exception of the parameter in the rate of reactions 7 and 8 for which the change that leads to the NALS configuration is accounted for (see S1 Text for a list of all changes in the parameters that lead from the ANLS to the NALS configuration).
It is demonstrated in Fig 4 that the most influential change in the shift from the ANLS to the NALS configuration is produced by the parameter that modifies R5, followed by the one modifying R9 and R10. In contrast, the change in the parameter that modifies R7 and R8 shifts the ANLS configuration away from the NALS one. This finding is in agreement with the substantially different influence of R5 in the fast dynamics of the ANLS and NALS cases, discussed previously; i.e., in the ANLS case R5 is fast enough that allows the coupling of the Glc and Lac pathways via the 7th reaction (so that the flow of Lac is directed from astrocytes to neurons), while in the NALS case R5 is slow enough, so that the coupling between the two pathways is mainly established via the 9th reaction (so that the flow of Lac is directed from neurons to astrocytes).
Assessment of CSP diagnostics
The conclusions reached on the basis of the various CSP diagnostic tools will be assessed here, by perturbing the rates of selected reaction rates. For simplicity, only the ANLS case will be considered; similar conclusions are reached in the NALS case. In particular, consider the rate R18, which does not contribute to the generation of the eight constraints, as shown in Table 1. In addition, the TPI results in Table 3 suggested that the rate R18 has a significant contribution to the generation of the characteristic time scale τ9. Finally, the II results in Table 5 suggested that R18 tends to decrease Laca,n. Therefore, a decrease of the rate constant of R18 is expected to increase τ9 and Laca,n. Indeed, as it is shown in the left part of Fig 5, a 20% decrease of R18 causes Lacn, Laca and Lacint to reach a constant value at t ≈ 1805 s, while in the unperturbed case this state is achieved faster, at t ≈ 1235 s. In addition, the levels of Laca,n are higher in the perturbed case, as expected. On the other hand, the increase of the rate R10 by 20% has negligible influence on the evolution of Laca,n, as shown at the right part of Fig 5. This is an expected result, since R10 has a negligible contribution to the eight constraints, to τ9 and to the slow evolution of the system; see Tables 1, 3 and 5.
A further demonstration of the usefulness of the CSP diagnostic tools emerges from the physical interpretation of the II results displayed in Table 5. In particular, consider the response of Laca,n and Glca,n when subjected to perturbations of R18, R9, R1f and R3f during neuronal activation. Such perturbations might affect the M constraints or the reduced model; see Eq (3). According to the API values in Table 1, R1f is the only rate from this set that influences the constraints, via the 4th and the 5th modes that relate to Glce and Glcint. According to the II values in Table 5 that assesses the influence of the reactions in the sow model, (i) R18 has the largest contribution towards decreasing Laca,n levels, while has negligible influence in the Glca,n, (ii) R9 exhibits a small influence in increasing Laca,n levels and a smaller influence in decreasing Glca,n levels, (iii) R1f is the major contributor towards increasing Glca,n, but has negligible influence in the evolution of Laca,n and (iv) R3f has a negligible contribution to the evolution of all variables.
The predicted influence of perturbations of R18, R9, R1f and R3f is demonstrated in Figs 6 and 7. Specifically, in the left part of Fig 6 it is shown that the decrease of R18 leads to increased Laca,n levels, but leaves Glca,n unaffected. On the right part of the same figure, it is shown that the increase of R9 leads to increased Laca,n levels and decreased Glca,n levels. Note that the response of Laca,n to a perturbation of R18 is more intense than that of Glca,n and Laca,n to a similar perturbation of R9, in accordance to the related II results displayed in Table 5. Also, in the left part of Fig 7 it is shown that the decrease of R1f leads to increased Glca,n levels, but leaves Laca,n levels unaffected. Finally, in the right part of the same figure it is shown that an increase of R3f has a negligible influence to both Glca,n and Laca,n, as predicted.
Investigating an exercise scenario
It was shown previously that a decisive factor on whether the ANLS or NALS configuration will develop is whether the 5th reaction contributes to the fast dynamics or not. In particular, it was shown that in the context of the ANLS (NALS) configuration the 5th reaction contributes (does not contribute) to the fast dynamics, allowing thus the major coupling between the Glc and the Lac paths to be established in astrocytes via the 7th reaction (in neurons via the 9th reaction). As a result, the Lac flow is directed from astrocytes to neurons in the ANLS case and in the opposite direction in the NALS case. These conclusions were reached for normal conditions (NC); among others, low Lac concentration in serum (Lacs = 1 mM).
However, under exercise conditions (EC) the Lac concentration in the serum reaches large values, about 15 − 25 mM (“all-out” maximal exertion) [123, 124]. A high Lac concentration in serum will certainly influence the concentrations along the Lac-path. Nonetheless, this factor will leave the concentrations of Glc unaffected, since the Glc-path is not coupled to the Lac-path. As a result, it is not expected to influence the degree to which the 5th reaction contributes to the fast dynamics; i.e., the conditions under which the Glc-path dynamics favor the development of ANLS or NALS. The only feature that might be affected is the flow direction of Lac between astrocytes and neurons, given the changes in the concentrations along the Lac-path that result by the high Lac concentration in the serum. These issues will be examined next.
In the exercise scenario considered here, the serum Lac concentration was increased to Lac = 20 mM (from Lac = 1 mM in NC). As expected, high serum Lac conditions, introduce significant changes to all Laci, but leave unaffected those of all Glci in both ANLS and NALS cases. This is displayed in the left panels of Fig 8, where the ratio of the concentrations computed under exercise conditions (EC) and under normal conditions (NC), yi,EC/yi,NC, is displayed. This ratio equals unity for Glci in all the 5 compartments, while it is slightly larger than unity for the Laci in all compartments, with the exception of the ratio for Lace which is around 6.
The increased concentrations of Laci modify the flow of lactate between astrocytes and neurons. As shown in the right panels of Fig 8, the high serum Lac concentration causes the rate of Lac-transport from astrocytes to interstitium (R15 = R15,f − R15,b > 0) to decrease and from interstitium to neurons (R16 = R16,f − R16,b > 0) to increase in the ANLS case. In contrast, in the NALS case the Lac-transport from neurons to interstitium (-R16 > 0) decreases and from interstitium to astrocytes (-R15 > 0) increases. Fig 8 shows that in both the ANLS and NALS cases the flow of lactate from basal lamina to intertitium (R14) increases. While in the ANLS case Lac flows from basal lamina to intertitium (R14 > 0) in both NC and EC, in the case of NALS the minor flow from intertitium to basal lamina (R14 < 0) in NC is significantly reversed in EC. A schematic presentation of the changes in the net rates R14, R15 and R16 is displayed in Fig 9. Clearly, the increased values of R16 in the ANLS case and of −R15 in the NALS case are in agreement with the nature of ANLS and NALS assumptions, respectively; since Lac is transported towards the neurons in the ANLS case and towards the astrocytes in the NALS case increases. However, the decreased values of R15 in the ANLS case and of −R16 in the NALS case are in disagreement; since the transport of Lac from astrocytes in the ANLS case and that from neurons in the NALS case decreases. In fact, as shown by the profile of R16 in the NALS case, at about t = 1300 s, Lac is transported to neurons (R16 > 0) and not from them, as the NALS hypothesis dictates.
In summary, it was shown that an increased serum Lac concentration resulted in an increased transport of Lac to neurons in the ANLS case, mainly due to the increased diffusion from basal lamina, accompanied by a decreased transport out of astrocytes. In contrast, an increased transport of Lac to astrocytes was recorded in the ANLS case, also due to the increased diffusion from basal lamina, accompanied by a decreased transport out of neurons. It can be shown that these features are in full agreement with the equilibrations that develop in the context of the 2nd and 3rd modes, as they are stated in Tables 1 and 2, which mainly involve the rates R16, R15 and R14 and to a lesser degree R13 in normal and exercise conditions and the rate R12 only in exercise conditions.
Discussion
Fast and slow dynamics influence the response of the system in both the ANLS and NALS frameworks in two distinct manners. The fast dynamics are responsible for the constraints that are established in phase space, within which the system evolves driven by processes responsible for the slow dynamics.
Constraints develop along both the Glucose and Lactate paths. Only one of these constraints couples the two paths; i.e., the Lactate to the Glucose path, while the latter is uncoupled to the former, see Fig 1. Along the Glucose path in both the ANLS and NALS cases the following constraints related to Glc-transport develop:
-
from endothelium to basal lamina (R2f − R2b) and from there to interstitium (R4f − R4b),
-
from serum to endothelium (R1f − R1b) and from there to basal lamina (R2f − R2b), with the contribution of a higher order correction from the Glc-transport in and out of the interstitium (R4f − R4b, R5f − R5b, R6f − R6b); the latter contribution based on the previous establishment of the constraint in (i),
-
from interstitium to neurons (R6f − R6b), from serum to endothelium (R1f − R1b) and only in the NALS case from astrocyte to interstitium (R5f − R5b), with the contribution of a higher order correction from the Glc-transport from basal lamina to interstitium (R4f − R4b) due to the constraints established previously in (i) and (ii), and
-
only in the ANLS case, from astrocytes to interstitium (R5f − R5b) and from there to neurons (R6f − R6b), with the glycolytic rate in astrocytes R7 contributing as a higher order correction.
As Fig 3 shows, these constraints involve all five compartments along the Glucose path. Along the Lactate path in both the ANLS and NALS cases the following constraints related to Lac-transport develop:
-
from basal lamina to interstitium (R14f − R14b), with the Lac-transport from interstitium via reactions 15 (R15f − R15b) and 16 (R16f − R16b) providing a higher order correction,
-
from astrocytes to interstitium (R15f − R15b) and from there to neurons (R16f − R16b),
-
from astrocytes to interstitium (R15f − R15b) and from there neurons (R16f − R16b), with a small contribution from reactions 7 and 9 from the Glucose path in the ANLS and NALS cases, respectively, and
-
from serum to endothelium (R11f − R11b) and from there to basal lamina (R12f − R12b), with only a small contribution from reaction 16 which is possible after the equilibria in (i), (ii) and (iii) are established.
As with the constraints related to the Glucose path, Fig 3 shows that these constraints involve all five compartments along the Lactate path.
Clearly, the major difference in the ANLS and NALS cases is the contribution to the fast dynamics of Glc-transport from astrocytes to interstitium (R5f − R5b). In the ANLS case this Glc-transport participates in both the fast dynamics and in the ensuing related constraints (large TPI and API in the 4th, 5th and 8th modes; see Table 1). However, in the NALS case this participation is greatly diminished (small TPI and API in the 4th and 5th modes, while the 8th mode is now a slow one; see Table 2). These features suggest that the strong coupling ot the Lac- to the Glc-path that is established in the ANLS case (mainly via the 7th reaction inside the asrtocytes, in the context of the 6th and 8th modes) diminishes significantly in the NALS case (where this coupling is now established mainly via the 9th reaction inside the neurons, in the context of only the 6th mode).
The finding that—due to the behavior of R5—the Lac- to the Glc-path coupling is mainly established inside the astrocytes in the ANLS case and inside the neurons in the NALS case determines the conditions promoting the ANLS or NALS development. In the ANLS case R5 participates in the equilibria of the 4th and 5th modes, which allow the Glc concentrations in both the astrocytes and neurons to be coupled directly with that in the serum. The formation of Lac from Glc is more intense in astrocytes (R7) than in neurons (R9), so that the typical for the ANLS case flow of Lac from astrocytes to neurons is generated. In the NALS case the participation of R5 to the equilibria of the 4th and 5th modes is significantly reduced, so that only the Glc concentration in neurons is coupled directly with that in the serum. Therefore, the formation of Lac from Glc in neurons (R9) dominates, so that the typical for the NALS case flow of Lac from neurons to astrocytes is generated.
The influence of the fast dynamics in determining the development of ANLS or NALS is fully reflected in the slow dynamics that characterize the evolution of the system within the established constraints. In particular, as shown in Table 5, in the ANLS case the Glca,n level is mainly determined by the rate of reactions in the first half of the Glc-path (R1, R2) and the glycolytic rates in astrocytes (R7), while the Laca,n level is mainly determined by rates in the second half of the Lac-path (R17 and R18) and by the rate R7. In essence, the glycolytic rate in astrocytes (R7) couples the Lac oxidative rates in astrocytes and neurons (R17 and R18) with the Lac inflow rate from the serum, via reactions at the start of the Glc-path (R1, R2). As a result, the flow along the Lac pathway in the ANLS case is directed from astrocytes to interstitium (R15 = R15,f − R15,b > 0) and from there to neurons (R16 = R16,f − R16,b > 0). In contrast, Table 5 shows that in the NALS case the same coupling is now established mainly via the glycolytic rate in neurons (R9), while the glycolytic rate in astrocytes (R7) exhibits a smaller influence. In this case the flow along the Lac pathway is directed from neurons to interstitium (R16 = R16,f − R16,b < 0) and from there to astrocytes (R15 = R15,f − R15,b < 0). In both cases, the direct coupling of the Lac oxidative rates in astrocytes and neurons (R17 and R18) with the Lac inflow rate from the serum via the reactions in the first of the Lac-path (R12 and R13) has no significant influence. Finally, in agreement to the model depicted schematically in Fig 1, Table 5 shows that in both the ANLS and NALS cases Glca,n is independent from Laca,n.
It was shown that the shift in the Glc to Lac activity along the Glc-path from the astrocytes to neurons, in the NALS case relative to the ANLS case is reflected in both the fast and slow dynamics. This shift is accompanied by a weakening of the Glca-consuming/producing reactions 3f, 3b, 5f, 5b, 7 and 8 and a strengthening of the Glcn-consuming reactions 9 and 10, as shown in Fig 10. In the same figure it is shown that in the NALS case the glycolytic rate in astrocytes (R7) decreases and the glycolytic rate in neurons (R9) increases, relative to the ANLS case, supporting the findings reported previously on the influence of these reactions in coupling the oxidative rates in astrocytes and neurons (R17 and R18) with the Lac inflow rate in the serum.
The large Lac concentrations in the serum considered in the context of the exercise conditions, were shown to affect the flow of Lac, but left unaffected the conditions under which the Glc concentrations in astrocytes and neurons are coupled directly with that in the serum, since the Glc-path is decoupled from the Lac-path in the model considered. In particular, it was shown that under exercise conditions the increased rate of Lac diffusion from basal lamina to interstitium modified the rate of Lac transport between astrocytes and neurons; i.e., in the ANLS case contributed to the increased transport to neurons and the decreased transport from the astrocytes and in the NALS case contributed to the increased transport to astrocytes and the decreased transport from neurons. It is noted here that although there were significant changes in the initial conditions and the values of some parameters in the cases of normal and exercise conditions, the basic dynamical features were not altered. Specifically, there were not changes recorded in the reactions contributing to the time scales and the emerging constraints, the variables related to the modes and the number of exhausted modes. Only moderate changes were noticed on the degree to which some reactions contributed to the constraints and to the evolution of the system within these constraints.
The findings reported here are supported by experimental results in the literature. The high glycolytic rate of astrocytes (R7) in the ANLS case has been detected experimentally [4, 16, 69, 135, 136], since neurons cannot afford to sustain high glycolytic rate (R9) [6, 16]. In addition, the dominance of the glycolytic uptake in neurons (R9) has also been detected in the NALS case during neuronal activation [13–15]. Here, the dominance of R7 over R9 in the ANLS case and of R9 over R7 in the NALS case, was identified both in the fast dynamics (i.e., in the generated constraints) and the slow dynamics (i.e., in the reactions that drive the system within these constraints).
The results reported in Tables 3, 4 and 5, regarding the TPI and II indices, show that the Lac oxidative metabolism in neurons (R18) has a larger impact than the oxidative metabolism in astrocytes (R17) in both the ANLS and NALS cases. This is consistent with the experimental finding that the Lac oxidative metabolism in neurons is more important than the respective metabolism in astrocytes [77, 135]. In addition, the identifications reported here in Table 5, show that the concentration of Laca,n is, among others, determined by the rates R7 and R9, which act towards increasing the rate of change of Laca,n. These rates have been reported to increase during neuronal activation in both ANLS (only R7) [16–18] and NALS (both R7 and R9) [19–21] cases, resulting in increasing Laca,n concentrations. This result is agreement with the elevated lactate levels that have been reported during neuronal activation [35, 137]. This feature has also been characterised as a moderate to slow process [12], which is also in agreement with our results, due to the aforementioned identifications in the slow modes of the system.
At this point it must be emphasized that the model employed here neglects metabolic processes which are currently considered important. Among them are (i) formation of pyruvate from glucose, (ii) oxidative phosphorylation in mitochondria, (iii) TCA cycle and (iv) consumption/production of ATP and NADH [8, 22, 88]. Clearly, it is possible that consideration of these processes might modify (i) the processes that are mainly responsible for the generation of the fast and slow dynamics and (ii) the conditions for which ANLS or NALS manifest. The degree to which these additional processes influence the response of the system, as it was examined here by analyzing algorithmically the fast and slow dynamics, is currently under investigation.
Finally, the ease by which the CSP tools lead to the determination of the factors that promote the development of ANLS or NALS must be highlighted. CSP allowed for these identifications by identifying the components of the model that are responsible for the fast and the slow dynamics, which are responsible for the development of the constraints in phase space (the fast) and for driving the system within these constraints (the slow).
Conclusions
A computational model describing the evolution of the brain lactate metabolism under neuronal activation was investigated using algorithmic tools of asymptotic analysis. The model introduced in Simpson et al. [9] and modified by Mangia et al. [10] was considered, which conforms to both the ANLS and NALS hypotheses. The analysis was carried out using the CSP method, which provides algorithmic tools for the analysis of multi-scale systems, so it is not hindered by the complexity and/or the size of the system under investigation and its tools can deliver systems-level understanding. CSP was employed in order to: (i) investigate the fast/slow dynamics of the system and the metabolic profiles of neurons and astrocytes during neuronal activation and (ii) compare the dynamics of the ANLS and NALS hypotheses, under normal and exercise conditions.
It was shown that the ANLS or NALS configuration develops depending on whether the Glc rate of transport from astrocytes to interstitium (R5) contributes or not to the generated constraints and the fast dynamics. In the ANLS case R5 participates in the equilibria of the 4th and 5th modes, which allow Glc in the astrocytes and neurons to be coupled directly with that in the serum. Lac formation from Glc is more intense in astrocytes (R7) than in neurons (R9), so that the typical for the ANLS case flow of Lac from astrocytes to neurons is generated. In the NALS case the participation of R5 to the equilibria of the 4th and 5th modes is significantly diminished, so that only Glc in neurons is coupled directly with that in the serum. Therefore, Lac formation from Glc is mainly taking place in neurons (R9), so that the typical for the NALS case flow of Lac from neurons to astrocytes is generated.
As demonstrated, CSP is not hindered by the size or the complexity of the mechanism under investigation, therefore more complicated and larger models of brain metabolism can be considered. In many Central Nervous System pathogenies (like Parkinson’s and Alzheimer’s diseases) or in metabolism-related diseases (e.g., diabetes), where the degeneration of the brain structures causes abnormalities in brain function, glucose and lactate metabolism is altered; see for example [7, 138, 139]. CSP could be effectively employed in order to identify optimal strategies for the control of the system towards desired outcomes, contributing thus to the development of therapeutic treatments.
Supporting information
S1 Text [pdf]
The brain lactate metabolism model.
S2 Text [pdf]
The computational singular perturbation methodology.
Zdroje
1. Howarth C, Gleeson P, Attwell D. Updated energy budgets for neural computation in the neocortex and cerebellum. Journal of Cerebral Blood Flow & Metabolism. 2012;32(7):1222–32. doi: 10.1038/jcbfm.2012.35
2. Falkowska A, Gutowska I, Goschorska M, Nowacki P, Chlubek D, Baranowska-Bosiacka I. Energy metabolism of the brain, including the cooperation between astrocytes and neurons, especially in the context of glycogen metabolism. International journal of molecular sciences. 2015;16(11):25959–81. doi: 10.3390/ijms161125939 26528968
3. Mergenthaler P, Lindauer U, Dienel GA, Meisel A. Sugar for the brain: the role of glucose in physiological and pathological brain function. Trends in neurosciences. 2013;36(10):587–97. doi: 10.1016/j.tins.2013.07.001 23968694
4. Pellerin L, Magistretti PJ. Glutamate uptake into astrocytes stimulates aerobic glycolysis: a mechanism coupling neuronal activity to glucose utilization. Proceedings of the National Academy of Sciences. 1994;91(22):10625–10629. doi: 10.1073/pnas.91.22.10625
5. Pellerin L, Magistretti PJ. Sweet sixteen for ANLS. Journal of Cerebral Blood Flow & Metabolism. 2012;32(7):1152–1166. doi: 10.1038/jcbfm.2011.149
6. Bélanger M, Allaman I, Magistretti PJ. Brain energy metabolism: focus on astrocyte-neuron metabolic cooperation. Cell metabolism. 2011;14(6):724–738. doi: 10.1016/j.cmet.2011.08.016 22152301
7. Newington JT, Harris RA, Cumming RC. Reevaluating metabolism in Alzheimer’s disease from the perspective of the astrocyte-neuron lactate shuttle model. Journal of neurodegenerative diseases. 2013;2013. doi: 10.1155/2013/234572 26316984
8. 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. Journal of Cerebral Blood Flow & Metabolism. 2010;30(3):586–602. doi: 10.1038/jcbfm.2009.232
9. 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(11):1766–91. doi: 10.1038/sj.jcbfm.9600521 17579656
10. 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. Journal of Neurochemistry. 2009;109:55–62. doi: 10.1111/j.1471-4159.2009.06003.x 19393009
11. Dienel GA. Brain lactate metabolism: the discoveries and the controversies. Journal of Cerebral Blood Flow & Metabolism. 2012;32(7):1107–1138. doi: 10.1038/jcbfm.2011.175
12. Chih CP, Roberts EL Jr. Energy substrates for neurons during neural activity: a critical review of the astrocyte-neuron lactate shuttle hypothesis. Journal of Cerebral Blood Flow & Metabolism. 2003;23(11):1263–1281. doi: 10.1097/01.WCB.0000081369.51727.6F
13. Patel AB, Lai JC, Chowdhury GM, Hyder F, Rothman DL, Shulman RG, Behar KL. Direct evidence for activity-dependent glucose phosphorylation in neurons with implications for the astrocyte-to-neuron lactate shuttle. Proceedings of the National Academy of Sciences. 2014:201403576. doi: 10.1073/pnas.1403576111
14. Lundgaard I, Li B, Xie L, Kang H, Sanggaard S, Haswell JD, Sun W, Goldman S, Blekot S, Nielsen M, Takano T. Direct neuronal glucose uptake heralds activity-dependent increases in cerebral metabolism. Nature communications. 2015;6:6807. doi: 10.1038/ncomms7807 25904018
15. Diaz-Garcia CM, Mongeon R, Lahmann C, Koveal D, Zucker H, Yellen G. Neuronal stimulation triggers neuronal glycolysis and not lactate uptake. Cell metabolism. 2017;26(2):361–74. doi: 10.1016/j.cmet.2017.06.021 28768175
16. Herrero-Mendez A, Almeida A, Fernández E, Maestre C, Moncada S, Bolaños JP. The bioenergetic and antioxidant status of neurons is controlled by continuous degradation of a key glycolytic enzyme by APC/C–Cdh1. Nature cell biology. 2009;11(6):747–752. doi: 10.1038/ncb1881 19448625
17. Patel JR, Brewer GJ. Age-related changes in neuronal glucose uptake in response to glutamate and β-amyloid. Journal of neuroscience research. 2003;72(4):527–536. doi: 10.1002/jnr.10602 12704814
18. Porras OH, Loaiza A, Barros LF. Glutamate mediates acute glucose transport inhibition in hippocampal neurons. Journal of Neuroscience. 2004;24(43):9669–9673. doi: 10.1523/JNEUROSCI.1882-04.2004 15509754
19. Erecińska M, Nelson D, Chance B. Depolarization-induced changes in cellular energy production. Proc. of the National Academy of Sciences. 1991;88(17):7600–7604. doi: 10.1073/pnas.88.17.7600
20. Kauppinen RA, Nicholls DG. Synaptosomal bioenergetics: the role of glycolysis, pyruvate oxidation and responses to hypoglycaemia. European journal of biochemistry. 1986;158(1):159–165. doi: 10.1111/j.1432-1033.1986.tb09733.x 2874024
21. Kauppinen RA, Taipale HT, Komulainen H. Interrelationships Between Glucose Metabolism, Energy State, and the Cytosolic Free Calcium Concentration in Cortical Synjaptosomes from the Guinea Pig. Journal of neurochemistry. 1989;53(3):766–771. doi: 10.1111/j.1471-4159.1989.tb11771.x 2503588
22. Jolivet R, Coggan JS, Allaman I, Magistretti PJ. Multi-timescale modeling of activity-dependent metabolic coupling in the neuron-glia-vasculature ensemble. PLoS computational biology. 2015 Feb 26;11(2):e1004036. doi: 10.1371/journal.pcbi.1004036 25719367
23. Juaristi I, Contreras P, González-Sánchez P, Pérez-Liébana I, González-Moreno L, Pardo B, del Arco A, Satrústegui J. The response to stimulation in neurons and astrocytes. Neurochemical research. 2019, https://doi.org/10.1007/s11064-019-02803-7. 31016552
24. Calì C, Tauffenberger A, Magistretti P. The strategic location of glycogen and lactate: from body energy reserve to brain plasticity. Frontiers in cellular neuroscience. 2019;13(82):1–7.
25. Wang Q, Hu Y, Wan J, Dong B, Sun J. Lactate: a novel signaling molecule in synaptic plasticity and drug addiction. BioEssays. 2019;41(8):1900008. doi: 10.1002/bies.201900008
26. Magistretti PJ, Allaman I. Lactate in the brain: from metabolic end-product to signalling molecule. Nature Reviews Neuroscience. 2018;19(4):235–249. doi: 10.1038/nrn.2018.19 29515192
27. Coggan JS, Keller D, Calì C, Lehväslaiho H, Markram H, Schürmann F, Magistretti PJ. Norepinephrine stimulates glycogenolysis in astrocytes to fuel neurons with lactate. PLOS Computational Biology. 2018;14(8):e1006392. doi: 10.1371/journal.pcbi.1006392 30161133
28. Barros LF, Joachim WD. Glucose and lactate supply to the synapse. Brain Research Reviews. 2010;63(1):149–159. doi: 10.1016/j.brainresrev.2009.10.002 19879896
29. Pellerin L, Bouzier-Sore AK, Aubert A, Serres S, Merle M, Costalat R, Magistretti PJ. Activity-dependent regulation of energy metabolism by astrocytes: an update. Glia. 2007;55(12):1251–1262. doi: 10.1002/glia.20528 17659524
30. Hyder F, Patel AB, Gjedde A, Rothman DL, Behar KL, Shulman RG. Neuronal–glial glucose oxidation and glutamatergic–GABAergic function. Journal of Cerebral Blood Flow & Metabolism. 2006;26(7):865–877. doi: 10.1038/sj.jcbfm.9600263
31. Jolivet R, Allaman I, Pellerin L, Magistretti PJ, Weber B. Comment on recent modeling studies of astrocyte-neuron metabolic interactions. Journal of Cerebral Blood Flow & Metabolism. 2010;30(12):1982–1986. doi: 10.1038/jcbfm.2010.132
32. Lundengård K, Cedersund G, Sten S, Leong F, Smedberg A, Elinder F, Engström M. Mechanistic mathematical modeling tests hypotheses of the neurovascular coupling in fMRI. PLOS Computational Biology. 2016;12(6):e1004971. doi: 10.1371/journal.pcbi.1004971 27310017
33. Rouach N, Koulakoff A, Abudara V, Willecke K, Giaume C. Astroglial metabolic networks sustain hippocampal synaptic transmission. American Association for the Advancement of Science. 2008;322(5907):1551–1555. doi: 10.1126/science.1164022
34. Kasischke KA, Vishwasrao HD, Fisher PJ, Zipfel WR, Webb WW. Neural activity triggers neuronal oxidative metabolism followed by astrocytic glycolysis. Science. 2004;305(5680):99–103. doi: 10.1126/science.1096485 15232110
35. Prichard J, Rothman D, Novotny E, Petroff O, Kuwabara T, Avison M, Howseman A, Hanstock C, Shulman R. Lactate rise detected by 1H NMR in human visual cortex during physiologic stimulation. Proceedings of the National Academy of Sciences of the United States of America. 1991;88(13):5829. doi: 10.1073/pnas.88.13.5829 2062861
36. Miyamoto K, Ishikura K, Kume K, Ohsawa M. Astrocyte-neuron lactate shuttle sensitizes nociceptive transmission in the spinal cord. Glia. 2019;67:27–36. doi: 10.1002/glia.23474 30430652
37. Mächler P, Wyss MT, Elsayed M, Stobart J, Gutierrez R, von Faber-Castell A, Kaelin V, Zuend M, San Martín A, Romero-Gómez I, Baeza-Lehnert F, Lengacher S, Schneider BL, Aebischer P, Magistretti PJ, Barros LF, Weber B. In vivo evidence for a lactate gradient from astrocytes to neurons. Cell Metabolism. 2016;23(1):94–102. doi: 10.1016/j.cmet.2015.10.010 26698914
38. Ruminot I, Schmälzle J, Leyton B, Barros LF, Deitmer JW. Tight coupling of astrocyte energy metabolism to synaptic activity revealed by genetically encoded FRET nanosensors in hippocampal tissue. Journal of Cerebral Blood Flow & Metabolism. 2019;39(3):513–523. doi: 10.1177/0271678X17737012
39. Liu L, MacKenzie KR, Putluri N, Maletić-Savatić M, Bellen HJ. The Glia-Neuron Lactate Shuttle and Elevated ROS Promote Lipid Synthesis in Neurons and Lipid Droplet Accumulation in Glia via APOE/D. Cell Metabolism. 2017;26(5):719–737. doi: 10.1016/j.cmet.2017.08.024 28965825
40. Cholet N, Pellerin L, Welker E, Lacombe P, Seylaz J, Magistretti PJ, Bonvento G. Local injection of antisense oligonucleotides targeted to the glial glutamate transporter GLAST decreases the metabolic response to somatosensory activation. Journal of Cerebral Blood Flow & Metabolism. 2001;21(4):404–412. doi: 10.1097/00004647-200104000-00009
41. Voutsinos-Porche B, Bonvento G, Tanaka K, Steiner P, Welker E, Chatton J-Y, Magistretti PJ, Pellerin L. Glial glutamate transporters mediate a functional metabolic crosstalk between neurons and astrocytes in the mouse developing cortex. Neuron. 2003;37(2):275–286. doi: 10.1016/s0896-6273(02)01170-4 12546822
42. Herard A-S, Dubois A, Escartin C, Tanaka K, Delzescaux T, Hantraye P, Bonvento G. Decreased metabolic response to visual stimulation in the superior colliculus of mice lacking the glial glutamate transporter GLT-1. European Journal of Neuroscience. 2005;22(7):1807–1811. doi: 10.1111/j.1460-9568.2005.04346.x 16197522
43. Chuquet J, Quilichini P, Nimchinsky EA, Buzsáki G. Predominant enhancement of glucose uptake in astrocytes versus neurons during activation of the somatosensory cortex. Journal of Neuroscience. 2010;30(45):15298–15303. doi: 10.1523/JNEUROSCI.0762-10.2010 21068334
44. Lin A-L, Fox PT, Hardies J, Duong TQ, Gao J-H. Nonlinear coupling between cerebral blood flow, oxygen consumption, and ATP production in human visual cortex. Proceedings of the National Academy of Sciences. 2010;107(418):8446–8451. doi: 10.1073/pnas.0909711107
45. Dienel GA. The “protected” glucose transport through the astrocytic endoplasmic reticulum is too slow to serve as a quantitatively-important highway for nutrient delivery. Journal of Neuroscience Research. 2019;97:854–862. doi: 10.1002/jnr.24432 31050047
46. Díaz-García CM, Yellen G. Neurons rely on glucose rather than astrocytic lactate during stimulation. Journal of Neuroscience Research. 2019;97:883–889. doi: 10.1002/jnr.24374 30575090
47. Dienel GA. Lack of appropriate stoichiometry: Strong evidence against an energetically important astrocyte–neuron lactate shuttle in brain. Journal of Neuroscience Research. 2017;95:2103–2125. doi: 10.1002/jnr.24015 28151548
48. Bak LK, Walls AB. CrossTalk opposing view: lack of evidence supporting an astrocyte-to-neuron lactate shuttle coupling neuronal activity to glucose utilisation in the brain. Journal of Physiology. 2018;596(3):351–353. doi: 10.1113/JP274945 29292507
49. Dienel GA. Lactate shuttling and lactate use as fuel after traumatic brain injury: metabolic considerations. Journal of Cerebral Blood Flow & Metabolism. 2014;34(11):1736–1748. doi: 10.1038/jcbfm.2014.153
50. Dienel GA. Fueling and imaging brain activation. JASN neuro. 2012;4(5):267–321.
51. Hertz L, Dienel GA. Lactate transport and transporters: general principles and functional roles in brain cells. Journal of neuroscience research. 2005;79(1-2):11–18. doi: 10.1002/jnr.20294 15586354
52. Hertz L. The astrocyte-neuron lactate shuttle: a challenge of a challenge. Journal of Cerebral Blood Flow & Metabolism. 2004;24(11):1241–1248. doi: 10.1097/00004647-200411000-00008
53. Hertz L, Dienel GA. Energy metabolism in the brain. International review of neurobiology. 2002;51:1–102. doi: 10.1016/s0074-7742(02)51003-5 12420357
54. Dienel GA, Hertz L. Glucose and lactate metabolism during brain activation. Journal of neuroscience research. 2001;66(5):824–838. doi: 10.1002/jnr.10079 11746408
55. Mangia S, Garreffa G, Bianciardi M, Giove F, Di Salle F, Maraviglia B. The aerobic brain: lactate decrease at the onset of neural activity. Neuroscience. 2003;118(1):7–10. doi: 10.1016/s0306-4522(02)00792-3 12676131
56. Maher F, Simpson IA. Modulation of expression of glucose transporters GLUT3 and GLUT1 by potassium and N-methyl-D-aspartate in cultured cerebellar granule neurons. Molecular and Cellular Neuroscience. 1994;5(4):369–375. doi: 10.1006/mcne.1994.1044 7804607
57. Castro MA, Pozo M, Cortés C, García MA, Concha II, Nualart F. Intracellular ascorbic acid inhibits transport of glucose by neurons, but not by astrocytes. Journal of neurochemistry. 2007;102(3):773–782. doi: 10.1111/j.1471-4159.2007.04631.x 17630983
58. Weisová P, Concannon CG, Devocelle M, Prehn JHM, Ward MW. Regulation of glucose transporter 3 surface expression by the AMP-activated protein kinase mediates tolerance to glutamate excitation in neurons. Journal of Neuroscience. 2009;29(9):2997–3008. doi: 10.1523/JNEUROSCI.0354-09.2009 19261894
59. Ferreira JM, Burnett AL, Rameau GA. Activity-dependent regulation of surface glucose transporter-3 Journal of Neuroscience. 2011;31(6):1991–1999. doi: 10.1523/JNEUROSCI.1850-09.2011 21307237
60. Almeida A, Almeida J, Bolaños JP, Moncada S. Different responses of astrocytes and neurons to nitric oxide: the role of glycolytically generated ATP in astrocyte protection. Proceedings of the National Academy of Sciences. 2001;98(26):15294–15299. doi: 10.1073/pnas.261560998
61. Dienel GA, Cruz NF. Astrocyte activation in working brain: energy supplied by minor substrates. Neurochemistry international. 2006;48(6-7):586–595. doi: 10.1016/j.neuint.2006.01.004 16513214
62. Gegg ME, Beltran B, Salas-Pino S, Bolanos JP, Clark JB, Moncada S, Heales SJR. Differential effect of nitric oxide on glutathione metabolism and mitochondrial function in astrocytes and neurones: implications for neuroprotection/neurodegeneration? Journal of neurochemistry. 2003;86(1):228–237. doi: 10.1046/j.1471-4159.2003.01821.x 12807442
63. Gjedde A, Marrett S. Glycolysis in neurons, not astrocytes, delays oxidative metabolism of human visual cortex during sustained checkerboard stimulation in vivo. Journal of Cerebral Blood Flow & Metabolism. 2001;21(12):1384–1392. doi: 10.1097/00004647-200112000-00002
64. Bak LK, Walls AB, 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. Journal of neurochemistry. 2009;109(1):87–93. doi: 10.1111/j.1471-4159.2009.05943.x 19393013
65. Cunnane S, Nugent S, Roy M, Courchesne-Loyer A, Croteau E, Tremblay S, Castellano A, Pifferi F, Bocti C, Paquet N and others. Brain fuel metabolism, aging, and Alzheimer’s disease. Nutrition. 2011;27(1):3–20. doi: 10.1016/j.nut.2010.07.021 21035308
66. Harris DL, Weston PJ, Harding JE. Lactate, rather than ketones, may provide alternative cerebral fuel in hypoglycaemic newborns. Archives of Disease in Childhood-Fetal and Neonatal Edition. 2015;100(2):F161–F164. doi: 10.1136/archdischild-2014-306435 25189167
67. Thudichum J. A. Treatise on the Chemical Constitution of the Brain. 1884. London: Bailliere, Tindall and Cox. 1962:262.
68. Himwich HE. Brain metabolism and cerebral disorders. Williams & Wilkins; 1951.
69. Itoh Y, Esaki T, Shimoji K, Cook M, Law MJ, Kaufman E, Sokoloff L. Dichloroacetate effects on glucose and lactate oxidation by neurons and astroglia in vitro and on glucose utilization by brain in vivo. Proceedings of the National Academy of Sciences. 2003;100(8):4879–84. doi: 10.1073/pnas.0831078100
70. Nakai T, Matsuo K, Kato C, Takehara Y, Isoda H, Moriya T, Okada T, Sakahara H. Post-stimulus response in hemodynamics observed by functional magnetic resonance imaging-Difference between the primary sensorimotor area and the supplementary motor area. Magnetic resonance imaging. 2000;18(10):1215–9. doi: 10.1016/s0730-725x(00)00217-4 11167041
71. Ances BM, Buerk DG, Greenberg JH, Detre JA. Temporal dynamics of the partial pressure of brain tissue oxygen during functional forepaw stimulation in rats. Neuroscience letters. 2001;306(1-2):106–10. doi: 10.1016/s0304-3940(01)01868-7 11403969
72. Fellows LK, Boutelle MG, Fillenz M. Physiological stimulation increases nonoxidative glucose metabolism in the brain of the freely moving rat. Journal of neurochemistry. 1993;60(4):1258–63. doi: 10.1111/j.1471-4159.1993.tb03285.x 8455025
73. Fray AE, Forsyth RJ, Boutelle MG, Fillenz M. The mechanisms controlling physiologically stimulated changes in rat brain glucose and lactate: a microdialysis study. The Journal of physiology. 1996;496(1):49–57. doi: 10.1113/jphysiol.1996.sp021664 8910195
74. Vega C, Martiel JL, Drouhault D, Burckhart MF, Coles JA. Uptake of locally applied deoxyglucose, glucose and lactate by axons and Schwann cells of rat vagus nerve. The Journal of physiology. 2003;546(2):551–64. doi: 10.1113/jphysiol.2002.029751 12527741
75. Kuhr WG, van den Berg CJ, Korf J. In vivo identification and quantitative evaluation of carrier-mediated transport of lactate at the cellular level in the striatum of conscious, freely moving rats. Journal of Cerebral Blood Flow & Metabolism. 1988;8(6):848–56. doi: 10.1038/jcbfm.1988.142
76. Hu Y, Wilson GS A temporary local energy pool coupled to neuronal activity: fluctuations of extracellular lactate levels in rat brain monitored with rapid-response enzyme-based sensor. Journal of neurochemistry. 1997;69(4):1484–90. doi: 10.1046/j.1471-4159.1997.69041484.x 9326277
77. Lebon V, Petersen KF, Cline GW, Shen J, Mason GF, Dufour S, Behar KL, Shulman GI, Rothman DL. Astroglial contribution to brain energy metabolism in humans revealed by 13C nuclear magnetic resonance spectroscopy: elucidation of the dominant pathway for neurotransmitter glutamate repletion and measurement of astrocytic oxidative metabolism. Journal of Neuroscience. 2002;22(5):1523–31. doi: 10.1523/JNEUROSCI.22-05-01523.2002 11880482
78. Buxton RB, Wong EC, Frank LR. Dynamics of blood flow and oxygenation changes during brain activation: the balloon model. Magnetic resonance in medicine. 1998;39(6):855–64. doi: 10.1002/mrm.1910390602 9621908
79. Gruetter R, Seaquist ER, Ugurbil K. A mathematical model of compartmentalized neurotransmitter metabolism in the human brain. American Journal of Physiology-Endocrinology And Metabolism. 2001;281(1):E100–12. doi: 10.1152/ajpendo.2001.281.1.E100 11404227
80. Aubert A, Costalat R, Valabrègue R. Modelling of the coupling between brain electrical activity and metabolism. Acta biotheoretica. 2001;49(4):301–26. doi: 10.1023/a:1014286728421 11804241
81. Aubert A, Costalat R. A model of the coupling between brain electrical activity, metabolism, and hemodynamics: application to the interpretation of functional neuroimaging. Neuroimage. 2002;17(3):1162–81. doi: 10.1006/nimg.2002.1224 12414257
82. Aubert A, Costalat R. Interaction between astrocytes and neurons studied using a mathematical model of compartmentalized energy metabolism. Journal of Cerebral Blood Flow & Metabolism. 2005; 25(11): 1476–90. doi: 10.1038/sj.jcbfm.9600144
83. Aubert A, Costalat R, Magistretti PJ, Pellerin L. Brain lactate kinetics: modeling evidence for neuronal lactate uptake upon activation. Proceedings of the National Academy of Sciences. 2005 Nov 8;102(45):16448–53. doi: 10.1073/pnas.0505427102
84. Aubert A, Pellerin L, Magistretti PJ, Costalat R. A coherent neurobiological framework for functional neuroimaging provided by a model integrating compartmentalized energy metabolism. Proceedings of the National Academy of Sciences. 2007 Mar 6;104(10):4188–93. doi: 10.1073/pnas.0605864104
85. Cloutier M, Bolger FB, Lowry JP, Wellstead P. An integrative dynamic model of brain energy metabolism using in vivo neurochemical measurements. Journal of computational neuroscience. 2009;27(3):391. doi: 10.1007/s10827-009-0152-8 19396534
86. Occhipinti R, Puchowicz MA, LaManna JC, Somersalo E, Calvetti D. Statistical analysis of metabolic pathways of brain metabolism at steady state. Annals of biomedical engineering. 2007;35(6):886–902. doi: 10.1007/s10439-007-9270-5 17385046
87. Occhipinti R, Somersalo E, Calvetti D. Astrocytes as the glucose shunt for glutamatergic neurons at high activity: an in silico study. Journal of neurophysiology. 2009;101(5):2528–38. doi: 10.1152/jn.90377.2008 18922953
88. Occhipinti R, Somersalo E, Calvetti D. Energetics of inhibition: insights with a computational model of the human GABAergic neuron–astrocyte cellular complex. Journal of Cerebral Blood Flow & Metabolism. 2010;30(11):1834–46. doi: 10.1038/jcbfm.2010.107
89. Calvetti D, Somersalo E. Menage a trois: the role of neurotransmitters in the energy metabolism of astrocytes, glutamatergic, and GABAergic neurons. Journal of Cerebral Blood Flow & Metabolism. 2012 Aug;32(8):1472–83. doi: 10.1038/jcbfm.2012.31
90. Blanchard S, Papadopoulo T, Bénar C-G, Voges N, Clerc M, Benali H, Warnking J, David O, Wendling F. Relationship between flow and metabolism in BOLD signals: insights from biophysical models. Brain topography, 2011 Nov 6;24(1):40–53 doi: 10.1007/s10548-010-0166-6 21057867
91. Winter F, Bludszuweit-Philipp C, Wolkenhauer O. Mathematical analysis of the influence of brain metabolism on the BOLD signal in Alzheimer’s disease. Journal of Cerebral Blood Flow & Metabolism, 2018;38(2):304–316 doi: 10.1177/0271678X17693024
92. Çakιr T, Alsan S, Saybaşιlι H, Akιn A, Ülgen K. Reconstruction and flux analysis of coupling between metabolic pathways of astrocytes and neurons: application to cerebral hypoxia. Theoretical Biology and Medical Modelling, 2007 Dec 10;4(1):48–66 doi: 10.1186/1742-4682-4-48
93. Massucci FA, DiNuzzo M, Giove F, Maraviglia B, Castillo IP, Marinari E, De Martino A. Energy metabolism and glutamate-glutamine cycle in the brain: a stoichiometric modeling perspective. BMC systems biology, 2013;7(1):103–117 doi: 10.1186/1752-0509-7-103 24112710
94. DiNuzzo M, Giove F, Maraviglia B, Mangia S. Computational flux balance analysis predicts that stimulation of energy metabolism in astrocytes and their metabolic interactions with neurons depend on uptake of K+ rather than glutamate. Neurochemical research, 2016 Sep 14;42(1):202–216 doi: 10.1007/s11064-016-2048-0 27628293
95. Calvetti D, Somersalo E. Dynamic activation model for a glutamatergic neurovascular unit. Journal of theoretical biology, 2010 Dec 19;274(1):12–29 doi: 10.1016/j.jtbi.2010.12.007 21176783
96. Somersalo E, Cheng Y, Calvetti D. The metabolism of neurons and astrocytes through mathematical models Annals of biomedical engineering, 2012 Sep 22;40(11):2328–2344 doi: 10.1007/s10439-012-0643-z 23001357
97. Frenklach M. Modeling. In: Gardiner WC, editor. Combustion Chemistry. Springer-Verlag, New York; 1984. p. 423–454.
98. Vajda S, Valko P, Turanyi T Principal component analysis of kinetic models. International Journal of Chemical Kinetics. 1985;17(1):55–81. doi: 10.1002/kin.550170107
99. Lam S, Goussis D. The CSP method for simplifying kinetics. International Journal of Chemical Kinetics. 1994;26(4):461–486. doi: 10.1002/kin.550260408
100. Lam S, Goussis D. Understanding complex chemical kinetics with computational singular perturbation. In: Symposium (International) on Combustion. 1989;22:931–941.
101. Hadjinicolaou M, Goussis DA. Asymptotic solution of stiff PDEs with the CSP method: the reaction diffusion equation. SIAM J. Sci. Comp., 1999;20:781–819. doi: 10.1137/S1064827596303995
102. Prager J, Najm HN, Valorani M, Goussis DA. Structure of n-heptane/air triple flames in partially-premixed mixing layers. Combustion and Flame. 2011;158(11):2128–44. doi: 10.1016/j.combustflame.2011.03.017
103. Goussis DA, Valorani M. An efficient iterative algorithm for the approximation of the fast and slow dynamics of stiff systems. Journal of Computational Physics. 2006;214(1):316–46. doi: 10.1016/j.jcp.2005.09.019
104. Neophytou MK, Goussis DA, Mastorakos E, Britter RE. The conceptual development of a simple scale-adaptive reactive pollutant dispersion model. Atmospheric Environment. 2005;39(15):2787–94. doi: 10.1016/j.atmosenv.2004.12.025
105. Kourdis PD, Goussis DA. Glycolysis in saccharomyces cerevisiae: algorithmic exploration of robustness and origin of oscillations. Mathematical biosciences. 2013;243(2):190–214. doi: 10.1016/j.mbs.2013.03.002 23517854
106. Surovtsova I, Simus N, Hübner K, Sahle S, Kummer U. Simplification of biochemical models: a general approach based on the analysis of the impact of individual species and reactions on the systems dynamics. BMC systems biology. 2012;6(1):14. doi: 10.1186/1752-0509-6-14 22390191
107. Samant A, Ogunnaike BA, Vlachos DG. A hybrid multiscale Monte Carlo algorithm (HyMSMC) to cope with disparity in time scales and species populations in intracellular networks. BMC bioinformatics. 2007;8(1):175. doi: 10.1186/1471-2105-8-175 17524148
108. Patsatzis DG, Goussis DA. A new Michaelis-Menten equation valid everywhere multi-scale dynamics prevails. Mathematical biosciences. 2019;315:108220. doi: 10.1016/j.mbs.2019.108220 31255632
109. Patsatzis DG, Maris DT, Goussis DA. Asymptotic analysis of a target-mediated drug disposition model: algorithmic and traditional approaches. Bulletin of mathematical biology, 2016;78(6):1121–1161. doi: 10.1007/s11538-016-0176-y 27271122
110. Michalaki LI, Goussis DA. Asymptotic analysis of a TMDD model: when a reaction contributes to the destruction of its product. Journal of mathematical biology. 2018:1–35.
111. Fenichel N. Geometric singular perturbation theory for ordinary differential equations. Journal of differential equations. 1979;31(1):53–98. doi: 10.1016/0022-0396(79)90152-9
112. Kaper TJ. An introduction to geometric methods and dynamical systems theory for singular perturbation problems. In: Cronin J, Robert J, O’Malley E (eds) Analyzing multiscale phenomena using singular perturbation methods. Proceedings of symposia in applied mathematics. 1999;56(1):85–131.
113. Kuehn C. Multiple time scale dynamics. Springer; 2015.
114. Zagaris A, Kaper HG, Kaper TJ. Analysis of the computational singular perturbation reduction method for chemical kinetics. Journal of Nonlinear Science. 2004;14(1):59–91. doi: 10.1007/s00332-003-0582-9
115. Zagaris A, Kaper HG, Kaper TJ. Fast and slow dynamics for the computational singular perturbation method. Multiscale Modeling & Simulation. 2004;2(4):613–638. doi: 10.1137/040603577
116. Zagaris A, Kaper HG, Kaper TJ. Two perspectives on reduction of ordinary differential equations. Mathematische Nachrichten. 2005;278(12-13):1629–1642. doi: 10.1002/mana.200410328
117. Kaper HG, Kaper TJ, Zagaris A. Geometry of the computational singular perturbation method. Mathematical modelling of natural phenomena. 2015;10(3):16–30. doi: 10.1051/mmnp/201510303
118. Hek G. Geometric singular perturbation theory in biological practice. Journal of mathematical biology. 2010;60(3):347–386. doi: 10.1007/s00285-009-0266-7 19347340
119. Auger P, de La Parra RB, Poggiale J-C, Sánchez E, Sanz L. Aggregation methods in dynamical systems and applications in population and community dynamics. Physics of Life Reviews. 2008;5(2):79–105. doi: 10.1016/j.plrev.2008.02.001
120. Kosiuk I, Szmolyan P. Geometric analysis of the Goldbeter minimal model for the embryonic cell cycle. Journal of mathematical biology. 2016;72(5):1337–1368. doi: 10.1007/s00285-015-0905-0 26100376
121. Popović N, Marr C, Swain PS. A geometric analysis of fast-slow models for stochastic gene expression. Journal of mathematical biology. 2016;72(1-2):87–122. doi: 10.1007/s00285-015-0876-1 25833185
122. Vanselow A, Wieczorek S, Feudel U. When very slow is too fast-collapse of a predator-prey system. Journal of theoretical biology. 2019;479:64–72. doi: 10.1016/j.jtbi.2019.07.008 31302207
123. Quistorff B, Secher NH, Van Lieshout JJ. Lactate fuels the human brain during exercise. The FASEB Journal. 2008;22(10):3443–3449. doi: 10.1096/fj.08-106104 18653766
124. Rasmussen P, Wyss MT, Lundby C. Cerebral glucose and lactate consumption during cerebral activation by physical activity in humans. The FASEB Journal. 2011;25(9):2865–2873. doi: 10.1096/fj.11-183822 21602451
125. Lam S, Goussis D. Conventional asymptotics and computational singular perturbation for simplified kinetics modelling. In: Reduced kinetic mechanisms and asymptotic approximations for methane-air flames; 1991. p. 227–242.
126. Tingas EA, Kyritsis DC, Goussis DA. Ignition delay control of DME/air and EtOH/air homogeneous autoignition with the use of various additives. Fuel. 2016;169:15–24. doi: 10.1016/j.fuel.2015.11.081
127. Goussis DA, Najm HN. Model reduction and physical understanding of slowly oscillating processes: the circadian cycle. Multiscale Modeling & Simulation. 2006;5(4):1297–1332. doi: 10.1137/060649768
128. Valorani M, Najm HN, Goussis DA. CSP analysis of a transient flame-vortex interaction: time scales and manifolds. Combustion and Flame. 2003;134(1):35–53. doi: 10.1016/S0010-2180(03)00067-1
129. Goussis DA, Skevis G. Nitrogen chemistry controlling steps in methane-air premixed flames. In: Bathe KJ, editor. Computational Fluid and Solid Mechanics. Elsevier, Amsterdam; 2005. p. 650–653.
130. Diamantis DJ, Mastorakos E, Goussis DA. H2/air autoignition: The nature and interaction of the developing explosive modes. Combustion Theory and Modelling. 2015;19(3):382–433. doi: 10.1080/13647830.2015.1027273
131. Goussis D, Lam S. A study of homogeneous methanol oxidation kinetics using CSP. In: Symposium (International) on Combustion. vol. 24. Elsevier; 1992. p. 113–120.
132. Goussis DA. Quasi steady state and partial equilibrium approximations: their relation and their validity. Combustion Theory and Modelling. 2012;16(5):869–926. doi: 10.1080/13647830.2012.680502
133. Kourdis PD, Palasantza AG, Goussis DA. Algorithmic asymptotic analysis of the NF-κB signaling system. Computers & Mathematics with Applications. 2013;65(10):1516–1534. doi: 10.1016/j.camwa.2012.11.004
134. Valorani M, Creta F, Donato F, Najm H, Goussis D. A CSP-based skeletal mechanism generation procedure: auto-ignition and premixed laminar flames in n-heptane/air mixtures. In: ECCOMAS CFD 2006: Proceedings of the European Conference on Computational Fluid Dynamics, Egmond aan Zee, The Netherlands, September 5-8, 2006. Delft University of Technology; European Community on Computational Methods in Applied Sciences (ECCOMAS); 2006.
135. Bouzier-Sore AK, Voisin P, Bouchaud V, Bezancon E, Franconi JM, Pellerin L. Competition between glucose and lactate as oxidative energy substrates in both neurons and astrocytes: a comparative NMR study. European Journal of Neuroscience. 2006;24(6):1687–1694. doi: 10.1111/j.1460-9568.2006.05056.x 17004932
136. Lovatt D, Sonnewald U, Waagepetersen HS, Schousboe A, He W, Lin JHC, et al. The transcriptome and metabolic gene signature of protoplasmic astrocytes in the adult murine cortex. Journal of Neuroscience. 2007;27(45):12255–12266. doi: 10.1523/JNEUROSCI.3404-07.2007 17989291
137. Ueki M, Linn F, Hossmann K-A. Functional activation of cerebral blood flow and metabolism before and after global ischemia of rat brain. Journal of Cerebral Blood Flow & Metabolism. 1988;8(4):486–494. doi: 10.1038/jcbfm.1988.89
138. Bohnen NI, Koeppe RA, Minoshima S, Giordani B, Albin RL, Frey KA, Kuhl DE. Cerebral glucose metabolic features of Parkinson disease and incident dementia: longitudinal study. Journal of Nuclear Medicine. 2011;52(6):848–855. doi: 10.2967/jnumed.111.089946 21571793
139. Action to Control Cardiovascular Risk in Diabetes Study Group. Effects of intensive glucose lowering in type 2 diabetes. New England journal of medicine. 2008;358(24):2545–2559. doi: 10.1056/NEJMoa0802743 18539917
Článek vyšel v časopise
PLOS One
2019 Číslo 12
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Methylsulfonylmethane increases osteogenesis and regulates the mineralization of the matrix by transglutaminase 2 in SHED cells
- Oregano powder reduces Streptococcus and increases SCFA concentration in a mixed bacterial culture assay
- The characteristic of patulous eustachian tube patients diagnosed by the JOS diagnostic criteria
- Parametric CAD modeling for open source scientific hardware: Comparing OpenSCAD and FreeCAD Python scripts