Enhanced Sensitivity to Rapid Input Fluctuations by Nonlinear Threshold Dynamics in Neocortical Pyramidal Neurons

The way in which single neurons transform input into output spike trains has fundamental consequences for network coding. Theories and modeling studies based on standard Integrate-and-Fire models implicitly assume that, in response to increasingly strong inputs, neurons modify their coding strategy by progressively reducing their selective sensitivity to rapid input fluctuations. Combining mathematical modeling with in vitro experiments, we demonstrate that, in L5 pyramidal neurons, the firing threshold dynamics adaptively adjust the effective timescale of somatic integration in order to preserve sensitivity to rapid signals over a broad range of input statistics. For that, a new Generalized Integrate-and-Fire model featuring nonlinear firing threshold dynamics and conductance-based adaptation is introduced that outperforms state-of-the-art neuron models in predicting the spiking activity of neurons responding to a variety of in vivo-like fluctuating currents. Our model allows for efficient parameter extraction and can be analytically mapped to a Generalized Linear Model in which both the input filter—describing somatic integration—and the spike-history filter—accounting for spike-frequency adaptation—dynamically adapt to the input statistics, as experimentally observed. Overall, our results provide new insights on the computational role of different biophysical processes known to underlie adaptive coding in single neurons and support previous theoretical findings indicating that the nonlinear dynamics of the firing threshold due to Na⁺-channel inactivation regulate the sensitivity to rapid input fluctuations.     

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.     

Fast negative feedback enables mammalian auditory nerve fibers to encode a wide dynamic range of sound intensities

Mammalian auditory nerve fibers (ANF) are remarkable for being able to encode a 40 dB, or hundred fold, range of sound pressure levels into their firing rate. Most of the fibers are very sensitive and raise their quiescent spike rate by a small amount for a faint sound at auditory threshold. Then as the sound intensity is increased, they slowly increase their spike rate, with some fibers going up as high as ～300 Hz. In this way mammals are able to combine sensitivity and wide dynamic range. They are also able to discern sounds embedded within background noise. ANF receive efferent feedback, which suggests that the fibers are readjusted according to the background noise in order to maximize the information content of their auditory spike trains. Inner hair cells activate currents in the unmyelinated distal dendrites of ANF where sound intensity is rate-coded into action potentials. We model this spike generator compartment as an attenuator that employs fast negative feedback. Input current induces rapid and proportional leak currents. This way ANF are able to have a linear frequency to input current (f-I) curve that has a wide dynamic range. The ANF spike generator remains very sensitive to threshold currents, but efferent feedback is able to lower its gain in response to noise.     

Firing rates of neurons with random excitation and inhibition

The expectation of the interspike interval for a Stein model neuron receiving Poisson excitation and inhibition is determined by solving a differential difference equation with both forward and backward differences. The method of solution relies on an asymptotic expansion at large initial hyperpolarizations. The asymptotic solution is continued to near threshold depolarization whereupon the boundary condition is employed along with recursion relations to obtain the complete solution. The dependency of the mean firing rate on excitation at fixed inhibition and on inhibition at fixed excitation is investigated as well as the threshold dependence at fixed input rates. The results are discussed in relation to those for intracellular current injection and synaptic input to real neurons.     

神经元放电阈值的可变性及其意义

神经元能够将不同时空模式的突触输入转化为时序精确的动作电位输出，这种灵活、可靠的信息编码方式是神经集群在动态环境或特定任务下产生所需活动模式的重要基础。动作电位的产生遵循全或无规律，只有当细胞膜电压达到放电阈值时，神经元才产生动作电位。放电阈值在细胞内和细胞间具有高度可变性，具体动态依赖于刺激输入和放电历史。特别是，放电阈值对动作电位起始前的膜电压变化十分敏感，这种状态依赖性产生的生物物理根源包括Na⁺失活和K⁺激活。在绝大多数神经元中，动作电位的触发位置是轴突起始端，这个位置处的阈值可变性是决定神经元对时空输入转化规律的关键因素。但是，电生理实验中动作电位的记录位置却通常是胞体或近端树突，此处的阈值可变性高于轴突起始端，而其产生的重要根源是轴突动作电位的反向传播。基于胞体测量的相关研究显示，放电阈值动态能够增强神经元的时间编码、特征选择、增益调控和同时侦测能力。本文首先介绍放电阈值的概念及量化方法，然后详细梳理近年来国内外关于放电阈值可变性及产生根源的研究进展，在此基础上归纳总结放电阈值可变性对神经元编码的重要性，最后对未来放电阈值的研究方向进行展望。     

