HippoUnit: A software tool for the automated testing and systematic comparison of detailed models of hippocampal neurons based on electrophysiological data

Anatomically and biophysically detailed data-driven neuronal models have become widely used tools for understanding and predicting the behavior and function of neurons. Due to the increasing availability of experimental data from anatomical and electrophysiological measurements as well as the growing number of computational and software tools that enable accurate neuronal modeling, there are now a large number of different models of many cell types available in the literature. These models were usually built to capture a few important or interesting properties of the given neuron type, and it is often unknown how they would behave outside their original context. In addition, there is currently no simple way of quantitatively comparing different models regarding how closely they match specific experimental observations. This limits the evaluation, re-use and further development of the existing models. Further, the development of new models could also be significantly facilitated by the ability to rapidly test the behavior of model candidates against the relevant collection of experimental data. We address these problems for the representative case of the CA1 pyramidal cell of the rat hippocampus by developing an open-source Python test suite, which makes it possible to automatically and systematically test multiple properties of models by making quantitative comparisons between the models and electrophysiological data. The tests cover various aspects of somatic behavior, and signal propagation and integration in apical dendrites. To demonstrate the utility of our approach, we applied our tests to compare the behavior of several different hippocampal CA1 pyramidal cell models from the ModelDB database against electrophysiological data available in the literature, and concluded that each of these models provides a good match to experimental results in some domains but not in others. We also show how we employed the test suite to aid the development of models within the European Human Brain Project (HBP), and describe the integration of the tests into the validation framework developed in the HBP, with the aim of facilitating more reproducible and transparent model building in the neuroscience community. Author summary Anatomically and biophysically detailed neuronal models are useful tools in neuroscience because they allow the prediction of the behavior and the function of the studied cell type under circumstances that are hard to investigate experimentally. However, most detailed biophysical models have been built to capture a few selected properties of the real neuron, and it is often unknown how they would behave under different circumstances, or whether they can be used to successfully answer different scientific questions. To help the modeling community develop better neural models, and make the process of model building more reproducible and transparent, we developed a test suite that enables the comparison of the behavior of models of neurons in the rat hippocampus and their evaluation against experimental data. Applying our tests to several models available in the literature, we show that each model is able to capture some of the important properties of the real neuron but fails to match experimental data in other domains. We also use the test suite in the model development workflow of the European Human Brain Project to aid the construction of better models of hippocampal neurons and networks.

[1]  Shawn D. Burton,et al.  NeuroElectro: a window to the world's neuron electrophysiology data , 2014, Front. Neuroinform..

[2]  I. Dittert,et al.  Temperature dependence of N-methyl-d-aspartate receptor channels and N-methyl-d-aspartate receptor excitatory postsynaptic currents , 2010, Neuroscience.

[3]  Andrew P Davison,et al.  The SONATA data format for efficient description of large-scale network models , 2020, PLoS computational biology.

[4]  Domenico Tegolo,et al.  A modeling study suggesting how a reduction in the context-dependent input on CA1 pyramidal neurons could generate schizophrenic behavior , 2011, Neural Networks.

[5]  Bartlett W. Mel,et al.  Distinguishing Linear vs. Non-Linear Integration in CA1 Radial Oblique Dendrites: It’s about Time , 2011, Front. Comput. Neurosci..

[6]  Ivan Raikov NineML – a description language for spiking neuron network modeling: the abstraction layer , 2010, BMC Neuroscience.

[7]  Szabolcs Káli,et al.  Data‐driven integration of hippocampal CA1 synaptic physiology in silico , 2020, Hippocampus.

[8]  Satish S. Nair,et al.  Neurons within the Same Network Independently Achieve Conserved Output by Differentially Balancing Variable Conductance Magnitudes , 2013, The Journal of Neuroscience.

[9]  M. Migliore,et al.  Functional significance of axonal Kv7 channels in hippocampal pyramidal neurons , 2008, Proceedings of the National Academy of Sciences.

[10]  A. Polsky,et al.  Submillisecond Precision of the Input-Output Transformation Function Mediated by Fast Sodium Dendritic Spikes in Basal Dendrites of CA1 Pyramidal Neurons , 2003, The Journal of Neuroscience.

[11]  Ivan Raikov,et al.  Interneuronal mechanisms of hippocampal theta oscillations in a full-scale model of the rodent CA1 circuit , 2016, eLife.

