Impact of adaptation currents on synchronization of coupled exponential integrate-and-fire neurons

The ability of spiking neurons to synchronize their activity in a network depends on the response behavior of these neurons as quantified by the phase response curve (PRC) and on coupling properties. The PRC characterizes the effects of transient inputs on spike timing and can be measured experimentally. Here we use the adaptive exponential integrate-and-fire (aEIF) neuron model to determine how subthreshold and spike-triggered slow adaptation currents shape the PRC. Based on that, we predict how synchrony and phase locked states of coupled neurons change in presence of synaptic delays and unequal coupling strengths. We find that increased subthreshold adaptation currents cause a transition of the PRC from only phase advances to phase advances and delays in response to excitatory perturbations. Increased spike-triggered adaptation currents on the other hand predominantly skew the PRC to the right. Both adaptation induced changes of the PRC are modulated by spike frequency, being more prominent at lower frequencies. Applying phase reduction theory, we show that subthreshold adaptation stabilizes synchrony for pairs of coupled excitatory neurons, while spike-triggered adaptation causes locking with a small phase difference, as long as synaptic heterogeneities are negligible. For inhibitory pairs synchrony is stable and robust against conduction delays, and adaptation can mediate bistability of in-phase and anti-phase locking. We further demonstrate that stable synchrony and bistable in/anti-phase locking of pairs carry over to synchronization and clustering of larger networks. The effects of adaptation in aEIF neurons on PRCs and network dynamics qualitatively reflect those of biophysical adaptation currents in detailed Hodgkin-Huxley-based neurons, which underscores the utility of the aEIF model for investigating the dynamical behavior of networks. Our results suggest neuronal spike frequency adaptation as a mechanism synchronizing low frequency oscillations in local excitatory networks, but indicate that inhibition rather than excitation generates coherent rhythms at higher frequencies.     

Balanced neurons: analysis of leaky integrate-and-fire neurons with reversal potentials

 A new technique is presented for analyzing leaky integrate-and-fire neurons that incorporates reversal potentials, which impose a biologically realistic lower bound to the membrane potential. The time distribution of the synaptic inputs is modeled as a Poisson process. The analysis is carried out in the Gaussian approximation, which comparison with numerical simulations confirms is most accurate in the limit of a large number of inputs. The hypothesis that the observed variability in the spike times of cortical neurons is caused by a balance of excitatory and inhibitory synaptic inputs is supported by the results for the coefficient of variation of the interspike intervals. Its value decreases with both increasing numbers and amplitude of inputs, and is consistently lower than 1.0 over a wide range of realistic parameter values. The dependence of the output spike rate upon the rate, number, and amplitude of the synaptic inputs, as well as upon the value of the inhibitory reversal potential, is given. Received: 15 February 2001 / Accepted in revised form: 27 March 2001     

Impulse pattern in bi-directionally coupled model neurons of different dynamics

The effects of bi-directional gap junction coupling of two model neurons with subthreshold oscillations have been examined when the individual neurons are operating at different dynamical states either in the tonic or bursting firing mode. Our simulations indicate that intermediate coupling strengths mostly lead to highly variable, often chaotic impulse patterns whereas transition to completely synchronized activity at high coupling strengths is generally going along with transitions to regular limit cycle activity. The synchronized activity pattern, however, can be completely different from the original pattern of the uncoupled neurons.     

Population density models of integrate-and-fire neurons with jumps: well-posedness

In this paper we study the well-posedness of different models of population of leaky integrate-and-fire neurons with a population density approach. The synaptic interaction between neurons is modeled by a potential jump at the reception of a spike. We study populations that are self excitatory or self inhibitory. We distinguish the cases where this interaction is instantaneous from the one where there is a repartition of conduction delays. In the case of a bounded density of delays both excitatory and inhibitory population models are shown to be well-posed. But without conduction delay the solution of the model of self excitatory neurons may blow up. We analyze the different behaviours of the model with jumps compared to its diffusion approximation.     

Spike-train spectra and network response functions for non-linear integrate-and-fire neurons

Reduced models have long been used as a tool for the analysis of the complex activity taking place in neurons and their coupled networks. Recent advances in experimental and theoretical techniques have further demonstrated the usefulness of this approach. Despite the often gross simplification of the underlying biophysical properties, reduced models can still present significant difficulties in their analysis, with the majority of exact and perturbative results available only for the leaky integrate-and-fire model. Here an elementary numerical scheme is demonstrated which can be used to calculate a number of biologically important properties of the general class of non-linear integrate-and-fire models. Exact results for the first-passage-time density and spike-train spectrum are derived, as well as the linear response properties and emergent states of recurrent networks. Given that the exponential integrate-fire model has recently been shown to agree closely with the experimentally measured response of pyramidal cells, the methodology presented here promises to provide a convenient tool to facilitate the analysis of cortical-network dynamics.     

