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

Simple model to incorporate statistical noise based on a modified hodgkin-huxley approach for external electrical field driven neural responses

Noise activity is known to affect neural networks, enhance the system response to weak external signals, and lead to stochastic resonance phenomenon that can effectively amplify signals in nonlinear systems. In most treatments, channel noise has been modeled based on multi-state Markov descriptions or the use stochastic differential equation models. Here we probe a computationally simple approach based on a minor modification of the traditional Hodgkin-Huxley approach to embed noise in neural response. Results obtained from numerous simulations with different excitation frequencies and noise amplitudes for the action potential firing show very good agreement with output obtained from well-established models. Furthermore, results from the Mann-Whitney U Test reveal a statistically insignificant difference. The distribution of the time interval between successive potential spikes obtained from this simple approach compared very well with the results of complicated Fox and Lu type methods at much reduced computational cost. This present method could also possibly be applied to the analysis of spatial variations and/or differences in characteristics of random incident electromagnetic signals.

A thermodynamical model of non-deterministic computation in cortical neural networks

Neuronal populations in the cerebral cortex engage in probabilistic coding, effectively encoding the state of the surrounding environment with high accuracy and extraordinary energy efficiency. A new approach models the inherently probabilistic nature of cortical neuron signaling outcomes as a thermodynamic process of non-deterministic computation. A mean field approach is used, with the trial Hamiltonian maximizing available free energy and minimizing the net quantity of entropy, compared with a reference Hamiltonian. Thermodynamic quantities are always conserved during the computation; free energy must be expended to produce information, and free energy is released during information compression, as correlations are identified between the encoding system and its surrounding environment. Due to the relationship between the Gibbs free energy equation and the Nernst equation, any increase in free energy is paired with a local decrease in membrane potential. As a result, this process of thermodynamic computation adjusts the likelihood of each neuron firing an action potential. This model shows that non-deterministic signaling outcomes can be achieved by noisy cortical neurons, through an energy-efficient computational process that involves optimally redistributing a Hamiltonian over some time evolution. Calculations demonstrate that the energy efficiency of the human brain is consistent with this model of non-deterministic computation, with net entropy production far too low to retain the assumptions of a classical system.

Resolving molecular contributions of ion channel noise to interspike interval variability through stochastic shielding

Molecular fluctuations can lead to macroscopically observable effects. The random gating of ion channels in the membrane of a nerve cell provides an important example. The contributions of independent noise sources to the variability of action potential timing have not previously been studied at the level of molecular transitions within a conductance-based model ion-state graph. Here we study a stochastic Langevin model for the Hodgkin-Huxley (HH) system based on a detailed representation of the underlying channel state Markov process, the "[Formula: see text]D model" introduced in (Pu and Thomas in Neural Computation 32(10):1775-1835, 2020). We show how to resolve the individual contributions that each transition in the ion channel graph makes to the variance of the interspike interval (ISI). We extend the mean return time (MRT) phase reduction developed in (Cao et al. in SIAM J Appl Math 80(1):422-447, 2020) to the second moment of the return time from an MRT isochron to itself. Because fixed-voltage spike detection triggers do not correspond to MRT isochrons, the inter-phase interval (IPI) variance only approximates the ISI variance. We find the IPI variance and ISI variance agree to within a few percent when both can be computed. Moreover, we prove rigorously, and show numerically, that our expression for the IPI variance is accurate in the small noise (large system size) regime; our theory is exact in the limit of small noise. By selectively including the noise associated with only those few transitions responsible for most of the ISI variance, our analysis extends the stochastic shielding (SS) paradigm (Schmandt and Galán in Phys Rev Lett 109(11):118101, 2012) from the stationary voltage clamp case to the current clamp case. We show numerically that the SS approximation has a high degree of accuracy even for larger, physiologically relevant noise levels. Finally, we demonstrate that the ISI variance is not an unambiguously defined quantity, but depends on the choice of voltage level set as the spike detection threshold. We find a small but significant increase in ISI variance, the higher the spike detection voltage, both for simulated stochastic HH data and for voltage traces recorded in in vitro experiments. In contrast, the IPI variance is invariant with respect to the choice of isochron used as a trigger for counting "spikes."

