Spiking resonances in models with the same slow resonant and fast amplifying currents but different subthreshold dynamic properties

The generation of spiking resonances in neurons (preferred spiking responses to oscillatory inputs) requires the interplay of the intrinsic ionic currents that operate at the subthreshold voltage level and the spiking mechanisms. Combinations of the same types of ionic currents in different parameter regimes may give rise to different types of nonlinearities in the voltage equation (e.g., parabolic- and cubic-like), generating subthreshold (membrane potential) oscillations patterns with different properties. These nonlinearities are not apparent in the model equations, but can be uncovered by plotting the voltage nullclines in the phase-plane diagram. We investigate the spiking resonant properties of conductance-based models that are biophysically equivalent at the subthreshold level (same ionic currents), but dynamically different (parabolic- and cubic-like voltage nullclines). As a case study we consider a model having a persistent sodium and a hyperpolarization-activated (h-) currents, which exhibits subthreshold resonance in the theta frequency band. We unfold the concept of spiking resonance into evoked and output spiking resonance. The former focuses on the input frequencies that are able to generate spikes, while the latter focuses on the output spiking frequencies regardless of the input frequency that generated these spikes. A cell can exhibit one or both types of resonances. We also measure spiking phasonance, which is an extension of subthreshold phasonance (zero-phase-shift response to oscillatory inputs) to the spiking regime. The subthreshold resonant properties of both types of models are communicated to the spiking regime for low enough input amplitudes as the voltage response for the subthreshold resonant frequency band raises above threshold. For higher input amplitudes evoked spiking resonance is no longer present in these models, but output spiking resonance is present primarily in the parabolic-like model due to a cycle skipping mechanism (involving mixed-mode oscillations), while the cubic-like model shows a better 1:1 entrainment. We use dynamical systems tools to explain the underlying mechanisms and the mechanistic differences between the resonance types. Our results demonstrate that the effective time scales that operate at the subthreshold regime to generate intrinsic subthreshold oscillations, mixed-mode oscillations and subthreshold resonance do not necessarily determine the existence of a preferred spiking response to oscillatory inputs in the same frequency band. The results discussed in this paper highlight both the complexity of the suprathreshold responses to oscillatory inputs in neurons having resonant and amplifying currents with different time scales and the fact that the identity of the participating ionic currents is not enough to predict the resulting patterns, but additional dynamic information, captured by the geometric properties of the phase-space diagram, is needed.     

Electrical Resonance in the θ Frequency Range in Olfactory Amygdala Neurons

The cortical amygdala receives direct olfactory inputs and is thought to participate in processing and learning of biologically relevant olfactory cues. As for other brain structures implicated in learning, the principal neurons of the anterior cortical nucleus (ACo) exhibit intrinsic subthreshold membrane potential oscillations in the θ-frequency range. Here we show that nearly 50% of ACo layer II neurons also display electrical resonance, consisting of selective responsiveness to stimuli of a preferential frequency (2–6 Hz). Their impedance profile resembles an electrical band-pass filter with a peak at the preferred frequency, in contrast to the low-pass filter properties of other neurons. Most ACo resonant neurons displayed frequency preference along the whole subthreshold voltage range. We used pharmacological tools to identify the voltage-dependent conductances implicated in resonance. A hyperpolarization-activated cationic current depending on HCN channels underlies resonance at resting and hyperpolarized potentials; notably, this current also participates in resonance at depolarized subthreshold voltages. KV7/KCNQ K⁺ channels also contribute to resonant behavior at depolarized potentials, but not in all resonant cells. Moreover, resonance was strongly attenuated after blockade of voltage-dependent persistent Na⁺ channels, suggesting an amplifying role. Remarkably, resonant neurons presented a higher firing probability for stimuli of the preferred frequency. To fully understand the mechanisms underlying resonance in these neurons, we developed a comprehensive conductance-based model including the aforementioned and leak conductances, as well as Hodgkin and Huxley-type channels. The model reproduces the resonant impedance profile and our pharmacological results, allowing a quantitative evaluation of the contribution of each conductance to resonance. It also replicates selective spiking at the resonant frequency and allows a prediction of the temperature-dependent shift in resonance frequency. Our results provide a complete characterization of the resonant behavior of olfactory amygdala neurons and shed light on a putative mechanism for network activity coordination in the intact brain.     