[12]  P. Dayan,et al.  A mathematical model explains saturating axon guidance responses to molecular gradients , 2016, eLife.

[13]  Henry Markram,et al.  BluePyOpt: Leveraging Open Source Software and Cloud Infrastructure to Optimise Model Parameters in Neuroscience , 2016, Front. Neuroinform..

[14]  J. Magee,et al.  Pathway Interactions and Synaptic Plasticity in the Dendritic Tuft Regions of CA1 Pyramidal Neurons , 2009, Neuron.

[15]  Romain Brette,et al.  Brian 2, an intuitive and efficient neural simulator , 2019, eLife.

[16]  Nicholas T. Carnevale,et al.  The NEURON Simulation Environment , 1997, Neural Computation.

[17]  Fiona E. N. LeBeau,et al.  Single-column thalamocortical network model exhibiting gamma oscillations, sleep spindles, and epileptogenic bursts. , 2005, Journal of neurophysiology.

[18]  J. Magee,et al.  Integrative Properties of Radial Oblique Dendrites in Hippocampal CA1 Pyramidal Neurons , 2006, Neuron.

[19]  E. Marder,et al.  Multiple models to capture the variability in biological neurons and networks , 2011, Nature Neuroscience.

[20]  J. Magee,et al.  Somatic EPSP amplitude is independent of synapse location in hippocampal pyramidal neurons , 2000, Nature Neuroscience.

[21]  Bartlett W. Mel,et al.  Pyramidal Neuron as Two-Layer Neural Network , 2003, Neuron.

[22]  Martin Styner,et al.  DTIPrep: quality control of diffusion-weighted images , 2014, Front. Neuroinform..

[23]  E. Marder,et al.  Variability, compensation and homeostasis in neuron and network function , 2006, Nature Reviews Neuroscience.

[24]  Tim Gollisch,et al.  Modeling Single-Neuron Dynamics and Computations: A Balance of Detail and Abstraction , 2006, Science.

[25]  Alexander Kozlov,et al.  The microcircuits of striatum in silico , 2020, Proceedings of the National Academy of Sciences.

[26]  Alois Knoll,et al.  The Human Brain Project: Creating a European Research Infrastructure to Decode the Human Brain , 2016, Neuron.

[27]  Addolorata Marasco,et al.  On the mechanisms underlying the depolarization block in the spiking dynamics of CA1 pyramidal neurons , 2012, Journal of Computational Neuroscience.

[28]  N. Spruston,et al.  Activity-dependent action potential invasion and calcium influx into hippocampal CA1 dendrites. , 1995, Science.

[29]  N. Spruston,et al.  Hippocampal Pyramidal Neurons Comprise Two Distinct Cell Types that Are Countermodulated by Metabotropic Receptors , 2012, Neuron.

[30]  Jeffrey M. Perkel,et al.  Why Jupyter is data scientists’ computational notebook of choice , 2018, Nature.

[31]  C. L. Rees,et al.  Hippocampome.org: a knowledge base of neuron types in the rodent hippocampus , 2015, eLife.

[32]  Sonja Grün,et al.  The Scientific Case for Brain Simulations , 2019, Neuron.

[33]  R. Nicoll,et al.  Analysis of excitatory synaptic action in pyramidal cells using whole‐cell recording from rat hippocampal slices. , 1990, The Journal of physiology.

[34]  Nace L. Golding,et al.  Compartmental Models Simulating a Dichotomy of Action Potential Backpropagation in Ca1 Pyramidal Neuron Dendrites , 2001, Journal of neurophysiology.

[35]  Pierre Yger,et al.  PyNN: A Common Interface for Neuronal Network Simulators , 2008, Front. Neuroinform..

[36]  Daniel Johnston,et al.  Intrinsic excitability of CA1 pyramidal neurones from the rat dorsal and ventral hippocampus , 2012, The Journal of physiology.

[37]  Szabolcs Káli,et al.  The physiological variability of channel density in hippocampal CA1 pyramidal cells and interneurons explored using a unified data-driven modeling workflow , 2018, PLoS Comput. Biol..

[38]  James G. King,et al.  The neocortical microcircuit collaboration portal: a resource for rat somatosensory cortex , 2015, Front. Neural Circuits.

[39]  C. Stevens,et al.  Voltage dependence of NMDA-activated macroscopic conductances predicted by single-channel kinetics , 1990, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[40]  Alon Korngreen,et al.  Is realistic neuronal modeling realistic? , 2016, Journal of neurophysiology.

