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

The dynamic brain: from spiking neurons to neural masses and cortical fields

The cortex is a complex system, characterized by its dynamics and architecture, which underlie many functions such as action, perception, learning, language, and cognition. Its structural architecture has been studied for more than a hundred years; however, its dynamics have been addressed much less thoroughly. In this paper, we review and integrate, in a unifying framework, a variety of computational approaches that have been used to characterize the dynamics of the cortex, as evidenced at different levels of measurement. Computational models at different space-time scales help us understand the fundamental mechanisms that underpin neural processes and relate these processes to neuroscience data. Modeling at the single neuron level is necessary because this is the level at which information is exchanged between the computing elements of the brain; the neurons. Mesoscopic models tell us how neural elements interact to yield emergent behavior at the level of microcolumns and cortical columns. Macroscopic models can inform us about whole brain dynamics and interactions between large-scale neural systems such as cortical regions, the thalamus, and brain stem. Each level of description relates uniquely to neuroscience data, from single-unit recordings, through local field potentials to functional magnetic resonance imaging (fMRI), electroencephalogram (EEG), and magnetoencephalogram (MEG). Models of the cortex can establish which types of large-scale neuronal networks can perform computations and characterize their emergent properties. Mean-field and related formulations of dynamics also play an essential and complementary role as forward models that can be inverted given empirical data. This makes dynamic models critical in integrating theory and experiments. We argue that elaborating principled and informed models is a prerequisite for grounding empirical neuroscience in a cogent theoretical framework, commensurate with the achievements in the physical sciences.

Mean field model of acetylcholine mediated dynamics in the thalamocortical system

A recent continuum model of the large scale electrical activity of the thalamocortical system is generalized to include cholinergic modulation. The model is examined analytically and numerically to determine the effect of acetylcholine (ACh) on its steady states, linear stability, spectrum, and temporal responses. Changing the ACh concentration moves the system between zones of one, three, and five steady states, showing that neuromodulation of synaptic strength is a possible mechanism by which multiple steady states emerge in the brain. The lowest firing rate steady state is always stable, and subsequent fixed points alternate between stable and unstable. Increasing ACh concentration changes the form of the spectrum. Increasing the tonic level of ACh concentration increases the magnitudes of the N100 and P200 in the evoked response potential (ERP), without changing the timing of these peaks. Driving the system with a pulse of cholinergic activity results in a transient increase in the firing rate of cortical neurons that lasts over 10s. Step-like increases in cortical ACh concentration cause increases in the firing rate of cortical neurons, with rapid responses due to fast acting nicotinic receptors and slower responses due to muscarinic receptor suppression of intracortical connections.

Narrow-band oscillations in probabilistic cellular automata

Dynamical properties of neural populations are studied using probabilistic cellular automata. Previous work demonstrated the emergence of critical behavior as the function of system noise and density of long-range axonal connections. Finite-size scaling theory identified critical properties, which were consistent with properties of a weak Ising universality class. The present work extends the studies to neural populations with excitatory and inhibitory interactions. It is shown that the populations can exhibit narrow-band oscillations when confined to a range of inhibition levels, with clear boundaries marking the parameter region of prominent oscillations. Phase diagrams have been constructed to characterize unimodal, bimodal, and quadromodal oscillatory states. The significance of these findings is discussed in the context of large-scale narrow-band oscillations in neural tissues, as observed in electroencephalographic and magnetoencephalographic measurements.

Computer modelling of epilepsy

Epilepsy is a complex set of disorders that can involve many areas of the cortex, as well as underlying deep-brain systems. The myriad manifestations of seizures, which can be as varied as déjà vu and olfactory hallucination, can therefore give researchers insights into regional functions and relations. Epilepsy is also complex genetically and pathophysiologically: it involves microscopic (on the scale of ion channels and synaptic proteins), macroscopic (on the scale of brain trauma and rewiring) and intermediate changes in a complex interplay of causality. It has long been recognized that computer modelling will be required to disentangle causality, to better understand seizure spread and to understand and eventually predict treatment efficacy. Over the past few years, substantial progress has been made in modelling epilepsy at levels ranging from the molecular to the socioeconomic. We review these efforts and connect them to the medical goals of understanding and treating the disorder.

Computational models of epileptic activity: a bridge between observation and pathophysiological interpretation

Epilepsy is a neurological disorder characterized by the recurrence of seizures. It affects 50 million people worldwide. Although a considerable number of new antiepileptic drugs with reduced side effects and toxicity have been introduced since the 1950s, 30% of patients remain pharmacoresistant. Although epilepsy research is making progress, advances in understanding drug resistance have been hampered by the complexity of the underlying neuronal systems responsible for epileptic activity. In such systems where short- or long-term plasticity plays a role, pathophysiological alterations may take place at subcellular (i.e., membrane ion channels and neurotransmitter receptors), cellular (neurons), tissular (networks of neurons) and regional (networks of networks of neurons) scales. In such a context, the demand for integrative approaches is high and neurocomputational models become recognized tools for tackling the complexity of epileptic phenomena. The purpose of this report is to provide an overview on computational modeling as a way of structuring and interpreting multimodal data recorded from the epileptic brain. Some examples are briefly described, which illustrate how computational models closely related with either experimental or clinical data can markedly advance our understanding of essential issues in epilepsy such as the transition from background to seizure activity. A commentary is also made on the potential use of such models in the study of therapeutic strategies such as rational drug design or electrical stimulations.

