Computational model prospective on the observation of proictal states in epileptic neuronal systems

Epileptic seizure suppression: A computational approach for identification and control using real data

Epilepsy affects millions of people worldwide every year and remains an open subject for research. Current development on this field has focused on obtaining computational models to better understand its triggering mechanisms, attain realistic descriptions and study seizure suppression. Controllers have been successfully applied to mitigate epileptiform activity in dynamic models written in state-space notation, whose applicability is, however, restricted to signatures that are accurately described by them. Alternatively, autoregressive modeling (AR), a typical data-driven tool related to system identification (SI), can be directly applied to signals to generate more realistic models, and since it is inherently convertible into state-space representation, it can thus be used for the artificial reconstruction and attenuation of seizures as well. Considering this, the first objective of this work is to propose an SI approach using AR models to describe real epileptiform activity. The second objective is to provide a strategy for reconstructing and mitigating such activity artificially, considering non-hybrid and hybrid controllers - designed from ictal and interictal events, respectively. The results show that AR models of relatively low order represent epileptiform activities fairly well and both controllers are effective in attenuating the undesired activity while simultaneously driving the signal to an interictal condition. These findings may lead to customized models based on each signal, brain region or patient, from which it is possible to better define shape, frequency and duration of external stimuli that are necessary to attenuate seizures.

Chaotic and stochastic dynamics of epileptiform-like activities in sclerotic hippocampus resected from patients with pharmacoresistant epilepsy

The types of epileptiform activity occurring in the sclerotic hippocampus with highest incidence are interictal-like events (II) and periodic ictal spiking (PIS). These activities are classified according to their event rates, but it is still unclear if these rate differences are consequences of underlying physiological mechanisms. Identifying new and more specific information related to these two activities may bring insights to a better understanding about the epileptogenic process and new diagnosis. We applied Poincaré map analysis and Recurrence Quantification Analysis (RQA) onto 35 in vitro electrophysiological signals recorded from slices of 12 hippocampal tissues surgically resected from patients with pharmacoresistant temporal lobe epilepsy. These analyzes showed that the II activity is related to chaotic dynamics, whereas the PIS activity is related to deterministic periodic dynamics. Additionally, it indicates that their different rates are consequence of different endogenous dynamics. Finally, by using two computational models we were able to simulate the transition between II and PIS activities. The RQA was applied to different periods of these simulations to compare the recurrences between artificial and real signals, showing that different ranges of regularity-chaoticity can be directly associated with the generation of PIS and II activities.

Temporal lobe epilepsy (TLE) is the most prevalent type of epilepsy in adults and hippocampal sclerosis is the major pathophysiological substrate of pharmaco-refractory TLE. Different patterns of epileptiform-like activity have been described in human hippocampal sclerosis, but the standard analysis applied to characterize the activities usually do not consider the nonlinear features that epileptiform patterns exhibit. Here, using Poincaré map and Recurrence Quantitative Analysis we characterized the most prevalent type of epileptiform-like activities—interictal-like events (II) and periodic ictal spiking (PIS), recorded in vitro from resected hippocampi of pharmacoresistant patients with TLE—according to their levels of stochasticity, chaoticity and determinism. The II activities showed to be more chaotic with complex rhythmicity than PIS activities. The nonlinear dynamic differences between II and PIS leads us to conjecture that they are expressions of different seizure susceptibility. We also identified that each hippocampal subfield expresses II and PIS activities in a specific and different way. Finally, from the modulation of internal parameters of two computational models, we show the conversion of one type of activity into the other, showing how specific neuron networks synchronize over time, leading to II and PIS activities and then into a generalized seizure.

Temporally Targeted Interactions With Pathologic Oscillations as Therapeutical Targets in Epilepsy and Beyond

