Parameter Estimation Methods for Single Neuron Models

With the advancement in computer technology, it has become possible to fit complex models to neuronal data. In this work, we test how two methods can estimate parameters of simple neuron models (passive soma) to more complex ones (neuron with one dendritic cylinder and two active conductances). The first method uses classical voltage traces resulting from current pulses injection (time domain), while the second uses measures of the neuron's response to sinusoidal stimuli (frequency domain). Both methods estimate correctly the parameters in all cases studied. However, the time-domain method is slower and more prone to estimation errors in the cable parameters than the frequency-domain method. Because with noisy data the goodness of fit does not distinguish between different solutions, we suggest that running the estimation procedure a large number of times might help find a good solution and can provide information about the interactions between parameters. Also, because the formulation used for the model's response in the frequency domain is analytical, one can derive a local sensitivity analysis for each parameter. This analysis indicates how well a parameter is likely to be estimated and helps choose an optimal stimulation protocol. Finally, the tests suggest a strategy for fitting single-cell models using the two methods examined.

[1]  W. R. Foster,et al.  Significance of conductances in Hodgkin-Huxley models. , 1993, Journal of neurophysiology.

[2]  S. Grillner,et al.  Voltage-clamp frequency domain analysis of NMDA-activated neurons. , 1993, The Journal of experimental biology.

[3]  D. Durand,et al.  The somatic shunt cable model for neurons. , 1984, Biophysical journal.

[4]  L. Moore,et al.  White noise analysis of cable properties of neuroblastoma cells and lamprey central neurons. , 1985, Journal of neurophysiology.

[5]  Paul B. Manis,et al.  The parameter identification problem for the somatic shunt model , 1992, Biological Cybernetics.

[6]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1990 .

[7]  L. E. Moore,et al.  Quantitative Analysis of Electrotonic Structure and Membrane Properties of NMDA-Activated Lamprey Spinal Neurons , 1995, Neural Computation.

[8]  John E. Dennis,et al.  An Adaptive Nonlinear Least-Squares Algorithm , 1977, TOMS.

[9]  Alexander Borst,et al.  The intrinsic electrophysiological characteristics of fly lobula plate tangential cells: I. Passive membrane properties , 1996, Journal of Computational Neuroscience.

[10]  K J Staley,et al.  Membrane properties of dentate gyrus granule cells: comparison of sharp microelectrode and whole-cell recordings. , 1992, Journal of neurophysiology.

[11]  M. Kawato Cable properties of a neuron model with non-uniform membrane resistivity. , 1984, Journal of theoretical biology.

[12]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[13]  Alan C. Hindmarsh,et al.  A Polyalgorithm for the Numerical Solution of Ordinary Differential Equations , 1975, TOMS.

[14]  C S Peskin,et al.  Recovery of cable properties through active and passive modeling of subthreshold membrane responses from laterodorsal tegmental neurons. , 1998, Journal of neurophysiology.

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

[16]  F. Dodge,et al.  Subthreshold Behavior and Phenomenological Impedance of the Squid Giant Axon , 1970, The Journal of general physiology.

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

[18]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1952, The Journal of physiology.

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

[20]  William H. Press,et al.  Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition , 1987 .

[21]  T. Valiante,et al.  White noise approach for estimating the passive electrical properties of neurons. , 1996, Journal of neurophysiology.

[22]  R. Burke,et al.  Electrotonic characteristics of alpha motoneurones of varying size , 1971, The Journal of physiology.

[23]  J. Bower,et al.  Exploring parameter space in detailed single neuron models: simulations of the mitral and granule cells of the olfactory bulb. , 1993, Journal of neurophysiology.

[24]  Joël Tabak,et al.  Simulation and Parameter Estimation Study of a Simple Neuronal Model of Rhythm Generation: Role of NMDA and Non-NMDA Receptors , 1998, Journal of Computational Neuroscience.

[25]  G. Major,et al.  The modelling of pyramidal neurones in the visual cortex , 1989 .

[26]  D. Durand,et al.  Estimation of electronic parameters of neurons using an inverse Fourier transform technique , 1992, IEEE Transactions on Biomedical Engineering.

[27]  D. Zipser,et al.  Identification models of the nervous system , 1992, Neuroscience.