Neural field theory of synaptic plasticity

A coupled neural field model for the standard consolidation theory

The standard consolidation theory states that short-term memories located in the hippocampus enable the consolidation of long-term memories in the neocortex. In other words, the neocortex slowly learns long-term memories with a transient support of the hippocampus that quickly learns unstable memories. However, it is not clear yet what could be the neurobiological mechanisms underlying these differences in learning rates and memory time-scales. Here, we propose a novel modeling approach of the standard consolidation theory, that focuses on its potential neurobiological mechanisms. In addition to synaptic plasticity and spike frequency adaptation, our model incorporates adult neurogenesis in the dentate gyrus as well as the difference in size between the neocortex and the hippocampus, that we associate with distance-dependent synaptic plasticity. We also take into account the interconnected spatial structure of the involved brain areas, by incorporating the above neurobiological mechanisms in a coupled neural field framework, where each area is represented by a separate neural field with intra- and inter-area connections. To our knowledge, this is the first attempt to apply neural fields to this process. Using numerical simulations and mathematical analysis, we explore the short-term and long-term dynamics of the model upon alternance of phases of hippocampal replay and retrieval cue of an external input. This external input is encodable as a memory pattern in the form of a multiple bump attractor pattern in the individual neural fields. In the model, hippocampal memory patterns become encoded first, before neocortical ones, because of the smaller distances between the bumps of the hippocampal memory patterns. As a result, retrieval of the input pattern in the neocortex at short time-scales necessitates the additional input delivered by the memory pattern of the hippocampus. Neocortical memory patterns progressively consolidate at longer times, up to a point where their retrieval does not need the support of the hippocampus anymore. At longer times, perturbation of the hippocampal neural fields by neurogenesis erases the hippocampus pattern, leading to a final state where the memory pattern is exclusively evoked in the neocortex. Therefore, the dynamics of our model successfully reproduces the main features of the standard consolidation theory. This suggests that neurogenesis in the hippocampus and distance-dependent synaptic plasticity coupled to synaptic depression and spike frequency adaptation, are indeed critical neurobiological processes in memory consolidation.

Method for quantifying arousal and consciousness in healthy states and severe brain injury via EEG-based measures of corticothalamic physiology

BACKGROUND

Characterization of normal arousal states has been achieved by fitting predictions of corticothalamic neural field theory (NFT) to electroencephalographic (EEG) spectra to yield relevant physiological parameters.

NEW METHOD

A prior fitting method is extended to distinguish conscious and unconscious states in healthy and brain injured subjects by identifying additional parameters and clusters in parameter space.

RESULTS

Fits of NFT predictions to EEG spectra are used to estimate neurophysiological parameters in healthy and brain injured subjects. Spectra are used from healthy subjects in wake and sleep and from patients with unresponsive wakefulness syndrome, in a minimally conscious state (MCS), and emerged from MCS. Subjects cluster into three groups in parameter space: conscious healthy (wake and REM), sleep, and brain injured. These are distinguished by the difference X-Y between corticocortical (X) and corticothalamic (Y) feedbacks, and by mean neural response rates α and β to incoming spikes. X-Y tracks consciousness in healthy individuals, with smaller values in wake/REM than sleep, but cannot distinguish between brain injuries. Parameters α and β differentiate deep sleep from wake/REM and brain injury.

COMPARISON WITH EXISTING METHODS

Other methods typically rely on laborious clinical assessment, manual EEG scoring, or evaluation of measures like Φ from integrated information theory, for which no efficient method exists. In contrast, the present method can be automated on a personal computer.

CONCLUSION

The method provides a means to quantify consciousness and arousal in healthy and brain injured subjects, but does not distinguish subtypes of brain injury.

Neural Field Continuum Limits and the Structure-Function Partitioning of Cognitive-Emotional Brain Networks

In The cognitive-emotional brain, Pessoa overlooks continuum effects on nonlinear brain network connectivity by eschewing neural field theories and physiologically derived constructs representative of neuronal plasticity. The absence of this content, which is so very important for understanding the dynamic structure-function embedding and partitioning of brains, diminishes the rich competitive and cooperative nature of neural networks and trivializes Pessoa's arguments, and similar arguments by other authors, on the phylogenetic and operational significance of an optimally integrated brain filled with variable-strength neural connections. Riemannian neuromanifolds, containing limit-imposing metaplastic Hebbian- and antiHebbian-type control variables, simulate scalable network behavior that is difficult to capture from the simpler graph-theoretic analysis preferred by Pessoa and other neuroscientists. Field theories suggest the partitioning and performance benefits of embedded cognitive-emotional networks that optimally evolve between exotic classical and quantum computational phases, where matrix singularities and condensations produce degenerate structure-function homogeneities unrealistic of healthy brains. Some network partitioning, as opposed to unconstrained embeddedness, is thus required for effective execution of cognitive-emotional network functions and, in our new era of neuroscience, should be considered a critical aspect of proper brain organization and operation.

