The Roles of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slow-wave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking (RS) cells with frequency adaptation and do not exhibit bursts in current-clamp experiments (in vitro). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings (gsyn) and external current inputs (I). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to the bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold (IL) and low-threshold (IT) ion channels play a role in increasing and decreasing the parameter conditions (gsyn and I) in which bistable dynamics occur, respectively. For high values of IL conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance IT increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.

The Role of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slowwave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking cells (RS) with frequency adaptation and do not exhibit bursts in current-clamp experiments ( in vitro ). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings ( g syn ) and external current inputs ( I ). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold ( I L ) and low-threshold ( I T ) ion channels play a role in increasing and decreasing the parameter conditions ( g syn and I ) in which bistable dynamics occur, respectively. For high values of I L conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance I T increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.

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.

Modeling Hippocampal CA1 Gabaergic Synapses of Audiogenic Rats

Wistar Audiogenic Rats (WARs) are genetically susceptible to sound-induced seizures that start in the brainstem and, in response to repetitive stimulation, spread to limbic areas, such as hippocampus. Analysis of the distribution of interevent intervals of GABAergic inhibitory postsynaptic currents (IPSCs) in CA1 pyramidal cells showed a monoexponential trend in Wistar rats, suggestive of a homogeneous population of synapses, but a biexponential trend in WARs. Based on this, we hypothesize that there are two populations of GABAergic synaptic release sites in CA1 pyramidal neurons from WARs. To address this hypothesis, we used a well-established neuronal computational model of a CA1 pyramidal neuron previously developed to replicate physiological properties of these cells. Our simulations replicated the biexponential trend only when we decreased the release frequency of synaptic currents by a factor of six in at least 40% of distal synapses. Our results suggest that almost half of the GABAergic synapses of WARs have a drastically reduced spontaneous release frequency. The computational model was able to reproduce the temporal dynamics of GABAergic inhibition that could underlie susceptibility to the spread of seizures.

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.

Firing properties of ventral medullary respiratory neurons in sino-aortic denervated rats

NEW FINDINGS

What is the central question of this study? After sino-aortic denervation (SAD), rats present normal levels of mean arterial pressure (MAP), high MAP variability and changes in breathing. However, mechanisms involved in SAD-induced respiratory changes and their impact on the modulation of sympathetic activity remain unclear. Herein, we characterized the firing frequency of medullary respiratory neurons after SAD. What is the main finding and its importance? Sino-aortic denervation-induced prolonged inspiration was associated with a reduced interburst frequency of pre-inspiratory/inspiratory neurons and an increased long-term variability of late inspiratory neurons, but no changes were observed in the ramp-inspiratory and post-inspiratory neurons. This imbalance in the respiratory network might contribute to the modulation of sympathetic activity after SAD.

ABSTRACT

In previous studies, we documented that after sino-aortic denervation (SAD) in rats there are significant changes in the breathing pattern, but no significant changes in sympathetic activity and mean arterial pressure compared with sham-operated rats. However, the neural mechanisms involved in the respiratory changes after SAD and the extent to which they might contribute to the observed normal sympathetic activity and mean arterial pressure remain unclear. Here, we hypothesized that after SAD, rats present with changes in the firing frequency of the ventral medullary inspiratory and post-inspiratory neurons. To test this hypothesis, male Wistar rats underwent SAD or sham surgery and 3 days later were surgically prepared for an in situ experiment. The duration of inspiration significantly increased in SAD rats. During inspiration, the total firing frequency of ramp-inspiratory, pre-inspiratory/inspiratory and late-inspiratory neurons was not different between groups. During post-inspiration, the total firing frequency of post-inspiratory neurons was also not different between groups. Furthermore, the data demonstrate a reduced interburst frequency of pre-inspiratory/inspiratory neurons and an increased long-term variability of late-inspiratory neurons in SAD compared with sham-operated rats. These findings indicate that the SAD-induced prolongation of inspiration was not accompanied by alterations in the total firing frequency of the ventral medullary respiratory neurons, but it was associated with changes in the long-term variability of late-inspiratory neurons. We suggest that the timing imbalance in the respiratory network in SAD rats might contribute to the modulation of presympathetic neurons after removal of baroreceptor afferents.

Dynamics of spontaneous activity in random networks with multiple neuron subtypes and synaptic noise : Spontaneous activity in networks with synaptic noise

Spontaneous cortical population activity exhibits a multitude of oscillatory patterns, which often display synchrony during slow-wave sleep or under certain anesthetics and stay asynchronous during quiet wakefulness. The mechanisms behind these cortical states and transitions among them are not completely understood. Here we study spontaneous population activity patterns in random networks of spiking neurons of mixed types modeled by Izhikevich equations. Neurons are coupled by conductance-based synapses subject to synaptic noise. We localize the population activity patterns on the parameter diagram spanned by the relative inhibitory synaptic strength and the magnitude of synaptic noise. In absence of noise, networks display transient activity patterns, either oscillatory or at constant level. The effect of noise is to turn transient patterns into persistent ones: for weak noise, all activity patterns are asynchronous non-oscillatory independently of synaptic strengths; for stronger noise, patterns have oscillatory and synchrony characteristics that depend on the relative inhibitory synaptic strength. In the region of parameter space where inhibitory synaptic strength exceeds the excitatory synaptic strength and for moderate noise magnitudes networks feature intermittent switches between oscillatory and quiescent states with characteristics similar to those of synchronous and asynchronous cortical states, respectively. We explain these oscillatory and quiescent patterns by combining a phenomenological global description of the network state with local descriptions of individual neurons in their partial phase spaces. Our results point to a bridge from events at the molecular scale of synapses to the cellular scale of individual neurons to the collective scale of neuronal populations.

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.

