Firing-rate models for neurons with a broad repertoire of spiking behaviors
© Heiberg et al; licensee BioMed Central Ltd. 2013
Published: 8 July 2013
Instantaneous firing rates are commonly used to describe either the compound spiking activity of neuron ensembles (population rate) or the trial-averaged response of individual neurons to multiple repetitions of the same stimulus. The dynamics of firing rates is often studied by means of population or firing-rate models. The main motivation for using such models rather than spiking neuron models is to reduce the dimensionality and complexity of the microscopic dynamics to allow analytical tractability, efficient simulation, and intuitive understanding. In most cases, population models have to be treated as purely phenomenological descriptions. Only under simplifying assumptions can they be derived or extracted from the single-neuron dynamics.
Previous studies in this field were mainly restricted to homogeneous populations of simple integrate-and-fire neurons which can, to some extent, mimic the rather stereotypical behaviour of regularly firing neurons (e.g., pyramidal cells in the neocortex). Electrophysiological studies have revealed a large diversity of neuron types with different dynamical behaviours, ranging from regular spiking to bursting cells, from accommodating to non-accommodating cells, from integrators to resonators, from depolarization- to hyperpolarization-activated cells, etc. So far, the effect of this neuron-type diversity on the dynamics of neural populations is unclear.
Here, we investigate the rate-response characteristics for a variety of neuron types in simulation experiments. To this end, we extend previous work on simple integrate-and-fire neurons [1, 2] to the multi-timescale adaptive threshold (MAT)  and the Izhikevich model . The model neurons are probed by spike trains modeled as realizations of inhomogeneous Poisson or Gamma point processes with sinusoidally modulated intensity. Average rates, modulation amplitudes and phases of the period-averaged spike responses are measured for a broad range of stimulus, neuron and synapse parameters. The stationary and nonstationary response properties are parameterized by means of nonlinear activation functions and linear filter kernels, respectively, serving as ingredients for simple linear-nonlinear firing rate models. We show that linear-nonlinear rate models can accurately predict the population responses to novel test stimuli for a broad range of parameters and neuron types.
Supported by the Research Council of Norway grant 178892/V30 eNeuro, the Helmholtz Alliance on Systems Biology, the Helmholtz Association in the Portfolio theme "Supercomputing and Modeling for the Human Brain", the Jülich Aachen Research Alliance (JARA), and EU Grant 269921 BrainScaleS. Simulations and data analysis were carried out using the NEST simulation tool (http://www.nest-initiative.org) and Python (http://www.python.org), using NOTUR computing resources.
- Nordlie E, Tetzlaff T, Einevoll GT: Rate Dynamics of Leaky Integrate-and-Fire Neurons with Strong Synapses. Front Comput Neurosci. 2010, 4: 149-PubMed CentralPubMedGoogle Scholar
- Heiberg T, Kriener B, Tetzlaff T, Casti A, Einevoll GT, Plesser HE: Firing-rate models for the retina-LGN connection derived from experiments. Neuroscience Meeting Planner. 2012, New Orleans, LA: Society for Neuroscience, 880.28.2012. OnlineGoogle Scholar
- Kobayashi R, Tsubo Y, Shinomoto S: Made-to-order spiking neuron model equipped with a multi-timescale adaptive threshold. Front Comput Neurosci. 2009, 3: 9PubMed CentralView ArticlePubMedGoogle Scholar
- Izhikevich EM: Simple Model of Spiking Neurons. 2003, IEEE Trans Neural Netw, 14: 1569-1572.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.