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.

[1]  Vladimir Parpura,et al.  A possible role of astrocytes in contextual memory retrieval: An analysis obtained using a quantitative framework , 2013, Front. Comput. Neurosci..

[2]  Pierre J Magistretti,et al.  Comment on Recent Modeling Studies of Astrocyte–Neuron Metabolic Interactions , 2010, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[3]  R A John Challiss,et al.  Determinants of Metabotropic Glutamate Receptor-5-mediated Ca2+ and Inositol 1,4,5-Trisphosphate Oscillation Frequency , 2002, The Journal of Biological Chemistry.

[4]  Rogier Min,et al.  The computational power of astrocyte mediated synaptic plasticity , 2012, Front. Comput. Neurosci..

[5]  Hugh D. Spence,et al.  Minimum information requested in the annotation of biochemical models (MIRIAM) , 2005, Nature Biotechnology.

[6]  Christof Koch,et al.  Worldwide initiatives to advance brain research , 2016, Nature Neuroscience.

[7]  Peter J. Hunter,et al.  The CellML Model Repository , 2008, Bioinform..

[8]  Marc-Oliver Gewaltig,et al.  Towards Reproducible Descriptions of Neuronal Network Models , 2009, PLoS Comput. Biol..

[9]  E. Ben-Jacob,et al.  Glutamate regulation of calcium and IP3 oscillating and pulsating dynamics in astrocytes , 2009, Journal of biological physics.

[10]  Kristofer E. Bouchard,et al.  High-Performance Computing in Neuroscience for Data-Driven Discovery, Integration, and Dissemination , 2016, Neuron.

[11]  H. Parri,et al.  Spontaneous astrocytic Ca2+ oscillations in situ drive NMDAR-mediated neuronal excitation , 2001, Nature Neuroscience.

[12]  C. Rose,et al.  Astrocyte calcium signals at Schaffer collateral to CA1 pyramidal cell synapses correlate with the number of activated synapses but not with synaptic strength , 2012, Hippocampus.

[13]  Olli Yli-Harja,et al.  Simulation tools for biochemical networks: evaluation of performance and usability , 2005, Bioinform..

[14]  Michael L. Hines,et al.  ModelDB - Making models publicly accessible to support computational neuroscience , 2003, Neuroinformatics.

[15]  Michael L. Hines,et al.  NeuroML: A Language for Describing Data Driven Models of Neurons and Networks with a High Degree of Biological Detail , 2010, PLoS Comput. Biol..

[16]  Richard H. Middleton,et al.  Spatial Quantification of Cytosolic Ca$^{2+}$ Accumulation in Nonexcitable Cells:An Analytical Study , 2014, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[17]  Tiina Manninen,et al.  Modeling Signal Transduction Leading to Synaptic Plasticity: Evaluation and Comparison of Five Models , 2011, EURASIP J. Bioinform. Syst. Biol..

[18]  R Heinrich,et al.  Intercellular Ca2+ wave propagation through gap-junctional Ca2+ diffusion: a theoretical study. , 2001, Biophysical journal.

[19]  Liam McDaid,et al.  Biophysically based computational models of astrocyte ~ neuron coupling and their functional significance , 2013, Front. Comput. Neurosci..

[20]  Yutaka Sakai,et al.  Similarity in Neuronal Firing Regimes across Mammalian Species , 2016, The Journal of Neuroscience.

[21]  Todd A Fiacco,et al.  Supporting Online Material Materials and Methods Som Text Figs. S1 to S4 Hippocampal Short-and Long-term Plasticity Are Not Modulated by Astrocyte Ca 2+ Signaling , 2022 .

[22]  Robert C. Cannon,et al.  LEMS: a language for expressing complex biological models in concise and hierarchical form and its use in underpinning NeuroML 2 , 2014, Front. Neuroinform..

[23]  Tiina Manninen,et al.  Computational Models of Astrocytes and Astrocyte–Neuron Interactions: Characterization, Reproducibility, and Future Perspectives , 2019, Springer Series in Computational Neuroscience.

