Spontaneous oscillations of dressed neurons: a new mechanism for epilepsy?

Astrocyte's self-repairing characteristics improve working memory in spiking neuronal networks

Astrocytes play a significant role in the working memory (WM) mechanism, yet their contribution to spiking neuron-astrocyte networks (SNAN) is underexplored. This study proposes a non-probabilistic SNAN incorporating a self-repairing (SR) mechanism through endocannabinoid pathways to facilitate WM function. Four experiments were conducted with different damaging patterns, replicating close-to-realistic synaptic impairments. Simulation results suggest that the SR process enhances WM performance by improving the consistency of neuronal firing. Moreover, the intercellular astrocytic [Ca]2+ transmission via gap junctions improves WM and SR processes. With increasing damage, WM and SR activities initially fail for non-matched samples and then for smaller and minimally overlapping matched samples. Simulation results also indicate that the inclusion of the SR mechanism in both random and continuous forms of damage improves the resilience of the WM by approximately 20%. This study highlights the importance of astrocytes in synaptically impaired networks.

Breathing cluster in complex neuron-astrocyte networks

Brain activities are featured by spatially distributed neural clusters of coherent firings and a spontaneous slow switching of the clusters between the coherent and incoherent states. Evidences from recent in vivo experiments suggest that astrocytes, a type of glial cell regarded previously as providing only structural and metabolic supports to neurons, participate actively in brain functions by regulating the neural firing activities, yet the underlying mechanism remains unknown. Here, introducing astrocyte as a reservoir of the glutamate released from the neuron synapses, we propose the model of the complex neuron-astrocyte network, and investigate the roles of astrocytes in regulating the cluster synchronization behaviors of networked chaotic neurons. It is found that a specific set of neurons on the network are synchronized and form a cluster, while the remaining neurons are kept as desynchronized. Moreover, during the course of network evolution, the cluster is switching between the synchrony and asynchrony states in an intermittent fashion, henceforth the phenomenon of "breathing cluster." By the method of symmetry-based analysis, we conduct a theoretical investigation on the synchronizability of the cluster. It is revealed that the contents of the cluster are determined by the network symmetry, while the breathing of the cluster is attributed to the interplay between the neural network and the astrocyte. The phenomenon of breathing cluster is demonstrated in different network models, including networks with different sizes, nodal dynamics, and coupling functions. The findings shed light on the cellular mechanism of astrocytes in regulating neural activities and give insights into the state-switching of the neocortex.

Analysis of Network Models with Neuron-Astrocyte Interactions

Neural networks, composed of many neurons and governed by complex interactions between them, are a widely accepted formalism for modeling and exploring global dynamics and emergent properties in brain systems. In the past decades, experimental evidence of computationally relevant neuron-astrocyte interactions, as well as the astrocytic modulation of global neural dynamics, have accumulated. These findings motivated advances in computational glioscience and inspired several models integrating mechanisms of neuron-astrocyte interactions into the standard neural network formalism. These models were developed to study, for example, synchronization, information transfer, synaptic plasticity, and hyperexcitability, as well as classification tasks and hardware implementations. We here focus on network models of at least two neurons interacting bidirectionally with at least two astrocytes that include explicitly modeled astrocytic calcium dynamics. In this study, we analyze the evolution of these models and the biophysical, biochemical, cellular, and network mechanisms used to construct them. Based on our analysis, we propose how to systematically describe and categorize interaction schemes between cells in neuron-astrocyte networks. We additionally study the models in view of the existing experimental data and present future perspectives. Our analysis is an important first step towards understanding astrocytic contribution to brain functions. However, more advances are needed to collect comprehensive data about astrocyte morphology and physiology in vivo and to better integrate them in data-driven computational models. Broadening the discussion about theoretical approaches and expanding the computational tools is necessary to better understand astrocytes' roles in brain functions.

Interplay between biochemical processes and network properties generates neuronal up and down states at the tripartite synapse

Neuronal up and down states have long been known to exist both in vitro and in vivo. A variety of functions and mechanisms have been proposed for their generation, but there has not been a clear connection between the functions and mechanisms. We explore the potential contribution of cellular-level biochemistry to the network-level mechanisms thought to underlie the generation of up and down states. We develop a neurochemical model of a single tripartite synapse, assumed to be within a network of similar tripartite synapses, to investigate possible function-mechanism links for the appearance of up and down states. We characterize the behavior of our model in different regions of parameter space and show that resource limitation at the tripartite synapse affects its ability to faithfully transmit input signals, leading to extinction-down states. Recovery of resources allows for "reignition" into up states. The tripartite synapse exhibits distinctive "regimes" of operation depending on whether ATP, neurotransmitter (glutamate), both, or neither, is limiting. Our model qualitatively matches the behavior of six disparate experimental systems, including both in vitro and in vivo models, without changing any model parameters except those related to the experimental conditions. We also explore the effects of varying different critical parameters within the model. Here we show that availability of energy, represented by ATP, and glutamate for neurotransmission at the cellular level are intimately related, and are capable of promoting state transitions at the network level as ignition and extinction phenomena. Our model is complementary to existing models of neuronal up and down states in that it focuses on cellular-level dynamics while still retaining essential network-level processes. Our model predicts the existence of a "final common pathway" of behavior at the tripartite synapse arising from scarcity of resources and may explain use dependence in the phenomenon of "local sleep." Ultimately, sleeplike behavior may be a fundamental property of networks of tripartite synapses.

Double resonance induced by group coupling with quenched disorder

Results show that the astrocytes can not only listen to the talk of large assemble of neurons but also give advice to the conversations and are significant sources of heterogeneous couplings as well. In the present work, we focus on such regulation character of astrocytes and explore the role of heterogeneous couplings among interacted neuron-astrocyte components in a signal response. We consider reduced dynamics in which the listening and advising processes of astrocytes are mapped into the form of group coupling, where the couplings are normally distributed. In both globally coupled overdamped bistable oscillators and an excitable FitzHugh-Nagumo (FHN) neuron model, we numerically and analytically demonstrate that two types of bell-shaped collective response curves can be obtained as the ensemble coupling strength or the heterogeneity of group coupling rise, respectively, which can be seen as a new type of double resonance. Furthermore, through the bifurcation analysis, we verify that these resonant signal responses stem from the competition between dispersion and aggregation induced by heterogeneous group and positive pairwise couplings, respectively. Our results contribute to a better understanding of the signal propagation in coupled systems with quenched disorder.