Self-organized neuronal oscillations rely on precisely orchestrated ensemble activity in reverberating neuronal networks. Chronic, non-malignant disorders of the brain are often coupled to pathological neuronal activity patterns. In addition to the characteristic behavioral symptoms, these disturbances are giving rise to both transient and persistent changes of various brain rhythms. Increasing evidence support the causal role of these "oscillopathies" in the phenotypic emergence of the disease symptoms, identifying neuronal network oscillations as potential therapeutic targets. While the kinetics of pharmacological therapy is not suitable to compensate the disease related fine-scale disturbances of network oscillations, external biophysical modalities (e.g., electrical stimulation) can alter spike timing in a temporally precise manner. These perturbations can warp rhythmic oscillatory patterns via resonance or entrainment. Properly timed phasic stimuli can even switch between the stable states of networks acting as multistable oscillators, substantially changing the emergent oscillatory patterns. Novel transcranial electric stimulation (TES) approaches offer more reliable neuronal control by allowing higher intensities with tolerable side-effect profiles. This precise temporal steerability combined with the non- or minimally invasive nature of these novel TES interventions make them promising therapeutic candidates for functional disorders of the brain. Here we review the key experimental findings and theoretical background concerning various pathological aspects of neuronal network activity leading to the generation of epileptic seizures. The conceptual and practical state of the art of temporally targeted brain stimulation is discussed focusing on the prevention and early termination of epileptic seizures.

Coherent false seizure prediction in epilepsy, coincidence or providence?

OBJECTIVE

Seizure forecasting using machine learning is possible, but the performance is far from ideal, as indicated by many false predictions and low specificity. Here, we examine false and missing alarms of two algorithms on long-term datasets to show that the limitations are less related to classifiers or features, but rather to intrinsic changes in the data.

METHODS

We evaluated two algorithms on three datasets by computing the correlation of false predictions and estimating the information transfer between both classification methods.

RESULTS

For 9 out of 12 individuals both methods showed a performance better than chance. For all individuals we observed a positive correlation in predictions. For individuals with strong correlation in false predictions we were able to boost the performance of one method by excluding test samples based on the results of the second method.

CONCLUSIONS

Substantially different algorithms exhibit a highly consistent performance and a strong coherency in false and missing alarms. Hence, changing the underlying hypothesis of a preictal state of fixed time length prior to each seizure to a proictal state is more helpful than further optimizing classifiers.

SIGNIFICANCE

The outcome is significant for the evaluation of seizure prediction algorithms on continuous data.

Active probing to highlight approaching transitions to ictal states in coupled neural mass models

The extraction of electrophysiological features that reliably forecast the occurrence of seizures is one of the most challenging goals in epilepsy research. Among possible approaches to tackle this problem is the use of active probing paradigms in which responses to stimuli are used to detect underlying system changes leading up to seizures. This work evaluates the theoretical and mechanistic underpinnings of this strategy using two coupled populations of the well-studied Wendling neural mass model. Different model settings are evaluated, shifting parameters (excitability, slow inhibition, or inter-population coupling gains) from normal towards ictal states while probing stimuli are applied every 2 seconds to the input of either one or both populations. The correlation between the extracted features and the ictogenic parameter shifting indicates if the impending transition to the ictal state may be identified in advance. Results show that not only can the response to the probing stimuli forecast seizures but this is true regardless of the altered ictogenic parameter. That is, similar feature changes are highlighted by probing stimuli responses in advance of the seizure including: increased response variance and lag-1 autocorrelation, decreased skewness, and increased mutual information between the outputs of both model subsets. These changes were mostly restricted to the stimulated population, showing a local effect of this perturbational approach. The transition latencies from normal activity to sustained discharges of spikes were not affected, suggesting that stimuli had no pro-ictal effects. However, stimuli were found to elicit interictal-like spikes just before the transition to the ictal state. Furthermore, the observed feature changes highlighted by probing the neuronal populations may reflect the phenomenon of critical slowing down, where increased recovery times from perturbations may signal the loss of a systems' resilience and are common hallmarks of an impending critical transition. These results provide more evidence that active probing approaches highlight information about underlying system changes involved in ictogenesis and may be able to play a role in assisting seizure forecasting methods which can be incorporated into early-warning systems that ultimately enable closing the loop for targeted seizure-controlling interventions.

