The Morris–Lecar neuron model embeds a leaky integrate-and-fire model

We show that the stochastic Morris–Lecar neuron, in a neighborhood of its stable point, can be approximated by a two-dimensional Ornstein–Uhlenbeck (OU) modulation of a constant circular motion. The associated radial OU process is an example of a leaky integrate-and-fire (LIF) model prior to firing. A new model constructed from a radial OU process together with a simple firing mechanism based on detailed Morris–Lecar firing statistics reproduces the Morris–Lecar Interspike Interval (ISI) distribution, and has the computational advantages of a LIF. The result justifies the large amount of attention paid to the LIF models.     

Local interneurons and information processing in the olfactory glomeruli of the moth Manduca sexta

Intracellular recordings were made from the major neurites of local interneurons in the moth antennal lobe. Antennal nerve stimulation evoked 3 patterns of postsynaptic activity: (i) a short-latency compound excitatory postsynaptic potential that, based on electrical stimulation of the antennal nerve and stimulation of the antenna with odors, represents a monosynaptic input from olfactory afferent axons (71 out of 86 neurons), (ii) a delayed activation of firing in response to both electrical- and odor-driven input (11 neurons), and (iii) a delayed membrane hyperpolarization in response to antennal nerve input (4 neurons).Simultaneous intracellular recordings from a local interneuron with short-latency responses and a projection (output) neuron revealed unidirectional synaptic interactions between these two cell types. In 20% of the 30 pairs studied, spontaneous and current-induced spiking activity in a local interneuron correlated with hyperpolarization and suppression of firing in a projection neuron. No evidence for recurrent or feedback inhibition of projection neurons was found. Furthermore, suppression of firing in an inhibitory local interneuron led to an increase in firing in the normally quiescent projection neuron, suggesting that a disinhibitory pathway may mediate excitation in projection neurons. This is the first direct evidence of an inhibitory role for local interneurons in olfactory information processing in insects. Through different types of multisynaptic interactions with projection neurons, local interneurons help to generate and shape the output from olfactory glomeruli in the antennal lobe.Abbreviations AL antennal lobe - EPSP excitatory postsynaptic potential - GABA -aminobutyric acid - IPSP inhibitory postsynaptic potential - LN local interneuron - MGC macroglomerular complex - OB olfactory bulb - PN projection neuron - TES N-tris[hydroxymethyl]methyl-2-aminoethane-sulfonic acid     

Synchronization and sensitivity enhancement of the Hodgkin-Huxley neurons due to inhibitory inputs

Recent experimental results imply that inhibitory postsynaptic potentials can play a functional role in realizing synchronization of neuronal firing in the brain. In order to examine the relation between inhibition and synchronous firing of neurons theoretically, we analyze possible effects of synchronization and sensitivity enhancement caused by inhibitory inputs to neurons with a biologically realistic model of the Hodgkin-Huxley equations. The result shows that, after an inhibitory spike, the firing probability of a single postsynaptic neuron exposed to random excitatory background activity oscillates with time. The oscillation of the firing probability can be related to synchronous firing of neurons receiving an inhibitory spike simultaneously. Further, we show that when an inhibitory spike input precedes an excitatory spike input, the presence of such preceding inhibition raises the firing probability peak of the neuron after the excitatory input. The result indicates that an inhibitory spike input can enhance the sensitivity of the postsynaptic neuron to the following excitatory spike input. Two neural network models based on these effects on postsynaptic neurons caused by inhibitory inputs are proposed to demonstrate possible mechanisms of detecting particular spatiotemporal spike patterns. Received: 15 April 1999 /Accepted in revised form: 25 November 1999     

Noise-induced divisive gain control in neuron models