Transition behavior of the seizure dynamics modulated by the astrocyte inositol triphosphate noise

Epilepsy is a neurological disorder with recurrent seizures, which convey complex dynamical characteristics including chaos and randomness. Until now, the underlying mechanism has not been fully elucidated, especially the bistable property beneath the epileptic random induction phenomena in certain conditions. Inspired by the recent finding that astrocyte GTPase-activating protein (G-protein)-coupled receptors could be involved in stochastic epileptic seizures, we proposed a neuron-astrocyte network model, incorporating the noise of the astrocytic second messenger, inositol triphosphate (IP3) that is modulated by G-protein-coupled receptor activation. Based on this model, we have statistically analyzed the transitions of epileptic seizures by performing repeatable simulation trials. Our simulation results show that the increase in the IP3 noise intensity induces depolarization-block epileptic seizures together with an increase in neuronal firing frequency, consistent with corresponding experiments. Meanwhile, the bistable states of the seizure dynamics were present under certain noise intensities, during which the neuronal firing pattern switches between regular sparse spiking and epileptic seizure states. This random presence of epileptic seizures is absent when the noise intensity continues to increase, accompanying with an increase in the epileptic depolarization block duration. The simulation results also shed light on the fact that calcium signals in astrocytes play significant roles in the pattern formations of the epileptic seizure. Our results provide a potential pathway for understanding the epileptic randomness in certain conditions.

Modeling Active Non-Markovian Oscillations

Modeling noisy oscillations of active systems is one of the current challenges in physics and biology. Because the physical mechanisms of such processes are often difficult to identify, we propose a linear stochastic model driven by a non-Markovian bistable noise that is capable of generating self-sustained periodic oscillation. We derive analytical predictions for most relevant dynamical and thermodynamic properties of the model. This minimal model turns out to describe accurately bistablelike oscillatory motion of hair bundles in bullfrog sacculus, extracted from experimental data. Based on and in agreement with these data, we estimate the power required to sustain such active oscillations to be of the order of 100  k_{B}T per oscillation cycle.

A New Real-Time Analog Circuit of Ca2+ Li-Rinzel Astrocyte Model Based on Analytical Method

Different biological models are used to study physical behaviors in neural networks. So far, various models of neural network components such as neurons, synapses, and astrocytes have been proposed. An astrocyte is one of the crucial parts introduced in multiple models. A model of astrocytes used as a good reference in various papers is the Li-Rinzel calcium model. This paper presents a real-time analog circuit of the Li-Rinzel calcium model based on common-mode (CM) in 180nm CMOS technology. To the best of our knowledge, this work is the first report to introduce a real-time analog Li-Rinzel model. The careful design of equations and low power consumption are essential features of this circuit. The real-time behavior is also crucial compared with the accelerated time circuits presented so far.

Controlling the spontaneous firing behavior of a neuron with astrocyte

Mounting evidence in recent years suggests that astrocytes, a sub-type of glial cells, not only serve metabolic and structural support for neurons and synapses but also play critical roles in the regulation of proper functioning of the nervous system. In this work, we investigate the effect of astrocytes on the spontaneous firing activity of a neuron through a combined model that includes a neuron-astrocyte pair. First, we show that an astrocyte may provide a kind of multistability in neuron dynamics by inducing different firing modes such as random and bursty spiking. Then, we identify the underlying mechanism of this behavior and search for the astrocytic factors that may have regulatory roles in different firing regimes. More specifically, we explore how an astrocyte can participate in the occurrence and control of spontaneous irregular spiking activity of a neuron in random spiking mode. Additionally, we systematically investigate the bursty firing regime dynamics of the neuron under the variation of biophysical facts related to the intracellular environment of the astrocyte. It is found that an astrocyte coupled to a neuron can provide a control mechanism for both spontaneous firing irregularity and burst firing statistics, i.e., burst regularity and size.

Neuron-Glia Interactions and Brain Circuits

Recent evidence suggests that glial cells take an active role in a number of brain functions that were previously attributed solely to neurons. For example, astrocytes, one type of glial cells, have been shown to promote coordinated activation of neuronal networks, modulate sensory-evoked neuronal network activity, and influence brain state transitions during development. This reinforces the idea that astrocytes not only provide the "housekeeping" for the neurons, but that they also play a vital role in supporting and expanding the functions of brain circuits and networks. Despite this accumulated knowledge, the field of computational neuroscience has mostly focused on modeling neuronal functions, ignoring the glial cells and the interactions they have with the neurons. In this chapter, we introduce the biology of neuron-glia interactions, summarize the existing computational models and tools, and emphasize the glial properties that may be important in modeling brain functions in the future.

Astrocytic modulation on neuronal electric mode selection induced by magnetic field effect

Astrocytes as well as electromagnetic induction have been primarily considered as main factors in regulating neuronal firing patterns in the recent decade. In this work, an improved neuron-astrocyte model in consideration of the modulation of astrocytes and the electromagnetic induction is employed to explore the extend to which both of the factors affect the firing modes of the neurons. The "alternation mode", defined as the alternative of neural normal spiking mode with the high-frequency bursting-like mode, clearly shows the functions of astrocytes on neurons. Moreover, the firing pattern of the neuron becomes more abnormal when astrocytes are hyper-excitable, the reason why the abnormal coupling of the astrocyte leads to the "alternation mode" of the neuron have been studied. In addition, the effect of electromagnetic induction manifests nonlinear characteristic towards neurons, complex firing modes of neurons are observed in the weaker field and a switching mode consists with quiescent and spiking mode appears when there is a higher stronger field. This approved model can reveal the normal or abnormal electric activities of neuron considered electromagnetic induction induced by the degree of excitability of the astrocyte. These results can provide potential understanding about the effects of astrocyte on neuronal activity when the coupling of electromagnetic field is considered.

Vibrational resonance in a neuron-astrocyte coupled model

