State-dependent effects of Na channel noise on neuronal burst generation

Dynamical response of Autaptic Izhikevich Neuron disturbed by Gaussian white noise

Using the improved memristive Izhikevich neuron model, the effects of autaptic connection as well as electromagnetic induction are studied on the dynamical behavior of neuronal spiking. Using bifurcation analysis for membrane potentials, the effects of autaptic and electromagnetic parameters on the mode transition in electrical activities of the neuron model are investigated. Furthermore, white Gaussian noise is considered in the neuron model, to evaluate the effect of electromagnetic disturbance on the firing pattern of the neuron using the coefficient of variation. The bifurcation diagram versus autaptic conductance and time delay has been extensively studied. The results show that the effects of autaptic connection as well as electromagnetic induction on the spiking behavior of neurons can be well demonstrated by using the Izhikevich model. The electrical activities of the Izhikevich neuron model become more complex when the effects of autaptic connection and electromagnetic induction are considered in the neuron model. Using the Izhikevich neuron model, the high variety of spiking/bursting patterns is represented in the bifurcation diagram of inter-spike interval versus autaptic or electromagnetic parameters. Noise can have distinct effects on the spiking activity of the neuron, for the subthreshold input current, increasing the intensity of the electromagnetic noise increases the regularity of the neuron spiking, but for the suprathreshold input current, the regularity of spiking decreases with noise.

Oscillations and variability in neuronal systems: interplay of autonomous transient dynamics and fast deterministic fluctuations

Neuronal systems are subject to rapid fluctuations both intrinsically and externally. These fluctuations can be disruptive or constructive. We investigate the dynamic mechanisms underlying the interactions between rapidly fluctuating signals and the intrinsic properties of the target cells to produce variable and/or coherent responses. We use linearized and non-linear conductance-based models and piecewise constant (PWC) inputs with short duration pieces. The amplitude distributions of the constant pieces consist of arbitrary permutations of a baseline PWC function. In each trial within a given protocol we use one of these permutations and each protocol consists of a subset of all possible permutations, which is the only source of uncertainty in the protocol. We show that sustained oscillatory behavior can be generated in response to various forms of PWC inputs independently of whether the stable equilibria of the corresponding unperturbed systems are foci or nodes. The oscillatory voltage responses are amplified by the model nonlinearities and attenuated for conductance-based PWC inputs as compared to current-based PWC inputs, consistent with previous theoretical and experimental work. In addition, the voltage responses to PWC inputs exhibited variability across trials, which is reminiscent of the variability generated by stochastic noise (e.g., Gaussian white noise). Our analysis demonstrates that both oscillations and variability are the result of the interaction between the PWC input and the target cell's autonomous transient dynamics with little to no contribution from the dynamics in vicinities of the steady-state, and do not require input stochasticity.

Ion channel noise shapes the electrical activity of endocrine cells

Endocrine cells in the pituitary gland typically display either spiking or bursting electrical activity, which is related to the level of hormone secretion. Recent work, which combines mathematical modelling with dynamic clamp experiments, suggests the difference is due to the presence or absence of a few large-conductance potassium channels. Since endocrine cells only contain a handful of these channels, it is likely that stochastic effects play an important role in the pattern of electrical activity. Here, for the first time, we explicitly determine the effect of such noise by studying a mathematical model that includes the realistic noisy opening and closing of ion channels. This allows us to investigate how noise affects the electrical activity, examine the origin of spiking and bursting, and determine which channel types are responsible for the greatest noise. Further, for the first time, we address the role of cell size in endocrine cell electrical activity, finding that larger cells typically display more bursting, while the smallest cells almost always only exhibit spiking behaviour.

Inverse stochastic resonance in networks of spiking neurons

