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.

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.

Background Alzheimer’s disease (AD) is the most common neurodegenerative disease ultimately manifesting as clinical dementia. Despite considerable effort and ample experimental data, the role of neuroinflammation related to systemic inflammation is still unsettled. While the implication of microglia is well recognized, the exact contribution of peripheral monocytes/macrophages is still largely unknown, especially concerning their role in the various stages of AD. Objectives AD develops over decades and its clinical manifestation is preceded by subjective memory complaints (SMC) and mild cognitive impairment (MCI); thus, the question arises how the peripheral innate immune response changes with the progression of the disease. Therefore, to further investigate the roles of monocytes/macrophages in the progression of AD we assessed their phenotypes and functions in patients at SMC, MCI and AD stages and compared them with cognitively healthy controls. We also conceptualised an idealised mathematical model to explain the functionality of monocytes/macrophages along the progression of the disease. Results We show that there are distinct phenotypic and functional changes in monocyte and macrophage populations as the disease progresses. Higher free radical production upon stimulation could already be observed for the monocytes of SMC patients. The most striking results show that activation of peripheral monocytes (hyperactivation) is the strongest in the MCI group, at the prodromal stage of the disease. Monocytes exhibit significantly increased chemotaxis, free radical production, and cytokine production in response to TLR2 and TLR4 stimulation. Conclusion Our data suggest that the peripheral innate immune system is activated during the progression from SMC through MCI to AD, with the highest levels of activation being in MCI subjects and the lowest in AD patients. Some of these parameters may be used as biomarkers, but more holistic immune studies are needed to find the best period of the disease for clinical intervention.

Canard-induced phenomena have been extensively studied in the last three decades, from both the mathematical and the application viewpoints. Canards in slow-fast systems with (at least) two slow variables, especially near folded-node singularities, give an essential generating mechanism for mixed-mode oscillations (MMOs) in the framework of smooth multiple timescale systems. There is a wealth of literature on such slow-fast dynamical systems and many models displaying canard-induced MMOs, particularly in neuroscience. In parallel, since the late 1990s several papers have shown that the canard phenomenon can be faithfully reproduced with piecewise-linear (PWL) systems in two dimensions, although very few results are available in the three-dimensional case. The present paper aims to bridge this gap by analyzing canonical PWL systems that display folded singularities, primary and secondary canards, with a similar control of the maximal winding number as in the smooth case. We also show that the singular phase portraits are compatible in both frameworks. Finally, we show using an example how to construct a (linear) global return and obtain robust PWL MMOs.Ministerio de Ciencia y Tecnología MTM2012-31821Junta de Andalucía P12-FQM-165