Recent findings have revealed that not only neurons but also astrocytes, a special type of glial cells, are major players of neuronal information processing. It is now widely accepted that they contribute to the regulation of their microenvironment by cross-talking with neurons via gliotransmitters. In this context, we here study the phenomenon of vibrational resonance in neurons by considering their interaction with astrocytes. Our analysis of a neuron-astrocyte pair reveals that intracellular dynamics of astrocytes can induce a double vibrational resonance effect in the weak signal detection performance of a neuron, exhibiting two distinct wells centred at different high-frequency driving amplitudes. We also identify the underlying mechanism of this behaviour, showing that the interaction of widely separated time scales of neurons, astrocytes and driving signals is the key factor for the emergence and control of double vibrational resonance. This article is part of the theme issue 'Vibrational and stochastic resonance in driven nonlinear systems (part 2)'.

Modeling Working Memory in a Spiking Neuron Network Accompanied by Astrocytes

We propose a novel biologically plausible computational model of working memory (WM) implemented by a spiking neuron network (SNN) interacting with a network of astrocytes. The SNN is modeled by synaptically coupled Izhikevich neurons with a non-specific architecture connection topology. Astrocytes generating calcium signals are connected by local gap junction diffusive couplings and interact with neurons via chemicals diffused in the extracellular space. Calcium elevations occur in response to the increased concentration of the neurotransmitter released by spiking neurons when a group of them fire coherently. In turn, gliotransmitters are released by activated astrocytes modulating the strength of the synaptic connections in the corresponding neuronal group. Input information is encoded as two-dimensional patterns of short applied current pulses stimulating neurons. The output is taken from frequencies of transient discharges of corresponding neurons. We show how a set of information patterns with quite significant overlapping areas can be uploaded into the neuron-astrocyte network and stored for several seconds. Information retrieval is organized by the application of a cue pattern representing one from the memory set distorted by noise. We found that successful retrieval with the level of the correlation between the recalled pattern and ideal pattern exceeding 90% is possible for the multi-item WM task. Having analyzed the dynamical mechanism of WM formation, we discovered that astrocytes operating at a time scale of a dozen of seconds can successfully store traces of neuronal activations corresponding to information patterns. In the retrieval stage, the astrocytic network selectively modulates synaptic connections in the SNN leading to successful recall. Information and dynamical characteristics of the proposed WM model agrees with classical concepts and other WM models.

Modulation of Astrocytes on Mode Selection of Neuron Firing Driven by Electromagnetic Induction

Both of astrocytes and electromagnetic induction are magnificent to modulate neuron firing by introducing feedback currents to membrane potential. An improved astro-neuron model considering both of the two factors is employed to investigate their different roles in modulation. The mixing mode, defined by combination of period bursting and depolarization blockage, characterizes the effect of astrocytes. Mixing mode and period bursting alternatively appear in parameter space with respect to the amplitude of feedback current on neuron from astrocyte modulation. However, magnetic flux obviously plays a role of neuron firing inhibition. It not only repels the mixing mode but also suppresses period bursting. The mixing mode becomes period bursting mode and even resting state when astrocytes are hyperexcitable. Abnormal activities of astrocytes are capable to induce depolarization blockage to compose the mixing mode together with bursting mode. But electromagnetic induction shows its strong ability of inhibition of neuron firing, which is also illustrated in the bifurcation diagram. Indeed, the combination of the two factors and appropriate choice of parameters show the great potential to control disorder of neuron firing like epilepsy.

Dynamic Transitions in Neuronal Network Firing Sustained by Abnormal Astrocyte Feedback

Astrocytes play a crucial role in neuronal firing activity. Their abnormal state may lead to the pathological transition of neuronal firing patterns and even induce seizures. However, there is still little evidence explaining how the astrocyte network modulates seizures caused by structural abnormalities, such as gliosis. To explore the role of gliosis of the astrocyte network in epileptic seizures, we first established a direct astrocyte feedback neuronal network model on the basis of the hippocampal CA3 neuron-astrocyte model to simulate the condition of gliosis when astrocyte processes swell and the feedback to neurons increases in an abnormal state. We analyzed the firing pattern transitions of the neuronal network when astrocyte feedback starts to change via increases in both astrocyte feedback intensity and the connection probability of astrocytes to neurons in the network. The results show that as the connection probability and astrocyte feedback intensity increase, neuronal firing transforms from a nonepileptic synchronous firing state to an asynchronous firing state, and when astrocyte feedback starts to become abnormal, seizure-like firing becomes more severe and synchronized; meanwhile, the synchronization area continues to expand and eventually transforms into long-term seizure-like synchronous firing. Therefore, our results prove that astrocyte feedback can regulate the firing of the neuronal network, and when the astrocyte network develops gliosis, there will be an increase in the induction rate of epileptic seizures.

Astrocyte-mediated spike-timing-dependent long-term depression modulates synaptic properties in the developing cortex

Astrocytes have been shown to modulate synaptic transmission and plasticity in specific cortical synapses, but our understanding of the underlying molecular and cellular mechanisms remains limited. Here we present a new biophysicochemical model of a somatosensory cortical layer 4 to layer 2/3 synapse to study the role of astrocytes in spike-timing-dependent long-term depression (t-LTD) in vivo. By applying the synapse model and electrophysiological data recorded from rodent somatosensory cortex, we show that a signal from a postsynaptic neuron, orchestrated by endocannabinoids, astrocytic calcium signaling, and presynaptic N-methyl-D-aspartate receptors coupled with calcineurin signaling, induces t-LTD which is sensitive to the temporal difference between post- and presynaptic firing. We predict for the first time the dynamics of astrocyte-mediated molecular mechanisms underlying t-LTD and link complex biochemical networks at presynaptic, postsynaptic, and astrocytic sites to the time window of t-LTD induction. During t-LTD a single astrocyte acts as a delay factor for fast neuronal activity and integrates fast neuronal sensory processing with slow non-neuronal processing to modulate synaptic properties in the brain. Our results suggest that astrocytes play a critical role in synaptic computation during postnatal development and are of paramount importance in guiding the development of brain circuit functions, learning and memory.

A Computational Model to Investigate GABA-Activated Astrocyte Modulation of Neuronal Excitation

