Neuronal model with distributed delay: analysis and simulation study for gamma distribution memory kernel

A single neuronal model incorporating distributed delay (memory)is proposed. The stochastic model has been formulated as a Stochastic Integro-Differential Equation (SIDE) which results in the underlying process being non-Markovian. A detailed analysis of the model when the distributed delay kernel has exponential form (weak delay) has been carried out. The selection of exponential kernel has enabled the transformation of the non-Markovian model to a Markovian model in an extended state space. For the study of First Passage Time (FPT) with exponential delay kernel, the model has been transformed to a system of coupled Stochastic Differential Equations (SDEs) in two-dimensional state space. Simulation studies of the SDEs provide insight into the effect of weak delay kernel on the Inter-Spike Interval(ISI) distribution. A measure based on Jensen–Shannon divergence is proposed which can be used to make a choice between two competing models viz. distributed delay model vis-á-vis LIF model. An interesting feature of the model is that the behavior of (CV(t))(ISI) (Coefficient of Variation) of the ISI distribution with respect to memory kernel time constant parameter η reveals that neuron can switch from a bursting state to non-bursting state as the noise intensity parameter changes. The membrane potential exhibits decaying auto-correlation structure with or without damped oscillatory behavior depending on the choice of parameters. This behavior is in agreement with empirically observed pattern of spike count in a fixed time window. The power spectral density derived from the auto-correlation function is found to exhibit single and double peaks. The model is also examined for the case of strong delay with memory kernel having the form of Gamma distribution. In contrast to fast decay of damped oscillations of the ISI distribution for the model with weak delay kernel, the decay of damped oscillations is found to be slower for the model with strong delay kernel.

[1]  Eugene M. Izhikevich,et al.  Resonate-and-fire neurons , 2001, Neural Networks.

[2]  A. Holden.  Models of the stochastic activity of neurones , 1976 .

[3]  Christof Koch,et al.  Computation and the single neuron , 1997, Nature.

[4]  K. Gopalsamy,et al.  Stability in asymmetric Hopfield nets with transmission delays , 1994 .

[5]  Wang,et al.  Global stability and local stability of Hopfield neural networks with delays. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[6]  R. Llinás.  The intrinsic electrophysiological properties of mammalian neurons: insights into central nervous system function. , 1988, Science.

[7]  M. Konishi,et al.  Neural map of interaural phase difference in the owl's brainstem. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Terry Williams,et al.  Probability and Statistics with Reliability, Queueing and Computer Science Applications , 1983 .

[9]  Bruno A. Olshausen,et al.  Book Review , 2003, Journal of Cognitive Neuroscience.

[10]  Y. Kuang.  Delay Differential Equations: With Applications in Population Dynamics , 2012 .

[11]  C. Koch,et al.  Methods in Neuronal Modeling: From Ions to Networks , 1998 .

[12]  Petr Lánský,et al.  Similarity of interspike interval distributions and information gain in a stationary neuronal firing , 2006, Biological Cybernetics.

[13]  Christof Koch,et al.  Do neurons have a voltage or a current threshold for action potential initiation? , 1995, Journal of Computational Neuroscience.

[14]  S. Ruan.  DELAY DIFFERENTIAL EQUATIONS IN SINGLE SPECIES DYNAMICS , 2006 .

[15]  P Lánský,et al.  On approximations of Stein's neuronal model. , 1984, Journal of theoretical biology.

[16]  R. Miura,et al.  Low-threshold calcium current and resonance in thalamic neurons: a model of frequency preference. , 1994, Journal of neurophysiology.

[17]  Jaime de la Rocha,et al.  Supplementary Information for the article ‘ Correlation between neural spike trains increases with firing rate ’ , 2007 .

[18]  Néstor Parga,et al.  Auto- and crosscorrelograms for the spike response of leaky integrate-and-fire neurons with slow synapses. , 2006, Physical review letters.

[19]  U. Heinemann,et al.  Dynamics of rat entorhinal cortex layer II and III cells: characteristics of membrane potential resonance at rest predict oscillation properties near threshold , 2004, The Journal of physiology.

[20]  L Schimansky-Geier,et al.  First passage time densities in resonate-and-fire models. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  Peter Dayan,et al.  Theoretical Neuroscience: Computational and Mathematical Modeling of Neural Systems , 2001 .

[22]  N. Macdonald.  Time lags in biological models , 1978 .

[23]  W. R. Buckland,et al.  Stochastic Models for Social Processes , 1967 .

[24]  Amir F. Atiya,et al.  How delays affect neural dynamics and learning , 1994, IEEE Trans. Neural Networks.

