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.

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.

Biochemical energy is the fundamental element that maintains both the adequate turnover of the biomolecular structures and the functional metabolic viability of unicellular organisms. The levels of ATP, ADP and AMP reflect roughly the energetic status of the cell, and a precise ratio relating them was proposed by Atkinson as the adenylate energy charge (AEC). Under growth-phase conditions, cells maintain the AEC within narrow physiological values, despite extremely large fluctuations in the adenine nucleotides concentration. Intensive experimental studies have shown that these AEC values are preserved in a wide variety of organisms, both eukaryotes and prokaryotes. Here, to understand some of the functional elements involved in the cellular energy status, we present a computational model conformed by some key essential parts of the adenylate energy system. Specifically, we have considered (I) the main synthesis process of ATP from ADP, (II) the main catalyzed phosphotransfer reaction for interconversion of ATP, ADP and AMP, (III) the enzymatic hydrolysis of ATP yielding ADP, and (IV) the enzymatic hydrolysis of ATP providing AMP. This leads to a dynamic metabolic model (with the form of a delayed differential system) in which the enzymatic rate equations and all the physiological kinetic parameters have been explicitly considered and experimentally tested in vitro. Our central hypothesis is that cells are characterized by changing energy dynamics (homeorhesis). The results show that the AEC presents stable transitions between steady states and periodic oscillations and, in agreement with experimental data these oscillations range within the narrow AEC window. Furthermore, the model shows sustained oscillations in the Gibbs free energy and in the total nucleotide pool. The present study provides a step forward towards the understanding of the fundamental principles and quantitative laws governing the adenylate energy system, which is a fundamental element for unveiling the dynamics of cellular life.

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.

We examine the properties of the transfer function F T = V m /V LFP between the intracellular membrane potential (V m ) and the local field potential (V LFP ) in cerebral cortex. We first show theoretically that, in the subthreshold regime, the frequency dependence of the extracellular medium and that of the membrane potential have a clear incidence on F T . The calculation of F T from experiments and the matching with theoretical expressions is possible for desynchronized states where individual current sources can be considered as independent. Using a mean-field approximation, we obtain a method to estimate the impedance of the extracellular medium without injecting currents. We examine the transfer function for bipolar (differential) LFPs and compare to simultaneous recordings of V m and V LFP during desynchronized states in rat barrel cortex in vivo. The experimentally derived F T matches the one derived theoretically, only if one assumes that the impedance of the extracellular medium is frequency-dependent, and varies as 1/ √ ω (Warburg impedance) for frequencies between 3 and 500 Hz. This constitutes indirect evidence that the extracellular medium is non-resistive, which has many possible consequences for modeling LFPs.