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.     

周期激励下FHN模型的两种多峰分布放电节律

文章揭示了外界周期脉冲激励下神经元系统产生的随机整数倍和混沌多峰放电节律的关系.随机节律统计直方图呈多峰分布、峰值指数衰减、不可预报且复杂度接近1;混沌节律统计直方图呈不同的多峰分布,峰值非指数衰减、有一定的可预报性且复杂度小于1.混沌节律在激励脉冲周期小于系统内在周期且刺激强度较大时产生,参数范围较小;而随机节律在激励脉冲周期大于系统内在周期且脉冲刺激强度小时,可与随机因素共同作用而产生,产生的参数范围较大.上述结果揭示了两类节律的动力学特性,为区分两类节律提供了实用指标.     

神经放电节律加周期分岔序列中的簇放电到峰放电的转迁

含快慢子系统的神经元数学模型仿真预期,神经放电节律经历加周期分岔序列,可以进一步表现激变,并通过逆倍周期分岔级联进入周期1峰放电。实验调节胞外钙离子浓度,观察到从周期1簇放电开始的带有随机节律的加周期分岔到簇内有多个峰的簇放电,再经激变转迁到峰放电节律的分岔序列,提供了这种分岔序列模式实验证据。实验所见之激变表现为簇放电节律的休止期消失,放电节律变为混沌峰放电和周期峰放电。作者利用随机Chay模型更加逼真地仿真再现了实验所见的分岔序列。该实验结果验证了以前的确定性数学模型的理论预期,并利用随机理论模型仿真了其在现实神经系统的表现;揭示了一类完整的神经放电节律的转换规律。     

交流外电场下映射神经元放电节律的分析

神经元不同的放电节律承载着不同的刺激信息。文章基于神经元映射模型,研究低频交流电场对神经元放电节律的影响。在外部刺激下映射模型表现出丰富的放电模式,包括周期簇放电、周期峰放电、交替放电和混沌放电。神经元对刺激频率和振幅的变化极为敏感,随着频率的增大,放电节律表现出从簇放电到峰放电和混沌放电的反向加周期分岔序列;在周期节律转迁过程中存在一种新的交替节律,其放电序列为两种周期放电模式的交替,峰峰间期序列具有整数倍特征。外电场的频率影响细胞内、外离子振荡周期,导致神经元放电与刺激信号同步,对放电节律的影响更为明显。研究结果揭示了交流外电场对神经元放电节律的作用规律,有助于探寻外电场对生物神经系统兴奋性的影响和神经系统疾病的致病机理。     

Six types of multistability in a neuronal model based on slow calcium current

Background

Multistability of oscillatory and silent regimes is a ubiquitous phenomenon exhibited by excitable systems such as neurons and cardiac cells. Multistability can play functional roles in short-term memory and maintaining posture. It seems to pose an evolutionary advantage for neurons which are part of multifunctional Central Pattern Generators to possess multistability. The mechanisms supporting multistability of bursting regimes are not well understood or classified.

Methodology/Principal Findings

Our study is focused on determining the bio-physical mechanisms underlying different types of co-existence of the oscillatory and silent regimes observed in a neuronal model. We develop a low-dimensional model typifying the dynamics of a single leech heart interneuron. We carry out a bifurcation analysis of the model and show that it possesses six different types of multistability of dynamical regimes. These types are the co-existence of 1) bursting and silence, 2) tonic spiking and silence, 3) tonic spiking and subthreshold oscillations, 4) bursting and subthreshold oscillations, 5) bursting, subthreshold oscillations and silence, and 6) bursting and tonic spiking. These first five types of multistability occur due to the presence of a separating regime that is either a saddle periodic orbit or a saddle equilibrium. We found that the parameter range wherein multistability is observed is limited by the parameter values at which the separating regimes emerge and terminate.

Conclusions

We developed a neuronal model which exhibits a rich variety of different types of multistability. We described a novel mechanism supporting the bistability of bursting and silence. This neuronal model provides a unique opportunity to study the dynamics of networks with neurons possessing different types of multistability.     

Coherence resonance for neuronal bursting with spike undershoot

