Morphology controls how hippocampal CA1 pyramidal neuron responds to uniform electric fields: a biophysical modeling study.

Responses of different neurons to electric field (EF) are highly variable, which depends on intrinsic properties of cell type. Here we use multi-compartmental biophysical models to investigate how morphologic features affect EF-induced responses in hippocampal CA1 pyramidal neurons. We find that the basic morphologies of neuronal elements, including diameter, length, bend, branch, and axon terminals, are all correlated with somatic depolarization through altering the current sources or sinks created by applied field. Varying them alters the EF threshold for triggering action potentials (APs), and then determines cell sensitivity to suprathreshold field. Introducing excitatory postsynaptic potential increases cell excitability and reduces morphology-dependent EF firing threshold. It is also shown that applying identical subthreshold EF results in distinct polarizations on cell membrane with different realistic morphologies. These findings shed light on the crucial role of morphologies in determining field-induced neural response from the point of view of biophysical models. The predictions are conducive to better understanding the variability in modulatory effects of EF stimulation at the cellular level, which could also aid the interpretations of how applied fields activate central nervous system neurons and affect relevant circuits.

Complexity of resting-state EEG activity in the patients with early-stage Parkinson's disease.

To investigate the abnormal brain activities in the early stage of Parkinson's disease (PD), the electroencephalogram (EEG) signals were recorded with 20 channels from non-dementia PD patients (18 patients, 8 females) and age matched healthy controls (18 subjects, 8 females) during the resting state. Two methods based on the ordinal patterns of the recorded series, i.e., permutation entropy (PE) and order index (OI), were introduced to characterize the complexity of the cortical activities for two groups. It was observed that the resting-state EEG of PD patients showed lower PE and higher OI than healthy controls, which indicated that the early-stage PD caused the reduced complexity of EEG. We further applied two methods to determine the complexity of EEG rhythms in five sub-bands. The results showed that the gamma, beta and alpha rhythms of PD patients were characterized by lower PE and higher OI, i.e., reduced complexity, than healthy subjects. No significant differences were observed in theta or delta rhythms between two groups. The findings suggested that PE and OI were promising methods to detect the abnormal changes in the dynamics of EEG signals associated with early-stage PD. Further, such changes in EEG complexity may be the early markers of the cortical or subcortical dysfunction caused by PD.

Metabolic Energy of Action Potentials Modulated by Spike Frequency Adaptation.

Spike frequency adaptation (SFA) exists in many types of neurons, which has been demonstrated to improve their abilities to process incoming information by synapses. The major carrier used by a neuron to convey synaptic signals is the sequences of action potentials (APs), which have to consume substantial metabolic energies to initiate and propagate. Here we use conductance-based models to investigate how SFA modulates the AP-related energy of neurons. The SFA is attributed to either calcium-activated K+ (IAHP) or voltage-activated K+ (IM) current. We observe that the activation of IAHP or IM increases the Na+ load used for depolarizing membrane, while produces few effects on the falling phase of AP. Then, the metabolic energy involved in Na+ current significantly increases from one AP to the next, while for K+ current it is less affected. As a consequence, the total energy cost by each AP gets larger as firing rate decays down. It is also shown that the minimum Na+ charge needed for the depolarization of each AP is unaffected during the course of SFA. This indicates that the activation of either adaptation current makes APs become less efficient to use Na+ influx for their depolarization. Further, our simulations demonstrate that the different biophysical properties of IM and IAHP result in distinct modulations of metabolic energy usage for APs. These investigations provide a fundamental link between adaptation currents and neuronal energetics, which could facilitate to interpret how SFA participates in neuronal information processing.

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.

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.

Investigation of EEG abnormalities in the early stage of Parkinson's disease.

The objective of the present study was to investigate brain activity abnormalities in the early stage of Parkinson's disease (PD). To achieve this goal, eyes-closed resting state electroencephalography (EEG) signals were recorded from 15 early-stage PD patients and 15 age-matched healthy controls. The AR Burg method and the wavelet packet entropy (WPE) method were used to characterize EEG signals in different frequency bands between the groups, respectively. In the case of the AR Burg method, an increase of relative powers in the δ- and θ-band, and a decrease of relative powers in the α- and ß-band were observed for patients compared with controls. For the WPE method, EEG signals from patients showed significant higher entropy over the global frequency domain. Furthermore, WPE in the γ-band of patients was higher than that of controls, while WPE in the δ-, Î¸-, α- and ß-band were all lower. All of these changes in EEG dynamics may represent early signs of cortical dysfunction, which have potential use as biomarkers of PD in the early stage. Our findings may be further used for early intervention and early diagnosis of PD.