A spatially extended model for macroscopic spike-wave discharges

Bifurcations and bursting in the Epileptor

The Epileptor is a phenomenological model for seizure activity that is used in a personalized large-scale brain modeling framework, the Virtual Epileptic Patient, with the aim of improving surgery outcomes for drug-resistant epileptic patients. Transitions between interictal and ictal states are modeled as bifurcations, enabling the definition of seizure classes in terms of onset/offset bifurcations. This establishes a taxonomy of seizures grounded in their essential underlying dynamics and the Epileptor replicates the activity of the most common class, as observed in patients with focal epilepsy, which is characterized by square-wave bursting properties. The Epileptor also encodes an additional mechanism to account for interictal spikes and spike and wave discharges. Here we use insights from a more generic model for square-wave bursting, based on the Unfolding Theory approach, to guide the bifurcation analysis of the Epileptor and gain a deeper understanding of the model and the role of its parameters. We show how the Epileptor's parameters can be modified to produce activities for other seizures classes of the taxonomy, as observed in patients, so that the large-scale brain models could be further personalized. Some of these classes have already been described in the literature in the Epileptor, others, predicted by the generic model, are new. Finally, we unveil how the interaction with the additional mechanism for spike and wave discharges alters the bifurcation structure of the main burster.

Dynamical Mechanism Analysis of Three Neuroregulatory Strategies on the Modulation of Seizures

This paper attempts to explore and compare the regulatory mechanisms of optogenetic stimulation (OS), deep brain stimulation (DBS) and electromagnetic induction on epilepsy. Based on the Wilson-Cowan model, we first demonstrate that the external input received by excitatory and inhibitory neural populations can induce rich dynamic bifurcation behaviors such as Hopf bifurcation, and make the system exhibit epileptic and normal states. Then, both OS and DBS are shown to be effective in controlling the epileptic state to a normal low-level state, and the stimulus parameters have a broad effective range. However, electromagnetic induction cannot directly control epilepsy to this desired state, even if it can significantly reduce the oscillation frequency of neural populations. One main difference worth noting is that the high spatiotemporal specificity of OS allows it to target inhibitory neuronal populations, whereas DBS and electromagnetic induction can only stimulate excitatory as well as inhibitory neuronal populations together. Next, the propagation behavior of epilepsy is explored under a typical three-node feedback loop structure. An increase in coupling strength accelerates and exacerbates epileptic activity in other brain regions. Finally, OS and DBS applied to the epileptic focus play similar positive roles in controlling the behavior of the area of seizure propagation, while electromagnetic induction still only achieves unsatisfactory effects. It is hoped that these dynamical results can provide insights into the treatment of epilepsy as well as other neurological disorders.

Mechanism analysis for excitatory interneurons dominating poly-spike wave and optimization of electrical stimulation

In addition to inhibitory interneurons, there exist excitatory interneurons (EINs) in the cortex, which mainly have excitatory projections to pyramidal neurons. In this study, we improve a thalamocortical model by introducing EIN, investigate the dominant role of EIN in generating spike and slow wave discharges (SWDs), and consider a non-rectangular pulse to control absence seizures. First, we display here that the improved model can reproduce typical SWDs of absence seizures. Moreover, we focus on the function of EIN by means of bifurcation analysis and find that EIN can induce transition behaviors under Hopf-type and fold limit cycle bifurcations. Specifically, the system has three stable solutions composing a tri-stable region. In this region, there are three attraction basins, which hints that external stimulation can drive the system trajectory from one basin to another, thereby eliminating abnormal oscillations. Furthermore, we compare the increasing ramp with rectangular pulse and optimize stimulation waveforms from the perspective of electrical charges input. The controlling role of the single increasing ramp to absence seizures is remarkable and the optimal stimulus parameters have been found theoretically. This work provides a computational model containing EIN and a theoretical basis for future physiological experiments and clinical research studies.

Regulatory Mechanism for Absence Seizures in Bidirectional Interactive Thalamocortical Model via Different Targeted Therapy Schemes

Recent clinical practice has found that the spike-wave discharge (SWD) scopes of absence seizures change from small cortical region to large thalamocortical networks, which has also been proved by theoretical simulation. The best biophysics explanation is that there are interactions between coupled cortico-thalamic and thalamocortical circuits. To agree with experiment results and describe the phenomena better, we constructed a coupled thalamocortical model with bidirectional channel (CTMBC) to account for the causes of absence seizures which are connected by the principle of two-way communication of neural pathways. By adjusting the coupling strength of bidirectional pathways, the spike-wave discharges are reproduced. Regulatory mechanism for absence seizures is further applied to CTMBC via four different targeted therapy schemes, such as deep brain stimulation (DBS), charge-balanced biphasic pulse (CBBP), coordinated reset stimulation (CRS) 1 : 0, and (CRS) 3 : 2. The new CTMBC model shows that neurodiversity in bidirectional interactive channel could supply theory reference for the bidirectional communication mode of thalamocortical networks and the hypothesis validation of pathogenesis.