Exact analytical results for integrate-and-fire neurons driven by excitatory shot noise

A neuron receives input from other neurons via electrical pulses, so-called spikes. The pulse-like nature of the input is frequently neglected in analytical studies; instead, the input is usually approximated to be Gaussian. Recent experimental studies have shown, however, that an assumption underlying this approximation is often not met: Individual presynaptic spikes can have a significant effect on a neuron’s dynamics. It is thus desirable to explicitly account for the pulse-like nature of neural input, i.e. consider neurons driven by a shot noise – a long-standing problem that is mathematically challenging. In this work, we exploit the fact that excitatory shot noise with exponentially distributed weights can be obtained as a limit case of dichotomous noise, a Markovian two-state process. This allows us to obtain novel exact expressions for the stationary voltage density and the moments of the interspike-interval density of general integrate-and-fire neurons driven by such an input. For the special case of leaky integrate-and-fire neurons, we also give expressions for the power spectrum and the linear response to a signal. We verify and illustrate our expressions by comparison to simulations of leaky-, quadratic- and exponential integrate-and-fire neurons.     

Traveling waves and pulses in a one-dimensional network of excitable integrate-and-fire neurons

 We study the existence and stability of traveling waves and pulses in a one-dimensional network of integrate-and-fire neurons with synaptic coupling. This provides a simple model of excitable neural tissue. We first derive a self-consistency condition for the existence of traveling waves, which generates a dispersion relation between velocity and wavelength. We use this to investigate how wave-propagation depends on various parameters that characterize neuronal interactions such as synaptic and axonal delays, and the passive membrane properties of dendritic cables. We also establish that excitable networks support the propagation of solitary pulses in the long-wavelength limit. We then derive a general condition for the (local) asymptotic stability of traveling waves in terms of the characteristic equation of the linearized firing time map, which takes the form of an integro-difference equation of infinite order. We use this to analyze the stability of solitary pulses in the long-wavelength limit. Solitary wave solutions are shown to come in pairs with the faster (slower) solution stable (unstable) in the case of zero axonal delays; for non-zero delays and fast synapses the stable wave can itself destabilize via a Hopf bifurcation. Received: 27 October 1998     

Fast inference of interactions in assemblies of stochastic integrate-and-fire neurons from spike recordings

We present two Bayesian procedures to infer the interactions and external currents in an assembly of stochastic integrate-and-fire neurons from the recording of their spiking activity. The first procedure is based on the exact calculation of the most likely time courses of the neuron membrane potentials conditioned by the recorded spikes, and is exact for a vanishing noise variance and for an instantaneous synaptic integration. The second procedure takes into account the presence of fluctuations around the most likely time courses of the potentials, and can deal with moderate noise levels. The running time of both procedures is proportional to the number S of spikes multiplied by the squared number N of neurons. The algorithms are validated on synthetic data generated by networks with known couplings and currents. We also reanalyze previously published recordings of the activity of the salamander retina (including from 32 to 40 neurons, and from 65,000 to 170,000 spikes). We study the dependence of the inferred interactions on the membrane leaking time; the differences and similarities with the classical cross-correlation analysis are discussed.     

The most likely voltage path and large deviations approximations for integrate-and-fire neurons

We develop theory and numerical methods for computing the most likely subthreshold voltage path of a noisy integrate-and-fire (IF) neuron, given observations of the neuron’s superthreshold spiking activity. This optimal voltage path satisfies a second-order ordinary differential (Euler-Lagrange) equation which may be solved analytically in a number of special cases, and which may be solved numerically in general via a simple “shooting” algorithm. Our results are applicable for both linear and nonlinear subthreshold dynamics, and in certain cases may be extended to correlated subthreshold noise sources. We also show how this optimal voltage may be used to obtain approximations to (1) the likelihood that an IF cell with a given set of parameters was responsible for the observed spike train; and (2) the instantaneous firing rate and interspike interval distribution of a given noisy IF cell. The latter probability approximations are based on the classical Freidlin-Wentzell theory of large deviations principles for stochastic differential equations. We close by comparing this most likely voltage path to the true observed subthreshold voltage trace in a case when intracellular voltage recordings are available in vitro. Action Editor: Peter Latham     

On the dynamics of a pair of coupled neurons subject to alternating input rates

We present a statistical analysis of the firing activity of two coupled neuronal units that interact according to a 'sending-receiving' model. The membrane potential's behavior of both units is described by the Stein equations under the additional assumption that the spikes released by the sending neuron constitute an extra excitation for the receiving one. We also assume the presence of an alternating behavior for the rates of inputs to the sending neuron. By means of ad hoc simulations, we obtain, and then discuss, some statistical results concerning the spike production times of the units within the subintervals of the alternating inputs, as well as the reaction times of the receiving neuron.     

