Cooperative behavior in the periodically modulated Wiener process: Noise-induced complexity in a model neutron

Low-dimensional firing-rate dynamics for populations of renewal-type spiking neurons

The macroscopic dynamics of large populations of neurons can be mathematically analyzed using low-dimensional firing-rate or neural-mass models. However, these models fail to capture spike synchronization effects and nonstationary responses of the population activity to rapidly changing stimuli. Here we derive low-dimensional firing-rate models for homogeneous populations of neurons modeled as time-dependent renewal processes. The class of renewal neurons includes integrate-and-fire models driven by white noise and has been frequently used to model neuronal refractoriness and spike synchronization dynamics. The derivation is based on an eigenmode expansion of the associated refractory density equation, which generalizes previous spectral methods for Fokker-Planck equations to arbitrary renewal models. We find a simple relation between the eigenvalues characterizing the timescales of the firing rate dynamics and the Laplace transform of the interspike interval density, for which explicit expressions are available for many renewal models. Retaining only the first eigenmode already yields a reliable low-dimensional approximation of the firing-rate dynamics that captures spike synchronization effects and fast transient dynamics at stimulus onset. We explicitly demonstrate the validity of our model for a large homogeneous population of Poisson neurons with absolute refractoriness and other renewal models that admit an explicit analytical calculation of the eigenvalues. The eigenmode expansion presented here provides a systematic framework for alternative firing-rate models in computational neuroscience based on spiking neuron dynamics with refractoriness.

DNA breathing dynamics under periodic forcing: Study of several distribution functions of relevant Brownian functionals

In this paper, we study DNA breathing dynamics in the presence of an external periodic force by proposing and inspecting several probability distribution functions (PDFs) of relevant Brownian functionals which specify the bubble lifetime, reactivity, and average size. We model the bubble dynamics process by an overdamped Langevin equation of broken base pairs on the Poland-Scheraga free energy landscape. Introducing an effective time-independent description for timescales larger than T[over ̃]=2π/ω (where ω is the frequency of external periodic force) and using an elegant backward Fokker-Planck method we derive closed form expressions of several PDFs associated with such stochastic processes. For instance, with an initial bubble size of x_{0}, we derive the following analytical expressions: (i) the PDF P(t_{f}|x_{0}) of the first passage time t_{f} which specifies the lifetime of the DNA breathing process, (ii) the PDF P(A|x_{0}) of the area A until the first passage time, and it provides much valuable information about the average bubble size and reactivity of the process, and (iii) the PDF P(M) associated with the maximum bubble size M of the breathing process before complete denaturation. Our analysis is limited to two limits: (a) large bubble size and (b) small bubble size. We further confirm our analytical predictions by computing the same PDFs with direct numerical simulations of the corresponding Langevin equations. We obtain very good agreement of our theoretical predictions with the numerically simulated results. Finally, several nontrivial scaling behaviors in the asymptotic limits for the above-mentioned PDFs are predicted, which can be verified further from experimental observation. Our main conclusion is that the large bubble dynamics is unaffected by the rapidly oscillating force, but the small bubble dynamics is significantly affected by the same periodic force.

New Type of Spectral Nonlinear Resonance Enhances Identification of Weak Signals

Some nonlinear systems possess innate capabilities of enhancing weak signal transmissions through a unique process called Stochastic Resonance (SR). However, existing SR mechanism suffers limited signal enhancement from inappropriate entraining signals. Here we propose a new and effective implementation, resulting in a new type of spectral resonance similar to SR but capable of achieving orders of magnitude higher signal enhancement than previously reported. By employing entraining frequency in the range of the weak signal, strong spectral resonances can be induced to facilitate nonlinear modulations and intermodulations, thereby strengthening the weak signal. The underlying physical mechanism governing the behavior of spectral resonances is examined, revealing the inherent advantages of the proposed spectral resonances over the existing implementation of SR. Wide range of parameters have been found for the optimal enhancement of any given weak signal and an analytical method is established to estimate these required parameters. A reliable algorithm is also developed for the identifications of weak signals using signal processing techniques. The present work can significantly improve existing SR performances and can have profound practical applications where SR is currently employed for its inherent technological advantages.

