Complementary responses to mean and variance modulations in the perfect integrate-and-fire model

In the perfect integrate-and-fire model (PIF), the membrane voltage is proportional to the integral of the input current since the time of the previous spike. It has been shown that the firing rate within a noise free ensemble of PIF neurons responds instantaneously to dynamic changes in the input current, whereas in the presence of white noise, model neurons preferentially pass low frequency modulations of the mean current. Here, we prove that when the input variance is perturbed while holding the mean current constant, the PIF responds preferentially to high frequency modulations. Moreover, the linear filters for mean and variance modulations are complementary, adding exactly to one. Since changes in the rate of Poisson distributed inputs lead to proportional changes in the mean and variance, these results imply that an ensemble of PIF neurons transmits a perfect replica of the time-varying input rate for Poisson distributed input. A more general argument shows that this property holds for any signal leading to proportional changes in the mean and variance of the input current.     

State-Dependent Effects of Na Channel Noise on Neuronal Burst Generation

We explore the effects of stochastic sodium (Na) channel activation on the variability and dynamics of spiking and bursting in a model neuron. The complete model segregates Hodgin-Huxley-type currents into two compartments, and undergoes applied current-dependent bifurcations between regimes of periodic bursting, chaotic bursting, and tonic spiking. Noise is added to simulate variable, finite sizes of the population of Na channels in the fast spiking compartment.During tonic firing, Na channel noise causes variability in interspike intervals (ISIs). The variance, as well as the sensitivity to noise, depend on the model's biophysical complexity. They are smallest in an isolated spiking compartment; increase significantly upon coupling to a passive compartment; and increase again when the second compartment also includes slow-acting currents. In this full model, sufficient noise can convert tonic firing into bursting.During bursting, the actions of Na channel noise are state-dependent. The higher the noise level, the greater the jitter in spike timing within bursts. The noise makes the burst durations of periodic regimes variable, while decreasing burst length duration and variance in a chaotic regime. Na channel noise blurs the sharp transitions of spike time and burst length seen at the bifurcations of the noise-free model. Close to such a bifurcation, the burst behaviors of previously periodic and chaotic regimes become essentially indistinguishable.We discuss biophysical mechanisms, dynamical interpretations and physiological implications. We suggest that noise associated with finite populations of Na channels could evoke very different effects on the intrinsic variability of spiking and bursting discharges, depending on a biological neuron's complexity and applied current-dependent state. We find that simulated channel noise in the model neuron qualitatively replicates the observed variability in burst length and interburst interval in an isolated biological bursting neuron.     

Noise-induced divisive gain control in neuron models

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

A note on neuronal firing and input variability.

The Ornstein-Uhlenbeck process with a constant forcing function has often been used as a model for the subthreshold membrane potential of a neuron. The mean, variance and coefficient of variation of the first passage time to a constant threshold are examined for this model in the limit of small synaptic noise and low thresholds. A comparison is made between the asymptotic results of Wan & Tuckwell, who used perturbation analysis, and several computationally simpler approximation methods. A generalization of Stein's method gives an overestimate of the mean interval while an approximation by a Wiener process with linear drift gives an underestimate of the mean interval. These bounds are simple to calculate and can be used as a prelude to a more detailed perturbation analysis.     

Increased Computational Accuracy in Multi-Compartmental Cable Models by a Novel Approach for Precise Point Process Localization

Compartmental models of dendrites are the most widely used tool for investigating their electrical behaviour. Traditional models assign a single potential to a compartment. This potential is associated with the membrane potential at the centre of the segment represented by the compartment. All input to that segment, independent of its location on the segment, is assumed to act at the centre of the segment with the potential of the compartment. By contrast, the compartmental model introduced in this article assigns a potential to each end of a segment, and takes into account the location of input to a segment on the model solution by partitioning the effect of this input between the axial currents at the proximal and distal boundaries of segments. For a given neuron, the new and traditional approaches to compartmental modelling use the same number of locations at which the membrane potential is to be determined, and lead to ordinary differential equations that are structurally identical. However, the solution achieved by the new approach gives an order of magnitude better accuracy and precision than that achieved by the latter in the presence of point process input.Action Editor: Alain Destexhe     

Transient Responses to Rapid Changes in Mean and Variance in Spiking Models

The mean input and variance of the total synaptic input to a neuron can vary independently, suggesting two distinct information channels. Here we examine the impact of rapidly varying signals, delivered via these two information conduits, on the temporal dynamics of neuronal firing rate responses. We examine the responses of model neurons to step functions in either the mean or the variance of the input current. Our results show that the temporal dynamics governing response onset depends on the choice of model. Specifically, the existence of a hard threshold introduces an instantaneous component into the response onset of a leaky-integrate-and-fire model that is not present in other models studied here. Other response features, for example a decaying oscillatory approach to a new steady-state firing rate, appear to be more universal among neuronal models. The decay time constant of this approach is a power-law function of noise magnitude over a wide range of input parameters. Understanding how specific model properties underlie these response features is important for understanding how neurons will respond to rapidly varying signals, as the temporal dynamics of the response onset and response decay to new steady-state determine what range of signal frequencies a population of neurons can respond to and faithfully encode.     

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.     

