Dynamics of large-scale brain activity in normal arousal states and epileptic seizures

Estimating the unmeasured membrane potential of neuronal populations from the EEG using a class of deterministic nonlinear filters

We present a model-based estimation method to reconstruct the unmeasured membrane potential of neuronal populations from a single-channel electroencephalographic (EEG) measurement. We consider a class of neural mass models that share a general structure, specifically the models by Stam et al (1999 Clin. Neurophysiol. 110 1801-13), Jansen and Rit (1995 Biol. Cybern. 73 357-66) and Wendling et al (2005 J. Clin. Neurophysiol. 22 343). Under idealized assumptions, we prove the global exponential convergence of our filter. Then, under more realistic assumptions, we investigate the robustness of our filter against model uncertainties and disturbances. Analytic proofs are provided for all results and our analyses are further illustrated via simulations.

A phenomenological model of seizure initiation suggests network structure may explain seizure frequency in idiopathic generalised epilepsy

We describe a phenomenological model of seizure initiation, consisting of a bistable switch between stable fixed point and stable limit-cycle attractors. We determine a quasi-analytic formula for the exit time problem for our model in the presence of noise. This formula--which we equate to seizure frequency--is then validated numerically, before we extend our study to explore the combined effects of noise and network structure on escape times. Here, we observe that weakly connected networks of 2, 3 and 4 nodes with equivalent first transitive components all have the same asymptotic escape times. We finally extend this work to larger networks, inferred from electroencephalographic recordings from 35 patients with idiopathic generalised epilepsies and 40 controls. Here, we find that network structure in patients correlates with smaller escape times relative to network structures from controls. These initial findings are suggestive that network structure may play an important role in seizure initiation and seizure frequency.

Corticothalamic dynamics: structure of parameter space, spectra, instabilities, and reduced model

Linear instabilities are analyzed in a physiologically based mean-field corticothalamic model and a reduced-parameter model derived from it. In both models, the stable zone corresponding to normal arousal states is bounded by a series of surfaces demarcating the onsets of instabilities. The stable zone is found to depend on delay and rate parameters, whose values have a simple relationship to the number of instabilities and dominant frequencies on the stable zone's boundary. The dominant frequencies of linear activity inside the stable zone are found to lie in clearly delineated regions, each corresponding to an instability surface on its boundary and having approximately the same dominant frequency. These regions are ordered in parameter space according to their dominant frequencies, and an instability associated with the intrathalamic loop is shown to have the highest frequency that can become unstable. This reveals an important role for the thalamus in controlling the stability and bandwidth of dynamics in the corticothalamic system as a whole. The reduced model is found to agree well with the full model in a wide region of parameter space and, thus, is a useful guide to the full model's dynamics.

Interrelating anatomical, effective, and functional brain connectivity using propagators and neural field theory

It is shown how to compute effective and functional connection matrices (eCMs and fCMs) from anatomical CMs (aCMs) and corresponding strength-of-connection matrices (sCMs) using propagator methods in which neural interactions play the role of scatterings. This analysis demonstrates how network effects dress the bare propagators (the sCMs) to yield effective propagators (the eCMs) that can be used to compute the covariances customarily used to define fCMs. The results incorporate excitatory and inhibitory connections, multiple structures and populations, asymmetries, time delays, and measurement effects. They can also be postprocessed in the same manner as experimental measurements for direct comparison with data and thereby give insights into the role of coarse-graining, thresholding, and other effects in determining the structure of CMs. The spatiotemporal results show how to generalize CMs to include time delays and how natural network modes give rise to long-range coherence at resonant frequencies. The results are demonstrated using tractable analytic cases via neural field theory of cortical and corticothalamic systems. These also demonstrate close connections between the structure of CMs and proximity to critical points of the system, highlight the importance of indirect links between brain regions and raise the possibility of imaging specific levels of indirect connectivity. Aside from the results presented explicitly here, the expression of the connections among aCMs, sCMs, eCMs, and fCMs in terms of propagators opens the way for propagator theory to be further applied to analysis of connectivity.

Modeling brain resonance phenomena using a neural mass model

Stimulation with rhythmic light flicker (photic driving) plays an important role in the diagnosis of schizophrenia, mood disorder, migraine, and epilepsy. In particular, the adjustment of spontaneous brain rhythms to the stimulus frequency (entrainment) is used to assess the functional flexibility of the brain. We aim to gain deeper understanding of the mechanisms underlying this technique and to predict the effects of stimulus frequency and intensity. For this purpose, a modified Jansen and Rit neural mass model (NMM) of a cortical circuit is used. This mean field model has been designed to strike a balance between mathematical simplicity and biological plausibility. We reproduced the entrainment phenomenon observed in EEG during a photic driving experiment. More generally, we demonstrate that such a single area model can already yield very complex dynamics, including chaos, for biologically plausible parameter ranges. We chart the entire parameter space by means of characteristic Lyapunov spectra and Kaplan-Yorke dimension as well as time series and power spectra. Rhythmic and chaotic brain states were found virtually next to each other, such that small parameter changes can give rise to switching from one to another. Strikingly, this characteristic pattern of unpredictability generated by the model was matched to the experimental data with reasonable accuracy. These findings confirm that the NMM is a useful model of brain dynamics during photic driving. In this context, it can be used to study the mechanisms of, for example, perception and epileptic seizure generation. In particular, it enabled us to make predictions regarding the stimulus amplitude in further experiments for improving the entrainment effect.

