Journal of Computational Neuroscience最新文献

The mechanisms underlying the generation of hippocampal epileptic seizures and interictal events and their interactions with the sleep-wake cycle are not yet fully understood. Indeed, medial temporal lobe epilepsy is associated with hippocampal abnormalities both at the neuronal (channelopathies, impaired potassium and chloride dynamics) and network level (neuronal and axonal loss, mossy fiber sprouting), with more frequent seizures during wakefulness compared with slow-wave sleep. In this article, starting from our previous computational modeling work of the hippocampal formation based on realistic topology and synaptic connectivity, we study the role of micro- and mesoscale pathological conditions of the epileptic hippocampus in the generation and maintenance of seizure-like theta and interictal oscillations. We show, through the simulations of hippocampal activity during slow-wave sleep and wakefulness that: (i) both mossy fiber sprouting and sclerosis account for seizure-like theta activity, (ii) but they have antagonist effects (seizure-like activity occurrence increases with sprouting but decreases with sclerosis), (iii) though impaired potassium and chloride dynamics have little influence on the generation of seizure-like activity, they do play a role on the generation of interictal patterns, and (iv) seizure-like activity and fast ripples are more likely to occur during wakefulness and interictal spikes during sleep.

Place cells develop spatially-tuned receptive fields during the early stages of novel environment exploration. The generative mechanism underlying these spatially-selective responses remains largely elusive, but has been associated with theta rhythmicity. An important factor implicating the transformation of silent cells to place cells is a spatially-uniform depolarization that is mediated by a persistent sodium current. This neuronal current is modulated by extracellular calcium concentration, which, in turn, is actively controlled by astrocytes. However, there is no established relationship between the neuronal depolarization and astrocytic activity. To consider this link, we designed a bioplausible computational model of a neuronal-astrocytic network, where astrocytes induced the transient emergence of place fields in silent cells, and accelerated the plasticity-induced consolidation of place cells. Interestingly, theta oscillations emerged naturally at the network level, resulting from the astrocytic modulation of subcellular neuronal properties. Our results suggest that astrocytes participate in spatial mapping and exploration, and further highlight the computational roles of these cells in the brain.

An important function of the brain is to predict which stimulus is likely to occur based on the perceived cues. The present research studied the branching behavior of a computational network model of populations of excitatory and inhibitory neurons, both analytically and through simulations. Results show how synaptic efficacy, retroactive inhibition and short-term synaptic depression determine the dynamics of selection between different branches predicting sequences of stimuli of different probabilities. Further results show that changes in the probability of the different predictions depend on variations of neuronal gain. Such variations allow the network to optimize the probability of its predictions to changing probabilities of the sequences without changing synaptic efficacy.

Temporal filters, the ability of postsynaptic neurons to preferentially select certain presynaptic input patterns over others, have been shown to be associated with the notion of information filtering and coding of sensory inputs. Short-term plasticity (depression and facilitation; STP) has been proposed to be an important player in the generation of temporal filters. We carry out a systematic modeling, analysis and computational study to understand how characteristic postsynaptic (low-, high- and band-pass) temporal filters are generated in response to periodic presynaptic spike trains in the presence STP. We investigate how the dynamic properties of these filters depend on the interplay of a hierarchy of processes, including the arrival of the presynaptic spikes, the activation of STP, its effect on the excitatory synaptic connection efficacy, and the response of the postsynaptic cell. These mechanisms involve the interplay of a collection of time scales that operate at the single-event level (roughly, during each presynaptic interspike-interval) and control the long-term development of the temporal filters over multiple presynaptic events. These time scales are generated at the levels of the presynaptic cell (captured by the presynaptic interspike-intervals), short-term depression and facilitation, synaptic dynamics and the post-synaptic cellular currents. We develop mathematical tools to link the single-event time scales with the time scales governing the long-term dynamics of the resulting temporal filters for a relatively simple model where depression and facilitation interact at the level of the synaptic efficacy change. We extend our results and tools to account for more complex models. These include multiple STP time scales and non-periodic presynaptic inputs. The results and ideas we develop have implications for the understanding of the generation of temporal filters in complex networks for which the simple feedforward network we investigate here is a building block.

Models of synaptic plasticity have been used to better understand neural development as well as learning and memory. One prominent classic model is the Bienenstock-Cooper-Munro (BCM) model that has been particularly successful in explaining plasticity of the visual cortex. Here, in an effort to include more biophysical detail in the BCM model, we incorporate 1) feedforward inhibition, and 2) the experimental observation that large synapses are relatively harder to potentiate than weak ones, while synaptic depression is proportional to the synaptic strength. These modifications change the outcome of unsupervised plasticity under the BCM model. The amount of feed-forward inhibition adds a parameter to BCM that turns out to determine the strength of competition. In the limit of strong inhibition the learning outcome is identical to standard BCM and the neuron becomes selective to one stimulus only (winner-take-all). For smaller values of inhibition, competition is weaker and the receptive fields are less selective. However, both BCM variants can yield realistic receptive fields.

Fibromyalgia (FM) is an unsolved central pain processing disturbance. We aim to provide a unifying model for FM pathogenesis based on a loop network involving thalamocortical regions, i.e., the ventroposterior lateral thalamus (VPL), the somatosensory cortex (SC), and the thalamic reticular nucleus (TRN). The dynamics of the loop have been described by three differential equations having neuron mean firing rates as variables and containing Hill functions to model mutual interactions among the loop elements. A computational analysis conducted with MATLAB has shown a transition from monostability to bistability of the loop behavior for a weakening of GABAergic transmission between TRN and VPL. This involves the appearance of a high-firing-rate steady state, which becomes dominant and is assumed to represent pathogenic pain processing giving rise to chronic pain. Our model is consistent with a bulk of literature evidence, such as neuroimaging and pharmacological data collected on FM patients, and with correlations between FM and immunoendocrine conditions, such as stress, perimenopause, chronic inflammation, obesity, and chronic dizziness. The model suggests that critical targets for FM treatment are to be found among immunoendocrine pathways leading to GABA/glutamate imbalance having an impact on the thalamocortical system.

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.