PubMed Health⌕ Search

Biomedical subjects

Nicolas Brunel

Publications and source records attributed to Nicolas Brunel.

15 recordsLinked to original sources

How noise affects the synchronization properties of recurrent networks of inhibitory neurons.

GABAergic interneurons play a major role in the emergence of various types of synchronous oscillatory patterns of activity in the central nervous system. Motivated by these experimental facts, modeling studies have investigated mechanisms for the emergence of coherent activity in networks of inhibitory neurons. However, most of these studies have focused either when the noise in the network is absent or weak or in the opposite situation when it is strong. Hence, a full picture of how noise affects the dynamics of such systems is still lacking. The aim of this letter is to provide a more comprehensive understanding of the mechanisms by which the asynchronous states in large, fully connected networks of inhibitory neurons are destabilized as a function of the noise level. Three types of single neuron models are considered: the leaky integrate-and-fire (LIF) model, the exponential integrate-and-fire (EIF), model and conductance-based models involving sodium and potassium Hodgkin-Huxley (HH) currents. We show that in all models, the instabilities of the asynchronous state can be classified in two classes. The first one consists of clustering instabilities, which exist in a restricted range of noise. These instabilities lead to synchronous patterns in which the population of neurons is broken into clusters of synchronously firing neurons. The irregularity of the firing patterns of the neurons is weak. The second class of instabilities, termed oscillatory firing rate instabilities, exists at any value of noise. They lead to cluster state at low noise. As the noise is increased, the instability occurs at larger coupling, and the pattern of firing that emerges becomes more irregular. In the regime of high noise and strong coupling, these instabilities lead to stochastic oscillations in which neurons fire in an approximately Poisson way with a common instantaneous probability of firing that oscillates in time.

Action Potentials↗

Contributions of intrinsic membrane dynamics to fast network oscillations with irregular neuronal discharges.

During fast oscillations in the local field potential (40-100 Hz gamma, 100-200 Hz sharp-wave ripples) single cortical neurons typically fire irregularly at rates that are much lower than the oscillation frequency. Recent computational studies have provided a mathematical description of such fast oscillations, using the leaky integrate-and-fire (LIF) neuron model. Here, we extend this theoretical framework to populations of more realistic Hodgkin-Huxley-type conductance-based neurons. In a noisy network of GABAergic neurons that are connected randomly and sparsely by chemical synapses, coherent oscillations emerge with a frequency that depends sensitively on the single cell's membrane dynamics. The population frequency can be predicted analytically from the synaptic time constants and the preferred phase of discharge during the oscillatory cycle of a single cell subjected to noisy sinusoidal input. The latter depends significantly on the single cell's membrane properties and can be understood in the context of the simplified exponential integrate-and-fire (EIF) neuron. We find that 200-Hz oscillations can be generated, provided the effective input conductance of single cells is large, so that the single neuron's phase shift is sufficiently small. In a two-population network of excitatory pyramidal cells and inhibitory neurons, recurrent excitation can either decrease or increase the population rhythmic frequency, depending on whether in a neuron the excitatory synaptic current follows or precedes the inhibitory synaptic current in an oscillatory cycle. Detailed single-cell properties have a substantial impact on population oscillations, even though rhythmicity does not originate from pacemaker neurons and is an emergent network phenomenon.

Action Potentials↗

Role of delays in shaping spatiotemporal dynamics of neuronal activity in large networks.

We study the effect of delays on the dynamics of large networks of neurons. We show that delays give rise to a wealth of bifurcations and to a rich phase diagram, which includes oscillatory bumps, traveling waves, lurching waves, standing waves arising via a period-doubling bifurcation, aperiodic regimes, and regimes of multistability. We study the existence and the stability of the various dynamical patterns analytically and numerically in a simplified rate model as a function of the interaction parameters. The results derived in that framework allow us to understand the origin of the diversity of dynamical states observed in large networks of spiking neurons.

Cell Communication↗

Dynamics of the instantaneous firing rate in response to changes in input statistics.

