Fitting experimental data to models that use morphological data from public databases

Ideally detailed neuron models should make use of morphological and electrophysiological data from the same cell. However, this rarely happens. Typically a modeler will choose a cell morphology from a public database, assign standard values for Ra, Cm, and other parameters and then do the modeling study. The assumption is that the model will produce results representative of what might be obtained experimentally. To test this assumption we developed models of CA1 hippocampal pyramidal neurons using 4 different morphologies obtained from 3 public databases. The multiple run fitter in NEURON was used to fit parameter values in each of the 4 morphological models to match experimental data recorded from 19 CA1 pyramidal cells. Fits with fixed standard parameter values produced results that were generally not representative of our experimental data. However, when parameter values were allowed to vary, excellent fits were obtained in almost all cases, but the fitted parameter values were very different among the 4 reconstructions and did not match standard values. The differences in fitted values can be explained by very different diameters, total lengths, membrane areas and volumes among the reconstructed cells, reflecting either cell heterogeneity or issues with the reconstruction data. The fitted values compensated for these differences to make the database cells and experimental cells more similar electrotonically. We conclude that models using fully reconstructed morphologies need to be calibrated with experimental data (even when morphological and electrophysiological data come from the same cell), model results should be generated with multiple reconstructions, morphological and experimental cells should come from the same strain of animal at the same age, and blind use of standard parameter values in models that use reconstruction data may not produce representative experimental results.

[1]  I. Módy,et al.  Characterization of synaptically elicited GABAB responses using patch‐clamp recordings in rat hippocampal slices. , 1993, The Journal of physiology.

[2]  J. Schiller,et al.  NMDA spikes in basal dendrites of cortical pyramidal neurons , 2000, Nature.

[3]  B Sakmann,et al.  Detailed passive cable models of whole-cell recorded CA3 pyramidal neurons in rat hippocampal slices , 1994, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[4]  Daniel Johnston,et al.  Voltage-gated ion channels in dendrites of hippocampal pyramidal neurons , 2006, Pflügers Archiv.

[5]  R. Chitwood,et al.  Passive electrotonic properties of rat hippocampal CA3 interneurones , 1999, The Journal of physiology.

[6]  D. Johnston,et al.  K+ channel regulation of signal propagation in dendrites of hippocampal pyramidal neurons , 1997, Nature.

[7]  L. Grover,et al.  Evidence for postsynaptic induction and expression of NMDA receptor independent LTP. , 1998, Journal of neurophysiology.

[8]  J. Kao,et al.  Compartmentalized and Binary Behavior of Terminal Dendrites in Hippocampal Pyramidal Neurons , 2001, Science.

[9]  R C Cannon,et al.  Panel on: “Distributed operating” systems and distributed databases , 1987, 1987 IEEE Third International Conference on Data Engineering.

[10]  D. DiFrancesco,et al.  Action of the hyperpolarization-activated current (Ih) blocker ZD 7288 in hippocampal CA1 neurons , 1997, Pflügers Archiv.

[11]  Michele Migliore,et al.  Role of an A-Type K+ Conductance in the Back-Propagation of Action Potentials in the Dendrites of Hippocampal Pyramidal Neurons , 1999, Journal of Computational Neuroscience.

[12]  A. Constanti,et al.  Mechanism of block by ZD 7288 of the hyperpolarization-activated inward rectifying current in guinea pig substantia nigra neurons in vitro. , 1995, Journal of neurophysiology.

[13]  W. Rall Time constants and electrotonic length of membrane cylinders and neurons. , 1969, Biophysical journal.

[14]  W. Holmes Modeling the effect of glutamate diffusion and uptake on NMDA and non-NMDA receptor saturation. , 1995, Biophysical journal.

[15]  J. Zhu,et al.  Maturation of layer 5 neocortical pyramidal neurons: amplifying salient layer 1 and layer 4 inputs by Ca2+ action potentials in adult rat tuft dendrites , 2000, The Journal of physiology.

[16]  Erik De Schutter,et al.  Effects of variability in anatomical reconstruction techniques on models of synaptic integration by dendrites: a comparison of three internet archives , 2004, The European journal of neuroscience.

[17]  T. Freund,et al.  Total number and distribution of inhibitory and excitatory synapses on hippocampal CA1 pyramidal cells , 2001, Neuroscience.

