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

Inverse stochastic resonance in adaptive small-world neural networks

Inverse stochastic resonance (ISR) is a counterintuitive phenomenon where noise reduces the oscillation frequency of an oscillator to a minimum occurring at an intermediate noise intensity, and sometimes even to the complete absence of oscillations. In neuroscience, ISR was first experimentally verified with cerebellar Purkinje neurons [Buchin et al., PLOS Comput. Biol. 12, e1005000 (2016)]. These experiments showed that ISR enables a locally optimal information transfer between the input and output spike train of neurons. Subsequent studies have further demonstrated the efficiency of information processing and transfer in neural networks with small-world network topology. We have conducted a numerical investigation into the impact of adaptivity on ISR in a small-world network of noisy FitzHugh-Nagumo (FHN) neurons, operating in a bi-metastable regime consisting of a metastable fixed point and a metastable limit cycle. Our results show that the degree of ISR is highly dependent on the value of the FHN model's timescale separation parameter ε. The network structure undergoes dynamic adaptation via mechanisms of either spike-time-dependent plasticity (STDP) with potentiation-/depression-domination parameter P or homeostatic structural plasticity (HSP) with rewiring frequency F. We demonstrate that both STDP and HSP amplify the effect of ISR when ε lies within the bi-stability region of FHN neurons. Specifically, at larger values of ε within the bi-stability regime, higher rewiring frequencies F are observed to enhance ISR at intermediate (weak) synaptic noise intensities, while values of P consistent with depression-domination (potentiation-domination) consistently enhance (deteriorate) ISR. Moreover, although STDP and HSP control parameters may jointly enhance ISR, P has a greater impact on improving ISR compared to F. Our findings inform future ISR enhancement strategies in noisy artificial neural circuits, aiming to optimize local information transfer between input and output spike trains in neuromorphic systems and prompt venues for experiments in neural networks.

Characterizing the information transmission of inverse stochastic resonance and noise-induced activity amplification in neuronal systems

Purkinje cells exhibit a reduction of the mean firing rate at intermediate-noise intensities, which is somewhat reminiscent of the response enhancement known as "stochastic resonance" (SR). Although the comparison with the stochastic resonance ends here, the current phenomenon has been given the name "inverse stochastic resonance" (ISR). Recent research has demonstrated that the ISR effect, like its close relative "nonstandard SR" [or, more correctly, noise-induced activity amplification (NIAA)], has been shown to stem from the weak-noise quenching of the initial distribution, in bistable regimes where the metastable state has a larger attraction basin than the global minimum. To understand the underlying mechanism of the ISR and NIAA phenomena, we study the probability distribution function of a one-dimensional system subjected to a bistable potential that has the property of symmetry, i.e., if we change the sign of one of its parameters, we can obtain both phenomena with the same properties in the depth of the wells and the width of their basins of attraction subjected to Gaussian white noise with variable intensity. Previous work has shown that one can theoretically determine the probability distribution function using the convex sum between the behavior at small and high noise intensities. To determine the probability distribution function more precisely, we resort to the "weighted ensemble Brownian dynamics simulation" model, which provides an accurate estimate of the probability distribution function for both low and high noise intensities and, most importantly, for the transition of both behaviors. In this way, on the one hand, we show that both phenomena emerge from a metastable system where, in the case of ISR, the global minimum of the system is in a state of lower activity, while in the case of NIAA, the global minimum is in a state of increased activity, the importance of which does not depend on the width of the basins of attraction. On the other hand, we see that quantifiers such as Fisher information, statistical complexity, and especially Shannon entropy fail to distinguish them, but they show the existence of the mentioned phenomena. Thus, noise management may well be a mechanism by which Purkinje cells find an efficient way to transmit information in the cerebral cortex.

Noise-tuned bursting in a Hedgehog burster

Noise can shape the firing behaviors of neurons. Here, we show that noise acting on the fast variable of the Hedgehog burster can tune the spike counts of bursts via the self-induced stochastic resonance (SISR) phenomenon. Using the distance matching condition, the critical transition positions on the slow manifolds can be predicted and the stochastic periodic orbits for various noise strengths are obtained. The critical transition positions on the slow manifold with non-monotonic potential differences exhibit a staircase-like dependence on the noise strength, which is also revealed by the stepwise change in the period of the stochastic periodic orbit. The noise-tuned bursting is more coherent within each step while displaying mixed-mode oscillations near the boundaries between the steps. When noise is large enough, noise-induced trapping of the slow variable can be observed, where the number of coexisting traps increases with the noise strength. It is argued that the robustness of SISR underlies the generality of the results discovered in this paper.

Efficient sampling and noisy decisions

Human decisions are based on finite information, which makes them inherently imprecise. But what determines the degree of such imprecision? Here, we develop an efficient coding framework for higher-level cognitive processes in which information is represented by a finite number of discrete samples. We characterize the sampling process that maximizes perceptual accuracy or fitness under the often-adopted assumption that full adaptation to an environmental distribution is possible, and show how the optimal process differs when detailed information about the current contextual distribution is costly. We tested this theory on a numerosity discrimination task, and found that humans efficiently adapt to contextual distributions, but in the way predicted by the model in which people must economize on environmental information. Thus, understanding decision behavior requires that we account for biological restrictions on information coding, challenging the often-adopted assumption of precise prior knowledge in higher-level decision systems.

Two paradigmatic scenarios for inverse stochastic resonance

