Noise from voltage-gated ion channels may influence neuronal dynamics in the entorhinal cortex

Network motifs in cellular neurophysiology

Concepts from network science and graph theory, including the framework of network motifs, have been frequently applied in studying neuronal networks and other biological complex systems. Network-based approaches can also be used to study the functions of individual neurons, where cellular elements such as ion channels and membrane voltage are conceptualized as nodes within a network, and their interactions are denoted by edges. Network motifs in this context provide functional building blocks that help to illuminate the principles of cellular neurophysiology. In this review we build a case that network motifs operating within neurons provide tools for defining the functional architecture of single-neuron physiology and neuronal adaptations. We highlight the presence of such computational motifs in the cellular mechanisms underlying action potential generation, neuronal oscillations, dendritic integration, and neuronal plasticity. Future work applying the network motifs perspective may help to decipher the functional complexities of neurons and their adaptation during health and disease.

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.

The enigmatic HCN channels: A cellular neurophysiology perspective

What physiological role does a slow hyperpolarization-activated ion channel with mixed cation selectivity play in the fast world of neuronal action potentials that are driven by depolarization? That puzzling question has piqued the curiosity of physiology enthusiasts about the hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, which are widely expressed across the body and especially in neurons. In this review, we emphasize the need to assess HCN channels from the perspective of how they respond to time-varying signals, while also accounting for their interactions with other co-expressing channels and receptors. First, we illustrate how the unique structural and functional characteristics of HCN channels allow them to mediate a slow negative feedback loop in the neurons that they express in. We present the several physiological implications of this negative feedback loop to neuronal response characteristics including neuronal gain, voltage sag and rebound, temporal summation, membrane potential resonance, inductive phase lead, spike triggered average, and coincidence detection. Next, we argue that the overall impact of HCN channels on neuronal physiology critically relies on their interactions with other co-expressing channels and receptors. Interactions with other channels allow HCN channels to mediate intrinsic oscillations, earning them the "pacemaker channel" moniker, and to regulate spike frequency adaptation, plateau potentials, neurotransmitter release from presynaptic terminals, and spike initiation at the axonal initial segment. We also explore the impact of spatially non-homogeneous subcellular distributions of HCN channels in different neuronal subtypes and their interactions with other channels and receptors. Finally, we discuss how plasticity in HCN channels is widely prevalent and can mediate different encoding, homeostatic, and neuroprotective functions in a neuron. In summary, we argue that HCN channels form an important class of channels that mediate a diversity of neuronal functions owing to their unique gating kinetics that made them a puzzle in the first place.

A dynamic spike threshold with correlated noise predicts observed patterns of negative interval correlations in neuronal spike trains

Negative correlations in the sequential evolution of interspike intervals (ISIs) are a signature of memory in neuronal spike-trains. They provide coding benefits including firing-rate stabilization, improved detectability of weak sensory signals, and enhanced transmission of information by improving signal-to-noise ratio. Primary electrosensory afferent spike-trains in weakly electric fish fall into two categories based on the pattern of ISI correlations: non-bursting units have negative correlations which remain negative but decay to zero with increasing lags (Type I ISI correlations), and bursting units have oscillatory (alternating sign) correlation which damp to zero with increasing lags (Type II ISI correlations). Here, we predict and match observed ISI correlations in these afferents using a stochastic dynamic threshold model. We determine the ISI correlation function as a function of an arbitrary discrete noise correlation function [Formula: see text], where k is a multiple of the mean ISI. The function permits forward and inverse calculations of the correlation function. Both types of correlation functions can be generated by adding colored noise to the spike threshold with Type I correlations generated with slow noise and Type II correlations generated with fast noise. A first-order autoregressive (AR) process with a single parameter is sufficient to predict and accurately match both types of afferent ISI correlation functions, with the type being determined by the sign of the AR parameter. The predicted and experimentally observed correlations are in geometric progression. The theory predicts that the limiting sum of ISI correlations is [Formula: see text] yielding a perfect DC-block in the power spectrum of the spike train. Observed ISI correlations from afferents have a limiting sum that is slightly larger at [Formula: see text] ([Formula: see text]). We conclude that the underlying process for generating ISIs may be a simple combination of low-order AR and moving average processes and discuss the results from the perspective of optimal coding.

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.

The voltage and spiking responses of subthreshold resonant neurons to structured and fluctuating inputs: persistence and loss of resonance and variability

We systematically investigate the response of neurons to oscillatory currents and synaptic-like inputs and we extend our investigation to non-structured synaptic-like spiking inputs with more realistic distributions of presynaptic spike times. We use two types of chirp-like inputs consisting of (i) a sequence of cycles with discretely increasing frequencies over time, and (ii) a sequence having the same cycles arranged in an arbitrary order. We develop and use a number of frequency-dependent voltage response metrics to capture the different aspects of the voltage response, including the standard impedance (Z) and the peak-to-trough amplitude envelope ([Formula: see text]) profiles. We show that Z-resonant cells (cells that exhibit subthreshold resonance in response to sinusoidal inputs) also show [Formula: see text]-resonance in response to sinusoidal inputs, but generally do not (or do it very mildly) in response to square-wave and synaptic-like inputs. In the latter cases the resonant response using Z is not predictive of the preferred frequencies at which the neurons spike when the input amplitude is increased above subthreshold levels. We also show that responses to conductance-based synaptic-like inputs are attenuated as compared to the response to current-based synaptic-like inputs, thus providing an explanation to previous experimental results. These response patterns were strongly dependent on the intrinsic properties of the participating neurons, in particular whether the unperturbed Z-resonant cells had a stable node or a focus. In addition, we show that variability emerges in response to chirp-like inputs with arbitrarily ordered patterns where all signals (trials) in a given protocol have the same frequency content and the only source of uncertainty is the subset of all possible permutations of cycles chosen for a given protocol. This variability is the result of the multiple different ways in which the autonomous transient dynamics is activated across cycles in each signal (different cycle orderings) and across trials. We extend our results to include high-rate Poisson distributed current- and conductance-based synaptic inputs and compare them with similar results using additive Gaussian white noise. We show that the responses to both Poisson-distributed synaptic inputs are attenuated with respect to the responses to Gaussian white noise. For cells that exhibit oscillatory responses to Gaussian white noise (band-pass filters), the response to conductance-based synaptic inputs are low-pass filters, while the response to current-based synaptic inputs may remain band-pass filters, consistent with experimental findings. Our results shed light on the mechanisms of communication of oscillatory activity among neurons in a network via subthreshold oscillations and resonance and the generation of network resonance.

Mycorrhizal and rhizobial interactions influence model grassland plant community structure and productivity