Inverse Stochastic Resonance (ISR) is a phenomenon in which the average spiking rate of a neuron exhibits a minimum with respect to noise. ISR has been studied in individual neurons, but here, we investigate ISR in scale-free networks, where the average spiking rate is calculated over the neuronal population. We use Hodgkin-Huxley model neurons with channel noise (i.e., stochastic gating variable dynamics), and the network connectivity is implemented via electrical or chemical connections (i.e., gap junctions or excitatory/inhibitory synapses). We find that the emergence of ISR depends on the interplay between each neuron's intrinsic dynamical structure, channel noise, and network inputs, where the latter in turn depend on network structure parameters. We observe that with weak gap junction or excitatory synaptic coupling, network heterogeneity and sparseness tend to favor the emergence of ISR. With inhibitory coupling, ISR is quite robust. We also identify dynamical mechanisms that underlie various features of this ISR behavior. Our results suggest possible ways of experimentally observing ISR in actual neuronal systems.

Noise, transient dynamics, and the generation of realistic interspike interval variation in square-wave burster neurons

First return maps of interspike intervals for biological neurons that generate repetitive bursts of impulses can display stereotyped structures (neuronal signatures). Such structures have been linked to the possibility of multicoding and multifunctionality in neural networks that produce and control rhythmical motor patterns. In some cases, isolating the neurons from their synaptic network reveals irregular, complex signatures that have been regarded as evidence of intrinsic, chaotic behavior. We show that incorporation of dynamical noise into minimal neuron models of square-wave bursting (either conductance-based or abstract) produces signatures akin to those observed in biological examples, without the need for fine tuning of parameters or ad hoc constructions for inducing chaotic activity. The form of the stochastic term is not strongly constrained and can approximate several possible sources of noise, e.g., random channel gating or synaptic bombardment. The cornerstone of this signature generation mechanism is the rich, transient, but deterministic dynamics inherent in the square-wave (saddle-node and homoclinic) mode of neuronal bursting. We show that noise causes the dynamics to populate a complex transient scaffolding or skeleton in state space, even for models that (without added noise) generate only periodic activity (whether in bursting or tonic spiking mode).

An Investigation of the Stochastic Hodgkin-Huxley Models Under Noisy Rate Functions

The effects of ion channel fluctuations on the transmembrane voltage activity are potentially profound in small-size excitable membrane patches. Different groups have extended Hodgkin-Huxley equations into stochastic differential equations to capture the effects of ion channel noise analytically (Fox & Lu, 1994 ; Linaro, Storace, & Giugliano, 2011 ; Güler, 2013 ). Studies have shown that the accuracy of spiking statistics by Fox and Lu's model does not match well with the corresponding statistics from the exact microscopic simulations. The models of both Linaro et al. and Güler, however, were found to produce highly accurate statistics. Here we extend the examination of these models to the case in which the rate functions for the opening and closing of gates are under the influence of noise. For that purpose, the usual rate functions are accompanied additively by Ornstein-Uhlenbeck–type stochastic angular variables. Moreover, we argue that the existence of such noise in the rate functions is a plausible physiological phenomenon for finite-size membranes. It is observed that the presence of noise in the rates is not effective on the degree of inaccuracies within the Fox and Lu model. Güler model's accuracy is found to remain high as in the case of noise free rates. But the performance of Linaro et al.’s model is seen to degrade seriously with the increasing strength of the introduced rate function noise. We attribute this failure of Linaro et al.’s model to the use of the covariance function of open channels at the steady state, in its derivation.

Stochastic Hodgkin-Huxley Equations with Colored Noise Terms in the Conductances

The excitability of cells is facilitated by voltage-gated ion channels. These channels accommodate a multiple number of gates individually. The possible impact of that gate multiplicity on the cell's function, specifically when the membrane area is of limited size, was investigated in the author's prior work (Güler, 2011 ). There, it was found that a nontrivially persistent correlation takes place between the transmembrane voltage fluctuations (also between the fluctuations in the gating variables) and the component of open channel fluctuations attributed to the gate multiplicity. This nontrivial phenomenon was found to be playing a major augmentative role for the elevation of excitability and spontaneous firing in small cells. In addition, the same phenomenon was found to be enhancing spike coherence significantly. Here we extend Fox and Lu's ( 1994 ) stochastic Hodgkin-Huxley equations by incorporating colored noise terms into the conductances there to obtain a formalism capable of capturing the addressed cross-correlations. Statistics of spike generation, spike coherence, firing efficiency, latency, and jitter from the articulated set of equations are found to be highly accurate in comparison with the corresponding statistics from the exact microscopic Markov simulations. This way, it is demonstrated vividly that our formulation overcomes the inherent inadequacy of the Fox and Lu equations. Finally, a recently proposed diffusion approximation method (Linaro, Storace, & Giugliano, 2011 ) is taken into consideration, and a discussion on its character is pursued.

