On approximations of Stein's neuronal model

Learning shapes cortical dynamics to enhance integration of relevant sensory input

Adaptive sensory behavior is thought to depend on processing in recurrent cortical circuits, but how dynamics in these circuits shapes the integration and transmission of sensory information is not well understood. Here, we study neural coding in recurrently connected networks of neurons driven by sensory input. We show analytically how information available in the network output varies with the alignment between feedforward input and the integrating modes of the circuit dynamics. In light of this theory, we analyzed neural population activity in the visual cortex of mice that learned to discriminate visual features. We found that over learning, slow patterns of network dynamics realigned to better integrate input relevant to the discrimination task. This realignment of network dynamics could be explained by changes in excitatory-inhibitory connectivity among neurons tuned to relevant features. These results suggest that learning tunes the temporal dynamics of cortical circuits to optimally integrate relevant sensory input.

Mean-reverting neuronal model based on two alternating patterns

A neuronal action potential model based on the generalised two-state Ornstein-Uhlenbeck process is studied. The model well describes all phases of a neuronal spike cycle, and the intrinsic parameters of the model have clear specification. Laplace transforms of a firing time are obtained explicitly. Formulae for the mean interspike intervals and their variances, as well as for the average duration of the relative refractory period are also obtained.

Multiplicative noise is beneficial for the transmission of sensory signals in simple neuron models

We study simple integrate-and-fire type models with multiplicative noise and consider the transmission of a weak and slow signal, i.e. a signal that evokes a small modulation of the instantaneous firing rate on time scales that are much larger than the membrane time scale and the mean interspike interval. The specific question of interest is whether and how the state-dependence of the noise can be optimized with respect to information transmission. First, in a simple model in which the noise intensity varies linearly with the state variable, we show analytically that multiplicative fluctuations may benefit the signal transfer and we elucidate the mechanism for this improvement. In a conductance-based integrate-and-fire model with synaptically filtered shot-noise input, we show by means of extended numerical simulations that also in a biophysically more relevant situation, multiplicative noise can enhance the signal-to-noise ratio. Our results shed light on a so far unexplored aspect of stochastic signal transmission in neural systems.

Voltage dependence of synaptic plasticity is essential for rate based learning with short stimuli

In computational neuroscience, synaptic plasticity rules are often formulated in terms of firing rates. The predominant description of in vivo neuronal activity, however, is the instantaneous rate (or spiking probability). In this article we resolve this discrepancy by showing that fluctuations of the membrane potential carry enough information to permit a precise estimate of the instantaneous rate in balanced networks. As a consequence, we find that rate based plasticity rules are not restricted to neuronal activity that is stable for hundreds of milliseconds to seconds, but can be carried over to situations in which it changes every few milliseconds. We illustrate this, by showing that a voltage-dependent realization of the classical BCM rule achieves input selectivity, even if stimulus duration is reduced to a few milliseconds each.

Dynamic Control of Synchronous Activity in Networks of Spiking Neurons

Oscillatory brain activity is believed to play a central role in neural coding. Accumulating evidence shows that features of these oscillations are highly dynamic: power, frequency and phase fluctuate alongside changes in behavior and task demands. The role and mechanism supporting this variability is however poorly understood. We here analyze a network of recurrently connected spiking neurons with time delay displaying stable synchronous dynamics. Using mean-field and stability analyses, we investigate the influence of dynamic inputs on the frequency of firing rate oscillations. We show that afferent noise, mimicking inputs to the neurons, causes smoothing of the system’s response function, displacing equilibria and altering the stability of oscillatory states. Our analysis further shows that these noise-induced changes cause a shift of the peak frequency of synchronous oscillations that scales with input intensity, leading the network towards critical states. We lastly discuss the extension of these principles to periodic stimulation, in which externally applied driving signals can trigger analogous phenomena. Our results reveal one possible mechanism involved in shaping oscillatory activity in the brain and associated control principles.

Stein's neuronal model with pooled renewal input

