Stochastic modeling of the neuronal activity in the subthalamic nucleus and model parameter identification from Parkinson patient data

Several stochastic models, with various degrees of complexity, have been proposed to model the neuronal activity from different parts of the human brain. In this article, we use a simple Ornstein–Uhlenbeck process (OUP) to model the spike activity recorded from the subthalamic nucleus of patients suffering from Parkinson’s disease at the time of implantation of the electrodes for deep brain stimulation. From the recorded data, which contains information about the spike times of a single neuron, we identify and extract the model parameters of the OUP. We then use these parameters to numerically simulate the inter-spike intervals and the voltage across the neuron membrane. We finally assess how well the proposed mathematical model fits to the measured data and compare it with other commonly adopted stochastic models. We show an excellent agreement between the computer-generated data according to the OUP model and the measured one, as well as the superiority of the OUP model when compared to the Poisson process model and the random walk model; thus, establishing the validity of the OUP as a simple yet biologically plausible model of the neuronal activity recorded from the subthalamic nucleus of Parkinson’s disease patients.     

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. Received: 8 November 1994/Accepted in revised form : 23 May 1995     

The parameters of the stochastic leaky integrate-and-fire neuronal model

Five parameters of one of the most common neuronal models, the diffusion leaky integrate-and-fire model, also known as the Ornstein-Uhlenbeck neuronal model, were estimated on the basis of intracellular recording. These parameters can be classified into two categories. Three of them (the membrane time constant, the resting potential and the firing threshold) characterize the neuron itself. The remaining two characterize the neuronal input. The intracellular data were collected during spontaneous firing, which in this case is characterized by a Poisson process of interspike intervals. Two methods for the estimation were applied, the regression method and the maximum-likelihood method. Both methods permit to estimate the input parameters and the membrane time constant in a short time window (a single interspike interval). We found that, at least in our example, the regression method gave more consistent results than the maximum-likelihood method. The estimates of the input parameters show the asymptotical normality, which can be further used for statistical testing, under the condition that the data are collected in different experimental situations. The model neuron, as deduced from the determined parameters, works in a subthreshold regimen. This result was confirmed by both applied methods. The subthreshold regimen for this model is characterized by the Poissonian firing. This is in a complete agreement with the observed interspike interval data. Action Editor: Nicolas Brunel     

A spatial stochastic neuronal model with Ornstein–Uhlenbeck input current

 We consider a spatial neuron model in which the membrane potential satisfies a linear cable equation with an input current which is a dynamical random process of the Ornstein–Uhlenbeck (OU) type. This form of current may represent an approximation to that resulting from the random opening and closing of ion channels on a neuron's surface or to randomly occurring synaptic input currents with exponential decay. We compare the results for the case of an OU input with those for a purely white-noise-driven cable model. The statistical properties, including mean, variance and covariance of the voltage response to an OU process input in the absence of a threshold are determined analytically. The mean and the variance are calculated as a function of time for various synaptic input locations and for values of the ratio of the time constant of decay of the input current to the time constant of decay of the membrane voltage in the physiological range for real neurons. The limiting case of a white-noise input current is obtained as the correlation time of the OU process approaches zero. The results obtained with an OU input current can be substantially different from those in the white-noise case. Using simulation of the terms in the series representation for the solution, we estimate the interspike interval distribution for various parameter values, and determine the effects of the introduction of correlation in the synaptic input stochastic process. Received: 5 March 2001 / Accepted in revised form: 7 August 2001     

Effects of active versus passive dendritic membranes on the transfer properties of a simulated neuron

