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.     

State-Dependent Effects of Na Channel Noise on Neuronal Burst Generation

We explore the effects of stochastic sodium (Na) channel activation on the variability and dynamics of spiking and bursting in a model neuron. The complete model segregates Hodgin-Huxley-type currents into two compartments, and undergoes applied current-dependent bifurcations between regimes of periodic bursting, chaotic bursting, and tonic spiking. Noise is added to simulate variable, finite sizes of the population of Na channels in the fast spiking compartment.During tonic firing, Na channel noise causes variability in interspike intervals (ISIs). The variance, as well as the sensitivity to noise, depend on the model's biophysical complexity. They are smallest in an isolated spiking compartment; increase significantly upon coupling to a passive compartment; and increase again when the second compartment also includes slow-acting currents. In this full model, sufficient noise can convert tonic firing into bursting.During bursting, the actions of Na channel noise are state-dependent. The higher the noise level, the greater the jitter in spike timing within bursts. The noise makes the burst durations of periodic regimes variable, while decreasing burst length duration and variance in a chaotic regime. Na channel noise blurs the sharp transitions of spike time and burst length seen at the bifurcations of the noise-free model. Close to such a bifurcation, the burst behaviors of previously periodic and chaotic regimes become essentially indistinguishable.We discuss biophysical mechanisms, dynamical interpretations and physiological implications. We suggest that noise associated with finite populations of Na channels could evoke very different effects on the intrinsic variability of spiking and bursting discharges, depending on a biological neuron's complexity and applied current-dependent state. We find that simulated channel noise in the model neuron qualitatively replicates the observed variability in burst length and interburst interval in an isolated biological bursting neuron.     

A model for the polarization of neurons by extrinsically applied electric fields.

A model is presented for the subthreshold polarization of a neuron by an applied electric field. It gives insight into how morphological features of a neuron affect its polarizability. The neuronal model consists of one or more extensively branched dendritic trees, a lumped somatic impedance, and a myelinated axon with nodes of Ranvier. The dendritic trees branch according to the 3/2-power rule of Rall, so that each tree has an equivalent cylinder representation. Equations for the membrane potential at the soma and at the nodes of Ranvier, given an arbitrary specified external potential, are derived. The solutions determine the contributions made by the dendritic tree and the axon to the net polarization at the soma. In the case of a spatially constant electric field, both the magnitude and sign of the polarization depend on simple combinations of parameters describing the neuron. One important combination is given by the ratio of internal resistances for longitudinal current spread along the dendritic tree trunk and along the axon. A second is given by the ratio between the DC space constant for the dendritic tree trunk and the distance between nodes of Ranvier in the axon. A third is given by the product of the electric field and the space constant for the trunk of the dendritic tree. When a neuron with a straight axon is subjected to a constant field, the membrane potential decays exponentially with distance from the soma. Thus, the soma seems to be a likely site for action potential initiation when the field is strong enough to elicit suprathreshold polarization. In a simple example, the way in which orientation of the various parts of the neuron affects its polarization is examined. When an axon with a bend is subjected to a spatially constant field, polarization is focused at the bend, and this is another likely site for action potential initiation.     

An ephaptic transmission model of CA3 pyramidal cells: an investigation into electric field effects

Extracellular electric fields existing throughout the living brain affect the neural coding and information processing via ephaptic transmission, independent of synapses. A two-compartment whole field effect model (WFEM) of pyramidal neurons embedded within a resistive array which simulates the extracellular medium i.e. ephapse is developed to study the effects of electric field on neuronal behaviors. We derive the two linearized filed effect models (LFEM-1 and LFEM-2) from WFEM at the stable resting state. Through matching these simplified models to the subthreshold membrane response in experiments of the resting pyramidal cells exposed to applied electric fields, we not only verify our proposed model’s validity but also found the key parameters which dominate subthreshold frequency response characteristic. Moreover, we find and give its underlying biophysical mechanism that the unsymmetrical properties of active ion channels results in the very different low-frequency response of somatic and dendritic compartments. Following, WFEM is used to investigate both direct-current (DC) and alternating-current field effect on the neural firing patterns by bifurcation analyses. We present that DC electric field could modulate neuronal excitability, with the positive field improving the excitability, the modest negative field suppressing the excitability, but interestingly, the larger negative field re-exciting the neuron back into spiking behavior. The neuron exposed to the sinusoidal electric field exhibits abundant firing patterns sensitive to the input frequency and intensity. In addition, the electrical properties of ephapse can modulate the efficacy of field effect. Our simulated results are qualitatively in line with the relevant experimental results and can explain some experimental phenomena. Furthermore, they are helpful to provide the predictions which can be tested in future experiments.     