Modifying spiking precision in conductance-based neuronal models

The temporal precision of a neuron's spiking can be characterized by calculating its "jitter," defined as the standard deviation of the timing of individual spikes in response to repeated presentations of a stimulus. Sub-millisecond jitters have been measured for neurons in a variety of experimental systems and appear to be functionally important in some instances. We have investigated how modifying a neuron's maximal conductances affects jitter using the leaky integrate-and-fire (LIF) model and an eight-conductance Hodgkin-Huxley type (HH8) model. We observed that jitter can be largely understood in the LIF model in terms of the neuron's filtering properties. In the HH8 model we found the role of individual conductances in determining jitter to be complicated and dependent on the model's spiking properties. Distinct behaviors were observed for populations with slow (<11.5 Hz) and fast (>11.5 Hz) spike rates and appear to be related to differences in a particular channel's activity at times just before spiking occurs.

Noise destroys feedback enhanced figure-ground segmentation but not feedforward figure-ground segmentation

Figure-ground (FG) segmentation is the separation of visual information into background and foreground objects. In the visual cortex, FG responses are observed in the late stimulus response period, when neurons fire in tonic mode, and are accompanied by a switch in cortical state. When such a switch does not occur, FG segmentation fails. Currently, it is not known what happens in the brain on such occasions. A biologically plausible feedforward spiking neuron model was previously devised that performed FG segmentation successfully. After incorporating feedback the FG signal was enhanced, which was accompanied by a change in spiking regime. In a feedforward model neurons respond in a bursting mode whereas in the feedback model neurons fired in tonic mode. It is known that bursts can overcome noise, while tonic firing appears to be much more sensitive to noise. In the present study, we try to elucidate how the presence of noise can impair FG segmentation, and to what extent the feedforward and feedback pathways can overcome noise. We show that noise specifically destroys the feedback enhanced FG segmentation and leaves the feedforward FG segmentation largely intact. Our results predict that noise produces failure in FG perception.

Simple, fast and accurate implementation of the diffusion approximation algorithm for stochastic ion channels with multiple states

BACKGROUND

The phenomena that emerge from the interaction of the stochastic opening and closing of ion channels (channel noise) with the non-linear neural dynamics are essential to our understanding of the operation of the nervous system. The effects that channel noise can have on neural dynamics are generally studied using numerical simulations of stochastic models. Algorithms based on discrete Markov Chains (MC) seem to be the most reliable and trustworthy, but even optimized algorithms come with a non-negligible computational cost. Diffusion Approximation (DA) methods use Stochastic Differential Equations (SDE) to approximate the behavior of a number of MCs, considerably speeding up simulation times. However, model comparisons have suggested that DA methods did not lead to the same results as in MC modeling in terms of channel noise statistics and effects on excitability. Recently, it was shown that the difference arose because MCs were modeled with coupled gating particles, while the DA was modeled using uncoupled gating particles. Implementations of DA with coupled particles, in the context of a specific kinetic scheme, yielded similar results to MC. However, it remained unclear how to generalize these implementations to different kinetic schemes, or whether they were faster than MC algorithms. Additionally, a steady state approximation was used for the stochastic terms, which, as we show here, can introduce significant inaccuracies.

MAIN CONTRIBUTIONS

We derived the SDE explicitly for any given ion channel kinetic scheme. The resulting generic equations were surprisingly simple and interpretable--allowing an easy, transparent and efficient DA implementation, avoiding unnecessary approximations. The algorithm was tested in a voltage clamp simulation and in two different current clamp simulations, yielding the same results as MC modeling. Also, the simulation efficiency of this DA method demonstrated considerable superiority over MC methods, except when short time steps or low channel numbers were used.

Neuromorphic control of stepping pattern generation: a dynamic model with analog circuit implementation