On two diffusion neuronal models with multiplicative noise: The mean first-passage time properties

Two diffusion processes with multiplicative noise, able to model the changes in the neuronal membrane depolarization between two consecutive spikes of a single neuron, are considered and compared. The processes have the same deterministic part but different stochastic components. The differences in the state-dependent variabilities, their asymptotic distributions, and the properties of the first-passage time across a constant threshold are investigated. Closed form expressions for the mean of the first-passage time of both processes are derived and applied to determine the role played by the parameters involved in the model. It is shown that for some values of the input parameters, the higher variability, given by the second moment, does not imply shorter mean first-passage time. The reason for that can be found in the complete shape of the stationary distribution of the two processes. Applications outside neuroscience are also mentioned.

Statistical structure of neural spiking under non-Poissonian or other non-white stimulation

Nerve cells in the brain generate sequences of action potentials with a complex statistics. Theoretical attempts to understand this statistics were largely limited to the case of a temporally uncorrelated input (Poissonian shot noise) from the neurons in the surrounding network. However, the stimulation from thousands of other neurons has various sorts of temporal structure. Firstly, input spike trains are temporally correlated because their firing rates can carry complex signals and because of cell-intrinsic properties like neural refractoriness, bursting, or adaptation. Secondly, at the connections between neurons, the synapses, usage-dependent changes in the synaptic weight (short-term plasticity) further shape the correlation structure of the effective input to the cell. From the theoretical side, it is poorly understood how these correlated stimuli, so-called colored noise, affect the spike train statistics. In particular, no standard method exists to solve the associated first-passage-time problem for the interspike-interval statistics with an arbitrarily colored noise. Assuming that input fluctuations are weaker than the mean neuronal drive, we derive simple formulas for the essential interspike-interval statistics for a canonical model of a tonically firing neuron subjected to arbitrarily correlated input from the network. We verify our theory by numerical simulations for three paradigmatic situations that lead to input correlations: (i) rate-coded naturalistic stimuli in presynaptic spike trains; (ii) presynaptic refractoriness or bursting; (iii) synaptic short-term plasticity. In all cases, we find severe effects on interval statistics. Our results provide a framework for the interpretation of firing statistics measured in vivo in the brain.

Addition of visual noise boosts evoked potential-based brain-computer interface

Although noise has a proven beneficial role in brain functions, there have not been any attempts on the dedication of stochastic resonance effect in neural engineering applications, especially in researches of brain-computer interfaces (BCIs). In our study, a steady-state motion visual evoked potential (SSMVEP)-based BCI with periodic visual stimulation plus moderate spatiotemporal noise can achieve better offline and online performance due to enhancement of periodic components in brain responses, which was accompanied by suppression of high harmonics. Offline results behaved with a bell-shaped resonance-like functionality and 7–36% online performance improvements can be achieved when identical visual noise was adopted for different stimulation frequencies. Using neural encoding modeling, these phenomena can be explained as noise-induced input-output synchronization in human sensory systems which commonly possess a low-pass property. Our work demonstrated that noise could boost BCIs in addressing human needs.

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.

Characteristic effects of stochastic oscillatory forcing on neural firing: analytical theory and comparison to paddlefish electroreceptor data

Stochastic signals with pronounced oscillatory components are frequently encountered in neural systems. Input currents to a neuron in the form of stochastic oscillations could be of exogenous origin, e.g. sensory input or synaptic input from a network rhythm. They shape spike firing statistics in a characteristic way, which we explore theoretically in this report. We consider a perfect integrate-and-fire neuron that is stimulated by a constant base current (to drive regular spontaneous firing), along with Gaussian narrow-band noise (a simple example of stochastic oscillations), and a broadband noise. We derive expressions for the nth-order interval distribution, its variance, and the serial correlation coefficients of the interspike intervals (ISIs) and confirm these analytical results by computer simulations. The theory is then applied to experimental data from electroreceptors of paddlefish, which have two distinct types of internal noisy oscillators, one forcing the other. The theory provides an analytical description of their afferent spiking statistics during spontaneous firing, and replicates a pronounced dependence of ISI serial correlation coefficients on the relative frequency of the driving oscillations, and furthermore allows extraction of certain parameters of the intrinsic oscillators embedded in these electroreceptors.