Neurostimulation stabilizes spiking neural networks by disrupting seizure-like oscillatory transitions

An improved understanding of the mechanisms underlying neuromodulatory approaches to mitigate seizure onset is needed to identify clinical targets for the treatment of epilepsy. Using a Wilson–Cowan-motivated network of inhibitory and excitatory populations, we examined the role played by intrinsic and extrinsic stimuli on the network’s predisposition to sudden transitions into oscillatory dynamics, similar to the transition to the seizure state. Our joint computational and mathematical analyses revealed that such stimuli, be they noisy or periodic in nature, exert a stabilizing influence on network responses, disrupting the development of such oscillations. Based on a combination of numerical simulations and mean-field analyses, our results suggest that high variance and/or high frequency stimulation waveforms can prevent multi-stability, a mathematical harbinger of sudden changes in network dynamics. By tuning the neurons’ responses to input, stimuli stabilize network dynamics away from these transitions. Furthermore, our research shows that such stabilization of neural activity occurs through a selective recruitment of inhibitory cells, providing a theoretical undergird for the known key role these cells play in both the healthy and diseased brain. Taken together, these findings provide new vistas on neuromodulatory approaches to stabilize neural microcircuit activity.

A taxonomy of seizure dynamotypes

Seizures are a disruption of normal brain activity present across a vast range of species and conditions. We introduce an organizing principle that leads to the first objective Taxonomy of Seizure Dynamics (TSD) based on bifurcation theory. The ‘dynamotype’ of a seizure is the dynamic composition that defines its observable characteristics, including how it starts, evolves and ends. Analyzing over 2000 focal-onset seizures from multiple centers, we find evidence of all 16 dynamotypes predicted in TSD. We demonstrate that patients’ dynamotypes evolve during their lifetime and display complex but systematic variations including hierarchy (certain types are more common), non-bijectivity (a patient may display multiple types) and pairing preference (multiple types may occur during one seizure). TSD provides a way to stratify patients in complement to present clinical classifications, a language to describe the most critical features of seizure dynamics, and a framework to guide future research focused on dynamical properties.

Epileptic seizures have been recognized for centuries. But it was only in the 1930s that it was realized that seizures are the result of out-of-control electrical activity in the brain. By placing electrodes on the scalp, doctors can identify when and where in the brain a seizure begins. But they cannot tell much about how the seizure behaves, that is, how it starts, stops or spreads to other areas. This makes it difficult to control and prevent seizures. It also helps explain why almost a third of patients with epilepsy continue to have seizures despite being on medication.

Saggio, Crisp et al. have now approached this problem from a new angle using methods adapted from physics and engineering. In these fields, “dynamics research” has been used with great success to predict and control the behavior of complex systems like electrical power grids. Saggio, Crisp et al. reasoned that applying the same approach to the brain would reveal the dynamics of seizures and that such information could then be used to categorize seizures into groups with similar properties. This would in effect create for seizures what the periodic table is for the elements.

Applying the dynamics research method to seizure data from more than a hundred patients from across the world revealed 16 types of seizure dynamics. These “dynamotypes” had distinct characteristics. Some were more common than others, and some tended to occur together. Individual patients showed different dynamotypes over time. By constructing a way to classify seizures based on the relationships between the dynamotypes, Saggio, Crisp et al. provide a new tool for clinicians and researchers studying epilepsy.

Previous clinical tools have focused on the physical symptoms of a seizure (referred to as the phenotype) or its potential genetic causes (genotype). The current approach complements these tools by adding the dynamotype: how seizures start, spread and stop in the brain. This approach has the potential to lead to new branches of research and better understanding and treatment of seizures.

Seizure prediction - ready for a new era