Characteristics of temporal fluctuations in the hyperpolarized state of the cortical slow oscillation

We present evidence for the hypothesis that transitions between the low- and high-firing states of the cortical slow oscillation correspond to neuronal phase transitions. By analyzing intracellular recordings of the membrane potential during the cortical slow oscillation in rats, we quantify the temporal fluctuations in power and the frequency centroid of the power spectrum in the period of time before "down" to "up" transitions. By taking appropriate averages over such events, we present these statistics as a function of time before transition. The results demonstrate an increase in fluctuation power and time scale broadly consistent with the slowing of systems close to phase transitions. The analysis is complicated and limited by the difficulty in identifying when transitions begin, and removing dc trends in membrane potential.

A dynamic extraversion model. The brain's response to a single dose of a stimulant drug

The aim of this paper is to present a mathematical dynamic modelling of the effect a stimulant drug has on different people which, at the same time, can be a useful tool for future brain studies. To this end, a dynamic model of the evolution of extraversion (considering its tonic and phasic aspects) has been constructed taking into account the unique personality trait theory and the general modelling methodology. This model consists of a delayed differential equation which, on one hand, considers that the active stimulus, a consequence of a single intake, is not constant; on the other hand, it contemplates that the state variable representing the phasic extraversion also represents the brain activation. The derivative of this state variable is calculated as the sum of the homeostatic control flow, the excitatory effect flow and the inhibitor effect flow. The solutions of this equation relate the tonic activation of an individual (that characterizes his or her personality) with his or her phasic activation level, whose evolution over time describes the organism's response to a single drug intake. These solutions quantitatively reproduce the predictions of current personality theories and anticipate vulnerability to drug misuse and addiction development.

Controlling limit-cycle behaviors of brain activity

The limit cycles of brain activity are studied using a compact continuum model that reproduces the main features of electroencephalographic signals, including bifurcations of fixed points and limit cycles in seizures. Frequencies and amplitudes are predicted analytically and related to physiology. Gaussian stimuli yield two distinct evoked responses in the linearly stable zone, consistent with experiment. Limit cycles can be initiated or suppressed by control signals or stimuli.

Physiology-based modeling of cortical auditory evoked potentials

Evoked potentials are the transient electrical responses caused by changes in the brain following stimuli. This work uses a physiology-based continuum model of neuronal activity in the human brain to calculate theoretical cortical auditory evoked potentials (CAEPs) from the model's linearized response. These are fitted to experimental data, allowing the fitted parameters to be related to brain physiology. This approach yields excellent fits to CAEP data, which can then be compared to fits of EEG spectra. It is shown that the differences between resting eyes-open EEG and standard CAEPs can be explained by changes in the physiology of populations of neurons in corticothalamic pathways, with notable similarities to certain aspects of slow-wave sleep. This pilot study demonstrates the ability of our model-based fitting method to provide information on the underlying physiology of the brain that is not available using standard methods.

Physiologically based calculation of steady-state evoked potentials and cortical wave velocities

Steady-state evoked potentials (SSEPs) elicited by sinusoidal stimuli are predicted from a physiologically-based model, including bielectrode and volume conduction effects. Comparison with visual SSEPs yields constraints on phase and latency of the retinothalamic transfer function that are consistent with experiment. Predictions of phase velocities measured as SSEPs cross the cortex are consistent with low values measured for slow waves in sleep, while resonant behavior induced by corticothalamic loops, especially near the alpha peak, contributes to wide scatter in waking-state phase velocity measurements comparable to effects from volume conduction. The common use of bielectrode derivations to compensate for volume conduction effects is examined and shown to be incomplete, tending to lead to underestimates of phase velocity, especially at low frequencies and near the alpha peak, due to incorrect elimination of true long-wavelength contributions to the SSEP.

Mean field model of acetylcholine mediated dynamics in the cerebral cortex

A recent continuum model of the large scale electrical activity of the cerebral cortex is generalized to include cholinergic modulation. In this model, dynamic modulation of synaptic strength acts over the time scales of nicotinic and muscarinic receptor action. The cortical model is analyzed to determine the effect of acetylcholine (ACh) on its steady states, linear stability, spectrum, and temporal responses to changes in subcortical input. ACh increases the firing rate in steady states of the system. Changing ACh concentration does not introduce oscillatory behavior into the system, but increases the overall spectral power. Model responses to pulses in subcortical input are affected by the tonic level of ACh concentration, with higher levels of ACh increasing the magnitude firing rate response of excitatory cortical neurons to pulses of subcortical input. Numerical simulations are used to explore the temporal dynamics of the model in response to changes in ACh concentration. Evidence is seen of a transition from a state in which intracortical inputs are emphasized to a state where thalamic afferents have enhanced influence. Perturbations in ACh concentration cause changes in the firing rate of cortical neurons, with rapid responses due to fast acting facilitatory effects of nicotinic receptors on subcortical afferents, and slower responses due to muscarinic suppression of intracortical connections. Together, these numerical simulations demonstrate that the actions of ACh could be a significant factor modulating early components of evoked response potentials.

