Ion Channel Density Regulates Switches between Regular and Fast Spiking in Soma but Not in Axons

The threshold firing frequency of a neuron is a characterizing feature of its dynamical behaviour, in turn determining its role in the oscillatory activity of the brain. Two main types of dynamics have been identified in brain neurons. Type 1 dynamics (regular spiking) shows a continuous relationship between frequency and stimulation current (f-I_stim) and, thus, an arbitrarily low frequency at threshold current; Type 2 (fast spiking) shows a discontinuous f-I_stim relationship and a minimum threshold frequency. In a previous study of a hippocampal neuron model, we demonstrated that its dynamics could be of both Type 1 and Type 2, depending on ion channel density. In the present study we analyse the effect of varying channel density on threshold firing frequency on two well-studied axon membranes, namely the frog myelinated axon and the squid giant axon. Moreover, we analyse the hippocampal neuron model in more detail. The models are all based on voltage-clamp studies, thus comprising experimentally measurable parameters. The choice of analysing effects of channel density modifications is due to their physiological and pharmacological relevance. We show, using bifurcation analysis, that both axon models display exclusively Type 2 dynamics, independently of ion channel density. Nevertheless, both models have a region in the channel-density plane characterized by an N-shaped steady-state current-voltage relationship (a prerequisite for Type 1 dynamics and associated with this type of dynamics in the hippocampal model). In summary, our results suggest that the hippocampal soma and the two axon membranes represent two distinct kinds of membranes; membranes with a channel-density dependent switching between Type 1 and 2 dynamics, and membranes with a channel-density independent dynamics. The difference between the two membrane types suggests functional differences, compatible with a more flexible role of the soma membrane than that of the axon membrane.     

Period doubling of calcium spike firing in a model of a Purkinje cell dendrite

Recordings from cerebellar Purkinje cell dendrites have revealed that in response to sustained current injection, the cell firing pattern can move from tonic firing of Ca²⁺ spikes to doublet firing and even to quadruplet firing or more complex firing. These firing patterns are not modified substantially if Na⁺ currents are blocked. We show that the experimental results can be viewed as a slow transition of the neuronal dynamics through a period-doubling bifurcation. To further support this conclusion and to understand the underlying mechanism that leads to doublet firing, we develop and study a simple, one-compartment model of Purkinje cell dendrite. The neuron can also exhibit quadruplet and chaotic firing patterns that are similar to the firing patterns that some of the Purkinje cells exhibit experimentally. The effects of parameters such as temperature, applied current, and potassium reversal potential in the model resemble their effects in experiments. The model dynamics involve three time scales. Ca²⁺- dependent K⁺ currents, with intermediate time scales, are responsible for the appearance of doublet firing, whereas a very slow hyperpolarizing current transfers the neuron from tonic to doublet firing. We use the fast-slow analysis to separate the effects of the three time scales. Fast-slow analysis of the neuronal dynamics, with the activation variable of the very slow, hyperpolarizing current considered as a parameter, reveals that the transitions occurs via a cascade of period-doubling bifurcations of the fast and intermediate subsystem as this slow variable increases. We carry out another analysis, with the Ca²⁺ concentration considered as a parameter, to investigate the conditions for the generation of doublet firing in systems with one effective variable with intermediate time scale, in which the rest state of the fast subsystem is terminated by a saddle-node bifurcation. We find that the scenario of period doubling in these systems can occur only if (1) the time scale of the intermediate variable (here, the decay rate of the calcium concentration) is slow enough in comparison with the interspike interval of the tonic firing at the transition but is not too slow and (2) there is a bistability of the fast subsystem of the spike-generating variables.     

Stochastically Gating Ion Channels Enable Patterned Spike Firing through Activity-Dependent Modulation of Spike Probability