Epilepsy is a common disorder characterized by recurrent seizures. An overwhelming majority of people with epilepsy regard the unpredictability of seizures as a major issue. More than 30 years of international effort have been devoted to the prediction of seizures, aiming to remove the burden of unpredictability and to couple novel, time-specific treatment to seizure prediction technology. A highly influential review published in 2007 concluded that insufficient evidence indicated that seizures could be predicted. Since then, several advances have been made, including successful prospective seizure prediction using intracranial EEG in a small number of people in a trial of a real-time seizure prediction device. In this Review, we examine advances in the field, including EEG databases, seizure prediction competitions, the prospective trial mentioned and advances in our understanding of the mechanisms of seizures. We argue that these advances, together with statistical evaluations, set the stage for a resurgence in efforts towards the development of seizure prediction methodologies. We propose new avenues of investigation involving a synergy between mechanisms, models, data, devices and algorithms and refine the existing guidelines for the development of seizure prediction technology to instigate development of a solution that removes the burden of the unpredictability of seizures.

The role that choice of model plays in predictions for epilepsy surgery

Mathematical modelling has been widely used to predict the effects of perturbations to brain networks. An important example is epilepsy surgery, where the perturbation in question is the removal of brain tissue in order to render the patient free of seizures. Different dynamical models have been proposed to represent transitions to ictal states in this context. However, our choice of which mathematical model to use to address this question relies on making assumptions regarding the mechanism that defines the transition from background to the seizure state. Since these mechanisms are unknown, it is important to understand how predictions from alternative dynamical descriptions compare. Herein we evaluate to what extent three different dynamical models provide consistent predictions for the effect of removing nodes from networks. We show that for small, directed, connected networks the three considered models provide consistent predictions. For larger networks, predictions are shown to be less consistent. However consistency is higher in networks that have sufficiently large differences in ictogenicity between nodes. We further demonstrate that heterogeneity in ictogenicity across nodes correlates with variability in the number of connections for each node.

Preventive and Abortive Strategies for Stimulation Based Control of Epilepsy: A Computational Model Study

Epilepsy is a condition in which periods of ongoing normal EEG activity alternate with periods of oscillatory behavior characteristic of epileptic seizures. The dynamics of the transitions between the two states are still unclear. Computational models provide a powerful tool to explore the underlying mechanisms of such transitions, with the purpose of eventually finding therapeutic interventions for this debilitating condition. In this study, the possibility to postpone seizures elicited by a decrease of inhibition is investigated by using external stimulation in a realistic bistable neuronal model consisting of two interconnected neuronal populations representing pyramidal cells and interneurons. In the simulations, seizures are induced by slowly decreasing the conductivity of GABA[Formula: see text] synaptic channels over time. Since the model is bistable, the system will change state from the initial steady state (SS) to the limit cycle (LS) state because of internal noise, when the inhibition falls below a certain threshold. Several state-independent stimulations paradigms are simulated. Their effectiveness is analyzed for various stimulation frequencies and intensities in combination with periodic and random stimulation sequences. The distributions of the time to first seizure in the presence of stimulation are compared with the situation without stimulation. In addition, stimulation protocols targeted to specific subsystems are applied with the objective of counteracting the baseline shift due to decreased inhibition in the system. Furthermore, an analytical model is used to investigate the effects of random noise. The relation between the strength of random noise stimulation, the control parameter of the system and the transitions between steady state and limit cycle are investigated. The study shows that it is possible to postpone epileptic activity by targeted stimulation in a realistic neuronal model featuring bistability and that it is possible to stop seizures by random noise in an analytical model.

Estimation of brain network ictogenicity predicts outcome from epilepsy surgery

Surgery is a valuable option for pharmacologically intractable epilepsy. However, significant post-operative improvements are not always attained. This is due in part to our incomplete understanding of the seizure generating (ictogenic) capabilities of brain networks. Here we introduce an in silico, model-based framework to study the effects of surgery within ictogenic brain networks. We find that factors conventionally determining the region of tissue to resect, such as the location of focal brain lesions or the presence of epileptiform rhythms, do not necessarily predict the best resection strategy. We validate our framework by analysing electrocorticogram (ECoG) recordings from patients who have undergone epilepsy surgery. We find that when post-operative outcome is good, model predictions for optimal strategies align better with the actual surgery undertaken than when post-operative outcome is poor. Crucially, this allows the prediction of optimal surgical strategies and the provision of quantitative prognoses for patients undergoing epilepsy surgery.