Arbuscular mycorrhizal (AM) fungi and rhizobium are likely important drivers of plant coexistence and grassland productivity due to complementary roles in supplying limiting nutrients. However, the interactive effects of mycorrhizal and rhizobial associations on plant community productivity and competitive dynamics remain unclear. To address this, we conducted a greenhouse experiment to determine the influences of these key microbial functional groups on communities comprising three plant species by comparing plant communities grown with or without each symbiont. We also utilized N-fertilization and clipping treatments to explore potential shifts in mycorrhizal and rhizobial benefits across abiotic and biotic conditions. Our research suggests AM fungi and rhizobium co-inoculation was strongly facilitative for plant community productivity and legume (Medicago sativa) growth and nodulation. Plant competitiveness shifted in the presence of AM fungi and rhizobium, favoring M. sativa over a neighboring C₄ grass (Andropogon gerardii) and C₃ forb (Ratibida pinnata). This may be due to rhizobial symbiosis as well as the relatively greater mycorrhizal growth response of M. sativa, compared to the other model plants. Clipping and N-fertilization altered relative costs and benefits of both symbioses, presumably by altering host-plant nitrogen and carbon dynamics, leading to a relative decrease in mycorrhizal responsiveness and proportional biomass of M. sativa relative to the total biomass of the entire plant community, with a concomitant relative increase in A. gerardii and R. pinnata proportional biomass. Our results demonstrate a strong influence of both microbial symbioses on host-plant competitiveness and community dynamics across clipping and N-fertilization treatments, suggesting the symbiotic rhizosphere community is critical for legume establishment in grasslands.

Stochasticity Versus Determinacy in Neurobiology: From Ion Channels to the Question of the "Free Will"

If one accepts that decisions are made by the brain and that neuronal mechanisms obey deterministic physical laws, it is hard to deny what some brain researchers postulate, such as "We do not do what we want, but we want what we do" and "We should stop talking about freedom. Our actions are determined by physical laws." This point of view has been substantially supported by spectacular neurophysiological experiments demonstrating action-related brain activity (readiness potentials, blood oxygen level-dependent signals) occurring up to several seconds before an individual becomes aware of his/her decision to perform the action. This report aims to counter the deterministic argument for the absence of free will by using experimental data, supplemented by computer simulations, to demonstrate that biological systems, specifically brain functions, are built on principle randomness, which is introduced already at the lowest level of neuronal information processing, the opening and closing of ion channels. Switching between open and closed states follows physiological laws but also makes use of randomness, which is apparently introduced by Brownian motion - principally unavoidable under all life-compatible conditions. Ion-channel stochasticity, manifested as noise, function is not smoothed out toward higher functional levels but can even be amplified by appropriate adjustment of the system's non-linearities. Examples shall be given to illustrate how stochasticity can propagate from ion channels to single neuron action potentials to neuronal network dynamics to the interactions between different brain nuclei up to the control of autonomic functions. It is proposed that this intrinsic stochasticity helps to keep the brain in a flexible state to explore diverse alternatives as a prerequisite of free decision-making.

Fast and Accurate Langevin Simulations of Stochastic Hodgkin-Huxley Dynamics

Fox and Lu introduced a Langevin framework for discrete-time stochastic models of randomly gated ion channels such as the Hodgkin-Huxley (HH) system. They derived a Fokker-Planck equation with state-dependent diffusion tensor D and suggested a Langevin formulation with noise coefficient matrix S such that SS⊤=D. Subsequently, several authors introduced a variety of Langevin equations for the HH system. In this article, we present a natural 14-dimensional dynamics for the HH system in which each directed edge in the ion channel state transition graph acts as an independent noise source, leading to a 14 × 28 noise coefficient matrix S. We show that (1) the corresponding 14D system of ordinary differential equations is consistent with the classical 4D representation of the HH system; (2) the 14D representation leads to a noise coefficient matrix S that can be obtained cheaply on each time step, without requiring a matrix decomposition; (3) sample trajectories of the 14D representation are pathwise equivalent to trajectories of Fox and Lu's system, as well as trajectories of several existing Langevin models; (4) our 14D representation (and those equivalent to it) gives the most accurate interspike interval distribution, not only with respect to moments but under both the L1 and L∞ metric-space norms; and (5) the 14D representation gives an approximation to exact Markov chain simulations that are as fast and as efficient as all equivalent models. Our approach goes beyond existing models, in that it supports a stochastic shielding decomposition that dramatically simplifies S with minimal loss of accuracy under both voltage- and current-clamp conditions.

How to correctly quantify neuronal phase-response curves from noisy recordings

At the level of individual neurons, various coding properties can be inferred from the input-output relationship of a cell. For small inputs, this relation is captured by the phase-response curve (PRC), which measures the effect of a small perturbation on the timing of the subsequent spike. Experimentally, however, an accurate experimental estimation of PRCs is challenging. Despite elaborate measurement efforts, experimental PRC estimates often cannot be related to those from modeling studies. In particular, experimental PRCs rarely resemble the characteristic theoretical PRC expected close to spike initiation, which is indicative of the underlying spike-onset bifurcation. Here, we show for conductance-based model neurons that the correspondence between theoretical and measured phase-response curve is lost when the stimuli used for the estimation are too large. In this case, the derived phase-response curve is distorted beyond recognition and takes on a generic shape that reflects the measurement protocol and masks the spike-onset bifurcation. We discuss how to identify appropriate stimulus strengths for perturbation and noise-stimulation methods, which permit to estimate PRCs that reliably reflect the spike-onset bifurcation - a task that is particularly difficult if a lower bound for the stimulus amplitude is dictated by prominent intrinsic neuronal noise.

Control of clustered action potential firing in a mathematical model of entorhinal cortex stellate cells

•

An SDE model of entorhinal cortex (EC) stellate cells is proposed.

•

Experimentally observed action potential clustering is investigated in the model.

•

Clusters are generated by subcritical-Hopf/homoclinic type bursting.

•

Potential mechanisms underlying changes in EC dynamics in dementia are presented.

The entorhinal cortex is a crucial component of our memory and spatial navigation systems and is one of the first areas to be affected in dementias featuring tau pathology, such as Alzheimer’s disease and frontotemporal dementia. Electrophysiological recordings from principle cells of medial entorhinal cortex (layer II stellate cells, mEC-SCs) demonstrate a number of key identifying properties including subthreshold oscillations in the theta (4–12 Hz) range and clustered action potential firing. These single cell properties are correlated with network activity such as grid firing and coupling between theta and gamma rhythms, suggesting they are important for spatial memory. As such, experimental models of dementia have revealed disruption of organised dorsoventral gradients in clustered action potential firing.

To better understand the mechanisms underpinning these different dynamics, we study a conductance based model of mEC-SCs. We demonstrate that the model, driven by extrinsic noise, can capture quantitative differences in clustered action potential firing patterns recorded from experimental models of tau pathology and healthy animals. The differential equation formulation of our model allows us to perform numerical bifurcation analyses in order to uncover the dynamic mechanisms underlying these patterns. We show that clustered dynamics can be understood as subcritical Hopf/homoclinic bursting in a fast-slow system where the slow sub-system is governed by activation of the persistent sodium current and inactivation of the slow A-type potassium current. In the full system, we demonstrate that clustered firing arises via flip bifurcations as conductance parameters are varied. Our model analyses confirm the experimentally suggested hypothesis that the breakdown of clustered dynamics in disease occurs via increases in AHP conductance.