We review and extend recent results on the instantaneous firing rate dynamics of simplified models of spiking neurons in response to noisy current inputs. It has been shown recently that the response of the instantaneous firing rate to small amplitude oscillations in the mean inputs depends in the large frequency limit f on the spike initiation dynamics. A particular simplified model, the exponential integrate-and-fire (EIF) model, has a response that decays as 1/f in the large frequency limit and describes very well the response of conductance-based models with a Hodgkin-Huxley type fast sodium current. Here, we show that the response of the EIF instantaneous firing rate also decays as 1/f in the case of an oscillation in the variance of the inputs for both white and colored noise. We then compute the initial transient response of the firing rate of the EIF model to a step change in its mean inputs and/or in the variance of its inputs. We show that in both cases the response speed is proportional to the neuron stationary firing rate and inversely proportional to a 'spike slope factor' Delta(T) that controls the sharpness of spike initiation: as 1/Delta(T) for a step change in mean inputs, and as 1/Delta(T) (2) for a step change in the variance in the inputs.

Action Potentials↗

Optimal information storage and the distribution of synaptic weights: perceptron versus Purkinje cell.

It is widely believed that synaptic modifications underlie learning and memory. However, few studies have examined what can be deduced about the learning process from the distribution of synaptic weights. We analyze the perceptron, a prototypical feedforward neural network, and obtain the optimal synaptic weight distribution for a perceptron with excitatory synapses. It contains more than 50% silent synapses, and this fraction increases with storage reliability: silent synapses are therefore a necessary byproduct of optimizing learning and reliability. Exploiting the classical analogy between the perceptron and the cerebellar Purkinje cell, we fitted the optimal weight distribution to that measured for granule cell-Purkinje cell synapses. The two distributions agreed well, suggesting that the Purkinje cell can learn up to 5 kilobytes of information, in the form of 40,000 input-output associations.

Afferent Pathways↗

How spike generation mechanisms determine the neuronal response to fluctuating inputs.

This study examines the ability of neurons to track temporally varying inputs, namely by investigating how the instantaneous firing rate of a neuron is modulated by a noisy input with a small sinusoidal component with frequency (f). Using numerical simulations of conductance-based neurons and analytical calculations of one-variable nonlinear integrate-and-fire neurons, we characterized the dependence of this modulation on f. For sufficiently high noise, the neuron acts as a low-pass filter. The modulation amplitude is approximately constant for frequencies up to a cutoff frequency, fc, after which it decays. The cutoff frequency increases almost linearly with the firing rate. For higher frequencies, the modulation amplitude decays as C/falpha, where the power alpha depends on the spike initiation mechanism. For conductance-based models, alpha = 1, and the prefactor C depends solely on the average firing rate and a spike "slope factor," which determines the sharpness of the spike initiation. These results are attributable to the fact that near threshold, the sodium activation variable can be approximated by an exponential function. Using this feature, we propose a simplified one-variable model, the "exponential integrate-and-fire neuron," as an approximation of a conductance-based model. We show that this model reproduces the dynamics of a simple conductance-based model extremely well. Our study shows how an intrinsic neuronal property (the characteristics of fast sodium channels) determines the speed with which neurons can track changes in input.

Action Potentials↗

Firing-rate resonance in a generalized integrate-and-fire neuron with subthreshold resonance.

Neurons that exhibit a peak at finite frequency in their membrane potential response to oscillatory inputs are widespread in the nervous system. However, the influence of this subthreshold resonance on spiking properties has not yet been thoroughly analyzed. To this end, generalized integrate-and-fire models are introduced that reproduce at the linear level the subthreshold behavior of any given conductance-based model. A detailed analysis is presented of the simplest resonant model of this kind that has two variables: the membrane potential and a supplementary voltage-gated resonant variable. The firing-rate modulation created by a noisy weak oscillatory drive, mimicking an in vivo environment, is computed numerically and analytically when the dynamics of the resonant variable is slow compared to that of the membrane potential. The results show that the firing-rate modulation is shaped by the subthreshold resonance. For weak noise, the firing-rate modulation has a minimum near the preferred subthreshold frequency. For higher noise, such as that prevailing in vivo, the firing-rate modulation peaks near the preferred subthreshold frequency.

Action Potentials↗

What determines the frequency of fast network oscillations with irregular neural discharges? I. Synaptic dynamics and excitation-inhibition balance.