Neuroscience aims to understand the enormously complex function of the normal and diseased brain. This, in turn, is the key to explaining human behavior and to developing novel diagnostic and therapeutic procedures. We develop and use models of mean activity in a single brain area, which provide a balance between tractability and plausibility. We use such a model to explain the resonance phenomenon in a photic driving experiment, which is routinely applied in the diagnosis of various diseases including epilepsy, migraine, schizophrenia and depression. Based on the model, we make predictions on the outcome of similar resonance experiments with periodic stimulation of the patients or participants. Our results are important for researchers and clinicians analyzing brain or behavioral data following periodic input.

Reduction theories elucidate the origins of complex biological rhythms generated by interacting delay-induced oscillations

Time delay is known to induce sustained oscillations in many biological systems such as electroencephalogram (EEG) activities and gene regulations. Furthermore, interactions among delay-induced oscillations can generate complex collective rhythms, which play important functional roles. However, due to their intrinsic infinite dimensionality, theoretical analysis of interacting delay-induced oscillations has been limited. Here, we show that the two primary methods for finite-dimensional limit cycles, namely, the center manifold reduction in the vicinity of the Hopf bifurcation and the phase reduction for weak interactions, can successfully be applied to interacting infinite-dimensional delay-induced oscillations. We systematically derive the complex Ginzburg-Landau equation and the phase equation without delay for general interaction networks. Based on the reduced low-dimensional equations, we demonstrate that diffusive (linearly attractive) coupling between a pair of delay-induced oscillations can exhibit nontrivial amplitude death and multimodal phase locking. Our analysis provides unique insights into experimentally observed EEG activities such as sudden transitions among different phase-locked states and occurrence of epileptic seizures.

Quantitative modelling of sleep dynamics

Arousal is largely controlled by the ascending arousal system of the hypothalamus and brainstem, which projects to the corticothalamic system responsible for electroencephalographic (EEG) signatures of sleep. Quantitative physiologically based modelling of brainstem dynamics theory is described here, using realistic parameters, and links to EEG are outlined. Verification against a wide range of experimental data is described, including arousal dynamics under normal conditions, sleep deprivation, stimuli, stimulants and jetlag, plus key features of wake and sleep EEGs.

The sleeping brain as a complex system

'Complexity science' is a rapidly developing research direction with applications in a multitude of fields that study complex systems consisting of a number of nonlinear elements with interesting dynamics and mutual interactions. This Theme Issue 'The complexity of sleep' aims at fostering the application of complexity science to sleep research, because the brain in its different sleep stages adopts different global states that express distinct activity patterns in large and complex networks of neural circuits. This introduction discusses the contributions collected in the present Theme Issue. We highlight the potential and challenges of a complex systems approach to develop an understanding of the brain in general and the sleeping brain in particular. Basically, we focus on two topics: the complex networks approach to understand the changes in the functional connectivity of the brain during sleep, and the complex dynamics of sleep, including sleep regulation. We hope that this Theme Issue will stimulate and intensify the interdisciplinary communication to advance our understanding of the complex dynamics of the brain that underlies sleep and consciousness.

The multiple time scales of sleep dynamics as a challenge for modelling the sleeping brain

A particular property of the sleeping brain is that it exhibits dynamics on very different time scales ranging from the typical sleep oscillations such as sleep spindles and slow waves that can be observed in electroencephalogram (EEG) segments of several seconds duration over the transitions between the different sleep stages on a time scale of minutes to the dynamical processes involved in sleep regulation with typical time constants in the range of hours. There is an increasing body of work on mathematical and computational models addressing these different dynamics, however, usually considering only processes on a single time scale. In this paper, we review and present a new analysis of the dynamics of human sleep EEG at the different time scales and relate the findings to recent modelling efforts pointing out both the achievements and remaining challenges.

Characterising the dynamics of EEG waveforms as the path through parameter space of a neural mass model: application to epilepsy seizure evolution

In this paper we propose that the dynamic evolution of EEG activity during epileptic seizures may be characterised as a path through parameter space of a neural mass model, reflecting gradual changes in underlying physiological mechanisms. Previous theoretical studies have shown how boundaries in parameter space of the model (so-called bifurcations) correspond to transitions in EEG waveforms between apparently normal, spike and wave and subsequently poly-spike and wave activity. In the present manuscript, we develop a multi-objective genetic algorithm that can estimate parameters of an underlying model from clinical data recordings. A standard approach to this problem is to transform both clinical data and model output into the frequency domain and then choose parameters that minimise the difference in their respective power spectra. Instead in the present manuscript, we estimate parameters in the time domain, their choice being determined according to the best fit obtained between the model output and specific features of the observed EEG waveform. This results in an approximate path through the bifurcation plane of the model obtained from clinical data. We present comparisons of such paths through parameter space from separate seizures from an individual subject, as well as between different subjects. Differences in the path reflect subtleties of variation in the dynamics of EEG, which at present appear indistinguishable using standard clinical techniques.