Intrinsic gain modulation and adaptive neural coding

In many cases, the computation of a neural system can be reduced to a receptive field, or a set of linear filters, and a thresholding function, or gain curve, which determines the firing probability; this is known as a linear/nonlinear model. In some forms of sensory adaptation, these linear filters and gain curve adjust very rapidly to changes in the variance of a randomly varying driving input. An apparently similar but previously unrelated issue is the observation of gain control by background noise in cortical neurons: the slope of the firing rate versus current (f-I) curve changes with the variance of background random input. Here, we show a direct correspondence between these two observations by relating variance-dependent changes in the gain of f-I curves to characteristics of the changing empirical linear/nonlinear model obtained by sampling. In the case that the underlying system is fixed, we derive relationships relating the change of the gain with respect to both mean and variance with the receptive fields derived from reverse correlation on a white noise stimulus. Using two conductance-based model neurons that display distinct gain modulation properties through a simple change in parameters, we show that coding properties of both these models quantitatively satisfy the predicted relationships. Our results describe how both variance-dependent gain modulation and adaptive neural computation result from intrinsic nonlinearity.     

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     

Effect of Voltage Sensitive Fluorescent Proteins on Neuronal Excitability

Fluorescent protein voltage sensors are recombinant proteins that are designed as genetically encoded cellular probes of membrane potential using mechanisms of voltage-dependent modulation of fluorescence. Several such proteins, including VSFP2.3 and VSFP3.1, were recently reported with reliable function in mammalian cells. They were designed as molecular fusions of the voltage sensor of Ciona intestinalis voltage sensor containing phosphatase with a fluorescence reporter domain. Expression of these proteins in cell membranes is accompanied by additional dynamic membrane capacitance, or “sensing capacitance”, with feedback effect on the native electro-responsiveness of targeted cells. We used recordings of sensing currents and fluorescence responses of VSFP2.3 and of VSFP3.1 to derive kinetic models of the voltage-dependent signaling of these proteins. Using computational neuron simulations, we quantitatively investigated the perturbing effects of sensing capacitance on the input/output relationship in two central neuron models, a cerebellar Purkinje and a layer 5 pyramidal neuron. Probe-induced sensing capacitance manifested as time shifts of action potentials and increased synaptic input thresholds for somatic action potential initiation with linear dependence on the membrane density of the probe. Whereas the fluorescence signal/noise grows with the square root of the surface density of the probe, the growth of sensing capacitance is linear. We analyzed the trade-off between minimization of sensing capacitance and signal/noise of the optical read-out depending on kinetic properties and cellular distribution of the probe. The simulation results suggest ways to reduce capacitive effects at a given level of signal/noise. Yet, the simulations indicate that significant improvement of existing probes will still be required to report action potentials in individual neurons in mammalian brain tissue in single trials.     

Linear n-compartment catenary models: Formulas to describe tracer amount in any compartment and identification of parameters from a concentration-time curve

Resolution of kinetic equations and parameter identification are discussed for n-compartment linear catenary models with elimination allowed from any compartment. For a given input, general formulas are derived to describe the tracer amount in any compartment as a function of the model parameters. Conversely, explicit procedures are given to identify the model parameters when the concentration-time curve is known in one arbitrary compartment, the tracer being injected into the same compartment. In this inverse problem, the solution is not unique: the model transfer rate constants can only be localized in a finite set of intervals.     

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 firstorder neurons as a function of the logarithm of stimulus concentration.     

The spatial properties of a model neuron increase its coding range

The coding properties of one-compartment and two-compartment model neurons are compared. The membrane depolarization in both models is described as a deterministic leaky integrator. Interspike intervals are identified with the periods between reset of the depolarization after firing and consecutive crossing of a fixed firing threshold. The two-point model has an input in the dendritic compartment and an output in the trigger-zone compartment. It is shown that the sensitivity threshold for the two-point model is shifted to the larger values of the input intensity with respect to the sensitivity threshold of its single-point counterpart. Further, its coding range is substantially larger than the coding range of the single-point model.     

Random currents through nerve membranes

The linear cable equation with uniform Poisson or white noise input current is employed as a model for the voltage across the membrane of a onedimensional nerve cylinder, which may sometimes represent the dendritic tree of a nerve cell. From the Green's function representation of the solutions, the mean, variance and covariance of the voltage are found. At large times, the voltage becomes asymptotically wide-sense stationary and we find the spectral density functions for various cable lengths and boundary conditions. For large frequencies the voltage exhibits “1/f ^3/2 noise”. Using the Fourier series representation of the voltage we study the moments of the firing times for the diffusion model with numerical techniques, employing a simplified threshold criterion. We also simulate the solution of the stochastic cable equation by two different methods in order to estimate the moments and density of the firing time.     

Firing-rate models capture essential response dynamics of LGN relay cells