Neural rate equations for bursting dynamics derived from conductance-based equations

A method of obtaining rate equations from conductance-based equations is developed and applied to fast-spiking and bursting neocortical neurons. It involves splitting systems of conductance-based equations into fast and slow subsystems, and averaging the effects of fast terms that drive the slowly varying quantities by showing that their average is closely proportional to the firing rate. The dependence of the firing rate on the injected current is then approximated in the analysis. The resulting behavior of the slow variables is then substituted back into the fast equations, with the further approximation of replacing the fast voltages in these terms by effective values. For bursting neurons the method yields two coupled limit-cycle oscillators: a self-exciting oscillator for the slow variables that commences limit-cycle oscillations at a critical current and modulates a fast spike-generating oscillator, thereby leading to slowly modulated bursts with a group of spikes in each burst. The dynamics of these coupled oscillators are then verified against those of the conductance-based equations. Finally, it is shown how to place the results in a form suitable for use in mean-field equations for neural population dynamics.

Visual gamma oscillations: waves, correlations, and other phenomena, including comparison with experimental data

Mean-field theory of brain dynamics is applied to explain the properties of gamma (> or approximately 30 Hz) oscillations of cortical activity often seen during vision experiments. It is shown that mm-scale patchy connections in the primary visual cortex can support collective gamma oscillations with the correct frequencies and spatial structure, even when driven by uncorrelated inputs. This occurs via resonances associated with the the periodic modulation of the network connections, rather than being due to single-cell properties alone. Near-resonant gamma waves are shown to obey the Schrödinger equation, which enables techniques and insights from quantum theory to be used in exploring these classical oscillations. Resulting predictions for gamma responses to stimuli account in a unified way for a wide range of experimental results, including why oscillations and zero-lag synchrony are associated, and variations in correlation functions with time delay, intercellular distance, and stimulus features. They also imply that gamma oscillations may enable a form of frequency multiplexing of neural signals. Most importantly, it is shown that correlations reproduce experimental results that show maximal correlations between cells that respond to related features, but little correlation with other cells, an effect that has been argued to be associated with segmentation of a scene into separate objects. Consistency with infill of missing contours and increase in response with length of bar-shaped stimuli are discussed. Background correlations expected in the absence of stimulation are also calculated and shown to be consistent in form with experimental measurements and similar to stimulus-induced correlations in structure. Finally, possible links of gamma instabilities to certain classes of photically induced seizures and visual hallucinations are discussed.

Spike-rate adaptation and neuronal bursting in a mean-field model of brain activity

Spike-rate adaptation is investigated within a mean-field model of brain activity. Two different mechanisms of negative feedback are considered; one involving modulation of the mean firing threshold, and the other, modulation of the mean synaptic strength. Adaptation to a constant stimulus is shown to take place for both mechanisms, and limit-cycle oscillations in the firing rate corresponding to bursts of neuronal activity are investigated. These oscillations are found to result from a Hopf bifurcation when the equilibrium lies between the local maximum and local minimum of a given nullcline. Oscillations with amplitudes significantly below the maximum firing rate are found over a narrow range of possible equilibriums.

Gap junctions mediate large-scale Turing structures in a mean-field cortex driven by subcortical noise

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.

Compact dynamical model of brain activity

A compact physiologically based mean-field formulation of brain dynamics is proposed to model observed brain activity and electroencephalographic (EEG) signals. In contrast to existing formulations, which are more detailed and complicated, our model is described by a single second-order delay differential equation that encapsulates salient aspects of the physiology. The model captures essential features of activity mediated by fast corticocortical connections and delayed feedbacks via extracortical pathways and external stimuli. In the linear regime, these features can be simply expressed by three coefficients derived from the properties of these physiological pathways and explicit nonlinear approximations are also derived. This compact model successfully reproduces the main features of experimental EEG's and the predictions of previous models, including resonance peaks in EEG spectra and nonlinear dynamics. As an illustration, key features of the dynamics of epileptic seizures are shown to be reproduced by the model. Due to its compact form, the model will facilitate insight into nonlinear brain dynamics via standard nonlinear techniques and will guide analysis and investigation of more complex models. It is thus a useful tool for analyzing complex brain activity, especially when it exhibits low-dimensional dynamics.

Pathological pattern formation and cortical propagation of epileptic seizures