Mean-field modeling of thalamocortical dynamics and a model-driven approach to EEG analysis

Higher brain function depends on task-dependent information flow between cortical regions. Converging lines of evidence suggest that interactions between cortical regions and the central thalamus play a key role in establishing the dynamic patterns of functional connectivity that normally support these processes. In patients with chronic disturbances of cognitive function due to severe brain injury, dysfunction of this circuitry likely plays a crucial role in pathogenesis. However, assaying thalamocortical interactions is challenging even in healthy subjects and more so in severely impaired patients. To approach this problem, we apply a dynamical-systems approach to motivate an analysis of the electroencephalogram (EEG). We begin with a model for a single thalamocortical module [Robinson PA, Rennie CJ, Rowe DL (2002) Phys Rev E Stat Nonlin Soft Matter Phys 65:041924; Robinson PA, Rennie CJ, Wright JJ, Bourke PD (1998) Phys Rev E Stat Nonlin Soft Matter Phys 58:3557-3571]. When two such modules interact via shared thalamic inhibition, multistable behavior emerges; each mode is characterized by a different pattern of coherence between cortical regions. This observation suggests that changing patterns of cortical coherence are a hallmark of normal thalamocortical dynamics. In a preliminary study, we test this idea by analyzing the EEG of a patient with chronic brain injury, who has a marked improvement in behavior and frontal brain metabolism in response to zolpidem. The analysis shows that following zolpidem administration, changing patterns of coherence are identified between the frontal lobes and between frontal and distant brain regions. These observations support the role of the central thalamus in the organization of patterns of cortical interactions and suggest how indexes of thalamocortical dynamics can be extracted from the EEG.

Dynamic causal modeling with neural fields

The aim of this paper is twofold: first, to introduce a neural field model motivated by a well-known neural mass model; second, to show how one can estimate model parameters pertaining to spatial (anatomical) properties of neuronal sources based on EEG or LFP spectra using Bayesian inference. Specifically, we consider neural field models of cortical activity as generative models in the context of dynamic causal modeling (DCM). This paper considers the simplest case of a single cortical source modeled by the spatiotemporal dynamics of hidden neuronal states on a bounded cortical surface or manifold. We build this model using multiple layers, corresponding to cortical lamina in the real cortical manifold. These layers correspond to the populations considered in classical (Jansen and Rit) neural mass models. This allows us to formulate a neural field model that can be reduced to a neural mass model using appropriate constraints on its spatial parameters. In turn, this enables one to compare and contrast the predicted responses from equivalent neural field and mass models respectively. We pursue this using empirical LFP data from a single electrode to show that the parameters controlling the spatial dynamics of cortical activity can be recovered, using DCM, even in the absence of explicit spatial information in observed data.

Transition to absence seizures and the role of GABA(A) receptors

Absence seizures appear to be initiated in a putative cortical ‘initiation site’ by the expression of medium-amplitude 5–9 Hz oscillations, which may in part be due to a decreased phasic GABA_A receptor function. These oscillations rapidly spread to other cortical areas and to the thalamus, leading to fully developed generalized spike and wave discharges. In thalamocortical neurons of genetic models, phasic GABA_A inhibition is either unchanged or increased, whereas tonic GABA_A inhibition is increased both in genetic and pharmacological models. This enhanced tonic inhibition is required for absence seizure generation, and in genetic models it results from a malfunction in the astrocytic GABA transporter GAT-1. Contradictory results from inbred and transgenic animals still do not allow us to draw firm conclusions on changes in phasic GABA_A inhibition in the GABAergic neurons of the nucleus reticularis thalami. Mathematical modelling may enhance our understanding of these competing hypotheses, by permitting investigations of their mechanistic aspects, hence enabling a greater understanding of the processes underlying seizure generation and evolution.

Neural field theory of synaptic plasticity

Plasticity is crucial to neural development, learning, and memory. In the common in vivo situation where postsynaptic neural activity results from multiple presynaptic inputs, it is shown that a widely used class of correlation-dependent and spike-timing dependent plasticity rules can be written in a form that can be incorporated into neural field theory, which enables their system-level dynamics to be investigated. It is shown that the resulting plasticity dynamics depends strongly on the stimulus spectrum via overall system frequency responses. In the case of perturbations that are approximately linear, explicit formulas are found for the dynamics in terms of stimulus spectra via system transfer functions. The resulting theory is applied to a simple model system to reveal how collective effects, especially resonances, can drastically modify system-level plasticity dynamics from that implied by single-neuron analyses. The simplified model illustrates the potential relevance of these effects in applications to brain stimulation, synaptic homeostasis, and epilepsy.

Biophysical mechanisms of multistability in resting-state cortical rhythms