Perturbations both trigger and delay seizures due to generic properties of slow-fast relaxation oscillators

The mechanisms underlying the emergence of seizures are one of the most important unresolved issues in epilepsy research. In this paper, we study how perturbations, exogenous or endogenous, may promote or delay seizure emergence. To this aim, due to the increasingly adopted view of epileptic dynamics in terms of slow-fast systems, we perform a theoretical analysis of the phase response of a generic relaxation oscillator. As relaxation oscillators are effectively bistable systems at the fast time scale, it is intuitive that perturbations of the non-seizing state with a suitable direction and amplitude may cause an immediate transition to seizure. By contrast, and perhaps less intuitively, smaller amplitude perturbations have been found to delay the spontaneous seizure initiation. By studying the isochrons of relaxation oscillators, we show that this is a generic phenomenon, with the size of such delay depending on the slow flow component. Therefore, depending on perturbation amplitudes, frequency and timing, a train of perturbations causes an occurrence increase, decrease or complete suppression of seizures. This dependence lends itself to analysis and mechanistic understanding through methods outlined in this paper. We illustrate this methodology by computing the isochrons, phase response curves and the response to perturbations in several epileptic models possessing different slow vector fields. While our theoretical results are applicable to any planar relaxation oscillator, in the motivating context of epilepsy they elucidate mechanisms of triggering and abating seizures, thus suggesting stimulation strategies with effects ranging from mere delaying to full suppression of seizures.

Despite its simplicity, the modelling of epileptic dynamics as a slow-fast transition between low and high activity states mediated by some slow feedback variable is a relatively novel albeit fruitful approach. This study is the first, to our knowledge, characterizing the response of such slow-fast models of epileptic brain to perturbations by computing its isochrons. Besides its numerical computation, we theoretically determine which factors shape the geometry of isochrons for planar slow-fast oscillators. As a consequence, we introduce a theoretical approach providing a clear understanding of the response of perturbations of slow-fast oscillators. Within the epilepsy context, this elucidates the origin of the contradictory role of interictal epileptiform discharges in the transition to seizure, manifested by both pro-convulsive and anti-convulsive effect depending on the amplitude, frequency and timing. More generally, this paper provides theoretical framework highlighting the role of the slow flow component on the response to perturbations in relaxation oscillators, pointing to the general phenomena in such slow-fast oscillators ubiquitous in biological systems.

Extracting the transition network of epileptic seizure onset

In presurgical monitoring, focal seizure onset is visually assessed from intracranial electroencephalogram (EEG), typically based on the selection of channels that show the strongest changes in amplitude and frequency. As epileptic seizure dynamics is increasingly considered to reflect changes in potentially distributed neural networks, it becomes important to also assess the interrelationships between channels. We propose a workflow to quantitatively extract the nodes and edges contributing to the seizure onset using an across-seizure scoring. We propose a quantification of the consistency of EEG channel contributions to seizure onset within a patient. The workflow is exemplified using recordings from patients with different degrees of seizure-onset consistency.

Spontaneous transitions to focal-onset epileptic seizures: A dynamical study

Given the complex temporal evolution of epileptic seizures, understanding their dynamic nature might be beneficial for clinical diagnosis and treatment. Yet, the mechanisms behind, for instance, the onset of seizures are still unknown. According to an existing classification, two basic types of dynamic onset patterns plus a number of more complex onset waveforms can be distinguished. Here, we introduce a basic three-variable model with two time scales to study potential mechanisms of spontaneous seizure onset. We expand the model to demonstrate how coupling of oscillators leads to more complex seizure onset waveforms. Finally, we test the response to pulse perturbation as a potential biomarker of interictal changes.

Divisive gain modulation enables flexible and rapid entrainment in a neocortical microcircuit model

Neocortical circuits exhibit a rich dynamic repertoire, and their ability to achieve entrainment (adjustment of their frequency to match the input frequency) is thought to support many cognitive functions and indicate functional flexibility. Although previous studies have explored the influence of various circuit properties on this phenomenon, the role of divisive gain modulation (or divisive inhibition) is unknown. This gain control mechanism is thought to be delivered mainly by the soma-targeting interneurons in neocortical microcircuits. In this study, we use a neural mass model of the neocortical microcircuit (extended Wilson-Cowan model) featuring both soma-targeting and dendrite-targeting interneuronal subpopulations to investigate the role of divisive gain modulation in entrainment. Our results demonstrate that the presence of divisive inhibition in the microcircuit, as delivered by the soma-targeting interneurons, enables its entrainment to a wider range of input frequencies. Divisive inhibition also promotes a faster entrainment, with the microcircuit needing less time to converge to the fully entrained state. We suggest that divisive inhibition, working alongside subtractive inhibition, allows for more adaptive oscillatory responses in neocortical circuits and, thus, supports healthy brain functioning.NEW & NOTEWORTHY We introduce a computational neocortical microcircuit model that features two inhibitory neural populations, with one providing subtractive and the other divisive inhibition to the excitatory population. We demonstrate that divisive inhibition widens the range of input frequencies to which the microcircuit can become entrained and diminishes the time needed to reach full entrainment. We suggest that divisive inhibition enables more adaptive oscillatory activity, with important implications for both normal and pathological brain function.