The input of Stein's model of a single neuron is usually described by using a Poisson process, which is assumed to represent the behaviour of spikes pooled from a large number of presynaptic spike trains. However, such a description of the input is not always appropriate as the variability cannot be separated from the intensity. Therefore, we create and study Stein's model with a more general input, a sum of equilibrium renewal processes. The mean and variance of the membrane potential are derived for this model. Using these formulas and numerical simulations, the model is analyzed to study the influence of the input variability on the properties of the membrane potential and the output spike trains. The generalized Stein's model is compared with the original Stein's model with Poissonian input using the relative difference of variances of membrane potential at steady state and the integral square error of output interspike intervals. Both of the criteria show large differences between the models for input with high variability.

Cooperative behavior in a jump diffusion model for a simple network of spiking neurons

The distribution of time intervals between successive spikes generated by a neuronal cell --the interspike intervals (ISI)-- may reveal interesting features of the underlying dynamics. In this study we analyze the ISI sequence --the spike train-- generated by a simple network of neurons whose output activity is modeled by a jump-diffusion process. We prove that, when specific ranges of the involved parameters are chosen, it is possible to observe multimodal ISI distributions which reveal that the modeled network fires with more than one single preferred time interval. Furthermore, the system exhibits resonance behavior, with modulation of the spike timings by the noise intensity. We also show that inhibition helps the signal transmission between the units of the simple network.

Diffusion approximation of neuronal models revisited

Leaky integrate-and-fire neuronal models with reversal potentials have a number of different diffusion approximations, each depending on the form of the amplitudes of the postsynaptic potentials. Probability distributions of the first-passage times of the membrane potential in the original model and its diffusion approximations are numerically compared in order to find which of the approximations is the most suitable one. The properties of the random amplitudes of postsynaptic potentials are discussed. It is shown on a simple example that the quality of the approximation depends directly on them.

A complex-valued firing-rate model that approximates the dynamics of spiking networks

Firing-rate models provide an attractive approach for studying large neural networks because they can be simulated rapidly and are amenable to mathematical analysis. Traditional firing-rate models assume a simple form in which the dynamics are governed by a single time constant. These models fail to replicate certain dynamic features of populations of spiking neurons, especially those involving synchronization. We present a complex-valued firing-rate model derived from an eigenfunction expansion of the Fokker-Planck equation and apply it to the linear, quadratic and exponential integrate-and-fire models. Despite being almost as simple as a traditional firing-rate description, this model can reproduce firing-rate dynamics due to partial synchronization of the action potentials in a spiking model, and it successfully predicts the transition to spike synchronization in networks of coupled excitatory and inhibitory neurons.

Estimation of Time-Dependent Input from Neuronal Membrane Potential

The set of firing rates of the presynaptic excitatory and inhibitory neurons constitutes the input signal to the postsynaptic neuron. Estimation of the time-varying input rates from intracellularly recorded membrane potential is investigated here. For that purpose, the membrane potential dynamics must be specified. We consider the Ornstein-Uhlenbeck stochastic process, one of the most common single-neuron models, with time-dependent mean and variance. Assuming the slow variation of these two moments, it is possible to formulate the estimation problem by using a state-space model. We develop an algorithm that estimates the paths of the mean and variance of the input current by using the empirical Bayes approach. Then the input firing rates are directly available from the moments. The proposed method is applied to three simulated data examples: constant signal, sinusoidally modulated signal, and constant signal with a jump. For the constant signal, the estimation performance of the method is comparable to that of the traditionally applied maximum likelihood method. Further, the proposed method accurately estimates both continuous and discontinuous time-variable signals. In the case of the signal with a jump, which does not satisfy the assumption of slow variability, the robustness of the method is verified. It can be concluded that the method provides reliable estimates of the total input firing rates, which are not experimentally measurable.

Synchronization in a noise-driven developing neural network

We use computer simulations to investigate the structural and dynamical properties of a developing neural network whose activity is driven by noise. Structurally, the constructed neural networks in our simulations exhibit the small-world properties that have been observed in several neural networks. The dynamical change of neuronal membrane potential is described by the Hodgkin-Huxley model, and two types of learning rules, including spike-timing-dependent plasticity (STDP) and inverse STDP, are considered to restructure the synaptic strength between neurons. Clustered synchronized firing (SF) of the network is observed when the network connectivity (number of connections/maximal connections) is about 0.75, in which the firing rate of neurons is only half of the network frequency. At the connectivity of 0.86, all neurons fire synchronously at the network frequency. The network SF frequency increases logarithmically with the culturing time of a growing network and decreases exponentially with the delay time in signal transmission. These conclusions are consistent with experimental observations. The phase diagrams of SF in a developing network are investigated for both learning rules.