The human alpha (8-12 Hz) rhythm is one of the most prominent, robust, and widely studied attributes of ongoing cortical activity. Contrary to the prevalent notion that it simply "waxes and wanes," spontaneous alpha activity bursts erratically between two distinct modes of activity. We now establish a mechanism for this multistable phenomenon in resting-state cortical recordings by characterizing the complex dynamics of a biophysical model of macroscopic corticothalamic activity. This is achieved by studying the predicted activity of cortical and thalamic neuronal populations in this model as a function of its dynamic stability and the role of nonspecific synaptic noise. We hence find that fluctuating noisy inputs into thalamic neurons elicit spontaneous bursts between low- and high-amplitude alpha oscillations when the system is near a particular type of dynamical instability, namely a subcritical Hopf bifurcation. When the postsynaptic potentials associated with these noisy inputs are modulated by cortical feedback, the SD of power within each of these modes scale in proportion to their mean, showing remarkable concordance with empirical data. Our state-dependent corticothalamic model hence exhibits multistability and scale-invariant fluctuations-key features of resting-state cortical activity and indeed, of human perception, cognition, and behavior-thus providing a unified account of these apparently divergent phenomena.

Synchronization in interdependent networks

We explore the synchronization behavior in interdependent systems, where the one-dimensional (1D) network (the intranetwork coupling strength J(I)) is ferromagnetically intercoupled (the strength J) to the Watts-Strogatz (WS) small-world network (the intranetwork coupling strength J(II)). In the absence of the internetwork coupling (J=0), the former network is well known not to exhibit the synchronized phase at any finite coupling strength, whereas the latter displays the mean-field transition. Through an analytic approach based on the mean-field approximation, it is found that for the weakly coupled 1D network (J(I)≪1) the increase of J suppresses synchrony, because the nonsynchronized 1D network becomes a heavier burden for the synchronization process of the WS network. As the coupling in the 1D network becomes stronger, it is revealed by the renormalization group (RG) argument that the synchronization is enhanced as J(I) is increased, implying that the more enhanced partial synchronization in the 1D network makes the burden lighter. Extensive numerical simulations confirm these expected behaviors, while exhibiting a reentrant behavior in the intermediate range of J(I). The nonmonotonic change of the critical value of J(II) is also compared with the result from the numerical RG calculation.

A spatially extended model for macroscopic spike-wave discharges

Spike-wave discharges are a distinctive feature of epileptic seizures. So far, they have not been reported in spatially extended neural field models. We study a space-independent version of the Amari neural field model with two competing inhibitory populations. We show that this competition leads to robust spike-wave dynamics if the inhibitory populations operate on different time-scales. The spike-wave oscillations present a fold/homoclinic type bursting. From this result we predict parameters of the extended Amari system where spike-wave oscillations produce a spatially homogeneous pattern. We propose this mechanism as a prototype of macroscopic epileptic spike-wave discharges. To our knowledge this is the first example of robust spike-wave patterns in a spatially extended neural field model.

Firing responses of bursting neurons with delayed feedback

Thalamic neurons, which play important roles in the genesis of rhythmic activities of the brain, show various bursting behaviors, particularly modulated by complex thalamocortical feedback via cortical neurons. As a first step to explore this complex neural system and focus on the effects of the feedback on the bursting behavior, a simple loop structure delayed in time and scaled by a coupling strength is added to a recent mean-field model of bursting neurons. Depending on the coupling strength and delay time, the modeled neurons show two distinct response patterns: one entrained to the unperturbed bursting frequency of the neurons and one entrained to the resonant frequency of the loop structure. Transitions between these two patterns are explored in the model's parameter space via extensive numerical simulations. It is found that at a fixed loop delay, there is a critical coupling strength at which the dominant response frequency switches from the unperturbed bursting frequency to the loop-induced one. Furthermore, alternating occurrence of these two response frequencies is observed when the delay varies at fixed coupling strength. The results demonstrate that bursting is coupled with feedback to yield new dynamics, which will provide insights into such effects in more complex neural systems.

Experimental observation of increased fluctuations in an order parameter before epochs of extended brain synchronization

The identification of epileptic seizure precursors has potential clinical relevance. It is conjectured that seizures may be represented by dynamical bifurcations and that an adequate order parameter to characterize brain dynamics is the phase difference in the oscillatory activity of neural systems. In this study, the critical point hypothesis that seizures, or more generally periods of widespread high synchronization, represent bifurcations is empirically tested by monitoring the growth of fluctuations in the putative order parameter of phase differences between magnetoencephalographic and electroencephalographic signals in nearby brain regions in patients with epilepsy and normal subjects during hyperventilation. Implications of the results with regard to epileptic phenomena are discussed.

Dynamics of coupled thalamocortical modules

We develop a model of thalamocortical dynamics using a shared population of thalamic neurons to couple distant cortical regions. Behavior of the model is determined as a function of the connection strengths with shared and unshared populations in the thalamus, either within a relay nucleus or the reticular nucleus. When the coupling is via the reticular nucleus, we locate solutions of the model where distant cortical regions maintain the same activity level, and regions where one region maintains an elevated activity level, suppressing activity in the other. We locate and investigate a region where both types of solutions exist and are stable, yielding a mechanism for spontaneous changes in global activity patterns. Power spectra and coherence are computed, and marked differences in the coherence are found between the two kinds of modes. When, on the other hand, the coupling is via a shared relay nuclei, the features seen with the reticular coupling are absent. These considerations suggest a role for the reticular nucleus in modulating long distance cortical communication.

