Membrane potential resonance frequency directly influences network frequency through electrical coupling

Distinct Mechanisms Underlie Electrical Coupling Resonance and Its Interaction with Membrane Potential Resonance

Neurons in oscillatory networks often exhibit membrane potential resonance, a peak impedance at a non-zero input frequency. In electrically coupled oscillatory networks, the coupling coefficient (the ratio of post- and prejunctional voltage responses) could also show resonance. Such coupling resonance may emerge from the interaction between the coupling current and resonance properties of the coupled neurons, but this relationship has not been clearly described. Additionally, it is unknown if the gap-junction mediated electrical coupling conductance may have frequency dependence. We examined these questions by recording a pair of electrically coupled neurons in the oscillatory pyloric network of the crab Cancer borealis. We performed dual current- and voltage-clamp recordings and quantified the frequency preference of the coupled neurons, the coupling coefficient, the electrical conductance, and the postjunctional neuronal response. We found that all components exhibit frequency selectivity, but with distinct preferred frequencies. Mathematical and computational analysis showed that membrane potential resonance of the postjunctional neuron was sufficient to give rise to resonance properties of the coupling coefficient, but not the coupling conductance. A distinct coupling conductance resonance frequency therefore emerges either from other circuit components or from the gating properties of the gap junctions. Finally, to explore the functional effect of the resonance of the coupling conductance, we examined its role in synchronizing neuronal the activities of electrically coupled bursting model neurons. Together, our findings elucidate factors that produce electrical coupling resonance and the function of this resonance in oscillatory networks.

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.

Frequency-preference response in covalent modification cycles under substrate sequestration conditions

Covalent modification cycles (CMCs) are basic units of signaling systems and their properties are well understood. However, their behavior has been mostly characterized in situations where the substrate is in excess over the modifying enzymes. Experimental data on protein abundance suggest that the enzymes and their target proteins are present in comparable concentrations, leading to substrate sequestration by the enzymes. In this enzyme-in-excess regime, CMCs have been shown to exhibit signal termination, the ability of the product to return to a stationary value lower than its peak in response to constant stimulation, while this stimulation is still active, with possible implications for the ability of systems to adapt to environmental inputs. We characterize the conditions leading to signal termination in CMCs in the enzyme-in-excess regime. We also demonstrate that this behavior leads to a preferred frequency response (band-pass filters) when the cycle is subjected to periodic stimulation, whereas the literature reports that CMCs investigated so far behave as low-pass filters. We characterize the relationship between signal termination and the preferred frequency response to periodic inputs and we explore the dynamic mechanism underlying these phenomena. Finally, we describe how the behavior of CMCs is reflected in similar types of responses in the cascades of which they are part. Evidence of protein abundance in vivo shows that enzymes and substrates are present in comparable concentrations, thus suggesting that signal termination and frequency-preference response to periodic inputs are also important dynamic features of cell signaling systems, which have been overlooked.

Asymmetrical voltage response in resonant neurons shaped by nonlinearities

The conventional impedance profile of a neuron can identify the presence of resonance and other properties of the neuronal response to oscillatory inputs, such as nonlinear response amplifications, but it cannot distinguish other nonlinear properties such as asymmetries in the shape of the voltage response envelope. Experimental observations have shown that the response of neurons to oscillatory inputs preferentially enhances either the upper or lower part of the voltage envelope in different frequency bands. These asymmetric voltage responses arise in a neuron model when it is submitted to high enough amplitude oscillatory currents of variable frequencies. We show how the nonlinearities associated to different ionic currents or present in the model as captured by its voltage equation lead to asymmetrical response and how high amplitude oscillatory currents emphasize this response. We propose a geometrical explanation for the phenomenon where asymmetries result not only from nonlinearities in their activation curves but also from nonlinearites captured by the nullclines in the phase-plane diagram and from the system's time-scale separation. In addition, we identify an unexpected frequency-dependent pattern which develops in the gating variables of these currents and is a product of strong nonlinearities in the system as we show by controlling such behavior by manipulating the activation curve parameters. The results reported in this paper shed light on the ionic mechanisms by which brain embedded neurons process oscillatory information.