How Sample Paths of Leaky Integrate-and-Fire Models Are Influenced by the Presence of a Firing Threshold

Neural membrane potential data are necessarily conditional on observation being prior to a firing time. In a stochastic leaky integrate-and-fire model, this corresponds to conditioning the process on not crossing a boundary. In the literature, simulation and estimation have almost always been done using unconditioned processes. In this letter, we determine the stochastic differential equations of a diffusion process conditioned to stay below a level S up to a fixed time t1 and of a diffusion process conditioned to cross the boundary for the first time at t1. This allows simulation of sample paths and identification of the corresponding mean process. Differences between the mean of free and conditioned processes are illustrated, as well as the role of noise in increasing these differences.

Survival probability and first-passage-time statistics of a Wiener process driven by an exponential time-dependent drift

The survival probability and the first-passage-time statistics are important quantities in different fields. The Wiener process is the simplest stochastic process with continuous variables, and important results can be explicitly found from it. The presence of a constant drift does not modify its simplicity; however, when the process has a time-dependent component the analysis becomes difficult. In this work we analyze the statistical properties of the Wiener process with an absorbing boundary, under the effect of an exponential time-dependent drift. Based on the backward Fokker-Planck formalism we set the time-inhomogeneous equation and conditions that rule the diffusion of the corresponding survival probability. We propose as the solution an expansion series in terms of the intensity of the exponential drift, resulting in a set of recurrence equations. We explicitly solve the expansion up to second order and comment on higher-order solutions. The first-passage-time density function arises naturally from the survival probability and preserves the proposed expansion. Explicit results, related properties, and limit behaviors are analyzed and extensively compared to numerical simulations.

Distribution of interspike intervals estimated from multiple spike trains observed in a short time window

Several nonparametric estimators of the probability distribution of interspike intervals are introduced. The methods are suitable for simultaneous spike trains observed in a time window of length comparable with the mean interspike interval. This reflects the situation in which a high number of input spike trains converge to a single cortical neuron that has to react in a relatively short time. The simulation study is performed to compare the estimators. For that purpose, several types of stationary point processes are considered as the models of neuronal activity. The methods permit one to estimate the distribution of interspike intervals even if practically none of them are observed. The Kaplan-Meier estimator seems to be the most flexible and reliable among all studied methods, but no direct conclusions as to how real neurons work can be deduced from it.

A first-order nonhomogeneous Markov model for the response of spiking neurons stimulated by small phase-continuous signals

We present a first-order nonhomogeneous Markov model for the interspike-interval density of a continuously stimulated spiking neuron. The model allows the conditional interspike-interval density and the stationary interspike-interval density to be expressed as products of two separate functions, one of which describes only the neuron characteristics and the other of which describes only the signal characteristics. The approximation shows particularly clearly that signal autocorrelations and cross-correlations arise as natural features of the interspike-interval density and are particularly clear for small signals and moderate noise. We show that this model simplifies the design of spiking neuron cross-correlation systems and describe a four-neuron mutual inhibition network that generates a cross-correlation output for two input signals.

Quasistatic approximation of the interspike interval distribution of neurons driven by time-dependent inputs

Variability in neural responses is a ubiquitous phenomenon in neurons, usually modeled with stochastic differential equations. In particular, stochastic integrate-and-fire models are widely used to simplify theoretical studies. The statistical properties of the generated spikes depend on the stimulating input current. Given that real sensory neurons are driven by time-dependent signals, here we study how the interspike interval distribution of integrate-and-fire neurons depends on the evolution of the stimulus in a quasistatic limit. We obtain a closed-form expression for this distribution, and we compare it to the one obtained with numerical simulations for several time-dependent currents. For slow inputs, the quasistatic distribution provides a very good description of the data. The results obtained for the integrate-and-fire model can be extended to other nonautonomous stochastic systems where the first passage time problem has an explicit solution.