Animals such as stick insects can adaptively walk on complex terrains by dynamically adjusting their stepping motion patterns. Inspired by the coupled Matsuoka and resonate-and-fire neuron models, we present a nonlinear oscillation model as the neuromorphic central pattern generator (CPG) for rhythmic stepping pattern generation. This dynamic model can also be used to actuate the motoneurons on a leg joint with adjustable driving frequencies and duty cycles by changing a few of the model parameters while operating such that different stepping patterns can be generated. A novel mixed-signal integrated circuit design of this dynamic model is subsequently implemented, which, although simplified, shares the equivalent output performance in terms of the adjustable frequency and duty cycle. Three identical CPG models being used to drive three joints can make an arthropod leg of three degrees of freedom. With appropriate initial circuit parameter settings, and thus suitable phase lags among joints, the leg is expected to walk on a complex terrain with adaptive steps. The adaptation is associated with the circuit parameters mediated both by the higher level nervous system and the lower level sensory signals. The model is realized using a 0.3- complementary metal-oxide-semiconductor process and the results are reported.

Temperature-dependent stochastic dynamics of the Huber-Braun neuron model

The response of a four-dimensional mammalian cold receptor model to different implementations of noise is studied across a wide temperature range. It is observed that for noisy activation kinetics, the parameter range decomposes into two regions in which the system reacts qualitatively completely different to small perturbations through noise, and these regions are separated by a homoclinic bifurcation. Noise implemented as an additional current yields a substantially different system response at low temperature values, while the response at high temperatures is comparable to activation-kinetic noise. We elucidate how this phenomenon can be understood in terms of state space dynamics and gives quantitative results on the statistics of interspike interval distributions across the relevant parameter range.

The what and where of adding channel noise to the Hodgkin-Huxley equations

Conductance-based equations for electrically active cells form one of the most widely studied mathematical frameworks in computational biology. This framework, as expressed through a set of differential equations by Hodgkin and Huxley, synthesizes the impact of ionic currents on a cell's voltage--and the highly nonlinear impact of that voltage back on the currents themselves--into the rapid push and pull of the action potential. Later studies confirmed that these cellular dynamics are orchestrated by individual ion channels, whose conformational changes regulate the conductance of each ionic current. Thus, kinetic equations familiar from physical chemistry are the natural setting for describing conductances; for small-to-moderate numbers of channels, these will predict fluctuations in conductances and stochasticity in the resulting action potentials. At first glance, the kinetic equations provide a far more complex (and higher-dimensional) description than the original Hodgkin-Huxley equations or their counterparts. This has prompted more than a decade of efforts to capture channel fluctuations with noise terms added to the equations of Hodgkin-Huxley type. Many of these approaches, while intuitively appealing, produce quantitative errors when compared to kinetic equations; others, as only very recently demonstrated, are both accurate and relatively simple. We review what works, what doesn't, and why, seeking to build a bridge to well-established results for the deterministic equations of Hodgkin-Huxley type as well as to more modern models of ion channel dynamics. As such, we hope that this review will speed emerging studies of how channel noise modulates electrophysiological dynamics and function. We supply user-friendly MATLAB simulation code of these stochastic versions of the Hodgkin-Huxley equations on the ModelDB website (accession number 138950) and http://www.amath.washington.edu/~etsb/tutorials.html.

Identification and continuity of the distributions of burst-length and interspike intervals in the stochastic Morris-Lecar neuron

Using the Morris-Lecar model neuron with a type II parameter set and K(+)-channel noise, we investigate the interspike interval distribution as increasing levels of applied current drive the model through a subcritical Hopf bifurcation. Our goal is to provide a quantitative description of the distributions associated with spiking as a function of applied current. The model generates bursty spiking behavior with sequences of random numbers of spikes (bursts) separated by interburst intervals of random length. This kind of spiking behavior is found in many places in the nervous system, most notably, perhaps, in stuttering inhibitory interneurons in cortex. Here we show several practical and inviting aspects of this model, combining analysis of the stochastic dynamics of the model with estimation based on simulations. We show that the parameter of the exponential tail of the interspike interval distribution is in fact continuous over the entire range of plausible applied current, regardless of the bifurcations in the phase portrait of the model. Further, we show that the spike sequence length, apparently studied for the first time here, has a geometric distribution whose associated parameter is continuous as a function of applied current over the entire input range. Hence, this model is applicable over a much wider range of applied current than has been thought.

