Journal of Computational Neuroscience最新文献

Networks of spiking neurons with adaption have been shown to be able to reproduce a wide range of neural activities, including the emergent population bursting and spike synchrony that underpin brain disorders and normal function. Exact mean-field models derived from spiking neural networks are extremely valuable, as such models can be used to determine how individual neurons and the network they reside within interact to produce macroscopic network behaviours. In the paper, we derive and analyze a set of exact mean-field equations for the neural network with spike frequency adaptation. Specifically, our model is a network of Izhikevich neurons, where each neuron is modeled by a two dimensional system consisting of a quadratic integrate and fire equation plus an equation which implements spike frequency adaptation. Previous work deriving a mean-field model for this type of network, relied on the assumption of sufficiently slow dynamics of the adaptation variable. However, this approximation did not succeed in establishing an exact correspondence between the macroscopic description and the realistic neural network, especially when the adaptation time constant was not large. The challenge lies in how to achieve a closed set of mean-field equations with the inclusion of the mean-field dynamics of the adaptation variable. We address this problem by using a Lorentzian ansatz combined with the moment closure approach to arrive at a mean-field system in the thermodynamic limit. The resulting macroscopic description is capable of qualitatively and quantitatively describing the collective dynamics of the neural network, including transition between states where the individual neurons exhibit asynchronous tonic firing and synchronous bursting. We extend the approach to a network of two populations of neurons and discuss the accuracy and efficacy of our mean-field approximations by examining all assumptions that are imposed during the derivation. Numerical bifurcation analysis of our mean-field models reveals bifurcations not previously observed in the models, including a novel mechanism for emergence of bursting in the network. We anticipate our results will provide a tractable and reliable tool to investigate the underlying mechanism of brain function and dysfunction from the perspective of computational neuroscience.

Understanding neural computation on the mechanistic level requires models of neurons and neuronal networks. To analyze such models one typically has to solve coupled ordinary differential equations (ODEs), which describe the dynamics of the underlying neural system. These ODEs are solved numerically with deterministic ODE solvers that yield single solutions with either no, or only a global scalar error indicator on precision. It can therefore be challenging to estimate the effect of numerical uncertainty on quantities of interest, such as spike-times and the number of spikes. To overcome this problem, we propose to use recently developed sampling-based probabilistic solvers, which are able to quantify such numerical uncertainties. They neither require detailed insights into the kinetics of the models, nor are they difficult to implement. We show that numerical uncertainty can affect the outcome of typical neuroscience simulations, e.g. jittering spikes by milliseconds or even adding or removing individual spikes from simulations altogether, and demonstrate that probabilistic solvers reveal these numerical uncertainties with only moderate computational overhead.

Place cells are hippocampal neurons encoding the position of an animal in space. Studies of place cells are essential to understanding the processing of information by neural networks of the brain. An important characteristic of place cell spike trains is phase precession. When an animal is running through the place field, the discharges of the place cells shift from the ascending phase of the theta rhythm through the minimum to the descending phase. The role of excitatory inputs to pyramidal neurons along the Schaffer collaterals and the perforant pathway in phase precession is described, but the role of local interneurons is poorly understood. Our goal is estimating of the contribution of field CA1 interneurons to the phase precession of place cells using mathematical methods. The CA1 field is chosen because it provides the largest set of experimental data required to build and verify the model. Our simulations discover optimal parameters of the excitatory and inhibitory inputs to the pyramidal neuron so that it generates a spike train with the effect of phase precession. The uniform inhibition of pyramidal neurons best explains the effect of phase precession. Among interneurons, axo-axonal neurons make the greatest contribution to the inhibition of pyramidal cells.

In the canonical view of visual processing the neural representation of complex objects emerges as visual information is integrated through a set of convergent, hierarchically organized processing stages, ending in the primate inferior temporal lobe. It seems reasonable to infer that visual perceptual categorization requires the integrity of anterior inferior temporal cortex (area TE). Many deep neural networks (DNNs) are structured to simulate the canonical view of hierarchical processing within the visual system. However, there are some discrepancies between DNNs and the primate brain. Here we evaluated the performance of a simulated hierarchical model of vision in discriminating the same categorization problems presented to monkeys with TE removals. The model was able to simulate the performance of monkeys with TE removals in the categorization task but performed poorly when challenged with visually degraded stimuli. We conclude that further development of the model is required to match the level of visual flexibility present in the monkey visual system.