Voltage distributions in extracellular brain recordings

Extracellular recordings of brain voltage signals have many uses, including the identification of spikes and the characterization of brain states via analysis of local field potential (LFP) or EEG recordings. Though the factors underlying the generation of these signals are time varying and complex, their analysis may be facilitated by an understanding of their statistical properties. To this end, we analyzed the voltage distributions of high-pass extracellular recordings from a variety of structures, including cortex, thalamus, and hippocampus, in monkeys, cats, and rodents. We additionally investigated LFP signals in these recordings as well as human EEG signals obtained during different sleep stages. In all cases, the distributions were accurately described by a Gaussian within ±1.5 standard deviations from zero. Outside these limits, voltages tended to be distributed exponentially, that is, they fell off linearly on log-linear frequency plots, with variable heights and slopes. A possible explanation for this is that sporadically and independently occurring events with individual Gaussian size distributions can sum to produce approximately exponential distributions. For the high-pass recordings, a second explanation results from a model of the noisy behavior of ion channels that produce action potentials via Hodgkin-Huxley kinetics. The distributions produced by this model, relative to the averaged potential, were also Gaussian with approximately exponential flanks. The model also predicted time-varying noise distributions during action potentials, which were observed in the extracellular spike signals. These findings suggest a principled method for detecting spikes in high-pass recordings and transient events in LFP and EEG signals.NEW & NOTEWORTHY We show that the voltage distributions in brain recordings, including high-pass extracellular recordings, the LFP, and human EEG, are accurately described by a Gaussian within ±1.5 standard deviations from zero, with heavy, exponential tails outside these limits. This offers a principled way of setting event detection thresholds in high-pass recordings. It also offers a means for identifying event-like, transient signals in LFP and EEG recordings which may correlate with other neural phenomena.

Firing statistics in the bistable regime of neurons with homoclinic spike generation

Neuronal voltage dynamics of regularly firing neurons typically has one stable attractor: either a fixed point (like in the subthreshold regime) or a limit cycle that defines the tonic firing of action potentials (in the suprathreshold regime). In two of the three spike onset bifurcation sequences that are known to give rise to all-or-none type action potentials, however, the resting-state fixed point and limit cycle spiking can coexist in an intermediate regime, resulting in bistable dynamics. Here, noise can induce switches between the attractors, i.e., between rest and spiking, and thus increase the variability of the spike train compared to neurons with only one stable attractor. Qualitative features of the resulting spike statistics depend on the spike onset bifurcations. This paper focuses on the creation of the spiking limit cycle via the saddle-homoclinic orbit (HOM) bifurcation and derives interspike interval (ISI) densities for a conductance-based neuron model in the bistable regime. The ISI densities of bistable homoclinic neurons are found to be unimodal yet distinct from the inverse Gaussian distribution associated with the saddle-node-on-invariant-cycle bifurcation. It is demonstrated that for the HOM bifurcation the transition between rest and spiking is mainly determined along the downstroke of the action potential-a dynamical feature that is not captured by the commonly used reset neuron models. The deduced spike statistics can help to identify HOM dynamics in experimental data.

The case for emulating insect brains using anatomical "wiring diagrams" equipped with biophysical models of neuronal activity

Developing whole-brain emulation (WBE) technology would provide immense benefits across neuroscience, biomedicine, artificial intelligence, and robotics. At this time, constructing a simulated human brain lacks feasibility due to limited experimental data and limited computational resources. However, I suggest that progress toward this goal might be accelerated by working toward an intermediate objective, namely insect brain emulation (IBE). More specifically, this would entail creating biologically realistic simulations of entire insect nervous systems along with more approximate simulations of non-neuronal insect physiology to make "virtual insects." I argue that this could be realistically achievable within the next 20 years. I propose that developing emulations of insect brains will galvanize the global community of scientists, businesspeople, and policymakers toward pursuing the loftier goal of emulating the human brain. By demonstrating that WBE is possible via IBE, simulating mammalian brains and eventually the human brain may no longer be viewed as too radically ambitious to deserve substantial funding and resources. Furthermore, IBE will facilitate dramatic advances in cognitive neuroscience, artificial intelligence, and robotics through studies performed using virtual insects.