[24]  Erik De Schutter Why Are Computational Neuroscience and Systems Biology So Separate? , 2008, PLoS Comput. Biol..

[25]  L. Venance,et al.  Control and Plasticity of Intercellular Calcium Waves in Astrocytes: A Modeling Approach , 2002, The Journal of Neuroscience.

[26]  Martin D. Haustein,et al.  Conditions and Constraints for Astrocyte Calcium Signaling in the Hippocampal Mossy Fiber Pathway , 2014, Neuron.

[27]  Nicolas P. Rougier,et al.  A long journey into reproducible computational neuroscience , 2015, Front. Comput. Neurosci..

[28]  Jeanette Kotaleski,et al.  Postsynaptic Signal Transduction Models for Long-Term Potentiation and Depression , 2010, Front. Comput. Neurosci..

[29]  Hiroaki Kitano,et al.  The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models , 2003, Bioinform..

[30]  Maxim Lavrentovich,et al.  A mathematical model of spontaneous calcium(II) oscillations in astrocytes. , 2008, Journal of theoretical biology.

[31]  W. Gibson,et al.  Origins of blood volume change due to glutamatergic synaptic activity at astrocytes abutting on arteriolar smooth muscle cells. , 2008, Journal of theoretical biology.

[32]  Angelo Di Garbo,et al.  Calcium signalling in astrocytes and modulation of neural activity , 2007, Biosyst..

[33]  Marja-Leena Linne,et al.  Comparison of Models for IP3 Receptor Kinetics Using Stochastic Simulations , 2013, PloS one.

[34]  Sonja Grün,et al.  Handling Metadata in a Neurophysiology Laboratory , 2016, Front. Neuroinform..

[35]  Edmund J. Crampin,et al.  Minimum Information About a Simulation Experiment (MIASE) , 2011, PLoS Comput. Biol..

[36]  Eduardo Soriano,et al.  Neuronal Activity Regulates Correlated Network Properties of Spontaneous Calcium Transients in Astrocytes In Situ , 2002, The Journal of Neuroscience.

[37]  Eshel Ben-Jacob,et al.  Computational quest for understanding the role of astrocyte signaling in synaptic transmission and plasticity , 2012, Front. Comput. Neurosci..

[38]  Alfonso Araque,et al.  Astrocyte calcium signal and gliotransmission in human brain tissue. , 2013, Cerebral cortex.

[39]  A Goldbeter,et al.  One-pool model for Ca2+ oscillations involving Ca2+ and inositol 1,4,5-trisphosphate as co-agonists for Ca2+ release. , 1993, Cell calcium.

[40]  H. Markram,et al.  The neural code between neocortical pyramidal neurons depends on neurotransmitter release probability. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[41]  Huda Akil,et al.  Neuroscience Training for the 21st Century , 2016, Neuron.

[42]  James A. Bednar,et al.  An automated and reproducible workflow for running and analyzing neural simulations using Lancet and IPython Notebook , 2013, Front. Neuroinform..

[43]  J. Keizer,et al.  A single-pool inositol 1,4,5-trisphosphate-receptor-based model for agonist-stimulated oscillations in Ca2+ concentration. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[44]  Eshel Ben-Jacob,et al.  Sparse short-distance connections enhance calcium wave propagation in a 3D model of astrocyte networks , 2014, Front. Comput. Neurosci..

[45]  Eshel Ben-Jacob,et al.  Astrocyte regulation of sleep circuits: experimental and modeling perspectives , 2012, Front. Comput. Neurosci..

[46]  A. Atri,et al.  A single-pool model for intracellular calcium oscillations and waves in the Xenopus laevis oocyte. , 1993, Biophysical journal.

[47]  S. Nadkarni,et al.  Spontaneous oscillations of dressed neurons: a new mechanism for epilepsy? , 2003, Physical review letters.

[48]  Vladimir Parpura,et al.  Data and model tango to aid the understanding of astrocyte-neuron signaling , 2014, Front. Comput. Neurosci..