Nitric oxide augments single persistent Na+ channel currents via the cGMP/PKG signaling pathway in Kenyon cells isolated from cricket mushroom bodies

The nitric oxide (NO)/cyclic GMP signaling pathway has been suggested to be important in the formation of olfactory memory in insects. However, the molecular targets of the NO signaling cascade in the central neurons associated with olfactory learning and memory have not been fully analyzed. In this study, we investigated the effects of NO donors on single voltage-dependent Na⁺ channels in intrinsic neurons, called Kenyon cells, in the mushroom bodies in the brain of the cricket. Step depolarization on cell-attached patch membranes induces single-channel currents with fast-activating and -inactivating brief openings at the beginning of the voltage steps followed by more persistently recurring brief openings all along the 150-ms pulses. Application of the NO donor S-nitrosoglutathione (GSNO) increased the number of channel openings of both types of single Na⁺ channels. This excitatory effect of GSNO on the activity of these Na⁺ channels was diminished by KT5823, an inhibitor of protein kinase G (PKG), indicating an involvement of PKG in the downstream pathway of NO. Application of KT5823 alone decreased the activity of the persistent Na⁺ channels without significant effects on the fast-inactivating Na⁺ channels. The membrane-permeable cGMP analog 8Br-cGMP increased the number of channel openings of both types of single Na⁺ channels, similar to the action of NO. Taken together, these results indicate that NO acts as a critical modulator of both fast-inactivating and persistent Na⁺ channels and that persistent Na⁺ channels are constantly upregulated by the endogenous cGMP/PKG signaling cascade. NEW & NOTEWORTHY This study clarified that nitric oxide (NO) increases the activity of both fast-inactivating and persistent Na⁺ channels via the cGMP/PKG signaling cascade in cricket Kenyon cells. The persistent Na⁺ channels are also found to be upregulated constantly by endogenous cGMP/PKG signaling. On the basis of the present results and the results of previous studies, we propose a hypothetical model explaining NO production and NO-dependent memory formation in cricket large Kenyon cells.

Degeneracy in the robust expression of spectral selectivity, subthreshold oscillations, and intrinsic excitability of entorhinal stellate cells

Biological heterogeneities are ubiquitous and play critical roles in the emergence of physiology at multiple scales. Although neurons in layer II (LII) of the medial entorhinal cortex (MEC) express heterogeneities in channel properties, the impact of such heterogeneities on the robustness of their cellular-scale physiology has not been assessed. Here, we performed a 55-parameter stochastic search spanning nine voltage- or calcium-activated channels to assess the impact of channel heterogeneities on the concomitant emergence of 10 in vitro electrophysiological characteristics of LII stellate cells (SCs). We generated 150,000 models and found a heterogeneous subpopulation of 449 valid models to robustly match all electrophysiological signatures. We employed this heterogeneous population to demonstrate the emergence of cellular-scale degeneracy in SCs, whereby disparate parametric combinations expressing weak pairwise correlations resulted in similar models. We then assessed the impact of virtually knocking out each channel from all valid models and demonstrate that the mapping between channels and measurements was many-to-many, a critical requirement for the expression of degeneracy. Finally, we quantitatively predict that the spike-triggered average of SCs should be endowed with theta-frequency spectral selectivity and coincidence detection capabilities in the fast gamma-band. We postulate this fast gamma-band coincidence detection as an instance of cellular-scale-efficient coding, whereby SC response characteristics match the dominant oscillatory signals in LII MEC. The heterogeneous population of valid SC models built here unveils the robust emergence of cellular-scale physiology despite significant channel heterogeneities, and forms an efficacious substrate for evaluating the impact of biological heterogeneities on entorhinal network function.

NEW & NOTEWORTHY We assessed the impact of heterogeneities in channel properties on the robustness of cellular-scale physiology of medial entorhinal cortical stellate neurons. We demonstrate that neuronal models with disparate channel combinations were endowed with similar physiological characteristics, as a consequence of the many-to-many mapping between channel properties and the physiological characteristics that they modulate. We predict that the spike-triggered average of stellate cells should be endowed with theta-frequency spectral selectivity and fast gamma-band coincidence detection capabilities.

Effects of channel blocking on information transmission and energy efficiency in squid giant axons

Action potentials are the information carriers of neural systems. The generation of action potentials involves the cooperative opening and closing of sodium and potassium channels. This process is metabolically expensive because the ions flowing through open channels need to be restored to maintain concentration gradients of these ions. Toxins like tetraethylammonium can block working ion channels, thus affecting the function and energy cost of neurons. In this paper, by computer simulation of the Hodgkin-Huxley neuron model, we studied the effects of channel blocking with toxins on the information transmission and energy efficiency in squid giant axons. We found that gradually blocking sodium channels will sequentially maximize the information transmission and energy efficiency of the axons, whereas moderate blocking of potassium channels will have little impact on the information transmission and will decrease the energy efficiency. Heavy blocking of potassium channels will cause self-sustained oscillation of membrane potentials. Simultaneously blocking sodium and potassium channels with the same ratio increases both information transmission and energy efficiency. Our results are in line with previous studies suggesting that information processing capacity and energy efficiency can be maximized by regulating the number of active ion channels, and this indicates a viable avenue for future experimentation.

Stochastic resonance improves vision in the severely impaired

We verified whether a stochastic resonance paradigm (SR), with random interference (“noise”) added in optimal amounts, improves the detection of sub-threshold visual information by subjects with retinal disorder and impaired vision as it does in the normally sighted. Six levels of dynamic, zero-mean Gaussian noise were added to each pixel of images (13 contrast levels) in which alphabet characters were displayed against a uniform gray background. Images were presented with contrast below the subjective threshold to 14 visually impaired subjects (age: 22–53 yrs.). The fraction of recognized letters varied between 0 and 0.3 at baseline and increased in all subjects when noise was added in optimal amounts; peak recognition ranged between 0.2 and 0.8 at noise sigmas between 6 and 30 grey scale values (GSV) and decreased in all subjects at noise levels with sigma above 30 GSV. The results replicate in the visually impaired the facilitation of visual information processing with images presented in SR paradigms that has been documented in sighted subjects. The effect was obtained with low-level image manipulation and application appears readily possible: it would enhance the efficiency of today vision-improving aids and help in the development of the visual prostheses hopefully available in the future.

BIOPHYSICAL PROPERTIES OF SUBTHRESHOLD RESONANCE OSCILLATIONS AND SUBTHRESHOLD MEMBRANE OSCILLATIONS IN NEURONS