Although the bursting patterns with spike undershoot are involved with the achievement of physiological or cognitive functions of brain with synaptic noise, noise induced-coherence resonance (CR) from resting state or subthreshold oscillations instead of bursting has been widely identified to play positive roles in information process. Instead, in the present paper, CR characterized by the increase firstly and then decease of peak value of power spectrum of spike trains is evoked from a bursting pattern with spike undershoot, which means that the minimal membrane potential within burst is lower than that of the subthreshold oscillations between bursts, while CR cannot be evoked from the bursting pattern without spike undershoot. With bifurcations and fast-slow variable dissection method, the bursting patterns with and without spike undershoot are classified into “Sub-Hopf/Fold” bursting and “Fold/Homoclinic” bursting, respectively. For the bursting with spike undershoot, the trajectory of the subthreshold oscillations is very close to that of the spikes within burst. Therefore, noise can induce more spikes from the subthreshold oscillations and modulate the bursting regularity, which leads to the appearance of CR. For the bursting pattern without spike undershoot, the trajectory of the quiescent state is not close to that of the spikes within burst, and noise cannot induce spikes from the quiescent state between bursts, which is cause for non-CR. The result provides a novel case of CR phenomenon and extends the scopes of CR concept, presents that noise can enhance rather than suppress information of the bursting patterns with spike undershoot, which are helpful for understanding the dynamics and the potential physiological or cognitive functions of the nerve fiber or brain neurons with such bursting patterns.     

神经起步点自发放电节律及节律转化的分岔规律

在神经起步点的实验中观察到了复杂多样的神经放电([Ca^2 ]o)节律模式,如周期簇放电、周期峰放电、混沌簇放电、混沌峰放电以及随机放电节律等。随着细胞外钙离子浓度的降低,神经放电节律从周期l簇放电,经过复杂的分岔过程(包括经倍周期分岔到混沌簇放电、混沌簇放电经激变到混沌峰放电、以及混沌峰放电经逆倍周期分岔到周期峰放电)转化为周期l峰放电。在神经放电理论模型——Chay模型中,调节与实验相关的参数(Ca^2 平衡电位),可以获得与实验相似的神经放电节律和节律转换规律。这表明复杂的神经放电节律之间存在着一定的分岔规律,它们是理解神经元信息编码的基础。     

时滞援助的两抑制性突触耦合的Chay 神经元的同步

一对抑制性突触耦合的混沌Chay神经元的同步模式被研究。结果表明当耦合强度超过临界值时,两抑制耦合的混沌Chay神经元能达到反相的同步。与此同时,两混沌的神经元变为周期而不是原来的混沌运动。然而,如果考虑耦合神经元信息的传导时滞,在有效的时滞下,两个耦合神经元的在相簇同步能增加。在相簇同步窗口的大小随着耦合强度的增加而增加。此结果对于我们理解神经元集群的运动是一个指导。     

Noise-induced impulse pattern modifications at different dynamical period-one situations in a computer model of temperature encoding

We used a minimal Hodgkin-Huxley type model of cold receptor discharges to examine how noise interferes with the non-linear dynamics of the ionic mechanisms of neuronal stimulus encoding. The model is based on the assumption that spike-generation depends on subthreshold oscillations. With physiologically plausible temperature scaling, it passes through different impulse patterns which, with addition of noise, are in excellent agreement with real experimental data. The interval distributions of purely deterministic simulations, however, exhibit considerable differences compared to the noisy simulations especially at the bifurcations of deterministically period-one discharges. We, therefore, analyzed the effects of noise in different situations of deterministically regular period-one discharges: (1) at high-temperatures near the transition to subthreshold oscillations and to burst discharges, and (2) at low-temperatures close to and more far away from the bifurcations to chaotic dynamics. The data suggest that addition of noise can considerably extend the dynamical behavior of the system with coexistence of different dynamical situations at deterministically fixed parameter constellations. Apart from well-described coexistence of spike-generating and subthreshold oscillations also mixtures of tonic and bursting patterns can be seen and even transitions to unstable period-one orbits seem to appear. The data indicate that cooperative effects between low- and high-dimensional dynamics have to be considered as qualitatively important factors in neuronal encoding.     

Complex bursting in pancreatic islets: a potential glycolytic mechanism