A method for detecting false bifurcations in dynamical systems: application to neural-field models

In this article, we present a method for tracking changes in curvature of limit cycle solutions that arise due to inflection points. In keeping with previous literature, we term these changes false bifurcations, as they appear to be bifurcations when considering a Poincaré section that is tangent to the solution, but in actual fact the deformation of the solution occurs smoothly as a parameter is varied. These types of solutions arise commonly in electroencephalogram models of absence seizures and correspond to the formation of spikes in these models. Tracking these transitions in parameter space allows regions to be defined corresponding to different types of spike and wave dynamics, that may be of use in clinical neuroscience as a means to classify different subtypes of the more general syndrome.

Axonal velocity distributions in neural field equations

By modelling the average activity of large neuronal populations, continuum mean field models (MFMs) have become an increasingly important theoretical tool for understanding the emergent activity of cortical tissue. In order to be computationally tractable, long-range propagation of activity in MFMs is often approximated with partial differential equations (PDEs). However, PDE approximations in current use correspond to underlying axonal velocity distributions incompatible with experimental measurements. In order to rectify this deficiency, we here introduce novel propagation PDEs that give rise to smooth unimodal distributions of axonal conduction velocities. We also argue that velocities estimated from fibre diameters in slice and from latency measurements, respectively, relate quite differently to such distributions, a significant point for any phenomenological description. Our PDEs are then successfully fit to fibre diameter data from human corpus callosum and rat subcortical white matter. This allows for the first time to simulate long-range conduction in the mammalian brain with realistic, convenient PDEs. Furthermore, the obtained results suggest that the propagation of activity in rat and human differs significantly beyond mere scaling. The dynamical consequences of our new formulation are investigated in the context of a well known neural field model. On the basis of Turing instability analyses, we conclude that pattern formation is more easily initiated using our more realistic propagator. By increasing characteristic conduction velocities, a smooth transition can occur from self-sustaining bulk oscillations to travelling waves of various wavelengths, which may influence axonal growth during development. Our analytic results are also corroborated numerically using simulations on a large spatial grid. Thus we provide here a comprehensive analysis of empirically constrained activity propagation in the context of MFMs, which will allow more realistic studies of mammalian brain activity in the future.

Neurophysiological changes with age probed by inverse modeling of EEG spectra

OBJECTIVE

To investigate age-associated changes in physiologically-based EEG spectral parameters in the healthy population.

METHODS

Eyes-closed EEG spectra of 1498 healthy subjects aged 6-86 years were fitted to a mean-field model of thalamocortical dynamics in a cross-sectional study. Parameters were synaptodendritic rates, cortical wave decay rates, connection strengths (gains), axonal delays for thalamocortical loops, and power normalizations. Age trends were approximated using smooth asymptotically linear functions with a single turning point. We also considered sex differences and relationships between model parameters and traditional quantitative EEG measures.

RESULTS

The cross-sectional data suggest that changes tend to be most rapid in childhood, generally leveling off at age 15-20 years. Most gains decrease in magnitude with age, as does power normalization. Axonal and dendritic delays decrease in childhood and then increase. Axonal delays and gains show small but significant sex differences.

CONCLUSIONS

Mean-field brain modeling allows interpretation of age-associated EEG trends in terms of physiological processes, including the growth and regression of white matter, influencing axonal delays, and the establishment and pruning of synaptic connections, influencing gains.

SIGNIFICANCE

This study demonstrates the feasibility of inverse modeling of EEG spectra as a noninvasive method for investigating large-scale corticothalamic dynamics, and provides a basis for future comparisons.

Comparative study of nonlinear properties of EEG signals of normal persons and epileptic patients

BACKGROUND

Investigation of the functioning of the brain in living systems has been a major effort amongst scientists and medical practitioners. Amongst the various disorder of the brain, epilepsy has drawn the most attention because this disorder can affect the quality of life of a person. In this paper we have reinvestigated the EEGs for normal and epileptic patients using surrogate analysis, probability distribution function and Hurst exponent.

RESULTS

Using random shuffled surrogate analysis, we have obtained some of the nonlinear features that was obtained by Andrzejak et al. [Phys Rev E 2001, 64:061907], for the epileptic patients during seizure. Probability distribution function shows that the activity of an epileptic brain is nongaussian in nature. Hurst exponent has been shown to be useful to characterize a normal and an epileptic brain and it shows that the epileptic brain is long term anticorrelated whereas, the normal brain is more or less stochastic. Among all the techniques, used here, Hurst exponent is found very useful for characterization different cases.

CONCLUSION

In this article, differences in characteristics for normal subjects with eyes open and closed, epileptic subjects during seizure and seizure free intervals have been shown mainly using Hurst exponent. The H shows that the brain activity of a normal man is uncorrelated in nature whereas, epileptic brain activity shows long range anticorrelation.

Bistability and non-Gaussian fluctuations in spontaneous cortical activity

