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.     

Bifurcation control of the Morris–Lecar neuron model via a dynamic state-feedback control

In this paper, we address the control problem of bifurcations in the Morris–Lecar (ML) neuron model. With the use of a dynamic state-feedback control, two Hopf bifurcation points in the ML neuron model with Type II excitability can be relocated to new desired locations simultaneously. Also, with the proposed control law, the neuronal excitability characteristics can be transformed from Type I excitability to Type II excitability by changing the type of bifurcation, in which the neuron goes from quiescence to periodic spiking from a saddle node on an invariant circle bifurcation to a Hopf bifurcation. Simulation results are provided.     

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.     

Weak electric fields detectability in a noisy neural network

We investigate the detectability of weak electric field in a noisy neural network based on Izhikevich neuron model systematically. The neural network is composed of excitatory and inhibitory neurons with similar ratio as that in the mammalian neocortex, and the axonal conduction delays between neurons are also considered. It is found that the noise intensity can modulate the detectability of weak electric field. Stochastic resonance (SR) phenomenon induced by white noise is observed when the weak electric field is added to the network. It is interesting that SR almost disappeared when the connections between neurons are cancelled, suggesting the amplification effects of the neural coupling on the synchronization of neuronal spiking. Furthermore, the network parameters, such as the connection probability, the synaptic coupling strength, the scale of neuron population and the neuron heterogeneity, can also affect the detectability of the weak electric field. Finally, the model sensitivity is studied in detail, and results show that the neural network model has an optimal region for the detectability of weak electric field signal.     

Synchronization of neuron population subject to steady DC electric field induced by magnetic stimulation

Electric fields, which are ubiquitous in the context of neurons, are induced either by external electromagnetic fields or by endogenous electric activities. Clinical evidences point out that magnetic stimulation can induce an electric field that modulates rhythmic activity of special brain tissue, which are associated with most brain functions, including normal and pathological physiological mechanisms. Recently, the studies about the relationship between clinical treatment for psychiatric disorders and magnetic stimulation have been investigated extensively. However, further development of these techniques is limited due to the lack of understanding of the underlying mechanisms supporting the interaction between the electric field induced by magnetic stimulus and brain tissue. In this paper, the effects of steady DC electric field induced by magnetic stimulation on the coherence of an interneuronal network are investigated. Different behaviors have been observed in the network with different topologies (i.e., random and small-world network, modular network). It is found that the coherence displays a peak or a plateau when the induced electric field varies between the parameter range we defined. The coherence of the neuronal systems depends extensively on the network structure and parameters. All these parameters play a key role in determining the range for the induced electric field to synchronize network activities. The presented results could have important implications for the scientific theoretical studies regarding the effects of magnetic stimulation on human brain.     

Dynamical analysis of periodic bursting in piece-wise linear planar neuron model

A piece-wise linear planar neuron model, namely, two-dimensional McKean model with periodic drive is investigated in this paper. Periodical bursting phenomenon can be observed in the numerical simulations. By assuming the formal solutions associated with different intervals of this non-autonomous system and introducing the generalized Jacobian matrix at the non-smooth boundaries, the bifurcation mechanism for the bursting solution induced by the slowly varying periodic drive is presented. It is shown that, the discontinuous Hopf bifurcation occurring at the non-smooth boundaries, i.e., the bifurcation taking place at the thresholds of the stimulation, leads the alternation between the rest state and spiking state. That is, different oscillation modes of this non-autonomous system convert periodically due to the non-smoothness of the vector field and the slow variation of the periodic drive as well.     

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.     

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.     

Relating ion channel expression, bifurcation structure, and diverse firing patterns in a model of an identified motor neuron

Neurons show diverse firing patterns. Even neurons belonging to a single chemical or morphological class, or the same identified neuron, can display different types of electrical activity. For example, motor neuron MN5, which innervates a flight muscle of adult Drosophila, can show distinct firing patterns under the same recording conditions. We developed a two-dimensional biophysical model and show that a core complement of just two voltage-gated channels is sufficient to generate firing pattern diversity. We propose Shab and DmNa _v to be two candidate genes that could encode these core currents, and find that changes in Shab channel expression in the model can reproduce activity resembling the main firing patterns observed in MN5 recordings. We use bifurcation analysis to describe the different transitions between rest and spiking states that result from variations in Shab channel expression, exposing a connection between ion channel expression, bifurcation structure, and firing patterns in models of membrane potential dynamics.     

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.     

Spontaneously emerging direction selectivity maps in visual cortex through STDP