When the local field potential of a cortical network displays coherent fast oscillations ( approximately 40-Hz gamma or approximately 200-Hz sharp-wave ripples), the spike trains of constituent neurons are typically irregular and sparse. The dichotomy between rhythmic local field and stochastic spike trains presents a challenge to the theory of brain rhythms in the framework of coupled oscillators. Previous studies have shown that when noise is large and recurrent inhibition is strong, a coherent network rhythm can be generated while single neurons fire intermittently at low rates compared to the frequency of the oscillation. However, these studies used too simplified synaptic kinetics to allow quantitative predictions of the population rhythmic frequency. Here we show how to derive quantitatively the coherent oscillation frequency for a randomly connected network of leaky integrate-and-fire neurons with realistic synaptic parameters. In a noise-dominated interneuronal network, the oscillation frequency depends much more on the shortest synaptic time constants (delay and rise time) than on the longer synaptic decay time, and approximately 200-Hz frequency can be realized with synaptic time constants taken from slice data. In a network composed of both interneurons and excitatory cells, the rhythmogenesis is a compromise between two scenarios: the fast purely interneuronal mechanism, and the slower feedback mechanism (relying on the excitatory-inhibitory loop). The properties of the rhythm are determined essentially by the ratio of time scales of excitatory and inhibitory currents and by the balance between the mean recurrent excitation and inhibition. Faster excitation than inhibition, or a higher excitation/inhibition ratio, favors the feedback loop and a much slower oscillation (typically in the gamma range).

Action Potentials↗

Retrospective and prospective persistent activity induced by Hebbian learning in a recurrent cortical network.

Recordings from cells in the associative cortex of monkeys performing visual working memory tasks link persistent neuronal activity, long-term memory and associative memory. In particular, delayed pair-associate tasks have revealed neuronal correlates of long-term memory of associations between stimuli. Here, a recurrent cortical network model with Hebbian plastic synapses is subjected to the pair-associate protocol. In a first stage, learning leads to the appearance of delay activity, representing individual images ('retrospective' activity). As learning proceeds, the same learning mechanism uses retrospective delay activity together with choice stimulus activity to potentiate synapses connecting neural populations representing associated images. As a result, the neural population corresponding to the pair-associate of the image presented is activated prior to its visual stimulation ('prospective' activity). The probability of appearance of prospective activity is governed by the strength of the inter-population connections, which in turn depends on the frequency of pairings during training. The time course of the transitions from retrospective to prospective activity during the delay period is found to depend on the fraction of slow, N-methyl-d-aspartate-like receptors at excitatory synapses. For fast recurrent excitation, transitions are abrupt; slow recurrent excitation renders transitions gradual. Both scenarios lead to a gradual rise of delay activity when averaged over many trials, because of the stochastic nature of the transitions. The model reproduces most of the neuro-physiological data obtained during such tasks, makes experimentally testable predictions and demonstrates how persistent activity (working memory) brings about the learning of long-term associations.

Animals↗

Dynamics and plasticity of stimulus-selective persistent activity in cortical network models.

Persistent neuronal activity is widespread in many areas of the cerebral cortex of monkeys performing cognitive tasks with a working memory component. Modeling studies have helped understanding of the conditions under which persistent activity can be sustained in cortical circuits. Here, we first review several basic models of persistent activity, including bistable models with excitation only and multistable models for working memory of a discrete set of pictures or objects with structured excitation and global inhibition. In many experiments, persistent activity has been shown to be subject to changes due to associative learning. In cortical network models, Hebbian learning shapes the synaptic structure and, in turn, the properties of persistent activity when pictures are associated together in the course of a task. It is shown how the theoretical models can reproduce basic experimental findings of neurophysiological recordings from inferior temporal and perirhinal cortices obtained using the following experimental protocols: (i) the pair-associate task; (ii) the pair-associate task with color switch; and (iii) the delay match to sample task with a fixed sequence of samples.

Animals↗

Firing rate of the noisy quadratic integrate-and-fire neuron.