Generalized neural field theory of cortical plasticity illustrated by an application to the linear phase of ocular dominance column formation in primary visual cortex

Physiologically based neural field theory (NFT) is extended to encompass cortical plasticity dynamics. An illustrative application is provided which treats the evolution of the connectivity of left- and right-eye visual stimuli to neuronal populations in the primary visual cortex (V1), and the initial, linear phase of formation of approximately one-dimensional (1D) ocular dominance columns (ODCs) that sets their transverse spatial scale. This links V1 activity, structure, and physiology within a single theory that already accounts for a range of other brain activity and connectivity phenomena, thereby enabling ODC formation and many other phenomena to be interrelated and cortical parameters to be constrained across multiple domains. The results accord with experimental ODC widths for realistic cortical parameters and are based directly on a unified description of the neuronal populations involved, their connection strengths, and the neuronal activity they support. Other key results include simple analytic approximations for ODC widths and the parameters of maximum growth rate, constraints on cortical excitatory and inhibitory gains, elucidation of the roles of specific poles of the V1 response function, and the fact that ODCs are not formed when input stimuli are fully correlated between eyes. This work provides a basis for further generalization of NFT to model other plasticity phenomena, thereby linking them to the range multiscale phenomena accounted for by NFT.

Spherical-harmonics mode decomposition of neural field equations

Large-scale neural networks can be described in the spatial continuous limit by neural field equations. For large-scale brain networks, the connectivity is typically translationally variant and imposes a large computational burden upon simulations. To reduce this burden, we take a semiquantitative approach and study the dynamics of neural fields described by a delayed integrodifferential equation. We decompose the connectivity into spatially variant and invariant contributions, which typically comprise the short- and long-range fiber systems, respectively. The neural fields are mapped on the two-dimensional spherical surface, which is choice consistent with routine mappings of cortical surfaces. Then, we perform mathematically a mode decomposition of the neural field equation into spherical harmonic basis functions. A spatial truncation of the leading orders at low wave number is consistent with the spatially coherent pattern formation of large-scale patterns observed in simulations and empirical brain imaging data and leads to a low-dimensional representation of the dynamics of the neural fields, bearing promise for an acceleration of the numerical simulations by orders of magnitude.

Next-generation neural mass and field modeling

The Wilson-Cowan population model of neural activity has greatly influenced our understanding of the mechanisms for the generation of brain rhythms and the emergence of structured brain activity. As well as the many insights that have been obtained from its mathematical analysis, it is now widely used in the computational neuroscience community for building large-scale in silico brain networks that can incorporate the increasing amount of knowledge from the Human Connectome Project. Here, we consider a neural population model in the spirit of that originally developed by Wilson and Cowan, albeit with the added advantage that it can account for the phenomena of event-related synchronization and desynchronization. This derived mean-field model provides a dynamic description for the evolution of synchrony, as measured by the Kuramoto order parameter, in a large population of quadratic integrate-and-fire model neurons. As in the original Wilson-Cowan framework, the population firing rate is at the heart of our new model; however, in a significant departure from the sigmoidal firing rate function approach, the population firing rate is now obtained as a real-valued function of the complex-valued population synchrony measure. To highlight the usefulness of this next-generation Wilson-Cowan style model, we deploy it in a number of neurobiological contexts, providing understanding of the changes in power spectra observed in electro- and magnetoencephalography neuroimaging studies of motor cortex during movement, insights into patterns of functional connectivity observed during rest and their disruption by transcranial magnetic stimulation, and to describe wave propagation across cortex.

Necessity of the sleep-wake cycle for synaptic homeostasis: system-level analysis of plasticity in the corticothalamic system

Neural field theory is used to study the system-level effects of plasticity in the corticothalamic system, where arousal states are represented parametrically by the connection strengths of the system, among other physiologically based parameters. It is found that the plasticity dynamics have no fixed points or closed cycles in the parameter space of the connection strengths, but parameter subregions exist where flows have opposite signs. Remarkably, these subregions coincide with previously identified regions that correspond to wake and slow-wave sleep, thus demonstrating state dependence of the sign of synaptic modification. We then show that a closed cycle in the parameter space is possible when the plasticity dynamics are driven by the ascending arousal system, which cycles the brain between sleep and wake to complete a closed loop that includes arcs through the opposite-flow subregions. Thus, it is concluded that both wake and sleep are necessary, and together are able to stabilize connection weights in the brain over the daily cycle, thereby providing quantitative realization of the synaptic homeostasis hypothesis.