Firing-rate models provide a practical tool for studying signal processing in the early visual system, permitting more thorough mathematical analysis than spike-based models. We show here that essential response properties of relay cells in the lateral geniculate nucleus (LGN) can be captured by surprisingly simple firing-rate models consisting of a low-pass filter and a nonlinear activation function. The starting point for our analysis are two spiking neuron models based on experimental data: a spike-response model fitted to data from macaque (Carandini et al. J. Vis., 20(14), 1–2011, 2007), and a model with conductance-based synapses and afterhyperpolarizing currents fitted to data from cat (Casti et al. J. Comput. Neurosci., 24(2), 235–252, 2008). We obtained the nonlinear activation function by stimulating the model neurons with stationary stochastic spike trains, while we characterized the linear filter by fitting a low-pass filter to responses to sinusoidally modulated stochastic spike trains. To account for the non-Poisson nature of retinal spike trains, we performed all analyses with spike trains with higher-order gamma statistics in addition to Poissonian spike trains. Interestingly, the properties of the low-pass filter depend only on the average input rate, but not on the modulation depth of sinusoidally modulated input. Thus, the response properties of our model are fully specified by just three parameters (low-frequency gain, cutoff frequency, and delay) for a given mean input rate and input regularity. This simple firing-rate model reproduces the response of spiking neurons to a step in input rate very well for Poissonian as well as for non-Poissonian input. We also found that the cutoff frequencies, and thus the filter time constants, of the rate-based model are unrelated to the membrane time constants of the underlying spiking models, in agreement with similar observations for simpler models.     

Axonal speeding: shaping synaptic potentials in small neurons by the axonal membrane compartment

The role of the axonal membrane compartment in synaptic integration is usually neglected. We show here that in interneurons of the cerebellar molecular layer, where dendrites are so short that the somatodendritic domain can be considered isopotential, the axonal membrane contributes a significant part of the cell input capacitance. We examine the impact of axonal membrane on synaptic integration by cutting the axon with two-photon illumination. We find that the axonal compartment acts as a sink for signals generated at fast conductance synapses, thus increasing the initial decay rate of corresponding synaptic potentials over the value predicted from the resistance-capacitance (RC) product of the cell membrane; signals generated at slower synapses are much less affected. This mechanism sharpens the spike firing precision of fast glutamatergic inputs without resorting to multisynaptic pathways.     

A review of the methods for signal estimation in stochastic diffusion leaky integrate-and-fire neuronal models

Parameters in diffusion neuronal models are divided into two groups; intrinsic and input parameters. Intrinsic parameters are related to the properties of the neuronal membrane and are assumed to be known throughout the paper. Input parameters characterize processes generated outside the neuron and methods for their estimation are reviewed here. Two examples of the diffusion neuronal model, which are based on the integrate-and-fire concept, are investigated--the Ornstein--Uhlenbeck model as the most common one and the Feller model as an illustration of state-dependent behavior in modeling the neuronal input. Two types of experimental data are assumed-intracellular describing the membrane trajectories and extracellular resulting in knowledge of the interspike intervals. The literature on estimation from the trajectories of the diffusion process is extensive and thus the stress in this review is set on the inference made from the interspike intervals.     

A Statistical Model for In Vivo Neuronal Dynamics

Single neuron models have a long tradition in computational neuroscience. Detailed biophysical models such as the Hodgkin-Huxley model as well as simplified neuron models such as the class of integrate-and-fire models relate the input current to the membrane potential of the neuron. Those types of models have been extensively fitted to in vitro data where the input current is controlled. Those models are however of little use when it comes to characterize intracellular in vivo recordings since the input to the neuron is not known. Here we propose a novel single neuron model that characterizes the statistical properties of in vivo recordings. More specifically, we propose a stochastic process where the subthreshold membrane potential follows a Gaussian process and the spike emission intensity depends nonlinearly on the membrane potential as well as the spiking history. We first show that the model has a rich dynamical repertoire since it can capture arbitrary subthreshold autocovariance functions, firing-rate adaptations as well as arbitrary shapes of the action potential. We then show that this model can be efficiently fitted to data without overfitting. We finally show that this model can be used to characterize and therefore precisely compare various intracellular in vivo recordings from different animals and experimental conditions.     

Stochastic Dynamic Model for Estimation of Rate Constants and Their Variances from Noisy and Heterogeneous PET Measurements

Tissue heterogeneity, radioactive decay and measurement noise are the main error sources in compartmental modeling used to estimate the physiologic rate constants of various radiopharmaceuticals from a dynamic PET study. We introduce a new approach to this problem by modeling the tissue heterogeneity with random rate constants in compartment models. In addition, the Poisson nature of the radioactive decay is included as a Poisson random variable in the measurement equations. The estimation problem will be carried out using the maximum likelihood estimation. With this approach, we do not only get accurate mean estimates for the rate constants, but also estimates for tissue heterogeneity within the region of interest and other possibly unknown model parameters, e.g. instrument noise variance, as well. We also avoid the problem of the optimal weighting of the data related to the conventionally used weighted least-squares method. The new approach was tested with simulated time–activity curves from the conventional three compartment – three rate constants model with normally distributed rate constants and with a noise mixture of Poisson and normally distributed random variables. Our simulation results showed that this new model gave accurate estimates for the mean of the rate constants, the measurement noise parameter and also for the tissue heterogeneity, i.e. for the variance of the rate constants within the region of interest.