A recent computational study of gain control via shunting inhibition has shown that the slope of the frequency-versus-input (f-I) characteristic of a neuron can be decreased by increasing the noise associated with the inhibitory input (Neural Comput. 13, 227-248). This novel noise-induced divisive gain control relies on the concommittant increase of the noise variance with the mean of the total inhibitory conductance. Here we investigate this effect using different neuronal models. The effect is shown to occur in the standard leaky integrate-and-fire (LIF) model with additive Gaussian white noise, and in the LIF with multiplicative noise acting on the inhibitory conductance. The noisy scaling of input currents is also shown to occur in the one-dimensional theta-neuron model, which has firing dynamics, as well as a large scale compartmental model of a pyramidal cell in the electrosensory lateral line lobe of a weakly electric fish. In this latter case, both the inhibition and the excitatory input have Poisson statistics; noise-induced divisive inhibition is thus seen in f-I curves for which the noise increases along with the input I. We discuss how the variation of the noise intensity along with inputs is constrained by the physiological context and the class of model used, and further provide a comparison of the divisive effect across models.     

Near Poisson-type firing produced by concurrent excitation and inhibition

The effect of inhibition on the firing variability is examined in this paper using the biologically-inspired temporal noisy-leaky integrator (TNLI) neuron model. The TNLI incorporates hyperpolarising inhibition with negative current pulses of controlled shapes and it also separates dendritic from somatic integration. The firing variability is observed by looking at the coefficient of variation (C(V)) (standard deviation/mean interspike interval) as a function of the mean interspike interval of firing (delta tM) and by comparing the results with the theoretical curve for random spike trains, as well as looking at the interspike interval (ISI) histogram distributions. The results show that with 80% inhibition, firing at high rates (up to 200 Hz) is nearly consistent with a Poisson-type variability, which complies with the analysis of cortical neuron firing recordings by Softky and Koch [1993, J. Neurosci. 13(1) 334-530]. We also demonstrate that the mechanism by which inhibition increases the C(V) values is by introducing more short intervals in the firing pattern as indicated by a small initial hump at the beginning of the ISI histogram distribution. The use of stochastic inputs and the separation of the dendritic and somatic integration which we model in TNLI, also affect the high firing, near Poisson-type (explained in the paper) variability produced. We have also found that partial dendritic reset increases slightly the firing variability especially at short ISIs.     

Multisynaptic activity in a pyramidal neuron model and neural code

The highly irregular firing of mammalian cortical pyramidal neurons is one of the most striking observation of the brain activity. This result affects greatly the discussion on the neural code, i.e. how the brain codes information transmitted along the different cortical stages. In fact it seems to be in favor of one of the two main hypotheses about this issue, named the rate code. But the supporters of the contrasting hypothesis, the temporal code, consider this evidence inconclusive. We discuss here a leaky integrate-and-fire model of a hippocampal pyramidal neuron intended to be biologically sound to investigate the genesis of the irregular pyramidal firing and to give useful information about the coding problem. To this aim, the complete set of excitatory and inhibitory synapses impinging on such a neuron has been taken into account. The firing activity of the neuron model has been studied by computer simulation both in basic conditions and allowing brief periods of over-stimulation in specific regions of its synaptic constellation. Our results show neuronal firing conditions similar to those observed in experimental investigations on pyramidal cortical neurons. In particular, the variation coefficient (CV) computed from the inter-spike intervals (ISIs) in our simulations for basic conditions is close to the unity as that computed from experimental data. Our simulation shows also different behaviors in firing sequences for different frequencies of stimulation.     

Spike tranfer in statistical neuron assemblies. II. Neuron--nonlinear spike source]

The mathematical model of the neuron function is known to rely on space summing of excitement. The spikes contribute to the inner state of the neuron the farther from cell soma the synapses are located. The difference between excitatory and inhibitory effect results in spike firing if only neural firing threshold is achieved. The values of spike flux have been estimated on the basis of the model of CA3 sector of the Hippocampus and were found to be 15 divided by 35 imp/s.     

Cortical network modeling: analytical methods for firing rates and some properties of networks of LIF neurons.