Frequency-dependent responses of neuronal models to oscillatory inputs in current versus voltage clamp

Action potential generation in neurons depends on a membrane potential threshold and therefore on how subthreshold inputs influence this voltage. In oscillatory networks, for example, many neuron types have been shown to produce membrane potential ([Formula: see text]) resonance: a maximum subthreshold response to oscillatory inputs at a nonzero frequency. Resonance is usually measured by recording [Formula: see text] in response to a sinusoidal current ([Formula: see text]), applied at different frequencies (f), an experimental setting known as current clamp (I-clamp). Several recent studies, however, use the voltage clamp (V-clamp) method to control [Formula: see text] with a sinusoidal input at different frequencies [[Formula: see text]] and measure the total membrane current ([Formula: see text]). The two methods obey systems of differential equations of different dimensionality, and while I-clamp provides a measure of electrical impedance [[Formula: see text]], V-clamp measures admittance [[Formula: see text]]. We analyze the relationship between these two measurement techniques. We show that, despite different dimensionality, in linear systems the two measures are equivalent: [Formula: see text]. However, nonlinear model neurons produce different values for Z and [Formula: see text]. In particular, nonlinearities in the voltage equation produce a much larger difference between these two quantities than those in equations of recovery variables that describe activation and inactivation kinetics. Neurons are inherently nonlinear, and notably, with ionic currents that amplify resonance, the voltage clamp technique severely underestimates the current clamp response. We demonstrate this difference experimentally using the PD neurons in the crab stomatogastric ganglion. These findings are instructive for researchers who explore cellular mechanisms of neuronal oscillations.

Membrane potential resonance in non-oscillatory neurons interacts with synaptic connectivity to produce network oscillations

Several neuron types have been shown to exhibit (subthreshold) membrane potential resonance (MPR), defined as the occurrence of a peak in their voltage amplitude response to oscillatory input currents at a preferred (resonant) frequency. MPR has been investigated both experimentally and theoretically. However, whether MPR is simply an epiphenomenon or it plays a functional role for the generation of neuronal network oscillations and how the latent time scales present in individual, non-oscillatory cells affect the properties of the oscillatory networks in which they are embedded are open questions. We address these issues by investigating a minimal network model consisting of (i) a non-oscillatory linear resonator (band-pass filter) with 2D dynamics, (ii) a passive cell (low-pass filter) with 1D linear dynamics, and (iii) nonlinear graded synaptic connections (excitatory or inhibitory) with instantaneous dynamics. We demonstrate that (i) the network oscillations crucially depend on the presence of MPR in the resonator, (ii) they are amplified by the network connectivity, (iii) they develop relaxation oscillations for high enough levels of mutual inhibition/excitation, and (iv) the network frequency monotonically depends on the resonators resonant frequency. We explain these phenomena using a reduced adapted version of the classical phase-plane analysis that helps uncovering the type of effective network nonlinearities that contribute to the generation of network oscillations. We extend our results to networks having cells with 2D dynamics. Our results have direct implications for network models of firing rate type and other biological oscillatory networks (e.g, biochemical, genetic).

Flexible resonance in prefrontal networks with strong feedback inhibition

Oscillations are ubiquitous features of brain dynamics that undergo task-related changes in synchrony, power, and frequency. The impact of those changes on target networks is poorly understood. In this work, we used a biophysically detailed model of prefrontal cortex (PFC) to explore the effects of varying the spike rate, synchrony, and waveform of strong oscillatory inputs on the behavior of cortical networks driven by them. Interacting populations of excitatory and inhibitory neurons with strong feedback inhibition are inhibition-based network oscillators that exhibit resonance (i.e., larger responses to preferred input frequencies). We quantified network responses in terms of mean firing rates and the population frequency of network oscillation; and characterized their behavior in terms of the natural response to asynchronous input and the resonant response to oscillatory inputs. We show that strong feedback inhibition causes the PFC to generate internal (natural) oscillations in the beta/gamma frequency range (>15 Hz) and to maximize principal cell spiking in response to external oscillations at slightly higher frequencies. Importantly, we found that the fastest oscillation frequency that can be relayed by the network maximizes local inhibition and is equal to a frequency even higher than that which maximizes the firing rate of excitatory cells; we call this phenomenon population frequency resonance. This form of resonance is shown to determine the optimal driving frequency for suppressing responses to asynchronous activity. Lastly, we demonstrate that the natural and resonant frequencies can be tuned by changes in neuronal excitability, the duration of feedback inhibition, and dynamic properties of the input. Our results predict that PFC networks are tuned for generating and selectively responding to beta- and gamma-rhythmic signals due to the natural and resonant properties of inhibition-based oscillators. They also suggest strategies for optimizing transcranial stimulation and using oscillatory networks in neuromorphic engineering.