The stochastic partial differential equations (SPDEs) stated by Steyn-Ross and co-workers constitute a model of mesoscopic electrical activity of the human cortex. A simplification in which spatial variation and stochastic input are neglected yields ordinary differential equations (ODEs), which are amenable to analysis by techniques of dynamical systems theory. Bifurcation diagrams are developed for the ODEs with increased subcortical excitation, showing that the model predicts oscillatory electrical activity in a large range of parameters. The full SPDEs with increased subcortical excitation produce travelling waves of electrical activity. These model results are compared with electrocortical data recorded at two subdural electrodes from a human subject undergoing a seizure. The model and observational results agree in two important respects during seizure: (i) the average frequency of maximum power, and (ii) the speed of spatial propagation of voltage peaks. This suggests that seizing activity on the human cortex may be understood as an example of pathological pattern formation. Included is a discussion of the applications and limitations of these results.

Patchy propagators, brain dynamics, and the generation of spatially structured gamma oscillations

Propagator theory of brain dynamics is generalized to incorporate a new class of patchy propagators that enable treatment of approximately periodic structures such as are seen in the visual cortex. Complex response fields are also incorporated to allow for features such as orientation preference and wave-number selectivity. The results are applied to the corticothalamic system associated with the primary visual cortex. It is found that this system can generate gamma ( > or = 30 Hz) oscillations during stimulation, whose properties are consistent with experimental findings on gamma frequency and bandwidth, and existence of fine-scale spatial structure. It is found that a potential resonance is associated with each reciprocal lattice vector corresponding to periodic modulations of the propagators. It is found that the lowest resonances are the most likely to give rise to noticeable spectral peaks and increases of correlation amplitude, length, and time, and that these aspects are prominent only if the system is close to marginal stability, in accord with previous measurements and discussions of cortical stability. These features also enable gamma resonances to be stimulus-evoked, with substantial resonance sharpening for relatively small changes in mean neural firing rate. The results also imply dependence of gamma frequency on stimulus features.

EEG MARKERS FOR COGNITIVE DECLINE IN ELDERLY SUBJECTS WITH SUBJECTIVE MEMORY COMPLAINTS

New treatments for Alzheimer's disease require early detection of cognitive decline. Most studies seeking to identify markers of early cognitive decline have focused on a limited number of measures. We sought to establish the profile of brain function measures which best define early neuropsychological decline. We compared subjects with subjective memory complaints to normative controls on a wide range of EEG derived measures, including a new measure of event-related spatio-temporal waves and biophysical modeling, which derives anatomical and physiological parameters based on subject's EEG measurements. Measures that distinguished the groups were then related to cognitive performance on a variety of learning and executive function tasks. The EEG measures include standard power measures, peak alpha frequency, EEG desynchronization to eyes-opening, and global phase synchrony. The most prominent differences in subjective memory complaint subjects were elevated alpha power and an increased number of spatio-temporal wave events. Higher alpha power and changes in wave activity related most strongly to a decline in verbal memory performance in subjects with subjective memory complaints, and also declines in maze performance and working memory reaction time. Interestingly, higher alpha power and wave activity were correlated with improved performance in reverse digit span in the subjective memory complaint group. The modeling results suggest that differences in the subjective memory complaint subjects were due to a decrease in cortical and thalamic inhibitory gains and slowed dendritic time-constants. The complementary profile that emerges from the variety of measures and analyses points to a nonlinear progression in electrophysiological changes from early neuropsychological decline to late-stage dementia, and electrophysiological changes in subjective memory complaint that vary in their relationships to a range of memory-related tasks.

Analysis of oscillatory patterns in the human sleep EEG using a novel detection algorithm

The different brain states during sleep are characterized by the occurrence of distinct oscillatory patterns such as spindles or delta waves. Using a new algorithm to detect oscillatory events in the electroencephalogram (EEG), we studied their properties and changes throughout the night. The present approach was based on the idea that the EEG may be described as a superposition of stochastically driven harmonic oscillators with damping and frequency varying in time. This idea was implemented by fitting autoregressive models to the EEG data. Oscillatory events were detected, whenever the damping of one or more frequencies was below a predefined threshold. Sleep EEG data of eight healthy young males were analyzed (four nights per subject). Oscillatory events occurred mainly in three frequency ranges, which correspond roughly to the classically defined delta (0-4.5 Hz), alpha (8-11.5 Hz) and sigma (11.5-16 Hz) bands. Their incidence showed small intra- but large inter-individual differences, in particular with respect to alpha events. The incidence and frequency of the events was characteristic for sleep stages and non-rapid eye movement (REM)-REM sleep cycles. The mean event frequency of delta and sigma (spindle) events decreased with the deepening of sleep. It was higher in the second half of the night compared with the first one for delta, alpha and sigma oscillations. The algorithm provides a general framework to detect and characterize oscillatory patterns in the EEG and similar signals.

A unifying explanation of primary generalized seizures through nonlinear brain modeling and bifurcation analysis