Subthreshold-level activities in neurons play a crucial role in neuronal oscillations. These small-amplitude oscillations have been suggested to be involved in synaptic plasticity and in determining the frequency of network oscillations. Subthreshold membrane oscillations (STOs) and subthreshold resonance oscillations (SROs) are the main constituents of subthreshold-level activities in neurons. In this study, a general theoretical framework for analyzing the mechanisms underlying STOs and SROs in neurons is presented. Results showed that the resting membrane potential and the hyperpolarization-activated potassium channel (h-channel) affect the subthreshold-level activities in stellate cells. The contribution of h-channel on resonance is attributed to its large time constant, which produces the time lag between I_h and the membrane potential. Conversely, the persistent sodium channels (Nap-channels) only play an amplifying role in these neurons.

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

Noise Enhances Action Potential Generation in Mouse Sensory Neurons via Stochastic Resonance

Noise can enhance perception of tactile and proprioceptive stimuli by stochastic resonance processes. However, the mechanisms underlying this general phenomenon remain to be characterized. Here we studied how externally applied noise influences action potential firing in mouse primary sensory neurons of dorsal root ganglia, modelling a basic process in sensory perception. Since noisy mechanical stimuli may cause stochastic fluctuations in receptor potential, we examined the effects of sub-threshold depolarizing current steps with superimposed random fluctuations. We performed whole cell patch clamp recordings in cultured neurons of mouse dorsal root ganglia. Noise was added either before and during the step, or during the depolarizing step only, to focus onto the specific effects of external noise on action potential generation. In both cases, step + noise stimuli triggered significantly more action potentials than steps alone. The normalized power norm had a clear peak at intermediate noise levels, demonstrating that the phenomenon is driven by stochastic resonance. Spikes evoked in step + noise trials occur earlier and show faster rise time as compared to the occasional ones elicited by steps alone. These data suggest that external noise enhances, via stochastic resonance, the recruitment of transient voltage-gated Na channels, responsible for action potential firing in response to rapid step-wise depolarizing currents.

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.

Cooperation and competition of gamma oscillation mechanisms

Oscillations of neuronal activity in different frequency ranges are thought to reflect important aspects of cortical network dynamics. Here we investigate how various mechanisms that contribute to oscillations in neuronal networks may interact. We focus on networks with inhibitory, excitatory, and electrical synapses, where the subnetwork of inhibitory interneurons alone can generate interneuron gamma (ING) oscillations and the interactions between interneurons and pyramidal cells allow for pyramidal-interneuron gamma (PING) oscillations. What type of oscillation will such a network generate? We find that ING and PING oscillations compete: The mechanism generating the higher oscillation frequency "wins"; it determines the frequency of the network oscillation and suppresses the other mechanism. For type I interneurons, the network oscillation frequency is equal to or slightly above the higher of the ING and PING frequencies in corresponding reduced networks that can generate only either of them; if the interneurons belong to the type II class, it is in between. In contrast to ING and PING, oscillations mediated by gap junctions and oscillations mediated by inhibitory synapses may cooperate or compete, depending on the type (I or II) of interneurons and the strengths of the electrical and chemical synapses. We support our computer simulations by a theoretical model that allows a full theoretical analysis of the main results. Our study suggests experimental approaches to deciding to what extent oscillatory activity in networks of interacting excitatory and inhibitory neurons is dominated by ING or PING oscillations and of which class the participating interneurons are.

Predicting the onset of period-doubling bifurcations in noisy cardiac systems

Biological, physical, and social systems often display qualitative changes in dynamics. Developing early warning signals to predict the onset of these transitions is an important goal. The current work is motivated by transitions of cardiac rhythms, where the appearance of alternating features in the timing of cardiac events is often a precursor to the initiation of serious cardiac arrhythmias. We treat embryonic chick cardiac cells with a potassium channel blocker, which leads to the initiation of alternating rhythms. We associate this transition with a mathematical instability, called a period-doubling bifurcation, in a model of the cardiac cells. Period-doubling bifurcations have been linked to the onset of abnormal alternating cardiac rhythms, which have been implicated in cardiac arrhythmias such as T-wave alternans and various tachycardias. Theory predicts that in the neighborhood of the transition, the system's dynamics slow down, leading to noise amplification and the manifestation of oscillations in the autocorrelation function. Examining the aggregates' interbeat intervals, we observe the oscillations in the autocorrelation function and noise amplification preceding the bifurcation. We analyze plots--termed return maps--that relate the current interbeat interval with the following interbeat interval. Based on these plots, we develop a quantitative measure using the slope of the return map to assess how close the system is to the bifurcation. Furthermore, the slope of the return map and the lag-1 autocorrelation coefficient are equal. Our results suggest that the slope and the lag-1 autocorrelation coefficient represent quantitative measures to predict the onset of abnormal alternating cardiac rhythms.

Quasicycles in the stochastic hybrid Morris-Lecar neural model

Intrinsic noise arising from the stochastic opening and closing of voltage-gated ion channels has been shown experimentally and mathematically to have important effects on a neuron's function. Study of classical neuron models with stochastic ion channels is becoming increasingly important, especially in understanding a cell's ability to produce subthreshold oscillations and to respond to weak periodic stimuli. While it is known that stochastic models can produce oscillations (quasicycles) in parameter regimes where the corresponding deterministic model has only a stable fixed point, little analytical work has been done to explore these connections within the context of channel noise. Using a stochastic hybrid Morris-Lecar (ML) model, we combine a system-size expansion in K(+) and a quasi-steady-state (QSS) approximation in persistent Na(+) in order to derive an effective Langevin equation that preserves the low-dimensional (planar) structure of the underlying deterministic ML model. (The QSS analysis exploits the fact that persistent Na(+) channels are fast.) By calculating the corresponding power spectrum, we determine analytically how noise significantly extends the parameter regime in which subthreshold oscillations occur.

Real-time Electrophysiology: Using Closed-loop Protocols to Probe Neuronal Dynamics and Beyond

Experimental neuroscience is witnessing an increased interest in the development and application of novel and often complex, closed-loop protocols, where the stimulus applied depends in real-time on the response of the system. Recent applications range from the implementation of virtual reality systems for studying motor responses both in mice and in zebrafish, to control of seizures following cortical stroke using optogenetics. A key advantage of closed-loop techniques resides in the capability of probing higher dimensional properties that are not directly accessible or that depend on multiple variables, such as neuronal excitability and reliability, while at the same time maximizing the experimental throughput. In this contribution and in the context of cellular electrophysiology, we describe how to apply a variety of closed-loop protocols to the study of the response properties of pyramidal cortical neurons, recorded intracellularly with the patch clamp technique in acute brain slices from the somatosensory cortex of juvenile rats. As no commercially available or open source software provides all the features required for efficiently performing the experiments described here, a new software toolbox called LCG was developed, whose modular structure maximizes reuse of computer code and facilitates the implementation of novel experimental paradigms. Stimulation waveforms are specified using a compact meta-description and full experimental protocols are described in text-based configuration files. Additionally, LCG has a command-line interface that is suited for repetition of trials and automation of experimental protocols.

