Computer science offers a large set of tools for prototyping, writing, running, testing, validating, sharing and reproducing results; however, computational science lags behind. In the best case, authors may provide their source code as a compressed archive and they may feel confident their research is reproducible. But this is not exactly true. James Buckheit and David Donoho proposed more than two decades ago that an article about computational results is advertising, not scholarship. The actual scholarship is the full software environment, code, and data that produced the result. This implies new workflows, in particular in peer-reviews. Existing journals have been slow to adapt: source codes are rarely requested and are hardly ever actually executed to check that they produce the results advertised in the article. ReScience is a peer-reviewed journal that targets computational research and encourages the explicit replication of already published research, promoting new and open-source implementations in order to ensure that the original research can be replicated from its description. To achieve this goal, the whole publishing chain is radically different from other traditional scientific journals. ReScience resides on GitHub where each new implementation of a computational study is made available together with comments, explanations, and software tests.

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.

More than a hundred biochemical species, activated by neurotransmitters binding to transmembrane receptors, are important in long-term potentiation (LTP) and long-term depression (LTD). To investigate which species and interactions are critical for synaptic plasticity, many computational postsynaptic signal transduction models have been developed. The models range from simple models with a single reversible reaction to detailed models with several hundred kinetic reactions. In this study, more than a hundred models are reviewed, and their features are compared and contrasted so that similarities and differences are more readily apparent. The models are classified according to the type of synaptic plasticity that is modeled (LTP or LTD) and whether they include diffusion or electrophysiological phenomena. Other characteristics that discriminate the models include the phase of synaptic plasticity modeled (induction, expression, or maintenance) and the simulation method used (deterministic or stochastic). We find that models are becoming increasingly sophisticated, by including stochastic properties, integrating with electrophysiological properties of entire neurons, or incorporating diffusion of signaling molecules. Simpler models continue to be developed because they are computationally efficient and allow theoretical analysis. The more complex models permit investigation of mechanisms underlying specific properties and experimental verification of model predictions. Nonetheless, it is difficult to fully comprehend the evolution of these models because (1) several models are not described in detail in the publications, (2) only a few models are provided in existing model databases, and (3) comparison to previous models is lacking. We conclude that the value of these models for understanding molecular mechanisms of synaptic plasticity is increasing and will be enhanced further with more complete descriptions and sharing of the published models.

No abstract

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.

The rate and temporal pattern of neural spiking each have the potential to influence computation. In the cerebellum, it has been hypothesized that the irregularity of interspike intervals in Purkinje cells affects their ability to transmit information to downstream neurons. Accordingly, during oculomotor behavior in mice and rhesus monkeys, mean irregularity of Purkinje cell spiking varied with mean eye velocity. However, moment-to-moment variations revealed a tight correlation between eye velocity and spike rate, with no additional information conveyed by spike irregularity. Moreover, when spike rate and irregularity were independently controlled using optogenetic stimulation, the eye movements elicited were well-described by a linear population rate code with 3–5 ms temporal precision. Biophysical and random-walk models identified biologically realistic parameter ranges that determine whether spike irregularity influences responses downstream. The results demonstrate cerebellar control of movements through a remarkably rapid rate code, with no evidence for an additional contribution of spike irregularity.

No abstract

BackgroundTo better understand the complex molecular level interactions seen in the pathogenesis of Alzheimer's disease, the results of the wet-lab and clinical studies can be complemented by mathematical models. Astrocytes are known to become reactive in Alzheimer's disease and their ionic equilibrium can be disturbed by interaction of the released and accumulated transmitters, such as serotonin, and peptides, including amyloid- peptides (A). We have here studied the effects of small amounts of A25–35 fragments on the transmitter-induced calcium signals in astrocytes by Fura-2AM fluorescence measurements and running simulations of the detected calcium signals.Methodology/Principal FindingsIntracellular calcium signals were measured in cultured rat cortical astrocytes following additions of serotonin and glutamate, or either of these transmitters together with A25–35. A25–35 increased the number of astrocytes responding to glutamate and exceedingly increased the magnitude of the serotonin-induced calcium signals. In addition to A25–35-induced effects, the contribution of intracellular calcium stores to calcium signaling was tested. When using higher stimulus frequency, the subsequent calcium peaks after the initial peak were of lower amplitude. This may indicate inadequate filling of the intracellular calcium stores between the stimuli. In order to reproduce the experimental findings, a stochastic computational model was introduced. The model takes into account the major mechanisms known to be involved in calcium signaling in astrocytes. Model simulations confirm the principal experimental findings and show the variability typical for experimental measurements.Conclusions/SignificanceNanomolar A25–35 alone does not cause persistent change in the basal level of calcium in astrocytes. However, even small amounts of A25–35, together with transmitters, can have substantial synergistic effects on intracellular calcium signals. Computational modeling further helps in understanding the mechanisms associated with intracellular calcium oscillations. Modeling the mechanisms is important, as astrocytes have an essential role in regulating the neuronal microenvironment of the central nervous system.