Ornstein-Uhlenbeck model neuron revisited

 The temporal patterns of action potentials fired by a two-point stochastic neuron model were investigated. In this model the membrane potential of the dendritic compartment follows the Orstein-Uhlenbeck process and is not affected by the spiking activity. The axonal compartment, corresponding to the spike initiation site, is described by a simplified RC circuit. Estimators of the mean and variance of the input, based on a sampling of the axonal membrane potential, were derived and applied to simulated data. The dependencies of the mean firing frequency and of the coefficient of variation and serial correlation of interspike intervals on the mean and variance of the input were also studied by computer simulation in both 1- and 2-point models. The main property distinguishing the 2-point model from the classical 1-point model is its ability to produce clusters of short (or long) intervals between spikes under conditions of constant stimulation, as often observed in real neurons. It is shown that the nearly linear frequency response of the neuron, starting with subthreshold values of the input, is accounted for by the variability of the input (noise), which indicates that noise can play a positive role in nervous systems. The linear response frequency with respect to noise of the models suggests that the neuron can function as a noise encoder. Received: 2 April 1993/Accepted in revised form: 15 September 1994     

Improved neuronal models for studying neural networks

Previous neuronal models used for the study of neural networks are considered. Equations are developed for a model which includes: 1) a normalized range of firing rates with decreased sensitivity at large excitatory or large inhibitory input levels, 2) a single rate constant for the increase in firing rate following step changes in the input, 3) one or more rate constants, as required to fit experimental data for the adaptation of firing rates to maintained inputs. Computed responses compare well with the types of neuronal responses observed experimentally. Depending on the parameters, overdamped increases and decreases, damped oscillatory or maintained oscillatory changes in firing rate are observed to step changes in the input. The integrodifferential equations describing the neuronal models can be represented by a set of first-order differential equations. Steady-state solutions for these equations can be obtained for constant inputs, as well as the stability of the solutions to small perturbations. The linear frequency response function is derived for sufficiently small time-varying inputs. The linear responses are also compared with the computed solutions for larger non-linear responses.     

Computation of spiking activity for a stochastic spatial neuron model: effects of spatial distribution of input on bimodality and CV of the ISI distribution