Systems analysis of the integrative activity of the neuron (1974)

This article is aimed at revising the traditional concept of neuronal activity based on pre-eminence of transmembrane potentials and "electric summation" on the neuron surface. It presents a historical survey of the emergence of the prevailing concept on propagation of potentials along conductive structures and reveals the psychological situation that determined the transfer of this concept to dendrites and the neuronal soma. Structural and biophysical properties of the neuron which do not permit information propagation along the neuronal membrane without crude distortion are critically discussed in detail.     

神经元放电阈值的可变性及其意义

神经元能够将不同时空模式的突触输入转化为时序精确的动作电位输出，这种灵活、可靠的信息编码方式是神经集群在动态环境或特定任务下产生所需活动模式的重要基础。动作电位的产生遵循全或无规律，只有当细胞膜电压达到放电阈值时，神经元才产生动作电位。放电阈值在细胞内和细胞间具有高度可变性，具体动态依赖于刺激输入和放电历史。特别是，放电阈值对动作电位起始前的膜电压变化十分敏感，这种状态依赖性产生的生物物理根源包括Na⁺失活和K⁺激活。在绝大多数神经元中，动作电位的触发位置是轴突起始端，这个位置处的阈值可变性是决定神经元对时空输入转化规律的关键因素。但是，电生理实验中动作电位的记录位置却通常是胞体或近端树突，此处的阈值可变性高于轴突起始端，而其产生的重要根源是轴突动作电位的反向传播。基于胞体测量的相关研究显示，放电阈值动态能够增强神经元的时间编码、特征选择、增益调控和同时侦测能力。本文首先介绍放电阈值的概念及量化方法，然后详细梳理近年来国内外关于放电阈值可变性及产生根源的研究进展，在此基础上归纳总结放电阈值可变性对神经元编码的重要性，最后对未来放电阈值的研究方向进行展望。     

Regulation of somatic firing dynamics by backpropagating dendritic spikes

Pyramidal cells of the apteronotid ELL have been shown to display a characteristic mechanism of burst discharge, which has been shown to play an important role in sensory coding. This form of bursting depends on a reciprocal dendro-somatic interaction, in which discharge of a somatic spike causes a dendritic spike, which in turn contributes a dendro-somatic current flow to create a depolarizing afterpotential (DAP) in the soma. We review here our recent work showing how the timing of this DAP influences the somatic firing dynamics, and how the degree of inactivation of dendritic Na⁺ currents can cause an increased delay between somatic and dendritic spikes. This ultimately allows the DAP to become more effective at increasing the excitability of the somatic spike generating mechanism. Further, this delay between dendritic and somatic spiking can be regulated by strongly hyperpolarizing GABA_B mediated dendritic inhibition, allowing the burst dynamics to fall under synaptic regulation. In contrast, a weaker, shunting inhibition due to GABA_A mediated dendritic inhibition can regulate the dendritic spike waveform to decrease the dendro-somatic current flow and the resulting DAP. We therefore show that the qualitative behaviour of an individual cell can depend on the degree of synaptic input, and the exact timing of events across the spatial extent of the neuron. Thus, our results serve to illustrate the complex dynamics that can be observed in cells with significant dendritic arborisation, a nearly ubiquitous adaptation amongst principal neurons.     

Impact of adaptation currents on synchronization of coupled exponential integrate-and-fire neurons