NFTsim: Theory and Simulation of Multiscale Neural Field Dynamics

A user ready, portable, documented software package, NFTsim, is presented to facilitate numerical simulations of a wide range of brain systems using continuum neural field modeling. NFTsim enables users to simulate key aspects of brain activity at multiple scales. At the microscopic scale, it incorporates characteristics of local interactions between cells, neurotransmitter effects, synaptodendritic delays and feedbacks. At the mesoscopic scale, it incorporates information about medium to large scale axonal ranges of fibers, which are essential to model dissipative wave transmission and to produce synchronous oscillations and associated cross-correlation patterns as observed in local field potential recordings of active tissue. At the scale of the whole brain, NFTsim allows for the inclusion of long range pathways, such as thalamocortical projections, when generating macroscopic activity fields. The multiscale nature of the neural activity produced by NFTsim has the potential to enable the modeling of resulting quantities measurable via various neuroimaging techniques. In this work, we give a comprehensive description of the design and implementation of the software. Due to its modularity and flexibility, NFTsim enables the systematic study of an unlimited number of neural systems with multiple neural populations under a unified framework and allows for direct comparison with analytic and experimental predictions. The code is written in C++ and bundled with Matlab routines for a rapid quantitative analysis and visualization of the outputs. The output of NFTsim is stored in plain text file enabling users to select from a broad range of tools for offline analysis. This software enables a wide and convenient use of powerful physiologically-based neural field approaches to brain modeling. NFTsim is distributed under the Apache 2.0 license.

The development and modelling of devices and paradigms for transcranial magnetic stimulation

Magnetic stimulation is a non-invasive neurostimulation technique that can evoke action potentials and modulate neural circuits through induced electric fields. Biophysical models of magnetic stimulation have become a major driver for technological developments and the understanding of the mechanisms of magnetic neurostimulation and neuromodulation. Major technological developments involve stimulation coils with different spatial characteristics and pulse sources to control the pulse waveform. While early technological developments were the result of manual design and invention processes, there is a trend in both stimulation coil and pulse source design to mathematically optimize parameters with the help of computational models. To date, macroscopically highly realistic spatial models of the brain, as well as peripheral targets, and user-friendly software packages enable researchers and practitioners to simulate the treatment-specific and induced electric field distribution in the brains of individual subjects and patients. Neuron models further introduce the microscopic level of neural activation to understand the influence of activation dynamics in response to different pulse shapes. A number of models that were designed for online calibration to extract otherwise covert information and biomarkers from the neural system recently form a third branch of modelling.

Toxic effect of acrylamide on the development of hippocampal neurons of weaning rats

Although numerous studies have examined the neurotoxicity of acrylamide in adult animals, the effects on neuronal development in the embryonic and lactational periods are largely unknown. Thus, we examined the toxicity of acrylamide on neuronal development in the hippocampus of fetal rats during pregnancy. Sprague-Dawley rats were mated with male rats at a 1:1 ratio. Rats were administered 0, 5, 10 or 20 mg/kg acrylamide intragastrically from embryonic days 6–21. The gait scores were examined in pregnant rats in each group to analyze maternal toxicity. Eight weaning rats from each group were also euthanized on postnatal day 21 for follow-up studies. Nissl staining was used to observe histological change in the hippocampus. Immunohistochemistry was conducted to observe the condition of neurites, including dendrites and axons. Western blot assay was used to measure the expression levels of the specific nerve axon membrane protein, growth associated protein 43, and the presynaptic vesicle membrane specific protein, synaptophysin. The gait scores of gravid rats significantly increased, suggesting that acrylamide induced maternal motor dysfunction. The number of neurons, as well as expression of growth associated protein 43 and synaptophysin, was reduced with increasing acrylamide dose in postnatal day 21 weaning rats. These data suggest that acrylamide exerts dose-dependent toxic effects on the growth and development of hippocampal neurons of weaning rats.

Continuous neural network with windowed Hebbian learning