Detecting dependencies between spike trains of pairs of neurons through copulas

The dynamics of a neuron are influenced by the connections with the network where it lies. Recorded spike trains exhibit patterns due to the interactions between neurons. However, the structure of the network is not known. A challenging task is to investigate it from the analysis of simultaneously recorded spike trains. We develop a non-parametric method based on copulas, that we apply to simulated data according to different bivariate Leaky Integrate and Fire models. The method discerns dependencies determined by the surrounding network, from those determined by direct interactions between the two neurons. Furthermore, the method recognizes the presence of delays in the spike propagation. This article is part of a Special Issue entitled "Neural Coding".

Neuronal model with distributed delay: analysis and simulation study for gamma distribution memory kernel

A single neuronal model incorporating distributed delay (memory)is proposed. The stochastic model has been formulated as a Stochastic Integro-Differential Equation (SIDE) which results in the underlying process being non-Markovian. A detailed analysis of the model when the distributed delay kernel has exponential form (weak delay) has been carried out. The selection of exponential kernel has enabled the transformation of the non-Markovian model to a Markovian model in an extended state space. For the study of First Passage Time (FPT) with exponential delay kernel, the model has been transformed to a system of coupled Stochastic Differential Equations (SDEs) in two-dimensional state space. Simulation studies of the SDEs provide insight into the effect of weak delay kernel on the Inter-Spike Interval(ISI) distribution. A measure based on Jensen-Shannon divergence is proposed which can be used to make a choice between two competing models viz. distributed delay model vis-á-vis LIF model. An interesting feature of the model is that the behavior of (CV(t))((ISI)) (Coefficient of Variation) of the ISI distribution with respect to memory kernel time constant parameter η reveals that neuron can switch from a bursting state to non-bursting state as the noise intensity parameter changes. The membrane potential exhibits decaying auto-correlation structure with or without damped oscillatory behavior depending on the choice of parameters. This behavior is in agreement with empirically observed pattern of spike count in a fixed time window. The power spectral density derived from the auto-correlation function is found to exhibit single and double peaks. The model is also examined for the case of strong delay with memory kernel having the form of Gamma distribution. In contrast to fast decay of damped oscillations of the ISI distribution for the model with weak delay kernel, the decay of damped oscillations is found to be slower for the model with strong delay kernel.

Finite post synaptic potentials cause a fast neuronal response

A generic property of the communication between neurons is the exchange of pulses at discrete time points, the action potentials. However, the prevalent theory of spiking neuronal networks of integrate-and-fire model neurons relies on two assumptions: the superposition of many afferent synaptic impulses is approximated by Gaussian white noise, equivalent to a vanishing magnitude of the synaptic impulses, and the transfer of time varying signals by neurons is assessable by linearization. Going beyond both approximations, we find that in the presence of synaptic impulses the response to transient inputs differs qualitatively from previous predictions. It is instantaneous rather than exhibiting low-pass characteristics, depends non-linearly on the amplitude of the impulse, is asymmetric for excitation and inhibition and is promoted by a characteristic level of synaptic background noise. These findings resolve contradictions between the earlier theory and experimental observations. Here we review the recent theoretical progress that enabled these insights. We explain why the membrane potential near threshold is sensitive to properties of the afferent noise and show how this shapes the neural response. A further extension of the theory to time evolution in discrete steps quantifies simulation artifacts and yields improved methods to cross check results.

Errors in estimation of the input signal for integrate-and-fire neuronal models

Estimation of the input parameters of stochastic (leaky) integrate-and-fire neuronal models is studied. It is shown that the presence of a firing threshold brings a systematic error to the estimation procedure. Analytical formulas for the bias are given for two models, the randomized random walk and the perfect integrator. For the third model considered, the leaky integrate-and-fire model, the study is performed by using Monte Carlo simulated trajectories. The bias is compared with other errors appearing during the estimation, and it is documented that the effect of the bias has to be taken into account in experimental studies.

Interspike interval statistics of a leaky integrate-and-fire neuron driven by Gaussian noise with large correlation times