Comparison of Langevin and Markov channel noise models for neuronal signal generation

The stochastic opening and closing of voltage-gated ion channels produce noise in neurons. The effect of this noise on the neuronal performance has been modeled using either an approximate or Langevin model based on stochastic differential equations or an exact model based on a Markov process model of channel gating. Yet whether the Langevin model accurately reproduces the channel noise produced by the Markov model remains unclear. Here we present a comparison between Langevin and Markov models of channel noise in neurons using single compartment Hodgkin-Huxley models containing either Na+ and K+, or only K+ voltage-gated ion channels. The performance of the Langevin and Markov models was quantified over a range of stimulus statistics, membrane areas, and channel numbers. We find that in comparison to the Markov model, the Langevin model underestimates the noise contributed by voltage-gated ion channels, overestimating information rates for both spiking and nonspiking membranes. Even with increasing numbers of channels, the difference between the two models persists. This suggests that the Langevin model may not be suitable for accurately simulating channel noise in neurons, even in simulations with large numbers of ion channels.

Noise-induced transitions in slow wave neuronal dynamics

Many neuronal systems exhibit slow random alternations and sudden switches in activity states. Models with noisy relaxation dynamics (oscillatory, excitable or bistable) account for these temporal, slow wave, patterns and the fluctuations within states. The noise-induced transitions in a relaxation dynamics are analogous to escape by a particle in a slowly changing double-well potential. In this formalism, we obtain semi-analytically the first and second order statistical properties: the distributions of the slow process at the transitions and the temporal correlations of successive switching events. We find that the temporal correlations can be used to help distinguish among biophysical mechanisms for the slow negative feedback, such as divisive or subtractive. We develop our results in the context of models for cellular pacemaker neurons; they also apply to mean-field models for spontaneously active networks with slow wave dynamics.

Variability of bursting patterns in a neuron model in the presence of noise

Spiking and bursting patterns of neurons are characterized by a high degree of variability. A single neuron can demonstrate endogenously various bursting patterns, changing in response to external disturbances due to synapses, or to intrinsic factors such as channel noise. We argue that in a model of the leech heart interneuron existing variations of bursting patterns are significantly enhanced by a small noise. In the absence of noise this model shows periodic bursting with fixed numbers of interspikes for most parameter values. As the parameter of activation kinetics of a slow potassium current is shifted to more hyperpolarized values of the membrane potential, the model undergoes a sequence of incremental spike adding transitions accumulating towards a periodic tonic spiking activity. Within a narrow parameter window around every spike adding transition, spike alteration of bursting is deterministically chaotic due to homoclinic bifurcations of a saddle periodic orbit. We have found that near these transitions the interneuron model becomes extremely sensitive to small random perturbations that cause a wide expansion and overlapping of the chaotic windows. The chaotic behavior is characterized by positive values of the largest Lyapunov exponent, and of the Shannon entropy of probability distribution of spike numbers per burst. The windows of chaotic dynamics resemble the Arnold tongues being plotted in the parameter plane, where the noise intensity serves as a second control parameter. We determine the critical noise intensities above which the interneuron model generates only irregular bursting within the overlapped windows.

Evaluation of stochastic differential equation approximation of ion channel gating models

Fox and Lu derived an algorithm based on stochastic differential equations for approximating the kinetics of ion channel gating that is simpler and faster than "exact" algorithms for simulating Markov process models of channel gating. However, the approximation may not be sufficiently accurate to predict statistics of action potential generation in some cases. The objective of this study was to develop a framework for analyzing the inaccuracies and determining their origin. Simulations of a patch of membrane with voltage-gated sodium and potassium channels were performed using an exact algorithm for the kinetics of channel gating and the approximate algorithm of Fox & Lu. The Fox & Lu algorithm assumes that channel gating particle dynamics have a stochastic term that is uncorrelated, zero-mean Gaussian noise, whereas the results of this study demonstrate that in many cases the stochastic term in the Fox & Lu algorithm should be correlated and non-Gaussian noise with a non-zero mean. The results indicate that: (i) the source of the inaccuracy is that the Fox & Lu algorithm does not adequately describe the combined behavior of the multiple activation particles in each sodium and potassium channel, and (ii) the accuracy does not improve with increasing numbers of channels.