Gamma-aminobutyric acid (GABA) is critical for proper neural network function and can activate astrocytes to induce neuronal excitability; however, the mechanism by which astrocytes transform inhibitory signaling to excitatory enhancement remains unclear. Computational modeling can be a powerful tool to provide further understanding of how GABA-activated astrocytes modulate neuronal excitation. In the present study, we implemented a biophysical neuronal network model to investigate the effects of astrocytes on excitatory pre- and postsynaptic terminals following exposure to increasing concentrations of external GABA. The model completely describes the effects of GABA on astrocytes and excitatory presynaptic terminals within the framework of glutamatergic gliotransmission according to neurophysiological findings. Utilizing this model, our results show that astrocytes can rapidly respond to incoming GABA by inducing Ca²⁺ oscillations and subsequent gliotransmitter glutamate release. Elevation in GABA concentrations not only naturally decreases neuronal spikes but also enhances astrocytic glutamate release, which leads to an increase in astrocyte-mediated presynaptic release and postsynaptic slow inward currents. Neuronal excitation induced by GABA-activated astrocytes partly counteracts the inhibitory effect of GABA. Overall, the model helps to increase knowledge regarding the involvement of astrocytes in neuronal regulation using simulated bath perfusion of GABA, which may be useful for exploring the effects of GABA-type antiepileptic drugs.

The Human Body as a Super Network: Digital Methods to Analyze the Propagation of Aging

Biological aging is a complex process involving multiple biological processes. These can be understood theoretically though considering them as individual networks-e.g., epigenetic networks, cell-cell networks (such as astroglial networks), and population genetics. Mathematical modeling allows the combination of such networks so that they may be studied in unison, to better understand how the so-called "seven pillars of aging" combine and to generate hypothesis for treating aging as a condition at relatively early biological ages. In this review, we consider how recent progression in mathematical modeling can be utilized to investigate aging, particularly in, but not exclusive to, the context of degenerative neuronal disease. We also consider how the latest techniques for generating biomarker models for disease prediction, such as longitudinal analysis and parenclitic analysis can be applied to as both biomarker platforms for aging, as well as to better understand the inescapable condition. This review is written by a highly diverse and multi-disciplinary team of scientists from across the globe and calls for greater collaboration between diverse fields of research.

Astrocyte-induced positive integrated information in neuron-astrocyte ensembles

Integrated information is a quantitative measure from information theory of how tightly all parts of a system are interconnected in terms of information exchange. In this study we show that astrocytes, playing an important role in regulation of information transmission between neurons, may contribute to a generation of positive integrated information in neuronal ensembles. Analytically and numerically we show that the presence of astrocytic regulation of neurotransmission may be essential for this information attribute in neuroastrocytic ensembles. Moreover, the proposed "spiking-bursting" mechanism of generating positive integrated information is shown to be generic and not limited to neuron-astrocyte networks and is given a complete analytic description.

Challenges in Reproducibility, Replicability, and Comparability of Computational Models and Tools for Neuronal and Glial Networks, Cells, and Subcellular Structures

The possibility to replicate and reproduce published research results is one of the biggest challenges in all areas of science. In computational neuroscience, there are thousands of models available. However, it is rarely possible to reimplement the models based on the information in the original publication, let alone rerun the models just because the model implementations have not been made publicly available. We evaluate and discuss the comparability of a versatile choice of simulation tools: tools for biochemical reactions and spiking neuronal networks, and relatively new tools for growth in cell cultures. The replicability and reproducibility issues are considered for computational models that are equally diverse, including the models for intracellular signal transduction of neurons and glial cells, in addition to single glial cells, neuron-glia interactions, and selected examples of spiking neuronal networks. We also address the comparability of the simulation results with one another to comprehend if the studied models can be used to answer similar research questions. In addition to presenting the challenges in reproducibility and replicability of published results in computational neuroscience, we highlight the need for developing recommendations and good practices for publishing simulation tools and computational models. Model validation and flexible model description must be an integral part of the tool used to simulate and develop computational models. Constant improvement on experimental techniques and recording protocols leads to increasing knowledge about the biophysical mechanisms in neural systems. This poses new challenges for computational neuroscience: extended or completely new computational methods and models may be required. Careful evaluation and categorization of the existing models and tools provide a foundation for these future needs, for constructing multiscale models or extending the models to incorporate additional or more detailed biophysical mechanisms. Improving the quality of publications in computational neuroscience, enabling progressive building of advanced computational models and tools, can be achieved only through adopting publishing standards which underline replicability and reproducibility of research results.

Computational Models for Calcium-Mediated Astrocyte Functions

The computational neuroscience field has heavily concentrated on the modeling of neuronal functions, largely ignoring other brain cells, including one type of glial cell, the astrocytes. Despite the short history of modeling astrocytic functions, we were delighted about the hundreds of models developed so far to study the role of astrocytes, most often in calcium dynamics, synchronization, information transfer, and plasticity in vitro, but also in vascular events, hyperexcitability, and homeostasis. Our goal here is to present the state-of-the-art in computational modeling of astrocytes in order to facilitate better understanding of the functions and dynamics of astrocytes in the brain. Due to the large number of models, we concentrated on a hundred models that include biophysical descriptions for calcium signaling and dynamics in astrocytes. We categorized the models into four groups: single astrocyte models, astrocyte network models, neuron-astrocyte synapse models, and neuron-astrocyte network models to ease their use in future modeling projects. We characterized the models based on which earlier models were used for building the models and which type of biological entities were described in the astrocyte models. Features of the models were compared and contrasted so that similarities and differences were more readily apparent. We discovered that most of the models were basically generated from a small set of previously published models with small variations. However, neither citations to all the previous models with similar core structure nor explanations of what was built on top of the previous models were provided, which made it possible, in some cases, to have the same models published several times without an explicit intention to make new predictions about the roles of astrocytes in brain functions. Furthermore, only a few of the models are available online which makes it difficult to reproduce the simulation results and further develop the models. Thus, we would like to emphasize that only via reproducible research are we able to build better computational models for astrocytes, which truly advance science. Our study is the first to characterize in detail the biophysical and biochemical mechanisms that have been modeled for astrocytes.

The role of neuron-glia interactions in the emergence of ultra-slow oscillations

Ultra-slow cortical oscillatory activity of 1-100 mHz has been recorded in human by electroencephalography and in dissociated cultures of cortical rat neurons, but the underlying mechanisms remain to be elucidated. This study presents a computational model of ultra-slow oscillatory activity based on the interaction between neurons and astrocytes. We predict that the frequency of these oscillations closely depends on activation of astrocytes in the network, which is reflected by oscillations of their intracellular calcium concentrations with periods between tens of seconds and minutes. An increase of intracellular calcium in astrocytes triggers the release of adenosine triphosphate from these cells which may alter transmission at nearby synapses by increasing or decreasing neurotransmitter release. These results provide theoretical support for the emerging awareness of astrocytes as active players in the regulation of neural activity and identify neuron-astrocyte interactions as a potential primary mechanism for the emergence of ultra-slow cortical oscillations.