In a simulated neuron with a dendritic tree, the relative effects of active and passive dendritic membranes on transfer properties were studied. The simulations were performed by means of a digital computer. The computations calculated the changes in transmembrane voltages of many compartments over time as a function of other biophysical variables. These variables were synaptic input intensity, critical firing threshold, rate of leakage of current across the membrane, and rate of longitudinal current spread between compartments. For both passive and active dendrites, the transfer properties of the soma studied for different rates of longitudinal current spread. With low rates of current spread, graded changes in firing threshold produced correspondingly graded changes in output discharge. With high rates of current spread, the neuron became a bistable operator where spiking was enhanced if the threshold was below a certain level and suppressed if the threshold was above that level. Since alterations in firing threshold were shown to have the same effect on firing rate as alterations in synaptic input intensity, the neuron can be said to change from graded to contrast-enhancing in its response to stimuli of different intensities. The presence or absence of dendritic spiking was found to have a significant effect on the integrative properties of the simulated neuron. In particular, contrast enhancement was considerably more pronounced in neurons with passive than with active dendrites in that somatic spike rates reached a higher maximum when dendrites were passive. With active dendrites, a less intense input was needed to initiate somatic spiking than with passive dendrites because a distal dendritic spike could easily propagate by means of longitudinal current spread to the soma. Once somatic spiking was initiated, though, spike rates tended to be lower with active than with passive dendrites because the soma recovered more slowly from its post-spike refractory period if it was also influenced by refractory periods in the dendrites. The experiment of comparing neurons with active and passive dendrites was repeated at a different, higher value of synaptic input. The same differences in transfer properties between the active and passive cases emerged as before. Spiking patterns in neurons with active dendrites were also affected by the time distribution of synaptic inputs. In a previous study, inputs had been random over both space and time, varying about a predetermined mean, whereas in the present study, inputs were random over space but uniform over time. When inputs were made uniform over time, spiking became more difficult to initiate and the transition from graded to bistable response became less sharp.     

A stochastic SIS epidemic with demography: initial stages and time to extinction

We study an open population stochastic epidemic model from the time of introduction of the disease, through a possible outbreak and to extinction. The model describes an SIS (susceptible–infective–susceptible) epidemic where all individuals, including infectious ones, reproduce at a given rate. An approximate expression for the outbreak probability is derived using a coupling argument. Further, we analyse the behaviour of the model close to quasi-stationarity, and the time to disease extinction, with the aid of a diffusion approximation. In this situation the number of susceptibles and infectives behaves as an Ornstein–Uhlenbeck process, centred around the stationary point, for an exponentially distributed time before going extinct.     

Influence of temporal correlation of synaptic input on the rate and variability of firing in neurons

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

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     

Relationship between applicability of current-based synapses and uniformity of firing patterns

The purpose of this paper is to identify situations in neural network modeling where current-based synapses are applicable. The applicability of current-based synapse model for studying post-transient behavior of neural networks is discussed in terms of average synaptic current strength induced by per spike during one firing cycle of a neuron (or briefly per spike synaptic current strength). It was found that current-based synapse models are applicable in both situations where both the interspike intervals of the neurons and the distribution of firing times of the neurons are uniform, and where the firing of all neurons is synchronized. If neither the interspike intervals nor the distribution of firing times of the neurons is uniform or the reversal potential is between the rest and threshold potentials, current-based synapse models may be oversimplified.     

A cross-interval spike train analysis: the correlation between spike generation and temporal integration of doublets

A stochastic spike train analysis technique is introduced to reveal the correlation between the firing of the next spike and the temporal integration period of two consecutive spikes (i.e., a doublet). Statistics of spike firing times between neurons are established to obtain the conditional probability of spike firing in relation to the integration period. The existence of a temporal integration period is deduced from the time interval between two consecutive spikes fired in a reference neuron as a precondition to the generation of the next spike in a compared neuron. This analysis can show whether the coupled spike firing in the compared neuron is correlated with the last or the second-to-last spike in the reference neuron. Analysis of simulated and experimentally recorded biological spike trains shows that the effects of excitatory and inhibitory temporal integration are extracted by this method without relying on any subthreshold potential recordings. The analysis also shows that, with temporal integration, a neuron driven by random firing patterns can produce fairly regular firing patterns under appropriate conditions. This regularity in firing can be enhanced by temporal integration of spikes in a chain of polysynaptically connected neurons. The bandpass filtering of spike firings by temporal integration is discussed. The results also reveal that signal transmission delays may be attributed not just to conduction and synaptic delays, but also to the delay time needed for temporal integration. Received: 3 March 1997 / Accepted in revised form: 6 November 1997     

On some computational results for single neurons' activity modeling

The classical Ornstein-Uhlenbeck diffusion neuronal model is generalized by inclusion of a time-dependent input whose strength exponentially decreases in time. The behavior of the membrane potential is consequently seen to be modeled by a process whose mean and covariance classify, it as Gaussian-Markov. The effect of the input on the neuron's firing characteristics is investigated by comparing the firing probability densities and distributions for such a process with the corresponding ones of the Ornstein-Uhlenbeck model. All numerical results are obtained by implementation of a recently developed computational method.     

Leader neurons in leaky integrate and fire neural network simulations