The brain is widely assumed to be a paradigmatic example of a complex, self-organizing system. As such, it should exhibit the classic hallmarks of nonlinearity, multistability, and "nondiffusivity" (large coherent fluctuations). Surprisingly, at least at the very large scale of neocortical dynamics, there is little empirical evidence to support this, and hence most computational and methodological frameworks for healthy brain activity have proceeded very reasonably from a purely linear and diffusive perspective. By studying the temporal fluctuations of power in human resting-state electroencephalograms, we show that, although these simple properties may hold true at some temporal scales, there is strong evidence for bistability and nondiffusivity in key brain rhythms. Bistability is manifest as nonclassic bursting between high- and low-amplitude modes in the alpha rhythm. Nondiffusivity is expressed through the irregular appearance of high amplitude "extremal" events in beta rhythm power fluctuations. The statistical robustness of these observations was confirmed through comparison with Gaussian-rendered phase-randomized surrogate data. Although there is a good conceptual framework for understanding bistability in cortical dynamics, the implications of the extremal events challenge existing frameworks for understanding large-scale brain systems.

Characterisation of cortical activity in response to deep brain stimulation of ventral-lateral nucleus: modelling and experiment

Motivated by its success as a therapeutic treatment in other neurological disorders, most notably Parkinson's disease, Deep Brain Stimulation (DBS) is currently being trialled in a number of patients with drug unresponsive epilepsies. However, the mechanisms by which DBS interferes with neuronal activity linked to the disorder are not well understood. Furthermore, there is a need to identify optimized values of parameters (for example in amplitude/frequency space) of the stimulation protocol with which one aims to achieve the desired outcome. In this paper we characterise the system response to stimulation, to gain an understanding of the role different brain regions play in generating the output observed in EEG. We perform a number of experiments in healthy rats, where the ventral-lateral thalamic nucleus is stimulated using a train of square-waves with different frequency and amplitudes. The response to stimulation in the motor cortex is recorded and the drive-response relationship over frequency/amplitude space is considered. Subsequently, we compare the experimental data with simulations of a mean-field model, finding good agreement between the output of the model and the experimental data--both in the time and frequency domains--when considering a transition to oscillatory activity in the cortex as the frequency of stimulation is increased. Overall, our study suggests that mean-field models can appropriately characterise the stimulus-response relationship of DBS in healthy animals. In this way, it constitutes a first step towards the goal of developing a closed-loop feedback control protocol for suppressing epileptic activity, by adaptively adjusting the stimulation protocol in response to EEG activity.

Transitions to spike-wave oscillations and epileptic dynamics in a human cortico-thalamic mean-field model

In this paper we present a detailed theoretical analysis of the onset of spike-wave activity in a model of human electroencephalogram (EEG) activity, relating this to clinical recordings from patients with absence seizures. We present a complete explanation of the transition from inter-ictal activity to spike and wave using a combination of bifurcation theory, numerical continuation and techniques for detecting the occurrence of inflection points in systems of delay differential equations (DDEs). We demonstrate that the initial transition to oscillatory behaviour occurs as a result of a Hopf bifurcation, whereas the addition of spikes arises as a result of an inflection point of the vector field. Strikingly these findings are consistent with EEG data recorded from patients with absence seizures and we present a discussion of the clinical significance of these results, suggesting potential new techniques for detection and anticipation of seizures.

Two-point heterogeneous connections in a continuum neural field model

We examine a novel heterogeneous connection scheme in a 1D continuum neural field model. Multiple two-point connections are added to a local connection function in order to model the "patchy" connections seen in, for example visual cortex. We use a numerical approach to solve the equations, choosing the locations of the two-point connections stochastically. We observe self-sustained persistent fluctuations of activity which can be classified into two types (one of which is similar to that seen in network models of discrete excitable neurons, the other being particular to this model). We study the effect of parameters such as system size and the range, number and strength of connections, on the probability that a particular realisation of the connections is able to exhibit persistent fluctuations.

Onset of polyspike complexes in a mean-field model of human electroencephalography and its application to absence epilepsy

In this paper, we introduce a modification of a mean-field model used to describe the brain's electrical activity as recorded via electroencephalography (EEG). The focus of the present study is to understand the mechanisms giving rise to the dynamics observed during absence epilepsy, one of the classical generalized syndromes. A systematic study of the data from a number of different subjects with absence epilepsy demonstrates a wide variety of dynamical phenomena in the recorded EEG. In addition to the classical spike and wave activity, there may be polyspike and wave, wave spike or even no discernible spike-wave onset during seizure events. The model we introduce is able to capture all of these different phenomena and we describe the bifurcations giving rise to these different types of seizure activity. We argue that such a model may provide a useful clinical tool for classifying different subclasses of absence epilepsy.

From working memory to epilepsy: dynamics of facilitation and inhibition in a cortical network

Persistent states are believed to be the correlate for short-term or working memory. Using a previously derived model for working memory, we show that disruption of the lateral inhibition can lead to a variety of pathological states. These states are analogs of reflex or pattern-sensitive epilepsy. Simulations, numerical bifurcation analysis, and fast-slow decomposition are used to explore the dynamics of this network.

Derivation and analysis of an ordinary differential equation mean-field model for studying clinically recorded epilepsy dynamics