We introduce an extension of the classical neural field equation where the dynamics of the synaptic kernel satisfies the standard Hebbian type of learning (synaptic plasticity). Here, a continuous network in which changes in the weight kernel occurs in a specified time window is considered. A novelty of this model is that it admits synaptic weight decrease as well as the usual weight increase resulting from correlated activity. The resulting equation leads to a delay-type rate model for which the existence and stability of solutions such as the rest state, bumps, and traveling fronts are investigated. Some relations between the length of the time window and the bump width is derived. In addition, the effect of the delay parameter on the stability of solutions is shown. Also numerical simulations for solutions and their stability are presented.

Possible mechanisms of synaptic plasticity modulation by extremely low-frequency magnetic fields

Understanding the biological mechanisms by which extremely low-frequency (ELF, < 300 Hz) magnetic fields (MFs) interact with human brain activity is an active field of research. Such knowledge is required by international agencies providing guidelines for general public and workers exposure to ELF MFs (such as ICNIRP, the International Commission on Non-Ionizing Radiation Protection). The identification of these interaction mechanisms is extremely challenging, since the effects of ELF MF exposure need to be monitored and understood at very different spatial (from micrometers to centimeters) and temporal (from milliseconds to minutes) scales. One possibility to overcome these issues is to develop biophysical models, based on the systems of mathematical equations describing the electric or metabolic activity of the brain tissue. Biophysical models of the brain activity offer the possibility to simulate how the brain tissue interacts with ELF MFs, in order to gain new insights into experimental data, and to test novel hypotheses regarding interaction mechanisms. This paper presents novel hypotheses regarding the effects of power line (60 Hz in North America) MFs on human brain activity, with arguments from biophysical models. We suggest a hypothetic chain of events that could bridge MF exposure with detectable effects on human neurophysiology. We also suggest novel directions of research in order to reach a convergence of biophysical models of brain activity and corresponding experimental data to identify interaction mechanisms.

Numerical modelling of plasticity induced by transcranial magnetic stimulation

We use neural field theory and spike-timing dependent plasticity to make a simple but biophysically reasonable model of long-term plasticity changes in the cortex due to transcranial magnetic stimulation (TMS). We show how common TMS protocols can be captured and studied within existing neural field theory. Specifically, we look at repetitive TMS protocols such as theta burst stimulation and paired-pulse protocols. Continuous repetitive protocols result mostly in depression, but intermittent repetitive protocols in potentiation. A paired pulse protocol results in depression at short ( < ∼ 10 ms) and long ( > ∼ 100 ms) interstimulus intervals, but potentiation for mid-range intervals. The model is sensitive to the choice of neural populations that are driven by the TMS pulses, and to the parameters that describe plasticity, which may aid interpretation of the high variability in existing experimental results. Driving excitatory populations results in greater plasticity changes than driving inhibitory populations. Modelling also shows the merit in optimizing a TMS protocol based on an individual's electroencephalogram. Moreover, the model can be used to make predictions about protocols that may lead to improvements in repetitive TMS outcomes.

The Virtual Brain: a simulator of primate brain network dynamics

We present The Virtual Brain (TVB), a neuroinformatics platform for full brain network simulations using biologically realistic connectivity. This simulation environment enables the model-based inference of neurophysiological mechanisms across different brain scales that underlie the generation of macroscopic neuroimaging signals including functional MRI (fMRI), EEG and MEG. Researchers from different backgrounds can benefit from an integrative software platform including a supporting framework for data management (generation, organization, storage, integration and sharing) and a simulation core written in Python. TVB allows the reproduction and evaluation of personalized configurations of the brain by using individual subject data. This personalization facilitates an exploration of the consequences of pathological changes in the system, permitting to investigate potential ways to counteract such unfavorable processes. The architecture of TVB supports interaction with MATLAB packages, for example, the well known Brain Connectivity Toolbox. TVB can be used in a client-server configuration, such that it can be remotely accessed through the Internet thanks to its web-based HTML5, JS, and WebGL graphical user interface. TVB is also accessible as a standalone cross-platform Python library and application, and users can interact with the scientific core through the scripting interface IDLE, enabling easy modeling, development and debugging of the scientific kernel. This second interface makes TVB extensible by combining it with other libraries and modules developed by the Python scientific community. In this article, we describe the theoretical background and foundations that led to the development of TVB, the architecture and features of its major software components as well as potential neuroscience applications.

Neural mass modeling of power-line magnetic fields effects on brain activity