We analytically investigate the interspike interval (ISI) density, the Fano factor, and the coefficient of variation of a leaky integrate-and-fire neuron model driven by exponentially correlated Gaussian noise with a large correlation time tau . We find a burstinglike behavior of the spike train, which is revealed by a dominant peak of the ISI density at small intraburst intervals and a slow power-law decay of long interburst intervals. The large, power-law distributed ISIs give rise to a coefficient of variation which diverges as square root [tau] . This leads to the paradoxical effect that ISI correlations, as expressed by the serial correlation coefficient, vanish for large correlation times. This is in contrast to findings of previous works on a simpler neuron model where the effect of noise correlations appeared in higher-order statistical measures.

Analysis of a stochastic neuronal model with excitatory inputs and state-dependent effects

We propose a stochastic model for the firing activity of a neuronal unit. It includes the decay effect of the membrane potential in absence of stimuli, and the occurrence of time-varying excitatory inputs governed by a Poisson process. The sample-paths of the membrane potential are piecewise exponentially decaying curves with jumps of random amplitudes occurring at the input times. An analysis of the probability distributions of the membrane potential and of the firing time is performed. In the special case of time-homogeneous stimuli the firing density is obtained in closed form, together with its mean and variance.

A review of the integrate-and-fire neuron model: I. Homogeneous synaptic input

The integrate-and-fire neuron model is one of the most widely used models for analyzing the behavior of neural systems. It describes the membrane potential of a neuron in terms of the synaptic inputs and the injected current that it receives. An action potential (spike) is generated when the membrane potential reaches a threshold, but the actual changes associated with the membrane voltage and conductances driving the action potential do not form part of the model. The synaptic inputs to the neuron are considered to be stochastic and are described as a temporally homogeneous Poisson process. Methods and results for both current synapses and conductance synapses are examined in the diffusion approximation, where the individual contributions to the postsynaptic potential are small. The focus of this review is upon the mathematical techniques that give the time distribution of output spikes, namely stochastic differential equations and the Fokker-Planck equation. The integrate-and-fire neuron model has become established as a canonical model for the description of spiking neurons because it is capable of being analyzed mathematically while at the same time being sufficiently complex to capture many of the essential features of neural processing. A number of variations of the model are discussed, together with the relationship with the Hodgkin-Huxley neuron model and the comparison with electrophysiological data. A brief overview is given of two issues in neural information processing that the integrate-and-fire neuron model has contributed to - the irregular nature of spiking in cortical neurons and neural gain modulation.

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.

Estimation of the input parameters in the Ornstein-Uhlenbeck neuronal model

The stochastic Ornstein-Uhlenbeck neuronal model is studied, and estimators of the model input parameters, depending on the firing regime of the process, are derived. Closed expressions for the Laplace transforms of the first two moments of the normalized first-passage time through a constant boundary in the suprathreshold regime are derived, which is used to define moment estimators. In the subthreshold regime, the exponentiality of the first-passage time is utilized to characterize the input parameters. In the threshold regime and for the Wiener process approximation, analytic expressions for the first-passage-time density are used to derive the maximum-likelihood estimators of the parameters. The methods are illustrated on simulated data under different conditions, including misspecification of the intrinsic parameters of the model. Finally, known approximations of the first-passage-time moments are improved.

Recurrent biological neural networks: the weak and noisy limit

A perturbative method is developed for calculating the effects of recurrent synaptic interactions between neurons embedded in a network. A series expansion is constructed that converges for networks with noisy membrane potential and weak synaptic connectivity. The terms of the series can be interpreted as loops of interactions between neurons, so the technique is called a loop expansion. A diagrammatic method is introduced that allows for construction of analytic expressions for the parameter dependencies of the spike-probability function and correlation functions. An analytic expression is obtained to predict the effect of the surrounding network on a neuron during an intracellular current injection. The analytic results are compared with simulations to test the range of their validity and significant effects of the the recurrent connections in network are accurately predicted by the loop expansion.

Duality of rate coding and temporal coding in multilayered feedforward networks

A functional role for precise spike timing has been proposed as an alternative hypothesis to rate coding. We show in this article that both the synchronous firing code and the population rate code can be used dually in a common framework of a single neural network model. Furthermore, these two coding mechanisms are bridged continuously by several modulatable model parameters, including shared connectivity, feedback strength, membrane leak rate, and neuron heterogeneity. The rates of change of these parameters are closely related to the response time and the timescale of learning.