Unified analysis of global and focal aspects of absence epilepsy via neural field theory of the corticothalamic system

Absence epilepsy is characterized by a sudden paroxysmal loss of consciousness accompanied by oscillatory activity propagating over many brain areas. Although primary generalized absence seizures are supported by the global corticothalamic system, converging experimental evidence supports a focal theory of absence epilepsy. Here a physiology-based corticothalamic model is investigated with spatial heterogeneity due to focal epilepsy to unify global and focal aspects of absence epilepsy. Numeric and analytic calculations are employed to investigate the emergent spatiotemporal dynamics as well as their underlying dynamical mechanisms. They can be categorized into three scenarios: suppressed epilepsy, focal seizures, or generalized seizures, as summarized from a phase diagram vs focal width and characteristic axon range. The corresponding temporal frequencies and spatial extents of cortical waves in generalized seizures and focal seizures agree well with experimental observations of global and focal aspects of absence epilepsy, respectively. The emergence of the spatiotemporal dynamics corresponding to focal seizures provides a biophysical explanation of the temporally higher frequency but spatially more localized cortical waves observed in genetic rat models that display characteristics of human absence epilepsy. Predictions are also presented for further experimental test.

Dynamic origin of spike and wave discharges in the brain

Spike and wave discharges are the main electrographic characteristic of a number of epileptic brain disorders including childhood absence epilepsy and photosensitive epilepsy. The basic dynamic mechanism that underlies the occurrence of these abnormal electrical patterns in the brain is not well understood. The current paper aims to provide a dynamic explanation for features and generation mechanism of spike and wave discharges in the brain. The main proposition of this study is that epileptic seizures could be interpreted as a resonance phenomenon rather than a limit cycle behavior. To shows this, a revised version of Jansen-Rit neural mass model is employed. The system can switch between monostable and bistable regimes, which are considered in this paper as wake and sleep states of the brain, respectively. In particular, it is shown that, in monostable region, the model can depict the alpha rhythm and alpha rhythm suppression due to mental activity. Frequency responses of the model near the bistable regime demonstrate that high amplitude harmonic excitation may lead to spike and wave like oscillations. Based on the computational results and the concept of stochastic resonance, a model for absence epilepsy is presented which can simulate spontaneous transitions between ictal and interictal states. Finally, it is shown that spike and wave discharges during epileptic seizures can be explained as a resonance phenomenon in a nonlinear system.

Predicting seizure by modeling synaptic plasticity based on EEG signals - a case study of inherited epilepsy

This paper explores the internal dynamical mechanisms of epileptic seizures through quantitative modeling based on full brain electroencephalogram (EEG) signals. Our goal is to provide seizure prediction and facilitate treatment for epileptic patients. Motivated by an earlier mathematical model with incorporated synaptic plasticity, we studied the nonlinear dynamics of inherited seizures through a differential equation model. First, driven by a set of clinical inherited electroencephalogram data recorded from a patient with diagnosed Glucose Transporter Deficiency, we developed a dynamic seizure model on a system of ordinary differential equations. The model was reduced in complexity after considering and removing redundancy of each EEG channel. Then we verified that the proposed model produces qualitatively relevant behavior which matches the basic experimental observations of inherited seizure, including synchronization index and frequency. Meanwhile, the rationality of the connectivity structure hypothesis in the modeling process was verified. Further, through varying the threshold condition and excitation strength of synaptic plasticity, we elucidated the effect of synaptic plasticity to our seizure model. Results suggest that synaptic plasticity has great effect on the duration of seizure activities, which support the plausibility of therapeutic interventions for seizure control.

A new description of epileptic seizures based on dynamic analysis of a thalamocortical model

Increasing evidence suggests that the brain dynamics can be interpreted from the viewpoint of nonlinear dynamical systems. The aim of this paper is to investigate the behavior of a thalamocortical model from this perspective. The model includes both cortical and sensory inputs that can affect the dynamic nature of the model. Driving response of the model subjected to various harmonic stimulations is considered to identify the effects of stimulus parameters on the cortical output. Detailed numerical studies including phase portraits, Poincare maps and bifurcation diagrams reveal a wide range of complex dynamics including period doubling and chaos in the output. Transition between different states can occur as the stimulation parameters are changed. In addition, the amplitude jump phenomena and hysteresis are shown to be possible as a result of the bending in the frequency response curve. These results suggest that the jump phenomenon due to the brain nonlinear resonance can be responsible for the transitions between ictal and interictal states.