Detailed numerical investigation of the dissipative stochastic mechanics based neuron model

Recently, a physical approach for the description of neuronal dynamics under the influence of ion channel noise was proposed in the realm of dissipative stochastic mechanics (Güler, Phys Rev E 76:041918, 2007). Led by the presence of a multiple number of gates in an ion channel, the approach establishes a viewpoint that ion channels are exposed to two kinds of noise: the intrinsic noise, associated with the stochasticity in the movement of gating particles between the inner and the outer faces of the membrane, and the topological noise, associated with the uncertainty in accessing the permissible topological states of open gates. Renormalizations of the membrane capacitance and of a membrane voltage dependent potential function were found to arise from the mutual interaction of the two noisy systems. The formalism therein was scrutinized using a special membrane with some tailored properties giving the Rose-Hindmarsh dynamics in the deterministic limit. In this paper, the resultant computational neuron model of the above approach is investigated in detail numerically for its dynamics using time-independent input currents. The following are the major findings obtained. The intrinsic noise gives rise to two significant coexisting effects: it initiates spiking activity even in some range of input currents for which the corresponding deterministic model is quiet and causes bursting in some other range of input currents for which the deterministic model fires tonically. The renormalization corrections are found to augment the above behavioral transitions from quiescence to spiking and from tonic firing to bursting, and, therefore, the bursting activity is found to take place in a wider range of input currents for larger values of the correction coefficients. Some findings concerning the diffusive behavior in the voltage space are also reported.

Dissipative stochastic mechanics for capturing neuronal dynamics under the influence of ion channel noise: formalism using a special membrane

Based on the idea conveyed in the author's prior study [Fluct. Noise Lett. 6, L147 (2006)], a physical approach for the description of neuronal dynamics under the influence of ion channel noise is developed in the realm of Nelson's stochastic mechanics when open to dissipative environments. The formalism therein is scrutinized using a special membrane with some tailored properties giving the Rose-Hindmarsh dynamics in the deterministic limit. Led by the presence of multiple number of gates in an ion channel, a dual viewpoint of channel noise is established. Then, stochastic mechanics is adopted to model those channel fluctuations emerging from the uncertainty in accessing the permissible topological states of open gates. A mutual interaction between the above fluctuations and the noise, emerging from the stochasticity in the movement of gating particles between the inner and the outer faces of the membrane, is portrayed within a system plus reservoir strategy. Induced by the interaction, renormalizations of the membrane capacitance and of a membrane voltage dependent potential are found to arise. Consequently, the equations of motion, for the expectation values of the variables and the pair correlation functions, are obtained in the collective membrane voltage space.

Interspike interval statistics in the stochastic Hodgkin-Huxley model: coexistence of gamma frequency bursts and highly irregular firing

When the classical Hodgkin-Huxley equations are simulated with Na- and K-channel noise and constant applied current, the distribution of interspike intervals is bimodal: one part is an exponential tail, as often assumed, while the other is a narrow gaussian peak centered at a short interspike interval value. The gaussian arises from bursts of spikes in the gamma-frequency range, the tail from the interburst intervals, giving overall an extraordinarily high coefficient of variation--up to 2.5 for 180,000 Na channels when I approximately 7 microA/cm(2). Since neurons with a bimodal ISI distribution are common, it may be a useful model for any neuron with class 2 firing. The underlying mechanism is due to a subcritical Hopf bifurcation, together with a switching region in phase-space where a fixed point is very close to a system limit cycle. This mechanism may be present in many different classes of neurons and may contribute to widely observed highly irregular neural spiking.

The rhythmic consequences of ion channel stochasticity