Are the input parameters of white noise driven integrate and fire neurons uniquely determined by rate and CV?

Integrate and fire (IF) neurons have found widespread applications in computational neuroscience. Particularly important are stochastic versions of these models where the driving consists of a synaptic input modeled as white Gaussian noise with mean mu and noise intensity D. Different IF models have been proposed, the firing statistics of which depends nontrivially on the input parameters mu and D. In order to compare these models among each other, one must first specify the correspondence between their parameters. This can be done by determining which set of parameters (mu,D) of each model is associated with a given set of basic firing statistics as, for instance, the firing rate and the coefficient of variation (CV) of the interspike interval (ISI). However, it is not clear a priori whether for a given firing rate and CV there is only one unique choice of input parameters for each model. Here we review the dependence of rate and CV on input parameters for the perfect, leaky, and quadratic IF neuron models and show analytically that indeed in these three models the firing rate and the CV uniquely determine the input parameters.

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).

Analytic expressions for rate and CV of a type I neuron driven by white gaussian noise

We study the one-dimensional normal form of a saddle-node system under the influence of additive gaussian white noise and a static "bias current" input parameter, a model that can be looked upon as the simplest version of a type I neuron with stochastic input. This is in contrast with the numerous studies devoted to the noise-driven leaky integrate-and-fire neuron. We focus on the firing rate and coefficient of variation (CV) of the interspike interval density, for which scaling relations with respect to the input parameter and noise intensity are derived. Quadrature formulas for rate and CV are numerically evaluated and compared to numerical simulations of the system and to various approximation formulas obtained in different limiting cases of the model. We also show that caution must be used to extend these results to the Theta neuron model with multiplicative gaussian white noise. The correspondence between the first passage time statistics for the saddle-node model and the Theta neuron model is obtained only in the Stratonovich interpretation of the stochastic Theta neuron model, while previous results have focused only on the Ito interpretation. The correct Stratonovich interpretation yields CVs that are still relatively high, although smaller than in the Ito interpretation; it also produces certain qualitative differences, especially at larger noise intensities. Our analysis provides useful relations for assessing the distance to threshold and the level of synaptic noise in real type I neurons from their firing statistics. We also briefly discuss the effect of finite boundaries (finite values of threshold and reset) on the firing statistics.

Ergodicity of spike trains: when does trial averaging make sense?

Neuronal information processing is often studied on the basis of spiking patterns. The relevant statistics such as firing rates calculated with the peri-stimulus time histogram are obtained by averaging spiking patterns over many experimental runs. However, animals should respond to one experimental stimulation in real situations, and what is available to the brain is not the trial statistics but the population statistics. Consequently, physiological ergodicity, namely, the consistency between trial averaging and population averaging, is implicitly assumed in the data analyses, although it does not trivially hold true. In this letter, we investigate how characteristics of noisy neural network models, such as single neuron properties, external stimuli, and synaptic inputs, affect the statistics of firing patterns. In particular, we show that how high membrane potential sensitivity to input fluctuations, inability of neurons to remember past inputs, external stimuli with large variability and temporally separated peaks, and relatively few contributions of synaptic inputs result in spike trains that are reproducible over many trials. The reproducibility of spike trains and synchronous firing are contrasted and related to the ergodicity issue. Several numerical calculations with neural network examples are carried out to support the theoretical results.

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

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

From subthreshold to firing-rate resonance

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

Input-output behaviour of a model neuron with alternating drift

The input-output behaviour of the Wiener neuronal model subject to alternating input is studied under the assumption that the effect of such an input is to make the drift itself of an alternating type. Firing densities and related statistics are obtained via simulations of the sample-paths of the process in the following three cases: the drift changes occur during random periods characterised by (i) exponential distribution, (ii) Erlang distribution with a preassigned shape parameter, and (iii) deterministic distribution. The obtained results are compared with those holding for the Wiener neuronal model subject to sinusoidal input.