Non-linear Membrane Properties in Entorhinal Cortical Stellate Cells Reduce Modulation of Input-Output Responses by Voltage Fluctuations

The presence of voltage fluctuations arising from synaptic activity is a critical component in models of gain control, neuronal output gating, and spike rate coding. The degree to which individual neuronal input-output functions are modulated by voltage fluctuations, however, is not well established across different cortical areas. Additionally, the extent and mechanisms of input-output modulation through fluctuations have been explored largely in simplified models of spike generation, and with limited consideration for the role of non-linear and voltage-dependent membrane properties. To address these issues, we studied fluctuation-based modulation of input-output responses in medial entorhinal cortical (MEC) stellate cells of rats, which express strong sub-threshold non-linear membrane properties. Using in vitro recordings, dynamic clamp and modeling, we show that the modulation of input-output responses by random voltage fluctuations in stellate cells is significantly limited. In stellate cells, a voltage-dependent increase in membrane resistance at sub-threshold voltages mediated by Na⁺ conductance activation limits the ability of fluctuations to elicit spikes. Similarly, in exponential leaky integrate-and-fire models using a shallow voltage-dependence for the exponential term that matches stellate cell membrane properties, a low degree of fluctuation-based modulation of input-output responses can be attained. These results demonstrate that fluctuation-based modulation of input-output responses is not a universal feature of neurons and can be significantly limited by subthreshold voltage-gated conductances.

The membrane voltage of neurons in vivo is dominated by noisy “background” fluctuations generated by network-based synaptic activity from nearby cells. It has been speculated that membrane voltage fluctuations in neurons play an important role in scaling the relationship between input amplitude and spike rate response. For this to be true, neuronal spike input-output behavior must be sensitive to physiological membrane voltage fluctuations. Using a combination of single cell recordings and modeling, we investigated the mechanisms through which voltage fluctuations modulate neuronal input-output responses. We find that neurons that express an increase in membrane input resistance with depolarization show low levels of noise-mediated modulation of input-output responses due, in part, to voltage trajectories that suppress the likelihood of generating a spike in response to random current input fluctuations. Hence, non-linear membrane properties arising from certain types of voltage-gated conductances limit noise-based modulation of neuronal input-output responses.

Low frequency electric current noise in glioma cell populations

Measuring the electrical activity of large and defined populations of cells is currently a major technical challenge to electrophysiology, especially in the picoampere-range. For this purpose, we developed and applied a bidirectional transducer based on a chip with interdigitated gold electrodes to record the electrical response of cultured glioma cells. Recent research determined that also non-neural brain glia cells are electrically active and excitable. Their transformed counterparts, e.g. glioma cells, were suggested to partially retain these electric features. Such electrophysiological studies however are usually performed on individual cells and are limited in their predictive power for the overall electrical activity of the multicellular tumour bulk. Our extremely low-noise measuring system allowed us to detect not only prominent electrical bursts of neuronal cells but also minute, yet constantly occurring and functional, membrane capacitive current oscillations across large populations of C6 glioma cells, which we termed electric current noise. At the same time, tumour cells of non-brain origin (HeLa) proved to be electrically quiescent in comparison. Finally, we determined that the glioma cell activity is primarily caused by the opening of voltage-gated Na⁺ and K⁺ ion channels and can be efficiently abolished using specific pharmacological inhibitors. Thus, we offer here a unique approach for studying electrophysiological properties of large cancer cell populations as an in vitro reference for tumour bulks in vivo.

Cell-attached single-channel recordings in intact prefrontal cortex pyramidal neurons reveal compartmentalized D1/D5 receptor modulation of the persistent sodium current

The persistent Na⁺ current (I_Nap) is believed to be an important target of dopamine modulation in prefrontal cortex (PFC) neurons. While past studies have tested the effects of dopamine on I_Nap, the results have been contradictory largely because of difficulties in measuring I_Nap using somatic whole-cell recordings. To circumvent these confounds we used the cell-attached patch-clamp technique to record single Na⁺ channels from the soma, proximal dendrite (PD) or proximal axon (PA) of intact prefrontal layer V pyramidal neurons. Under baseline conditions, numerous well resolved Na⁺ channel openings were recorded that exhibited an extrapolated reversal potential of 73 mV, a slope conductance of 14–19 pS and were blocked by tetrodotoxin (TTX). While similar in most respects, the propensity to exhibit prolonged bursts lasting >40 ms was many fold greater in the axon than the soma or dendrite. Bath application of the D1/D5 receptor agonist SKF81297 shifted the ensemble current activation curve leftward and increased the number of late events recorded from the PD but not the soma or PA. However, the greatest effect was on prolonged bursting where the D1/D5 receptor agonist increased their occurrence 3 fold in the PD and nearly 7 fold in the soma, but not at all in the PA. As a result, D1/D5 receptor activation equalized the probability of prolonged burst occurrence across the proximal axosomatodendritic region. Therefore, D1/D5 receptor modulation appears to be targeted mainly to Na⁺ channels in the PD/soma and not the PA. By circumventing the pitfalls of previous attempts to study the D1/D5 receptor modulation of I_Nap, we demonstrate conclusively that D1/D5 receptor activation can increase the I_Nap generated proximally, however questions still remain as to how D1/D5 receptor modulates Na⁺ currents in the more distal initial segment where most of the I_Nap is normally generated.

Asymptotic phase for stochastic oscillators

Oscillations and noise are ubiquitous in physical and biological systems. When oscillations arise from a deterministic limit cycle, entrainment and synchronization may be analyzed in terms of the asymptotic phase function. In the presence of noise, the asymptotic phase is no longer well defined. We introduce a new definition of asymptotic phase in terms of the slowest decaying modes of the Kolmogorov backward operator. Our stochastic asymptotic phase is well defined for noisy oscillators, even when the oscillations are noise dependent. It reduces to the classical asymptotic phase in the limit of vanishing noise. The phase can be obtained either by solving an eigenvalue problem, or by empirical observation of an oscillating density's approach to its steady state.

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.

Volterra dendritic stimulus processors and biophysical spike generators with intrinsic noise sources

We consider a class of neural circuit models with internal noise sources arising in sensory systems. The basic neuron model in these circuits consists of a dendritic stimulus processor (DSP) cascaded with a biophysical spike generator (BSG). The dendritic stimulus processor is modeled as a set of nonlinear operators that are assumed to have a Volterra series representation. Biophysical point neuron models, such as the Hodgkin-Huxley neuron, are used to model the spike generator. We address the question of how intrinsic noise sources affect the precision in encoding and decoding of sensory stimuli and the functional identification of its sensory circuits. We investigate two intrinsic noise sources arising (i) in the active dendritic trees underlying the DSPs, and (ii) in the ion channels of the BSGs. Noise in dendritic stimulus processing arises from a combined effect of variability in synaptic transmission and dendritic interactions. Channel noise arises in the BSGs due to the fluctuation of the number of the active ion channels. Using a stochastic differential equations formalism we show that encoding with a neuron model consisting of a nonlinear DSP cascaded with a BSG with intrinsic noise sources can be treated as generalized sampling with noisy measurements. For single-input multi-output neural circuit models with feedforward, feedback and cross-feedback DSPs cascaded with BSGs we theoretically analyze the effect of noise sources on stimulus decoding. Building on a key duality property, the effect of noise parameters on the precision of the functional identification of the complete neural circuit with DSP/BSG neuron models is given. We demonstrate through extensive simulations the effects of noise on encoding stimuli with circuits that include neuron models that are akin to those commonly seen in sensory systems, e.g., complex cells in V1.