In this paper we describe how an ordinary differential equation model of corticothalamic interactions may be obtained from a more general system of delay differential equations. We demonstrate that transitions to epileptic dynamics via changes in system parameters are qualitatively the same as in the original model with delay, as well as demonstrating that the onset of epileptic activity may arise due to regions of bistability. Hence, the model presents in one unique framework, two competing theories for the genesis of epileptiform activity. Similarities between model transitions and clinical data are presented and we argue that statistics obtained from, and a parameter estimation of this model may be a potential means of classifying and predicting the onset and offset of seizure activity.

Mean-field modeling of the basal ganglia-thalamocortical system. I Firing rates in healthy and parkinsonian states

Parkinsonism leads to various electrophysiological changes in the basal ganglia-thalamocortical system (BGTCS), often including elevated discharge rates of the subthalamic nucleus (STN) and the output nuclei, and reduced activity of the globus pallidus external (GPe) segment. These rate changes have been explained qualitatively in terms of the direct/indirect pathway model, involving projections of distinct striatal populations to the output nuclei and GPe. Although these populations partly overlap, evidence suggests dopamine depletion differentially affects cortico-striato-pallidal connection strengths to the two pallidal segments. Dopamine loss may also decrease the striatal signal-to-noise ratio, reducing both corticostriatal coupling and striatal firing thresholds. Additionally, nigrostriatal degeneration may cause secondary changes including weakened lateral inhibition in the GPe, and mesocortical dopamine loss may decrease intracortical excitation and especially inhibition. Here a mean-field model of the BGTCS is presented with structure and parameter estimates closely based on physiology and anatomy. Changes in model rates due to the possible effects of dopamine loss listed above are compared with experiment. Our results suggest that a stronger indirect pathway, possibly combined with a weakened direct pathway, is compatible with empirical evidence. However, altered corticostriatal connection strengths are probably not solely responsible for substantially increased STN activity often found. A lower STN firing threshold, weaker intracortical inhibition, and stronger striato-GPe inhibition help explain the relatively large increase in STN rate. Reduced GPe-GPe inhibition and a lower GPe firing threshold can account for the comparatively small decrease in GPe rate frequently observed. Changes in cortex, GPe, and STN help normalize the cortical rate, also in accord with experiments. The model integrates the basal ganglia into a unified framework along with an existing thalamocortical model that already accounts for a wide range of electrophysiological phenomena. A companion paper discusses the dynamics and oscillations of this combined system.

Cortical local and long-range synchronization interplay in human absence seizure initiation

Brain activity relies on transient, fluctuating interactions between segregated neuronal populations. Synchronization within a single and between distributed neuronal clusters reflects the dynamics of these cooperative patterns. Thus absence epilepsy can be used as a model for integrated, large-scale investigation of the emergence of pathological collective dynamics in the brain. Indeed, spike-wave discharges (SWD) of an absence seizure are thought to reflect abnormal cortical hypersynchronization. In this paper, we address two questions: how and where do SWD arise in the human brain? Therefore, we explored the spatio-temporal dynamics of interactions within and between widely distributed cortical sites using magneto-encephalographic recordings of spontaneous absence seizures. We then extracted, from their time-frequency analysis, local synchronization of cortical sources and long-range synchronization linking distant sites. Our analyses revealed a reproducible sequence of 1) long-range desynchronization, 2) increased local synchronization and 3) increased long-range synchronization. Although both local and long-range synchronization displayed different spatio-temporal profiles, their cortical projection within an initiation time window overlap and reveal a multifocal fronto-central network. These observations contradict the classical view of sudden generalized synchronous activities in absence epilepsy. Furthermore, they suggest that brain states transition may rely on multi-scale processes involving both local and distant interactions.

Mean-field modeling of the basal ganglia-thalamocortical system. II Dynamics of parkinsonian oscillations

Neuronal correlates of Parkinson's disease (PD) include a shift to lower frequencies in the electroencephalogram (EEG) and enhanced synchronized oscillations at 3-7 and 7-30 Hz in the basal ganglia, thalamus, and cortex. This study describes the dynamics of a recent physiologically based mean-field model of the basal ganglia-thalamocortical system, and shows how it accounts for many key electrophysiological correlates of PD. Its detailed functional connectivity comprises partially segregated direct and indirect pathways through two populations of striatal neurons, a hyperdirect pathway involving a corticosubthalamic projection, thalamostriatal feedback, and local inhibition in striatum and external pallidum (GPe). In a companion paper, realistic steady-state firing rates were obtained for the healthy state, and after dopamine loss modeled by weaker direct and stronger indirect pathways, reduced intrapallidal inhibition, lower firing thresholds of the GPe and subthalamic nucleus (STN), a stronger projection from striatum to GPe, and weaker cortical interactions. Here it is shown that oscillations around 5 and 20 Hz can arise with a strong indirect pathway, which also causes increased synchronization throughout the basal ganglia. Furthermore, increased theta power with progressive nigrostriatal degeneration is correlated with reduced alpha power and peak frequency, in agreement with empirical results. Unlike the hyperdirect pathway, the indirect pathway sustains oscillations with phase relationships that coincide with those found experimentally. Alterations in the responses of basal ganglia to transient stimuli accord with experimental observations. Reduced cortical gains due to both nigrostriatal and mesocortical dopamine loss lead to slower changes in cortical activity and may be related to bradykinesia. Finally, increased EEG power found in some studies may be partly explained by a lower effective GPe firing threshold, reduced GPe-GPe inhibition, and/or weaker intracortical connections in parkinsonian patients. Strict separation of the direct and indirect pathways is not necessary to obtain these results.

