The metabolism of neurons and astrocytes through mathematical models

Lactate's behavioral switch in the brain: An in-silico model

Emerging evidence emphasizes lactate's involvement in both physiological processes (energy metabolism, memory, etc.) and disease (traumatic brain injury, epilepsy, etc.). Furthermore, the usefulness of mathematical modeling in deciphering underlying dynamics of the brain to investigate lactate roles and mechanisms of action has been well established. Here, we analyze a novel mathematical model of brain lactate exchanges between four compartments: neurons, astrocytes, capillaries, and extracellular space. A system of four ordinary differential equations is proposed to explain interactions between these compartments. We first optimize and analyze the model's parameters under normal, resting state conditions, and then use it to simulate changes linked to elevated arterial lactate. Our results show that even though increased arterial lactate results in increased uptake by astrocytes and release to the extracellular space, it cannot strongly recover the initial drop in neuronal lactate concentration. Also, we show that the direction of lactate transport between the compartments is influenced by the maximum astrocyte production rate and the transport rate between astrocytes and extracellular space.

Virtual Populations for Quantitative Systems Pharmacology Models

Quantitative systems pharmacology (QSP) places an emphasis on dynamic systems modeling, incorporating considerations from systems biology modeling and pharmacodynamics. The goal of QSP is often to quantitatively predict the effects of clinical therapeutics, their combinations, and their doses on clinical biomarkers and endpoints. In order to achieve this goal, strategies for incorporating clinical data into model calibration are critical. Virtual population (VPop) approaches facilitate model calibration while faced with challenges encountered in QSP model application, including modeling a breadth of clinical therapies, biomarkers, endpoints, utilizing data of varying structure and source, capturing observed clinical variability, and simulating with models that may require more substantial computational time and resources than often found in pharmacometrics applications. VPops are frequently developed in a process that may involve parameterization of isolated pathway models, integration into a larger QSP model, incorporation of clinical data, calibration, and quantitative validation that the model with the accompanying, calibrated VPop is suitable to address the intended question or help with the intended decision. Here, we introduce previous strategies for developing VPops in the context of a variety of therapeutic and safety areas: metabolic disorders, drug-induced liver injury, autoimmune diseases, and cancer. We introduce methodological considerations, prior work for sensitivity analysis and VPop algorithm design, and potential areas for future advancement. Finally, we give a more detailed application example of a VPop calibration algorithm that illustrates recent progress and many of the methodological considerations. In conclusion, although methodologies have varied, VPop strategies have been successfully applied to give valid clinical insights and predictions with the assistance of carefully defined and designed calibration and validation strategies. While a uniform VPop approach for all potential QSP applications may be challenging given the heterogeneity in use considerations, we anticipate continued innovation will help to drive VPop application for more challenging cases of greater scale while developing new rigorous methodologies and metrics.

MITODYN: An Open Source Software for Quantitative Modeling of Mitochondrial and Cellular Energy Metabolic Flux Dynamics in Health and Disease

Mitochondrial respiratory chain (RC) transforms the reductive power of NADH or FADH2 oxidation into a proton gradient between the matrix and cytosolic sides of the inner mitochondrial membrane, that ATP synthase uses to generate ATP. This process constitutes a bridge between carbohydrates' central metabolism and ATP-consuming cellular functions. Moreover, the RC is responsible for a large part of reactive oxygen species (ROS) generation that play signaling and oxidizing roles in cells. Mathematical methods and computational analysis are required to understand and predict the possible behavior of this metabolic system. Here we propose a software tool that helps to analyze individual steps of respiratory electron transport in their dynamics, thus deepening understanding of the mechanism of energy transformation and ROS generation in the RC. This software's core is a kinetic model of the RC represented by a system of ordinary differential equations (ODEs). This model enables the analysis of complex dynamic behavior of the RC, including multistationarity and oscillations. The proposed RC modeling method can be applied to study respiration and ROS generation in various organisms and naturally extended to explore carbohydrates' metabolism and linked metabolic processes.

Brain energetics, mitochondria, and traumatic brain injury