Information coding by ensembles of resonant neurons

In the present paper, we propose a novel neural procedure for signal processing and coding based on the subthreshold oscillations and resonance of the neural membrane potential that could be used by real neurons to perform frequency spectra analysis and information coding of incoming signals. Taking into account the biophysical properties of the neural membranes, we note that the subthreshold resonant behaviour they exhibit can be used to analyse incoming signals and represent them in the frequency domain. We study the reliability of the representation of signals depending on the biophysical parameters of the neurons, the fault-tolerance of this coding scheme and its robustness against noise and in the presence of spikes. The principal characteristics of our system are the use of the physical phenomenon of neural resonance (rarely considered in the literature for signal coding); it fits well with the biophysical parameters of most neurons that exhibit subthreshold oscillations; it is compatible with experimental data; and it can be easily integrated in a more general model of information processing and coding that includes communication between neurons based on spikes.     

Long-term potentiation in rat hippocampal neurons is accompanied by spatially widespread changes in intrinsic oscillatory dynamics and excitability

Oscillations in neural activity are a prominent feature of many brain states. Individual hippocampal neurons exhibit intrinsic membrane potential oscillations and intrinsic resonance in the theta frequency range. We found that the subthreshold resonance frequency of CA1 pyramidal neurons was location dependent, varying more than 3-fold between the soma and the distal dendrites. Furthermore, activity- and NMDA-receptor-dependent long-term plasticity increased this resonance frequency through changes in h channel properties. The increase in resonance frequency and an associated reduction in excitability were nearly identical in the soma and the first 300 mum of the apical dendrites. These spatially widespread changes accompanying long-term synaptic potentiation also reduced the neuron's ability to elicit spikes evoked through a nonpotentiated synaptic pathway. Our results suggest that the frequency response of these neurons depends on the dendritic location of their inputs and that activity can regulate their response dynamics within an oscillating neural network.     

Resonances and noise in a stochastic Hindmarsh-Rose model of thalamic neurons

Thalamic neurons exhibit subthreshold resonance when stimulated with small sine wave signals of varying frequency and stochastic resonance when noise is added to these signals. We study a stochastic Hindmarsh-Rose model using Monte-Carlo simulations to investigate how noise, in conjunction with subthreshold resonance, leads to a preferred frequency in the firing pattern. The resulting stochastic resonance (SR) exhibits a preferred firing frequency that is approximately exponential in its dependence on the noise amplitude. In similar experiments, frequency dependent SR is found in the reliability of detection of alpha-function inputs under noise, which are more realistic inputs for neurons. A mathematical analysis of the equations reveals that the frequency preference arises from the dynamics of the slow variable. Noise can then transfer the resonance over the firing threshold because of the proximity of the fast subsystem to a Hopf bifurcation point. Our results may have implications for the behavior of thalamic neurons in a network, with noise switching the membrane potential between different resonance modes.     

Membrane Resonance and Stochastic Resonance Modulate Firing Patterns of Thalamocortical Neurons

We examined the interactions of subthreshold membrane resonance and stochastic resonance using whole-cell patch clamp recordings in thalamocortical neurons of rat brain slices, as well as with a Hodgkin-Huxley-type mathematical model of thalamocortical neurons. The neurons exhibited the subthreshold resonance when stimulated with small amplitude sine wave currents of varying frequency, and stochastic resonance when noise was added to sine wave inputs. Stochastic resonance was manifest as a maximum in signal-to-noise ratio of output response to subthreshold periodic input combined with noise. Stochastic resonance in conjunction with subthreshold resonance resulted in action potential patterns that showed frequency selectivity for periodic inputs. Stochastic resonance was maximal near subthreshold resonance frequency and a high noise level was required for detection of high frequency signals. We speculate that combined membrane and stochastic resonances have physiological utility in coupling synaptic activity to preferred firing frequency and in network synchronization under noise.     