The ability of spiking neurons to synchronize their activity in a network depends on the response behavior of these neurons as quantified by the phase response curve (PRC) and on coupling properties. The PRC characterizes the effects of transient inputs on spike timing and can be measured experimentally. Here we use the adaptive exponential integrate-and-fire (aEIF) neuron model to determine how subthreshold and spike-triggered slow adaptation currents shape the PRC. Based on that, we predict how synchrony and phase locked states of coupled neurons change in presence of synaptic delays and unequal coupling strengths. We find that increased subthreshold adaptation currents cause a transition of the PRC from only phase advances to phase advances and delays in response to excitatory perturbations. Increased spike-triggered adaptation currents on the other hand predominantly skew the PRC to the right. Both adaptation induced changes of the PRC are modulated by spike frequency, being more prominent at lower frequencies. Applying phase reduction theory, we show that subthreshold adaptation stabilizes synchrony for pairs of coupled excitatory neurons, while spike-triggered adaptation causes locking with a small phase difference, as long as synaptic heterogeneities are negligible. For inhibitory pairs synchrony is stable and robust against conduction delays, and adaptation can mediate bistability of in-phase and anti-phase locking. We further demonstrate that stable synchrony and bistable in/anti-phase locking of pairs carry over to synchronization and clustering of larger networks. The effects of adaptation in aEIF neurons on PRCs and network dynamics qualitatively reflect those of biophysical adaptation currents in detailed Hodgkin-Huxley-based neurons, which underscores the utility of the aEIF model for investigating the dynamical behavior of networks. Our results suggest neuronal spike frequency adaptation as a mechanism synchronizing low frequency oscillations in local excitatory networks, but indicate that inhibition rather than excitation generates coherent rhythms at higher frequencies.     

High Prevalence of Multistability of Rest States and Bursting in a Database of a Model Neuron

Flexibility in neuronal circuits has its roots in the dynamical richness of their neurons. Depending on their membrane properties single neurons can produce a plethora of activity regimes including silence, spiking and bursting. What is less appreciated is that these regimes can coexist with each other so that a transient stimulus can cause persistent change in the activity of a given neuron. Such multistability of the neuronal dynamics has been shown in a variety of neurons under different modulatory conditions. It can play either a functional role or present a substrate for dynamical diseases. We considered a database of an isolated leech heart interneuron model that can display silent, tonic spiking and bursting regimes. We analyzed only the cases of endogenous bursters producing functional half-center oscillators (HCOs). Using a one parameter (the leak conductance ()) bifurcation analysis, we extended the database to include silent regimes (stationary states) and systematically classified cases for the coexistence of silent and bursting regimes. We showed that different cases could exhibit two stable depolarized stationary states and two hyperpolarized stationary states in addition to various spiking and bursting regimes. We analyzed all cases of endogenous bursters and found that 18% of the cases were multistable, exhibiting coexistences of stationary states and bursting. Moreover, 91% of the cases exhibited multistability in some range of . We also explored HCOs built of multistable neuron cases with coexisting stationary states and a bursting regime. In 96% of cases analyzed, the HCOs resumed normal alternating bursting after one of the neurons was reset to a stationary state, proving themselves robust against this perturbation.     

Biophysical basis for three distinct dynamical mechanisms of action potential initiation