The aim of this paper is to explain critical features of the human primary generalized epilepsies by investigating the dynamical bifurcations of a nonlinear model of the brain's mean field dynamics. The model treats the cortex as a medium for the propagation of waves of electrical activity, incorporating key physiological processes such as propagation delays, membrane physiology, and corticothalamic feedback. Previous analyses have demonstrated its descriptive validity in a wide range of healthy states and yielded specific predictions with regards to seizure phenomena. We show that mapping the structure of the nonlinear bifurcation set predicts a number of crucial dynamic processes, including the onset of periodic and chaotic dynamics as well as multistability. Quantitative study of electrophysiological data supports the validity of these predictions. Hence, we argue that the core electrophysiological and cognitive differences between tonic-clonic and absence seizures are predicted and interrelated by the global bifurcation diagram of the model's dynamics. The present study is the first to present a unifying explanation of these generalized seizures using the bifurcation analysis of a dynamical model of the brain.

Multiscale brain modelling

A central difficulty of brain modelling is to span the range of spatio-temporal scales from synapses to the whole brain. This paper overviews results from a recent model of the generation of brain electrical activity that incorporates both basic microscopic neurophysiology and large-scale brain anatomy to predict brain electrical activity at scales from a few tenths of a millimetre to the whole brain. This model incorporates synaptic and dendritic dynamics, nonlinearity of the firing response, axonal propagation and corticocortical and corticothalamic pathways. Its relatively few parameters measure quantities such as synaptic strengths, corticothalamic delays, synaptic and dendritic time constants, and axonal ranges, and are all constrained by independent physiological measurements. It reproduces quantitative forms of electroencephalograms seen in various states of arousal, evoked response potentials, coherence functions, seizure dynamics and other phenomena. Fitting model predictions to experimental data enables underlying physiological parameters to be inferred, giving a new non-invasive window into brain function that complements slower, but finer-resolution, techniques such as fMRI. Because the parameters measure physiological quantities relating to multiple scales, and probe deep structures such as the thalamus, this will permit the testing of a range of hypotheses about vigilance, cognition, drug action and brain function. In addition, referencing to a standardized database of subjects adds strength and specificity to characterizations obtained.

Propagator theory of brain dynamics

A physiologically based continuum model of brain dynamics is extended to incorporate arbitrary numbers of structures and neural populations, multiple outgoing fields of activity from a single population of neurons to various targets, improved treatment of converging or diverging projections and mesoscopic structure, and generalized connections to quantities observable via electroencephalography and other methods. The results are applied to study the corticothalamic system, predicting an intracortical resonance that leads to enhancements of electroencephalographic activity in the gamma (>30 Hz) range. This resonance involves feedback loops incorporating slow, short-range inhibitory fibers.

Modelling event-related responses in the brain

The aim of this work was to investigate the mechanisms that shape evoked electroencephalographic (EEG) and magneto-encephalographic (MEG) responses. We used a neuronally plausible model to characterise the dependency of response components on the models parameters. This generative model was a neural mass model of hierarchically arranged areas using three kinds of inter-area connections (forward, backward and lateral). We investigated how responses, at each level of a cortical hierarchy, depended on the strength of connections or coupling. Our strategy was to systematically add connections and examine the responses of each successive architecture. We did this in the context of deterministic responses and then with stochastic spontaneous activity. Our aim was to show, in a simple way, how event-related dynamics depend on extrinsic connectivity. To emphasise the importance of nonlinear interactions, we tried to disambiguate the components of event-related potentials (ERPs) or event-related fields (ERFs) that can be explained by a linear superposition of trial-specific responses and those engendered nonlinearly (e.g., by phase-resetting). Our key conclusions were; (i) when forward connections, mediating bottom-up or extrinsic inputs, are sufficiently strong, nonlinear mechanisms cause a saturation of excitatory interneuron responses. This endows the system with an inherent stability that precludes nondissipative population dynamics. (ii) The duration of evoked transients increases with the hierarchical depth or level of processing. (iii) When backward connections are added, evoked transients become more protracted, exhibiting damped oscillations. These are formally identical to late or endogenous components seen empirically. This suggests that late components are mediated by reentrant dynamics within cortical hierarchies. (iv) Bilateral connections produce similar effects to backward connections but can also mediate zero-lag phase-locking among areas. (v) Finally, with spontaneous activity, ERPs/ERFs can arise from two distinct mechanisms: For low levels of (stimulus related and ongoing) activity, the systems response conforms to a quasi-linear superposition of separable responses to the fixed and stochastic inputs. This is consistent with classical assumptions that motivate trial averaging to suppress spontaneous activity and disclose the ERP/ERF. However, when activity is sufficiently high, there are nonlinear interactions between the fixed and stochastic inputs. This interaction is expressed as a phase-resetting and represents a qualitatively different explanation for the ERP/ERF.

Stimulant drug action in attention deficit hyperactivity disorder (ADHD): inference of neurophysiological mechanisms via quantitative modelling

OBJECTIVE