It is still an open question as to whether, and how, direction-selective neuronal responses in primary visual cortex are generated by feedforward thalamocortical or recurrent intracortical connections, or a combination of both. Here we present an investigation that concentrates on and, only for the sake of simplicity, restricts itself to intracortical circuits, in particular, with respect to the developmental aspects of direction selectivity through spike-timing-dependent synaptic plasticity. We show that directional responses can emerge in a recurrent network model of visual cortex with spiking neurons that integrate inputs mainly from a particular direction, thus giving rise to an asymmetrically shaped receptive field. A moving stimulus that enters the receptive field from this (preferred) direction will activate a neuron most strongly because of the increased number and/or strength of inputs from this direction and since delayed isotropic inhibition will neither overlap with, nor cancel excitation, as would be the case for other stimulus directions. It is demonstrated how direction-selective responses result from spatial asymmetries in the distribution of synaptic contacts or weights of inputs delivered to a neuron by slowly conducting intracortical axonal delay lines. By means of spike-timing-dependent synaptic plasticity with an asymmetric learning window this kind of coupling asymmetry develops naturally in a recurrent network of stochastically spiking neurons in a scenario where the neurons are activated by unidirectionally moving bar stimuli and even when only intrinsic spontaneous activity drives the learning process. We also present simulation results to show the ability of this model to produce direction preference maps similar to experimental findings     

Generation of Very Slow Neuronal Rhythms and Chaos Near the Hopf Bifurcation in Single Neuron Models

We have presented a new generation mechanism of slow spiking or repetitive discharges with extraordinarily long inter-spike intervals using the modified Hodgkin-Huxley equations (Doi and Kumagai, 2001). This generation process of slow firing is completely different from that of the well-known potassium A-current in that the steady-state current-voltage relation of the neuronal model is monotonic rather than the N-shaped one of the A-current. In this paper, we extend the previous results and show that the very slow spiking generically appears in both the three-dimensional Hodgkin-Huxley equations and the three dimensional Bonhoeffer-van der Pol (or FitzHugh-Nagumo) equations. The generation of repetitive discharges or the destabilization of the unique equilibrium point (resting potential) is a simple Hopf bifurcation. We also show that the generation of slow spiking does not depend on the stability of the Hopf bifurcation: supercritical or subcritical. The dynamics of slow spiking is investigated in detail and we demonstrate that the phenomenology of slow spiking can be categorized into two types according to the type of the corresponding bifurcation of a fast subsystem: Hopf or saddle-node bifurcation.     

In vivo measurement of cortical impedance spectrum in monkeys: implications for signal propagation

To combine insights obtained from electric field potentials (LFPs) and neuronal spiking activity (MUA) we need a better understanding of the relative spatial summation of these indices of neuronal activity. Compared to MUA, the LFP has greater spatial coherence, resulting in lower spatial specificity and lower stimulus selectivity. A differential propagation of low- and high-frequency electric signals supposedly underlies this phenomenon, which could result from cortical tissue specifically attenuating higher frequencies, i.e., from a frequency-dependent impedance spectrum. Here we directly measure the cortical impedance spectrum in vivo in monkey primary visual cortex. Our results show that impedance is independent of frequency, is homogeneous and tangentially isotropic within gray matter, and can be theoretically predicted assuming a pure-resistive conductor. We propose that the spatial summation of LFP and MUA is determined by the size of these signals' generators and the nature of neural events underlying them, rather than by biophysical properties of gray matter.     

Brain signal predictions from multi-scale networks using a linearized framework

Simulations of neural activity at different levels of detail are ubiquitous in modern neurosciences, aiding the interpretation of experimental data and underlying neural mechanisms at the level of cells and circuits. Extracellular measurements of brain signals reflecting transmembrane currents throughout the neural tissue remain commonplace. The lower frequencies (≲ 300Hz) of measured signals generally stem from synaptic activity driven by recurrent interactions among neural populations and computational models should also incorporate accurate predictions of such signals. Due to limited computational resources, large-scale neuronal network models (≳ 10⁶ neurons or so) often require reducing the level of biophysical detail and account mainly for times of action potentials (‘spikes’) or spike rates. Corresponding extracellular signal predictions have thus poorly accounted for their biophysical origin.Here we propose a computational framework for predicting spatiotemporal filter kernels for such extracellular signals stemming from synaptic activity, accounting for the biophysics of neurons, populations, and recurrent connections. Signals are obtained by convolving population spike rates by appropriate kernels for each connection pathway and summing the contributions. Our main results are that kernels derived via linearized synapse and membrane dynamics, distributions of cells, conduction delay, and volume conductor model allow for accurately capturing the spatiotemporal dynamics of ground truth extracellular signals from conductance-based multicompartment neuron networks. One particular observation is that changes in the effective membrane time constants caused by persistent synapse activation must be accounted for.The work also constitutes a major advance in computational efficiency of accurate, biophysics-based signal predictions from large-scale spike and rate-based neuron network models drastically reducing signal prediction times compared to biophysically detailed network models. This work also provides insight into how experimentally recorded low-frequency extracellular signals of neuronal activity may be approximately linearly dependent on spiking activity. A new software tool LFPykernels serves as a reference implementation of the framework.     

Computational Modeling of Seizure Dynamics Using Coupled Neuronal Networks: Factors Shaping Epileptiform Activity