Interplay of activation kinetics and the derivative conductance determines resonance properties of neurons

In a neuron with hyperpolarization activated current (I_{h}), the correct input frequency leads to an enhancement of the output response. This behavior is known as resonance and is well described by the neuronal impedance. In a simple neuron model we derive equations for the neuron's resonance and we link its frequency and existence with the biophysical properties of I_{h}. For a small voltage change, the component of the ratio of current change to voltage change (dI/dV) due to the voltage-dependent conductance change (dg/dV) is known as derivative conductance (G_{h}^{Der}). We show that both G_{h}^{Der} and the current activation kinetics (characterized by the activation time constant τ_{h}) are mainly responsible for controlling the frequency and existence of resonance. The increment of both factors (G_{h}^{Der} and τ_{h}) greatly contributes to the appearance of resonance. We also demonstrate that resonance is voltage dependent due to the voltage dependence of G_{h}^{Der}. Our results have important implications and can be used to predict and explain resonance properties of neurons with the I_{h} current.

Gap junction plasticity as a mechanism to regulate network-wide oscillations

Cortical oscillations are thought to be involved in many cognitive functions and processes. Several mechanisms have been proposed to regulate oscillations. One prominent but understudied mechanism is gap junction coupling. Gap junctions are ubiquitous in cortex between GABAergic interneurons. Moreover, recent experiments indicate their strength can be modified in an activity-dependent manner, similar to chemical synapses. We hypothesized that activity-dependent gap junction plasticity acts as a mechanism to regulate oscillations in the cortex. We developed a computational model of gap junction plasticity in a recurrent cortical network based on recent experimental findings. We showed that gap junction plasticity can serve as a homeostatic mechanism for oscillations by maintaining a tight balance between two network states: asynchronous irregular activity and synchronized oscillations. This homeostatic mechanism allows for robust communication between neuronal assemblies through two different mechanisms: transient oscillations and frequency modulation. This implies a direct functional role for gap junction plasticity in information transmission in cortex.

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.

Mechanisms of generation of membrane potential resonance in a neuron with multiple resonant ionic currents

Neuronal membrane potential resonance (MPR) is associated with subthreshold and network oscillations. A number of voltage-gated ionic currents can contribute to the generation or amplification of MPR, but how the interaction of these currents with linear currents contributes to MPR is not well understood. We explored this in the pacemaker PD neurons of the crab pyloric network. The PD neuron MPR is sensitive to blockers of H- (I_H) and calcium-currents (I_Ca). We used the impedance profile of the biological PD neuron, measured in voltage clamp, to constrain parameter values of a conductance-based model using a genetic algorithm and obtained many optimal parameter combinations. Unlike most cases of MPR, in these optimal models, the values of resonant- (f_res) and phasonant- (f_{ϕ = 0}) frequencies were almost identical. Taking advantage of this fact, we linked the peak phase of ionic currents to their amplitude, in order to provide a mechanistic explanation the dependence of MPR on the I_Ca gating variable time constants. Additionally, we found that distinct pairwise correlations between I_Ca parameters contributed to the maintenance of f_res and resonance power (Q_Z). Measurements of the PD neuron MPR at more hyperpolarized voltages resulted in a reduction of f_res but no change in Q_Z. Constraining the optimal models using these data unmasked a positive correlation between the maximal conductances of I_H and I_Ca. Thus, although I_H is not necessary for MPR in this neuron type, it contributes indirectly by constraining the parameters of I_Ca.