The transformation of synaptic input into patterns of spike output is a fundamental operation that is determined by the particular complement of ion channels that a neuron expresses. Although it is well established that individual ion channel proteins make stochastic transitions between conducting and non-conducting states, most models of synaptic integration are deterministic, and relatively little is known about the functional consequences of interactions between stochastically gating ion channels. Here, we show that a model of stellate neurons from layer II of the medial entorhinal cortex implemented with either stochastic or deterministically gating ion channels can reproduce the resting membrane properties of stellate neurons, but only the stochastic version of the model can fully account for perithreshold membrane potential fluctuations and clustered patterns of spike output that are recorded from stellate neurons during depolarized states. We demonstrate that the stochastic model implements an example of a general mechanism for patterning of neuronal output through activity-dependent changes in the probability of spike firing. Unlike deterministic mechanisms that generate spike patterns through slow changes in the state of model parameters, this general stochastic mechanism does not require retention of information beyond the duration of a single spike and its associated afterhyperpolarization. Instead, clustered patterns of spikes emerge in the stochastic model of stellate neurons as a result of a transient increase in firing probability driven by activation of HCN channels during recovery from the spike afterhyperpolarization. Using this model, we infer conditions in which stochastic ion channel gating may influence firing patterns in vivo and predict consequences of modifications of HCN channel function for in vivo firing patterns.     

Co-existent activity patterns in inhibitory neuronal networks with short-term synaptic depression

A network of two neurons mutually coupled through inhibitory synapses that display short-term synaptic depression is considered. We show that synaptic depression expands the number of possible activity patterns that the network can display and allows for co-existence of different patterns. Specifically, the network supports different types of n-m anti-phase firing patterns, where one neuron fires n spikes followed by the other neuron firing m spikes. When maximal synaptic conductances are identical, n-n anti-phase firing patterns are obtained and there are conductance intervals over which different pairs of these solutions co-exist. The multitude of n-m anti-phase patterns and their co-existence are not found when the synapses are non-depressing. Geometric singular perturbation methods for dynamical systems are applied to the original eight-dimensional model system to derive a set of one-dimensional conditions for the existence and co-existence of different anti-phase solutions. The generality and validity of these conditions are demonstrated through numerical simulations utilizing the Hodgkin-Huxley and Morris-Lecar neuronal models.     

Convergent neuromodulation onto a network neuron can have divergent effects at the network level

Different neuromodulators often target the same ion channel. When such modulators act on different neuron types, this convergent action can enable a rhythmic network to produce distinct outputs. Less clear are the functional consequences when two neuromodulators influence the same ion channel in the same neuron. We examine the consequences of this seeming redundancy using a mathematical model of the crab gastric mill (chewing) network. This network is activated in vitro by the projection neuron MCN1, which elicits a half-center bursting oscillation between the reciprocally-inhibitory neurons LG and Int1. We focus on two neuropeptides which modulate this network, including a MCN1 neurotransmitter and the hormone crustacean cardioactive peptide (CCAP). Both activate the same voltage-gated current (I _MI ) in the LG neuron. However, I _MI-MCN1 , resulting from MCN1 released neuropeptide, has phasic dynamics in its maximal conductance due to LG presynaptic inhibition of MCN1, while I _MI-CCAP retains the same maximal conductance in both phases of the gastric mill rhythm. Separation of time scales allows us to produce a 2D model from which phase plane analysis shows that, as in the biological system, I _MI-MCN1 and I _MI-CCAP primarily influence the durations of opposing phases of this rhythm. Furthermore, I _MI-MCN1 influences the rhythmic output in a manner similar to the Int1-to-LG synapse, whereas I _MI-CCAP has an influence similar to the LG-to-Int1 synapse. These results show that distinct neuromodulators which target the same voltage-gated ion channel in the same network neuron can nevertheless produce distinct effects at the network level, providing divergent neuromodulator actions on network activity.     

The potassium A-current,low firing rates and rebound excitation in Hodgkin-Huxley models

It is widely believed, following the work of Connor and Stevens (1971,J. Physiol. Lond. 214, 31–53) that the ability to fire action potentials over a wide frequency range, especially down to very low rates, is due to the transient, potassium A-current (I _A). Using a reduction of the classical Hodgkin-Huxley model, we study the effects ofI _A on steady firing rate, especially in the near-threshold regime for the onset of firing. A minimum firing rate of zero corresponds to a homoclinic bifurcation of periodic solutions at a critical level of stimulating current. It requires that the membrane's steady-state current-voltage relation be N-shaped rather than monotonic. For experimentally based genericI _A parameters, the model does not fire at arbitrarily low rates, although it can for the more atypicalI _A parameters given by Connor and Stevens for the crab axon. When theI _A inactivation rate is slow, we find that the transient potassium current can mediate more complex firing patterns, such as periodic bursting in some parameter regimes. The number of spikes per burst increases asg _A decreases and as inactivation rate decreases. We also study howI _A affects properties of transient voltage responses, such as threshold and firing latency for anodal break excitation. We provide mathematical explanations for several of these dynamic behaviors using bifurcation theory and averaging methods.     

