Input-output relation and energy efficiency in the neuron with different spike threshold dynamics.

Neuron encodes and transmits information through generating sequences of output spikes, which is a high energy-consuming process. The spike is initiated when membrane depolarization reaches a threshold voltage. In many neurons, threshold is dynamic and depends on the rate of membrane depolarization (dV/dt) preceding a spike. Identifying the metabolic energy involved in neural coding and their relationship to threshold dynamic is critical to understanding neuronal function and evolution. Here, we use a modified Morris-Lecar model to investigate neuronal input-output property and energy efficiency associated with different spike threshold dynamics. We find that the neurons with dynamic threshold sensitive to dV/dt generate discontinuous frequency-current curve and type II phase response curve (PRC) through Hopf bifurcation, and weak noise could prohibit spiking when bifurcation just occurs. The threshold that is insensitive to dV/dt, instead, results in a continuous frequency-current curve, a type I PRC and a saddle-node on invariant circle bifurcation, and simultaneously weak noise cannot inhibit spiking. It is also shown that the bifurcation, frequency-current curve and PRC type associated with different threshold dynamics arise from the distinct subthreshold interactions of membrane currents. Further, we observe that the energy consumption of the neuron is related to its firing characteristics. The depolarization of spike threshold improves neuronal energy efficiency by reducing the overlap of Na(+) and K(+) currents during an action potential. The high energy efficiency is achieved at more depolarized spike threshold and high stimulus current. These results provide a fundamental biophysical connection that links spike threshold dynamics, input-output relation, energetics and spike initiation, which could contribute to uncover neural encoding mechanism.

Biophysical Insights into How Spike Threshold Depends on the Rate of Membrane Potential Depolarization in Type I and Type II Neurons.

Dynamic spike threshold plays a critical role in neuronal input-output relations. In many neurons, the threshold potential depends on the rate of membrane potential depolarization (dV/dt) preceding a spike. There are two basic classes of neural excitability, i.e., Type I and Type II, according to input-output properties. Although the dynamical and biophysical basis of their spike initiation has been established, the spike threshold dynamic for each cell type has not been well described. Here, we use a biophysical model to investigate how spike threshold depends on dV/dt in two types of neuron. It is observed that Type II spike threshold is more depolarized and more sensitive to dV/dt than Type I. With phase plane analysis, we show that each threshold dynamic arises from the different separatrix and K+ current kinetics. By analyzing subthreshold properties of membrane currents, we find the activation of hyperpolarizing current prior to spike initiation is a major factor that regulates the threshold dynamics. The outward K+ current in Type I neuron does not activate at the perithresholds, which makes its spike threshold insensitive to dV/dt. The Type II K+ current activates prior to spike initiation and there is a large net hyperpolarizing current at the perithresholds, which results in a depolarized threshold as well as a pronounced threshold dynamic. These predictions are further attested in several other functionally equivalent cases of neural excitability. Our study provides a fundamental description about how intrinsic biophysical properties contribute to the threshold dynamics in Type I and Type II neurons, which could decipher their significant functions in neural coding.

Spike-frequency adaptation of a two-compartment neuron modulated by extracellular electric fields.

Spike-frequency adaptation has been shown to play an important role in neural coding. Based on a reduced two-compartment model, here we investigate how two common adaptation currents, i.e., voltage-sensitive potassium current (I(M)) and calcium-sensitive potassium current (I(AHP)), modulate neuronal responses to extracellular electric fields. It is shown that two adaptation mechanisms lead to distinct effects on the dynamical behavior of the neuron to electric fields. These effects depend on a neuronal morphological parameter that characterizes the ratio of soma area to total membrane area and internal coupling conductance. In the case of I(AHP) current, changing the morphological parameter switches spike initiation dynamics between saddle-node on invariant cycle bifurcation and supercritical Hopf bifurcation, whereas it only switches between subcritical and supercritical Hopf bifurcations for I(M) current. Unlike the morphological parameter, internal coupling conductance is unable to alter the bifurcation scenario for both adaptation currents. We also find that the electric field threshold for triggering neuronal steady-state firing is determined by two parameters, especially by the morphological parameter. Furthermore, the neuron with I(AHP) current generates mixed-mode oscillations through the canard phenomenon for some small values of the morphological parameter. All these results suggest that morphological properties play a critical role in field-induced effects on neuronal dynamics, which could qualitatively alter the outcome of adaptation by modulating internal current between soma and dendrite. The findings are readily testable in experiments, which could help to reveal the mechanisms underlying how the neuron responds to electric field stimulus.

Neuronal spike initiation modulated by extracellular electric fields.

