this paper we propose that the irregular ISI arises as a consequence of a specific problem that cortical neurons must solve: the problem of dynamic range or gain control. Cortical neurons receive 3000--10,000 synaptic contacts, 85% of which are asymmetric and hence presumably excitatory (Peters, 1987; Braitenberg and Schuz, 1991). More than half of these contacts are thought to arise from neurons within a 100--200 #m radius of the target cell, reflecting the stereotypical columnar organization of neocortex. Because neurons within a cortical column typically share similar physiological properties, the conditions that excite one neuron are likely to excite a considerable fraction of its afferent input as well (Mountcastle, 1978; Peters and Sethares, 1991), creating a scenario in which saturation of the neuron's firing rate could easily occur. This problem is exacerbated by the fact that EPSPs from individual axons appear to exert substantial impact on the membrane potential (Mason et al., 1991; Otmakhov Received Sept. 15, 1997; revised Feb. 25, 1998; accepted March 3, 1998.

Several formulations of correlation-based Hebbian learning are reviewed. On the presynaptic side, activity is described either by a firing rate or by presynaptic spike arrival. The state of the postsynaptic neuron can be described by its membrane potential, its firing rate, or the timing of backpropagating action potentials (BPAPs). It is shown that all of the above formulations can be derived from the point of view of an expansion. In the absence of BPAPs potentials, it is natural to correlate presynaptic spikes with the postsynaptic membrane potential. Time windows of spike time dependent plasticity arise naturally, if the timing of postsynaptic spikes is available at the site of the synapse as it is the case in the presence of BPAPs. With an appropriate choice of parameters, Hebbian synaptic plasticity has intrinsic normalization properties that stabilizes postsynaptic firing rates and leads to subtractive weight normalization.

We demonstrate that single-variable integrate-and-fire models can quantitatively capture the dynamics of a physiologically-detailed model for fast-spiking cortical neurons. Through a systematic set of approximations, we reduce the conductance based model to two variants of integrate-and-fire models. In the first variant (non-linear integrate-and-fire model), parameters depend on the instantaneous membrane potential whereas in the second variant, they depend on the time elapsed since the last spike (Spike Response Model). The direct reduction links features of the simple models to biophysical features of the full conductance based model. To quantitatively

We analyze the effect of noise in integrate-and-fire neurons driven by timedependent input, and compare the diffusion approximation for the membrane potential to escape noise. It is shown that for time-dependent sub-threshold input, diffusive noise can be replaced by escape noise with a hazard function that has a Gaussian dependence upon the distance between the (noise-free) membrane voltage and threshold. The approximation is improved if we add to the hazard function a probability current proportional to the derivative of the voltage. Stochastic resonance in response to periodic input occurs in both noise models and exhibits similar characteristics.

We propose a biophysical mechanism for the high interspike interval variability observed in cortical spike trains. The key lies in the nonlinear dynamics of cortical spike generation, which are consistent with type I membranes where saddle-node dynamics underlie excitability (Rinzel & Ermentrout, 1989). We present a canonical model for type I membranes, the θ-neuron. The θ-neuron is a phase model whose dynamics reflect salient features of type I membranes. This model generates spike trains with coefficient of variation (CV) above 0.6 when brought to firing by noisy inputs. This happens because the timing of spikes for a type I excitable cell is exquisitely sensitive to the amplitude of the suprathreshold stimulus pulses. A noisy input current, giving random amplitude “kicks” to the cell, evokes highly irregular firing across a wide range of firing rates; an intrinsically oscillating cell gives regular spike trains. We corroborate the results with simulations of the Morris-Lecar (M-L) neural model with random synaptic inputs: type I M-L yields high CVs. When this model is modified to have type II dynamics (periodicity arises via a Hopf bifurcation), however, it gives regular spike trains (CV below 0.3). Our results suggest that the high CV values such as those observed in cortical spike trains are an intrinsic characteristic of type I membranes driven to firing by “random” inputs. In contrast, neural oscillators or neurons exhibiting type II excitability should produce regular spike trains.

you have obtained prior permission, you may not download an entire issue of a journal or multiple copies of articles, and you may use content in the JSTOR archive only for your personal, non-commercial use. Please contact the publisher regarding any further use of this work. Publisher contact information may be obtained at

ABSTRACT The spike trains that transmit information between neurons are stochastic. We used the theory of random point processes and simulation methods to investigate the influence of temporal correlation of synaptic input current on firing statistics. The theory accounts for two sources for temporal correlation: synchrony between spikes in presynaptic input trains and the unitary synaptic current time course. Simulations show that slow temporal correlation of synaptic input leads to high variability in firing. In a leaky integrate-and-fire neuron model with spike afterhyperpolarization the theory accurately predicts the firing rate when the spike threshold is higher than two standard deviations of the membrane potential fluctuations. For lower thresholds the spike afterhyperpolarization reduces the firing rate below the theory’s predicted level when the synaptic correlation decays rapidly. If the synaptic correlation decays slower than the spike afterhyperpolarization, spike bursts can occur during single broad peaks of input fluctuations, increasing the firing rate over the prediction. Spike bursts lead to a coefficient of variation for the interspike intervals that can exceed one, suggesting an explanation of high coefficient of variation for interspike intervals observed in vivo.

A novel approach to cortical modeling was introduced by Knight et al. (1996). In their presentation cortical dynamics is formulated in terms of in- teracting populations of neurons, a perspective that is in part motivated by modern cortical imaging (For a review see Sirovich and Kaplan (2002)). The approach

Abstract. We model spiking neurons in locus coeruleus (LC), a brain nucleus involved in modulating cognitive performance, and compare with recent experimental data. Extracellular recordings from LC of monkeys performing target detection and selective attention tasks show varying responses dependent on stimuli and performance accuracy. From membrane voltage and ion channel equations, we derive a phase oscillator model for LC neurons. Average spiking probabilities of a pool of cells over many trials are then computed via a probability density formulation. These show that: (1) Post-stimulus response is elevated in populations with lower spike rates; (2) Responses decay exponentially due to noise and variable pre-stimulus spike rates; and (3) Shorter stimuli preferentially cause depressed post-activation spiking. These results allow us to propose mechanisms for the different LC responses observed across behavioral and task conditions, and to make explicit the role of baseline firing rates and the duration of task-related inputs in determining LC response.

Expressions for the moments of the membrane potential in Stein's model with synaptic (excitatory and inhibitory) reversal potentials are presented. The mean and variance of the voltage trajectories are used to approximate the mean and standard deviation of the interspike intervals under certain conditions (Stein, 1967) that yield a small coefficient of variation. The limiting case of excitatory input only and negligible decay places a lower bound on the coefficient of variation of the interspike intervals.