Distal Spike Initiation Zone Location Estimation by Morphological Simulation of Ionic Current Filtering Demonstrated in a Novel Model of an Identified Drosophila Motoneuron

Studying ion channel currents generated distally from the recording site is difficult because of artifacts caused by poor space clamp and membrane filtering. A computational model can quantify artifact parameters for correction by simulating the currents only if their exact anatomical location is known. We propose that the same artifacts that confound current recordings can help pinpoint the source of those currents by providing a signature of the neuron’s morphology. This method can improve the recording quality of currents initiated at the spike initiation zone (SIZ) that are often distal to the soma in invertebrate neurons. Drosophila being a valuable tool for characterizing ion currents, we estimated the SIZ location and quantified artifacts in an identified motoneuron, aCC/MN1-Ib, by constructing a novel multicompartmental model. Initial simulation of the measured biophysical channel properties in an isopotential Hodgkin-Huxley type neuron model partially replicated firing characteristics. Adding a second distal compartment, which contained spike-generating Na⁺ and K⁺ currents, was sufficient to simulate aCC’s in vivo activity signature. Matching this signature using a reconstructed morphology predicted that the SIZ is on aCC’s primary axon, 70 μm after the most distal dendritic branching point. From SIZ to soma, we observed and quantified selective morphological filtering of fast activating currents. Non-inactivating K⁺ currents are filtered ∼3 times less and despite their large magnitude at the soma they could be as distal as Na⁺ currents. The peak of transient component (NaT) of the voltage-activated Na⁺ current is also filtered more than the magnitude of slower persistent component (NaP), which can contribute to seizures. The corrected NaP/NaT ratio explains the previously observed discrepancy when the same channel is expressed in different cells. In summary, we used an in vivo signature to estimate ion channel location and recording artifacts, which can be applied to other neurons.     

Control of firing patterns by two transient potassium currents: leading spike,latency, bistability

Transient potassium currents distinctively affect firing properties, particularly in regulating the latency before repetitive firing. Pyramidal cells of the dorsal cochlear nucleus (DCN) have two transient potassium currents, I _Kif and I _Kis, fast and slowly inactivating, respectively, and they exhibit firing patterns with dramatically variable latencies. They show immediate repetitive firing, or only after a long latency with or without a leading spike, the so-called pauser and buildup patterns. We consider a conductance-based, ten-variable, single-compartment model for the DCN pyramidal cells (Kanold and Manis 2001). We develop and analyze a reduced three-variable integrate-and-fire model (KM-LIF) which captures the qualitative firing features. We apply dynamical systems methods to explain the underlying biophysical and mathematical mechanisms for the firing behaviors, including the characteristic firing patterns, the latency phase, the onset of repetitive firing, and some discontinuities in the timing of latency duration (e.i. first spike latency and first inter spike interval). Moreover, we obtain new insights associated with the leading spike by phase plane analysis. We further demonstrate the effects of possible heterogeneity of I _Kis. The latency before repetitive firing can be controlled to cover a large range by tuning of the relative amounts of I _Kif and I _Kis. Finally, we find for the full system robust bistability when enough I _Kis is present.     

Recovery from slow inactivation of Shab K+ channels

We have recently examined slow inactivation of Shab channels. Here we extend our characterization of Shab slow inactivation by presenting the properties of recovery from inactivation. The observations support our proposal that Shab reaches the same inactivated state either from open or closed states and suggest that closed and open state inactivation share the same mechanism. Regarding the latter, we also show that external K⁺ and TEA slow down recovery from inactivation in agreement with the hypothesis that the mechanism of Shab inactivation qualitatively differs from C-type inactivation.     

Non-monotonic effects of GABAergic synaptic inputs on neuronal firing