Stimulus-induced transitions between spike-wave discharges and spindles with the modulation of thalamic reticular nucleus

It is believed that thalamic reticular nucleus (TRN) controls spindles and spike-wave discharges (SWD) in seizure or sleeping processes. The dynamical mechanisms of spatiotemporal evolutions between these two types of activity, however, are not well understood. In light of this, we first use a single-compartment thalamocortical neural field model to investigate the effects of TRN on occurrence of SWD and its transition. Results show that the increasing inhibition from TRN to specific relay nuclei (SRN) can lead to the transition of system from SWD to slow-wave oscillation. Specially, it is shown that stimulations applied in the cortical neuronal populations can also initiate the SWD and slow-wave oscillation from the resting states under the typical inhibitory intensity from TRN to SRN. Then, we expand into a 3-compartment coupled thalamocortical model network in linear and circular structures, respectively, to explore the spatiotemporal evolutions of wave states in different compartments. The main results are: (i) for the open-ended model network, SWD induced by stimulus in the first compartment can be transformed into sleep-like slow UP-DOWN and spindle states as it propagates into the downstream compartments; (ii) for the close-ended model network, weak stimulations performed in the first compartment can result in the consistent experimentally observed spindle oscillations in all three compartments; in contrast, stronger periodic single-pulse stimulations applied in the first compartment can induce periodic transitions between SWD and spindle oscillations. Detailed investigations reveal that multi-attractor coexistence mechanism composed of SWD, spindles and background state underlies these state evolutions. What's more, in order to demonstrate the state evolution stability with respect to the topological structures of neural network, we further expand the 3-compartment coupled network into 10-compartment coupled one, with linear and circular structures, and nearest-neighbor (NN) coupled network as well as its realization of small-world (SW) topology via random rewiring, respectively. Interestingly, for the cases of linear and circular connetivities, qualitatively similar results were obtained in addition to the more irregularity of firings. However, SWD can be eventually transformed into the consistent low-amplitude oscillations for both NN and SW networks. In particular, SWD evolves into the slow spindling oscillations and background tonic oscillations within the NN and SW network, respectively. Our modeling and simulation studies highlight the effect of network topology in the evolutions of SWD and spindling oscillations, which provides new insights into the mechanisms of cortical seizures development.

The pacemaker role of thalamic reticular nucleus in controlling spike-wave discharges and spindles

Absence epilepsy, characterized by 2-4 Hz spike-wave discharges (SWDs), can be caused by pathological interactions within the thalamocortical system. Cortical spindling oscillations are also demonstrated to involve the oscillatory thalamocortical rhythms generated by the synaptic circuitry of the thalamus and cortex. This implies that SWDs and spindling oscillations can share the common thalamocortical mechanism. Additionally, the thalamic reticular nucleus (RE) is hypothesized to regulate the onsets and propagations of both the epileptic SWDs and sleep spindles. Based on the proposed single-compartment thalamocortical neural field model, we firstly investigate the stimulation effect of RE on the initiations, terminations, and transitions of SWDs. It is shown that the activations and deactivations of RE triggered by single-pulse stimuli can drive the cortical subsystem to behave as the experimentally observed onsets and self-abatements of SWDs, as well as the transitions from 2-spike and wave discharges (2-SWDs) to SWDs. In particular, with increasing inhibition from RE to the specific relay nucleus (TC), rich transition behaviors in cortex can be obtained through the upstream projection path, RE→TC→Cortex. Although some of the complex dynamical patterns can be expected from the earlier single compartment thalamocortical model, the effect of brain network topology on the emergence of SWDs and spindles, as well as the transitions between them, has not been fully investigated. We thereby develop a spatially extended 3-compartment coupled network model with open-/closed-end connective configurations, to investigate the spatiotemporal effect of RE on the SWDs and spindles. Results show that the degrees of activations of RE₁ can induce the rich spatiotemporal evolution properties including the propagations from SWDs to spindles within different compartments and the transitions between them, through the RE₁→TC₁→Cortex₁ and Cortex₁→Cortex₂→Cortex₃ projecting paths, respectively. Overall, those results imply that RE possesses the pacemaker function in controlling SWDs and spindling oscillations, which computationally provide causal support for the involvement of RE in absence seizures and sleep spindles.

Prediction of the Seizure Suppression Effect by Electrical Stimulation via a Computational Modeling Approach

