A stochastic model for interacting neurons in the olfactory bulb.

We focus on interacting neurons organized in a block-layered network devoted to the information processing from the sensory system to the brain. Specifically, we consider the firing activity of olfactory sensory neurons, periglomerular, granule and mitral cells in the context of the neuronal activity of the olfactory bulb. We propose and investigate a stochastic model of a layered and modular network to describe the dynamic behavior of each prototypical neuron, taking into account both its role (excitatory/inhibitory) and its location within the network. We adopt specific Gauss-Markov processes suitable to provide reliable estimates of the firing activity of the different neurons, given their linkages. Furthermore, we study the impact of selective excitation/inhibition on the information transmission by means of simulations and numerical estimates obtained through a Volterra integral approach.

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.

Linked Gauss-Diffusion processes for modeling a finite-size neuronal network.

A Leaky Integrate-and-Fire (LIF) model with stochastic current-based linkages is considered to describe the firing activity of neurons interacting in a (2×2)-size feed-forward network. In the subthreshold regime and under the assumption that no more than one spike is exchanged between coupled neurons, the stochastic evolution of the neuronal membrane voltage is subject to random jumps due to interactions in the network. Linked Gauss-Diffusion processes are proposed to describe this dynamics and to provide estimates of the firing probability density of each neuron. To this end, an iterated integral equation-based approach is applied to evaluate numerically the first passage time density of such processes through the firing threshold. Asymptotic approximations of the firing densities of surrounding neurons are used to obtain closed-form expressions for the mean of the involved processes and to simplify the numerical procedure. An extension of the model to an (N×N)-size network is also given. Histograms of firing times obtained by simulations of the LIF dynamics and numerical firings estimates are compared.

On a stochastic leaky integrate-and-fire neuronal model.

The leaky integrate-and-fire neuronal model proposed in Stevens and Zador (1998), in which time constant and resting potential are postulated to be time dependent, is revisited within a stochastic framework in which the membrane potential is mathematically described as a gauss-diffusion process. The first-passage-time probability density, miming in such a context the firing probability density, is evaluated by either the Volterra integral equation of Buonocore, Nobile, and Ricciardi ( 1987 ) or, when possible, by the asymptotics of Giorno, Nobile, and Ricciardi (1990). The model examined here represents an extension of the classic leaky integrate-and-fire one based on the Ornstein-Uhlenbeck process in that it is in principle compatible with the inclusion of some other physiological characteristics such as relative refractoriness. It also allows finer tuning possibilities in view of its accounting for certain qualitative as well as quantitative features, such as the behavior of the time course of the membrane potential prior to firings and the computation of experimentally measurable statistical descriptors of the firing time: mean, median, coefficient of variation, and skewness. Finally, implementations of this model are provided in connection with certain experimental evidence discussed in the literature.

On a pulsating Brownian motor and its characterization.

As a model of Brownian motor we consider the jump diffusion motion of a particle in the presence of an asymmetric periodic potential with a unique minimum and subject to half-period space shifts at the instants of occurrence of two Poisson processes. The relevant quantities, i.e., probability current, effective driving force, stall force, power and efficiency of the motor are explicitly calculated as averages of certain functions of the random variable representing the particle position.

On Myosin II dynamics in the presence of external loads.

We address the controversial hot question concerning the validity of the loose coupling versus the lever-arm theories in the actomyosin dynamics by re-interpreting and extending the phenomenological washboard potential model proposed by some of us in a previous paper. In this new model a Brownian motion harnessing thermal energy is assumed to co-exist with the deterministic swing of the lever-arm, to yield an excellent fit of the set of data obtained by some of us on the sliding of Myosin II heads on immobilized actin filaments under various load conditions. Our theoretical arguments are complemented by accurate numerical simulations, and the robustness of the model is tested via different choices of parameters and potential profiles.

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

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

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.

On a non-Markov neuronal model and its approximations.

Single neuron's activity modeling is considered with reference to some earlier contributions in which a non-Markov Gaussian process is assumed to describe the time course of the neuron's membrane potential. After re-formulating the problem in a rigorous framework and pinpointing the limits of validity of such a model, the available results on the firing probability density are compared with those obtained by us by means of an ad hoc numerical algorithm implemented for the leaky integrator diffusion firing model and with some data constructed by a simulation procedure of non-Markov Gaussian processes with pre-assigned covariances. Throughout this paper, the notion of 'correlation time' plays a fundamental role for the neuronal coding process modeling.