Neural mass models are an appropriate framework to study brain activity, combining a high degree of biological realism while being mathematically tractable. These models have been used, with a certain success, to simulate brain electric (electroencephalography, EEG) and metabolic (functional magnetic resonance imaging, fMRI) activity. However, concrete applications of neural mass models have remained limited to date. Motivated by experimental results obtained in humans, we propose in this paper a neural mass model designed to study the interaction between power-line magnetic fields (MFs) (60 Hz in North America) and brain activity. The model includes pyramidal cells; dendrite-projecting, slow GABAergic neurons; soma-projecting, fast GABAergic neurons; and glutamatergic interneurons. A simple phenomenological model of interaction between the induced electric field and neuron membranes is also considered, along with a model of post-synaptic calcium concentration and associated changes in synaptic weights Simulated EEG signals are produced in a simple protocol, both in the absence and presence of a 60 Hz MF. These results are discussed based on results obtained previously in humans. Notably, results highlight that (1) EEG alpha (8-12 Hz) power can be modulated by weak membrane depolarizations induced by the exposure; (2) the level of input noise has a significant impact on EEG power modulation; and (3) the threshold value in MF flux density resulting in a significant effect on the EEG depends on the type of neuronal populations modulated by the MF exposure. Results obtained from the model shed new light on the effects of power-line MFs on brain activity, and will provide guidance in future human experiments. This may represent a valuable contribution to international regulation agencies setting guidelines on MF values to which the general public and workers can be exposed.

Inferring network properties of cortical neurons with synaptic coupling and parameter dispersion

Computational models at different space-time scales allow us to understand the fundamental mechanisms that govern neural processes and relate uniquely these processes to neuroscience data. In this work, we propose a novel neurocomputational unit (a mesoscopic model which tell us about the interaction between local cortical nodes in a large scale neural mass model) of bursters that qualitatively captures the complex dynamics exhibited by a full network of parabolic bursting neurons. We observe that the temporal dynamics and fluctuation of mean synaptic action term exhibits a high degree of correlation with the spike/burst activity of our population. With heterogeneity in the applied drive and mean synaptic coupling derived from fast excitatory synapse approximations we observe long term behavior in our population dynamics such as partial oscillations, incoherence, and synchrony. In order to understand the origin of multistability at the population level as a function of mean synaptic coupling and heterogeneity in the firing rate threshold we employ a simple generative model for parabolic bursting recently proposed by Ghosh et al. (2009). Further, we use here a mean coupling formulated for fast spiking neurons for our analysis of generic model. Stability analysis of this mean field network allow us to identify all the relevant network states found in the detailed biophysical model. We derive here analytically several boundary solutions, a result which holds for any number of spikes per burst. These findings illustrate the role of oscillations occurring at slow time scales (bursts) on the global behavior of the network.

Neural field theory of plasticity in the cerebral cortex

A generalized timing-dependent plasticity rule is incorporated into a recent neural field theory to explore synaptic plasticity in the cerebral cortex, with both excitatory and inhibitory populations included. Analysis in the time and frequency domains reveals that cortical network behavior gives rise to a saddle-node bifurcation and resonant frequencies, including a gamma-band resonance. These system resonances constrain cortical synaptic dynamics and divide it into four classes, which depend on the type of synaptic plasticity window. Depending on the dynamical class, synaptic strengths can either have a stable fixed point, or can diverge in the absence of a separate saturation mechanism. Parameter exploration shows that time-asymmetric plasticity windows, which are signatures of spike-timing dependent plasticity, enable the richest variety of synaptic dynamics to occur. In particular, we predict a zone in parameter space which may allow brains to attain the marginal stability phenomena observed experimentally, although additional regulatory mechanisms may be required to maintain these parameters.

Complementarity of spike- and rate-based dynamics of neural systems

Relationships between spiking-neuron and rate-based approaches to the dynamics of neural assemblies are explored by analyzing a model system that can be treated by both methods, with the rate-based method further averaged over multiple neurons to give a neural-field approach. The system consists of a chain of neurons, each with simple spiking dynamics that has a known rate-based equivalent. The neurons are linked by propagating activity that is described in terms of a spatial interaction strength with temporal delays that reflect distances between neurons; feedback via a separate delay loop is also included because such loops also exist in real brains. These interactions are described using a spatiotemporal coupling function that can carry either spikes or rates to provide coupling between neurons. Numerical simulation of corresponding spike- and rate-based methods with these compatible couplings then allows direct comparison between the dynamics arising from these approaches. The rate-based dynamics can reproduce two different forms of oscillation that are present in the spike-based model: spiking rates of individual neurons and network-induced modulations of spiking rate that occur if network interactions are sufficiently strong. Depending on conditions either mode of oscillation can dominate the spike-based dynamics and in some situations, particularly when the ratio of the frequencies of these two modes is integer or half-integer, the two can both be present and interact with each other.