Interspike interval statistics in the Ornstein-Uhlenbeck neuronal model with signal-dependent noise

The stochastic leaky integrate-and-fire (LIF) continuous model is studied under the condition that the amplitude of noise is a function of the input signal. The coefficient of variation (CV) of interspike intervals (ISIs) is investigated for different types of dependencies between the noise and the signal. Finally, we present the CV and the ISI density resulting from the special choice of parameters of the input that gave rise to a contra-intuitive behavior of the transfer function in Lánský and Sacerdote [Phys. Lett. A 285 (2001) 132].

Properties of the extra-positional signal in hippocampal place cell discharge derived from the overdispersion in location-specific firing

There is a good deal of evidence that in the rodent, an internal model of the external world is encoded by hippocampal pyramidal cells, called 'place cells'. During free exploration, the activity of place cells is higher within a small part of the space, called the firing field, and virtually silent elsewhere. We have previously shown that the spiking activity during passes through the firing field is characterized not only by the high firing rate, but also by its very high variability ('overdispersion'). This overdispersion indicates that place cells carry information in addition to position. Here we demonstrate by simulations of an integrate-and-fire neuronal model that while a rat is foraging in an open space this additional information may arise from a process that alternatingly modulates the inputs to place cells by about 10% with a mean period of about 1 s. We propose that the overdispersion reflects switches of the rats attention between different spatial reference frames of the environment. This predicts that the overdispersion will not be observed in rats that use only room-based cues for navigation. We show that while place cell firing is overdispersed in rats during foraging in an open arena, the firing is less overdispersed during the same behavior in the same environment, when the rats have been trained to use only room-based and not arena-based cues to navigate.

Shot noise in the leaky integrate-and-fire neuron

We study the influence of noise on the transmission of temporal information by a leaky integrate-and-fire neuron using the theory of shot noise. The model includes a finite number of synapses and has a membrane potential variance de facto modulated by the input signal. The phenomenon of stochastic resonance in spiking neurons is analytically exhibited using an inhomogeneous Poisson process model of the spike trains, and links with the traditional Ornstein-Uhlenbeck process obtained by a diffusion approximation are given. It is shown that the modulated membrane potential variance inherent to the model gives better signal processing capabilities than the diffusion approximation.

Calculation of interspike intervals for integrate-and-fire neurons with poisson distribution of synaptic inputs

We present a new technique for calculating the interspike intervals of integrate-and-fire neurons. There are two new components to this technique. First, the probability density of the summed potential is calculated by integrating over the distribution of arrival times of the afferent post-synaptic potentials (PSPs), rather than using conventional stochastic differential equation techniques. A general formulation of this technique is given in terms of the probability distribution of the inputs and the time course of the postsynaptic response. The expressions are evaluated in the gaussian approximation, which gives results that become more accurate for large numbers of small-amplitude PSPs. Second, the probability density of output spikes, which are generated when the potential reaches threshold, is given in terms of an integral involving a conditional probability density. This expression is a generalization of the renewal equation, but it holds for both leaky neurons and situations in which there is no time-translational invariance. The conditional probability density of the potential is calculated using the same technique of integrating over the distribution of arrival times of the afferent PSPs. For inputs with a Poisson distribution, the known analytic solutions for both the perfect integrator model and the Stein model (which incorporates membrane potential leakage) in the diffusion limit are obtained. The interspike interval distribution may also be calculated numerically for models that incorporate both membrane potential leakage and a finite rise time of the postsynaptic response. Plots of the relationship between input and output firing rates, as well as the coefficient of variation, are given, and inputs with varying rates and amplitudes, including inhibitory inputs, are analyzed. The results indicate that neurons functioning near their critical threshold, where the inputs are just sufficient to cause firing, display a large variability in their spike timings.

Input parameters in a one-dimensional neuronal model with reversal potentials

An equation for a stochastic neuronal model describing the initiation of action potentials is studied for the case in which the depolarization of the membrane potential is restricted by the reversal potentials. It is assumed that the values of the membrane potential can be continuously recorded between consecutive spikes. Under this assumption, the estimators of the model parameters are derived and the methods for testing the model are proposed. The objective of the methods presented in this contribution is to provide neuroscientists with quantitative means in order to estimate parameters of stochastic neuronal models.