In this paper, we identified factors that can affect seizure suppression via electrical stimulation by an integrative study based on experimental and computational approach. Preferentially, we analyzed the characteristics of seizure-like events (SLEs) using our previous in vitro experimental data. The results were analyzed in two groups classified according to the size of the effective region, in which the SLE was able to be completely suppressed by local stimulation. However, no significant differences were found between these two groups in terms of signal features or propagation characteristics (i.e., propagation delays, frequency spectrum, and phase synchrony). Thus, we further investigated important factors using a computational model that was capable of evaluating specific influences on effective region size. In the proposed model, signal transmission between neurons was based on two different mechanisms: synaptic transmission and the electrical field effect. We were able to induce SLEs having similar characteristics with differentially weighted adjustments for the two transmission methods in various noise environments. Although the SLEs had similar characteristics, their suppression effects differed. First of all, the suppression effect occurred only locally where directly received the stimulation effect in the high noise environment, but it occurred in the entire network in the low noise environment. Interestingly, in the same noise environment, the suppression effect was different depending on SLE propagation mechanism; only a local suppression effect was observed when the influence of the electrical field transmission was very weak, whereas a global effect was observed with a stronger electrical field effect. These results indicate that neuronal activities synchronized by a strong electrical field effect respond more sensitively to partial changes in the entire network. In addition, the proposed model was able to predict that stimulation of a seizure focus region is more effective for suppression. In conclusion, we confirmed the possibility of a computational model as a simulation tool to analyze the efficacy of deep brain stimulation (DBS) and investigated the key factors that determine the size of an effective region in seizure suppression via electrical stimulation.

Understanding Epileptiform After-Discharges as Rhythmic Oscillatory Transients

Electro-cortical activity in patients with epilepsy may show abnormal rhythmic transients in response to stimulation. Even when using the same stimulation parameters in the same patient, wide variability in the duration of transient response has been reported. These transients have long been considered important for the mapping of the excitability levels in the epileptic brain but their dynamic mechanism is still not well understood. To investigate the occurrence of abnormal transients dynamically, we use a thalamo-cortical neural population model of epileptic spike-wave activity and study the interaction between slow and fast subsystems. In a reduced version of the thalamo-cortical model, slow wave oscillations arise from a fold of cycles (FoC) bifurcation. This marks the onset of a region of bistability between a high amplitude oscillatory rhythm and the background state. In vicinity of the bistability in parameter space, the model has excitable dynamics, showing prolonged rhythmic transients in response to suprathreshold pulse stimulation. We analyse the state space geometry of the bistable and excitable states, and find that the rhythmic transient arises when the impending FoC bifurcation deforms the state space and creates an area of locally reduced attraction to the fixed point. This area essentially allows trajectories to dwell there before escaping to the stable steady state, thus creating rhythmic transients. In the full thalamo-cortical model, we find a similar FoC bifurcation structure. Based on the analysis, we propose an explanation of why stimulation induced epileptiform activity may vary between trials, and predict how the variability could be related to ongoing oscillatory background activity. We compare our dynamic mechanism with other mechanisms (such as a slow parameter change) to generate excitable transients, and we discuss the proposed excitability mechanism in the context of stimulation responses in the epileptic cortex.

Within brain area tractography suggests local modularity using high resolution connectomics

Previous structural brain connectivity studies have mainly focussed on the macroscopic scale of around 1,000 or fewer brain areas (network nodes). However, it has recently been demonstrated that high resolution structural connectomes of around 50,000 nodes can be generated reproducibly. In this study, we infer high resolution brain connectivity matrices using diffusion imaging data from the Human Connectome Project. With such high resolution we are able to analyse networks within brain areas in a single subject. We show that the global network has a scale invariant topological organisation, which means there is a hierarchical organisation of the modular architecture. Specifically, modules within brain areas are spatially localised. We find that long range connections terminate between specific modules, whilst short range connections via highly curved association fibers terminate within modules. We suggest that spatial locations of white matter modules overlap with cytoarchitecturally distinct grey matter areas and may serve as the structural basis for function specialisation within brain areas. Future studies might elucidate how brain diseases change this modular architecture within brain areas.

Controllability Analysis of the Neural Mass Model with Dynamic Parameters

The development of control technology for the brain is of potential significance to the prevention and treatment of neuropsychiatric disorders and the improvement of humans' mental health. A controllability analysis of the brain is necessary to ensure the feasibility of the brain control. In this letter, we investigate the influences of dynamical parameters on the controllability in the neural mass model by using controllability indices as quantitative indicators. The indices are obtained by computing Lie brackets and condition numbers of the system model. We show how controllability changes with important parameters of our dynamical (neuronal) model. Our results suggest that the underlying dynamical parameters have certain ranges with better controllability. We hope it can play potential roles in therapy for brain nervous disorder disease.

Stimulus-induced Epileptic Spike-Wave Discharges in Thalamocortical Model with Disinhibition