We obtain computational results for a new extended spatial neuron model in which the neuronal electrical depolarization from resting level satisfies a cable partial differential equation and the synaptic input current is also a function of space and time, obeying a first order linear partial differential equation driven by a two-parameter random process. The model is first described explicitly with the inclusion of all biophysical parameters. Simplified equations are obtained with dimensionless space and time variables. A standard parameter set is described, based mainly on values appropriate for cortical pyramidal cells. When the noise is small and the mean voltage crosses threshold, a formula is derived for the expected time to spike. A simulation algorithm, involving one-dimensional random processes is given and used to obtain moments and distributions of the interspike interval (ISI). The parameters used are those for a near balanced state and there is great sensitivity of the firing rate around the balance point. This sensitivity may be related to genetically induced pathological brain properties (Rett's syndrome). The simulation procedure is employed to find the ISI distribution for some simple patterns of synaptic input with various relative strengths for excitation and inhibition. With excitation only, the ISI distribution is unimodal of exponential type and with a large coefficient of variation. As inhibition near the soma grows, two striking effects emerge. The ISI distribution shifts first to bimodal and then to unimodal with an approximately Gaussian shape with a concentration at large intervals. At the same time the coefficient of variation of the ISI drops dramatically to less than 1/5 of its value without inhibition.     

On the mechanisms underlying the depolarization block in the spiking dynamics of CA1 pyramidal neurons

Under sustained input current of increasing strength neurons eventually stop firing, entering a depolarization block. This is a robust effect that is not usually explored in experiments or explicitly implemented or tested in models. However, the range of current strength needed for a depolarization block could be easily reached with a random background activity of only a few hundred excitatory synapses. Depolarization block may thus be an important property of neurons that should be better characterized in experiments and explicitly taken into account in models at all implementation scales. Here we analyze the spiking dynamics of CA1 pyramidal neuron models using the same set of ionic currents on both an accurate morphological reconstruction and on its reduction to a single-compartment. The results show the specific ion channel properties and kinetics that are needed to reproduce the experimental findings, and how their interplay can drastically modulate the neuronal dynamics and the input current range leading to a depolarization block. We suggest that this can be one of the rate-limiting mechanisms protecting a CA1 neuron from excessive spiking activity.     

Temperature-sensitive properties of rat suprachiasmatic nucleus neurons

The hypothalamic suprachiasmatic nucleus (SCN) contains a heterogeneous population of neurons, some of which are temperature sensitive in their firing rate activity. Neuronal thermosensitivity may provide cues that synchronize the circadian clock. In addition, through synaptic inhibition on nearby cells, thermosensitive neurons may provide temperature compensation to other SCN neurons, enabling postsynaptic neurons to maintain a constant firing rate despite changes in temperature. To identify mechanisms of neuronal thermosensitivity, whole cell patch recordings monitored resting and transient potentials of SCN neurons in rat hypothalamic tissue slices during changes in temperature. Firing rate temperature sensitivity is not due to thermally dependent changes in the resting membrane potential, action potential threshold, or amplitude of the fast afterhyperpolarizing potential (AHP). The primary mechanism of neuronal thermosensitivity resides in the depolarizing prepotential, which is the slow depolarization that occurs prior to the membrane potential reaching threshold. In thermosensitive neurons, warming increases the prepotential's rate of depolarization, such that threshold is reached sooner. This shortens the interspike interval and increases the firing rate. In some SCN neurons, the slow component of the AHP provides an additional mechanism for thermosensitivity. In these neurons, warming causes the slow AHP to begin at a more depolarized level, and this, in turn, shortens the interspike interval to increase firing rate.     

Inactivation mode of sodium channels defines the different maximal firing rates of conventional versus atypical midbrain dopamine neurons

Two subpopulations of midbrain dopamine (DA) neurons are known to have different dynamic firing ranges in vitro that correspond to distinct projection targets: the originally identified conventional DA neurons project to the dorsal striatum and the lateral shell of the nucleus accumbens, whereas an atypical DA population with higher maximum firing frequencies projects to prefrontal regions and other limbic regions including the medial shell of nucleus accumbens. Using a computational model, we show that previously identified differences in biophysical properties do not fully account for the larger dynamic range of the atypical population and predict that the major difference is that originally identified conventional cells have larger occupancy of voltage-gated sodium channels in a long-term inactivated state that recovers slowly; stronger sodium and potassium conductances during action potential firing are also predicted for the conventional compared to the atypical DA population. These differences in sodium channel gating imply that longer intervals between spikes are required in the conventional population for full recovery from long-term inactivation induced by the preceding spike, hence the lower maximum frequency. These same differences can also change the bifurcation structure to account for distinct modes of entry into depolarization block: abrupt versus gradual. The model predicted that in cells that have entered depolarization block, it is much more likely that an additional depolarization can evoke an action potential in conventional DA population. New experiments comparing lateral to medial shell projecting neurons confirmed this model prediction, with implications for differential synaptic integration in the two populations.     

Reliability of spike and burst firing in thalamocortical relay cells

The reliability and precision of the timing of spikes in a spike train is an important aspect of neuronal coding. We investigated reliability in thalamocortical relay (TCR) cells in the acute slice and also in a Morris-Lecar model with several extensions. A frozen Gaussian noise current, superimposed on a DC current, was injected into the TCR cell soma. The neuron responded with spike trains that showed trial-to-trial variability, due to amongst others slow changes in its internal state and the experimental setup. The DC current allowed to bring the neuron in different states, characterized by a well defined membrane voltage (between ?80 and ?50 mV) and by a specific firing regime that on depolarization gradually shifted from a predominantly bursting regime to a tonic spiking regime. The filtered frozen white noise generated a spike pattern output with a broad spike interval distribution. The coincidence factor and the Hunter and Milton measure were used as reliability measures of the output spike train. In the experimental TCR cell as well as the Morris-Lecar model cell the reliability depends on the shape (steepness) of the current input versus spike frequency output curve. The model also allowed to study the contribution of three relevant ionic membrane currents to reliability: a T-type calcium current, a cation selective h-current and a calcium dependent potassium current in order to allow bursting, investigate the consequences of a more complex current-frequency relation and produce realistic firing rates. The reliability of the output of the TCR cell increases with depolarization. In hyperpolarized states bursts are more reliable than single spikes. The analytically derived relations were capable to predict several of the experimentally recorded spike features.     

Spike latency and jitter of neuronal membrane patches with stochastic Hodgkin-Huxley channels

Ion channel stochasticity can influence the voltage dynamics of neuronal membrane, with stronger effects for smaller patches of membrane because of the correspondingly smaller number of channels. We examine this question with respect to first spike statistics in response to a periodic input of membrane patches including stochastic Hodgkin-Huxley channels, comparing these responses to spontaneous firing. Without noise, firing threshold of the model depends on frequency—a sinusoidal stimulus is subthreshold for low and high frequencies and suprathreshold for intermediate frequencies. When channel noise is added, a stimulus in the lower range of subthreshold frequencies can influence spike output, while high subthreshold frequencies remain subthreshold. Both input frequency and channel noise strength influence spike timing. Specifically, spike latency and jitter have distinct minima as a function of input frequency, showing a resonance like behavior. With either no input, or low frequency subthreshold input, or input in the low or high suprathreshold frequency range, channel noise reduces latency and jitter, with the strongest impact for the lowest input frequencies. In contrast, for an intermediate range of suprathreshold frequencies, where an optimal input gives a minimum latency, the noise effect reverses, and spike latency and jitter increase with channel noise. Thus, a resonant minimum of the spike response as a function of frequency becomes more pronounced with less noise. Spike latency and jitter also depend on the initial phase of the input, resulting in minimal latencies at an optimal phase, and depend on the membrane time constant, with a longer time constant broadening frequency tuning for minimal latency and jitter. Taken together, these results suggest how stochasticity of ion channels may influence spike timing and thus coding for neurons with functionally localized concentrations of channels, such as in “hot spots” of dendrites, spines or axons.     

Modeling the Afferent Dynamics of the Baroreflex Control System

In this study we develop a modeling framework for predicting baroreceptor firing rate as a function of blood pressure. We test models within this framework both quantitatively and qualitatively using data from rats. The models describe three components: arterial wall deformation, stimulation of mechanoreceptors located in the BR nerve-endings, and modulation of the action potential frequency. The three sub-systems are modeled individually following well-established biological principles. The first submodel, predicting arterial wall deformation, uses blood pressure as an input and outputs circumferential strain. The mechanoreceptor stimulation model, uses circumferential strain as an input, predicting receptor deformation as an output. Finally, the neural model takes receptor deformation as an input predicting the BR firing rate as an output. Our results show that nonlinear dependence of firing rate on pressure can be accounted for by taking into account the nonlinear elastic properties of the artery wall. This was observed when testing the models using multiple experiments with a single set of parameters. We find that to model the response to a square pressure stimulus, giving rise to post-excitatory depression, it is necessary to include an integrate-and-fire model, which allows the firing rate to cease when the stimulus falls below a given threshold. We show that our modeling framework in combination with sensitivity analysis and parameter estimation can be used to test and compare models. Finally, we demonstrate that our preferred model can exhibit all known dynamics and that it is advantageous to combine qualitative and quantitative analysis methods.     

Method of evaluating the synaptic input to a single mesencephalic neuron]