In this paper, we highlight the topological properties of leader neurons whose existence is an experimental fact. Several experimental studies show the existence of leader neurons in population bursts of activity in 2D living neural networks (Eytan and Marom, J Neurosci 26(33):8465–8476, 2006; Eckmann et al., New J Phys 10(015011), 2008). A leader neuron is defined as a neuron which fires at the beginning of a burst (respectively network spike) more often than we expect by chance considering its mean firing rate. This means that leader neurons have some burst triggering power beyond a chance-level statistical effect. In this study, we characterize these leader neuron properties. This naturally leads us to simulate neural 2D networks. To build our simulations, we choose the leaky integrate and fire (lIF) neuron model (Gerstner and Kistler 2002; Cessac, J Math Biol 56(3):311–345, 2008), which allows fast simulations (Izhikevich, IEEE Trans Neural Netw 15(5):1063–1070, 2004; Gerstner and Naud, Science 326:379–380, 2009). The dynamics of our lIF model has got stable leader neurons in the burst population that we simulate. These leader neurons are excitatory neurons and have a low membrane potential firing threshold. Except for these two first properties, the conditions required for a neuron to be a leader neuron are difficult to identify and seem to depend on several parameters involved in the simulations themselves. However, a detailed linear analysis shows a trend of the properties required for a neuron to be a leader neuron. Our main finding is: A leader neuron sends signals to many excitatory neurons as well as to few inhibitory neurons and a leader neuron receives only signals from few other excitatory neurons. Our linear analysis exhibits five essential properties of leader neurons each with different relative importance. This means that considering a given neural network with a fixed mean number of connections per neuron, our analysis gives us a way of predicting which neuron is a good leader neuron and which is not. Our prediction formula correctly assesses leadership for at least ninety percent of neurons.     

Sensitive dependence of the coefficient of variation of interspike intervals on the lower boundary of membrane potential for the leaky integrate-and-fire neuron model

After the report of Softky and Koch [Softky, W.R., Koch, C., 1993. The highly irregular firing of cortical cells is inconsistent with temporal integration of random EPSPs. J. Neurosci. 13, 334-350], leaky integrate-and-fire models have been investigated to explain high coefficient of variation (CV) of interspike intervals (ISIs) at high firing rates observed in the cortex. The purpose of this paper is to study the effect of the position of a lower boundary of membrane potential on the possible value of CV of ISIs based on the diffusional leaky integrate-and-fire models with and without reversal potentials. Our result shows that the irregularity of ISIs for the diffusional leaky integrate-and-fire neuron significantly changes by imposing a lower boundary of membrane potential, which suggests the importance of the position of the lower boundary as well as that of the firing threshold when we study the statistical properties of leaky integrate-and-fire neuron models. It is worth pointing out that the mean-CV plot of ISIs for the diffusional leaky integrate-and-fire neuron with reversal potentials shows a close similarity to the experimental result obtained in Softky and Koch [Softky, W.R., Koch, C., 1993. The highly irregular firing of cortical cells is inconsistent with temporal integration of random EPSPs. J. Neurosci. 13, 334-350].     

Dynamical features of higher-order correlation events: impact on cortical cells

Cortical neurons receive signals from thousands of other neurons. The statistical properties of the input spike trains substantially shape the output response properties of each neuron. Experimental and theoretical investigations have mostly focused on the second order statistical features of the input spike trains (mean firing rates and pairwise correlations). Little is known of how higher order correlations affect the integration and firing behavior of a cell independently of the second order statistics. To address this issue, we simulated the dynamics of a population of 5000 neurons, controlling both their second order and higher-order correlation properties to reflect physiological data. We then used these ensemble dynamics as the input stage to morphologically reconstructed cortical cells (layer 5 pyramidal, layer 4 spiny stellate cell), and to an integrate and fire neuron. Our results show that changes done solely to the higher-order correlation properties of the network’s dynamics significantly affect the response properties of a target neuron, both in terms of output rate and spike timing. Moreover, the neuronal morphology and voltage dependent mechanisms of the target neuron considerably modulate the quantitative aspects of these effects. Finally, we show how these results affect sparseness of neuronal representations, tuning properties, and feature selectivity of cortical cells. An erratum to this article can be found at     

A stochastic model and a functional central limit theorem for information processing in large systems of neurons