Fluctuation-induced memory retrieval in a pulsed neural network storing sparse patterns with hierarchical correlations

An associative memory in a pulsed neural network composed of the FitzHugh-Nagumo models storing sparse patterns with hierarchical correlations is investigated. The memory patterns composed of 0/1 digits are represented by the synchronous periodic firings of the neurons. It is found that the target pattern and the OR pattern are retrieved individually by controlling the intensity of fluctuations in the system.

Theoretical analysis of array-enhanced stochastic resonance in the diffusively coupled FitzHugh-Nagumo equation

The array-enhanced stochastic resonance (AESR) in the diffusively coupled FitzHugh-Nagumo equation is investigated. The two properties of AESR, namely, the scaling of the optimal noise intensity and the enhancement of the maximum value of the correlation coefficient as a function of the coupling strength, are analyzed theoretically. By transforming the dynamics of N elements into that of the mean and the deviation from it, it is found that AESR is caused by the correlation between them. A low-dimensional model that reproduces the above properties is constructed.

Stochastic dynamics in a two-dimensional oscillator near a saddle-node bifurcation

We study the oscillator equations describing a particular class of nonlinear amplifier, exemplified in this work by a two-junction superconducting quantum interference device. This class of dynamic system is described by a potential energy function that can admit minima (corresponding to stable solutions of the dynamic equations), or "running states" wherein the system is biased so that the potential minima disappear and the solutions display spontaneous oscillations. Just beyond the onset of the spontaneous oscillations, the system is known to show significantly enhanced sensitivity to very weak magnetic signals. The global phase space structure allows us to apply a center manifold technique to approximate analytically the oscillatory behavior just past the (saddle-node) bifurcation and compute the oscillation period, which obeys standard scaling laws. In this regime, the dynamics can be represented by an "integrate-fire" model drawn from the computational neuroscience repertoire; in fact, we obtain an "interspike interval" probability density function and an associated power spectral density (computed via Renewal theory) that agree very well with the results obtained via numerical simulations. Notably, driving the system with one or more time sinusoids produces a noise-lowering injection locking effect and/or heterodyning.

Noninvasive control of stochastic resonance

External feedback can enhance (or depress) the response of a noisy bistable system to monochromatic signals, significantly magnifying its natural stochastic resonance. We compare and contrast a variety of such feedback strategies, using both numerical simulations and analog electronic experiments. These noninvasive control techniques are especially valuable for noisy bistable systems that are difficult or impossible to modify internally.

Stochastic resonance in neuron models: endogenous stimulation revisited

The paradigm of stochastic resonance (SR)-the idea that signal detection and transmission may benefit from noise-has met with great interest in both physics and the neurosciences. We investigate here the consequences of reducing the dynamics of a periodically driven neuron to a renewal process (stimulation with reset or endogenous stimulation). This greatly simplifies the mathematical analysis, but we show that stochastic resonance as reported earlier occurs in this model only as a consequence of the reduced dynamics.

Effect of spatial extension on noise-enhanced phase locking in a leaky integrate-and-fire model of a neuron

Signal transmission enhanced by noise has been recently investigated in detail on the single compartment, also referred to as single point, leaky integrate-and-fire model neuron under a subthreshold stimulation. In this paper we study how this phenomenon is influenced by taking into account the spatial characteristics of the neuron. A stochastic two-point leaky integrate-and-fire model, comprising a dendritic compartment and trigger zone, under periodic stimulation is studied. A method of how to measure synchronization between the signal and the output in both, experiments and models, is proposed. This method is based on a distance between the exact periodic spiking, as expected for sufficiently strong and noiseless stimulation, and neuronal activity evoked by a subthreshold signal corrupted by noise. It is shown that qualitatively the same phenomenon, phase-locking enhanced by the noise, as found in the spatially unstructured neuron is produced by the spatially complex neuron. However, quantitatively there are significant differences. Namely, the two-point model neuron is more robust against the noise and therefore its amplitude has to be higher to enhance the signal. Further, it is found that the range of the critical levels of noise is larger for the two-point model than for the single-point one. Finally, the enhancing effect at the optimal noise is more efficient in the single-point model and thus the firing patterns at their optimal noise levels are different in both models.