GABA is generally known as the principal inhibitory neurotransmitter in the nervous system, usually acting by hyperpolarizing membrane potential. However, GABAergic currents sometimes exhibit non-inhibitory effects, depending on the brain region, developmental stage or pathological condition. Here, we investigate the diverse effects of GABA on the firing rate of several single neuron models, using both analytical calculations and numerical simulations. We find that GABAergic synaptic conductance and output firing rate exhibit three qualitatively different regimes as a function of GABA reversal potential, E_GABA: monotonically decreasing for sufficiently low E_GABA (inhibitory), monotonically increasing for E_GABA above firing threshold (excitatory); and a non-monotonic region for intermediate values of E_GABA. In the non-monotonic regime, small GABA conductances have an excitatory effect while large GABA conductances show an inhibitory effect. We provide a phase diagram of different GABAergic effects as a function of GABA reversal potential and glutamate conductance. We find that noisy inputs increase the range of E_GABA for which the non-monotonic effect can be observed. We also construct a micro-circuit model of striatum to explain observed effects of GABAergic fast spiking interneurons on spiny projection neurons, including non-monotonicity, as well as the heterogeneity of the effects. Our work provides a mechanistic explanation of paradoxical effects of GABAergic synaptic inputs, with implications for understanding the effects of GABA in neural computation and development.     

A Computational Model of Neuro-Glio-Vascular Loop Interactions

We present a computational, biophysical model of neuron-astrocyte-vessel interaction. Unlike other cells, neurons convey “hunger” signals to the vascular network via an intervening layer of glial cells (astrocytes); vessels dilate and release glucose which fuels neuronal firing. Existing computational models focus on only parts of this loop (neuron→astrocyte→vessel→neuron), whereas the proposed model describes the entire loop. Neuronal firing causes release of a neurotransmitter like glutamate which triggers release of vasodilator by astrocytes via a cascade of biochemical events. Vasodilators released from astrocytic endfeet cause blood vessels to dilate and release glucose into the interstitium, part of which is taken up by the astrocyticendfeet. Glucose is converted into lactate in the astrocyte and transported into the neuron. Glucose from the interstitium and lactate (produced from glucose) influx from astrocyte are converted into ATP in the neuron. Neuronal ATP is used to drive the Na⁺/K⁺ATPase pumps, which maintain ionic gradients necessary for neuronal firing. When placed in the metabolic loop, the neuron exhibits sustained firing only when the stimulation current is more than a minimum threshold. For various combinations of initial neuronal [ATP] and external current, the neuron exhibits a variety of firing patterns including sustained firing, firing after an initial pause, burst firing etc. Neurovascular interactions under conditions of constricted vessels are also studied. Most models of cerebral circulation describe neurovascular interactions exclusively in the “forward” neuron→vessel direction. The proposed model indicates possibility of “reverse” influence also, with vasomotion rhythms influencing neural firing patterns. Another idea that emerges out of the proposed work is that brain''s computations may be more comprehensively understood in terms of neuro-glial-vascular dynamics and not in terms of neural firing alone.     

K+-Dependent Selectivity and External Ca2+ Block of Shab K+ Channels

Potassium channels allow the selective flux of K⁺ excluding the smaller, and more abundant in the extracellular solution, Na⁺ ions. Here we show that Shab is a typical K⁺ channel that excludes Na⁺ under bi-ionic, Na_o/K_i or Na_o/Rb_i, conditions. However, when internal K⁺ is replaced by Cs⁺ (Na_o/Cs_i), stable inward Na⁺ and outward Cs⁺ currents are observed. These currents show that Shab selectivity is not accounted for by protein structural elements alone, as implicit in the snug-fit model of selectivity. Additionally, here we report the block of Shab channels by external Ca²⁺ ions, and compare the effect that internal K⁺ replacement exerts on both Ca²⁺ and TEA block. Our observations indicate that Ca²⁺ blocks the channels at a site located near the external TEA binding site, and that this pore region changes conformation under conditions that allow Na⁺ permeation. In contrast, the latter ion conditions do not significantly affect the binding of quinidine to the pore central cavity. Based on our observations and the structural information derived from the NaK bacterial channel, we hypothesize that Ca²⁺ is probably coordinated by main chain carbonyls of the pore´s first K⁺-binding site.     

Effects of ion channel activity on development of dorsal root ganglion neurons