Fluctuation Scaling of Neuronal Firing and Bursting in Spontaneously Active Brain Circuits

We employed high-density microelectrode arrays to investigate spontaneous firing patterns of neurons in brain circuits of the primary somatosensory cortex (S1) in mice. We recorded from over 150 neurons for 10min in each of eight different experiments, identified their location in S1, sorted their action potentials (spikes), and computed their power spectra and inter-spike interval (ISI) statistics. Of all persistently active neurons, 92% fired with a single dominant frequency - regularly firing neurons (RNs) - from 1 to 8Hz while 8% fired in burst with two dominant frequencies - bursting neurons (BNs) - corresponding to the inter-burst (2-6Hz) and intra-burst intervals (20-160Hz). RNs were predominantly located in layers 2/3 and 5/6 while BNs localized to layers 4 and 5. Across neurons, the standard deviation of ISI was a power law of its mean, a property known as fluctuation scaling, with a power law exponent of 1 for RNs and 1.25 for BNs. The power law implies that firing and bursting patterns are scale invariant: the firing pattern of a given RN or BN resembles that of another RN or BN, respectively, after a time contraction or dilation. An explanation for this scale invariance is discussed in the context of previous computational studies as well as its potential role in information processing.

Data-driven inference for stationary jump-diffusion processes with application to membrane voltage fluctuations in pyramidal neurons

The emergent activity of biological systems can often be represented as low-dimensional, Langevin-type stochastic differential equations. In certain systems, however, large and abrupt events occur and violate the assumptions of this approach. We address this situation here by providing a novel method that reconstructs a jump-diffusion stochastic process based solely on the statistics of the original data. Our method assumes that these data are stationary, that diffusive noise is additive, and that jumps are Poisson. We use threshold-crossing of the increments to detect jumps in the time series. This is followed by an iterative scheme that compensates for the presence of diffusive fluctuations that are falsely detected as jumps. Our approach is based on probabilistic calculations associated with these fluctuations and on the use of the Fokker-Planck and the differential Chapman-Kolmogorov equations. After some validation cases, we apply this method to recordings of membrane noise in pyramidal neurons of the electrosensory lateral line lobe of weakly electric fish. These recordings display large, jump-like depolarization events that occur at random times, the biophysics of which is unknown. We find that some pyramidal cells increase their jump rate and noise intensity as the membrane potential approaches spike threshold, while their drift function and jump amplitude distribution remain unchanged. As our method is fully data-driven, it provides a valuable means to further investigate the functional role of these jump-like events without relying on unconstrained biophysical models.

Double inverse stochastic resonance with dynamic synapses

We investigate the behavior of a model neuron that receives a biophysically realistic noisy postsynaptic current based on uncorrelated spiking activity from a large number of afferents. We show that, with static synapses, such noise can give rise to inverse stochastic resonance (ISR) as a function of the presynaptic firing rate. We compare this to the case with dynamic synapses that feature short-term synaptic plasticity and show that the interval of presynaptic firing rate over which ISR exists can be extended or diminished. We consider both short-term depression and facilitation. Interestingly, we find that a double inverse stochastic resonance (DISR), with two distinct wells centered at different presynaptic firing rates, can appear.

Effects of ion channel noise on neural circuits: an application to the respiratory pattern generator to investigate breathing variability

Neural activity generally displays irregular firing patterns even in circuits with apparently regular outputs, such as motor pattern generators, in which the output frequency fluctuates randomly around a mean value. This "circuit noise" is inherited from the random firing of single neurons, which emerges from stochastic ion channel gating (channel noise), spontaneous neurotransmitter release, and its diffusion and binding to synaptic receptors. Here we demonstrate how to expand conductance-based network models that are originally deterministic to include realistic, physiological noise, focusing on stochastic ion channel gating. We illustrate this procedure with a well-established conductance-based model of the respiratory pattern generator, which allows us to investigate how channel noise affects neural dynamics at the circuit level and, in particular, to understand the relationship between the respiratory pattern and its breath-to-breath variability. We show that as the channel number increases, the duration of inspiration and expiration varies, and so does the coefficient of variation of the breath-to-breath interval, which attains a minimum when the mean duration of expiration slightly exceeds that of inspiration. For small channel numbers, the variability of the expiratory phase dominates over that of the inspiratory phase, and vice versa for large channel numbers. Among the four different cell types in the respiratory pattern generator, pacemaker cells exhibit the highest sensitivity to channel noise. The model shows that suppressing input from the pons leads to longer inspiratory phases, a reduction in breathing frequency, and larger breath-to-breath variability, whereas enhanced input from the raphe nucleus increases breathing frequency without changing its pattern.