A simple stochastic model of spatially complex neurons

A method for studying the coding properties of a multicompartmental integrate-and-fire neuron of arbitrary geometry is presented. Depolarization at each compartment evolves like a leaky integrator with an after-firing reset imposed only at the trigger zone. The frequency of firing at the steady-state regime is related to the properties of the multidimensional input. The decreasing variability of subthreshold depolarization from the dendritic tree to the trigger zone is shown for an input that is corrupted by a white noise. The role of a Poissonian noise is also investigated. The proposed method gives an estimate of the mean interspike interval that can be used to study the input output transfer function of the system. Both types of the stochastic inputs result in broadening the transfer function with respect to the deterministic case.

The variability of basal action potential firing is positively correlated with bursting in hypothalamic oxytocin neurones

Magnocellular oxytocin neurones are proposed as a suitable system for studying the mechanisms involved in the regulation of neuronal bursting activity. They display high frequency (50 sp./s) bursts of spikes (approximately every 300 s), in response to specific stimuli, which are superimposed on a variable level of basal activity and are tightly co-ordinated as a result of network interactions. The relationship between the strength of the bursting activity (as quantified by burst amplitude and interburst interval) and the characteristics of the interburst basal activity were assessed. During control conditions, mean basal activity and variability of firing increased just before bursts. During experimental conditions leading to burst facilitation, burst amplitude increased and interburst interval decreased while a sustained increase in mean firing rate occurred. Variability of firing (measured by both the standard deviation of firing rate, and the index of dispersion which corrected this standard deviation for differences in mean firing rate), increased demonstrating an increase in spike clustering greater than expected as a result of increased basal activity. When bursting was restrained (i.e. interburst interval increased), mean basal activity increased substantially, but index of dispersion decreased. A narrowing of the interspike interval distribution occurred, indicating increased regularity of firing. The aspect of basal activity most strongly correlated with bursting was variability of firing rate. The strongest correlate of burst amplitude was the standard deviation of mean firing rate, whereas the strongest and most consistent correlate of interburst interval was the index of dispersion. In conclusion, bursting behaviour is most strongly related to the irregularity rather than the level of basal activity.

Coherence and stochastic resonance in a two-state system

The subject of our study is a two-state dynamics driven by Gaussian white noise and a weak harmonic signal. The system resulting from a piecewise linear FizHugh-Nagumo model in the case of perfect time scale separation between fast and slow variables shows either bistable, excitable, or oscillatory behavior. Its output spectra as well as the spectral power amplification of the signal can be calculated for arbitrary noise strength and frequency, allowing characterization of the coherence resonance in the bistable and excitable regimes as well as quantification of nonadiabatic resonances with respect to the external signal in all regimes.

Stochastic resonance in one-dimensional diffusion with one reflecting and one absorbing end point

An analysis of the nonmonotonic dependence of the mean-free-passage time on the frequency of a periodic signal [stochastic resonance (SR)] for diffusion on a segment with one absorbing and one reflecting end point shows that SR exists only for some restricted values of parameters. SR always exists if the periodic telegraph signal is replaced by a random one. The latter case is considered in detail.

Statistical analysis of stochastic resonance in a simple setting

A subthreshold signal may be detected if noise is added to the data. We study a simple model, consisting of a constant signal to which at uniformly spaced times independent and identically distributed noise variables with known distribution are added. A detector records the times at which the noisy signal exceeds a threshold. There is an optimal noise level, called stochastic resonance. We explore the detectability of the signal in a system with one or more detectors, with different thresholds. We use a statistical detectability measure, the asymptotic variance of the best estimator of the signal from the thresholded data, or equivalently, the Fisher information in the data. In particular, we determine optimal configurations of detectors, varying the distances between the thresholds and the signal, as well as the noise level. The approach generalizes to nonconstant signals.