Inverse stochastic resonance comprises a nonlinear response of an oscillatory system to noise where the frequency of noise-perturbed oscillations becomes minimal at an intermediate noise level. We demonstrate two generic scenarios for inverse stochastic resonance by considering a paradigmatic model of two adaptively coupled stochastic active rotators whose local dynamics is close to a bifurcation threshold. In the first scenario, shown for the two rotators in the excitable regime, inverse stochastic resonance emerges due to a biased switching between the oscillatory and the quasi-stationary metastable states derived from the attractors of the noiseless system. In the second scenario, illustrated for the rotators in the oscillatory regime, inverse stochastic resonance arises due to a trapping effect associated with a noise-enhanced stability of an unstable fixed point. The details of the mechanisms behind the resonant effect are explained in terms of slow-fast analysis of the corresponding noiseless systems.

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.

Bistability in Hodgkin-Huxley-type equations

We study how initial conditions of the Hodgkin-Huxley model affect the dynamics of simulated neurons. We systematically vary the amplitudes of depolarization currents in order to bring neuron dynamics to stable equilibrium. Our results demonstrate that simulated neurons can have spontaneous spiking or a silent state, depending on the initial conditions. We propose the methodology to study the circumstances under which Purkinje cells transit between hyperpolarized quiescent state (down state) and a depolarized spiking state (up state). We show that results derived using the Hodgkin-Huxley methodology should be carefully analyzed before suggesting a direct relevance to neuroprosthetic implants.

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.

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.

Accurate Langevin approaches to simulate Markovian channel dynamics

The stochasticity of ion-channels dynamic is significant for physiological processes on neuronal cell membranes. Microscopic simulations of the ion-channel gating with Markov chains can be considered to be an accurate standard. However, such Markovian simulations are computationally demanding for membrane areas of physiologically relevant sizes, which makes the noise-approximating or Langevin equation methods advantageous in many cases. In this review, we discuss the Langevin-like approaches, including the channel-based and simplified subunit-based stochastic differential equations proposed by Fox and Lu, and the effective Langevin approaches in which colored noise is added to deterministic differential equations. In the framework of Fox and Lu's classical models, several variants of numerical algorithms, which have been recently developed to improve accuracy as well as efficiency, are also discussed. Through the comparison of different simulation algorithms of ion-channel noise with the standard Markovian simulation, we aim to reveal the extent to which the existing Langevin-like methods approximate results using Markovian methods. Open questions for future studies are also discussed.

A review of methods for identifying stochastic resonance in simulations of single neuron models

Stochastic resonance (SR) is said to be observed when the presence of noise in a nonlinear system enables an output signal from the system to better represent some feature of an input signal than it does in the absence of noise. The effect has been observed in models of individual neurons, and in experiments performed on real neural systems. Despite the ubiquity of biophysical sources of stochastic noise in the nervous system, however, it has not yet been established whether neuronal computation mechanisms involved in performance of specific functions such as perception or learning might exploit such noise as an integral component, such that removal of the noise would diminish performance of these functions. In this paper we revisit the methods used to demonstrate stochastic resonance in models of single neurons. This includes a previously unreported observation in a multicompartmental model of a CA1-pyramidal cell. We also discuss, as a contrast to these classical studies, a form of 'stochastic facilitation', known as inverse stochastic resonance. We draw on the reviewed examples to argue why new approaches to studying 'stochastic facilitation' in neural systems need to be developed.

Modeling the influence of short term depression in vesicle release and stochastic calcium channel gating on auditory nerve spontaneous firing statistics

We propose several modifications to an existing computational model of stochastic vesicle release in inner hair cell ribbon synapses, with the aim of producing simulated auditory nerve fiber spiking data that more closely matches empirical data. Specifically, we studied the inter-spike-interval (ISI) distribution, and long and short term ISI correlations in spontaneous spiking in post-synaptic auditory nerve fibers. We introduced short term plasticity to the pre-synaptic release probability, in a manner analogous to standard stochastic models of cortical short term synaptic depression. This modification resulted in a similar distribution of vesicle release intervals to that estimated from empirical data. We also introduced a biophysical stochastic model of calcium channel opening and closing, but showed that this model is insufficient for generating a match with empirically observed spike correlations. However, by combining a phenomenological model of channel noise and our short term depression model, we generated short and long term correlations in auditory nerve spontaneous activity that qualitatively match empirical data.

Stochastic information transfer from cochlear implant electrodes to auditory nerve fibers

Cochlear implants, also called bionic ears, are implanted neural prostheses that can restore lost human hearing function by direct electrical stimulation of auditory nerve fibers. Previously, an information-theoretic framework for numerically estimating the optimal number of electrodes in cochlear implants has been devised. This approach relies on a model of stochastic action potential generation and a discrete memoryless channel model of the interface between the array of electrodes and the auditory nerve fibers. Using these models, the stochastic information transfer from cochlear implant electrodes to auditory nerve fibers is estimated from the mutual information between channel inputs (the locations of electrodes) and channel outputs (the set of electrode-activated nerve fibers). Here we describe a revised model of the channel output in the framework that avoids the side effects caused by an "ambiguity state" in the original model and also makes fewer assumptions about perceptual processing in the brain. A detailed comparison of how different assumptions on fibers and current spread modes impact on the information transfer in the original model and in the revised model is presented. We also mathematically derive an upper bound on the mutual information in the revised model, which becomes tighter as the number of electrodes increases. We found that the revised model leads to a significantly larger maximum mutual information and corresponding number of electrodes compared with the original model and conclude that the assumptions made in this part of the modeling framework are crucial to the model's overall utility.