Journal of Computational Neuroscience最新文献

We characterise cortical dynamics using partial differential equations (PDEs), analysing various connectivity patterns within the cortical sheet. This exploration yields diverse dynamics, encompassing wave equations and limit cycle activity. We presume balanced equations between excitatory and inhibitory neuronal units, reflecting the ubiquitous oscillatory patterns observed in electrophysiological measurements. Our derived dynamics comprise lowest-order wave equations (i.e., the Klein-Gordon model), limit cycle waves, higher-order PDE formulations, and transitions between limit cycles and near-zero states. Furthermore, we delve into the symmetries of the models using the Lagrangian formalism, distinguishing between continuous and discontinuous symmetries. These symmetries allow for mathematical expediency in the analysis of the model and could also be useful in studying the effect of symmetrical input from distributed cortical regions. Overall, our ability to derive multiple constraints on the fields - and predictions of the model - stems largely from the underlying assumption that the brain operates at a critical state. This assumption, in turn, drives the dynamics towards oscillatory or semi-conservative behaviour. Within this critical state, we can leverage results from the physics literature, which serve as analogues for neural fields, and implicit construct validity. Comparisons between our model predictions and electrophysiological findings from the literature - such as spectral power distribution across frequencies, wave propagation speed, epileptic seizure generation, and pattern formation over the cortical surface - demonstrate a close match. This study underscores the importance of utilizing symmetry preserving PDE formulations for further mechanistic insights into cortical activity.

Activity in layer 2/3 of the mouse primary visual cortex has been shown to depend both on visual input and the mouse's locomotion. Moreover, this activity is altered by a mismatch between the observed visual flow and the predicted visual flow from locomotion. Here, I present a simple computational model that explains previously reported recordings from layer 2/3 neurons in mice. In my model, layer 2/3 encodes the velocity difference between the estimate from visual flow and the prediction from locomotion using a neural population code. Moreover, I describe a hypothesized mechanism for how the brain may carry out computations of variables encoded in population codes. This mechanism may point to a general principle for computing any mathematical function in the brain.

The hippocampal formation is critical for episodic memory, with area Cornu Ammonis 3 (CA3) a necessary substrate for auto-associative pattern completion. Recent theoretical and experimental evidence suggests that the formation and retrieval of cell assemblies enable these functions. Yet, how cell assemblies are formed and retrieved in a full-scale spiking neural network (SNN) of CA3 that incorporates the observed diversity of neurons and connections within this circuit is not well understood. Here, we demonstrate that a data-driven SNN model quantitatively reflecting the neuron type-specific population sizes, intrinsic electrophysiology, connectivity statistics, synaptic signaling, and long-term plasticity of the mouse CA3 is capable of robust auto-association and pattern completion via cell assemblies. Our results show that a broad range of assembly sizes could successfully and systematically retrieve patterns from heavily incomplete or corrupted cues after a limited number of presentations. Furthermore, performance was robust with respect to partial overlap of assemblies through shared cells, substantially enhancing memory capacity. These novel findings provide computational evidence that the specific biological properties of the CA3 circuit produce an effective neural substrate for associative learning in the mammalian brain.

We demonstrate a model of chirp-velocity sensitivity in the inferior colliculus (IC) that retains the tuning to amplitude modulation (AM) that was established in earlier models. The mechanism of velocity sensitivity is sequence detection by octopus cells of the posteroventral cochlear nucleus, which have been proposed in physiological studies to respond preferentially to the order of arrival of cross-frequency inputs of different amplitudes. Model architecture is based on coincidence detection of a combination of excitatory and inhibitory inputs. Chirp-sensitivity of the IC output is largely controlled by the strength and timing of the chirp-sensitive octopus-cell inhibitory input. AM tuning is controlled by inhibition and excitation that are tuned to the same frequency. We present several example neurons that demonstrate the feasibility of the model in simulating realistic chirp-sensitivity and AM tuning for a wide range of characteristic frequencies. Additionally, we explore the systematic impact of varying parameters on model responses. The proposed model can be used to assess the contribution of IC chirp-velocity sensitivity to responses to complex sounds, such as speech.

Firing rate models for describing the mean-field activities of neuronal ensembles can be used effectively to study network function and dynamics, including synchronization and rhythmicity of excitatory-inhibitory populations. However, traditional Wilson-Cowan-like models, even when extended to include an explicit dynamic synaptic activation variable, are found unable to capture some dynamics such as Interneuronal Network Gamma oscillations (ING). Use of an explicit delay is helpful in simulations at the expense of complicating mathematical analysis. We resolve this issue by introducing a dynamic variable, u, that acts as an effective delay in the negative feedback loop between firing rate (r) and synaptic gating of inhibition (s). In effect, u endows synaptic activation with second order dynamics. With linear stability analysis, numerical branch-tracking and simulations, we show that our r-u-s rate model captures some key qualitative features of spiking network models for ING. We also propose an alternative formulation, a v-u-s model, in which mean membrane potential v satisfies an averaged current-balance equation. Furthermore, we extend the framework to E-I networks. With our six-variable v-u-s model, we demonstrate in firing rate models the transition from Pyramidal-Interneuronal Network Gamma (PING) to ING by increasing the external drive to the inhibitory population without adjusting synaptic weights. Having PING and ING available in a single network, without invoking synaptic blockers, is plausible and natural for explaining the emergence and transition of two different types of gamma oscillations.