We review current thinking about, and draw connections between, brain energetics and metabolism, and between mitochondria and traumatic brain injury. Energy is fundamental to proper brain function. Its creation in a useful form for neurons and glia, and consistently in response to the brain's high energy needs, is critical for physiological pathways. Dysfunction in the mechanisms of energy production is at the center of neurological and neuropsychiatric pathologies. We examine the connections between energetics and mitochondria - the organelle responsible for almost all the energy production in the cell - and how secondary pathologies in traumatic brain injury result from energetic dysfunction. This paper interweaves these topics, a necessity since they are closely coupled, and identifies where there exist a lack of understanding and of data. In addition to summarizing current thinking in these disciplines, our goal is to suggest a framework for the mathematical modeling of mechanisms and pathways based on optimal energetic decisions.

Contributions of synaptic and astrocyte physiology to the anaesthetised encephalogram revealed using a computational model

BACKGROUND

General anaesthesia is known to enhance inhibitory synaptic transmission to produce characteristic effects on the EEG and reduction in brain metabolism secondary to reduced neuronal activity. Evidence suggests that anaesthesia might have a direct effect on synaptic metabolic processes, and this relates to anaesthesia sensitivity. We explored elements of synaptic transmission looking for possible contributions to the anaesthetised EEG and how it may modulate anaesthesia sensitivity.

METHODS

We developed a Hodgkin-Huxley-type neural network computer simulation capable of mimicking anaesthetic prolongation of gamma-aminobutyric acid (GABA)ergic inhibitory postsynaptic potentials (IPSPs), and capable of altering postsynaptic ion homeostasis and neurotransmitter recycling. We examined their interactions on simulated electrocorticography (sECoG), and compared these with published anaesthesia EEG spectra.

RESULTS

The sECoG spectra from the model were comparable with published normal awake EEG spectra. Prolongation of IPSP duration in the model caused inhibition of high frequencies and saturation of low frequencies with a peak in keeping with current evidence. IPSP prolongation alone was unable to reproduce alpha rhythms or the generalised increase in EEG power found with anaesthesia. Adding inhibition of postsynaptic ion homeostasis to IPSP prolongation helped retain alpha rhythms, increased sECoG power, and antagonised the slow-wave saturation peak in a dose-dependent fashion that appeared dependent on the postsynaptic membrane potential, providing a plausible mechanism for how metabolic changes can modulate anaesthesia sensitivity.

CONCLUSIONS

Our model suggests how metabolic processes can modulate anaesthesia and produce non-receptor dependent drug sensitivity.

Computational singular perturbation analysis of brain lactate metabolism

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.

A computational model integrating brain electrophysiology and metabolism highlights the key role of extracellular potassium and oxygen

The human brain is a small organ which uses a disproportionate amount of the total metabolic energy production in the body. While it is well understood that the most significant energy sink is the maintenance of the neuronal membrane potential during the brain signaling activity, the role of astrocytes in the energy balance continues to be the topic of a lot of research. A key function of astrocytes, besides clearing glutamate from the synaptic clefts, is the potassium clearing after neuronal activation. Extracellular potassium plays a significant role in triggering neuronal firing, and elevated concentration of potassium may lead to abnormal firing patterns, e.g., seizures, thus emphasizing the importance of the glial K⁺ buffering role. The predictive mathematical model proposed in this paper elucidates the role of glial potassium clearing in brain energy metabolism, integrating a detailed model of the ion dynamics which regulates neuronal firing with a four compartment metabolic model. Because of the very different characteristic time scales of electrophysiology and metabolism, care must be taken when coupling the two models to ensure that the predictions, e.g., neuronal firing frequencies and the oxygen-glucose index (OGI) of the brain during activation and rest, are in agreement with empirical observations. The temporal multi-scale nature of the problem requires the design of new computational tools to ensure a stable and accurate numerical treatment. The model predictions for different protocols, including combinations of elevated activation and ischemic episodes, are in good agreement with experimental observations reported in the literature.

QSP Toolbox: Computational Implementation of Integrated Workflow Components for Deploying Multi-Scale Mechanistic Models