Self-Consistent Scheme for Spike-Train Power Spectra in Heterogeneous Sparse Networks

Recurrent networks of spiking neurons can be in an asynchronous state characterized by low or absent cross-correlations and spike statistics which resemble those of cortical neurons. Although spatial correlations are negligible in this state, neurons can show pronounced temporal correlations in their spike trains that can be quantified by the autocorrelation function or the spike-train power spectrum. Depending on cellular and network parameters, correlations display diverse patterns (ranging from simple refractory-period effects and stochastic oscillations to slow fluctuations) and it is generally not well-understood how these dependencies come about. Previous work has explored how the single-cell correlations in a homogeneous network (excitatory and inhibitory integrate-and-fire neurons with nearly balanced mean recurrent input) can be determined numerically from an iterative single-neuron simulation. Such a scheme is based on the fact that every neuron is driven by the network noise (i.e., the input currents from all its presynaptic partners) but also contributes to the network noise, leading to a self-consistency condition for the input and output spectra. Here we first extend this scheme to homogeneous networks with strong recurrent inhibition and a synaptic filter, in which instabilities of the previous scheme are avoided by an averaging procedure. We then extend the scheme to heterogeneous networks in which (i) different neural subpopulations (e.g., excitatory and inhibitory neurons) have different cellular or connectivity parameters; (ii) the number and strength of the input connections are random (Erdős-Rényi topology) and thus different among neurons. In all heterogeneous cases, neurons are lumped in different classes each of which is represented by a single neuron in the iterative scheme; in addition, we make a Gaussian approximation of the input current to the neuron. These approximations seem to be justified over a broad range of parameters as indicated by comparison with simulation results of large recurrent networks. Our method can help to elucidate how network heterogeneity shapes the asynchronous state in recurrent neural networks.

Non-Decaying postsynaptics potentials and delayed spikes in hippocampal pyramidal neurons generated by a zero slope conductance created by the persistent Na+ current

The negative slope conductance created by the persistent sodium current (I_NaP) prolongs the decay phase of excitatory postsynaptic potentials (EPSPs). In a recent study, we demonstrated that this effect was due to an increase of the membrane time constant. When the negative slope conductance opposes completely the positive slope conductances of the other currents it creates a zero slope conductance region. In this region the membrane time constant is infinite and the decay phase of the EPSPs is virtually absent. Here we show that non-decaying EPSPs are present in CA1 hippocampal pyramidal cells in the zero slope conductance region, in the suprathreshold range of membrane potential. Na⁺ channel block with tetrodotoxin abolishes the non-decaying EPSPs. Interestingly, the non-decaying EPSPs are observed only in response to artificial excitatory postsynaptic currents (aEPSCs) of small amplitude, and not in response to aEPSCs of big amplitude. We also observed concomitantly delayed spikes with long latencies and high variability only in response to small amplitude aEPSCs. Our results showed that in CA1 pyramidal neurons I_NaP creates non-decaying EPSPs and delayed spikes in the subthreshold range of membrane potentials, which could potentiate synaptic integration of synaptic potentials coming from distal regions of the dendritic tree.

Mechanisms of Self-Sustained Oscillatory States in Hierarchical Modular Networks with Mixtures of Electrophysiological Cell Types

In a network with a mixture of different electrophysiological types of neurons linked by excitatory and inhibitory connections, temporal evolution leads through repeated epochs of intensive global activity separated by intervals with low activity level. This behavior mimics “up” and “down” states, experimentally observed in cortical tissues in absence of external stimuli. We interpret global dynamical features in terms of individual dynamics of the neurons. In particular, we observe that the crucial role both in interruption and in resumption of global activity is played by distributions of the membrane recovery variable within the network. We also demonstrate that the behavior of neurons is more influenced by their presynaptic environment in the network than by their formal types, assigned in accordance with their response to constant current.

Sustained oscillations, irregular firing, and chaotic dynamics in hierarchical modular networks with mixtures of electrophysiological cell types

The cerebral cortex exhibits neural activity even in the absence of external stimuli. This self-sustained activity is characterized by irregular firing of individual neurons and population oscillations with a broad frequency range. Questions that arise in this context, are: What are the mechanisms responsible for the existence of neuronal spiking activity in the cortex without external input? Do these mechanisms depend on the structural organization of the cortical connections? Do they depend on intrinsic characteristics of the cortical neurons? To approach the answers to these questions, we have used computer simulations of cortical network models. Our networks have hierarchical modular architecture and are composed of combinations of neuron models that reproduce the firing behavior of the five main cortical electrophysiological cell classes: regular spiking (RS), chattering (CH), intrinsically bursting (IB), low threshold spiking (LTS), and fast spiking (FS). The population of excitatory neurons is built of RS cells (always present) and either CH or IB cells. Inhibitory neurons belong to the same class, either LTS or FS. Long-lived self-sustained activity states in our network simulations display irregular single neuron firing and oscillatory activity similar to experimentally measured ones. The duration of self-sustained activity strongly depends on the initial conditions, suggesting a transient chaotic regime. Extensive analysis of the self-sustained activity states showed that their lifetime expectancy increases with the number of network modules and is favored when the network is composed of excitatory neurons of the RS and CH classes combined with inhibitory neurons of the LTS class. These results indicate that the existence and properties of the self-sustained cortical activity states depend on both the topology of the network and the neuronal mixture that comprises the network.