Mean discharge frequency locking in the response of a noisy neuron model to subthreshold periodic stimulation

Leaky integrate-and-fire neuron models display stochastic resonance-like behavior when stimulated by subthreshold periodic signal and noise. Previous works have shown that matching between the time scales of the noise induced discharges and the modulation period can account for this phenomenon at low modulation amplitudes, but not large subthreshold modulation amplitude. In order to examine the discharge patterns of the model in this regime, we introduce a method for the computation of the power spectral density of the discharge train. Using this method, we clarify the role of the distribution of the input phase at discharge times. Finally, we argue that for large subthreshold inputs, mean discharge frequency locking accounts for the enhanced response.

Stochastic resonance in thermally activated reactions: Application to biological ion channels

At the molecular level many thermally activated reactions can be viewed as Poisson trains of events whose instantaneous rates are defined by the reaction activation barrier height and an effective collision frequency. When the barrier height depends on an external parameter, variation in this parameter induces variation in the event rate. Extending our previous work, we offer a detailed theoretical analysis of signal transduction properties of these reactions considering the external parameter as an input signal and the train of resulting events as an output signal. The addition of noise to the system input facilitates signal transduction in two ways. First, for a linear relationship between the barrier height and the external parameter the output signal power grows exponentially with the mean square fluctuation of the noise. Second, for noise of a sufficiently high bandwidth, its addition increases output signal quality measured as the signal-to-noise ratio (SNR). The output SNR reaches a maximum at optimal noise intensity defined by the reaction sensitivity to the external parameter, reaction initial rate, and the noise bandwidth. We apply this theory to ion channels of excitable biological membranes. Based on classical results of Hodgkin and Huxley we show that open/closed transitions of voltage-gated ion channels can be treated as thermally activated reactions whose activation barriers change linearly with applied transmembrane voltage. As an experimental example we discuss our recent results obtained with polypeptide alamethicin incorporated into planar lipid bilayers.(c) 1998 American Institute of Physics.

Minireview of stochastic resonance

We present an introductory overview of the subject of stochastic resonance. As researchers' interest in the phenomenon has spread from physics to biology, new questions both fundamental and practical have emerged. After reviewing some key aspects of the subject, we describe a promising candidate for exploring the possible beneficial effects of random noise in sensory transduction. (c) 1998 American Institute of Physics.

Stochastic resonance in mammalian neuronal networks

We present stochastic resonance observed in the dynamics of neuronal networks from mammalian brain. Both sinusoidal signals and random noise were superimposed into an applied electric field. As the amplitude of the noise component was increased, an optimization (increase then decrease) in the signal-to-noise ratio of the network response to the sinusoidal signal was observed. The relationship between the measures used to characterize the dynamics is discussed. Finally, a computational model of these neuronal networks that includes the neuronal interactions with the electric field is presented to illustrate the physics behind the essential features of the experiment. (c) 1998 American Institute of Physics.

Signal transduction across alamethicin ion channels in the presence of noise

We have studied voltage-dependent ion channels of alamethicin reconstituted into an artificial planar lipid bilayer membrane from the point of view of electric signal transduction. Signal transduction properties of these channels are highly sensitive to the external electric noise. Specifically, addition of bandwidth-restricted "white" noise of 10-20 mV (r.m.s.) to a small sine wave input signal increases the output signal by approximately 20-40 dB conserving, and even slightly increasing, the signal-to-noise ratio at the system output. We have developed a small-signal adiabatic theory of stochastic resonance for a threshold-free system of voltage-dependent ion channels. This theory describes our main experimental findings giving good qualitative understanding of the underlying mechanism. It predicts the right value of the output signal-to-noise ratio and provides a reliable estimate for the noise intensity corresponding to its maximum. Our results suggest that the alamethicin channel in a lipid bilayer is a good model system for studies of mechanisms of primary electrical signal processing in biology showing an important feature of signal transduction improvement by a fluctuating environment.