Gap Junctions as Common Cause of High-Frequency Oscillations and Epileptic Seizures in a Computational Cascade of Neuronal Mass and Compartmental Modeling

High frequency oscillations (HFO) appear to be a promising marker for delineating the seizure onset zone (SOZ) in patients with localization related epilepsy. It remains, however, a purely observational phenomenon and no common mechanism has been proposed to relate HFOs and seizure generation. In this work we show that a cascade of two computational models, one on detailed compartmental scale and a second one on neural mass scale can explain both the autonomous generation of HFOs and the presence of epileptic seizures as emergent properties. To this end we introduce axonal–axonal gap junctions on a microscopic level and explore their impact on the higher level neural mass model (NMM). We show that the addition of gap junctions can generate HFOs and simultaneously shift the operational point of the NMM from a steady state network into bistable behavior that can autonomously generate epileptic seizures. The epileptic properties of the system, or the probability to generate epileptic type of activity, increases gradually with the increase of the density of axonal–axonal gap junctions. We further demonstrate that ad hoc HFO detectors used in previous studies are applicable to our simulated data.

Role of multiple-scale modeling of epilepsy in seizure forecasting

Over the past three decades, a number of seizure prediction, or forecasting, methods have been developed. Although major achievements were accomplished regarding the statistical evaluation of proposed algorithms, it is recognized that further progress is still necessary for clinical application in patients. The lack of physiological motivation can partly explain this limitation. Therefore, a natural question is raised: can computational models of epilepsy be used to improve these methods? Here, we review the literature on the multiple-scale neural modeling of epilepsy and the use of such models to infer physiologic changes underlying epilepsy and epileptic seizures. The authors argue how these methods can be applied to advance the state-of-the-art in seizure forecasting.

A critical role for network structure in seizure onset: a computational modeling approach

Recent clinical work has implicated network structure as critically important in the initiation of seizures in people with idiopathic generalized epilepsies. In line with this idea, functional networks derived from the electroencephalogram (EEG) at rest have been shown to be significantly different in people with generalized epilepsy compared to controls. In particular, the mean node degree of networks from the epilepsy cohort was found to be statistically significantly higher than those of controls. However, the mechanisms by which these network differences can support recurrent transitions into seizures remain unclear. In this study, we use a computational model of the transition into seizure dynamics to explore the dynamic consequences of these differences in functional networks. We demonstrate that networks with higher mean node degree are more prone to generating seizure dynamics in the model and therefore suggest a mechanism by which increased mean node degree of brain networks can cause heightened ictogenicity.

MULTIPLE OSCILLATORY STATES IN MODELS OF COLLECTIVE NEURONAL DYNAMICS

In our previous studies, we showed that the both realistic and analytical computational models of neural dynamics can display multiple sustained states (attractors) for the same values of model parameters. Some of these states can represent normal activity while other, of oscillatory nature, may represent epileptic types of activity. We also showed that a simplified, analytical model can mimic this type of behavior and can be used instead of the realistic model for large scale simulations. The primary objective of the present work is to further explore the phenomenon of multiple stable states, co-existing in the same operational model, or phase space, in systems consisting of large number of interconnected basic units. As a second goal, we aim to specify the optimal method for state control of the system based on inducing state transitions using appropriate external stimulus. We use here interconnected model units that represent the behavior of neuronal populations as an effective dynamic system. The model unit is an analytical model (S. Kalitzin et al., Epilepsy Behav. 22 (2011) S102–S109) and does not correspond directly to realistic neuronal processes (excitatory–inhibitory synaptic interactions, action potential generation). For certain parameter choices however it displays bistable dynamics imitating the behavior of realistic neural mass models. To analyze the collective behavior of the system we applied phase synchronization analysis (PSA), principal component analysis (PCA) and stability analysis using Lyapunov exponent (LE) estimation. We obtained a large variety of stable states with different dynamic characteristics, oscillatory modes and phase relations between the units. These states can be initiated by appropriate initial conditions; transitions between them can be induced stochastically by fluctuating variables (noise) or by specific inputs. We propose a method for optimal reactive control, allowing forced transitions from one state (attractor) into another.