NEW & NOTEWORTHY

A major source of noise in neuronal circuits is the "flickering" of ion currents passing through the neurons' membranes (channel noise), which cannot be suppressed experimentally. Computational simulations are therefore the best way to investigate the effects of this physiological noise by manipulating its level at will. We investigate the role of noise in the respiratory pattern generator and show that endogenous, breath-to-breath variability is tightly linked to the respiratory pattern.

Stochastic and deterministic dynamics of intrinsically irregular firing in cortical inhibitory interneurons

Most cortical neurons fire regularly when excited by a constant stimulus. In contrast, irregular-spiking (IS) interneurons are remarkable for the intrinsic variability of their spike timing, which can synchronize amongst IS cells via specific gap junctions. Here, we have studied the biophysical mechanisms of this irregular spiking in mice, and how IS cells fire in the context of synchronous network oscillations. Using patch-clamp recordings, artificial dynamic conductance injection, pharmacological analysis and computational modeling, we show that spike time irregularity is generated by a nonlinear dynamical interaction of voltage-dependent sodium and fast-inactivating potassium channels just below spike threshold, amplifying channel noise. This active irregularity may help IS cells synchronize with each other at gamma range frequencies, while resisting synchronization to lower input frequencies.

DOI:http://dx.doi.org/10.7554/eLife.16475.001

Neurons send information to other neurons in the brain by generating fast electrical pulses called action potentials (or spikes). When stimulated by input signals of a constant size, neurons generally respond with regular patterns of spiking leading to rhythmical brain activity. However, neurons known as irregular spiking interneurons are unique: the relationship between the input they receive and whether or not they produce a spike appears to be random. The molecular mechanism behind this phenomenon is not clear.

Mendonça et al. set out to investigate whether irregular spiking is truly random, or whether there is some degree of predictability. The experiments used genetically modified mice in which irregular spiking interneurons were specifically labeled with a fluorescent protein, which made them easier to find to record their electrical activity. Sophisticated statistical analyses showed that these neurons are not firing at random. Instead, there is a pattern to the timings of the spikes they produce.

It was previously known that electrical spikes in neurons are generated by sodium ions and potassium ions moving across the membrane that surrounds each cell. Proteins called ion channels provide routes for these ions to pass through the membrane. Mendonça et al. show that compared to other neurons, irregular spiking interneurons have larger numbers of a specific type of potassium ion channel. Mimicking the effect of increasing the number of these potassium ion channels in the interneurons made the firing pattern of these neurons more irregular, while decreasing the number of these channels made the firing patterns more regular and predictable.

A computer model of an irregular spiking interneuron showed that the activity of these potassium ion channels and a type of sodium ion channel plays a key role in producing irregular electrical spiking. Further analysis showed that irregular spiking interneurons can synchronize their activity with fast, but not slow, rhythms in brain activity.

The findings of Mendonça et al. suggest that irregular spiking interneurons can disrupt slow regular electrical activity in the brain. Rhythms in brain activity vary depending on whether we are awake or asleep, and are altered in diseases such as epilepsy and schizophrenia. Now that we have a better understanding of how irregular spiking interneurons work, it should be possible to find out how they coordinate their activity with each other, and what effect they have on animal behavior.

DOI:http://dx.doi.org/10.7554/eLife.16475.002

Fluid limit theorems for stochastic hybrid systems with application to neuron models

In this paper we establish limit theorems for a class of stochastic hybrid systems (continuous deterministic dynamics coupled with jump Markov processes) in the fluid limit (small jumps at high frequency), thus extending known results for jump Markov processes. We prove a functional law of large numbers with exponential convergence speed, derive a diffusion approximation, and establish a functional central limit theorem. We apply these results to neuron models with stochastic ion channels, as the number of channels goes to infinity, estimating the convergence to the deterministic model. In terms of neural coding, we apply our central limit theorems to numerically estimate the impact of channel noise both on frequency and spike timing coding.