The circuitry of cortical networks involves interacting populations of excitatory (E) and inhibitory (I) neurons whose relationships are now known to a large extent. Inputs to E- and I-cells may have their origins in remote or local cortical areas. We consider a rudimentary model involving E- and I-cells. One of our goals is to test an analytic approach to finding firing rates in neural networks without using a diffusion approximation and to this end we consider in detail networks of excitatory neurons with leaky integrate-and-fire (LIF) dynamics. A simple measure of synchronization, denoted by S(q), where q is between 0 and 100 is introduced. Fully connected E-networks have a large tendency to become dominated by synchronously firing groups of cells, except when inputs are relatively weak. We observed random or asynchronous firing in such networks with diverse sets of parameter values. When such firing patterns were found, the analytical approach was often able to accurately predict average neuronal firing rates. We also considered several properties of E-E networks, distinguishing several kinds of firing pattern. Included were those with silences before or after periods of intense activity or with periodic synchronization. We investigated the occurrence of synchronized firing with respect to changes in the internal excitatory postsynaptic potential (EPSP) magnitude in a network of 100 neurons with fixed values of the remaining parameters. When the internal EPSP size was less than a certain value, synchronization was absent. The amount of synchronization then increased slowly as the EPSP amplitude increased until at a particular EPSP size the amount of synchronization abruptly increased, with S(5) attaining the maximum value of 100%. We also found network frequency transfer characteristics for various network sizes and found a linear dependence of firing frequency over wide ranges of the external afferent frequency, with non-linear effects at lower input frequencies. The theory may also be applied to sparsely connected networks, whose firing behaviour was found to change abruptly as the probability of a connection passed through a critical value. The analytical method was also found to be useful for a feed-forward excitatory network and a network of excitatory and inhibitory neurons.     

Attentional modulation of firing rate and synchrony in a model cortical network

The response of a neuron in the visual cortex to stimuli of different contrast placed in its receptive field is commonly characterized using the contrast response curve. When attention is directed into the receptive field of a V4 neuron, its contrast response curve is shifted to lower contrast values (Reynolds et al., 2000). The neuron will thus be able to respond to weaker stimuli than it responded to without attention. Attention also increases the coherence between neurons responding to the same stimulus (Fries et al., 2001). We studied how the firing rate and synchrony of a densely interconnected cortical network varied with contrast and how they were modulated by attention. The changes in contrast and attention were modeled as changes in driving current to the network neurons. We found that an increased driving current to the excitatory neurons increased the overall firing rate of the network, whereas variation of the driving current to inhibitory neurons modulated the synchrony of the network. We explain the synchrony modulation in terms of a locking phenomenon during which the ratio of excitatory to inhibitory firing rates is approximately constant for a range of driving current values. We explored the hypothesis that contrast is represented primarily as a drive to the excitatory neurons, whereas attention corresponds to a reduction in driving current to the inhibitory neurons. Using this hypothesis, the model reproduces the following experimental observations: (1) the firing rate of the excitatory neurons increases with contrast; (2) for high contrast stimuli, the firing rate saturates and the network synchronizes; (3) attention shifts the contrast response curve to lower contrast values; (4) attention leads to stronger synchronization that starts at a lower value of the contrast compared with the attend-away condition. In addition, it predicts that attention increases the delay between the inhibitory and excitatory synchronous volleys produced by the network, allowing the stimulus to recruit more downstream neurons. Action Editor: David Golomb     

Firing rate of a retinal neuron are not predictable from interspike interval statistics.