We calculate the firing rate of the quadratic integrate-and-fire neuron in response to a colored noise input current. Such an input current is a good approximation to the noise due to the random bombardment of spikes, with the correlation time of the noise corresponding to the decay time of the synapses. The key parameter that determines the firing rate is the ratio of the correlation time of the colored noise, tau(s), to the neuronal time constant, tau(m). We calculate the firing rate exactly in two limits: when the ratio, tau(s)/tau(m), goes to zero (white noise) and when it goes to infinity. The correction to the short correlation time limit is O(tau(s)/tau(m)), which is qualita tively different from that of the leaky integrate-and-fire neuron, where the correction is O( radical tau(s)/tau(m)). The difference is due to the different boundary conditions of the probability density function of the membrane potential of the neuron at firing threshold. The correction to the long correlation time limit is O(tau(m)/tau(s)). By combining the short and long correlation time limits, we derive an expression that provides a good approximation to the firing rate over the whole range of tau(s)/tau(m) in the suprathreshold regime-that is, in a regime in which the average current is sufficient to make the cell fire. In the subthreshold regime, the expression breaks down somewhat when tau(s) becomes large compared to tau(m).

Action Potentials↗

From subthreshold to firing-rate resonance.

Many types of neurons exhibit subthreshold resonance. However, little is known about whether this frequency preference influences spike emission. Here, the link between subthreshold resonance and firing rate is examined in the framework of conductance-based models. A classification of the subthreshold properties of a general class of neurons is first provided. In particular, a class of neurons is identified in which the input impedance exhibits a suppression at a nonzero low frequency as well as a peak at higher frequency. The analysis is then extended to the effect of subthreshold resonance on the dynamics of the firing rate. The considered input current comprises a background noise term, mimicking the massive synaptic bombardment in vivo. Of interest is the modulatory effect an additional weak oscillating current has on the instantaneous firing rate. When the noise is weak and firing regular, the frequency most preferentially modulated is the firing rate itself. Conversely, when the noise is strong and firing irregular, the modulation is strongest at the subthreshold resonance frequency. These results are demonstrated for two specific conductance-based models and for a generalization of the integrate-and-fire model that captures subthreshold resonance. They suggest that resonant neurons are able to communicate their frequency preference to postsynaptic targets when the level of noise is comparable to that prevailing in vivo.

Algorithms↗

Dynamics of the firing probability of noisy integrate-and-fire neurons.

Cortical neurons in vivo undergo a continuous bombardment due to synaptic activity, which acts as a major source of noise. Here, we investigate the effects of the noise filtering by synapses with various levels of realism on integrate-and-fire neuron dynamics. The noise input is modeled by white (for instantaneous synapses) or colored (for synapses with a finite relaxation time) noise. Analytical results for the modulation of firing probability in response to an oscillatory input current are obtained by expanding a Fokker-Planck equation for small parameters of the problem - when both the amplitude of the modulation is small compared to the background firing rate and the synaptic time constant is small compared to the membrane time constant. We report here the detailed calculations showing that if a synaptic decay time constant is included in the synaptic current model, the firing-rate modulation of the neuron due to an oscillatory input remains finite in the high-frequency limit with no phase lag. In addition, we characterize the low-frequency behavior and the behavior of the high-frequency limit for intermediate decay times. We also characterize the effects of introducing a rise time to the synaptic currents and the presence of several synaptic receptors with different kinetics. In both cases, we determine, using numerical simulations, an effective decay time constant that describes the neuronal response completely.

Action Potentials↗

A continuous attractor network model without recurrent excitation: maintenance and integration in the head direction cell system.

Motivated by experimental observations of the head direction system, we study a three population network model that operates as a continuous attractor network. This network is able to store in a short-term memory an angular variable (the head direction) as a spatial profile of activity across neurons in the absence of selective external inputs, and to accurately update this variable on the basis of angular velocity inputs. The network is composed of one excitatory population and two inhibitory populations, with inter-connections between populations but no connections within the neurons of a same population. In particular, there are no excitatory-to-excitatory connections. Angular velocity signals are represented as inputs in one inhibitory population (clockwise turns) or the other (counterclockwise turns). The system is studied using a combination of analytical and numerical methods. Analysis of a simplified model composed of threshold-linear neurons gives the conditions on the connectivity for (i) the emergence of the spatially selective profile, (ii) reliable integration of angular velocity inputs, and (iii) the range of angular velocities that can be accurately integrated by the model. Numerical simulations allow us to study the proposed scenario in a large network of spiking neurons and compare their dynamics with that of head direction cells recorded in the rat limbic system. In particular, we show that the directional representation encoded by the attractor network can be rapidly updated by external cues, consistent with the very short update latencies observed experimentally by Zugaro et al. (2003) in thalamic head direction cells.

Action Potentials↗