[18]  L. Grover,et al.  Evidence for involvement of group II/III metabotropic glutamate receptors in NMDA receptor-independent long-term potentiation in area CA1 of rat hippocampus. , 1999, Journal of neurophysiology.

[19]  G. Ascoli,et al.  Quantitative morphometry of hippocampal pyramidal cells: Differences between anatomical classes and reconstructing laboratories , 2004, The Journal of comparative neurology.

[20]  W Rall,et al.  Matching dendritic neuron models to experimental data. , 1992, Physiological reviews.

[21]  D. Jaeger Accurate Reconstruction of Neuronal Morphology , 2000 .

[22]  I Segev,et al.  Untangling dendrites with quantitative models. , 2000, Science.

[23]  S. Hoffman,et al.  Funding for malaria genome sequencing , 1997, Nature.

[24]  D. Amaral,et al.  A quantitative analysis of the dendritic organization of pyramidal cells in the rat hippocampus , 1995, The Journal of comparative neurology.

[25]  Bartlett W. Mel,et al.  Dendrites: bug or feature? , 2003, Current Opinion in Neurobiology.

[26]  Arnold R. Kriegstein,et al.  Whole cell recording from neurons in slices of reptilian and mammalian cerebral cortex , 1989, Journal of Neuroscience Methods.

[27]  J. Magee Dendritic Hyperpolarization-Activated Currents Modify the Integrative Properties of Hippocampal CA1 Pyramidal Neurons , 1998, The Journal of Neuroscience.

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

[29]  Gongyu Y. Shen,et al.  Computational analysis of action potential initiation in mitral cell soma and dendrites based on dual patch recordings. , 1999, Journal of neurophysiology.

[30]  T. Sejnowski,et al.  [Letters to nature] , 1996, Nature.

[31]  Erik De Schutter,et al.  Computational neuroscience : realistic modeling for experimentalists , 2000 .

[32]  D DiFrancesco,et al.  Properties of the hyperpolarization-activated current in rat hippocampal CA1 pyramidal cells. , 1993, Journal of neurophysiology.

[33]  L. Nowak,et al.  The role of divalent cations in the N‐methyl‐D‐aspartate responses of mouse central neurones in culture. , 1988, The Journal of physiology.

[34]  D A Turner,et al.  Dendritic properties of hippocampal CA1 pyramidal neurons in the rat: Intracellular staining in vivo and in vitro , 1998, The Journal of comparative neurology.

[35]  Arjen van Ooyen,et al.  The effect of dendritic topology on firing patterns in model neurons , 2002, Network.

[36]  K. Harris,et al.  Three-dimensional structure of dendritic spines and synapses in rat hippocampus (CA1) at postnatal day 15 and adult ages: implications for the maturation of synaptic physiology and long-term potentiation. , 1992, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[37]  M. Häusser,et al.  Propagation of action potentials in dendrites depends on dendritic morphology. , 2001, Journal of neurophysiology.

[38]  M. Häusser,et al.  Compartmental models of rat cerebellar Purkinje cells based on simultaneous somatic and dendritic patch‐clamp recordings , 2001, The Journal of physiology.

[39]  Andreas T. Schaefer,et al.  Coincidence detection in pyramidal neurons is tuned by their dendritic branching pattern. , 2003, Journal of neurophysiology.

[40]  William R Holmes,et al.  Analysis and comparison of morphological reconstructions of hippocampal field CA1 pyramidal cells , 2005, Hippocampus.

[41]  D. Vasilyev,et al.  Postnatal Development of the Hyperpolarization-Activated Excitatory Current Ih in Mouse Hippocampal Pyramidal Neurons , 2002, The Journal of Neuroscience.

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

[43]  Richard Durbin,et al.  The computing neuron , 1989 .

[44]  D. Jaffe,et al.  Passive normalization of synaptic integration influenced by dendritic architecture. , 1999, Journal of neurophysiology.

[45]  N. Spruston,et al.  Determinants of Voltage Attenuation in Neocortical Pyramidal Neuron Dendrites , 1998, The Journal of Neuroscience.

[46]  Stuart D. Washington,et al.  Effects of dendritic morphology on CA3 pyramidal cell electrophysiology: a simulation study , 2002, Brain Research.

[47]  R. C Cannon,et al.  An on-line archive of reconstructed hippocampal neurons , 1998, Journal of Neuroscience Methods.