Epileptic absence seizure characterized by the typical 2–4 Hz spike-wave discharges (SWD) are known to arise due to the physiologically abnormal interactions within the thalamocortical network. By introducing a second inhibitory neuronal population in the cortical system, here we propose a modified thalamocortical field model to mathematically describe the occurrences and transitions of SWD under the mutual functions between cortex and thalamus, as well as the disinhibitory modulations of SWD mediated by the two different inhibitory interneuronal populations. We first show that stimulation can induce the recurrent seizures of SWD in the modified model. Also, we demonstrate the existence of various types of firing states including the SWD. Moreover, we can identify the bistable parametric regions where the SWD can be both induced and terminated by stimulation perturbations applied in the background resting state. Interestingly, in the absence of stimulation disinhibitory functions between the two different interneuronal populations can also both initiate and abate the SWD, which suggests that the mechanism of disinhibition is comparable to the effect of stimulation in initiating and terminating the epileptic SWD. Hopefully, the obtained results can provide theoretical evidences in exploring dynamical mechanism of epileptic seizures.

Disinhibition-Induced Delayed Onset of Epileptic Spike-Wave Discharges in a Five Variable Model of Cortex and Thalamus

Based on a modified neural field network model composed of cortex and thalamus, we here propose a computational framework to investigate the onset control of absence seizure, which is characterized by the spike-wave discharges. Firstly, we briefly demonstrate the existence of various transition types in Taylor's model by increasing the thalamic input. Furthermore, after the disinhibitory function is reasonably introduced into the Taylor's model, we can observe the occurrence of various transition states of firing patterns with different dominant frequencies as the thalamic input is varied under different disinhibitory effects onto the pyramidal neural population. Interestingly, it is found that the onset of spike-wave discharges can be delayed as the disinhibitory input is considered. More importantly, we explore bifurcation mechanism of firing transitions as some key parameters are changed. And also, we observe other dynamical states, such as simple oscillations and saturated discharges with different spatial scales, which are consistent with previous theoretical or experimental findings.

Predicting Surgery Targets in Temporal Lobe Epilepsy through Structural Connectome Based Simulations

Temporal lobe epilepsy (TLE) is a prevalent neurological disorder resulting in disruptive seizures. In the case of drug resistant epilepsy resective surgery is often considered. This is a procedure hampered by unpredictable success rates, with many patients continuing to have seizures even after surgery. In this study we apply a computational model of epilepsy to patient specific structural connectivity derived from diffusion tensor imaging (DTI) of 22 individuals with left TLE and 39 healthy controls. We validate the model by examining patient-control differences in simulated seizure onset time and network location. We then investigate the potential of the model for surgery prediction by performing in silico surgical resections, removing nodes from patient networks and comparing seizure likelihood post-surgery to pre-surgery simulations. We find that, first, patients tend to transit from non-epileptic to epileptic states more often than controls in the model. Second, regions in the left hemisphere (particularly within temporal and subcortical regions) that are known to be involved in TLE are the most frequent starting points for seizures in patients in the model. In addition, our analysis also implicates regions in the contralateral and frontal locations which may play a role in seizure spreading or surgery resistance. Finally, the model predicts that patient-specific surgery (resection areas chosen on an individual, model-prompted, basis and not following a predefined procedure) may lead to better outcomes than the currently used routine clinical procedure. Taken together this work provides a first step towards patient specific computational modelling of epilepsy surgery in order to inform treatment strategies in individuals.

Temporal lobe epilepsy (TLE) is a disorder characterised by unpredictable seizures, where surgical removal of brain tissue is often the final treatment option. In roughly 30% of cases surgery procedures are unsuccessful at preventing future seizures. This paper shows the application of a computational model which uses patient derived brain connectivity to predict the success rates of surgery in people with TLE. We consider the brains of 22 patients as networks, with brain regions as nodes and the white matter connections between them as edges. The brain network is unique to each subject and produced from brain imaging scans of 22 patients and 39 controls. Seizures are simulated before and after surgery, where surgery in the model is the removal of nodes from the network. The model successfully identifies regions known to be involved in TLE, and its predicted success rates for surgery are close to the results found in reality. The model additionally provides patient specific recommendations for surgical procedures, which in simulations show improved results compared to standard surgery in every case. This is a first step towards designing personalised surgery procedures in order to improve surgery success rates.

Disinhibition-induced transitions between absence and tonic-clonic epileptic seizures