Based on a reduced two-compartment model, the dynamical and biophysical mechanism underlying the spike initiation of the neuron to extracellular electric fields is investigated in this paper. With stability and phase plane analysis, we first investigate in detail the dynamical properties of neuronal spike initiation induced by geometric parameter and internal coupling conductance. The geometric parameter is the ratio between soma area and total membrane area, which describes the proportion of area occupied by somatic chamber. It is found that varying it could qualitatively alter the bifurcation structures of equilibrium as well as neuronal phase portraits, which remain unchanged when varying internal coupling conductance. By analyzing the activating properties of somatic membrane currents at subthreshold potentials, we explore the relevant biophysical basis of spike initiation dynamics induced by these two parameters. It is observed that increasing geometric parameter could greatly decrease the intensity of the internal current flowing from soma to dendrite, which switches spike initiation dynamics from Hopf bifurcation to SNIC bifurcation; increasing internal coupling conductance could lead to the increase of this outward internal current, whereas the increasing range is so small that it could not qualitatively alter the spike initiation dynamics. These results highlight that neuronal geometric parameter is a crucial factor in determining the spike initiation dynamics to electric fields. The finding is useful to interpret the functional significance of neuronal biophysical properties in their encoding dynamics, which could contribute to uncovering how neuron encodes electric field signals.

Exploring how extracellular electric field modulates neuron activity through dynamical analysis of a two-compartment neuron model.

To investigate how extracellular electric field modulates neuron activity, a reduced two-compartment neuron model in the presence of electric field is introduced in this study. Depending on neuronal geometric and internal coupling parameters, the behaviors of the model have been studied extensively. The neuron model can exist in quiescent state or repetitive spiking state in response to electric field stimulus. Negative electric field mainly acts as inhibitory stimulus to the neuron, positive weak electric field could modulate spiking frequency and spike timing when the neuron is already active, and positive electric fields with sufficient intensity could directly trigger neuronal spiking in the absence of other stimulations. By bifurcation analysis, it is observed that there is saddle-node on invariant circle bifurcation, supercritical Hopf bifurcation and subcritical Hopf bifurcation appearing in the obtained two parameter bifurcation diagrams. The bifurcation structures and electric field thresholds for triggering neuron firing are determined by neuronal geometric and coupling parameters. The model predicts that the neurons with a nonsymmetric morphology between soma and dendrite, are more sensitive to electric field stimulus than those with the spherical structure. These findings suggest that neuronal geometric features play a crucial role in electric field effects on the polarization of neuronal compartments. Moreover, by determining the electric field threshold of our biophysical model, we could accurately distinguish between suprathreshold and subthreshold electric fields. Our study highlights the effects of extracellular electric field on neuronal activity from the biophysical modeling point of view. These insights into the dynamical mechanism of electric field may contribute to the investigation and development of electromagnetic therapies, and the model in our study could be further extended to a neuronal network in which the effects of electric fields on network activity may be investigated.

Effects of extremely low-frequency magnetic fields on the response of a conductance-based neuron model.

To provide insights into the modulation of neuronal activity by extremely low-frequency (ELF) magnetic field (MF), we present a conductance-based neuron model and introduce ELF sinusoidal MF as an additive voltage input. By analyzing spike times and spiking frequency, it is observed that neuron with distinct spiking patterns exhibits different response properties in the presence of MF exposure. For tonic spiking neuron, the perturbations of MF exposure on spike times is maximized at the harmonics of neuronal intrinsic spiking frequency, while it is maximized at the harmonics of bursting frequency for burst spiking neuron. As MF intensity increases, the perturbations also increase. Compared with tonic spiking, bursting dynamics are less sensitive to the perturbations of ELF MF exposure. Further, ELF MF exposure is more prone to perturb neuronal spike times relative to spiking frequency. Our finding suggests that the resonance may be one of the neural mechanisms underlying the modulatory effects of the low-intensity ELF MFs on neuronal activities. The results highlight the impacts of ELF MFs exposure on neuronal activity from the single cell level, and demonstrate various factors including ELF MF properties and neuronal spiking characteristics could determine the outcome of exposure. These insights into the mechanism of MF exposure may be relevant for the design of multi-intensity magnetic stimulus protocols, and may even contribute to the interpretation of MF effects on the central nervous systems.

Propagation of spiking regularity and double coherence resonance in feedforward networks.

We investigate the propagation of spiking regularity in noisy feedforward networks (FFNs) based on FitzHugh-Nagumo neuron model systematically. It is found that noise could modulate the transmission of firing rate and spiking regularity. Noise-induced synchronization and synfire-enhanced coherence resonance are also observed when signals propagate in noisy multilayer networks. It is interesting that double coherence resonance (DCR) with the combination of synaptic input correlation and noise intensity is finally attained after the processing layer by layer in FFNs. Furthermore, inhibitory connections also play essential roles in shaping DCR phenomena. Several properties of the neuronal network such as noise intensity, correlation of synaptic inputs, and inhibitory connections can serve as control parameters in modulating both rate coding and the order of temporal coding.