[41]  Szabolcs Káli,et al.  Distinct properties of two major excitatory inputs to hippocampal pyramidal cells: a computational study , 2005, The European journal of neuroscience.

[42]  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..

[43]  H. Markram,et al.  Data-driven integration of hippocampal CA1 synapse physiology in silico , 2019, bioRxiv.

[44]  Gerkin Richard,et al.  NeuroUnit: Validation Tests for Neuroscience Models , 2013 .

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

[46]  Henry Markram,et al.  Models of Neocortical Layer 5b Pyramidal Cells Capturing a Wide Range of Dendritic and Perisomatic Active Properties , 2011, PLoS Comput. Biol..

[47]  A. Korngreen,et al.  A Quantitative Description of Dendritic Conductances and Its Application to Dendritic Excitation in Layer 5 Pyramidal Neurons , 2014, The Journal of Neuroscience.

[48]  Calvin J. Schneider,et al.  Toward a full-scale computational model of the rat dentate gyrus , 2012, Front. Neural Circuits.

[49]  Jeffrey C Magee,et al.  Sleep deprivation‐induced alterations in excitatory synaptic transmission in the CA1 region of the rat hippocampus , 2006, The Journal of physiology.

[50]  J. Magee,et al.  State-Dependent Dendritic Computation in Hippocampal CA1 Pyramidal Neurons , 2006, The Journal of Neuroscience.

[51]  Perry L. Miller,et al.  Twenty years of ModelDB and beyond: building essential modeling tools for the future of neuroscience , 2016, Journal of Computational Neuroscience.

[52]  E. Marder,et al.  Failure of averaging in the construction of a conductance-based neuron model. , 2002, Journal of neurophysiology.

[53]  E. Marder,et al.  Robust circuit rhythms in small circuits arise from variable circuit components and mechanisms , 2015, Current Opinion in Neurobiology.

[54]  Bartlett W. Mel,et al.  Arithmetic of Subthreshold Synaptic Summation in a Model CA1 Pyramidal Cell , 2003, Neuron.

[55]  Szabolcs Káli,et al.  A flexible, interactive software tool for fitting the parameters of neuronal models , 2014, Front. Neuroinform..

[56]  Richard C. Gerkin,et al.  Collaborative infrastructure for test-driven scientific model validation , 2014, ICSE Companion.

[57]  N. Spruston,et al.  Synapse Distribution Suggests a Two-Stage Model of Dendritic Integration in CA1 Pyramidal Neurons , 2009, Neuron.

[58]  Attila Losonczy,et al.  Associative pairing enhances action potential back‐propagation in radial oblique branches of CA1 pyramidal neurons , 2007, The Journal of physiology.

[59]  Henry Markram,et al.  A Novel Multiple Objective Optimization Framework for Constraining Conductance-Based Neuron Models by Experimental Data , 2007, Front. Neurosci..

[60]  R. Traub,et al.  A model of a CA3 hippocampal pyramidal neuron incorporating voltage-clamp data on intrinsic conductances. , 1991, Journal of neurophysiology.

[61]  James M. Bower,et al.  A Comparative Survey of Automated Parameter-Search Methods for Compartmental Neural Models , 1999, Journal of Computational Neuroscience.

[62]  James M. Bower,et al.  The 40-year history of modeling active dendrites in cerebellar Purkinje cells: emergence of the first single cell “community model” , 2015, Front. Comput. Neurosci..

[63]  I. Soltesz,et al.  Modulation of network behaviour by changes in variance in interneuronal properties , 2002, The Journal of physiology.

[64]  Nelson Spruston,et al.  Factors mediating powerful voltage attenuation along CA1 pyramidal neuron dendrites , 2005, The Journal of physiology.

[65]  Daniel Johnston,et al.  Mapping the electrophysiological and morphological properties of CA1 pyramidal neurons along the longitudinal hippocampal axis , 2016, Hippocampus.

[66]  N Spruston,et al.  Resting and active properties of pyramidal neurons in subiculum and CA1 of rat hippocampus. , 2000, Journal of neurophysiology.

[67]  Eduardo Ros,et al.  Modeling the Cerebellar Microcircuit: New Strategies for a Long-Standing Issue , 2016, Front. Cell. Neurosci..