Many neuron types exhibit membrane potential resonance (MPR) in which the neuron produces the largest response to oscillatory input at some preferred (resonant) frequency and, in many systems, the network frequency is correlated with neuronal MPR. MPR is captured by a peak in the impedance vs. frequency curve (Z-profile), which is shaped by the dynamics of voltage-gated ionic currents. Although neuron types can express variable levels of ionic currents, they may have a stable resonant frequency. We used the PD neuron of the crab pyloric network to understand how MPR emerges from the interplay of the biophysical properties of multiple ionic currents, each capable of generating resonance. We show the contribution of an inactivating current at the resonant frequency in terms of interacting time constants. We measured the Z-profile of the PD neuron and explored possible combinations of model parameters that fit this experimentally measured profile. We found that the Z-profile constrains and defines correlations among parameters associated with ionic currents. Furthermore, the resonant frequency and amplitude are sensitive to different parameter sets and can be preserved by co-varying pairs of parameters along their correlation lines. Furthermore, although a resonant current may be present in a neuron, it may not directly contribute to MPR, but constrain the properties of other currents that generate MPR. Finally, constraining model parameters further to those that modify their MPR properties to changes in voltage range produces maximal conductance correlations.

Resonance modulation, annihilation and generation of anti-resonance and anti-phasonance in 3D neuronal systems: interplay of resonant and amplifying currents with slow dynamics

Subthreshold (membrane potential) resonance and phasonance (preferred amplitude and zero-phase responses to oscillatory inputs) in single neurons arise from the interaction between positive and negative feedback effects provided by relatively fast amplifying currents and slower resonant currents. In 2D neuronal systems, amplifying currents are required to be slave to voltage (instantaneously fast) for these phenomena to occur. In higher dimensional systems, additional currents operating at various effective time scales may modulate and annihilate existing resonances and generate antiresonance (minimum amplitude response) and antiphasonance (zero-phase response with phase monotonic properties opposite to phasonance). We use mathematical modeling, numerical simulations and dynamical systems tools to investigate the mechanisms underlying these phenomena in 3D linear models, which are obtained as the linearization of biophysical (conductance-based) models. We characterize the parameter regimes for which the system exhibits the various types of behavior mentioned above in the rather general case in which the underlying 2D system exhibits resonance. We consider two cases: (i) the interplay of two resonant gating variables, and (ii) the interplay of one resonant and one amplifying gating variables. Increasing levels of an amplifying current cause (i) a response amplification if the amplifying current is faster than the resonant current, (ii) resonance and phasonance attenuation and annihilation if the amplifying and resonant currents have identical dynamics, and (iii) antiresonance and antiphasonance if the amplifying current is slower than the resonant current. We investigate the underlying mechanisms by extending the envelope-plane diagram approach developed in previous work (for 2D systems) to three dimensions to include the additional gating variable, and constructing the corresponding envelope curves in these envelope-space diagrams. We find that antiresonance and antiphasonance emerge as the result of an asymptotic boundary layer problem in the frequency domain created by the different balances between the intrinsic time constants of the cell and the input frequency f as it changes. For large enough values of f the envelope curves are quasi-2D and the impedance profile decreases with the input frequency. In contrast, for f ≪ 1 the dynamics are quasi-1D and the impedance profile increases above the limiting value in the other regime. Antiresonance is created because the continuity of the solution requires the impedance profile to connect the portions belonging to the two regimes. If in doing so the phase profile crosses the zero value, then antiphasonance is also generated.

Bifurcation transitions in gap-junction-coupled neurons

Here we investigate transitions occurring in the dynamical states of pairs of distinct neurons electrically coupled, with one neuron tonic and the other bursting. Depending on the dynamics of the individual neurons, and for strong enough coupling, they synchronize either in a tonic or a bursting regime, or initially tonic transitioning to bursting via a period doubling cascade. Certain intrinsic properties of the individual neurons such as minimum firing rates are carried over into the dynamics of the coupled neurons affecting their ultimate synchronous state.