To infer the neural mechanisms underlying tonic transitions in the electroencephalogram (EEG) in 11 adolescents diagnosed with attention deficit hyperactivity disorder (ADHD) before and after treatment with stimulant medication.

METHODS

A biophysical model was used to analyse electroencephalographic (EEG) measures of tonic brain activity at multiple scalp sites before and after treatment with medication.

RESULTS

It was observed that stimulants had the affect of significantly reducing the parameter controlling activation in the intrathalamic pathway involving the thalamic reticular nucleus (TRN) and the parameter controlling excitatory cortical activity. The effect of stimulant medication was also found to be preferentially localized within subcortical nuclei projecting towards frontal and central scalp sites.

CONCLUSIONS

It is suggested that the action of stimulant medication occurs via suppression of the locus coeruleus, which in turn reduces stimulation of the TRN, and improves cortical arousal. The effects localized to frontal and central sites are consistent with the occurrence of frontal delta-theta EEG abnormalities in ADHD, and existing theories of hypoarousal.

SIGNIFICANCE

To our knowledge, this is the first study where a detailed biophysical model of the brain has been used to estimate changes in neurophysiological parameters underlying the effects of stimulant medication in ADHD.

Modeling the minimally conscious state: measurements of brain function and therapeutic possibilities

The minimally conscious state (MCS) defines a functional level of recovery following severe brain injuries. Patients in MCS demonstrate unequivocal evidence of response to their environment yet fail to recover the ability to communicate. Drawing on recent functional brain-imaging studies, pathological data, and neurophysiological investigations, models of brain function in MCS are proposed. MCS models are compared and contrasted with models of the vegetative state (VS), a condition characterized by wakeful appearance and unconsciousness. VS reflects a total loss of cognitive function and failure to recover basic aspects of the normal physiologic brain state associated with wakefulness. MCS may represent a recovery of the minimal dynamic architecture required to organize behavioral sets and respond to sensory stimuli. Several pathophysiological mechanisms that might limit further recovery in MCS patients are considered. Implications for future research directions and possible therapeutic strategies are reviewed.

Analysis of the electroencephalographic activity associated with thalamic tumors

A physiologically based model of corticothalamic dynamics is used to investigate the electroencephalographic (EEG) activity associated with tumors of the thalamus. Tumor activity is modeled by introducing localized two-dimensional spatial non-uniformities into the model parameters, and calculating the resulting activity via the coupling of spatial eigenmodes. The model is able to reproduce various qualitative features typical of waking eyes-closed EEGs in the presence of a thalamic tumor, such as the appearance of abnormal peaks at theta ( approximately 3Hz) and spindle ( approximately 12Hz) frequencies, the attenuation of normal eyes-closed background rhythms, and the onset of epileptic activity, as well as the relatively normal EEGs often observed. The results indicate that the abnormal activity at theta and spindle frequencies arises when a small portion of the brain is forced into an over-inhibited state due to the tumor, in which there is an increase in the firing of (inhibitory) thalamic reticular neurons. The effect is heightened when there is a concurrent decrease in the firing of (excitatory) thalamic relay neurons, which are in any case inhibited by the reticular ones. This is likely due to a decrease in the responsiveness of the peritumoral region to cholinergic inputs from the brainstem, and a corresponding depolarization of thalamic reticular neurons, and hyperpolarization of thalamic relay neurons, similar to the mechanism active during slow-wave sleep. The results indicate that disruption of normal thalamic activity is essential to generate these spectral peaks. Furthermore, the present work indicates that high-voltage and epileptiform EEGs are caused by a tumor-induced local over-excitation of the thalamus, which propagates to the cortex. Experimental findings relating to local over-inhibition and over-excitation are discussed. It is also confirmed that increasing the size of the tumor leads to greater abnormalities in the observable EEG. The usefulness of EEG for localizing the tumor is investigated.

Unifying and interpreting the spectral wavenumber content of EEGs, ECoGs, and ERPs

A biological model of corticothalamic dynamics is used to investigate the spatial power spectrum (wavenumber spectrum) of electrical activity in the brain. The model provides a single framework for unifying different aspects of activity. Comparisons of the predicted spectra with published electrocorticographic, electroencephalographic, and evoked response potential data enable physiology and anatomy to be inferred, producing results which are complementary to those obtained from comparisons in the frequency domain; the inferred quantities are consistent with, and complementary to, direct physiological and anatomical measurements. We also use the model to quantify the interdependence of the wavenumber and frequency domains, and deduce that further experiments that cover large wavenumber and frequency ranges simultaneously would greatly increase our knowledge of brain function. We conclude that both the frequency and wavenumber domains should be studied in order to build the fullest picture of brain dynamics: the two domains are both complementary and interdependent.

Estimation of multiscale neurophysiologic parameters by electroencephalographic means