Epileptic seizure dynamics span multiple scales in space and time. Understanding seizure mechanisms requires identifying the relations between seizure components within and across these scales, together with the analysis of their dynamical repertoire. Mathematical models have been developed to reproduce seizure dynamics across scales ranging from the single neuron to the neural population. In this study, we develop a network model of spiking neurons and systematically investigate the conditions, under which the network displays the emergent dynamic behaviors known from the Epileptor, which is a well-investigated abstract model of epileptic neural activity. This approach allows us to study the biophysical parameters and variables leading to epileptiform discharges at cellular and network levels. Our network model is composed of two neuronal populations, characterized by fast excitatory bursting neurons and regular spiking inhibitory neurons, embedded in a common extracellular environment represented by a slow variable. By systematically analyzing the parameter landscape offered by the simulation framework, we reproduce typical sequences of neural activity observed during status epilepticus. We find that exogenous fluctuations from extracellular environment and electro-tonic couplings play a major role in the progression of the seizure, which supports previous studies and further validates our model. We also investigate the influence of chemical synaptic coupling in the generation of spontaneous seizure-like events. Our results argue towards a temporal shift of typical spike waves with fast discharges as synaptic strengths are varied. We demonstrate that spike waves, including interictal spikes, are generated primarily by inhibitory neurons, whereas fast discharges during the wave part are due to excitatory neurons. Simulated traces are compared with in vivo experimental data from rodents at different stages of the disorder. We draw the conclusion that slow variations of global excitability, due to exogenous fluctuations from extracellular environment, and gap junction communication push the system into paroxysmal regimes. We discuss potential mechanisms underlying such machinery and the relevance of our approach, supporting previous detailed modeling studies and reflecting on the limitations of our methodology.     

An extended model for a spiking neuron class

This paper proposes an extension to the model of a spiking neuron for information processing in artificial neural networks, developing a new approach for the dynamic threshold of the integrate-and-fire neuron. This new approach invokes characteristics of biological neurons such as the behavior of chemical synapses and the receptor field. We demonstrate how such a digital model of spiking neurons can solve complex nonlinear classification with a single neuron, performing experiments for the classical XOR problem. Compared with rate-coded networks and the classical integrate-and-fire model, the trained network demonstrated faster information processing, requiring fewer neurons and shorter learning periods. The extended model validates all the logic functions of biological neurons when such functions are necessary for the proper flow of binary codes through a neural network.     

Dynamics of spontaneous activity in random networks with multiple neuron subtypes and synaptic noise

During slow-wave sleep, brain electrical activity is dominated by the slow (< 1 Hz) electroencephalogram (EEG) oscillations characterized by the periodic transitions between active (or Up) and silent (or Down) states in the membrane voltage of the cortical and thalamic neurons. Sleep slow oscillation is believed to play critical role in consolidation of recent memories. Past computational studies, based on the Hodgkin-Huxley type neuronal models, revealed possible intracellular and network mechanisms of the neuronal activity during sleep, however, they failed to explore the large-scale cortical network dynamics depending on collective behavior in the large populations of neurons. In this new study, we developed a novel class of reduced discrete time spiking neuron models for large-scale network simulations of wake and sleep dynamics. In addition to the spiking mechanism, the new model implemented nonlinearities capturing effects of the leak current, the Ca²⁺ dependent K⁺ current and the persistent Na⁺ current that were found to be critical for transitions between Up and Down states of the slow oscillation. We applied the new model to study large-scale two-dimensional cortical network activity during slow-wave sleep. Our study explained traveling wave dynamics and characteristic synchronization properties of transitions between Up and Down states of the slow oscillation as observed in vivo in recordings from cats. We further predict a critical role of synaptic noise and slow adaptive currents for spike sequence replay as found during sleep related memory consolidation.     

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.     

Multiplication and stimulus invariance in a looming-sensitive neuron.

Multiplicative operations and invariance of neuronal responses are thought to play important roles in the processing of neural information in many sensory systems. Yet the biophysical mechanisms that underlie both multiplication and invariance of neuronal responses in vivo, either at the single cell or at the network level, remain to a large extent unknown. Recent work on an identified neuron in the locust visual system (the LGMD neuron) that responds well to objects looming on a collision course towards the animal suggests that this cell represents a good model to investigate the biophysical basis of multiplication and invariance at the single neuron level. Experimental and theoretical results are consistent with multiplication being implemented by subtraction of two logarithmic terms followed by exponentiation via active membrane conductances, according to a x 1/b = exp(log(a) - log(b)). Invariance appears to be in part due to non-linear integration of synaptic inputs within the dendritic tree of this neuron.     

Firing statistics of inhibitory neuron with delayed feedback. II: Non-Markovian behavior

The instantaneous state of a neural network consists of both the degree of excitation of each neuron the network is composed of and positions of impulses in communication lines between the neurons. In neurophysiological experiments, the neuronal firing moments are registered, but not the state of communication lines. But future spiking moments depend essentially on the past positions of impulses in the lines. This suggests, that the sequence of intervals between firing moments (inter-spike intervals, ISIs) in the network could be non-Markovian.