The paper deals with information transmission in large systems of neurons. We model the membrane potential in a single neuron belonging to a cell tissue by a non time-homogeneous Cox-Ingersoll-Ross type diffusion; in terms of its time-varying expectation, this stochastic process can convey deterministic signals. We model the spike train emitted by this neuron as a Poisson point process compensated by the occupation time of the membrane potential process beyond the excitation threshold. In a large system of neurons 1≤i≤N processing independently the same deterministic signal, we prove a functional central limit theorem for the pooled spike train collected from the N neurons. This pooled spike train allows to recover the deterministic signal, up to some shape transformation which is explicit.     

Short-window spectral analysis using AMVAR and multitaper methods: a comparison

We compare two popular methods for estimating the power spectrum from short data windows, namely the adaptive multivariate autoregressive (AMVAR) method and the multitaper method. By analyzing a simulated signal (embedded in a background Ornstein–Uhlenbeck noise process) we demonstrate that the AMVAR method performs better at detecting short bursts of oscillations compared to the multitaper method. However, both methods are immune to jitter in the temporal location of the signal. We also show that coherence can still be detected in noisy bivariate time series data by the AMVAR method even if the individual power spectra fail to show any peaks. Finally, using data from two monkeys performing a visuomotor pattern discrimination task, we demonstrate that the AMVAR method is better able to determine the termination of the beta oscillations when compared to the multitaper method.     

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.     

Library-based numerical reduction of the Hodgkin–Huxley neuron for network simulation

We present an efficient library-based numerical method for simulating the Hodgkin–Huxley (HH) neuronal networks. The key components in our numerical method involve (i) a pre-computed high resolution data library which contains typical neuronal trajectories (i.e., the time-courses of membrane potential and gating variables) during the interval of an action potential (spike), thus allowing us to avoid resolving the spikes in detail and to use large numerical time steps for evolving the HH neuron equations; (ii) an algorithm of spike-spike corrections within the groups of strongly coupled neurons to account for spike-spike interactions in a single large time step. By using the library method, we can evolve the HH networks using time steps one order of magnitude larger than the typical time steps used for resolving the trajectories without the library, while achieving comparable resolution in statistical quantifications of the network activity, such as average firing rate, interspike interval distribution, power spectra of voltage traces. Moreover, our large time steps using the library method can break the stability requirement of standard methods (such as Runge–Kutta (RK) methods) for the original dynamics. We compare our library-based method with RK methods, and find that our method can capture very well phase-locked, synchronous, and chaotic dynamics of HH neuronal networks. It is important to point out that, in essence, our library-based HH neuron solver can be viewed as a numerical reduction of the HH neuron to an integrate-and-fire (I&F) neuronal representation that does not sacrifice the gating dynamics (as normally done in the analytical reduction to an I&F neuron).     

A review of the integrate-and-fire neuron model: II. Inhomogeneous synaptic input and network properties

The integrate-and-fire neuron model describes the state of a neuron in terms of its membrane potential, which is determined by the synaptic inputs and the injected current that the neuron receives. When the membrane potential reaches a threshold, an action potential (spike) is generated. This review considers the model in which the synaptic input varies periodically and is described by an inhomogeneous Poisson process, with both current and conductance synapses. The focus is on the mathematical methods that allow the output spike distribution to be analyzed, including first passage time methods and the Fokker–Planck equation. Recent interest in the response of neurons to periodic input has in part arisen from the study of stochastic resonance, which is the noise-induced enhancement of the signal-to-noise ratio. Networks of integrate-and-fire neurons behave in a wide variety of ways and have been used to model a variety of neural, physiological, and psychological phenomena. The properties of the integrate-and-fire neuron model with synaptic input described as a temporally homogeneous Poisson process are reviewed in an accompanying paper (Burkitt in Biol Cybern, 2006).     

Cooperative mechanism for improving the discriminating ability in the chemoreceptor neuron Binomial case

The discriminating ability (selectivity) of the chemoreceptor neuron is compared with that of its receptor proteins. The process of neuronal triggering is expected to be cooperative and threshold type in a sense that the neuron fires a spike if and only if the number of receptor proteins which are bound with odor molecules is above a definite threshold. The binomial distribution is utilized to estimate the firing probability if a definite odor is applied. It is established that a chemoreceptor neuron can have a much higher selectivity than its individual receptor proteins, provided that the chemical stimuli are presented at low concentrations. A possibility for the above mechanism to be valid in other sensory systems is discussed. Received: 20 July 1998 / Accepted in revised form: 30 April 1999