It is shown that new model-based electroencephalographic (EEG) methods can quantify neurophysiologic parameters that underlie EEG generation in ways that are complementary to and consistent with standard physiologic techniques. This is done by isolating parameter ranges that give good matches between model predictions and a variety of experimental EEG-related phenomena simultaneously. Resulting constraints range from the submicrometer synaptic level to length scales of tens of centimeters, and from timescales of around 1 ms to 1 s or more, and are found to be consistent with independent physiologic and anatomic measures. In the process, a new method of obtaining model parameters from the data is developed, including a Monte Carlo implementation for use when not all input data are available. Overall, the approaches used are complementary to other methods, constraining allowable parameter ranges in different ways and leading to much tighter constraints overall. EEG methods often provide the most restrictive individual constraints. This approach opens a new, noninvasive window on quantitative brain analysis, with the ability to monitor temporal changes, and the potential to map spatial variations. Unlike traditional phenomenologic quantitative EEG measures, the methods proposed here are based explicitly on physiology and anatomy.

Controlled test for predictive power of Lyapunov exponents: their inability to predict epileptic seizures

Lyapunov exponents are a set of fundamental dynamical invariants characterizing a system's sensitive dependence on initial conditions. For more than a decade, it has been claimed that the exponents computed from electroencephalogram (EEG) or electrocorticogram (ECoG) signals can be used for prediction of epileptic seizures minutes or even tens of minutes in advance. The purpose of this paper is to examine the predictive power of Lyapunov exponents. Three approaches are employed. (1) We present qualitative arguments suggesting that the Lyapunov exponents generally are not useful for seizure prediction. (2) We construct a two-dimensional, nonstationary chaotic map with a parameter slowly varying in a range containing a crisis, and test whether this critical event can be predicted by monitoring the evolution of finite-time Lyapunov exponents. This can thus be regarded as a "control test" for the claimed predictive power of the exponents for seizure. We find that two major obstacles arise in this application: statistical fluctuations of the Lyapunov exponents due to finite time computation and noise from the time series. We show that increasing the amount of data in a moving window will not improve the exponents' detective power for characteristic system changes, and that the presence of small noise can ruin completely the predictive power of the exponents. (3) We report negative results obtained from ECoG signals recorded from patients with epilepsy. All these indicate firmly that, the use of Lyapunov exponents for seizure prediction is practically impossible as the brain dynamical system generating the ECoG signals is more complicated than low-dimensional chaotic systems, and is noisy.

Oscillatory phase transition and pulse propagation in noisy integrate-and-fire neurons

We study nonlocally coupled noisy integrate-and-fire neurons with the Fokker-Planck equation. A propagating pulse state and a wavy state appear as a phase transition from an asynchronous state. We also find a solution in which traveling pulses are emitted periodically from a pacemaker region.

Spatially uniform and nonuniform analyses of electroencephalographic dynamics,with application to the topography of the alpha rhythm

Corticothalamic dynamics are investigated using a model in which spatial nonuniformities are incorporated via the coupling of spatial eigenmodes. Comparison of spectra generated using the nonuniform analysis with those generated using a uniform one demonstrates that, for most frequencies, local activity is only weakly dependent on activity elsewhere in the cortex; however, dispersion of low-wave-number activity ensures that distant dynamics influence local dynamics at low frequencies (below approximately 2 Hz ), and at the alpha frequency (approximately 10 Hz ), where propagating signals are inherently weakly damped, and wavelengths are large. When certain model parameters have similar spatial profiles, as is expected from physiology, the low-frequency discrepancies tend to cancel, and the uniform analysis with local parameter values is an adequate approximation to the full nonuniform one across the whole spectrum, at least for large-scale nonuniformities. After comparing the uniform and nonuniform analyses, we consider one possible application of the nonuniform analysis: studying the phenomenon of occipital alpha dominance, whereby the alpha frequency and power are greater at the back of the head (occipitally) than at the front. In order to infer realistic nonuniformities in the model parameters, the uniform version of the model is first fitted to data recorded from 98 normal subjects in a waking, eyes-closed state. This yields a set of parameters at each of five electrode sites along the midline. The inferred parameter nonuniformities are consistent with anatomical and physiological constraints. Introducing these spatial profiles into the full nonuniform model then quantitatively reproduces observed site-dependent variations in the alpha power and frequency. The results confirm that the frequency shift is mainly due to a decrease in the corticothalamic propagation delay, but indicate that the delay nonuniformity cannot account for the observed occipital increase in alpha power; the occipital alpha dominance is due to decreased cortical gains and increased thalamic gains in occipital regions compared to frontal ones.

A disturbance of nonlinear interdependence in scalp EEG of subjects with first episode schizophrenia

