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.

Adaptive control of Parkinson's state based on a nonlinear computational model with unknown parameters.

The objective here is to explore the use of adaptive input-output feedback linearization method to achieve an improved deep brain stimulation (DBS) algorithm for closed-loop control of Parkinson's state. The control law is based on a highly nonlinear computational model of Parkinson's disease (PD) with unknown parameters. The restoration of thalamic relay reliability is formulated as the desired outcome of the adaptive control methodology, and the DBS waveform is the control input. The control input is adjusted in real time according to estimates of unknown parameters as well as the feedback signal. Simulation results show that the proposed adaptive control algorithm succeeds in restoring the relay reliability of the thalamus, and at the same time achieves accurate estimation of unknown parameters. Our findings point to the potential value of adaptive control approach that could be used to regulate DBS waveform in more effective treatment of PD.

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.

Closed-loop control of the thalamocortical relay neuron's Parkinsonian state based on slow variable.

A novel closed-loop control strategy is proposed to control Parkinsonian state based on a computational model. By modeling thalamocortical relay neurons under external electric field, a slow variable feedback control is applied to restore its relay functionality. Qualitative and quantitative analysis demonstrates the performance of feedback controller based on slow variable is more efficient compared with traditional feedback control based on fast variable. These findings point to the potential value of model-based design of feedback controllers for Parkinson's disease.

Characterizing electrical signals evoked by acupuncture through complex network mapping: a new perspective on acupuncture.

The electrical signals are obtained in spinal dorsal root after different manipulations of acupuncture (MA) being taken at the 'Zusanli' point of the experiment rats. After combining the analysis of the data generated from neuronal network model and that evoked by acupuncture, it is found that features of neuronal chaotic rate time series induced by periodic stimuli can be characterized by complex network approach. The features of signals evoked by MA 'nb' 'nx' (twisting) and MA 'tb' 'tx' (lifting and thrusting) are shown to be different according to the topologies of the mapped networks. This study provides us a new perspective on the analysis of acupuncture and may give potential helps on clinical treatment.

[Influence of acupuncture of Zusanli (ST 36) on connectivity of brain functional network in healthy subjects].

OBJECTIVE: To observe the effect of acupuncture of Zusanli (ST 36) on electroencephalogram (EEG) so as to probe into its law in regulating the interconnectivity of brain functional network. METHODS: A total of 9 healthy young volunteer students (6 male, 3 female) participated in the present study. They were asked to take a dorsal position on a test-bed. EEG signals were acquired from 22 surface scalp electrodes (Fp1, Fp2, F7, F3, F2, F4, F8, A1, T3, C3, C2, C4, T4, A2, T5, P3, P2, P4, T6, O2, O1 and O2) fixed on the subject's head. Acupuncture stimulation was applied to the right Zusanli (ST 36) by manipulating the filiform needle with uniform reducing-reinforcing method and at a frequency of about 50 cycles/min for 2 min. Then the stimulation was stopped for 10 min, and repeated once again (needle-twirling frequency: 150 and 200 cycles/min), 3 times altogether. The acquired EEG data were analyzed by using coherence estimation method, average path length, average clustering coefficient, and the average degree of the articulation points (nodes) for analyzing the synchronization of EEG signals before, during and after acupuncture. RESULTS: In comparison with pre-acupuncture, the coherence amplitude values of EEG-delta (1-4 Hz) and y (31-47 Hz) waves were increased significantly after acupuncture of ST 36. No significant changes were found in the amplitude values of EEG-theta (5-8 Hz), -alpha (9-13 Hz) and-beta (14-30 Hz) waves after acupuncture stimulation. During and after acupuncture, the synchronism values of EEG-delta waves of different leads and numbers of interconnectivity between every two brain functional regions in majority of the 9 volunteers were increased clearly. In all volunteers, the degree values of all nodes except A1 and A2, the average clustering coefficients along with the increase of the threshold (r), and the average path lengths of the brain functional network of EEG-delta waves during and after acupuncture were also increased evidently (the latter two items, P < 0.05), suggesting an increase of the information exchange and functional connectivity of different brain regions. CONCLUSION: Acupuncture of Zusanli (ST 36) can increase the amplitude and synchronization of EEG-delta waves of different leads, and potentiate the functional interconnectivity of brain functional network.

Vibrational resonance in feedforward network.

This paper investigates vibrational resonance in multi-layer feedforward network (FFN) based on FitzHugh-Nagumo (FHN) neuron model. High-frequency stimuli can improve the input-output linearity of firing rates, especially for the inputs with low firing rate. For FFN network, it is found that high-frequency disturbances play important roles in enhancing the propagation of weak signal through layers. Synfire-enhanced phenomenon of signal propagation is also observed in multi-layers network, when the signal transmission is affected by high-frequency disturbances. Network connections are found to be important for the propagation of weak signal. Besides that, the characteristics of high-frequency stimuli such as heterogeneity and frequency can also modulate the propagation of neural code through layers.

Bifurcations in Morris-Lecar model exposed to DC electric field.

As an important neuron model, the Morris-Lecar (ML) equations can exhibit classes I and II excitabilities with appropriate system parameters. In this paper, the effects of external DC electric field on the neuro-computational properties of ML model are investigated using bifurcation analysis. We obtain the bifurcation diagram in two dimensional parameter space of externally applied DC current and trans-membrane potential induced by external DC electric field. The bifurcation sets partition the two dimensional parameter space about the qualitatively different behaviors of the ML model. Thus the neuron's information encodes the stimulus information, and vice versa, which is significant in neural control. Furthermore, we identify the electric field as a key parameter to control the transitions among four different excitability and spiking properties.