Systematic analysis of the contributions of stochastic voltage gated channels to neuronal noise

Electrical signaling in neurons is mediated by the opening and closing of large numbers of individual ion channels. The ion channels' state transitions are stochastic and introduce fluctuations in the macroscopic current through ion channel populations. This creates an unavoidable source of intrinsic electrical noise for the neuron, leading to fluctuations in the membrane potential and spontaneous spikes. While this effect is well known, the impact of channel noise on single neuron dynamics remains poorly understood. Most results are based on numerical simulations. There is no agreement, even in theoretical studies, on which ion channel type is the dominant noise source, nor how inclusion of additional ion channel types affects voltage noise. Here we describe a framework to calculate voltage noise directly from an arbitrary set of ion channel models, and discuss how this can be use to estimate spontaneous spike rates.

Frequency-dependent signal processing in apical dendrites of hippocampal CA1 pyramidal cells

Depending on an animal's behavioral state, hippocampal CA1 pyramidal cells receive distinct patterns of excitatory and inhibitory synaptic inputs. The time-dependent changes in the frequencies of these inputs and the nonuniform distribution of voltage-gated channels lead to dynamic fluctuations in membrane conductance. In this study, using a whole-cell patch-clamp method, we attempted to record and analyze the frequency dependencies of membrane responsiveness in Wistar rat hippocampal CA1 pyramidal cells following noise current injection directly into dendrites and somata under pharmacological blockade of all synaptic inputs. To estimate the frequency-dependent properties of membrane potential, membrane impedance was determined from the voltage response divided by the input current in the frequency domain. The cell membrane of most neurons showed low-pass filtering properties in all regions. In particular, the properties were strongly expressed in the somata or proximal dendrites. Moreover, the data revealed nonuniform distribution of dendritic impedance, which was high in the intermediate segment of the apical dendritic shaft (∼220-260μm from the soma). The low-pass filtering properties in the apical dendrites were more enhanced by membrane depolarization than those in the somata. Coherence spectral analysis revealed high coherence between the input signal and the output voltage response in the theta-gamma frequency range, and large lags emerged in the distal dendrites in the gamma frequency range. Our results suggest that apical dendrites of hippocampal CA1 pyramidal cells integrate synaptic inputs according to the frequency components of the input signal along the dendritic segments receiving the inputs.

Membrane potential-dependent integration of synaptic inputs in entorhinal stellate neurons

Stellate cells (SCs) of the medial entorhinal cortex exhibit robust spontaneous membrane-potential oscillations (MPOs) in the theta (4-12 Hz) frequency band as well as theta-frequency resonance in their membrane impedance spectra. Past experimental and modeling work suggests that these features may contribute to the phase-locking of SCs to the entorhinal theta rhythm and may be important for forming the hexagonally tiled grid cell place fields exhibited by these neurons in vivo. Among the major biophysical mechanisms contributing to MPOs is a population of persistent (non-inactivating or slowly inactivating) sodium channels. The resulting persistent sodium conductance (GNaP ) gives rise to an apparent increase in input resistance as the cell approaches threshold. In this study, we used dynamic clamp to test the hypothesis that this increased input resistance gives rise to voltage-dependent, and thus MPO phase-dependent, changes in the amplitude of excitatory and inhibitory post-synaptic potential (PSP) amplitudes. We find that PSP amplitude depends on membrane potential, exhibiting a 5-10% increase in amplitude per mV depolarization. The effect is larger than-and sums quasi-linearly with-the effect of the synaptic driving force, V - Esyn . Given that input-driven MPOs 10 mV in amplitude are commonly observed in MEC stellate cells in vivo, this voltage- and phase-dependent synaptic gain is large enough to modulate PSP amplitude by over 50% during theta-frequency MPOs. Phase-dependent synaptic gain may therefore impact the phase locking and phase precession of grid cells in vivo to ongoing network oscillations. © 2014 Wiley Periodicals, Inc.

Command-line cellular electrophysiology for conventional and real-time closed-loop experiments

BACKGROUND

Current software tools for electrophysiological experiments are limited in flexibility and rarely offer adequate support for advanced techniques such as dynamic clamp and hybrid experiments, which are therefore limited to laboratories with a significant expertise in neuroinformatics.

NEW METHOD

We have developed lcg, a software suite based on a command-line interface (CLI) that allows performing both standard and advanced electrophysiological experiments. Stimulation protocols for classical voltage and current clamp experiments are defined by a concise and flexible meta description that allows representing complex waveforms as a piece-wise parametric decomposition of elementary sub-waveforms, abstracting the stimulation hardware. To perform complex experiments lcg provides a set of elementary building blocks that can be interconnected to yield a large variety of experimental paradigms.

RESULTS

We present various cellular electrophysiological experiments in which lcg has been employed, ranging from the automated application of current clamp protocols for characterizing basic electrophysiological properties of neurons, to dynamic clamp, response clamp, and hybrid experiments. We finally show how the scripting capabilities behind a CLI are suited for integrating experimental trials into complex workflows, where actual experiment, online data analysis and computational modeling seamlessly integrate.

COMPARISON WITH EXISTING METHODS

We compare lcg with two open source toolboxes, RTXI and RELACS.

CONCLUSIONS

We believe that lcg will greatly contribute to the standardization and reproducibility of both simple and complex experiments. Additionally, on the long run the increased efficiency due to a CLI will prove a great benefit for the experimental community.

A novel model incorporating two variability sources for describing motor evoked potentials

OBJECTIVE

Motor evoked potentials (MEPs) play a pivotal role in transcranial magnetic stimulation (TMS), e.g., for determining the motor threshold and probing cortical excitability. Sampled across the range of stimulation strengths, MEPs outline an input-output (IO) curve, which is often used to characterize the corticospinal tract. More detailed understanding of the signal generation and variability of MEPs would provide insight into the underlying physiology and aid correct statistical treatment of MEP data.

METHODS

A novel regression model is tested using measured IO data of twelve subjects. The model splits MEP variability into two independent contributions, acting on both sides of a strong sigmoidal nonlinearity that represents neural recruitment. Traditional sigmoidal regression with a single variability source after the nonlinearity is used for comparison.

RESULTS