The mathematical model of the spike activity of a neuron with synaptic input from many other neurons [1], describes adequately the firing of 5 from 7 neurons in the tegmentum of mesencephalic cat and changes of their activity evoked by glutamate iontophoresis. For these 5 neurons the estimates of the PSPs' average frequency of the threshold depolarization and of the constant decay of the EPSP were received. For different neurons the values of these parameters are 4--100 KHz, 100--800 average unitary EPSPs and 4--30 msec correspondingly. The stationary value of the average membrane potential (SVAMP) in all 5 neurons was removed significantly from the resting potential toward the threshold potential. SWAMP could be changed by the glutamate iontophoresis in such a degree to overlap the threshold potential.     

A Theoretical Analysis of Neuronal Variability

A simple neuronal model is assumed in which, after a refractory period, excitatory and inhibitory exponentially decaying inputs of constant size occur at random intervals and sum until a threshold is reached. The distribution of time intervals between successive neuronal firings (interresponse time histogram), the firing rate as a function of input frequency, the variability in the time course of depolarization from trial to trial, and the strength-duration curve are derived for this model. The predictions are compared with data from the literature and good qualitative agreement is found. All parameters are experimentally measurable and a direct test of the theory is possible with present techniques. The assumptions of the model are relaxed and the effects of such experimentally found phenomena as relative refractory and supernormal periods, adaptation, potentiation, and rhythmic slow potentials are discussed. Implications for gross behavior studies are considered briefly.     