Fluid limit theorems for stochastic hybrid systems with application to neuron models

In this paper we establish limit theorems for a class of stochastic hybrid systems (continuous deterministic dynamics coupled with jump Markov processes) in the fluid limit (small jumps at high frequency), thus extending known results for jump Markov processes. We prove a functional law of large numbers with exponential convergence speed, derive a diffusion approximation, and establish a functional central limit theorem. We apply these results to neuron models with stochastic ion channels, as the number of channels goes to infinity, estimating the convergence to the deterministic model. In terms of neural coding, we apply our central limit theorems to numerically estimate the impact of channel noise both on frequency and spike timing coding.

Distribution of transition times in a stochastic model of excitable cell: Insights into the cell-intrinsic mechanisms of randomness in neuronal interspike intervals

In this paper, we develop an analytical approach to studying random patterns of activity in excitable cells. Our analytical approach uses a two-state stochastic model of excitable system based on the electrophysiological properties of refractoriness and restitution, which characterize cell recovery after excitation. By applying the notion of probability density flux, we derive the distributions of transition times between states and the distribution of interspike interval (ISI) durations for a constant applied stimulus. The derived ISI distribution is unimodal and, provided that the time spent in the excited state is constant, can be approximated by a Rayleigh peak followed by an exponential tail. We then explore the role of the model parameters in determining the shape of the derived distributions and the ISI coefficient of variation. Finally, we use our analytical results to study simulation results from the stochastic Morris-Lecar neuron and from a three-state extension of the proposed stochastic model, which is capable of reproducing multimodal ISI histograms.

Channel based generating function approach to the stochastic Hodgkin-Huxley neuronal system

Internal and external fluctuations, such as channel noise and synaptic noise, contribute to the generation of spontaneous action potentials in neurons. Many different Langevin approaches have been proposed to speed up the computation but with waning accuracy especially at small channel numbers. We apply a generating function approach to the master equation for the ion channel dynamics and further propose two accelerating algorithms, with an accuracy close to the Gillespie algorithm but with much higher efficiency, opening the door for expedited simulation of noisy action potential propagating along axons or other types of noisy signal transduction.

Stochastic sensitivity analysis of noise-induced suppression of firing and giant variability of spiking in a Hodgkin-Huxley neuron model

We study the stochastic dynamics of a Hodgkin-Huxley neuron model in a regime of coexistent stable equilibrium and a limit cycle. In this regime, noise may suppress periodic firing by switching the neuron randomly to a quiescent state. We show that at a critical value of the injected current, the mean firing rate depends weakly on noise intensity, while the neuron exhibits giant variability of the interspike intervals and spike count. To reveal the dynamical origin of this noise-induced effect, we develop the stochastic sensitivity analysis and use the Mahalanobis metric for this four-dimensional stochastic dynamical system. We show that the critical point of giant variability corresponds to the matching of the Mahalanobis distances from attractors (stable equilibrium and limit cycle) to a three-dimensional surface separating their basins of attraction.

Channel-noise-induced critical slowing in the subthreshold Hodgkin-Huxley neuron

The dynamics of a spiking neuron approaching threshold is investigated in the framework of Markov-chain models describing the random state-transitions of the underlying ion-channel proteins. We characterize subthreshold channel-noise-induced transmembrane potential fluctuations in both type-I (integrator) and type-II (resonator) parametrizations of the classic conductance-based Hodgkin-Huxley equations. As each neuron approaches spiking threshold from below, numerical simulations of stochastic trajectories demonstrate pronounced growth in amplitude simultaneous with decay in frequency of membrane voltage fluctuations induced by ion-channel state transitions. To explore this progression of fluctuation statistics, we approximate the exact Markov treatment with a 12-variable channel-based stochastic differential equation (SDE) and its Ornstein-Uhlenbeck (OU) linearization and show excellent agreement between Markov and SDE numerical simulations. Predictions of the OU theory with respect to membrane potential fluctuation variance, autocorrelation, correlation time, and spectral density are also in agreement and illustrate the close connection between the eigenvalue structure of the associated deterministic bifurcations and the observed behavior of the noisy Markov traces on close approach to threshold for both integrator and resonator point-neuron varieties.