Interspike interval densities of resonate and fire neurons

The subthreshold dynamics of a neuron can follow one of the two patterns: resonant neurons generate intrinsic subthreshold membrane potential oscillations, whereas in nonresonant neurons these oscillations are not observed. Here, we investigate how these subthreshold behaviors affect the suprathreshold response. Both types of neurons are described by a resonate and fire model, with the stable fixpoint being either a focus or a node. Using analytic expression for a linear oscillator model with threshold and reset, we calculate the multimodal interspike interval densities. We show that a change in model parameters induces qualitative changes in the interspike interval densities.     

Noise enhances subthreshold oscillations in injured primary sensory neurons

Noise can play a constructive role in the detection of weak signals in various kinds of peripheral receptors and neurons. What the mechanism underlying the effect of noise is remains unclear. Here, the perforated patch-clamp technique was used on isolated cells from chronic compression of the dorsal root ganglion (DRG) model. Our data provided new insight indicating that, under conditions without external signals, noise can enhance subthreshold oscillations, which was observed in a certain type of neurons with high-frequency (20-100 Hz) intrinsic resonance from injured DRG neurons. The occurrence of subthreshold oscillation considerably decreased the threshold potential for generating repetitive firing. The above effects of noise can be abolished by blocking the persistent sodium current (I(Na, P)). Utilizing a mathematical neuron model we further simulated the effect of noise on subthreshold oscillation and firing, and also found that noise can enhance the electrical activity through autonomous stochastic resonance. Accordingly, we propose a new concept of the effects of noise on neural intrinsic activity, which suggests that noise may be an important factor for modulating the excitability of neurons and generation of chronic pain signals.     

Coherence resonance for neuronal bursting with spike undershoot

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

Squid Axon Membrane Response to White Noise Stimulation

The current from a white noise generator was applied as a stimulus to a space-clamped squid axon in double sucrose gap. The membrane current and the voltage response of the membrane were then amplified, recorded on magnetic tape, and the stimulus was cross-correlated with the response. With subthreshold stimuli, a cross-correlation function resembling that obtained from a resonant parallel circuit is obtained. As the intensity of the input noise is increased, the cross-correlation function resembles that obtained from a less damped oscillatory circuit. When the noise intensity is further increased so that an appreciable frequency of action potentials is observed, an additional component appears in the experimental cross-correlogram. The subthreshold cross-correlogram is analyzed theoretically in terms of the linearized Hodgkin-Huxley equations. The subthreshold axon approximates a parallel resonant circuit. The circuit parameters are temperature dependent, with resonant frequency varying from approximately 100 Hz at 10°C to approximately 250 Hz at 20°C. The Q₁₀ of the resonant frequency is equal to 1.9. These values are in agreement with values found previously for subthreshold oscillations following a single action potential.     

History-dependent excitability as a single-cell substrate of transient memory for information discrimination

Neurons react differently to incoming stimuli depending upon their previous history of stimulation. This property can be considered as a single-cell substrate for transient memory, or context-dependent information processing: depending upon the current context that the neuron "sees" through the subset of the network impinging on it in the immediate past, the same synaptic event can evoke a postsynaptic spike or just a subthreshold depolarization. We propose a formal definition of History-Dependent Excitability (HDE) as a measure of the propensity to firing in any moment in time, linking the subthreshold history-dependent dynamics with spike generation. This definition allows the quantitative assessment of the intrinsic memory for different single-neuron dynamics and input statistics. We illustrate the concept of HDE by considering two general dynamical mechanisms: the passive behavior of an Integrate and Fire (IF) neuron, and the inductive behavior of a Generalized Integrate and Fire (GIF) neuron with subthreshold damped oscillations. This framework allows us to characterize the sensitivity of different model neurons to the detailed temporal structure of incoming stimuli. While a neuron with intrinsic oscillations discriminates equally well between input trains with the same or different frequency, a passive neuron discriminates better between inputs with different frequencies. This suggests that passive neurons are better suited to rate-based computation, while neurons with subthreshold oscillations are advantageous in a temporal coding scheme. We also address the influence of intrinsic properties in single-cell processing as a function of input statistics, and show that intrinsic oscillations enhance discrimination sensitivity at high input rates. Finally, we discuss how the recognition of these cell-specific discrimination properties might further our understanding of neuronal network computations and their relationships to the distribution and functional connectivity of different neuronal types.     