Studies of mouse dorsal root ganglion neurons in vitro demonstrate that ion channel function and regulation can influence a wide range of developmental processes. The work suggests that much as exposure to different trophic factors, the pattern of impulse activity a neuron experiences can have significant structural and functional effects during development. Studies concerning effects of ion channel activity on growth cone motility, axon fasciculation, synaptic plasticity, myelination, and intracellular signaling pathways regulating gene expression are presented in the context of changes in endogenous firing patterns during development. © 1998 John Wiley & Sons, Inc. J Neurobiol 37: 158–170, 1998     

Modeling NaV1.1/SCN1A sodium channel mutations in a microcircuit with realistic ion concentration dynamics suggests differential GABAergic mechanisms leading to hyperexcitability in epilepsy and hemiplegic migraine

Loss of function mutations of SCN1A, the gene coding for the voltage-gated sodium channel Na_V1.1, cause different types of epilepsy, whereas gain of function mutations cause sporadic and familial hemiplegic migraine type 3 (FHM-3). However, it is not clear yet how these opposite effects can induce paroxysmal pathological activities involving neuronal networks’ hyperexcitability that are specific of epilepsy (seizures) or migraine (cortical spreading depolarization, CSD). To better understand differential mechanisms leading to the initiation of these pathological activities, we used a two-neuron conductance-based model of interconnected GABAergic and pyramidal glutamatergic neurons, in which we incorporated ionic concentration dynamics in both neurons. We modeled FHM-3 mutations by increasing the persistent sodium current in the interneuron and epileptogenic mutations by decreasing the sodium conductance in the interneuron. Therefore, we studied both FHM-3 and epileptogenic mutations within the same framework, modifying only two parameters. In our model, the key effect of gain of function FHM-3 mutations is ion fluxes modification at each action potential (in particular the larger activation of voltage-gated potassium channels induced by the Na_V1.1 gain of function), and the resulting CSD-triggering extracellular potassium accumulation, which is not caused only by modifications of firing frequency. Loss of function epileptogenic mutations, on the other hand, increase GABAergic neurons’ susceptibility to depolarization block, without major modifications of firing frequency before it. Our modeling results connect qualitatively to experimental data: potassium accumulation in the case of FHM-3 mutations and facilitated depolarization block of the GABAergic neuron in the case of epileptogenic mutations. Both these effects can lead to pyramidal neuron hyperexcitability, inducing in the migraine condition depolarization block of both the GABAergic and the pyramidal neuron. Overall, our findings suggest different mechanisms of network hyperexcitability for migraine and epileptogenic Na_V1.1 mutations, implying that the modifications of firing frequency may not be the only relevant pathological mechanism.     

The Influence of Cold Temperature on Cellular Excitability of Hippocampal Networks

The hippocampus plays an important role in short term memory, learning and spatial navigation. A characteristic feature of the hippocampal region is its expression of different electrical population rhythms and activities during different brain states. Physiological fluctuations in brain temperature affect the activity patterns in hippocampus, but the underlying cellular mechanisms are poorly understood. In this work, we investigated the thermal modulation of hippocampal activity at the cellular network level. Primary cell cultures of mouse E17 hippocampus displayed robust network activation upon light cooling of the extracellular solution from baseline physiological temperatures. The activity generated was dependent on action potential firing and excitatory glutamatergic synaptic transmission. Involvement of thermosensitive channels from the transient receptor potential (TRP) family in network activation by temperature changes was ruled out, whereas pharmacological and immunochemical experiments strongly pointed towards the involvement of temperature-sensitive two-pore-domain potassium channels (K_2P), TREK/TRAAK family. In hippocampal slices we could show an increase in evoked and spontaneous synaptic activity produced by mild cooling in the physiological range that was prevented by chloroform, a K_2P channel opener. We propose that cold-induced closure of background TREK/TRAAK family channels increases the excitability of some hippocampal neurons, acting as a temperature-sensitive gate of network activation. Our findings in the hippocampus open the possibility that small temperature variations in the brain in vivo, associated with metabolism or blood flow oscillations, act as a switch mechanism of neuronal activity and determination of firing patterns through regulation of thermosensitive background potassium channel activity.     

Impact of slow K<Superscript>+</Superscript> currents on spike generation can be described by an adaptive threshold model