The electrical activity of insulin-secreting pancreatic islets of Langerhans is characterized by bursts of action potentials. Most often this bursting is periodic, but in some cases it is modulated by an underlying slower rhythm. We suggest that the modulatory rhythm for this complex bursting pattern is due to oscillations in glycolysis, while the bursting itself is generated by some other slow process. To demonstrate this hypothesis, we couple a minimal model of glycolytic oscillations to a minimal model for activity-dependent bursting in islets. We show that the combined model can reproduce several complex bursting patterns from mouse islets published in the literature, and we illustrate how these complex oscillations are produced through the use of a fast/slow analysis.     

Bursting, beating, and chaos in an excitable membrane model.

We have studied periodic as well as aperiodic behavior in the self-sustained oscillations exhibited by the Hodgkin-Huxley type model of Chay, T. R., and J. Keizer (Biophys. J., 1983, 42:181-190) for the pancreatic beta-cell. Numerical solutions reveal a variety of patterns as the glucose-dependent parameter kCa is varied. These include regimes of periodic beating (continuous spiking) and bursting modes and, in the transition between these modes, aperiodic responses. Such aperiodic behavior for a nonrandom system has been called deterministic chaos and is characterized by distinguishing features found in previous studies of chaos in nonbiophysical systems and here identified for an (endogenously active) excitable membrane model. To parallel the successful analysis of chaos in other physical/chemical contexts we introduce a simplified, but quantitative, one-variable, discrete-time representation of the dynamics. It describes the evolution of intracellular calcium (which activates a potassium conductance) from one spike upstroke to the next and exhibits the various modes of behavior.     

Analysis of Chaotic Resonance in Izhikevich Neuron Model

In stochastic resonance (SR), the presence of noise helps a nonlinear system amplify a weak (sub-threshold) signal. Chaotic resonance (CR) is a phenomenon similar to SR but without stochastic noise, which has been observed in neural systems. However, no study to date has investigated and compared the characteristics and performance of the signal responses of a spiking neural system in some chaotic states in CR. In this paper, we focus on the Izhikevich neuron model, which can reproduce major spike patterns that have been experimentally observed. We examine and classify the chaotic characteristics of this model by using Lyapunov exponents with a saltation matrix and Poincaré section methods in order to address the measurement challenge posed by the state-dependent jump in the resetting process. We found the existence of two distinctive states, a chaotic state involving primarily turbulent movement and an intermittent chaotic state. In order to assess the signal responses of CR in these classified states, we introduced an extended Izhikevich neuron model by considering weak periodic signals, and defined the cycle histogram of neuron spikes as well as the corresponding mutual correlation and information. Through computer simulations, we confirmed that both chaotic states in CR can sensitively respond to weak signals. Moreover, we found that the intermittent chaotic state exhibited a prompter response than the chaotic state with primarily turbulent movement.     

Recurrent inhibitory dynamics: the role of state-dependent distributions of conduction delay times

We have formulated and analysed a dynamic model for recurrent inhibition that takes into account the state dependence of the delayed feedback signal (due to the variation in threshold of fibres with their size) and the distribution of these delays (due to the distribution of fibre diameters in the feedback pathway). Using a combination of analytic and numerical tools, we have analysed the behaviour of this model. Depending on the parameter values chosen, as well as the initial preparation of the system, there may be a spectrum of post-synaptic firing dynamics ranging from stable constant values through periodic bursting (limit cycle) behaviour and chaotic firing as well as bistable behaviours. Using detailed parameter estimation for a physiologically motivated example (the CA3-basket cell-mossy fibre system in the hippocampus), we present some of these numerical behaviours. The numerical results corroborate the results of the analytic characterization of the solutions. Namely, for some parameter values the model has a single stable steady state while for the others there is a bistability in which the eventual behaviour depends on the magnitude of stimulation (the initial function).     

Dynamics and bifurcations of the adaptive exponential integrate-and-fire model

Recently, several two-dimensional spiking neuron models have been introduced, with the aim of reproducing the diversity of electrophysiological features displayed by real neurons while keeping a simple model, for simulation and analysis purposes. Among these models, the adaptive integrate-and-fire model is physiologically relevant in that its parameters can be easily related to physiological quantities. The interaction of the differential equations with the reset results in a rich and complex dynamical structure. We relate the subthreshold features of the model to the dynamical properties of the differential system and the spike patterns to the properties of a Poincaré map defined by the sequence of spikes. We find a complex bifurcation structure which has a direct interpretation in terms of spike trains. For some parameter values, spike patterns are chaotic.     