Cortical Resonance Frequencies Emerge from Network Size and Connectivity

Neural oscillations occur within a wide frequency range with different brain regions exhibiting resonance-like characteristics at specific points in the spectrum. At the microscopic scale, single neurons possess intrinsic oscillatory properties, such that is not yet known whether cortical resonance is consequential to neural oscillations or an emergent property of the networks that interconnect them. Using a network model of loosely-coupled Wilson-Cowan oscillators to simulate a patch of cortical sheet, we demonstrate that the size of the activated network is inversely related to its resonance frequency. Further analysis of the parameter space indicated that the number of excitatory and inhibitory connections, as well as the average transmission delay between units, determined the resonance frequency. The model predicted that if an activated network within the visual cortex increased in size, the resonance frequency of the network would decrease. We tested this prediction experimentally using the steady-state visual evoked potential where we stimulated the visual cortex with different size stimuli at a range of driving frequencies. We demonstrate that the frequency corresponding to peak steady-state response inversely correlated with the size of the network. We conclude that although individual neurons possess resonance properties, oscillatory activity at the macroscopic level is strongly influenced by network interactions, and that the steady-state response can be used to investigate functional networks.     

Internal noise-driven circadian oscillator in Drosophila

An internal noise-driven oscillator was studied in a two-variable Drosophila model, where both positive feedback and negative feedback are crucial to the circadian oscillations. It is shown that internal noise could sustain reliable oscillations for the parameter which produces a stable steady state in the deterministic system. The noise-sustained oscillations are interpreted by using phase plane analysis. The period of such oscillations fluctuates slightly around the period of deterministic oscillations and the coherence of oscillations becomes the best at an optimal internal noise intensity, indicating the occurrence of intrinsic coherence resonance. In addition, in the oscillatory region, the coherence of noisy circadian oscillations is suppressed by the internal noise, but the period is hardly affected, demonstrating the robustness of the Drosophila model for circadian rhythms to the intrinsic noise.     

Interplay of Intrinsic and Synaptic Conductances in the Generation of High-Frequency Oscillations in Interneuronal Networks with Irregular Spiking

High-frequency oscillations (above 30 Hz) have been observed in sensory and higher-order brain areas, and are believed to constitute a general hallmark of functional neuronal activation. Fast inhibition in interneuronal networks has been suggested as a general mechanism for the generation of high-frequency oscillations. Certain classes of interneurons exhibit subthreshold oscillations, but the effect of this intrinsic neuronal property on the population rhythm is not completely understood. We study the influence of intrinsic damped subthreshold oscillations in the emergence of collective high-frequency oscillations, and elucidate the dynamical mechanisms that underlie this phenomenon. We simulate neuronal networks composed of either Integrate-and-Fire (IF) or Generalized Integrate-and-Fire (GIF) neurons. The IF model displays purely passive subthreshold dynamics, while the GIF model exhibits subthreshold damped oscillations. Individual neurons receive inhibitory synaptic currents mediated by spiking activity in their neighbors as well as noisy synaptic bombardment, and fire irregularly at a lower rate than population frequency. We identify three factors that affect the influence of single-neuron properties on synchronization mediated by inhibition: i) the firing rate response to the noisy background input, ii) the membrane potential distribution, and iii) the shape of Inhibitory Post-Synaptic Potentials (IPSPs). For hyperpolarizing inhibition, the GIF IPSP profile (factor iii)) exhibits post-inhibitory rebound, which induces a coherent spike-mediated depolarization across cells that greatly facilitates synchronous oscillations. This effect dominates the network dynamics, hence GIF networks display stronger oscillations than IF networks. However, the restorative current in the GIF neuron lowers firing rates and narrows the membrane potential distribution (factors i) and ii), respectively), which tend to decrease synchrony. If inhibition is shunting instead of hyperpolarizing, post-inhibitory rebound is not elicited and factors i) and ii) dominate, yielding lower synchrony in GIF networks than in IF networks.     