A neuron that is stimulated by rectangular current injections initially responds with a high firing rate, followed by a decrease in the firing rate. This phenomenon is called spike-frequency adaptation and is usually mediated by slow K⁺ currents, such as the M-type K⁺ current (I _M ) or the Ca²⁺-activated K⁺ current (I _AHP ). It is not clear how the detailed biophysical mechanisms regulate spike generation in a cortical neuron. In this study, we investigated the impact of slow K⁺ currents on spike generation mechanism by reducing a detailed conductance-based neuron model. We showed that the detailed model can be reduced to a multi-timescale adaptive threshold model, and derived the formulae that describe the relationship between slow K⁺ current parameters and reduced model parameters. Our analysis of the reduced model suggests that slow K⁺ currents have a differential effect on the noise tolerance in neural coding.     

Ion channel density and threshold dynamics of repetitive firing in a cortical neuron model

Modifying the density and distribution of ion channels in a neuron (by natural up- and down-regulation, by pharmacological intervention or by spontaneous mutations) changes its activity pattern. In the present investigation, we analyze how the impulse patterns are regulated by the density of voltage-gated channels in a model neuron, based on voltage clamp measurements of hippocampal interneurons. At least three distinct oscillatory patterns, associated with three distinct regions in the Na-K channel density plane, were found. A stability analysis showed that the different regions are characterized by saddle-node, double-orbit, and Hopf bifurcation threshold dynamics, respectively. Single strongly graded action potentials occur in an area outside the oscillatory regions, but less graded action potentials occur together with repetitive firing over a considerable range of channel densities. The presently found relationship between channel densities and oscillatory behavior may be relevance for understanding principal spiking patterns of cortical neurons (regular firing and fast spiking). It may also be of relevance for understanding the action of pharmacological compounds on brain oscillatory activity.     

The Interaction between the Drosophila EAG Potassium Channel and the Protein Kinase CaMKII Involves an Extensive Interface at the Active Site of the Kinase

The Drosophila EAG (dEAG) potassium channel is the founding member of the superfamily of KNCH channels, which are involved in cardiac repolarization, neuronal excitability and cellular proliferation. In flies, dEAG is involved in regulation of neuron firing and assembles with CaMKII to form a complex implicated in memory formation. We have characterized the interaction between the kinase domain of CaMKII and a 53-residue fragment of the dEAG channel that includes a canonical CaMKII recognition sequence. Crystal structures together with biochemical/biophysical analysis show a substrate–kinase complex with an unusually tight and extensive interface that appears to be strengthened by phosphorylation of the channel fragment. Electrophysiological recordings show that catalytically active CaMKII is required to observe active dEAG channels. A previously identified phosphorylation site in the recognition sequence is not the substrate for this crucial kinase activity, but rather contributes importantly to the tight interaction of the kinase with the channel. The available data suggest that the dEAG channel is a docking platform for the kinase and that phosphorylation of the channel's kinase recognition sequence modulates the strength of the interaction between the channel and the kinase.     

Mechanisms of Gain Control by Voltage-Gated Channels in Intrinsically-Firing Neurons

Gain modulation is a key feature of neural information processing, but underlying mechanisms remain unclear. In single neurons, gain can be measured as the slope of the current-frequency (input-output) relationship over any given range of inputs. While much work has focused on the control of basal firing rates and spike rate adaptation, gain control has been relatively unstudied. Of the limited studies on gain control, some have examined the roles of synaptic noise and passive somatic currents, but the roles of voltage-gated channels present ubiquitously in neurons have been less explored. Here, we systematically examined the relationship between gain and voltage-gated ion channels in a conductance-based, tonically-active, model neuron. Changes in expression (conductance density) of voltage-gated channels increased (Ca²⁺ channel), reduced (K⁺ channels), or produced little effect (h-type channel) on gain. We found that the gain-controlling ability of channels increased exponentially with the steepness of their activation within the dynamic voltage window (voltage range associated with firing). For depolarization-activated channels, this produced a greater channel current per action potential at higher firing rates. This allowed these channels to modulate gain by contributing to firing preferentially at states of higher excitation. A finer analysis of the current-voltage relationship during tonic firing identified narrow voltage windows at which the gain-modulating channels exerted their effects. As a proof of concept, we show that h-type channels can be tuned to modulate gain by changing the steepness of their activation within the dynamic voltage window. These results show how the impact of an ion channel on gain can be predicted from the relationship between channel kinetics and the membrane potential during firing. This is potentially relevant to understanding input-output scaling in a wide class of neurons found throughout the brain and other nervous systems.