[25]  Carson C. Chow,et al.  Noise shaping in populations of coupled model neurons. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[26]  R. Llinás,et al.  In vitro neurons in mammalian cortical layer 4 exhibit intrinsic oscillatory activity in the 10- to 50-Hz frequency range. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[27]  T. Verechtchaguina,et al.  First passage time densities in non-Markovian models with subthreshold oscillations , 2006 .

[28]  Robert Zwanzig,et al.  Hydrodynamic Theory of the Velocity Correlation Function , 1970 .

[29]  L. Ricciardi,et al.  The Ornstein-Uhlenbeck process as a model for neuronal activity , 1979, Biological Cybernetics.

[30]  Petr Lánský,et al.  The parameters of the stochastic leaky integrate-and-fire neuronal model , 2006, Journal of Computational Neuroscience.

[31]  B. Scott Jackson,et al.  Including Long-Range Dependence in Integrate-and-Fire Models of the High Interspike-Interval Variability of Cortical Neurons , 2004, Neural Computation.

[32]  D. Bartholomew,et al.  Stochastic Models for Social Processes. , 1970 .

[33]  Subrahmanyan Chandrasekhar,et al.  Brownian Motion, Dynamical Friction, and Stellar Dynamics , 1949 .

[34]  Xiaofeng Liao,et al.  Bifurcation analysis on a two-neuron system with distributed delays in the frequency domain , 2004, Neural Networks.

[35]  K. Vahala.  Handbook of stochastic methods for physics, chemistry and the natural sciences , 1986, IEEE Journal of Quantum Electronics.

[36]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[37]  Fabrizio Gabbiani,et al.  Principles of spike train analysis , 1996 .

[38]  Hagai Bergman,et al.  The neuronal refractory period causes a short-term peak in the autocorrelation function , 2001, Journal of Neuroscience Methods.

[39]  R. Llinás,et al.  Dendritic calcium conductances generate high-frequency oscillation in thalamocortical neurons. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[40]  J. Goodisman.  Hellmann-Feynman Theorem in Thomas-Fermi and Related Theories , 1970 .

[41]  H. Tuckwell.  Introduction to Theoretical Neurobiology: Linear Cable Theory and Dendritic Structure , 1988 .

[42]  J J Hopfield,et al.  Neural computation by concentrating information in time. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[43]  Jawaharlal Karmeshu,et al.  Entropy Measures, Maximum Entropy Principle and Emerging Applications , 2003 .

[44]  B. Mandelbrot,et al.  RANDOM WALK MODELS FOR THE SPIKE ACTIVITY OF A SINGLE NEURON. , 1964, Biophysical journal.

[45]  J Rinzel,et al.  Influence of temporal correlation of synaptic input on the rate and variability of firing in neurons. , 2000, Biophysical journal.

[46]  P Lánský,et al.  First-passage-time problem for simulated stochastic diffusion processes. , 1994, Computers in biology and medicine.

[47]  Hal L. Smith,et al.  An introduction to delay differential equations with applications to the life sciences / Hal Smith , 2010 .

[48]  Kishor S. Trivedi.  Probability and Statistics with Reliability, Queuing, and Computer Science Applications , 1984 .

[49]  Robert G. Turcott,et al.  Long-duration correlation in the sequence of action potentials in an insect visual interneuron , 1995 .

[50]  José Carlos Príncipe,et al.  The gamma model--A new neural model for temporal processing , 1992, Neural Networks.

[51]  H. Haken.  Synergetics: Introduction and Advanced Topics , 2004 .

[52]  R. Stein.  A THEORETICAL ANALYSIS OF NEURONAL VARIABILITY. , 1965, Biophysical journal.

[53]  Henry C. Tuckwell,et al.  Introduction to theoretical neurobiology , 1988 .

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

[55]  J. Bélair.  Stability in a model of a delayed neural network , 1993 .

[56]  K. Gopalsamy,et al.  Convergence under dynamical thresholds with delays , 1997, IEEE Trans. Neural Networks.

[57]  Christof Koch,et al.  Detecting and Estimating Signals in Noisy Cable Structures, I: Neuronal Noise Sources , 1999, Neural Computation.

[58]  C. Zucca,et al.  A Monte Carlo Method for the Simulation of First Passage Times of Diffusion Processes , 2001 .

[59]  S. Levin.  Lectu re Notes in Biomathematics , 1983 .

[60]  Eduardo R. Caianiello,et al.  Decision equation for binary systems. Application to neuronal behavior , 2004, Kybernetik.