Computation in a Single Neuron: Hodgkin and Huxley Revisited

A spiking neuron computes by transforming a complex dynamical input into a train of action potentials, or spikes. The computation performed by the neuron can be formulated as dimensional reduction, or feature detection, followed by a nonlinear decision function over the low-dimensional space. Generalizations of the reverse correlation technique with white noise input provide a numerical strategy for extracting the relevant low-dimensional features from experimental data, and information theory can be used to evaluate the quality of the low-dimensional approximation. We apply these methods to analyze the simplest biophysically realistic model neuron, the Hodgkin-Huxley (HH) model, using this system to illustrate the general methodological issues. We focus on the features in the stimulus that trigger a spike, explicitly eliminating the effects of interactions between spikes. One can approximate this triggering feature space as a two-dimensional linear subspace in the high-dimensional space of input histories, capturing in this way a substantial fraction of the mutual information between inputs and spike time. We find that an even better approximation, however, is to describe the relevant subspace as two dimensional but curved; in this way, we can capture 90 of the mutual information even at high time resolution. Our analysis provides a new understanding of the computational properties of the HH model. While it is common to approximate neural behavior as integrate and fire, the HH model is not an integrator nor is it well described by a single threshold.

[1]  R. Reid,et al.  Predicting Every Spike A Model for the Responses of Visual Neurons , 2001, Neuron.

[2]  H. K. Hartline,et al.  THE RECEPTIVE FIELDS OF OPTIC NERVE FIBERS , 1940 .

[3]  D. Bray Protein molecules as computational elements in living cells , 1995, Nature.

[4]  F ROSENBLATT,et al.  The perceptron: a probabilistic model for information storage and organization in the brain. , 1958, Psychological review.

[5]  H. Barlow Summation and inhibition in the frog's retina , 1953, The Journal of physiology.

[6]  K. Sen,et al.  Spectral-temporal Receptive Fields of Nonlinear Auditory Neurons Obtained Using Natural Sounds , 2022 .

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

[8]  S. Yoshizawa,et al.  An Active Pulse Transmission Line Simulating Nerve Axon , 1962, Proceedings of the IRE.

[9]  Steven W. Zucker,et al.  Logical/Linear Operators for Image Curves , 1995, IEEE Trans. Pattern Anal. Mach. Intell..

[10]  Wulfram Gerstner,et al.  Reduction of the Hodgkin-Huxley Equations to a Single-Variable Threshold Model , 1997, Neural Computation.

[11]  R. FitzHugh Impulses and Physiological States in Theoretical Models of Nerve Membrane. , 1961, Biophysical journal.

[12]  G B Stanley,et al.  Reconstruction of Natural Scenes from Ensemble Responses in the Lateral Geniculate Nucleus , 1999, The Journal of Neuroscience.

[13]  R. Keynes The ionic channels in excitable membranes. , 1975, Ciba Foundation symposium.

[14]  Erkki Oja,et al.  Signal Separation by Nonlinear Hebbian Learning , 1995 .

[15]  Eve Marder,et al.  Reduction of conductance-based neuron models , 1992, Biological Cybernetics.

[16]  P Kuyper,et al.  Triggered correlation. , 1968, IEEE transactions on bio-medical engineering.

[17]  D. Hubel,et al.  Receptive fields, binocular interaction and functional architecture in the cat's visual cortex , 1962, The Journal of physiology.

[18]  W. Bialek,et al.  Features and dimensions: Motion estimation in fly vision , 2005, q-bio/0505003.

[19]  Adrienne L. Fairhall,et al.  What Can a Single Neuron Compute? , 2000, NIPS.

[20]  Amos Storkey,et al.  Advances in Neural Information Processing Systems 20 , 2007 .

[21]  William Bialek,et al.  Synergy in a Neural Code , 2000, Neural Computation.

[22]  William Bialek,et al.  Adaptive Rescaling Maximizes Information Transmission , 2000, Neuron.

[23]  Bernhard E. Boser,et al.  A training algorithm for optimal margin classifiers , 1992, COLT '92.

[24]  Naftali Tishby,et al.  The information bottleneck method , 2000, ArXiv.

[25]  T. Sejnowski,et al.  Comparison of current-driven and conductance-driven neocortical model neurons with Hodgkin-Huxley voltage-gated channels. , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[26]  Michael J. Berry,et al.  The Neural Code of the Retina , 1999, Neuron.

[27]  C. Koch,et al.  Encoding of visual information by LGN bursts. , 1999, Journal of neurophysiology.

[28]  James A. Simmons,et al.  A possible neuronal basis for representation of acoustic scenes in auditory cortex of the big brown bat , 1993, Nature.

[29]  Michael J. Berry,et al.  The structure and precision of retinal spike trains. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[30]  Idan Segev,et al.  Ion Channel Stochasticity May Be Critical in Determining the Reliability and Precision of Spike Timing , 1998, Neural Computation.

[31]  G D Lewen,et al.  Neural coding of naturalistic motion stimuli , 2001, Network.

[32]  Isabelle Guyon,et al.  Automatic Capacity Tuning of Very Large VC-Dimension Classifiers , 1992, NIPS.

[33]  R. Reid,et al.  Temporal Coding of Visual Information in the Thalamus , 2000, The Journal of Neuroscience.

[34]  L. Abbott,et al.  Model neurons: From Hodgkin-Huxley to hopfield , 1990 .

[35]  Adrienne L. Fairhall,et al.  Efficiency and ambiguity in an adaptive neural code , 2001, Nature.

[36]  Thomas M. Cover,et al.  Elements of Information Theory , 2005 .

[37]  M. Diamond,et al.  The Role of Spike Timing in the Coding of Stimulus Location in Rat Somatosensory Cortex , 2001, Neuron.

[38]  William Bialek,et al.  Spikes: Exploring the Neural Code , 1996 .

[39]  W. Bialek,et al.  Universal Statistical Behavior of Neural Spike Trains , 1998 .

[40]  S T Roweis,et al.  Nonlinear dimensionality reduction by locally linear embedding. , 2000, Science.

[41]  Vivien A. Casagrande,et al.  Biophysics of Computation: Information Processing in Single Neurons , 1999 .

[42]  William Bialek,et al.  Maximally Informative Dimensions: Analyzing Neural Responses to Natural Signals , 2002, NIPS.

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

[44]  Christof Koch,et al.  Biophysics of Computation: Information Processing in Single Neurons (Computational Neuroscience Series) , 1998 .

[45]  G D Lewen,et al.  Reproducibility and Variability in Neural Spike Trains , 1997, Science.

[46]  S. W. Kuffler Discharge patterns and functional organization of mammalian retina. , 1953, Journal of neurophysiology.

[47]  T. Sejnowski,et al.  Reliability of spike timing in neocortical neurons. , 1995, Science.

[48]  Adrienne L. Fairhall,et al.  What Causes a Neuron to Spike? , 2003, Neural Computation.

[49]  William Bialek,et al.  Real-time performance of a movement-sensitive neuron in the blowfly visual system: coding and information transfer in short spike sequences , 1988, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[50]  William Bialek,et al.  Analyzing Neural Responses to Natural Signals: Maximally Informative Dimensions , 2002, Neural Computation.

[51]  William Bialek,et al.  Statistical properties of spike trains: universal and stimulus-dependent aspects. , 1999, Physical review. E, Statistical, nonlinear, and soft matter physics.

[52]  H. Barlow,et al.  Retinal ganglion cells responding selectively to direction and speed of image motion in the rabbit , 1964, The Journal of physiology.

[53]  A. A. Mullin,et al.  Principles of neurodynamics , 1962 .