Characterizing spiking in noisy type II neurons

Understanding the dynamics of noisy neurons remains an important challenge in neuroscience. Here, we describe a simple probabilistic model that accurately describes the firing behavior in a large class (type II) of neurons. To demonstrate the usefulness of this model, we show how it accurately predicts the interspike interval (ISI) distributions, bursting patterns and mean firing rates found by: (1) simulations of the classic Hodgkin-Huxley model with channel noise, (2) experimental data from squid giant axon with a noisy input current and (3) experimental data on noisy firing from a neuron within the suprachiasmatic nucleus (SCN). This simple model has 6 parameters, however, in some cases, two of these parameters are coupled and only 5 parameters account for much of the known behavior. From these parameters, many properties of spiking can be found through simple calculation. Thus, we show how the complex effects of noise can be understood through a simple and general probabilistic model.

The ISI distribution of the stochastic Hodgkin-Huxley neuron

The simulation of ion-channel noise has an important role in computational neuroscience. In recent years several approximate methods of carrying out this simulation have been published, based on stochastic differential equations, and all giving slightly different results. The obvious, and essential, question is: which method is the most accurate and which is most computationally efficient? Here we make a contribution to the answer. We compare interspike interval histograms from simulated data using four different approximate stochastic differential equation (SDE) models of the stochastic Hodgkin-Huxley neuron, as well as the exact Markov chain model simulated by the Gillespie algorithm. One of the recent SDE models is the same as the Kurtz approximation first published in 1978. All the models considered give similar ISI histograms over a wide range of deterministic and stochastic input. Three features of these histograms are an initial peak, followed by one or more bumps, and then an exponential tail. We explore how these features depend on deterministic input and on level of channel noise, and explain the results using the stochastic dynamics of the model. We conclude with a rough ranking of the four SDE models with respect to the similarity of their ISI histograms to the histogram of the exact Markov chain model.

Channel-noise-induced stochastic facilitation in an auditory brainstem neuron model

Neuronal membrane potentials fluctuate stochastically due to conductance changes caused by random transitions between the open and closed states of ion channels. Although it has previously been shown that channel noise can nontrivially affect neuronal dynamics, it is unknown whether ion-channel noise is strong enough to act as a noise source for hypothesized noise-enhanced information processing in real neuronal systems, i.e., "stochastic facilitation". Here we demonstrate that biophysical models of channel noise can give rise to two kinds of recently discovered stochastic facilitation effects in a Hodgkin-Huxley-like model of auditory brainstem neurons. The first, known as slope-based stochastic resonance (SBSR), enables phasic neurons to emit action potentials that can encode the slope of inputs that vary slowly relative to key time constants in the model. The second, known as inverse stochastic resonance (ISR), occurs in tonically firing neurons when small levels of noise inhibit tonic firing and replace it with burstlike dynamics. Consistent with previous work, we conclude that channel noise can provide significant variability in firing dynamics, even for large numbers of channels. Moreover, our results show that possible associated computational benefits may occur due to channel noise in neurons of the auditory brainstem. This holds whether the firing dynamics in the model are phasic (SBSR can occur due to channel noise) or tonic (ISR can occur due to channel noise).

Stochastic sensitivity analysis of the noise-induced excitability in a model of a hair bundle

We study effect of weak noise on the dynamics of a hair bundle model near the excitability threshold and near a subcritical Hopf bifurcation. We analyze numerically noise-induced structural changes in the probability density and the power spectral density of the model. In particular, we show that weak noise can induce oscillations with two distinct frequencies in both excitable and limit-cycle regimes. We then applied a recently developed technique of stochastic sensitivity functions which allows us to estimate threshold values of noise intensity corresponding to these transitions.

Origin of intrinsic irregular firing in cortical interneurons