On the relationship between synaptic input and spike output jitter in individual neurons

What is the relationship between the temporal jitter in the arrival times of individual synaptic inputs to a neuron and the resultant jitter in its output spike? We report that the rise time of firing rates of cells in striate and extrastriate visual cortex in the macaque monkey remain equally sharp at different stages of processing. Furthermore, as observed by others, multiunit recordings from single units in the primate frontal lobe reveal a strong peak in their cross-correlation in the 10-150 msec range with very small temporal jitter (on the order of 1 msec). We explain these results using numerical models to study the relationship between the temporal jitter in excitatory and inhibitory synaptic input and the variability in the spike output timing in integrate-and-fire units and in a biophysically and anatomically detailed model of a cortical pyramidal cell. We conclude that under physiological circumstances, the standard deviation in the output jitter is linearly related to the standard deviation in the input jitter, with a constant of less than one. Thus, the timing jitter in successive layers of such neurons will converge to a small value dictated by the jitter in axonal propagation times.

On the comparison of Feller and Ornstein-Uhlenbeck models for neural activity

Diffusion processes have been extensively used to describe membrane potential behavior. In this approach the interspike interval has a theoretical counterpart in the first-passage-time of the diffusion model employed. Since the mathematical complexity of the first-passage-time problem increases with attempts to make the models more realistic it seems useful to compare the features of different models in order to highlight their relative performance. In this paper we compare the Feller and Ornstein-Uhlenbeck models under three different criteria derived from the level of information available about their parameters. We conclude that the Feller model is preferable when complete knowledge of the characterizing parameters is assumed. On the other hand, when only limited information about the parameters is available, such as the mean firing time and the histogram shape, no advantage arises from using this more complex model.

Stochastic model neuron without resetting of dendritic potential: application to the olfactory system

A two-dimensional neuronal model, in which the membrane potential of the dendrite evolves independently from that at the trigger zone of the axon, is proposed and studied. In classical one-dimensional neuronal models the dendritic and axonal potentials cannot be distinguished, and thus they are reset to resting level after firing of an action potential, whereas in the present model the dendritic potential is not reset. The trigger zone is modelled by a simplified leaky integrator (RC circuit) and the dendritic compartment can be described by any of the classical one-dimensional neuronal models. The new model simulates observed features of the firing dynamics which are not displayed by classical models, namely positive correlation between interspike intervals and endogenous bursting. It gives a more natural account of features already accounted for in previous models, such as the absence of an upper limit for the coefficient of variation of intervals (i.e. irregular firing). It allows the first- and second-order neurons of the olfactory system to be described with the same basic assumptions, which was not the case in one-point models. Nevertheless it keeps the main qualitative properties found previously, such as the existence of three regimens of firing with increasing stimulus concentration and the sigmoid shape of the firing frequency of first-order neurons as a function of the logarithm of stimulus concentration.

A diffusion model of a neuron and neural nets

In the paper a diffusion model of a neuron is treated. A new, less restrictive than usually, condition of applicability of a diffusion model is presented. As a result the point-process-to-point-process model of a neuron is obtained, which produces an output signal of the same kind as the accepted input signals.

Coding of odor intensity

A model for coding of odor intensity in the first two neuronal layers of olfactory systems is proposed. First, the occupation and activation by odorant molecules of receptor proteins of different types borne by the first order neurons are described as birth and death processes. The occupation (birth) rate depends on the concentration of the odorant, whereas the probability of activation of an occupied receptor depends on the type of the odorant. Second, the spike generation mechanism proposed for the first order neuron depends on the level of the generator potential evoked by the activated receptors and on a time-decaying threshold which is reset to infinity after each spike. The various resulting stochastic regimes of firing activity at different concentrations are described. Third, each second order neuron is influenced by excitation coming from numerous first order neurons, lateral inhibition from other second order neurons, and self-inhibition. All these incoming signals are integrated at the second order neuron. The firing activity of the first and second order neurons is modeled by a first passage time scheme. For both types of neuron the shapes of the curves predicted by the model for the mean firing frequency as a function of stimulus concentration are shown to be in accordance with available experimental results.

Simulation of a diffusion process with randomly distributed jumps in neuronal context