Reproducibility and Comparability of Computational Models for Astrocyte Calcium Excitability

The scientific community across all disciplines faces the same challenges of ensuring accessibility, reproducibility, and efficient comparability of scientific results. Computational neuroscience is a rapidly developing field, where reproducibility and comparability of research results have gained increasing interest over the past years. As the number of computational models of brain functions is increasing, we chose to address reproducibility using four previously published computational models of astrocyte excitability as an example. Although not conventionally taken into account when modeling neuronal systems, astrocytes have been shown to take part in a variety of in vitro and in vivo phenomena including synaptic transmission. Two of the selected astrocyte models describe spontaneous calcium excitability, and the other two neurotransmitter-evoked calcium excitability. We specifically addressed how well the original simulation results can be reproduced with a reimplementation of the models. Additionally, we studied how well the selected models can be reused and whether they are comparable in other stimulation conditions and research settings. Unexpectedly, we found out that three of the model publications did not give all the necessary information required to reimplement the models. In addition, we were able to reproduce the original results of only one of the models completely based on the information given in the original publications and in the errata. We actually found errors in the equations provided by two of the model publications; after modifying the equations accordingly, the original results were reproduced more accurately. Even though the selected models were developed to describe the same biological event, namely astrocyte calcium excitability, the models behaved quite differently compared to one another. Our findings on a specific set of published astrocyte models stress the importance of proper validation of the models against experimental wet-lab data from astrocytes as well as the careful review process of models. A variety of aspects of model development could be improved, including the presentation of models in publications and databases. Specifically, all necessary mathematical equations, as well as parameter values, initial values of variables, and stimuli used should be given precisely for successful reproduction of scientific results.

From in silico astrocyte cell models to neuron-astrocyte network models: A review

The idea that astrocytes may be active partners in synaptic information processing has recently emerged from abundant experimental reports. Because of their spatial proximity to neurons and their bidirectional communication with them, astrocytes are now considered as an important third element of the synapse. Astrocytes integrate and process synaptic information and by doing so generate cytosolic calcium signals that are believed to reflect neuronal transmitter release. Moreover, they regulate neuronal information transmission by releasing gliotransmitters into the synaptic cleft affecting both pre- and postsynaptic receptors. Concurrent with the first experimental reports of the astrocytic impact on neural network dynamics, computational models describing astrocytic functions have been developed. In this review, we give an overview over the published computational models of astrocytic functions, from single-cell dynamics to the tripartite synapse level and network models of astrocytes and neurons.

Dynamic transition of neuronal firing induced by abnormal astrocytic glutamate oscillation

The gliotransmitter glutamate released from astrocytes can modulate neuronal firing by activating neuronal N-methyl-D-aspartic acid (NMDA) receptors. This enables astrocytic glutamate(AG) to be involved in neuronal physiological and pathological functions. Based on empirical results and classical neuron-glial "tripartite synapse" model, we propose a practical model to describe extracellular AG oscillation, in which the fluctuation of AG depends on the threshold of calcium concentration, and the effect of AG degradation is considered as well. We predict the seizure-like discharges under the dysfunction of AG degradation duration. Consistent with our prediction, the suppression of AG uptake by astrocytic transporters, which operates by modulating the AG degradation process, can account for the emergence of epilepsy.

Astrocytes as new targets to improve cognitive functions

Astrocytes are now viewed as key elements of brain wiring as well as neuronal communication. Indeed, they not only bridge the gap between metabolic supplies by blood vessels and neurons, but also allow fine control of neurotransmission by providing appropriate signaling molecules and insulation through a tight enwrapping of synapses. Recognition that astroglia is essential to neuronal communication is nevertheless fairly recent and the large body of evidence dissecting such role has focused on the synaptic level by identifying neuro- and gliotransmitters uptaken and released at synaptic or extrasynaptic sites. Yet, more integrated research deciphering the impact of astroglial functions on neuronal network activity have led to the reasonable assumption that the role of astrocytes in supervising synaptic activity translates in influencing neuronal processing and cognitive functions. Several investigations using recent genetic tools now support this notion by showing that inactivating or boosting astroglial function directly affects cognitive abilities. Accordingly, brain diseases resulting in impaired cognitive functions have seen their physiopathological mechanisms revisited in light of this primary protagonist of brain processing. We here provide a review of the current knowledge on the role of astrocytes in cognition and in several brain diseases including neurodegenerative disorders, psychiatric illnesses, as well as other conditions such as epilepsy. Potential astroglial therapeutic targets are also discussed.

Information transmission in a neuron-astrocyte coupled model

A coupled model containing two neurons and one astrocyte is constructed by integrating Hodgkin-Huxley neuronal model and Li-Rinzel calcium model. Based on this hybrid model, information transmission between neurons is studied numerically. Our results show that when the successive spikes are produced in neuron 1 (N1), the bursting-like spikes (BLSs) occur in two neurons simultaneously during the spikes being transferred to neuron 2 (N2). The existence of the astrocyte and a higher expression level of mGluRs facilitate the occurrence of BLSs, but the rate of occurrence is not sensitive to the parameters. Furthermore, time delay τ occurs during the information transmission, and τ is almost independent of the effect of the astrocyte. Additionally, we found that low coupling strength may result in the distortion of the information, and this distortion is also proven to be almost independent of the astrocyte.

A possible role of astrocytes in contextual memory retrieval: An analysis obtained using a quantitative framework

The hippocampus is central to our understanding of memory formation and retrieval. Its CA1 region is known for encoding contextual memory. Here, using a computational approach, which embeds existing physiological data, we propose a particular role of astrocytes in contextual memory retrieval. We provide a quantitative framework under which the astrocyte modulates the firing of a context-associated CA1 pyramidal neurons, resulting in a prominent tuning of neurons to a delta rhythm. Using the very framework, we further studied astrocytic function in the modulation of neuronal firing under pathological conditions, i.e., during astrocytic induction of epileptiform discharge in CA1 pyramidal neurons. Thus, we provide a quantitative framework that would aid understanding of the Schaffer collateral-CA1 tripartite synapse in health and disease.