Ion channels are membrane-spanning proteins with central pores through which ions cross neuronal membranes. The pores through each ion channel flicker between open and closed states, starting and stopping the flow of ions and the electrical current they carry. Hence the current flickers on and off, varying widely on very short time scales. Recent evidence suggests that this noisy current is a source of rhythmic behaviors in neurons. In this review, we begin by providing an illustrative model that links the stochastic flicker of ion channels to neuronal rhythms. The author explores recent experimental work that shows channel flicker is necessary for at least one rhythm that characterizes a class of cortical neurons in vitro. Finally, the author highlights a number of novel studies that link ion channel stochasticity to neuronal rhythmic behaviors in other interesting ways.

Perinatal maturation of the respiratory rhythm generator in mammals: from experimental results to computational simulation

The survival of neonatal mammals requires a correct function of the respiratory rhythm generator (RRG), and therefore, the processes that control its prenatal maturation are of vital importance. In humans, lambs and rodents, foetal breathing movements (FBMs) occur early during gestation, are episodic, sensitive to bioamines, central hypoxia and inputs from CNS upper structures, and evolve with developmental age. In vitro, the foetal rodent RRG studied in preparations where the upper CNS structures are lacking continuously produces a rhythmic command, which is sensitive to hypoxia and bioaminergic inputs. The rhythm is slow with variable periods 4 days before birth. It becomes faster 2 days before birth, similar to the postnatal rhythm. Compelling evidence suggests that a region of the RRG called the preBötzinger complex (PBC) contains respiratory pacemaker neurones which play a primary role in perinatal rhythmogenesis. Although the RRG functions during early gestation, no pacemakers are found in the putative PBC area and its electrical stimulation and lesion do not affect the early foetal rhythm. To know whether the early foetal and perinatal rhythms originate from either pacemaker neurones or network connection properties, and to know which maturational processes might explain the appearance of PBC pacemakers and the rhythm increase during perinatal development, we computationally modelled maturing RRG. Our model shows that both network noise and persistent sodium conductance are crucial for rhythmogenesis and that a slight increase in the persistent sodium conductance can solve the pacemaker versus network dilemma in a noisy network.

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. Pontryagin's averaging technique is used to locate the periodic orbits in the phase space.

Whole cell stochastic model reproduces the irregularities found in the membrane potential of bursting neurons

Irregular intrinsic behavior of neurons seems ubiquitous in the nervous system. Even in circuits specialized to provide periodic and reliable patterns to control the repetitive activity of muscles, such as the pyloric central pattern generator (CPG) of the crustacean stomatogastric ganglion (STG), many bursting motor neurons present irregular activity when deprived from synaptic inputs. Moreover, many authors attribute to these irregularities the role of providing flexibility and adaptation capabilities to oscillatory neural networks such as CPGs. These irregular behaviors, related to nonlinear and chaotic properties of the cells, pose serious challenges to developing deterministic Hodgkin-Huxley-type (HH-type) conductance models. Only a few deterministic HH-type models based on experimental conductance values were able to show such nonlinear properties, but most of these models are based on slow oscillatory dynamics of the cytosolic calcium concentration that were never found experimentally in STG neurons. Based on an up-to-date single-compartment deterministic HH-type model of a STG neuron, we developed a stochastic HH-type model based on the microscopic Markovian states that an ion channel can achieve. We used tools from nonlinear analysis to show that the stochastic model is able to express the same kind of irregularities, sensitivity to initial conditions, and low dimensional dynamics found in the neurons isolated from the STG. Without including any nonrealistic dynamics in our whole cell stochastic model, we show that the nontrivial dynamics of the membrane potential naturally emerge from the interplay between the microscopic probabilistic character of the ion channels and the nonlinear interactions among these elements. Moreover, the experimental irregular behavior is reproduced by the stochastic model for the same parameters for which the membrane potential of the original deterministic model exhibits periodic oscillations.

Transition between tonic spiking and bursting in a neuron model via the blue-sky catastrophe

We study a continuous and reversible transition between periodic tonic spiking and bursting activities in a neuron model. It is described as the blue-sky catastrophe, which is a homoclinic bifurcation of a saddle-node periodic orbit of codimension one. This transition constitutes a biophysically plausible mechanism for the regulation of burst duration that increases with no bound like 1/square root alpha-alpha0 as the transition value alpha0 is approached.