Frequency sensitivity in Hodgkin–Huxley systems

 The frequency sensitivity of weak periodic signal detection has been studied via numerical simulations for both a single neuron and a neuronal network. The dependence of the critical amplitude of the signal upon its frequency and a resonance between the intrinsic oscillations of a neuron and the signal could account for the frequency sensitivity. In the presence of both a subthreshold periodic signal and noise, the signal-to-noise ratio (SNR) of the output of either a single neuron or a neuronal network present the typical characteristics of stochastic resonance. In particular, there exists a frequency-sensitive range of 30–100 Hz, and for signals with frequencies within this range the SNRs have large values. This implies that the system under consideration (a single neuron or a neuronal network) is more sensitive to the detection of periodic signals, and the frequency sensitivity may be of a functional significance to signal processing. Received: 26 October 1999 / Accepted in revised form: 25 July 2000     

The generation of phase differences and frequency changes in a network model of inferior olive subthreshold oscillations

It is commonly accepted that the Inferior Olive (IO) provides a timing signal to the cerebellum. Stable subthreshold oscillations in the IO can facilitate accurate timing by phase-locking spikes to the peaks of the oscillation. Several theoretical models accounting for the synchronized subthreshold oscillations have been proposed, however, two experimental observations remain an enigma. The first is the observation of frequent alterations in the frequency of the oscillations. The second is the observation of constant phase differences between simultaneously recorded neurons. In order to account for these two observations we constructed a canonical network model based on anatomical and physiological data from the IO. The constructed network is characterized by clustering of neurons with similar conductance densities, and by electrical coupling between neurons. Neurons inside a cluster are densely connected with weak strengths, while neurons belonging to different clusters are sparsely connected with stronger connections. We found that this type of network can robustly display stable subthreshold oscillations. The overall frequency of the network changes with the strength of the inter-cluster connections, and phase differences occur between neurons of different clusters. Moreover, the phase differences provide a mechanistic explanation for the experimentally observed propagating waves of activity in the IO. We conclude that the architecture of the network of electrically coupled neurons in combination with modulation of the inter-cluster coupling strengths can account for the experimentally observed frequency changes and the phase differences.     

Voltage dependence of subthreshold resonance frequency in layer II of medial entorhinal cortex

The resonance properties of individual neurons in entorhinal cortex (EC) may contribute to their functional properties in awake, behaving rats. Models propose that entorhinal grid cells could arise from shifts in the intrinsic frequency of neurons caused by changes in membrane potential owing to depolarizing input from neurons coding velocity. To test for potential changes in intrinsic frequency, we measured the resonance properties of neurons at different membrane potentials in neurons in medial and lateral EC. In medial entorhinal neurons, the resonant frequency of individual neurons decreased in a linear manner as the membrane potential was depolarized between -70 and -55 mV. At more hyperpolarized membrane potentials, cells asymptotically approached a maximum resonance frequency. Consistent with the previous studies, near resting potential, the cells of the medial EC possessed a decreasing gradient of resonance frequency along the dorsal to ventral axis, and cells of the lateral EC lacked resonant properties, regardless of membrane potential or position along the medial to lateral axis within lateral EC. Application of 10 μM ZD7288, the H-channel blocker, abolished all resonant properties in MEC cells, and resulted in physiological properties very similar to lateral EC cells. These results on resonant properties show a clear change in frequency response with depolarization that could contribute to the generation of grid cell firing properties in the medial EC.     