Electrophysiological experiments have long revealed the existence of two-way transitions between absence and tonic-clonic epileptic seizures in the cerebral cortex. Based on a modified spatially-extended Taylor & Baier neural field model, we here propose a computational framework to mathematically describe the transition dynamics between these epileptic seizures. We first demonstrate the existence of various transition types that are induced by disinhibitory functions between two inhibitory variables in an isolated Taylor & Baier model. Moreover, we show that these disinhibition-induced transitions can lead to stable tonic-clonic oscillations as well as periodic spike with slow-wave discharges, which are the hallmark of absence seizures. We also observe fascinating dynamical states, such as periodic 2-spike with slow-wave discharges, tonic death, bursting oscillations, as well as saturated firing. Most importantly, we identify paths that represent physiologically plausible transitions between absence and tonic-clonic seizures in the modified spatially-extended Taylor & Baier model.

Optimal control based seizure abatement using patient derived connectivity

Epilepsy is a neurological disorder in which patients have recurrent seizures. Seizures occur in conjunction with abnormal electrical brain activity which can be recorded by the electroencephalogram (EEG). Often, this abnormal brain activity consists of high amplitude regular spike-wave oscillations as opposed to low amplitude irregular oscillations in the non-seizure state. Active brain stimulation has been proposed as a method to terminate seizures prematurely, however, a general and widely-applicable approach to optimal stimulation protocols is still lacking. In this study we use a computational model of epileptic spike-wave dynamics to evaluate the effectiveness of a pseudospectral method to simulated seizure abatement. We incorporate brain connectivity derived from magnetic resonance imaging of a subject with idiopathic generalized epilepsy. We find that the pseudospectral method can successfully generate time-varying stimuli that abate simulated seizures, even when including heterogeneous patient specific brain connectivity. The strength of the stimulus required varies in different brain areas. Our results suggest that seizure abatement, modeled as an optimal control problem and solved with the pseudospectral method, offers an attractive approach to treatment for in vivo stimulation techniques. Further, if optimal brain stimulation protocols are to be experimentally successful, then the heterogeneity of cortical connectivity should be accounted for in the development of those protocols and thus more spatially localized solutions may be preferable.

A computational study of stimulus driven epileptic seizure abatement

Active brain stimulation to abate epileptic seizures has shown mixed success. In spike-wave (SW) seizures, where the seizure and background state were proposed to coexist, single-pulse stimulations have been suggested to be able to terminate the seizure prematurely. However, several factors can impact success in such a bistable setting. The factors contributing to this have not been fully investigated on a theoretical and mechanistic basis. Our aim is to elucidate mechanisms that influence the success of single-pulse stimulation in noise-induced SW seizures. In this work, we study a neural population model of SW seizures that allows the reconstruction of the basin of attraction of the background activity as a four dimensional geometric object. For the deterministic (noise-free) case, we show how the success of response to stimuli depends on the amplitude and phase of the SW cycle, in addition to the direction of the stimulus in state space. In the case of spontaneous noise-induced seizures, the basin becomes probabilistic introducing some degree of uncertainty to the stimulation outcome while maintaining qualitative features of the noise-free case. Additionally, due to the different time scales involved in SW generation, there is substantial variation between SW cycles, implying that there may not be a fixed set of optimal stimulation parameters for SW seizures. In contrast, the model suggests an adaptive approach to find optimal stimulation parameters patient-specifically, based on real-time estimation of the position in state space. We discuss how the modelling work can be exploited to rationally design a successful stimulation protocol for the abatement of SW seizures using real-time SW detection.

Structural connectivity based whole brain modelling in epilepsy

Epilepsy is a neurological condition characterised by the recurrence of seizures. During seizures multiple brain areas can behave abnormally. Rather than considering each abnormal area in isolation, one can consider them as an interconnected functional 'network'. Recently, there has been a shift in emphasis to consider epilepsy as a disorder involving more widespread functional brain networks than perhaps was previously thought. The basis for these functional networks is proposed to be the static structural brain network established through the connectivity of the white matter. Additionally, it has also been argued that time varying aspects of epilepsy are of crucial importance and as such computational models of these dynamical properties have recently advanced. We describe how dynamic computer models can be combined with static human in vivo connectivity obtained through diffusion weighted magnetic resonance imaging. We predict that in future the use of these two methods in concert will lead to predictions for optimal surgery and brain stimulation sites for epilepsy and other neurological disorders.

Modern concepts of seizure modeling

Seizures are complex phenomena spanning multiple spatial and temporal scales, from ion dynamics to communication between brain regions, from milliseconds (spikes) to days (interseizure intervals). Because of the existence of such multiple scales, the experimental evaluation of the mechanisms underlying the initiation, propagation, and termination of epileptic seizures is a difficult problem. Theoretical models and numerical simulations provide new tools to investigate seizure mechanisms at multiple scales. In this chapter, we review different theoretical approaches and their contributions to our understanding of seizure mechanisms.

Travelling waves in a neural field model with refractoriness