Coexistence of Tonic Spiking Oscillations in a Leech Neuron Model

The leech neuron model studied here has a remarkable dynamical plasticity. It exhibits a wide range of activities including various types of tonic spiking and bursting. In this study we apply methods of the qualitative theory of dynamical systems and the bifurcation theory to analyze the dynamics of the leech neuron model with emphasis on tonic spiking regimes. We show that the model can demonstrate bi-stability, such that two modes of tonic spiking coexist. Under a certain parameter regime, both tonic spiking modes are represented by the periodic attractors. As a bifurcation parameter is varied, one of the attractors becomes chaotic through a cascade of period-doubling bifurcations, while the other remains periodic. Thus, the system can demonstrate co-existence of a periodic tonic spiking with either periodic or chaotic tonic spiking. Pontryagins averaging technique is used to locate the periodic orbits in the phase space.     

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.     

Pollen coupling of forest trees: forming synchronized and periodic reproduction out of chaos

Many of the tree species in mature forests show masting; their reproductive activity has a large variance between years and is often synchronized between different individuals. In this paper, we analyse a globally coupled map model in which trees accumulate photosynthate every year, produce flowers when the energy reserve level exceeds a threshold, and set seeds and fruits at a rate limited by pollen availability. Without pollen limitation, the trees in the forest show independent chaotic fluctuation. Coupling of trees via pollen exchange results in reproduction being synchronized partially or completely over the forest. The whole forest shows diverse dynamical behaviors determined by the values of two essential parameters; the depletion coefficient k and the coupling strength beta. We find perfectly synchronized periodic reproduction, synchronized reproduction with a chaotic time series, clustering phenomena, and chaotic reproduction of trees without synchronization over individuals. There are many parameter windows in which synchronized reproduction of trees shows a stable periodic fluctuation. For perfectly synchronized forests, we can calculate all the Lyapunov exponents analytically. They show that synchronized reproduction of all the trees in the forest can only occur when trees flower at low (but positive) levels in a significant fraction of years, resulting in small fruit sets due to outcrossed pollen limitation. This is consistent with the observation that the distinction between mast years and non-mast years is often not clear cut.     

Pattern recall in networks of chaotic neurons

This research investigates the potential utility of chaotic dynamics in neural information processing. A novel chaotic spiking neural network model is presented which is composed of non-linear dynamic state (NDS) neurons. The activity of each NDS neuron is driven by a set of non-linear equations coupled with a threshold based spike output mechanism. If time-delayed self-connections are enabled then the network stabilises to a periodic pattern of activation. Previous publications of this work have demonstrated that the chaotic dynamics which drive the network activity ensure that an extremely large number of such periodic patterns can be generated by this network. This paper presents a major extension to this model which enables the network to recall a pattern of activity from a selection of previously stabilised patterns.     

Switching from simple to complex oscillations in calcium signaling

We present a new model for calcium oscillations based on experiments in hepatocytes. The model considers feedback inhibition on the initial agonist receptor complex by calcium and activated phospholipase C, as well as receptor type-dependent self-enhanced behavior of the activated G(alpha) subunit. It is able to show simple periodic oscillations and periodic bursting, and it is the first model to display chaotic bursting in response to agonist stimulations. Moreover, our model offers a possible explanation for the differences in dynamic behavior observed in response to different agonists in hepatocytes.     

神经放电加周期分岔中由随机自共振引起一类新节律

当改变实验性神经起步点细胞外[Ca^2 ]时,放电节律表现出从周期1节律转换为周期4节律的加周期分岔序列。其中,周期n节律转换为周期n 1节律的过程中(n=1,2,3)存在一种新的具有交替特征的节律,该新节律为周期n簇与周期n 1簇放电的交替,并且周期n 1簇的时间间隔序列呈现出整数倍特征。确定性神经放电理论模型(chay模型)只能模拟周期n节律直接到周期n 1节律的加周期分岔序列;而随机chay模型可以模拟实验中的加周期分岔过程和新节律。进一步,新节律被确认是经随机自共振机制产生的。这不仅解释了实验现象,也将随机自共振的产生区间从以前认识到的Hopf分岔点附近扩大到加周期分岔点附近,同时扩大了噪声在神经放电和神经编码中起重要作用的参数区间。