Cortical spike trains are highly irregular both during ongoing, spontaneous activity and when driven at high firing rates. There is uncertainty about the source of this irregularity, ranging from intrinsic noise sources in neurons to collective effects in large-scale cortical networks. Cortical interneurons display highly irregular spike times (coefficient of variation of the interspike intervals >1) in response to dc-current injection in vitro. This is in marked contrast to cortical pyramidal cells, which spike highly irregularly in vivo, but regularly in vitro. We show with in vitro recordings and computational models that this is due to the fast activation kinetics of interneuronal K(+) currents. This explanation holds over a wide parameter range and with Gaussian white, power-law, and Ornstein-Uhlenbeck noise. The intrinsically irregular spiking of interneurons could contribute to the irregularity of the cortical network.

Channel-based Langevin approach for the stochastic Hodgkin-Huxley neuron

Stochasticity in ion channel gating is the major source of intrinsic neuronal noise, which can induce many important effects in neuronal dynamics. Several numerical implementations of the Langevin approach have been proposed to approximate the Markovian dynamics of the Hodgkin-Huxley neuronal model. In this work an improved channel-based Langevin approach is proposed by introducing a truncation procedure to limit the state fractions in the range of [0, 1]. The truncated fractions are put back into the state fractions in the next time step for channel noise calculation. Our simulations show that the bounded Langevin approaches combined with the restored process give better approximations to the statistics of action potentials with the Markovian method. As a result, in our approach the channel state fractions are disturbed by two terms of noise: an uncorrelated Gaussian noise and a time-correlated noise obtained from the truncated fractions. We suggest that the restoration of truncated fractions is a critical process for a bounded Langevin method.

Conductance-based neuron models and the slow dynamics of excitability

In recent experiments, synaptically isolated neurons from rat cortical culture, were stimulated with periodic extracellular fixed-amplitude current pulses for extended durations of days. The neuron’s response depended on its own history, as well as on the history of the input, and was classified into several modes. Interestingly, in one of the modes the neuron behaved intermittently, exhibiting irregular firing patterns changing in a complex and variable manner over the entire range of experimental timescales, from seconds to days. With the aim of developing a minimal biophysical explanation for these results, we propose a general scheme, that, given a few assumptions (mainly, a timescale separation in kinetics) closely describes the response of deterministic conductance-based neuron models under pulse stimulation, using a discrete time piecewise linear mapping, which is amenable to detailed mathematical analysis. Using this method we reproduce the basic modes exhibited by the neuron experimentally, as well as the mean response in each mode. Specifically, we derive precise closed-form input-output expressions for the transient timescale and firing rates, which are expressed in terms of experimentally measurable variables, and conform with the experimental results. However, the mathematical analysis shows that the resulting firing patterns in these deterministic models are always regular and repeatable (i.e., no chaos), in contrast to the irregular and variable behavior displayed by the neuron in certain regimes. This fact, and the sensitive near-threshold dynamics of the model, indicate that intrinsic ion channel noise has a significant impact on the neuronal response, and may help reproduce the experimentally observed variability, as we also demonstrate numerically. In a companion paper, we extend our analysis to stochastic conductance-based models, and show how these can be used to reproduce the details of the observed irregular and variable neuronal response.

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.

Stochastic differential equation models for ion channel noise in Hodgkin-Huxley neurons

The random transitions of ion channels between conducting and nonconducting states generate a source of internal fluctuations in a neuron, known as channel noise. The standard method for modeling the states of ion channels nonlinearly couples continuous-time Markov chains to a differential equation for voltage. Beginning with the work of R. F. Fox and Y.-N. Lu [Phys. Rev. E 49, 3421 (1994)], there have been attempts to generate simpler models that use stochastic differential equation (SDEs) to approximate the stochastic spiking activity produced by Markov chain models. Recent numerical investigations, however, have raised doubts that SDE models can capture the stochastic dynamics of Markov chain models.We analyze three SDE models that have been proposed as approximations to the Markov chain model: one that describes the states of the ion channels and two that describe the states of the ion channel subunits. We show that the former channel-based approach can capture the distribution of channel noise and its effects on spiking in a Hodgkin-Huxley neuron model to a degree not previously demonstrated, but the latter two subunit-based approaches cannot. Our analysis provides intuitive and mathematical explanations for why this is the case. The temporal correlation in the channel noise is determined by the combinatorics of bundling subunits into channels, but the subunit-based approaches do not correctly account for this structure. Our study confirms and elucidates the findings of previous numerical investigations of subunit-based SDE models. Moreover, it presents evidence that Markov chain models of the nonlinear, stochastic dynamics of neural membranes can be accurately approximated by SDEs. This finding opens a door to future modeling work using SDE techniques to further illuminate the effects of ion channel fluctuations on electrically active cells.