Transduction of graded synaptic input into trains of all-or-none action potentials (spikes) is a crucial step in neural coding. Hodgkin identified three classes of neurons with qualitatively different analog-to-digital transduction properties. Despite widespread use of this classification scheme, a generalizable explanation of its biophysical basis has not been described. We recorded from spinal sensory neurons representing each class and reproduced their transduction properties in a minimal model. With phase plane and bifurcation analysis, each class of excitability was shown to derive from distinct spike initiating dynamics. Excitability could be converted between all three classes by varying single parameters; moreover, several parameters, when varied one at a time, had functionally equivalent effects on excitability. From this, we conclude that the spike-initiating dynamics associated with each of Hodgkin's classes represent different outcomes in a nonlinear competition between oppositely directed, kinetically mismatched currents. Class 1 excitability occurs through a saddle node on invariant circle bifurcation when net current at perithreshold potentials is inward (depolarizing) at steady state. Class 2 excitability occurs through a Hopf bifurcation when, despite net current being outward (hyperpolarizing) at steady state, spike initiation occurs because inward current activates faster than outward current. Class 3 excitability occurs through a quasi-separatrix crossing when fast-activating inward current overpowers slow-activating outward current during a stimulus transient, although slow-activating outward current dominates during constant stimulation. Experiments confirmed that different classes of spinal lamina I neurons express the subthreshold currents predicted by our simulations and, further, that those currents are necessary for the excitability in each cell class. Thus, our results demonstrate that all three classes of excitability arise from a continuum in the direction and magnitude of subthreshold currents. Through detailed analysis of the spike-initiating process, we have explained a fundamental link between biophysical properties and qualitative differences in how neurons encode sensory input.     

Intrinsic and network rhythmogenesis in a reduced traub model for CA3 neurons

We have developed a two-compartment, eight-variable model of a CA3 pyramidal cell as a reduction of a complex 19-compartment cable model [Traub et al, 1991]. Our reduced model segregates the fast currents for sodium spiking into a proximal, soma-like, compartment and the slower calcium and calcium-mediated currents into a dendrite-like compartment. In each model periodic bursting gives way to repetitive soma spiking as somatic injected current increases. Steady dendritic stimulation can produce periodic bursting of significantly higher frequency (8–20 Hz) than can steady somatic input (<8 Hz). Bursting in our model occurs only for an intermediate range of electronic coupling conductance. It depends on the segregation of channel types and on the coupling current that flows back-and-forth between compartments. When the soma and dendrite are tightly coupled electrically, our model reduces to a single compartment and does not burst. Network simulations with our model using excitatory AMPA and NMDA synapses (without inhibition) give results similar to those obtained with the complex cable model [Traub et al, 1991; Traub et al, 1992]. Brief stimulation of a single cell in a resting network produces multiple synchronized population bursts, with fast AMPA synapses providing the dominant synchronizing mechanism. The number of bursts increases with the level of maximal NMDA conductance. For high enough maximal NMDA conductance synchronized bursting repeats indefinitely. We find that two factors can cause the cells to desynchronize when AMPA synapses are blocked: heterogeneity of properties amongst cells and intrinsically chaotic burst dynamics. But even when cells are identical, they may synchronize only approximately rather than exactly. Since our model has a limited number of parameters and variables, we have studied its cellular and network dynamics computationally with relative ease and over wide parameter ranges. Thereby, we identify some qualitative features that parallel or are distinguished from those of other neuronal systems; e.g., we discuss how bursting here differs from that in some classical models.     

Ion concentration dynamics as a mechanism for neuronal bursting

We describe a simple conductance-based model neuron that includes intra- and extracellular ion concentration dynamics and show that this model exhibits periodic bursting. The bursting arises as the fast-spiking behavior of the neuron is modulated by the slow oscillatory behavior in the ion concentration variables and vice versa. By separating these time scales and studying the bifurcation structure of the neuron, we catalog several qualitatively different bursting profiles that are strikingly similar to those seen in experimental preparations. Our work suggests that ion concentration dynamics may play an important role in modulating neuronal excitability in real biological systems.     

Dynamics and bifurcations of the adaptive exponential integrate-and-fire model

Recently, several two-dimensional spiking neuron models have been introduced, with the aim of reproducing the diversity of electrophysiological features displayed by real neurons while keeping a simple model, for simulation and analysis purposes. Among these models, the adaptive integrate-and-fire model is physiologically relevant in that its parameters can be easily related to physiological quantities. The interaction of the differential equations with the reset results in a rich and complex dynamical structure. We relate the subthreshold features of the model to the dynamical properties of the differential system and the spike patterns to the properties of a Poincaré map defined by the sequence of spikes. We find a complex bifurcation structure which has a direct interpretation in terms of spike trains. For some parameter values, spike patterns are chaotic.     

Effects of active versus passive dendritic membranes on the transfer properties of a simulated neuron

