Journal of Computational Neuroscience最新文献

The standard protocol for studying the spiking properties of single neurons is the application of current steps while monitoring the voltage response. Although this is informative, the jump in applied current is artificial. A more physiological input is where the applied current is ramped up, reflecting chemosensory input. Unsurprisingly, neurons can respond differently to the two protocols, since ion channel activation and inactivation are affected differently. Understanding the effects of current ramps, and changes in their slopes, is facilitated by mathematical models. However, techniques for analyzing current ramps are under-developed. In this article, we demonstrate how current ramps can be analyzed in single neuron models. The primary issue is the presence of gating variables that activate on slow time scales and are therefore far from equilibrium throughout the ramp. The use of an appropriate fast-slow analysis technique allows one to fully understand the neural response to ramps of different slopes. This study is motivated by data from olfactory bulb dopamine neurons, where both fast ramp (tens of milliseconds) and slow ramp (tens of seconds) protocols are used to understand the spiking profiles of the cells. The slow ramps generate experimental bifurcation diagrams with the applied current as a bifurcation parameter, thereby establishing asymptotic spiking activity patterns. The faster ramps elicit purely transient behavior that is of relevance to most physiological inputs, which are short in duration. The two protocols together provide a broader understanding of the neuron's spiking profile and the role that slowly activating ion channels can play.

The external segment of globus pallidus (GPe) is a network of oscillatory neurons connected by inhibitory synapses. We studied the intrinsic dynamic and the response to a shared brief inhibitory stimulus in a model GPe network. Individual neurons were simulated using a phase resetting model based on measurements from mouse GPe neurons studied in slices. The neurons showed a broad heterogeneity in their firing rates and in the shapes and sizes of their phase resetting curves. Connectivity in the network was set to match experimental measurements. We generated statistically equivalent neuron heterogeneity in a small-world model, in which 99% of connections were made with near neighbors and 1% at random, and in a model with entirely random connectivity. In both networks, the resting activity was slowed and made more irregular by the local inhibition, but it did not show any periodic pattern. Cross-correlations among neuron pairs were limited to directly connected neurons. When stimulated by a shared inhibitory input, the individual neuron responses separated into two groups: one with a short and stereotyped period of inhibition followed by a transient increase in firing probability, and the other responding with a sustained inhibition. Despite differences in firing rate, the responses of the first group of neurons were of fixed duration and were synchronized across cells.

Multilevel Monte Carlo (MLMC) methods aim to speed up computation of statistics from dynamical simulations. MLMC is easy to implement and is sometimes very effective, but its efficacy may depend on the underlying dynamics. We apply MLMC to networks of spiking neurons and assess its effectiveness on prototypical models of cortical circuitry under different conditions. We find that MLMC can be very efficient for computing reliable features, i.e., features of network dynamics that are reproducible upon repeated presentation of the same external forcing. In contrast, MLMC is less effective for complex, internally generated activity. Qualitative explanations are given using concepts from random dynamical systems theory.

Recurrent neural networks of spiking neurons can exhibit long lasting and even persistent activity. Such networks are often not robust and exhibit spike and firing rate statistics that are inconsistent with experimental observations. In order to overcome this problem most previous models had to assume that recurrent connections are dominated by slower NMDA type excitatory receptors. Usually, the single neurons within these networks are very simple leaky integrate and fire neurons or other low dimensional model neurons. However real neurons are much more complex, and exhibit a plethora of active conductances which are recruited both at the sub and supra threshold regimes. Here we show that by including a small number of additional active conductances we can produce recurrent networks that are both more robust and exhibit firing-rate statistics that are more consistent with experimental results. We show that this holds both for bi-stable recurrent networks, which are thought to underlie working memory and for slowly decaying networks which might underlie the estimation of interval timing. We also show that by including these conductances, such networks can be trained to using a simple learning rule to predict temporal intervals that are an order of magnitude larger than those that can be trained in networks of leaky integrate and fire neurons.

The majority of seizures recorded in humans and experimental animal models can be described by a generic phenomenological mathematical model, the Epileptor. In this model, seizure-like events (SLEs) are driven by a slow variable and occur via saddle node (SN) and homoclinic bifurcations at seizure onset and offset, respectively. Here we investigated SLEs at the single cell level using a biophysically relevant neuron model including a slow/fast system of four equations. The two equations for the slow subsystem describe ion concentration variations and the two equations of the fast subsystem delineate the electrophysiological activities of the neuron. Using extracellular K⁺ as a slow variable, we report that SLEs with SN/homoclinic bifurcations can readily occur at the single cell level when extracellular K⁺ reaches a critical value. In patients and experimental models, seizures can also evolve into sustained ictal activity (SIA) and depolarization block (DB), activities which are also parts of the dynamic repertoire of the Epileptor. Increasing extracellular concentration of K⁺ in the model to values found during experimental status epilepticus and DB, we show that SIA and DB can also occur at the single cell level. Thus, seizures, SIA, and DB, which have been first identified as network events, can exist in a unified framework of a biophysical model at the single neuron level and exhibit similar dynamics as observed in the Epileptor.Author Summary: Epilepsy is a neurological disorder characterized by the occurrence of seizures. Seizures have been characterized in patients in experimental models at both macroscopic and microscopic scales using electrophysiological recordings. Experimental works allowed the establishment of a detailed taxonomy of seizures, which can be described by mathematical models. We can distinguish two main types of models. Phenomenological (generic) models have few parameters and variables and permit detailed dynamical studies often capturing a majority of activities observed in experimental conditions. But they also have abstract parameters, making biological interpretation difficult. Biophysical models, on the other hand, use a large number of variables and parameters due to the complexity of the biological systems they represent. Because of the multiplicity of solutions, it is difficult to extract general dynamical rules. In the present work, we integrate both approaches and reduce a detailed biophysical model to sufficiently low-dimensional equations, and thus maintaining the advantages of a generic model. We propose, at the single cell level, a unified framework of different pathological activities that are seizures, depolarization block, and sustained ictal activity.

Transcranial magnetic stimulation (TMS) is an effective method to treat neurophysiological disorders by modulating the electrical activities of neurons. Neurons can exhibit complex nonlinear behaviors underlying the external stimuli. Currently, we do not know how stimulation interacts with endogenous neural activity. In this paper, the effects of magnetic field on spiking neuron, bursting neuron and bistable neuron are studied based on the Hodgkin-Huxley (HH) neuron model. The results show that the neurons in three different states can exhibit different dynamic responses under magnetic field stimulation. The magnetic field stimulation could increase or decrease the firing frequencies of spiking neuron, bursting neuron and bistable neuron. The transitions between different firing patterns of neurons can be promoted by changing the parameters of the magnetic field. Magnetic field stimulation has a minimal impact on the firing temporal sequence sequences in bursting neuron than that in spiking neuron and bistable neuron. These results provided an insight into the impact of neuronal states on neuronal dynamic responses under brain stimulation and show that subtle changes in external conditions and stimuli can cause complex neuronal responses. This study can help us understand the state-dependent coding mechanism of neurons under electromagnetic stimulation.