The influence of depolarization block on seizure-like activity in networks of excitatory and inhibitory neurons

The inhibitory restraint necessary to suppress aberrant activity can fail when inhibitory neurons cease to generate action potentials as they enter depolarization block. We investigate possible bifurcation structures that arise at the onset of seizure-like activity resulting from depolarization block in inhibitory neurons. Networks of conductance-based excitatory and inhibitory neurons are simulated to characterize different types of transitions to the seizure state, and a mean field model is developed to verify the generality of the observed phenomena of excitatory-inhibitory dynamics. Specifically, the inhibitory population’s activation function in the Wilson-Cowan model is modified to be non-monotonic to reflect that inhibitory neurons enter depolarization block given strong input. We find that a physiological state and a seizure state can coexist, where the seizure state is characterized by high excitatory and low inhibitory firing rate. Bifurcation analysis of the mean field model reveals that a transition to the seizure state may occur via a saddle-node bifurcation or a homoclinic bifurcation. We explain the hysteresis observed in network simulations using these two bifurcation types. We also demonstrate that extracellular potassium concentration affects the depolarization block threshold; the consequent changes in bifurcation structure enable the network to produce the tonic to clonic phase transition observed in biological epileptic networks.     

Effect of an exponentially decaying threshold on the firing statistics of a stochastic integrate-and-fire neuron

We study a white-noise driven integrate-and-fire (IF) neuron with a time-dependent threshold. We give analytical expressions for mean and variance of the interspike interval assuming that the modification of the threshold value is small. It is shown that the variability of the interval can become both smaller or larger than in the case of constant threshold depending on the decay rate of threshold. We also show that the relative variability is minimal for a certain finite decay rate of the threshold. Furthermore, for slow threshold decay the leaky IF model shows a minimum in the coefficient of variation whenever the firing rate of the neuron matches the decay rate of the threshold. This novel effect can be seen if the firing rate is changed by varying the noise intensity or the mean input current.     

Surface EMG-force modelling for the biceps brachii and its experimental evaluation during isometric isotonic contractions

The aim of this study was to evaluate a surface electromyography (sEMG) signal and force model for the biceps brachii muscle during isotonic isometric contractions for an experimental set-up as well as for a simulation. The proposed model includes a new rate coding scheme and a new analytical formulation of the muscle force generation. The proposed rate coding scheme supposes varying minimum and peak firing frequencies according to motor unit (MU) type (I or II). Practically, the proposed analytical mechanogram allows us to tune the force contribution of each active MU according to its type and instantaneous firing rate. A subsequent sensitivity analysis using a Monte Carlo simulation allows deducing optimised input parameter ranges that guarantee a realistic behaviour of the proposed model according to two existing criteria and an additional one. In fact, this proposed new criterion evaluates the force generation efficiency according to neural intent. Experiments and simulations, at varying force levels and using the optimised parameter ranges, were performed to evaluate the proposed model. As a result, our study showed that the proposed sEMG–force modelling can emulate the biceps brachii behaviour during isotonic isometric contractions.     

Effect of skin surface temperature on localization error and two-point threshold

The average localization error and two-point threshold were determined from 9 subjects at 5 skin temperatures on the volar forearm. The temperatures were physiological zero (31 degrees--32 degrees C), 19 degrees C, 25 degrees C, 38 degrees C, and 44 degrees C. Localization error was not significantly affected by skin surface temperature. This was verified on palm and forehead locations. The two-point threshold was increased by skin surface temperatures above physiological zero and decreased by skin surface temperatures below physiological zero. These results are consistent with a previous report of similar changes in punctate pressure sensitivity with varied skin surface temperatures. Such cross-model effects are considered in terms of cutaneous sensory mechanisms combining broad receptor specificity and subsequent coding of information by patterns of neural impulses.