An integrate-and-fire model for synchronized bursting in a network of cultured cortical neurons

It has been suggested that spontaneous synchronous neuronal activity is an essential step in the formation of functional networks in the central nervous system. The key features of this type of activity consist of bursts of action potentials with associated spikes of elevated cytoplasmic calcium. These features are also observed in networks of rat cortical neurons that have been formed in culture. Experimental studies of these cultured networks have led to several hypotheses for the mechanisms underlying the observed synchronized oscillations. In this paper, bursting integrate-and-fire type mathematical models for regular spiking (RS) and intrinsic bursting (IB) neurons are introduced and incorporated through a small-world connection scheme into a two-dimensional excitatory network similar to those in the cultured network. This computer model exhibits spontaneous synchronous activity through mechanisms similar to those hypothesized for the cultured experimental networks. Traces of the membrane potential and cytoplasmic calcium from the model closely match those obtained from experiments. We also consider the impact on network behavior of the IB neurons, the geometry and the small world connection scheme. Action Editor: David Golomb     

Different types of noise in leaky integrate-and-fire model of neuronal dynamics with discrete periodical input

Different variants of stochastic leaky integrate-and-fire model for the membrane depolarisation of neurons are investigated. The model is driven by a constant input and equidistant pulses of fixed amplitude. These two types of signal are considered under the influence of three types of noise: white noise, jitter on interpulse distance, and noise in the amplitude of pulses. The results of computational experiments demonstrate the enhancement of the signal by noise in subthreshold regime and deterioration of the signal if it is sufficiently strong to carry the information in absence of noise. Our study holds mainly to central neurons that process discrete pulses although an application in sensory system is also available.     

Experimental and computational determination of tRNA dynamics

As the molecular representation of the genetic code, tRNA plays a central role in the translational machinery where it interacts with several proteins and other RNAs during the course of protein synthesis. These interactions exploit the dynamic flexibility of tRNA. In this minireview, we discuss the effects of modified bases, ions, and proteins on tRNA structure and dynamics and the challenges of observing its motions over the cycle of translation.     

Predicting the synaptic information efficacy in cortical layer 5 pyramidal neurons using a minimal integrate-and-fire model

Synaptic information efficacy (SIE) is a statistical measure to quantify the efficacy of a synapse. It measures how much information is gained, on the average, about the output spike train of a postsynaptic neuron if the input spike train is known. It is a particularly appropriate measure for assessing the input–output relationship of neurons receiving dynamic stimuli. Here, we compare the SIE of simulated synaptic inputs measured experimentally in layer 5 cortical pyramidal neurons in vitro with the SIE computed from a minimal model constructed to fit the recorded data. We show that even with a simple model that is far from perfect in predicting the precise timing of the output spikes of the real neuron, the SIE can still be accurately predicted. This arises from the ability of the model to predict output spikes influenced by the input more accurately than those driven by the background current. This indicates that in this context, some spikes may be more important than others. Lastly we demonstrate another aspect where using mutual information could be beneficial in evaluating the quality of a model, by measuring the mutual information between the model’s output and the neuron’s output. The SIE, thus, could be a useful tool for assessing the quality of models of single neurons in preserving input–output relationship, a property that becomes crucial when we start connecting these reduced models to construct complex realistic neuronal networks.     

Modeling brain dynamics using computational neurogenetic approach

The paper introduces a novel computational approach to brain dynamics modeling that integrates dynamic gene–protein regulatory networks with a neural network model. Interaction of genes and proteins in neurons affects the dynamics of the whole neural network. Through tuning the gene–protein interaction network and the initial gene/protein expression values, different states of the neural network dynamics can be achieved. A generic computational neurogenetic model is introduced that implements this approach. It is illustrated by means of a simple neurogenetic model of a spiking neural network of the generation of local field potential. Our approach allows for investigation of how deleted or mutated genes can alter the dynamics of a model neural network. We conclude with the proposal how to extend this approach to model cognitive neurodynamics.

Sawfishes stealth revealed using computational fluid dynamics

Detailed computational fluid dynamics simulations for the rostrum of three species of sawfish (Pristidae) revealed that negligible turbulent flow is generated from all rostra during lateral swipe prey manipulation and swimming. These results suggest that sawfishes are effective stealth hunters that may not be detected by their teleost prey's lateral line sensory system during pursuits. Moreover, during lateral swipes, the rostra were found to induce little velocity into the surrounding fluid. Consistent with previous data of sawfish feeding behaviour, these data indicate that the rostrum is therefore unlikely to be used to stir up the bottom to uncover benthic prey. Whilst swimming with the rostrum inclined at a small angle to the horizontal, the coefficient of drag of the rostrum is relatively low and the coefficient of lift is zero.