The distribution of MEP amplitudes varied across different stimulation strengths, violating statistical assumptions in traditional regression models. In contrast to the conventional regression model, the dual variability source model better described the IO characteristics including phenomena such as changing distribution spread and skewness along the IO curve.

CONCLUSIONS

MEP variability is best described by two sources that most likely separate variability in the initial excitation process from effects occurring later on. The new model enables more accurate and sensitive estimation of the IO curve characteristics, enhancing its power as a detection tool, and may apply to other brain stimulation modalities. Furthermore, it extracts new information from the IO data concerning the neural variability-information that has previously been treated as noise.

Optimal size of stochastic Hodgkin-Huxley neuronal systems for maximal energy efficiency in coding pulse signals

The generation and conduction of action potentials (APs) represents a fundamental means of communication in the nervous system and is a metabolically expensive process. In this paper, we investigate the energy efficiency of neural systems in transferring pulse signals with APs. By analytically solving a bistable neuron model that mimics the AP generation with a particle crossing the barrier of a double well, we find the optimal number of ion channels that maximizes the energy efficiency of a neuron. We also investigate the energy efficiency of a neuron population in which the input pulse signals are represented with synchronized spikes and read out with a downstream coincidence detector neuron. We find an optimal number of neurons in neuron population, as well as the number of ion channels in each neuron that maximizes the energy efficiency. The energy efficiency also depends on the characters of the input signals, e.g., the pulse strength and the interpulse intervals. These results are confirmed by computer simulation of the stochastic Hodgkin-Huxley model with a detailed description of the ion channel random gating. We argue that the tradeoff between signal transmission reliability and energy cost may influence the size of the neural systems when energy use is constrained.

Consequences of converting graded to action potentials upon neural information coding and energy efficiency

Information is encoded in neural circuits using both graded and action potentials, converting between them within single neurons and successive processing layers. This conversion is accompanied by information loss and a drop in energy efficiency. We investigate the biophysical causes of this loss of information and efficiency by comparing spiking neuron models, containing stochastic voltage-gated Na(+) and K(+) channels, with generator potential and graded potential models lacking voltage-gated Na(+) channels. We identify three causes of information loss in the generator potential that are the by-product of action potential generation: (1) the voltage-gated Na(+) channels necessary for action potential generation increase intrinsic noise and (2) introduce non-linearities, and (3) the finite duration of the action potential creates a 'footprint' in the generator potential that obscures incoming signals. These three processes reduce information rates by ∼50% in generator potentials, to ∼3 times that of spike trains. Both generator potentials and graded potentials consume almost an order of magnitude less energy per second than spike trains. Because of the lower information rates of generator potentials they are substantially less energy efficient than graded potentials. However, both are an order of magnitude more efficient than spike trains due to the higher energy costs and low information content of spikes, emphasizing that there is a two-fold cost of converting analogue to digital; information loss and cost inflation.

How to build a grid cell

Neurons in the medial entorhinal cortex fire action potentials at regular spatial intervals, creating a striking grid-like pattern of spike rates spanning the whole environment of a navigating animal. This remarkable spatial code may represent a neural map for path integration. Recent advances using patch-clamp recordings from entorhinal cortex neurons in vitro and in vivo have revealed how the microcircuitry in the medial entorhinal cortex may contribute to grid cell firing patterns, and how grid cells may transform synaptic inputs into spike output during firing field crossings. These new findings provide key insights into the ingredients necessary to build a grid cell.

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.

Active decorrelation in the basal ganglia

The cytoarchitecturally-homogeneous appearance of the globus pallidus, subthalamic nucleus and substantia nigra has long been said to imply a high degree of afferent convergence and sharing of inputs by nearby neurons. Moreover, axon collaterals of neurons in the external segment of the globus pallidus and the substantia nigra pars reticulata arborize locally and make inhibitory synapses on other cells of the same type. These features suggest that the connectivity of the basal ganglia may impose spike-time correlations among the cells, and it has been puzzling that experimental studies have failed to demonstrate such correlations. One possible solution arises from studies of firing patterns in basal ganglia cells, which reveal that they are nearly all pacemaker cells. Their high rate of firing does not depend on synaptic excitation, but they fire irregularly because a dense barrage of synaptic inputs normally perturbs the timing of their autonomous activity. Theoretical and computational studies show that the responses of repetitively-firing neurons to shared input or mutual synaptic coupling often defy classical intuitions about temporal synaptic integration. The patterns of spike-timing among such neurons depend on the ionic mechanism of pacemaking, the level of background uncorrelated cellular and synaptic noise, and the firing rates of the neurons, as well as the properties of their synaptic connections. Application of these concepts to the basal ganglia circuitry suggests that the connectivity and physiology of these nuclei may be configured to prevent the establishment of permanent spike-timing relationships between neurons. The development of highly synchronous oscillatory patterns of activity in Parkinson's disease may result from the loss of pacemaking by some basal ganglia neurons, and accompanying breakdown of the mechanisms responsible for active decorrelation.

Comparison of models for IP3 receptor kinetics using stochastic simulations

Inositol 1,4,5-trisphosphate receptor (IP₃R) is a ubiquitous intracellular calcium (Ca²⁺) channel which has a major role in controlling Ca²⁺ levels in neurons. A variety of computational models have been developed to describe the kinetic function of IP₃R under different conditions. In the field of computational neuroscience, it is of great interest to apply the existing models of IP₃R when modeling local Ca²⁺ transients in dendrites or overall Ca²⁺ dynamics in large neuronal models. The goal of this study was to evaluate existing IP₃R models, based on electrophysiological data. This was done in order to be able to suggest suitable models for neuronal modeling. Altogether four models (Othmer and Tang, 1993; Dawson etal., 2003; Fraiman and Dawson, 2004; Doi etal., 2005) were selected for a more detailed comparison. The selection was based on the computational efficiency of the models and the type of experimental data that was used in developing the model. The kinetics of all four models were simulated by stochastic means, using the simulation software STEPS, which implements the Gillespie stochastic simulation algorithm. The results show major differences in the statistical properties of model functionality. Of the four compared models, the one by Fraiman and Dawson (2004) proved most satisfactory in producing the specific features of experimental findings reported in literature. To our knowledge, the present study is the first detailed evaluation of IP₃R models using stochastic simulation methods, thus providing an important setting for constructing a new, realistic model of IP₃R channel kinetics for compartmental modeling of neuronal functions. We conclude that the kinetics of IP₃R with different concentrations of Ca²⁺ and IP₃ should be more carefully addressed when new models for IP₃R are developed.

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.

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.

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.

Transient sodium current at subthreshold voltages: activation by EPSP waveforms

