Abstract. We present a mean-field model of the cortex that attempts to describe the gross changes in brain electrical activity for the cycles of natural sleep. We incorporate within the model two major sleep modulatory effects: slow changes in both synaptic efficiency and in neuron resting voltage caused by the ∼90-min cycling in acetylcholine, together with even slower changes in resting voltage caused by gradual elimination during sleep of somnogens (fatigue agents) such as adenosine. We argue that the change from slow-wave sleep (SWS) to rapid-eye-movement (REM) sleep can be understood as a firstorder phase transition from a low-firing, coherent state to a high-firing, desychronized cortical state. We show that the model predictions for changes in EEG power, spectral distribution, and correlation time at the SWS-to-REM transition are consistent not only with those observed in clinical recordings of a sleeping human subject, but also with the on-cortex EEG patterns recently reported by Destexhe et al. [J. Neurosci. 19(11), (1999) [4595][4596][4597][4598][4599][4600][4601][4602][4603][4604][4605][4606][4607][4608] for the sleeping cat.

One of the grand puzzles in neuroscience is establishing the link between cognition and the disparate patterns of spontaneous and task-induced brain activity that can be measured clinically using a wide range of detection modalities such as scalp electrodes and imaging tomography. High-level brain function is not a single-neuron property, yet emerges as a cooperative phenomenon of multiply-interacting populations of neurons. Therefore a fruitful modeling approach is to picture the cerebral cortex as a continuum characterized by parameters that have been averaged over a small volume of cortical tissue. Such mean-field cortical models have been used to investigate gross patterns of brain behavior such as anesthesia, the cycles of natural sleep, memory and erasure in slow-wave sleep, and epilepsy. There is persuasive and accumulating evidence that direct gap-junction connections between inhibitory neurons promote synchronous oscillatory behavior both locally and across distances of some centimeters, but, to date, continuum models have ignored gap-junction connectivity. In this paper we employ simple mean-field arguments to derive an expression for D2, the diffusive coupling strength arising from gap-junction connections between inhibitory neurons. Using recent neurophysiological measurements reported by Fukuda [J. Neurosci. 26, 3434 (2006)], we estimate an upper limit of D2 approximately 0.6cm2. We apply a linear stability analysis to a standard mean-field cortical model, augmented with gap-junction diffusion, and find this value for the diffusive coupling strength to be close to the critical value required to destabilize the homogeneous steady state. Computer simulations demonstrate that larger values of D2 cause the noise-driven model cortex to spontaneously crystalize into random mazelike Turing structures: centimeter-scale spatial patterns in which regions of high-firing activity are intermixed with regions of low-firing activity. These structures are consistent with the spatial variations in brain activity patterns detected with the BOLD (blood oxygen-level-dependent) signal detected with magnetic resonance imaging, and may provide a natural substrate for synchronous gamma-band rhythms observed across separated EEG (electroencephalogram) electrodes.

The Wilson–Cowan neural field equations describe the dynamical behavior of a 1-D continuum of excitatory and inhibitory cortical neural aggregates, using a pair of coupled integro-differential equations. Here we use bifurcation theory and small-noise linear stochastics to study the range of a phase transitions—sudden qualitative changes in the state of a dynamical system emerging from a bifurcation—accessible to the Wilson–Cowan network. Specifically, we examine saddle-node, Hopf, Turing, and Turing–Hopf instabilities. We introduce stochasticity by adding small-amplitude spatio-temporal white noise, and analyze the resulting subthreshold fluctuations using an Ornstein–Uhlenbeck linearization. This analysis predicts divergent changes in correlation and spectral characteristics of neural activity during close approach to bifurcation from below. We validate these theoretical predictions using numerical simulations. The results demonstrate the role of noise in the emergence of critically slowed precursors in both space and time, and suggest that these early-warning signals are a universal feature of a neural system close to bifurcation. In particular, these precursor signals are likely to have neurobiological significance as early warnings of impending state change in the cortex. We support this claim with an analysis of the in vitro local field potentials recorded from slices of mouse-brain tissue. We show that in the period leading up to emergence of spontaneous seizure-like events, the mouse field potentials show a characteristic spectral focusing toward lower frequencies concomitant with a growth in fluctuation variance, consistent with critical slowing near a bifurcation point. This observation of biological criticality has clear implications regarding the feasibility of seizure prediction.

The calcium dependent plasticity (CaDP) approach to the modeling of synaptic weight change is applied using a neural field approach to realistic repetitive transcranial magnetic stimulation (rTMS) protocols. A spatially-symmetric nonlinear neural field model consisting of populations of excitatory and inhibitory neurons is used. The plasticity between excitatory cell populations is then evaluated using a CaDP approach that incorporates metaplasticity. The direction and size of the plasticity (potentiation or depression) depends on both the amplitude of stimulation and duration of the protocol. The breaks in the inhibitory theta-burst stimulation protocol are crucial to ensuring that the stimulation bursts are potentiating in nature. Tuning the parameters of a spike-timing dependent plasticity (STDP) window with a Monte Carlo approach to maximize agreement between STDP predictions and the CaDP results reproduces a realistically-shaped window with two regions of depression in agreement with the existing literature. Developing understanding of how TMS interacts with cells at a network level may be important for future investigation.

No abstract