Modeling brain activation patterns for the default and cognitive states

We argue that spatial patterns of cortical activation observed with EEG, MEG and fMRI might arise from spontaneous self-organisation of interacting populations of excitatory and inhibitory neurons. We examine the dynamical behavior of a mean-field cortical model that includes chemical and electrical (gap-junction) synapses, focusing on two limiting cases: the "slow-soma" limit with slow voltage feedback from soma to dendrite, and the "fast-soma" limit in which the feedback action of soma voltage onto dendrite reversal potentials is instantaneous. For slow soma-dendrite feedback, we find a low-frequency (approximately 1 Hz) dynamic Hopf instability, and a stationary Turing instability that catalyzes formation of patterned distributions of cortical firing-rate activity with pattern wavelength approximately 2 cm. Turing instability can only be triggered when gap-junction diffusion between inhibitory neurons is strong, but patterning is destroyed if the tonic level of subcortical excitation is raised sufficiently. Interaction between the Hopf and Turing instabilities may describe the non-cognitive background or "default" state of the brain, as observed by BOLD imaging. In the fast-soma limit, the model predicts a high-frequency Hopf (approximately 35 Hz) instability, and a traveling-wave gamma-band instability that manifests as a 2-D standing-wave pattern oscillating in place at approximately 30 Hz. Small levels of inhibitory diffusion enhance and broaden the definition of the gamma antinodal regions by suppressing higher-frequency spatial modes, but gamma emergence is not contingent on the presence of inhibitory gap junctions; higher levels of diffusion suppress gamma activity. Fast-soma instabilities are enhanced by increased subcortical stimulation. Prompt soma-dendrite feedback may be an essential component of the genesis and large-scale cortical synchrony of gamma activity observed at the point of cognition.

Modeling brain dynamics using computational neurogenetic approach

The paper introduces a novel computational approach to brain dynamics modeling that integrates dynamic gene-protein regulatory networks with a neural network model. Interaction of genes and proteins in neurons affects the dynamics of the whole neural network. Through tuning the gene-protein interaction network and the initial gene/protein expression values, different states of the neural network dynamics can be achieved. A generic computational neurogenetic model is introduced that implements this approach. It is illustrated by means of a simple neurogenetic model of a spiking neural network of the generation of local field potential. Our approach allows for investigation of how deleted or mutated genes can alter the dynamics of a model neural network. We conclude with the proposal how to extend this approach to model cognitive neurodynamics.

Modeling distributed axonal delays in mean-field brain dynamics

The range of conduction delays between connected neuronal populations is often modeled as a single discrete delay, assumed to be an effective value averaging over all fiber velocities. This paper shows the effects of distributed delays on signal propagation. A distribution acts as a linear filter, imposing an upper frequency cutoff that is inversely proportional to the delay width. Distributed thalamocortical and corticothalamic delays are incorporated into a physiologically based mean-field model of the cortex and thalamus to illustrate their effects on the electroencephalogram (EEG). The power spectrum is acutely sensitive to the width of the thalamocortical delay distribution, and more so than the corticothalamic distribution, because all input signals must travel along the thalamocortical pathway. This imposes a cutoff frequency above which the spectrum is overly damped. The positions of spectral peaks in the resting EEG depend primarily on the distribution mean, with only weak dependences on distribution width. Increasing distribution width increases the stability of fixed point solutions. A single discrete delay successfully approximates a distribution for frequencies below a cutoff that is inversely proportional to the delay width, provided that other model parameters are moderately adjusted. A pair of discrete delays together having the same mean, variance, and skewness as the distribution approximates the distribution over the same frequency range without needing parameter adjustment. Delay distributions with large fractional widths are well approximated by low-order differential equations.

A study on the possible usefulness of detrended fluctuation analysis of the electroencephalogram background activity in Alzheimer's disease

We studied the EEG background activity of Alzheimer's disease (AD) patients with detrended fluctuation analysis (DFA). DFA provides an estimation of the scaling information and long-range correlations in time series. We recorded the EEG in 11 AD patients and 11 age-matched controls. Our results showed two scaling regions in all subjects' channels (for limited time scales from 0.01 to 0.04 s and from 0.08 to 0.43 s, respectively), with a clear bend when their corresponding slopes (alpha(1) and alpha(2)) were different. No significant differences between groups were found with alpha(1). However, alpha(2) values were significantly lower in control subjects at electrodes T5, T6, and O1 (p < 0.01, Student's t-test). These findings suggest that the scaling behavior of the EEG is sensitive to AD. Although alpha(2) values allowed us to separate AD patients and controls, accuracies were lower than with spectral analysis. However, a forward stepwise linear discriminant analysis with a leave-one-out cross-validation procedure showed that the combined use of DFA and spectral analysis could improve the diagnostic accuracy of each individual technique. Thus, although spectral analysis outperforms DFA, the combined use of both techniques may increase the insight into brain dysfunction in AD.