It has been proposed that schizophrenia arises through a disturbance of coupling between large-scale cortical systems. This "disconnection hypothesis" is tested by applying a measure of dynamical interdependence to scalp EEG data. EEG data were collected from 40 subjects with a first episode of schizophrenia and 40 matched healthy controls. An algorithm for the detection of dynamical interdependence was applied to six pairs of bipolar electrodes in each subject. The topographic organization of the interdependence was calculated and served as the principle measure of cortical integration. The rate of occurrence of dynamical interdependence did not statistically differ between subject groups at any of the sites. However, the topography across the scalp was significantly different between the two groups. Specifically, nonlinear interdependence tended to occur in larger concurrent "clusters" across the scalp in schizophrenia than in the healthy subjects. This disturbance was reflected most strongly in left intrahemispheric coupling and did not differ significantly according to symptomatology. Medication dose and subject arousal were not observed to be confounding factors. The study of dynamical interdependence in scalp EEG data does not support a straightforward interpretation of the disconnection hypothesis-that there is a decrease in the strength of functional coupling between adjacent cortical regions. Rather, it suggests a dysregulation in the organization of dynamical interactions across supraregional brain systems.

Nonuniform corticothalamic continuum model of electroencephalographic spectra with application to split-alpha peaks

Recent theoretical work has successfully predicted electroencephalographic spectra from physiology using a model corticothalamic system with spatially uniform parameters. The present work incorporates parameter nonuniformities into this model via the coupling they induce between spatial eigenmodes. Splitting of the spectral alpha peak, an effect seen in a small percentage of the normal population, is investigated as an illustrative special case. It is confirmed that weak splitting can arise from mode structure if the peak is sufficiently sharp, even for uniform parameters. However, it is further demonstrated that greater splitting can result from nonuniformities, and it is argued that this mechanism for split alpha is better able to account quantitatively for this effect than previously suggested alternatives of pacemakers or purely cortical resonances. On introducing nonuniformities in corticothalamic loop time delays, we find that the alpha frequency also varies as one moves from the front to the back of the head, in accord with observations, and that analogous (but less distinct) variations are seen in the beta peak. Analysis shows realistic variations of around +/-10 ms relative to the mean loop delay of approximately 80 ms can account for observed splittings of about 1 Hz. It is also suggested that subjects who display clear alpha splitting form the tail of a distribution of magnitude of cortical inhomogeneity, rather than a separate population.

Neurophysical modeling of brain dynamics

A recent neurophysical model of brain electrical activity is outlined and applied to EEG phenomena. It incorporates single-neuron physiology and the large-scale anatomy of corticocortical and corticothalamic pathways, including synaptic strengths, dendritic propagation, nonlinear firing responses, and axonal conduction. Small perturbations from steady states account for observed EEGs as functions of arousal. Evoked response potentials (ERPs), correlation, and coherence functions are also reproduced. Feedback via thalamic nuclei is critical in determining the forms of these quantities, the transition between sleep and waking, and stability against seizures. Many disorders correspond to significant changes in EEGs, which can potentially be quantified in terms of the underlying physiology using this theory. In the nonlinear regime, limit cycles are often seen, including a regime in which they have the characteristic petit mal 3 Hz spike-and-wave form.

Wave-number spectrum of electrocorticographic signals

A physiologically based continuum model of corticothalamic electrodynamics is generalized and used to derive the theoretical form of the electrocorticographic (ECoG) wave-number spectrum. A one-dimensional projection of the spectrum is derived, as is the azimuthally averaged two-dimensional spectrum for isotropic and anisotropic cortices. The predicted spectra are found to consist of a low-k plateau followed by three regions of power-law decrease, which result from filtering of the electrical activity through physical structures at different scales in the cortex. The magnitude of the maximum theoretical power-law exponent is larger for the two-dimensional (2D) spectrum than for its 1D counterpart. The predicted spectra agree well with experimental data obtained from 1D and 2D recording arrays on the cortical surface, enabling the structures in the brain that are important in determining spatial cortical dynamics to be identified. The cortical dispersion relation predicted by our model is also investigated, providing insight into the relationships between temporal and spatial brain dynamics.

Interpretation of scaling properties of electroencephalographic fluctuations via spectral analysis and underlying physiology

Detrended fluctuation analysis has recently demonstrated the existence of two approximate temporal scaling regimes in locally detrended human electroencephalographic (EEG) fluctuations, and has suggested a connection between the location of the breakpoint between regimes and the alpha resonance near 10 Hz. It is shown here that these scalings can be explained in terms of the filtering of the underlying power spectrum implied by the detrending process. Using a recent physiologically based model of EEG generation, the main features of the scalings, and deviations from them, are related to the underlying physiology of dendritic propagation and muscle electrical activity, and it is concluded that the effects of such physiological features are usually clearer in spectra.

Wave-number spectrum of electroencephalographic signals

A recently developed, physiologically based continuum model of corticothalamic electrodynamics is used to derive the theoretical form of the electroencephalographic wave-number spectrum and its projection onto a one-dimensional recording array. The projected spectrum is found to consist of a plateau followed by regions of power-law decrease with various exponents, which are dependent on both model parameters and temporal frequency. The theoretical spectrum is compared with experimental results obtained in other studies, showing good agreement. The model provides a framework for understanding the nature of the spatial power spectrum by linking the underlying physiology with the large-scale dynamics of the brain.