In a simulated neuron with a dendritic tree, the relative effects of active and passive dendritic membranes on transfer properties were studied. The simulations were performed by means of a digital computer. The computations calculated the changes in transmembrane voltages of many compartments over time as a function of other biophysical variables. These variables were synaptic input intensity, critical firing threshold, rate of leakage of current across the membrane, and rate of longitudinal current spread between compartments. For both passive and active dendrites, the transfer properties of the soma studied for different rates of longitudinal current spread. With low rates of current spread, graded changes in firing threshold produced correspondingly graded changes in output discharge. With high rates of current spread, the neuron became a bistable operator where spiking was enhanced if the threshold was below a certain level and suppressed if the threshold was above that level. Since alterations in firing threshold were shown to have the same effect on firing rate as alterations in synaptic input intensity, the neuron can be said to change from graded to contrast-enhancing in its response to stimuli of different intensities. The presence or absence of dendritic spiking was found to have a significant effect on the integrative properties of the simulated neuron. In particular, contrast enhancement was considerably more pronounced in neurons with passive than with active dendrites in that somatic spike rates reached a higher maximum when dendrites were passive. With active dendrites, a less intense input was needed to initiate somatic spiking than with passive dendrites because a distal dendritic spike could easily propagate by means of longitudinal current spread to the soma. Once somatic spiking was initiated, though, spike rates tended to be lower with active than with passive dendrites because the soma recovered more slowly from its post-spike refractory period if it was also influenced by refractory periods in the dendrites. The experiment of comparing neurons with active and passive dendrites was repeated at a different, higher value of synaptic input. The same differences in transfer properties between the active and passive cases emerged as before. Spiking patterns in neurons with active dendrites were also affected by the time distribution of synaptic inputs. In a previous study, inputs had been random over both space and time, varying about a predetermined mean, whereas in the present study, inputs were random over space but uniform over time. When inputs were made uniform over time, spiking became more difficult to initiate and the transition from graded to bistable response became less sharp.     

Sharpness of Spike Initiation in Neurons Explained by Compartmentalization

In cortical neurons, spikes are initiated in the axon initial segment. Seen at the soma, they appear surprisingly sharp. A standard explanation is that the current coming from the axon becomes sharp as the spike is actively backpropagated to the soma. However, sharp initiation of spikes is also seen in the input–output properties of neurons, and not only in the somatic shape of spikes; for example, cortical neurons can transmit high frequency signals. An alternative hypothesis is that Na channels cooperate, but it is not currently supported by direct experimental evidence. I propose a simple explanation based on the compartmentalization of spike initiation. When Na channels are placed in the axon, the soma acts as a current sink for the Na current. I show that there is a critical distance to the soma above which an instability occurs, so that Na channels open abruptly rather than gradually as a function of somatic voltage.     

Cholinergic modulation of neuronal spike reactions to dendritic and somatic application of excitatory acids

Acetylcholine effects on neuronal firing responses evoked by somatic or dendritic applications of excitatory amino acids were studied in slices of guinea-pig parietal cortex. Excitatory reactions initiated by dendritic activation were enhanced by acetylcholine wherever it was iontophoretically applied: either to soma or dendrites. The effect consisted in shortening spike response latencies and increasing response intensity and duration. The modified responses were recorded within 1-min interval after acetylcholine microinjections at a distance within 300 microns of the soma. Parameters of responses to somatic applications of excitatory amino acids were not significantly changed by acetylcholine. The results suggest that acetylcholine improves dendritic propagation rather than membrane excitability.     

Identification of the cable parameters in the somatic shunt model

 We show that the first five moments of the soma potential and soma current uniquely and stably determine the soma conductance and capacitance and the dendritic electrotonic length, conductance, and capacitance in the so-called somatic shunt model of the passive behavior of a neuron. We test our resulting input admittance algorithm on synthetic data and demonstrate the regularizing effect of knowledge of the ratio of soma to dendrite surface areas. Received: 9 June 1999 / Accepted in revised form: 24 January 2000