At one level of abstraction neural tissue can be regarded as a medium for turning local synaptic activity into output signals that propagate over large distances via axons to generate further synaptic activity that can cause reverberant activity in networks that possess a mixture of excitatory and inhibitory connections. This output is often taken to be a firing rate, and the mathematical form for the evolution equation of activity depends upon a spatial convolution of this rate with a fixed anatomical connectivity pattern. Such formulations often neglect the metabolic processes that would ultimately limit synaptic activity. Here we reinstate such a process, in the spirit of an original prescription by Wilson and Cowan (Biophys J 12:1–24, 1972), using a term that multiplies the usual spatial convolution with a moving time average of local activity over some refractory time-scale. This modulation can substantially affect network behaviour, and in particular give rise to periodic travelling waves in a purely excitatory network (with exponentially decaying anatomical connectivity), which in the absence of refractoriness would only support travelling fronts. We construct these solutions numerically as stationary periodic solutions in a co-moving frame (of both an equivalent delay differential model as well as the original delay integro-differential model). Continuation methods are used to obtain the dispersion curve for periodic travelling waves (speed as a function of period), and found to be reminiscent of those for spatially extended models of excitable tissue. A kinematic analysis (based on the dispersion curve) predicts the onset of wave instabilities, which are confirmed numerically.

Towards a large-scale model of patient-specific epileptic spike-wave discharges

Clinical electroencephalographic (EEG) recordings of the transition into generalised epileptic seizures show a sudden onset of spike-wave dynamics from a low-amplitude irregular background. In addition, non-trivial and variable spatio-temporal dynamics are widely reported in combined EEG/fMRI studies on the scale of the whole cortex. It is unknown whether these characteristics can be accounted for in a large-scale mathematical model with fixed heterogeneous long-range connectivities. Here, we develop a modelling framework with which to investigate such EEG features. We show that a neural field model composed of a few coupled compartments can serve as a low-dimensional prototype for the transition between irregular background dynamics and spike-wave activity. This prototype then serves as a node in a large-scale network with long-range connectivities derived from human diffusion-tensor imaging data. We examine multivariate properties in 42 clinical EEG seizure recordings from 10 patients diagnosed with typical absence epilepsy and 50 simulated seizures from the large-scale model using 10 DTI connectivity sets from humans. The model can reproduce the clinical feature of stereotypy where seizures are more similar within a patient than between patients, essentially creating a patient-specific fingerprint. We propose the approach as a feasible technique for the investigation of patient-specific large-scale epileptic features in space and time.

The importance of modeling epileptic seizure dynamics as spatio-temporal patterns

The occurrence of seizures is the common feature across the spectrum of epileptic disorders. We describe how the use of mechanistic neural population models leads to novel insight into the dynamic mechanisms underlying two important types of epileptic seizures. We specifically stress the need for a spatio-temporal description of the rhythms to deal with the complexity of the pathophenotype. Adapted to functional and structural patient data, the macroscopic models may allow a patient-specific description of seizures and prediction of treatment outcome.

Phase space approach for modeling of epileptic dynamics

Epileptic electroencephalography recordings can be described in terms of four prototypic wave forms: fast sinusoidal oscillations, large slow waves, fast spiking, and spike waves. On the macroscopic level, these wave forms have been modeled by different mechanistic models which share canonical features. Here we derive a minimal model of excitatory and inhibitory processes with features common to all previous models. We can infer that at least three interacting processes are required to support the prototypic epileptic dynamics. Based on a separation of time scales we analyze the model in terms of interacting manifolds in phase space. This allows qualitative reverse engineering of all epileptic wave forms and transitions between them. We propose this method as a complement to traditional approaches to modeling epileptiform rhythms.

Scaling effects and spatio-temporal multilevel dynamics in epileptic seizures

Epileptic seizures are one of the most well-known dysfunctions of the nervous system. During a seizure, a highly synchronized behavior of neural activity is observed that can cause symptoms ranging from mild sensual malfunctions to the complete loss of body control. In this paper, we aim to contribute towards a better understanding of the dynamical systems phenomena that cause seizures. Based on data analysis and modelling, seizure dynamics can be identified to possess multiple spatial scales and on each spatial scale also multiple time scales. At each scale, we reach several novel insights. On the smallest spatial scale we consider single model neurons and investigate early-warning signs of spiking. This introduces the theory of critical transitions to excitable systems. For clusters of neurons (or neuronal regions) we use patient data and find oscillatory behavior and new scaling laws near the seizure onset. These scalings lead to substantiate the conjecture obtained from mean-field models that a Hopf bifurcation could be involved near seizure onset. On the largest spatial scale we introduce a measure based on phase-locking intervals and wavelets into seizure modelling. It is used to resolve synchronization between different regions in the brain and identifies time-shifted scaling laws at different wavelet scales. We also compare our wavelet-based multiscale approach with maximum linear cross-correlation and mean-phase coherence measures.