Quantitative systems pharmacology (QSP) modeling has become increasingly important in pharmaceutical research and development, and is a powerful tool to gain mechanistic insights into the complex dynamics of biological systems in response to drug treatment. However, even once a suitable mathematical framework to describe the pathophysiology and mechanisms of interest is established, final model calibration and the exploration of variability can be challenging and time consuming. QSP models are often formulated as multi-scale, multi-compartment nonlinear systems of ordinary differential equations. Commonly accepted modeling strategies, workflows, and tools have promise to greatly improve the efficiency of QSP methods and improve productivity. In this paper, we present the QSP Toolbox, a set of functions, structure array conventions, and class definitions that computationally implement critical elements of QSP workflows including data integration, model calibration, and variability exploration. We present the application of the toolbox to an ordinary differential equations-based model for antibody drug conjugates. As opposed to a single stepwise reference model calibration, the toolbox also facilitates simultaneous parameter optimization and variation across multiple in vitro, in vivo, and clinical assays to more comprehensively generate alternate mechanistic hypotheses that are in quantitative agreement with available data. The toolbox also includes scripts for developing and applying virtual populations to mechanistic exploration of biomarkers and efficacy. We anticipate that the QSP Toolbox will be a useful resource that will facilitate implementation, evaluation, and sharing of new methodologies in a common framework that will greatly benefit the community.

A spatially distributed computational model of brain cellular metabolism

This paper develops a three-dimensional spatially distributed model of brain cellular metabolism and investigates how the locus of the synaptic activity in reference to the capillaries and diffusion affects the behavior of the model, a type of analysis which is impossible to carry out in spatially lumped models, which are shown to be consistent spatially averaged approximations of the distributed model. To avoid a geometrically detailed modeling of the complex structure of the tissue consisting of different cell types and the extracellular space, the distributed model is based on a novel multi-domain formulation of reaction-diffusion equations, accounting also for separate mitochondria. The model reduction relating the spatially distributed model and lower dimensional reduced models, including the well-mixed spatially lumped compartment model, is carefully explained. We illustrate the effects of losing the spatial resolution with a computed example which is based on a reduced one-dimensional distributed radial model, and look into how the model behaves when the locus of the synaptic activity in reference to the capillaries is changed. By averaging the fluxes and concentrations in the distributed radial model to correspond to quantities in a lumped model, and further by estimating the parameters in the lumped, we conclude that varying the locus of the synaptic activity in reference to the capillaries alters significantly the lumped model parameters. This observation seems to be consequential for parameter estimation procedures from data when the spatial resolution is insufficient to determine the locus of the activity, indicating that the model uncertainty is an inherent feature of lumped models.

Modelling biochemical features of mitochondrial neuropathology

BACKGROUND

The neuropathology of mitochondrial disease is well characterised. However, pathophysiological mechanisms at the level of biochemistry and cell biology are less clear. Progress in this area has been hampered by the limited accessibility of neurologically relevant material for analysis.

SCOPE OF REVIEW

Here we discuss the recent development of a variety of model systems that have greatly extended our capacity to understand the biochemical features associated with mitochondrial neuropathology. These include animal and cell based models, with mutations in both nuclear and mitochondrial DNA encoded genes, which aim to recapitulate the neuropathology and cellular biochemistry of mitochondrial diseases.

MAJOR CONCLUSIONS

Analysis of neurological tissue and cells from these models suggests that although there is no unifying mode of pathogenesis, dysfunction of the oxidative phosphorylation (OXPHOS) system is often central. This can be associated with altered reactive oxygen species (ROS) generation, disruption of the mitochondrial membrane potential (ΔΨm) and inadequate ATP synthesis. Thus, other cellular processes such as calcium (Ca(2+)) homeostasis, cellular signaling and mitochondrial morphology could be altered, ultimately compromising viability of neuronal cells.

GENERAL SIGNIFICANCE