Tetrodotoxin (TTX)-sensitive sodium channels carry large transient currents during action potentials and also "persistent" sodium current, a noninactivating TTX-sensitive current present at subthreshold voltages. We examined gating of subthreshold sodium current in dissociated cerebellar Purkinje neurons and hippocampal CA1 neurons, studied at 37°C with near-physiological ionic conditions. Unexpectedly, in both cell types small voltage steps at subthreshold voltages activated a substantial component of transient sodium current as well as persistent current. Subthreshold EPSP-like waveforms also activated a large component of transient sodium current, but IPSP-like waveforms engaged primarily persistent sodium current with only a small additional transient component. Activation of transient as well as persistent sodium current at subthreshold voltages produces amplification of EPSPs that is sensitive to the rate of depolarization and can help account for the dependence of spike threshold on depolarization rate, as previously observed in vivo.

Contribution of near-threshold currents to intrinsic oscillatory activity in rat medial entorhinal cortex layer II stellate cells

The temporal lobe is well known for its oscillatory activity associated with exploration, navigation, and learning. Intrinsic membrane potential oscillations (MPOs) and resonance of stellate cells (SCs) in layer II of the entorhinal cortex are thought to contribute to network oscillations and thereby to the encoding of spatial information. Generation of both MPOs and resonance relies on the expression of specific voltage-dependent ion currents such as the hyperpolarization-activated cation current (I(H)), the persistent sodium current (I(NaP)), and the noninactivating muscarine-modulated potassium current (I(M)). However, the differential contributions of these currents remain a matter of debate. We therefore examined how they modify neuronal excitability near threshold and generation of near-threshold MPOs and resonance in vitro. We found that resonance mainly relied on I(H) and was reduced by I(H) blockers and modulated by cAMP and an I(M) enhancer but that neither of the currents exhibited full control over MPOs in these cells. As previously reported, I(H) controlled a theta-frequency component of MPOs such that blockade of I(H) resulted in fewer regular oscillations that retained low-frequency components and high peak amplitude. However, pharmacological inhibition and augmentation of I(M) also affected MPO frequencies and amplitudes. In contrast to other cell types, inhibition of I(NaP) did not result in suppression of MPOs but only in a moderation of their properties. We reproduced the experimentally observed effects in a single-compartment stochastic model of SCs, providing further insight into the interactions between different ionic conductances.

Spontaneous excitation patterns computed for axons with injury-like impairments of sodium channels and Na/K pumps

In injured neurons, “leaky” voltage-gated sodium channels (Nav) underlie dysfunctional excitability that ranges from spontaneous subthreshold oscillations (STO), to ectopic (sometimes paroxysmal) excitation, to depolarizing block. In recombinant systems, mechanical injury to Nav1.6-rich membranes causes cytoplasmic Na⁺-loading and “Nav-CLS”, i.e., coupled left-(hyperpolarizing)-shift of Nav activation and availability. Metabolic injury of hippocampal neurons (epileptic discharge) results in comparable impairment: left-shifted activation and availability and hence left-shifted I_Na-window. A recent computation study revealed that CLS-based I_Na-window left-shift dissipates ion gradients and impairs excitability. Here, via dynamical analyses, we focus on sustained excitability patterns in mildly damaged nodes, in particular with more realistic Gaussian-distributed Nav-CLS to mimic “smeared” injury intensity. Since our interest is axons that might survive injury, pumps (sine qua non for live axons) are included. In some simulations, pump efficacy and system volumes are varied. Impacts of current noise inputs are also characterized. The diverse modes of spontaneous rhythmic activity evident in these scenarios are studied using bifurcation analysis. For “mild CLS injury”, a prominent feature is slow pump/leak-mediated E_Ion oscillations. These slow oscillations yield dynamic firing thresholds that underlie complex voltage STO and bursting behaviors. Thus, Nav-CLS, a biophysically justified mode of injury, in parallel with functioning pumps, robustly engenders an emergent slow process that triggers a plethora of pathological excitability patterns. This minimalist “device” could have physiological analogs. At first nodes of Ranvier and at nociceptors, e.g., localized lipid-tuning that modulated Nav midpoints could produce Nav-CLS, as could co-expression of appropriately differing Nav isoforms.

Nerve cells damaged by trauma, stroke, epilepsy, inflammatory conditions etc, have chronically leaky sodium channels that eventually kill. The usual job of sodium channels is to make brief voltage signals –action potentials– for long distance propagation. After sodium channels open to generate action potentials, sodium pumps work harder to re-establish the intracellular/extracellular sodium imbalance that is, literally, the neuron's battery for firing action potentials. Wherever tissue damage renders membranes overly fluid, we hypothesize, sodium channels become chronically leaky. Our experimental findings justify this. In fluidized membranes, sodium channel voltage sensors respond too easily, letting channels spend too much time open. Channels leak, pumps respond. By mathematical modeling, we show that in damaged channel-rich membranes the continual pump/leak counterplay would trigger the kinds of bizarre intermittent action potential bursts typical of injured neurons. Arising ectopically from injury regions, such neuropathic firing is unrelated to events in the external world. Drugs that can silence these deleterious electrical barrages without blocking healthy action potentials are needed. If fluidized membranes house the problematic leaky sodium channels, then drug side effects could be diminished by using drugs that accumulate most avidly into fluidized membranes, and that bind their targets with highest affinity there.

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.

Distinct developmental patterns in the expression of transient, persistent, and resurgent Na+ currents in entorhinal cortex layer-II neurons

Sub- and near-threshold voltage-dependent Na+ currents (VDSCs) are of major importance in determining the electrical properties of medial entorhinal cortex (mEC) layer-II neurons. Developmental changes in the ability of mEC layer-II stellate cells (SCs) to generate Na+ -dependent, subthreshold electrical events have been reported between P14 and P18. In this study we examined the modifications occurring in the various components of VDSCs during postnatal development of mEC SCs. The transient, resurgent, and persistent Na+ currents (I(NaT), I(NaR), and I(NaP), respectively) showed distinct patterns of developmental expression in the time window considered (P5 to P24-27). All three currents prominently and steeply increased in absolute amplitude and conductance from P5 to at least P16. However, capacitive charge accumulation, an index of membrane surface area, also markedly increased in the same time window, and in the case of I(NaT) the specific conductance per unit of accumulated capacitive charge remained relatively constant. By contrast, specific I(NaR) and I(NaP) conductances showed a significant tendency to increase, especially from P5 to P18. Neither I(NaR) nor I(NaP) represented a constant fraction of the total Na+ current at all developmental ages. Indeed, detectable levels of I(NaR) and I(NaP) were present in only ~20% and ~70%, respectively, of the cells on P5, and were observed in all cells only from P10 onwards. Moreover, the average I(NaR)-to-I(NaT) conductance ratio increased steadily from ~0.004 (P5) up to a plateau level of ~0.05 (P22+), whereas the I(NaP)-to-I(NaT) conductance ratio increased only from ~0.009 on P5 to ~0.02 on P22+. The relative increase in conductance ratio from P5 to P22 was significantly greater for I(NaR) than for I(NaP), indicating that I(NaR) expression starts later than that of I(NaP). These findings show that in mEC layer-II SCs the single functional components of the VDSC are regulated differentially from each other as far as their developmental expression is concerned.