Summary:Purpose: The occurrence of abnormal dynamics in a physiological system can become manifest as a sudden qualitative change in the behavior of characteristic physiologic variables. We assume that this is what happens in the brain with regard to epilepsy. We consider that neuronal networks involved in epilepsy possess multistable dynamics (i.e., they may display several dynamic states). To illustrate this concept, we may assume, for simplicity, that at least two states are possible: an interictal one characterized by a normal, apparently random, steady-state of ongoing activity, and another one that is characterized by the paroxysmal occurrence of a synchronous oscillations (seizure).Methods: By using the terminology of the mathematics of nonlinear systems, we can say that such a bistable system has two attractors, to which the trajectories describing the system's output converge, depending on initial conditions and on the system's parameters. In phase-space, the basins of attraction corresponding to the two states are separated by what is called a "separatrix." We propose, schematically, that the transition between the normal ongoing and the seizure activity can take place according to three basic models:Model I: In certain epileptic brains (e.g., in absence seizures of idiopathic primary generalized epilepsies), the distance between "normal steady-state" and "paroxysmal" attractors is very small in contrast to that of a normal brain (possibly due to genetic and/or developmental factors). In the former, discrete random fluctuations of some variables can be sufficient for the occurrence of a transition to the paroxysmal state. In this case, such seizures are not predictable.Model II and model III: In other kinds of epileptic brains (e.g., limbic cortex epilepsies), the distance between "normal steadystate" and "paroxysmal" attractors is, in general, rather large, such that random fluctuations, of themselves, are commonly not capable of triggering a seizure. However, in these brains, neuronal networks have abnormal features characterized by unstable parameters that are very vulnerable to the influence of endogenous (model II) and/or exogenous (model III) factors. In these cases, these critical parameters may gradually change with time, in such a way that the attractor can deform either gradually or suddenly, with the consequence that the distance between the basin of attraction of the normal state and the separatrix tends to zero. This can lead, eventually, to a transition to a seizure.Results: The changes of the system's dynamics preceding a seizure in these models either may be detectable in the EEG and thus the route to the seizure may be predictable, or may be unobservable by using only measurements of the dynamical state. It is thinkable, however, that in some cases, changes in the excitability state of the underlying networks may be uncovered by using appropriate stimuli configurations before changes in the dynamics of the ongoing EEG activity are evident. A typical example of model III that we discuss...

No abstract

No abstract

No abstract

In this overview, we consider epilepsies as dynamical diseases of brain systems since they are manifestations of the property of neuronal networks to display multistable dynamics. To illustrate this concept we may assume that at least two states of the epileptic brain are possible: the interictal state characterized by a normal, apparently random, steady-state electroencephalography (EEG) ongoing activity, and the ictal state, that is characterized by paroxysmal occurrence of synchronous oscillations and is generally called, in neurology, a seizure. The transition between these two states can either occur: 1) as a continuous sequence of phases, like in some cases of mesial temporal lobe epilepsy (MTLE); or 2) as a sudden leap, like in most cases of absence seizures. In the mathematical terminology of nonlinear systems, we can say that in the first case the system's attractor gradually deforms from an interictal to an ictal attractor. The causes for such a deformation can be either endogenous or external. In this type of ictal transition, the seizure possibly may be anticipated in its early, preclinical phases. In the second case, where a sharp critical transition takes place, we can assume that the system has at least two simultaneous interictal and ictal attractors all the time. To which attractor the trajectories converge, depends on the initial conditions and the system's parameters. An essential question in this scenario is how the transition between the normal ongoing and the seizure activity takes place. Such a transition can occur either due to the influence of external or endogenous factors or due to a random perturbation and, thus, it will be unpredictable. These dynamical changes may not be detectable from the analysis of the ongoing EEG, but they may be observable only by measuring the system's response to externally administered stimuli. In the special cases of reflex epilepsy, the leap between the normal ongoing attractor and the ictal attractor is caused by a well-defined external perturbation. Examples from these different scenarios are presented and discussed.

In this paper, we investigate the dynamical scenarios of transitions between normal and paroxysmal state in epilepsy. We assume that some epileptic neural network are bistable i.e., they feature two operational states, ictal and interictal that co-exist. The transitions between these two states may occur according to a Poisson process, a random walk process or as a result of deterministic time-dependent mechanisms. We analyze data from animal models of absence epilepsy, human epilepsies and in vitro models. The distributions of durations of ictal and interictal epochs are fitted with a gamma distribution. On the basis of qualitative features of the fits, we identify the dynamical processes that may have generated the underlying data. The analysis showed that the following hold. 1) The dynamics of ictal epochs differ from those of interictal states. 2) Seizure initiation can be accounted for by a random walk process while seizure termination is often mediated by deterministic mechanisms. 3) In certain cases, the transitions between ictal and interictal states can be modeled by a Poisson process operating in a bistable network. These results imply that exact prediction of seizure occurrence is not possible but termination of an ictal state by appropriate counter stimulation might be feasible.

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.

Steady state visual evoked potentials (SSVEPs) are steady state oscillatory potentials elicited in the electroencephalogram (EEG) by flicker stimulation. The frequency of these responses maches the frequency of the stimulation and of its harmonics and subharmonics. In this study, we investigated the origin of the harmonic and subharmonic components of SSVEPs, which are not well understood. We applied both sine and square wave visual stimulation at 5 and 15 Hz to human subjects and analyzed the properties of the fundamental responses and harmonically related components. In order to interpret the results, we used the well-established neural mass model that consists of interacting populations of excitatory and inhibitory cortical neurons. In our study, this model provided a simple explanation for the origin of SSVEP spectra, and showed that their harmonic and subharmonic components are a natural consequence of the nonlinear properties of neuronal populations and the resonant properties of the modeled network. The model also predicted multiples of subharmonic responses, which were subsequently confirmed using experimental data.