Mechanisms of neuronal dysfunction in mitochondrial disease are only just beginning to be characterised, are system dependent and complex, and not merely driven by energy deficiency. The diversity of pathogenic mechanisms emphasises the need for characterisation in a wide range of models, as different therapeutic strategies are likely to be needed for different diseases. This article is part of a Special Issue entitled Frontiers of Mitochondrial Research.

Unraveling the complex metabolic nature of astrocytes

Since the initial description of astrocytes by neuroanatomists of the nineteenth century, a critical metabolic role for these cells has been suggested in the central nervous system. Nonetheless, it took several technological and conceptual advances over many years before we could start to understand how they fulfill such a role. One of the important and early recognized metabolic function of astrocytes concerns the reuptake and recycling of the neurotransmitter glutamate. But the description of this initial property will be followed by several others including an implication in the supply of energetic substrates to neurons. Indeed, despite the fact that like most eukaryotic non-proliferative cells, astrocytes rely on oxidative metabolism for energy production, they exhibit a prominent aerobic glycolysis capacity. Moreover, this unusual metabolic feature was found to be modulated by glutamatergic activity constituting the initial step of the neurometabolic coupling mechanism. Several approaches, including biochemical measurements in cultured cells, genetic screening, dynamic cell imaging, nuclear magnetic resonance spectroscopy and mathematical modeling, have provided further insights into the intrinsic characteristics giving rise to these key features of astrocytes. This review will provide an account of the different results obtained over several decades that contributed to unravel the complex metabolic nature of astrocytes that make this cell type unique.

Energy metabolism and glutamate-glutamine cycle in the brain: a stoichiometric modeling perspective

Background

The energetics of cerebral activity critically relies on the functional and metabolic interactions between neurons and astrocytes. Important open questions include the relation between neuronal versus astrocytic energy demand, glucose uptake and intercellular lactate transfer, as well as their dependence on the level of activity.

Results

We have developed a large-scale, constraint-based network model of the metabolic partnership between astrocytes and glutamatergic neurons that allows for a quantitative appraisal of the extent to which stoichiometry alone drives the energetics of the system. We find that the velocity of the glutamate-glutamine cycle (V_cyc) explains part of the uncoupling between glucose and oxygen utilization at increasing V_cyc levels. Thus, we are able to characterize different activation states in terms of the tissue oxygen-glucose index (OGI). Calculations show that glucose is taken up and metabolized according to cellular energy requirements, and that partitioning of the sugar between different cell types is not significantly affected by V_cyc. Furthermore, both the direction and magnitude of the lactate shuttle between neurons and astrocytes turn out to depend on the relative cell glucose uptake while being roughly independent of V_cyc.

Conclusions

These findings suggest that, in absence of ad hoc activity-related constraints on neuronal and astrocytic metabolism, the glutamate-glutamine cycle does not control the relative energy demand of neurons and astrocytes, and hence their glucose uptake and lactate exchange.

Quantitative in silico Analysis of Neurotransmitter Pathways Under Steady State Conditions

The modeling of glutamate/GABA-glutamine cycling in the brain tissue involving astrocytes, glutamatergic and GABAergic neurons leads to a complex compartmentalized metabolic network that comprises neurotransmitter synthesis, shuttling, and degradation. Without advanced computational tools, it is difficult to quantitatively track possible scenarios and identify viable ones. In this article, we follow a sampling-based computational paradigm to analyze the biochemical network in a multi-compartment system modeling astrocytes, glutamatergic, and GABAergic neurons, and address some questions about the details of transmitter cycling, with particular emphasis on the ammonia shuttling between astrocytes and neurons, and the synthesis of transmitter GABA. More specifically, we consider the joint action of the alanine-lactate shuttle, the branched chain amino acid shuttle, and the glutamine-glutamate cycle, as well as the role of glutamate dehydrogenase (GDH) activity. When imposing a minimal amount of bound constraints on reaction and transport fluxes, a preferred stoichiometric steady state equilibrium requires an unrealistically high reductive GDH activity in neurons, indicating the need for additional bound constants which were included in subsequent computer simulations. The statistical flux balance analysis also suggests a stoichiometrically viable role for leucine transport as an alternative to glutamine for replenishing the glutamate pool in neurons.