[49]  Leone Fronzoni,et al.  The influence of the astrocyte field on neuronal dynamics and synchronization , 2009, Journal of biological physics.

[50]  Jacky L. Snoep,et al.  Reproducible computational biology experiments with SED-ML - The Simulation Experiment Description Markup Language , 2011, BMC Systems Biology.

[51]  Geneviève Dupont,et al.  A model for Ca2+ oscillations stimulated by the type 5 metabotropic glutamate receptor: an unusual mechanism based on repetitive, reversible phosphorylation of the receptor. , 2011, Biochimie.

[52]  D E Postnov,et al.  Dynamical patterns of calcium signaling in a functional model of neuron–astrocyte networks , 2009, Journal of biological physics.

[53]  Catherine M Lloyd,et al.  CellML: its future, present and past. , 2004, Progress in biophysics and molecular biology.

[54]  Jacky L. Snoep,et al.  BioModels Database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems , 2005, Nucleic Acids Res..

[55]  J Riera,et al.  Modeling the spontaneous Ca2+ oscillations in astrocytes: Inconsistencies and usefulness. , 2011, Journal of integrative neuroscience.

[56]  Marja-Leena Linne,et al.  Astrocyte-neuron interactions: from experimental research-based models to translational medicine. , 2014, Progress in molecular biology and translational science.

[57]  Terrence J. Sejnowski,et al.  Computational models of neuron-astrocyte interaction in epilepsy , 2012, Front. Comput. Neurosci..

[58]  James G. King,et al.  Reconstruction and Simulation of Neocortical Microcircuitry , 2015, Cell.

[59]  B. J. Roth,et al.  A mathematical model of agonist-induced propagation of calcium waves in astrocytes. , 1995, Cell calcium.

[60]  H. Parri,et al.  The role of Ca2+ in the generation of spontaneous astrocytic Ca2+ oscillations , 2003, Neuroscience.

[61]  Nicholas T. Carnevale,et al.  ModelDB: A Database to Support Computational Neuroscience , 2004, Journal of Computational Neuroscience.

[62]  P. Jung,et al.  Stochastic properties of Ca(2+) release of inositol 1,4,5-trisphosphate receptor clusters. , 2002, Biophysical journal.

[63]  Geneviève Dupont,et al.  Spatiotemporal organization of Ca2+ dynamics: A modeling‐based approach , 2010, HFSP journal.

[64]  A Goldbeter,et al.  Bursting, chaos and birhythmicity originating from self-modulation of the inositol 1,4,5-trisphosphate signal in a model for intracellular Ca2+ oscillations , 1999, Bulletin of mathematical biology.

[65]  Liam McDaid,et al.  Self-repair in a bidirectionally coupled astrocyte-neuron (AN) system based on retrograde signaling , 2012, Front. Comput. Neurosci..

[66]  J. Rinzel,et al.  Equations for InsP3 receptor-mediated [Ca2+]i oscillations derived from a detailed kinetic model: a Hodgkin-Huxley like formalism. , 1994, Journal of theoretical biology.

[67]  T. Ozaki,et al.  Quantifying the uncertainty of spontaneous Ca2+ oscillations in astrocytes: particulars of Alzheimer's disease. , 2011, Biophysical journal.

[68]  Andrew P. Davison,et al.  Learning from the Past: Approaches for Reproducibility in Computational Neuroscience , 2013 .

[69]  N. Brunel,et al.  Astrocytes: Orchestrating synaptic plasticity? , 2015, Neuroscience.

[70]  Kenneth D. Harris,et al.  Data Sharing for Computational Neuroscience , 2008, Neuroinformatics.

[71]  Michael L. Hines,et al.  Interoperability of Neuroscience Modeling Software: Current Status and Future Directions , 2007, Neuroinformatics.

[72]  Nicolas Liaudet,et al.  Astrocyte Ca2+ signalling: an unexpected complexity , 2014, Nature Reviews Neuroscience.