Divide and conquer: functional segregation of synaptic inputs by astrocytic microdomains could alleviate paroxysmal activity following brain trauma

Traumatic brain injury often leads to epileptic seizures. Among other factors, homeostatic synaptic plasticity (HSP) mediates posttraumatic epileptogenesis through unbalanced synaptic scaling, partially compensating for the trauma-incurred loss of neural excitability. HSP is mediated in part by tumor necrosis factor alpha (TNFα), which is released locally from reactive astrocytes early after trauma in response to chronic neuronal inactivity. During this early period, TNFα is likely to be constrained to its glial sources; however, the contribution of glia-mediated spatially localized HSP to post-traumatic epileptogenesis remains poorly understood. We used computational model to investigate the reorganization of collective neural activity early after trauma. Trauma and synaptic scaling transformed asynchronous spiking into paroxysmal discharges. The rate of paroxysms could be reduced by functional segregation of synaptic input into astrocytic microdomains. Thus, we propose that trauma-triggered reactive gliosis could exert both beneficial and deleterious effects on neural activity.

A computational model of neuro-glio-vascular loop interactions

We present a computational, biophysical model of neuron-astrocyte-vessel interaction. Unlike other cells, neurons convey “hunger” signals to the vascular network via an intervening layer of glial cells (astrocytes); vessels dilate and release glucose which fuels neuronal firing. Existing computational models focus on only parts of this loop (neuron→astrocyte→vessel→neuron), whereas the proposed model describes the entire loop. Neuronal firing causes release of a neurotransmitter like glutamate which triggers release of vasodilator by astrocytes via a cascade of biochemical events. Vasodilators released from astrocytic endfeet cause blood vessels to dilate and release glucose into the interstitium, part of which is taken up by the astrocyticendfeet. Glucose is converted into lactate in the astrocyte and transported into the neuron. Glucose from the interstitium and lactate (produced from glucose) influx from astrocyte are converted into ATP in the neuron. Neuronal ATP is used to drive the Na⁺/K⁺ATPase pumps, which maintain ionic gradients necessary for neuronal firing. When placed in the metabolic loop, the neuron exhibits sustained firing only when the stimulation current is more than a minimum threshold. For various combinations of initial neuronal [ATP] and external current, the neuron exhibits a variety of firing patterns including sustained firing, firing after an initial pause, burst firing etc. Neurovascular interactions under conditions of constricted vessels are also studied. Most models of cerebral circulation describe neurovascular interactions exclusively in the “forward” neuron→vessel direction. The proposed model indicates possibility of “reverse” influence also, with vasomotion rhythms influencing neural firing patterns. Another idea that emerges out of the proposed work is that brain's computations may be more comprehensively understood in terms of neuro-glial-vascular dynamics and not in terms of neural firing alone.

Computational model of neuron-astrocyte interactions during focal seizure generation

Empirical research in the last decade revealed that astrocytes can respond to neurotransmitters with Ca²⁺ elevations and generate feedback signals to neurons which modulate synaptic transmission and neuronal excitability. This discovery changed our basic understanding of brain function and provided new perspectives for how astrocytes can participate not only to information processing, but also to the genesis of brain disorders, such as epilepsy. Epilepsy is a neurological disorder characterized by recurrent seizures that can arise focally at restricted areas and propagate throughout the brain. Studies in brain slice models suggest that astrocytes contribute to epileptiform activity by increasing neuronal excitability through a Ca²⁺-dependent release of glutamate. The underlying mechanism remains, however, unclear. In this study, we implemented a parsimonious network model of neurons and astrocytes. The model consists of excitatory and inhibitory neurons described by Izhikevich's neuron dynamics. The experimentally observed Ca²⁺ change in astrocytes in response to neuronal activity was modeled with linear equations. We considered that glutamate is released from astrocytes above certain intracellular Ca²⁺ concentrations thus providing a non-linear positive feedback signal to neurons. Propagating seizure-like ictal discharges (IDs) were reliably evoked in our computational model by repeatedly exciting a small area of the network, which replicates experimental results in a slice model of focal ID in entorhinal cortex. We found that the threshold of focal ID generation was lowered when an excitatory feedback-loop between astrocytes and neurons was included. Simulations show that astrocytes can contribute to ID generation by directly affecting the excitatory/inhibitory balance of the neuronal network. Our model can be used to obtain mechanistic insights into the distinct contributions of the different signaling pathways to the generation and propagation of focal IDs.

Computational models of neuron-astrocyte interaction in epilepsy

Astrocytes actively shape the dynamics of neurons and neuronal ensembles by affecting several aspects critical to neuronal function, such as regulating synaptic plasticity, modulating neuronal excitability, and maintaining extracellular ion balance. These pathways for astrocyte-neuron interaction can also enhance the information-processing capabilities of brains, but in other circumstances may lead the brain on the road to pathological ruin. In this article, we review the existing computational models of astrocytic involvement in epileptogenesis, focusing on their relevance to existing physiological data.

A mathematical model of the tripartite synapse: astrocyte-induced synaptic plasticity

In this paper, we present a biologically detailed mathematical model of tripartite synapses, where astrocytes modulate short-term synaptic plasticity. The model consists of a pre-synaptic bouton, a post-synaptic dendritic spine-head, a synaptic cleft and a peri-synaptic astrocyte controlling Ca(2 + ) dynamics inside the synaptic bouton. This in turn controls glutamate release dynamics in the cleft. As a consequence of this, glutamate concentration in the cleft has been modeled, in which glutamate reuptake by astrocytes has also been incorporated. Finally, dendritic spine-head dynamics has been modeled. As an application, this model clearly shows synaptic potentiation in the hippocampal region, i.e., astrocyte Ca(2 + ) mediates synaptic plasticity, which is in conformity with the majority of the recent findings (Perea and Araque (Science 317, 1083-1086, 2007); Henneberger et al. (Nature 463, 232-236, 2010); Navarrete et al. (PLoS Biol. 10, e1001259, 2012)).

A mathematical model for astrocytes mediated LTP at single hippocampal synapses