In stochastic neuronal models, an interspike interval corresponds to the time interval during which the process imitating the membrane potential reaches a threshold from an initial depolarization. For neurons with an extensive dendritic structure, a stochastic process combining diffusion and discontinuous development of its trajectory is considered a good description of the membrane potential. Due to a lack of analytical solutions of the threshold passage distribution for such a process, a method for computer simulation is introduced here. For the diffusion Ornstein-Uhlenbeck process with exponentially distributed moments of constant jumps a program is given. The relation between the simulation step, accuracy of simulation and amount of computing time required is discussed.

Effects of afterhyperpolarization on neuronal firing

Stein's model for a neuron is studied. This model is modified to take into account the effects of afterhyperpolarization on the neuronal firing. The relative refractory phase, following the absolute one, is modelled by a time-increasing amplitude of postsynaptic potentials and it is also incorporated into the model. Besides the simulation of the model, some theoretical results and approximation methods are derived. Afterhyperpolarization tends to preserve the linearity of the frequency transfer characteristic and it has a limited effect on the moments of the interspike intervals in general. The main effects are seen at high firing rates and in the removal of short intervals in the interspike interval histogram.

The thalamocortical contribution to epilepsy

The experimental literature has dealt intensively with the cortical contribution to epilepsy. Possibly because of the direction of technological advance, much less attention has been paid to the role of other structures. A model which emphasizes the role of some of those non-cortical structures, specifically that of thalamocortical modulation of cortical excitability, is developed. Some aspects of the petit mal seizure, a seizure type considered by some investigators to involve thalamocortical mechanisms, are predicted by the model. Although the thalamocortical mechanisms under study are not the only mechanisms underlying seizures, a full understanding of the phenomenology of epilepsy needs to take into account the role of subcortical modification of cortical activities in addition to other mechanisms. Gloor has described two types of epileptogenesis: type I characteristic of non-convulsive seizures and type II characteristics of convulsions. There is disagreement as to whether or not the two mechanisms represent qualitatively different phenomena. Utilizing the thalamocortical model, it can be shown that the two types of epileptogenesis are qualitatively different. Furthermore, the thalamocortical model leads to a possible explanation of clinically different profiles of antiepileptic efficacy of medications.

Generalized Stein's model for anatomically complex neurons

A neuron with a large dendritic structure is considered. The number of synapses located on the dendrites is substantially higher than on the soma. The synaptic input effect on the neuronal excitability decreases with distance between a synapse ending and the trigger zone. Two areas are distinguished in accordance with the effect of synaptic input--dendritic and somatic. The dendritic area, when compared to the soma, is characterized by much higher intensity of its activation but the amplitudes of synaptically evoked changes of the membrane potential at the trigger zone are in general small. This situation is suitable for a diffusion approximation. However, on the soma, especially in the proximity of the trigger zone, the membrane potential changes are a large fraction of the threshold depolarization. The membrane potential at the trigger zone is modelled by a one-dimensional stochastic process. The diffusion Ornstein-Uhlenbeck process serves as a basis of the model; however, at the moments of somatic synapses activation its voltage changes in jumps. Their sizes represent the amplitudes of the evoked postsynaptic potentials. The unimodal histograms of interspike intervals can be explained by the model. The values of the coefficient of variation greater than one are connected with substantial inhibition.

A stochastic model for the membrane potential of a stimulated neuron

We present a simple model describing the transition between the prefiring, firing and postfiring phases of a single neuron in a large neural net. Using typical values for the physiological parameters that enter the model, we find average interspike times that are close to those reported in experimental measurements.

The effect of a random initial value in neural first-passage-time models

The effect of a random initial value is examined in several stochastic integrate-and-fire neural models with a constant threshold and a constant input. The three models considered are approximations of Stein's model, namely: (1) a leaky integrator with deterministic trajectories, (2) a Wiener process with drift, and (3) an Ornstein-Uhlenbeck process. For model 1, different distributions for the initial value lead to commonly observed interspike interval distributions. For model 2, a discrete and a uniform distribution for the initial value are examined along with some parameter estimation procedures. For model 3, with a truncated normal distribution for the initial value, the coefficient of variation is shown to be greater than 1, and as the threshold becomes large the first-passage-time distribution approaches an exponential distribution. The relationships among the models and between them and previous models are also discussed, along with the robustness of the model assumptions and methods of their verification. The effects of a random initial value are found to be most pronounced at high firing rates.