Multimodal transition and stochastic antiresonance in squid giant axons

The experimental data of Takahashi [Physica D 43, 318 (1990)], on the response of squid giant axons stimulated by periodic sequence of short current pulses is interpreted within the Hodgkin-Huxley model. The minimum of the firing rate as a function of the stimulus amplitude I0 in the high-frequency regime is due to the multimodal transition. Below this singular point only odd multiples of the driving period remain and the system is sensitive to noise. The coefficient of variation has a maximum and the firing rate has a minimum as a function of the noise intensity, which is an indication of the stochastic coherence antiresonance. The model calculations reproduce the frequency of occurrence of the most common modes in the vicinity of the transition. A linear relation of output frequency vs I0 above the transition is also confirmed.

Intrinsic variability of latency to first-spike

The assessment of the variability of neuronal spike timing is fundamental to gain understanding of latency coding. Based on recent mathematical results, we investigate theoretically the impact of channel noise on latency variability. For large numbers of ion channels, we derive the asymptotic distribution of latency, together with an explicit expression for its variance. Consequences in terms of information processing are studied with Fisher information in the Morris-Lecar model. A competition between sensitivity to input and precision is responsible for favoring two distinct regimes of latencies.

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.

Correlation analysis a tool for comparing relaxation-type models to experimental data

We describe a new technique for comparing mathematical models to the biological systems that are described. This technique is appropriate for systems that produce relaxation oscillations or bursting oscillations, and takes advantage of noise that is inherent to all biological systems. Both types of oscillations are composed of active phases of activity followed by silent phases, repeating periodically. The presence of noise adds variability to the durations of the different phases. The central idea of the technique is that the active phase duration may be correlated with either/both the previous or next silent phase duration, and the resulting correlation pattern provides information about the dynamic structure of the system. Correlation patterns can easily be determined by making scatter plots and applying correlation analysis to the cluster of data points. This could be done both with experimental data and with model simulation data. If the model correlation pattern is in general agreement with the experimental data, then this adds support for the validity of the model. Otherwise, the model must be corrected. While this tool is only one test of many required to validate a mathematical model, it is easy to implement and is noninvasive.

Discrimination with Spike Times and ISI Distributions

We study the discrimination capability of spike time sequences using the Chernoff distance as a metric. We assume that spike sequences are generated by renewal processes and study how the Chernoff distance depends on the shape of interspike interval (ISI) distribution. First, we consider a lower bound to the Chernoff distance because it has a simple closed form. Then we consider specific models of ISI distributions such as the gamma, inverse gaussian (IG), exponential with refractory period (ER), and that of the leaky integrate-and-fire (LIF) neuron. We found that the discrimination capability of spike times strongly depends on high-order moments of ISI and that it is higher when the spike time sequence has a larger skewness and a smaller kurtosis. High variability in terms of coefficient of variation (CV) does not necessarily mean that the spike times have less discrimination capability. Spike sequences generated by the gamma distribution have the minimum discrimination capability for a given mean and variance of ISI. We used series expansions to calculate the mean and variance of ISIs for LIF neurons as a function of the mean input level and the input noise variance. Spike sequences from an LIF neuron are more capable of discrimination than those of IG and gamma distributions when the stationary voltage level is close to the neuron's threshold value of the neuron.

Irregular firing of isolated cortical interneurons in vitro driven by intrinsic stochastic mechanisms

Pharmacologically isolated GABAergic irregular spiking and stuttering interneurons in the mouse visual cortex display highly irregular spike times, with high coefficients of variation approximately 0.9-3, in response to a depolarizing, constant current input. This is in marked contrast to cortical pyramidal cells, which spike quite regularly in response to the same current injection. We applied time-series analysis methods to show that the irregular behavior of the interneurons was not a consequence of low-dimensional, deterministic processes. These methods were also applied to the Hindmarsh and Rose neuronal model to confirm that the methods are adequate for the types of data under investigation. This result has important consequences for the origin of fluctuations observed in the cortex in vivo.