Computational singular perturbation analysis of brain lactate metabolism

• Published in: PloS one Vol. 14; no. 12; p. e0226094
• Main Authors: Patsatzis, Dimitris G, Tingas, Efstathios-Al, Goussis, Dimitris A, Sarathy, S Mani
• Format: Journal Article
• Language: English
• Published: United States Public Library of Science 17.12.2019; Public Library of Science (PLoS)
• Subjects: Activation; Algorithms; Applied mathematics; Astrocytes; Astrocytes - cytology; Astrocytes - metabolism; Biology and Life Sciences; Blood flow; Brain; Brain - metabolism; Brain research; Computational neuroscience; Computer simulation; Energy; Humans; Hypotheses; Lactic acid; Lactic Acid - metabolism; Mathematical models; Metabolism; Metabolites; Models, Theoretical; Neurogenesis; Neurons; Neurons - cytology; Neurons - metabolism; Perturbation methods; Physical Sciences; Research and Analysis Methods; Singular perturbation; Spacecraft components; Validity; Saudi Arabia
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0226094

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.