We derive a next generation neural mass model of a population of quadratic-integrate-and-fire neurons, with slow adaptation, and conductance-based AMPAR, GABAR and nonlinear NMDAR synapses. We show that the Lorentzian ansatz assumption can be satisfied by introducing a piece-wise polynomial approximation of the nonlinear voltage-dependent magnesium block of NMDAR current. We study the dynamics of the resulting system for two example cases of excitatory cortical neurons and inhibitory striatal neurons. Bifurcation diagrams are presented comparing the different dynamical regimes as compared to the case of linear NMDAR currents, along with sample comparison simulation time series demonstrating different possible oscillatory solutions. The omission of the nonlinearity of NMDAR currents results in a shift in the range (and possible disappearance) of the constant high firing rate regime, along with a modulation in the amplitude and frequency power spectrum of oscillations. Moreover, nonlinear NMDAR action is seen to be state-dependent and can have opposite effects depending on the type of neurons involved and the level of input firing rate received. The presented model can serve as a computationally efficient building block in whole brain network models for investigating the differential modulation of different types of synapses under neuromodulatory influence or receptor specific malfunction.

Theta burst stimulation (TBS) is a form of repetitive transcranial magnetic stimulation (rTMS) with unknown underlying mechanisms and highly variable responses across subjects. To investigate these issues, we developed a simple computational model. Our model consisted of two neurons linked by an excitatory synapse that incorporates two mechanisms: short-term plasticity (STP) and spike-timing-dependent plasticity (STDP). We applied a variable-amplitude current through I-clamp with a TBS time pattern to the pre- and post-synaptic neurons, simulating synaptic plasticity. We analyzed the results and provided an explanation for the effects of TBS, as well as the variability of responses to it. Our findings suggest that the interplay of STP and STDP mechanisms determines the direction of plasticity, which selectively affects synapses in extended neurons and underlies functional effects. Our model describes how the timing, number, and intensity of pulses delivered to neurons during rTMS contribute to induced plasticity. This not only successfully explains the different effects of intermittent TBS (iTBS) and continuous TBS (cTBS), but also predicts the results of other protocols such as 10 Hz rTMS. We propose that the variability in responses to TBS can be attributed to the variable span of neuronal thresholds across individuals and sessions. Our model suggests a biologically plausible mechanism for the diverse responses to TBS protocols and aligns with experimental data on iTBS and cTBS outcomes. This model could potentially aid in improving TBS and rTMS protocols and customizing treatments for patients, brain areas, and brain disorders.

Replicating neural responses observed in biological systems using artificial neural networks holds significant promise in the fields of medicine and engineering. In this study, we employ ultra-fast artificial neurons based on antiferromagnetic (AFM) spin Hall oscillators to emulate the biological withdrawal reflex responsible for self-preservation against noxious stimuli, such as pain or temperature. As a result of utilizing the dynamics of AFM neurons, we are able to construct an artificial neural network that can mimic the functionality and organization of the biological neural network responsible for this reflex. The unique features of AFM neurons, such as inhibition that stems from an effective AFM inertia, allow for the creation of biologically realistic neural network components, like the interneurons in the spinal cord and antagonist motor neurons. To showcase the effectiveness of AFM neuron modeling, we conduct simulations of various scenarios that define the withdrawal reflex, including responses to both weak and strong sensory stimuli, as well as voluntary suppression of the reflex.

The computational resources of a neuromorphic network model introduced earlier were investigated in the first paper of this series. It was argued that a form of ubiquitous spontaneous local convolution enabled logical gate-like neural motifs to form into hierarchical feed-forward structures of the Hubel-Wiesel type. Here we investigate concomitant data-like structures and their dynamic rôle in memory formation, retrieval, and replay. The mechanisms give rise to the need for general inhibitory sculpting, and the simulation of the replay of episodic memories, well known in humans and recently observed in rats. Other consequences include explanations of such findings as the directional flows of neural waves in memory formation and retrieval, visual anomalies and memory deficits in schizophrenia, and the operation of GABA agonist drugs in suppressing episodic memories. We put forward the hypothesis that all neural logical operations and feature extractions are of the convolutional hierarchical type described here and in the earlier paper, and exemplified by the Hubel-Wiesel model of the visual cortex, but that in more general cases the precise geometric layering might be obscured and so far undetected.