The intervals between successive action potentials (impulses, or "spikes") produced the maintained firing of a neuron (ISIs) are often treated as if they were independent on each other; that is, an impulse train is considered as a stationary renewal process. If this is so, the variability of the mean rate of firing impulses in a sequence of temporal windows should be predictable from the distribution of ISIs. This was found not to be the case for the maintained firing of retinal ganglion cells in goldfish. Although some evident nonstationarity sometimes resulted in greater variability of the observed rate distributions than those predicted (for relatively long temporal windows), as a general rule the observed rate distributions were considerable less dispersed than would be predicted by sampling of the ISI distributions. This was taken as evidence of long-term serial dependency between successive ISIs; however, two standard test for dependency (autocorrelations and serial correlograms failed to to reveal structure of sufficiently long duration to account for the effect noted.     

Stochastic automaton models for interaction of excitatory and inhibitory impulse sequences in neurons

The mechanisms by which excitatory and inhibitory input impulse sequences interact in changing the spike probability in neurons are examined in the two mathematical neuron models; one is a real-time neuron model which is close to physiological reality, and the other a stochastic automaton model for the temporal pattern discrimination proposed in the previous paper (Tsukada et al., 1976), which is developed in this paper as neuron models for interaction of excitatory and inhibitory input impulse sequences. The interval distributions of the output spike train from these models tend to be multimodal and are compared with those used for experimental data, reported by Bishop et al. (1964) for geniculate neuron activity and Poisson process deleting model analyzed by Ten Hoopen et al. (1966). Special attention, moreover, should be paid to how different forms of inhibitory input are transformed into the output interval distributions through these neuron models. These results exhibit a clear correlation between inhibitory input form and output interval distribution. More detailed information on this mechanism is obtained from the computations of recurrence-time under the stationary condition to go from active state to itself for the first time, each of which is influenced by the inhibitory input forms. In addition to these facts, some resultant characteristics on interval histogram and serial correlation are discussed in relation to physiological data from the literature.     

Modelling Feedback Excitation,Pacemaker Properties and Sensory Switching of Electrically Coupled Brainstem Neurons Controlling Rhythmic Activity

What cellular and network properties allow reliable neuronal rhythm generation or firing that can be started and stopped by brief synaptic inputs? We investigate rhythmic activity in an electrically-coupled population of brainstem neurons driving swimming locomotion in young frog tadpoles, and how activity is switched on and off by brief sensory stimulation. We build a computational model of 30 electrically-coupled conditional pacemaker neurons on one side of the tadpole hindbrain and spinal cord. Based on experimental estimates for neuron properties, population sizes, synapse strengths and connections, we show that: long-lasting, mutual, glutamatergic excitation between the neurons allows the network to sustain rhythmic pacemaker firing at swimming frequencies following brief synaptic excitation; activity persists but rhythm breaks down without electrical coupling; NMDA voltage-dependency doubles the range of synaptic feedback strengths generating sustained rhythm. The network can be switched on and off at short latency by brief synaptic excitation and inhibition. We demonstrate that a population of generic Hodgkin-Huxley type neurons coupled by glutamatergic excitatory feedback can generate sustained asynchronous firing switched on and off synaptically. We conclude that networks of neurons with NMDAR mediated feedback excitation can generate self-sustained activity following brief synaptic excitation. The frequency of activity is limited by the kinetics of the neuron membrane channels and can be stopped by brief inhibitory input. Network activity can be rhythmic at lower frequencies if the neurons are electrically coupled. Our key finding is that excitatory synaptic feedback within a population of neurons can produce switchable, stable, sustained firing without synaptic inhibition.     

Synaptic depression increases the selectivity of a neuron to its preferred pattern and binarizes the neural code

The preferred pattern of a neuron is defined here by the set of features detected by its excitatory inputs. It is shown that the Leaky integrate-and-fire (LIF) model of a neuron has a poor selectivity to its preferred pattern. Its response is determined by the total current injected by input spike trains. Thus, a few inputs with a high activity (an incomplete pattern) can elicit the same response as many inputs (a complete pattern) with a weak activity. A theoretical model of depressing synapse with linear recovery is proposed which eliminates this problem. Using this model, the time-averaged current injected in the soma by a spike train becomes independent on its frequency. The neural code thus becomes binary, and the response strength of the target neuron depends only on the number of active inputs. Simulations show that a biological model of strong synaptic depression has effects similar to those of the ideal linear model. The best selectivity is obtained with long somatic decay time constants (>50 ms) and with depression recovery time constants larger or equal to the somatic decay time constant. Thus, by eliminating information carried in the input firing rate, a neuron can improve its pattern recognition performance.     

A model for the variability of interspike intervals during sustained firing of a retinal neuron.

The statistics of the variability of interspike intervals of ganglion cells in the retina of goldfish are modeled by assuming the noise in an integrate-and-fire mechanism is proportional to the reciprocal of a normally distributed variable. This model meets the constraint that the coefficient of variation of the interspike. This does not change when the mean firing rate of the neuron changes. Alternative sources of variability of interspike intervals are discussed.     

Statistical Analysis of Membrane Potential Fluctuations: Relation with Presynaptic Spike Train

In a study of integration at the single neuron level, the relationships between the postsynaptic membrane potential and the presynaptic spike train were analyzed. Fluctuations in membrane potential of neurons in the visceral ganglion of Aplysia were measured and described by histograms. The histogram estimates the probability density function of the membrane potential. Comparisons were made among histograms when there was no synaptic input, and when there was a single input in which variations were made in the PSP (postsynaptic potential) sign, i.e. excitatory or inhibitory, and arrival statistics, e.g. slow or fast, regular, Poisson-like, or patterned. This was examined in cells where the membrane potential was constant and in cells in which there was spontaneous pacemaker activity. The form of the histogram depended on whether the neuron was spontaneously quiescent or a pacemaker, or whether it received presynaptic input and, if it did, on the sign and temporal characteristics of such input. From such histograms the mean firing rate of output spike trains can be predicted; additional information of a temporal nature is required, however, to predict features of the interval structure of the output train. Suggestions are made concerning the way the nervous system might utilize the information summarized in the membrane potential histogram.     

Recurrent networks of submucous neurons controlling intestinal secretion: a modeling study

Secretomotor neurons, immunoreactive for vasoactive intestinal peptide (VIP), are important in controlling chloride secretion in the small intestine. These neurons form functional synapses with other submucosal VIP neurons and transmit via slow excitatory postsynaptic potentials (EPSPs). Thus they form a recurrent network with positive feedback. Intrinsic sensory neurons within the submucosa are also likely to form recurrent networks with positive feedback, provide substantial output to VIP neurons, and receive input from VIP neurons. If positive feedback within recurrent networks is sufficiently large, then neurons in the network respond to even small stimuli by firing at their maximum possible rate, even after the stimulus is removed. However, it is not clear whether such a mechanism operates within the recurrent networks of submucous neurons. We investigated this question by performing computer simulations of realistic models of VIP and intrinsic sensory neuron networks. In the expected range of electrophysiological properties, we found that activity in the VIP neuron network decayed slowly after cessation of a stimulus, indicating that positive feedback is not strong enough to support the uncontrolled firing state. The addition of intrinsic sensory neurons produced a low stable firing rate consistent with the common finding that basal secretory activity is, in part, neurogenic. Changing electrophysiological properties enables these recurrent networks to support the uncontrolled firing state, which may have implications with hypersecretion in the presence of enterotoxins such as cholera-toxin.     

Neural spike generator entrainment by regularly spaced synaptic input

 It has been known for 30 years that the output of a repetitively firing neuron or pacemaker can be synchronized (locked) to regularly spaced inhibitory or excitatory postsynaptic input potentials. Conditions for stable locking have been determined mathematically, demonstrated in computer simulation, and locking has been observed in vivo. We have developed a neural spike generator circuit model which exhibits stable locking to externally derived simulated inhibitory or excitatory post-synaptic inputs. Conditions for stable 1 : 1 lock, in which pacemaker output frequency matches that of the periodic input, are derived. These take the form of expressions for stable delay and convergence factor which incorporate known or measurable parameters of the circuit model. The expressions have been evaluated and shown to compare satisfactorily with experimental observations of locking by our circuit model. Received: 28 March 1996 / Accepted in revised form: 18 February 1997     

Determinants of synaptic integration and heterogeneity in rebound firing explored with data-driven models of deep cerebellar nucleus cells

Significant inroads have been made to understand cerebellar cortical processing but neural coding at the output stage of the cerebellum in the deep cerebellar nuclei (DCN) remains poorly understood. The DCN are unlikely to just present a relay nucleus because Purkinje cell inhibition has to be turned into an excitatory output signal, and DCN neurons exhibit complex intrinsic properties. In particular, DCN neurons exhibit a range of rebound spiking properties following hyperpolarizing current injection, raising the question how this could contribute to signal processing in behaving animals. Computer modeling presents an ideal tool to investigate how intrinsic voltage-gated conductances in DCN neurons could generate the heterogeneous firing behavior observed, and what input conditions could result in rebound responses. To enable such an investigation we built a compartmental DCN neuron model with a full dendritic morphology and appropriate active conductances. We generated a good match of our simulations with DCN current clamp data we recorded in acute slices, including the heterogeneity in the rebound responses. We then examined how inhibitory and excitatory synaptic input interacted with these intrinsic conductances to control DCN firing. We found that the output spiking of the model reflected the ongoing balance of excitatory and inhibitory input rates and that changing the level of inhibition performed an additive operation. Rebound firing following strong Purkinje cell input bursts was also possible, but only if the chloride reversal potential was more negative than −70 mV to allow de-inactivation of rebound currents. Fast rebound bursts due to T-type calcium current and slow rebounds due to persistent sodium current could be differentially regulated by synaptic input, and the pattern of these rebounds was further influenced by HCN current. Our findings suggest that active properties of DCN neurons could play a crucial role for signal processing in the cerebellum.     

A massively connected subthalamic nucleus leads to the generation of widespread pulses.

A composite model of the subthalamic nucleus is developed from physiological and anatomical considerations. First, study of a geometric model of the anatomical arrangements of projection neurons within the nucleus indicates that they form a massively connected network. Second, given the excitatory nature of these neurons, their threshold and peak firing rates, a simple model of neuron responses reveals that large regions of this highly interconnected nucleus can respond to excitatory input in the form of a wide-spread uniform pulse. Such widespread pulses of activity may act as a braking signal that resets the major basal ganglia output nuclei.     

Effect of steroid milieu on gonadotropin-releasing hormone-1 neuron firing pattern and luteinizing hormone levels in male mice

GnRH neuronal function is regulated by gonadal hormone feedback. In males, testosterone can act directly or be converted to either dihydrotestosterone (DHT) or estradiol (E2). We examined central steroid feedback by recording firing of green fluorescent protein (GFP)-identified GnRH neurons in brain slices from male mice that were intact, castrated, or castrated and treated with implants containing DHT, E2, or E2 + DHT. Castration increased LH levels. DHT or E2 alone partially suppressed LH, whereas E2 + DHT reduced LH to intact levels. Despite the inhibitory actions on LH, the combination of E2 + DHT increased GnRH neuron activity relative to other treatments, reflected in mean firing rate, amplitude of peaks in firing rate, and area under the curve of firing rate vs. time. Cluster8 was used to identify peaks in firing activity that may be correlated with hormone release. Castration increased the frequency of peaks in firing rate. Treatment with DHT failed to reduce frequency of these peaks. In contrast, treatment with E2 reduced peak frequency to intact levels. The frequency of peaks in firing rate was intermediate in animals treated with E2 + DHT, perhaps suggesting the activating effects of this combination partially counteracts the inhibitory actions of E2. These data indicate that E2 mediates central negative feedback in males primarily by affecting the pattern of GnRH neuron activity, and that androgens combined with estrogens have a central activating effect on GnRH neurons. The negative feedback induced by E2 + DHT to restore LH to intact levels may mask an excitatory central effect of this combination.