Ambiguous sensory information can lead to spontaneous alternations between perceptual states, recently shown to extend to tactile perception. The authors recently proposed a simplified form of tactile rivalry which evokes two competing percepts for a fixed difference in input amplitudes across antiphase, pulsatile stimulation of the left and right fingers. This study addresses the need for a tactile rivalry model that captures the dynamics of perceptual alternations and that incorporates the structure of the somatosensory system. The model features hierarchical processing with two stages. The first and the second stages of model could be located at the secondary somatosensory cortex (area S2), or in higher areas driven by S2. The model captures dynamical features specific to the tactile rivalry percepts and produces general characteristics of perceptual rivalry: input strength dependence of dominance times (Levelt's proposition II), short-tailed skewness of dominance time distributions and the ratio of distribution moments. The presented modelling work leads to experimentally testable predictions. The same hierarchical model could generalise to account for percept formation, competition and alternations for bistable stimuli that involve pulsatile inputs from the visual and auditory domains.

Parkinson's disease (PD) and animal models of PD feature enhanced oscillations in several frequency bands in the basal ganglia (BG). Past research has emphasized the enhancement of 13-30 Hz beta oscillations. Recently, however, oscillations in the delta band (0.5-4 Hz) have been identified as a robust predictor of dopamine loss and motor dysfunction in several BG regions in mouse models of PD. In particular, delta oscillations in the substantia nigra pars reticulata (SNr) were shown to lead oscillations in motor cortex (M1) and persist under M1 lesion, but it is not clear where these oscillations are initially generated. In this paper, we use a computational model to study how delta oscillations may arise in the SNr due to projections from the globus pallidus externa (GPe). We propose a network architecture that incorporates inhibition in SNr from oscillating GPe neurons and other SNr neurons. In our simulations, this configuration yields firing patterns in model SNr neurons that match those measured in vivo. In particular, we see the spontaneous emergence of near-antiphase active-predicting and inactive-predicting neural populations in the SNr, which persist under the inclusion of STN inputs based on experimental recordings. These results demonstrate how delta oscillations can propagate through BG nuclei despite imperfect oscillatory synchrony in the source site, narrowing down potential targets for the source of delta oscillations in PD models and giving new insight into the dynamics of SNr oscillations.

Efforts on selective neural stimulation have concentrated on segregating axons based on their size and geometry. Nonetheless, axons of the white matter or peripheral nerves may also differ in their electrophysiological properties. The primary objective of this study was to investigate the possibility of selective activation of axons by leveraging an assumed level of diversity in passive (C_m & G_leak) and active membrane properties (K_temp & G_namax). First, the stimulus waveforms with hyperpolarizing (HPP) and depolarizing pre-pulsing (DPP) were tested on selectivity in a local membrane model. The default value of membrane capacitance (C_m) was found to play a critical role in sensitivity of the chronaxie time (Chr) and rheobase (Rhe) to variations of all the four membrane parameters. Decreasing the default value of C_m, and thus the passive time constant of the membrane, amplified the sensitivity to the active parameters, K_temp and G_Namax, on Chr. The HPP waveform could selectively activate neurons even if they were diversified by membrane leakage (G_leak) only, and produced higher selectivity than DPP when parameters are varied in pairs. Selectivity measures were larger when the passive parameters (C_m & G_leak) were varied together, compared to the active parameters. Second, this novel mechanism of selectivity was investigated with non-rectangular waveforms for the stimulating phase (and HPP) in the same local membrane model. Simulation results suggest that Kt² is the most selective waveform followed by Linear and Gaussian waveforms. Traditional rectangular pulse was among the least selective of all. Finally, a compartmental axon model confirmed the main findings of the local model that Kt² is the most selective, but rank ordered the other waveforms differently. These results suggest a potentially novel mechanism of stimulation selectivity, leveraging electrophysiological variations in membrane properties, that can lead to various neural prosthetic applications.

Pyramidal neurons display a variety of active conductivities and complex morphologies that support nonlinear dendritic computation. Given growing interest in understanding the ability of pyramidal neurons to classify real-world data, in our study we applied both a detailed pyramidal neuron model and the perceptron learning algorithm to classify real-world ECG data. We used Gray coding to generate spike patterns from ECG signals as well as investigated the classification performance of the pyramidal neuron's subcellular regions. Compared with the equivalent single-layer perceptron, the pyramidal neuron performed poorly due to a weight constraint. A proposed mirroring approach for inputs, however, significantly boosted the classification performance of the neuron. We thus conclude that pyramidal neurons can classify real-world data and that the mirroring approach affects performance in a way similar to non-constrained learning.

In this paper we propose a neurogeometrical model of the behaviour of cells of the arm area of the primary motor cortex (M1). We will mathematically express as a fiber bundle the hypercolumnar organization of this cortical area, first modelled by Georgopoulos (Georgopoulos et al., 1982; Georgopoulos, 2015). On this structure, we will consider the selective tuning of M1 neurons of kinematic variables of positions and directions of movement. We will then extend this model to encode the notion of fragments introduced by Hatsopoulos et al. (2007) which describes the selectivity of neurons to movement direction varying in time. This leads to consider a higher dimensional geometrical structure where fragments are represented as integral curves. A comparison with the curves obtained through numerical simulations and experimental data will be presented. Moreover, neural activity shows coherent behaviours represented in terms of movement trajectories pointing to a specific pattern of movement decomposition Kadmon Harpaz et al. (2019). Here, we will recover this pattern through a spectral clustering algorithm in the subriemannian structure we introduced, and compare our results with the neurophysiological one of Kadmon Harpaz et al. (2019).