DYNAMICS OF COLLECTIVE MULTI-STABILITY IN MODELS OF MULTI-UNIT NEURONAL SYSTEMS

In this study, we investigate the correspondence between dynamic patterns of behavior in two types of computational models of neuronal activity. The first model type is the realistic neuronal model; the second model type is the phenomenological or analytical model. In the simplest model set-up of two interconnected units, we define a parameter space for both types of systems where their behavior is similar. Next we expand the analytical model to two sets of 90 fully interconnected units with some overlap, which can display multi-stable behavior. This system can be in three classes of states: (i) a class consisting of a single resting state, where all units of a set are in steady state, (ii) a class consisting of multiple preserving states, where subsets of the units of a set participate in limit cycle, and (iii) a class consisting of a single saturated state, where all units of a set are recruited in a global limit cycle. In the third and final part of the work, we demonstrate that phase synchronization of units can be detected by a single output unit.

Connectivity and phase coherence in neural network models of interconnected Z(4)-bi-stable units

A phenomenological neural network model with bi-stable oscillatory units is used to model up- and down-states. These states have been observed in vivo in biological neuronal systems and feature oscillatory, limit cycle type of behavior in the up-states. A network is formed by a set of interconnected units. Two different types of network layouts are considered in this work: networks with hierarchical connections and hubs and networks with random connections. The phase coherence between the different units is analyzed and compared to the connectivity distance between nodes. In addition the connectivity degree of a node is associated to the average phase coherence with all other units. The results show that we may be able to identify the set of hubs in a network based on the phase coherence estimates between the different nodes. If the network is very dense or randomly connected, the underlying network structure, however, can not be derived uniquely from the phase coherence.

REACTIVE CONTROL OF EPILEPTIFORM DISCHARGES IN REALISTIC COMPUTATIONAL NEURONAL MODELS WITH BISTABILITY

We aim to derive fully autonomous seizure suppression paradigms based on reactive control of neuronal dynamics. A previously derived computational model of seizure generation describing collective degrees of freedom and featuring bistable dynamics is used. A novel technique for real-time control of epileptogenicity is introduced. The reactive control reduces practically all seizures in the model. The study indicates which parameters provide the maximal seizure reduction with minimal intervention. An adaptive scheme is proposed that optimizes the stimulation parameters in nonstationary situations.

Quantification of spontaneous and evoked HFO's in SEEG recording and prospective for pre-surgical diagnostics. Case study

High frequency oscillations (HFO) in stereo electroencephalographic (SEEG) signals have been recently the focus of attention as biomarkers that can have potential predictive power for the spatial location and possibly the timing of the onset of epileptic seizures. In this work we present a case study where we compare two quantitative paradigms for automated detection of biomarkers, one based on spontaneous SEEG recordings of HFOs and the other using activity induced by direct electrical stimulation (relative Phase Clustering Index algorithm). We compare the performance of these automated methods with manually detected HFO ripples by a trained EEG analyst and explore their potential diagnostic relevance. Intracranial recordings from patients undergoing pre-surgical evaluation are processed with a combination of morphological filtering and the analysis of the auto-correlation function. The results were compared to those obtained by visual inspection and to results from an active paradigm involving stimulation with 20 Hz trains of biphasic pulses. The quantity of HFOs, estimated automatically, or "rippleness", was found to correspond to the findings of a trained EEG analyst. The relative phase clustering index (rPCI) obtained using periodic stimulation appeared to be associated with the closeness to the seizure onset zone (SOZ) detected from ictal epochs. The HFO estimates were also indicative for the SOZ but with less specificity.