Many contemporary studies have shown that astrocytes play a significant role in modulating both short and long form of synaptic plasticity. There are very few experimental models which elucidate the role of astrocyte over Long-term Potentiation (LTP). Recently, Perea and Araque (Science 317:1083-1086, 2007) demonstrated a role of astrocytes in induction of LTP at single hippocampal synapses. They suggested a purely pre-synaptic basis for induction of this N-methyl-D-Aspartate (NMDA) Receptor-independent LTP. Also, the mechanisms underlying this pre-synaptic induction were not investigated. Here, in this article, we propose a mathematical model for astrocyte modulated LTP which successfully imitates the experimental findings of Perea and Araque (Science 317:1083-1086, 2007). Our study suggests the role of retrograde messengers, possibly Nitric Oxide (NO), for this pre-synaptically modulated LTP.

Bifurcation mechanisms of regular and chaotic network signaling in brain astrocytes

Bifurcation mechanisms underlying calcium oscillations in the network of astrocytes are investigated. Network model includes the dynamics of intracellular calcium concentration and intercellular diffusion of inositol 1,4,5-trisphosphate through gap junctions. Bifurcation analysis of underlying nonlinear dynamical system is presented. Parameter regions and principle bifurcation boundaries have been delineated and described. We show how variations of the diffusion rate can lead to generation of network calcium oscillations in originally nonoscillating cells. Different scenarios of regular activity and its transitions to chaotic dynamics have been obtained. Then, the bifurcations have been associated with statistical characteristics of calcium signals showing that different bifurcation scenarios yield qualitative changes in experimentally measurable quantities of the astrocyte activity, e.g., statistics of calcium spikes.

A neuron-astrocyte transistor-like model for neuromorphic dressed neurons

Experimental evidences on the role of the synaptic glia as an active partner together with the bold synapse in neuronal signaling and dynamics of neural tissue strongly suggest to investigate on a more realistic neuron-glia model for better understanding human brain processing. Among the glial cells, the astrocytes play a crucial role in the tripartite synapsis, i.e. the dressed neuron. A well-known two-way astrocyte-neuron interaction can be found in the literature, completely revising the purely supportive role for the glia. The aim of this study is to provide a computationally efficient model for neuron-glia interaction. The neuron-glia interactions were simulated by implementing the Li-Rinzel model for an astrocyte and the Izhikevich model for a neuron. Assuming the dressed neuron dynamics similar to the nonlinear input-output characteristics of a bipolar junction transistor, we derived our computationally efficient model. This model may represent the fundamental computational unit for the development of real-time artificial neuron-glia networks opening new perspectives in pattern recognition systems and in brain neurophysiology.

A Model of Indispensability of a Large Glial Layer in Cerebrovascular Circulation

We formulate the problem of oxygen delivery to neural tissue as a problem of association. Input to a pool of neurons in one brain area must be matched in space and time with metabolic inputs from the vascular network via the glial network. We thus have a model in which neural, glial, and vascular layers are connected bidirectionally, in that order. Connections between neuro-glial and glial-vascular stages are trained by an unsupervised learning mechanism such that input to the neural layer is sustained by the precisely patterned delivery of metabolic inputs from the vascular layer via the glial layer. Simulations show that the capacity of such a system to sustain patterns is weak when the glial layer is absent. Capacity is higher when a glial layer is present and increases with the layer size. The proposed formulation of neurovascular interactions raises many intriguing questions about the role of glial cells in cerebral circulation.

Dynamical Mean Field Model of a Neural-Glial Mass

Our goal is to model the behavior of an ensemble of interacting neurons and astrocytes (the neural-glial mass). For this, a model describing N tripartite synapses is proposed. Each tripartite synapse consists of presynaptic and postsynaptic nerve terminals, as well as the synaptically associated astrocytic microdomain, and is described by a system of 13 stochastic differential equations. Then, by applying the dynamical mean field approximation (DMA) (Hasegawa, 2003a , 2003b ) the system of 13N equations is reduced to 13(13 + 2) = 195 deterministic differential equations for the means and the second-order moments of local and global variables. Simulations are carried out for studying the response of the neural-glial mass to external inputs applied to either the presynaptic terminals or the astrocytes. Three cases were considered: the astrocytes influence only the presynaptic terminal, only the postsynaptic terminal, or both the presynaptic and postsynaptic terminals. As a result, a wide range of responses varying from singles spikes to train of spikes was evoked on presynaptic and postsynaptic terminals. The experimentally observed phenomenon of spontaneous activity in astrocytes was replicated on the neural-glial mass. The model predicts that astrocytes can have a strong and activity-dependent influence on synaptic transmission. Finally, simulations show that the dynamics of astrocytes influences the synchronization ratio between neurons, predicting a peak in the synchronization for specific values of the astrocytes’ parameters.

Simulation of spontaneous Ca2+ oscillations in astrocytes mediated by voltage-gated calcium channels

The purpose of this computational study was to investigate the possible role of voltage-gated Ca(2+) channels in spontaneous Ca(2+) oscillations of astrocytes. By incorporating different types of voltage-gated Ca(2+) channels and a previous model, this study reproduced typical Ca(2+) oscillations in silico. Our model could mimic the oscillatory phenomenon under a wide range of experimental conditions, including resting membrane potential (-75 to -60 mV), extracellular Ca(2+) concentration (0.1 to 1500 muM), temperature (20 to 37 degrees C), and blocking specific Ca(2+) channels. By varying the experimental conditions, the amplitude and duration of Ca(2+) oscillations changed slightly (both <25%), while the frequency changed significantly ( approximately 400%). This indicates that spontaneous Ca(2+) oscillations in astrocytes might be an all-or-none process, which might be frequency-encoded in signaling. Moreover, the properties of Ca(2+) oscillations were found to be related to the dynamics of Ca(2+) influx, and not only to a constant influx. Therefore, calcium channels dynamics should be used in studying Ca(2+) oscillations. This work provides a platform to explore the still unclear mechanism of spontaneous Ca(2+) oscillations in astrocytes.

Dynamical patterns of calcium signaling in a functional model of neuron-astrocyte networks

We propose a functional mathematical model for neuron-astrocyte networks. The model incorporates elements of the tripartite synapse and the spatial branching structure of coupled astrocytes. We consider glutamate-induced calcium signaling as a specific mode of excitability and transmission in astrocytic-neuronal networks. We reproduce local and global dynamical patterns observed experimentally.

Epilepsy, regulation of brain energy metabolism and neurotransmission