Spike latency and jitter of neuronal membrane patches with stochastic Hodgkin-Huxley channels

Ion channel stochasticity can influence the voltage dynamics of neuronal membrane, with stronger effects for smaller patches of membrane because of the correspondingly smaller number of channels. We examine this question with respect to first spike statistics in response to a periodic input of membrane patches including stochastic Hodgkin-Huxley channels, comparing these responses to spontaneous firing. Without noise, firing threshold of the model depends on frequency—a sinusoidal stimulus is subthreshold for low and high frequencies and suprathreshold for intermediate frequencies. When channel noise is added, a stimulus in the lower range of subthreshold frequencies can influence spike output, while high subthreshold frequencies remain subthreshold. Both input frequency and channel noise strength influence spike timing. Specifically, spike latency and jitter have distinct minima as a function of input frequency, showing a resonance like behavior. With either no input, or low frequency subthreshold input, or input in the low or high suprathreshold frequency range, channel noise reduces latency and jitter, with the strongest impact for the lowest input frequencies. In contrast, for an intermediate range of suprathreshold frequencies, where an optimal input gives a minimum latency, the noise effect reverses, and spike latency and jitter increase with channel noise. Thus, a resonant minimum of the spike response as a function of frequency becomes more pronounced with less noise. Spike latency and jitter also depend on the initial phase of the input, resulting in minimal latencies at an optimal phase, and depend on the membrane time constant, with a longer time constant broadening frequency tuning for minimal latency and jitter. Taken together, these results suggest how stochasticity of ion channels may influence spike timing and thus coding for neurons with functionally localized concentrations of channels, such as in “hot spots” of dendrites, spines or axons.     

Internal noise-sustained circadian rhythms in a Drosophila model

Circadian rhythmic processes, mainly regulated by gene expression at the molecular level, have inherent stochasticity. Their robustness or resistance to internal noise has been extensively investigated by most of the previous studies. This work focuses on the constructive roles of internal noise in a reduced Drosophila model, which incorporates negative and positive feedback loops, each with a time delay. It is shown that internal noise sustains reliable oscillations with periods close to 24 h in a region of parameter space, where the deterministic kinetics would evolve to a stable steady state. The amplitudes of noise-sustained oscillations are significantly affected by the variation of internal noise level, and the best performance of the oscillations could be found at an optimal noise intensity, indicating the occurrence of intrinsic coherence resonance. In the oscillatory region of the deterministic model, the coherence of noisy circadian oscillations is suppressed by internal noise, while the period remains nearly constant over a large range of noise intensity, demonstrating robustness of the Drosophila model for circadian rhythms to intrinsic noise. In addition, the effects of time delay in the positive feedback on the oscillations are also investigated. It is found that the time delay could efficiently tune the performance of the noise-sustained oscillations. These results might aid understanding of the exploitation of intracellular noise in biochemical and genetic regulatory systems.     

How could a chirp be more effective than a louder clock--resonant transfer of energy between subthreshold excitation pulses and excitable tissues

By definition, electrical resonance of excitable biomembranes reflects their remarkable ability to select input signals with preferential frequencies as ones that will pass through the chain of information exchange within a coordinated organism. This feature is known to be caused and regulated by the voltage-dependent kinetics of transmembrane ion channels. Knowledge about the resonant behavior of excitable cells may prove crucial toward better underpinning those essential properties of neurons which make them to function as a highly coordinated network. In this work we present novel experimental evidence which strongly supports the paradigm of selective communication between excitable cells, by no more than efficient transfer of energy which takes place when the frequency of subthreshold presynaptic pulses matches the natural, resonant frequency of a target excitable tissue. As a possible direct application, our novel experimental results prove the functional importance of bursting activity of pre-synaptic neurons as an effective physiological mechanism for discriminatory communication between neurons.