Seizures are the result of a sudden and temporary synchronization of neuronal activity, the reason for which is not clearly understood. Astrocytes participate in the control of neurotransmitter storage and neurotransmission efficacy. They provide fuel to neurons, which need a high level of energy to sustain normal and pathological neuronal activities, such as during epilepsy. Various genetic or induced animal models have been developed and used to study epileptogenic mechanisms. Methionine sulfoximine induces both seizures and the accumulation of brain glycogen, which might be considered as a putative energy store to neurons in various animals. Animals subjected to methionine sulfoximine develop seizures similar to the most striking form of human epilepsy, with a long pre-convulsive period of several hours, a long convulsive period during up to 48 hours and a post convulsive period during which they recover normal behavior. The accumulation of brain glycogen has been demonstrated in both the cortex and cerebellum as early as the pre-convulsive period, indicating that this accumulation is not a consequence of seizures. The accumulation results from an activation of gluconeogenesis specifically localized to astrocytes, both in vivo and in vitro. Both seizures and brain glycogen accumulation vary when using different inbred strains of mice. C57BL/6J is the most "resistant" strain to methionine sulfoximine, while CBA/J is the most "sensitive" one. The present review describes the data obtained on methionine sulfoximine dependent seizures and brain glycogen in the light of neurotransmission, highlighting the relevance of brain glycogen content in epilepsies.

Signal propagation in cortical networks: a digital signal processing approach

This work reports a digital signal processing approach to representing and modeling transmission and combination of signals in cortical networks. The signal dynamics is modeled in terms of diffusion, which allows the information processing undergone between any pair of nodes to be fully characterized in terms of a finite impulse response (FIR) filter. Diffusion without and with time decay are investigated. All filters underlying the cat and macaque cortical organization are found to be of low-pass nature, allowing the cortical signal processing to be summarized in terms of the respective cutoff frequencies (a high cutoff frequency meaning little alteration of signals through their intermixing). Several findings are reported and discussed, including the fact that the incorporation of temporal activity decay tends to provide more diversified cutoff frequencies. Different filtering intensity is observed for each community in those networks. In addition, the brain regions involved in object recognition tend to present the highest cutoff frequencies for both the cat and macaque networks.

The influence of the astrocyte field on neuronal dynamics and synchronization

Astrocytes can sense local synaptic release of glutamate by metabotropic glutamate receptors. Receptor activation in turn can mediate transient increases of astrocytic intracellular calcium concentration through inositol 1,4,5-trisphosphate production. Notably, the perturbation of calcium concentration can propagate to other adjacent astrocytes. Astrocytic calcium signaling can therefore be linked to synaptic information transfer between neurons. On the other hand, astrocytes can also modulate neuronal activity by feeding back onto synaptic terminals in a fashion that depends on their intracellular calcium concentration. Thus, astrocytes can also be active partners in neuronal network activity. The aim of our study is to provide a computationally simple network model of mutual neuron-astrocyte interactions, in order to investigate the possible roles of astrocytes in neuronal network dynamics. In particular, we focus on the information entropy of neuronal firing of the whole network, considering how it could be affected by neuron-glial interactions.

Dynamics of a minimal neural model consisting of an astrocyte, a neuron, and an interneuron

In this paper, a biophysical neural network model consisting of a pyramidal neuron, an interneuron, and the astrocyte is studied. The corresponding dynamical properties are mainly investigated by using numerical simulations. The results show that the presence of the adenosine triphosphate and of the interneuron impacts the overall neural activity. It is shown that the fluxes of calcium through the cellular membrane strongly affect the modulation of the neural activity arising from the astrocyte.

On the role of synchrony for neuron-astrocyte interactions and perceptual conscious processing

Recent research on brain correlates of cognitive processes revealed the occurrence of global synchronization during conscious processing of sensory stimuli. In spite of technological progress in brain imaging, an explanation of the computational role of synchrony is still a highly controversial issue. In this study, we depart from an analysis of the usage of blood-oxygen-level-dependent functional magnetic resonance imaging for the study of cognitive processing, leading to the identification of evoked local field potentials as the vehicle for sensory patterns that compose conscious episodes. Assuming the "astrocentric hypothesis" formulated by James M. Robertson (astrocytes being the final stage of conscious processing), we propose that the role of global synchrony in perceptual conscious processing is to induce the transfer of information patterns embodied in local field potentials to astrocytic calcium waves, further suggesting that these waves are responsible for the "binding" of spatially distributed patterns into unitary conscious episodes.

Spontaneous calcium signals induced by gap junctions in a network model of astrocytes

The dynamics of a network model of astrocytes coupled by gap junctions is investigated. Calcium dynamics of the single cell is described by the biophysical model comprising the set of three nonlinear differential equations. Intercellular dynamics is provided by the diffusion of inositol 1,4,5-trisphosphate (IP3) through gap junctions between neighboring astrocytes. It is found that the diffusion induces the appearance of spontaneous activity patterns in the network. Stability of the network steady state is analyzed. It is proved that the increase of the diffusion coefficient above a certain critical value yields the generation of low-amplitude subthreshold oscillatory signals in a certain frequency range. It is shown that such spontaneous oscillations can facilitate calcium pulse generation and provide a certain time scale in astrocyte signaling.

The influence of sodium and potassium dynamics on excitability, seizures, and the stability of persistent states. II. Network and glial dynamics

In these companion papers, we study how the interrelated dynamics of sodium and potassium affect the excitability of neurons, the occurrence of seizures, and the stability of persistent states of activity. We seek to study these dynamics with respect to the following compartments: neurons, glia, and extracellular space. We are particularly interested in the slower time-scale dynamics that determine overall excitability, and set the stage for transient episodes of persistent oscillations, working memory, or seizures. In this second of two companion papers, we present an ionic current network model composed of populations of Hodgkin-Huxley type excitatory and inhibitory neurons embedded within extracellular space and glia, in order to investigate the role of micro-environmental ionic dynamics on the stability of persistent activity. We show that these networks reproduce seizure-like activity if glial cells fail to maintain the proper micro-environmental conditions surrounding neurons, and produce several experimentally testable predictions. Our work suggests that the stability of persistent states to perturbation is set by glial activity, and that how the response to such perturbations decays or grows may be a critical factor in a variety of disparate transient phenomena such as working memory, burst firing in neonatal brain or spinal cord, up states, seizures, and cortical oscillations.