Optimal control based seizure abatement using patient derived connectivity

Virtual brain twins: from basic neuroscience to clinical use

Virtual brain twins are personalized, generative and adaptive brain models based on data from an individual's brain for scientific and clinical use. After a description of the key elements of virtual brain twins, we present the standard model for personalized whole-brain network models. The personalization is accomplished using a subject's brain imaging data by three means: (1) assemble cortical and subcortical areas in the subject-specific brain space; (2) directly map connectivity into the brain models, which can be generalized to other parameters; and (3) estimate relevant parameters through model inversion, typically using probabilistic machine learning. We present the use of personalized whole-brain network models in healthy ageing and five clinical diseases: epilepsy, Alzheimer's disease, multiple sclerosis, Parkinson's disease and psychiatric disorders. Specifically, we introduce spatial masks for relevant parameters and demonstrate their use based on the physiological and pathophysiological hypotheses. Finally, we pinpoint the key challenges and future directions.

Spike-and-wave discharges of absence seizures in a sleep waves-constrained corticothalamic model

AIMS

Recurrent network activity in corticothalamic circuits generates physiological and pathological EEG waves. Many computer models have simulated spike-and-wave discharges (SWDs), the EEG hallmark of absence seizures (ASs). However, these models either provided detailed simulated activity only in a selected territory (i.e., cortical or thalamic) or did not test whether their corticothalamic networks could reproduce the physiological activities that are generated by these circuits.

METHODS

Using a biophysical large-scale corticothalamic model that reproduces the full extent of EEG sleep waves, including sleep spindles, delta, and slow (<1 Hz) waves, here we investigated how single abnormalities in voltage- or transmitter-gated channels in the neocortex or thalamus led to SWDs.

RESULTS

We found that a selective increase in the tonic γ-aminobutyric acid type A receptor (GABA-A) inhibition of first-order thalamocortical (TC) neurons or a selective decrease in cortical phasic GABA-A inhibition is sufficient to generate ~4 Hz SWDs (as in humans) that invariably start in neocortical territories. Decreasing the leak conductance of higher-order TC neurons leads to ~7 Hz SWDs (as in rodent models) while maintaining sleep spindles at 7-14 Hz.

CONCLUSION

By challenging key features of current mechanistic views, this simulated ictal corticothalamic activity provides novel understanding of ASs and makes key testable predictions.

Personalized dynamic network models of the human brain as a future tool for planning and optimizing epilepsy therapy

Epilepsy is a common neurological disorder, with one third of patients not responding to currently available antiepileptic drugs. The proportion of pharmacoresistant epilepsies has remained unchanged for many decades. To cure epilepsy and control seizures requires a paradigm shift in the development of new approaches to epilepsy diagnosis and treatment. Contemporary medicine has benefited from the exponential growth of computational modeling, and the application of network dynamics theory to understanding and treating human brain disorders. In epilepsy, the introduction of these approaches has led to personalized epileptic network modeling that can explore the patient's seizure genesis and predict the functional impact of resection on its individual network's propensity to seize. The application of the dynamic systems approach to neurostimulation therapy of epilepsy allows designing stimulation strategies that consider the patient's seizure dynamics and long-term fluctuations in the stability of their epileptic networks. In this article, we review, in a nontechnical fashion suitable for a broad neuroscientific audience, recent progress in personalized dynamic brain network modeling that is shaping the future approach to the diagnosis and treatment of epilepsy.

Closed-loop optogenetic control of the dynamics of neural activity in non-human primates

Electrical neurostimulation is effective in the treatment of neurological disorders, but associated recording artefacts generally limit its applications to open-loop stimuli. Real-time and continuous closed-loop control of brain activity can, however, be achieved by pairing concurrent electrical recordings and optogenetics. Here we show that closed-loop optogenetic stimulation with excitatory opsins enables the precise manipulation of neural dynamics in brain slices from transgenic mice and in anaesthetized non-human primates. The approach generates oscillations in quiescent tissue, enhances or suppresses endogenous patterns in active tissue and modulates seizure-like bursts elicited by the convulsant 4-aminopyridine. A nonlinear model of the phase-dependent effects of optical stimulation reproduced the modulation of cycles of local-field potentials associated with seizure oscillations, as evidenced by the systematic changes in the variability and entropy of the phase-space trajectories of seizures, which correlated with changes in their duration and intensity. We also show that closed-loop optogenetic neurostimulation could be delivered using intracortical optrodes incorporating light-emitting diodes. Closed-loop optogenetic approaches may be translatable to therapeutic applications in humans.

Suppression of seizure in childhood absence epilepsy using robust control of deep brain stimulation: a simulation study

Deep brain stimulation (DBS) is a promising technique to relieve the symptoms in patients with intractable seizures. Although the DBS therapy for seizure suppression dates back more than 40 years, determining stimulation parameters is a significant challenge to the success of this technique. One solution to this challenge with application in a real DBS system is to design a closed-loop control system to regulate the stimulation intensity using computational models of epilepsy automatically. The main goal of the current study is to develop a robust control technique based on adaptive fuzzy terminal sliding mode control (AFTSMC) for eliminating the oscillatory spiking behavior in childhood absence epilepsy (CAE) dynamical model consisting of cortical, thalamic relay, and reticular nuclei neurons. To this end, the membrane voltage dynamics of the three coupled neurons are considered as a three-input three-output nonlinear state delay system. A fuzzy logic system is developed to estimate the unknown nonlinear dynamics of the current and delayed states of the model embedded in the control input. Chattering-free control input (continuous DBS pulses) without any singularity problem is the superiority of the proposed control method. To guarantee the bounded stability of the closed-loop system in a finite time, the upper bounds of the external disturbance and minimum estimation errors are updated online with adaptive laws without any offline tuning phase. Simulation results are provided to show the robustness of AFTSMC in the presence of uncertainty and external disturbances.

Towards network-guided neuromodulation for epilepsy

Epilepsy is well-recognized as a disorder of brain networks. There is a growing body of research to identify critical nodes within dynamic epileptic networks with the aim to target therapies that halt the onset and propagation of seizures. In parallel, intracranial neuromodulation, including deep brain stimulation and responsive neurostimulation, are well-established and expanding as therapies to reduce seizures in adults with focal-onset epilepsy; and there is emerging evidence for their efficacy in children and generalized-onset seizure disorders. The convergence of these advancing fields is driving an era of 'network-guided neuromodulation' for epilepsy. In this review, we distil the current literature on network mechanisms underlying neurostimulation for epilepsy. We discuss the modulation of key 'propagation points' in the epileptogenic network, focusing primarily on thalamic nuclei targeted in current clinical practice. These include (i) the anterior nucleus of thalamus, now a clinically approved and targeted site for open loop stimulation, and increasingly targeted for responsive neurostimulation; and (ii) the centromedian nucleus of the thalamus, a target for both deep brain stimulation and responsive neurostimulation in generalized-onset epilepsies. We discuss briefly the networks associated with other emerging neuromodulation targets, such as the pulvinar of the thalamus, piriform cortex, septal area, subthalamic nucleus, cerebellum and others. We report synergistic findings garnered from multiple modalities of investigation that have revealed structural and functional networks associated with these propagation points - including scalp and invasive EEG, and diffusion and functional MRI. We also report on intracranial recordings from implanted devices which provide us data on the dynamic networks we are aiming to modulate. Finally, we review the continuing evolution of network-guided neuromodulation for epilepsy to accelerate progress towards two translational goals: (i) to use pre-surgical network analyses to determine patient candidacy for neurostimulation for epilepsy by providing network biomarkers that predict efficacy; and (ii) to deliver precise, personalized and effective antiepileptic stimulation to prevent and arrest seizure propagation through mapping and modulation of each patients' individual epileptogenic networks.

Prevalence and scalable control of localized networks

The ability to control network dynamics is essential for ensuring desirable functionality of many technological, biological, and social systems. Such systems often consist of a large number of network elements, and controlling large-scale networks remains challenging because the computation and communication requirements increase prohibitively fast with network size. Here, we introduce a notion of network locality that can be exploited to make the control of networks scalable, even when the dynamics are nonlinear. We show that network locality is captured by an information metric and is almost universally observed across real and model networks. In localized networks, the optimal control actions and system responses are both shown to be necessarily concentrated in small neighborhoods induced by the information metric. This allows us to develop localized algorithms for determining network controllability and optimizing the placement of driver nodes. This also allows us to develop a localized algorithm for designing local feedback controllers that approach the performance of the corresponding best global controllers, while incurring a computational cost orders-of-magnitude lower. We validate the locality, performance, and efficiency of the algorithms in Kuramoto oscillator networks, as well as three large empirical networks: synchronization dynamics in the Eastern US power grid, epidemic spreading mediated by the global air-transportation network, and Alzheimer's disease dynamics in a human brain network. Taken together, our results establish that large networks can be controlled with computation and communication costs comparable to those for small networks.

Role of coupling distances in a coupled thalamocortical network for regulation of epilepsy

The recent theoretical modeling of coupled cortical thalamic network is an important advance toward the spatiotemporal dynamics of the brain. However, the diversity of coupling distances is ignored, and the better choice of deep brain stimulation (DBS) parameters to control epilepsy is still a challenge so far. A modeling object of this paper is to establish a coupled cortical thalamic model with uncertain coupling distances including nine combinations. Based on the pathways formed by pyramidal neuronal population (PY), thalamic reticular nucleus (RE) and thalamic relay nucleus (TC), we simulate the spike-wave discharges (SWD) at 2-4Hz which are the main manifestations of absence episodes. It is demonstrated that combination (1/3, 1/9) between the left and right ventricles is the optimal coupling distance of the proposed model by analyzing the percentage of SWD. A stimulating object of this paper is to find an optimum parameter range of DBS. One of the important results is that the number of SWD is inversely proportional to the amplitude, another one is that the number of SWD shows a U-shaped trend with the change of frequency. The present study has laidtheoryfoundationforthebrainplasticity, which will provide an important theoretical basis and direction for the treatment of absence epilepsy in the future. In brief, hopefully our simulation results will provide some help to patients.

Drug-resistant focal epilepsy in children is associated with increased modal controllability of the whole brain and epileptogenic regions

Network control theory provides a framework by which neurophysiological dynamics of the brain can be modelled as a function of the structural connectome constructed from diffusion MRI. Average controllability describes the ability of a region to drive the brain to easy-to-reach neurophysiological states whilst modal controllability describes the ability of a region to drive the brain to difficult-to-reach states. In this study, we identify increases in mean average and modal controllability in children with drug-resistant epilepsy compared to healthy controls. Using simulations, we purport that these changes may be a result of increased thalamocortical connectivity. At the node level, we demonstrate decreased modal controllability in the thalamus and posterior cingulate regions. In those undergoing resective surgery, we also demonstrate increased modal controllability of the resected parcels, a finding specific to patients who were rendered seizure free following surgery. Changes in controllability are a manifestation of brain network dysfunction in epilepsy and may be a useful construct to understand the pathophysiology of this archetypical network disease. Understanding the mechanisms underlying these controllability changes may also facilitate the design of network-focussed interventions that seek to normalise network structure and function.

Control theory illustrates the energy efficiency in the dynamic reconfiguration of functional connectivity

The brain’s functional connectivity fluctuates over time instead of remaining steady in a stationary mode even during the resting state. This fluctuation establishes the dynamical functional connectivity that transitions in a non-random order between multiple modes. Yet it remains unexplored how the transition facilitates the entire brain network as a dynamical system and what utility this mechanism for dynamic reconfiguration can bring over the widely used graph theoretical measurements. To address these questions, we propose to conduct an energetic analysis of functional brain networks using resting-state fMRI and behavioral measurements from the Human Connectome Project. Through comparing the state transition energy under distinct adjacent matrices, we justify that dynamic functional connectivity leads to 60% less energy cost to support the resting state dynamics than static connectivity when driving the transition through default mode network. Moreover, we demonstrate that combining graph theoretical measurements and our energy-based control measurements as the feature vector can provide complementary prediction power for the behavioral scores (Combination vs. Control: t = 9.41, p = 1.64e−13; Combination vs. Graph: t = 4.92, p = 3.81e−6). Our approach integrates statistical inference and dynamical system inspection towards understanding brain networks.

A framework that allows for the statistical investigation of the dynamic aspect of functional connectivity derived from resting-state fMRI is developed that is shown to complementarily predict individual differences in behavioral measurements compared to existing approaches.

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.

The Critical Modulatory Role of Spiny Stellate Cells in Seizure Onset Based on Dynamic Analysis of a Neural Mass Model

Growing evidence suggests that excitatory neurons in the brain play a significant role in seizure generation. Nonetheless, spiny stellate cells are cortical excitatory non-pyramidal neurons in the brain, whose basic role in seizure occurrence is not well understood. In the present research, we study the critical role of spiny stellate cells or the excitatory interneurons (EI), for the first time, in epileptic seizure generation using an extended neural mass model inspired by a thalamocortical model originally introduced by another research group. Applying bifurcation analysis on this modified model, we investigated the rich dynamics corresponding to the epileptic seizure onset and transition between interictal and ictal states caused by EI connectivity to other cell types. Our results indicate that the transition between interictal and ictal states (preictal signal) corresponds to a supercritical Hopf bifurcation, and thus, the extended model suggests that before seizure onset, the amplitude and frequency of neural activities gradually increase. Moreover, we showed that (1) the altered function of GABAergic and glutamatergic receptors of EI can cause seizure, and (2) the pathway between the thalamic relay nucleus and EI facilitates the transition from interictal to ictal activity by decreasing the preictal period. Thereafter, we considered both sensory and cortical periodic inputs to study model responses to various harmonic stimulations. Bifurcation analysis of the model, in this case, suggests that the initial state of the model might be the main cause for the transition between interictal and ictal states as the stimulus frequency changes. The extended thalamocortical model shows also that the amplitude jump phenomenon and non-linear resonance behavior result from the preictal state of the modified model. These results can be considered as a step forward to a deeper understanding of the mechanisms underlying the transition from normal activities to epileptic activities.

Control of epileptic seizures by electrical stimulation: a model-based study

High frequency electrical stimulation of brain is commonly used in research experiments and clinical trials as a modern tool for control of epileptic seizures. However, the mechanistic basis by which periodic external stimuli alter the brain state is not well understood. This study provides a computational insight into the mechanism of seizure suppression by high frequency stimulation (HFS). In particular, a modified version of the Jansen-Rit neural mass model is employed, in which EEG signals can be considered as the input. The proposed model reproduces seizure-like activity in the output during the ictal period of the input signal. By applying a control signal to the model, a wide range of stimulation amplitudes and frequencies are systematically explored. Simulation results reveal that HFS can effectively suppress the seizure-like activity. Our results suggest that HFS has the ability of shifting the operating state of neural populations away from a critical condition. Furthermore, a closed-loop control strategy is proposed in this paper. The main objective has been to considerably reduce the control effort needed for blocking abnormal activity of the brain. Such an energy reduction could be of practical importance, to reduce possible side effects and increase battery life for implanted neurostimulators.

Falling off a limit cycle using phase-agnostic stimuli: Applications to clinical oscillopathies

For over a century, physiological studies have shown that precisely timed pulses can switch off a biological oscillator. This empiric finding has shaped our mechanistic understanding of how perturbations start, stop, and reset biological oscillators and has led to treatments that suppress pathological oscillations using electrical pulses given within specified therapeutic phase windows. Here, we present evidence, using numerical simulations of models of epileptic seizures and reentrant tachycardia, that the phase window can be opened to the entire cycle using novel complex stimulus waveforms. Our results reveal that the trajectories are displaced by such phase-agnostic stimuli off the oscillator's limit cycle and corralled into a region where oscillation is suppressed, irrespective of the phase at which the stimulus was applied. Our findings suggest the need for broadening theoretical understanding of how complex perturbing waveforms interact with biological oscillators to access their arrhythmic states. In clinical practice, oscillopathies may be treated more effectively with non-traditional stimulus waveforms that obviate the need for phase specificity.

Dynamic Transitions of Epilepsy Waveforms Induced by Astrocyte Dysfunction and Electrical Stimulation

Experimental studies have shown that astrocytes participate in epilepsy through inducing the release of glutamate. Meanwhile, considering the disinhibition circuit among inhibitory neuronal populations with different time scales and the feedforward inhibition connection from thalamic relay nucleus to cortical inhibitory neuronal population, here, we propose a modified thalamocortical field model to systematically investigate the mechanism of epilepsy. Firstly, our results show that rich firing activities can be induced by astrocyte dysfunction, including high or low saturated state, high- or low-frequency clonic, spike-wave discharge (SWD), and tonic. More importantly, with the enhancement of feedforward inhibition connection, SWD and tonic oscillations will disappear. In other words, all these pathological waveforms can be suppressed or eliminated. Then, we explore the control effects after different external stimulations applying to thalamic neuronal population. We find that single-pulse stimulation can not only suppress but also induce pathological firing patterns, such as SWD, tonic, and clonic oscillations. And we further verify that deep brain stimulation can control absence epilepsy by regulating the amplitude and pulse width of stimulation. In addition, based on our modified model, 3 : 2 coordinated reset stimulation strategies with different intensities are compared and a more effective and safer stimulation mode is proposed. Our conclusions are expected to give more theoretical insights into the treatment of epilepsy.

Control analysis of electrical stimulation for epilepsy waveforms in a thalamocortical network

Physiological experiments and computational models both show that the thalamic reticular nucleus (RE) participates in inducing various firing patterns of cortex. Absence seizure, featured by 2-4 Hz spike-wave discharges (SWD) oscillation, is a high incidence of disease in children. Lots of electrophysiological experiments have verified the correlation between absence seizures and RE, however, the dynamical mechanisms are not well understood. Based on previous Taylor model, we firstly study the effects of external input and self-inhibition of RE on epilepsy transition. We show that increasing external input and self-inhibition of RE can lead the system from epileptic state to normal state, and vice versa. Next, we explore two stimulus strategies added in RE and various transition behaviors can be induced, such as high saturated state to clonic. Meanwhile, as the intensity of stimulation increasing, they can not only suppress the SWD, but also produce tonic-clonic oscillation. Finally, the control of DBS on single neuron cluster and two neuron clusters are compared and we find stimulating RE and TC simultaneously is a superior mode to stimulate anyone of RE or TC. It is hoped that the results we obtained will have an enlightenment on clinical treatment.

Network structure of cascading neural systems predicts stimulus propagation and recovery

OBJECTIVE

Many neural systems display spontaneous, spatiotemporal patterns of neural activity that are crucial for information processing. While these cascading patterns presumably arise from the underlying network of synaptic connections between neurons, the precise contribution of the network's local and global connectivity to these patterns and information processing remains largely unknown.

APPROACH

Here, we demonstrate how network structure supports information processing through network dynamics in empirical and simulated spiking neurons using mathematical tools from linear systems theory, network control theory, and information theory.

MAIN RESULTS

In particular, we show that activity, and the information that it contains, travels through cycles in real and simulated networks.

SIGNIFICANCE

Broadly, our results demonstrate how cascading neural networks could contribute to cognitive faculties that require lasting activation of neuronal patterns, such as working memory or attention.

White Matter Network Architecture Guides Direct Electrical Stimulation through Optimal State Transitions

Optimizing direct electrical stimulation for the treatment of neurological disease remains difficult due to an incomplete understanding of its physical propagation through brain tissue. Here, we use network control theory to predict how stimulation spreads through white matter to influence spatially distributed dynamics. We test the theory’s predictions using a unique dataset comprising diffusion weighted imaging and electrocorticography in epilepsy patients undergoing grid stimulation. We find statistically significant shared variance between the predicted activity state transitions and the observed activity state transitions. We then use an optimal control framework to posit testable hypotheses regarding which brain states and structural properties will efficiently improve memory encoding when stimulated. Our work quantifies the role that white matter architecture plays in guiding the dynamics of direct electrical stimulation and offers empirical support for the utility of network control theory in explaining the brain’s response to stimulation.

Stiso et al. report evidence that network control theory can explain the propagation of electrical stimulation through the human brain and quantify how white matter connectivity is crucial for driving spatially distributed changes in activity. Furthermore, they use network control theory to predict stimulation outcome in specific cognitive contexts.

Control of brain network dynamics across diverse scales of space and time

The human brain is composed of distinct regions that are each associated with particular functions and distinct propensities for the control of neural dynamics. However, the relation between these functions and control profiles is poorly understood, as is the variation in this relation across diverse scales of space and time. Here we probe the relation between control and dynamics in brain networks constructed from diffusion tensor imaging data in a large community sample of young adults. Specifically, we probe the control properties of each brain region and investigate their relationship with dynamics across various spatial scales using the Laplacian eigenspectrum. In addition, through analysis of regional modal controllability and partitioning of modes, we determine whether the associated dynamics are fast or slow, as well as whether they are alternating or monotone. We find that brain regions that facilitate the control of energetically easy transitions are associated with activity on short length scales and slow timescales. Conversely, brain regions that facilitate control of difficult transitions are associated with activity on long length scales and fast timescales. Built on linear dynamical models, our results offer parsimonious explanations for the activity propagation and network control profiles supported by regions of differing neuroanatomical structure.

Virtual resection predicts surgical outcome for drug-resistant epilepsy

Patients with drug-resistant epilepsy often require surgery to become seizure-free. While laser ablation and implantable stimulation devices have lowered the morbidity of these procedures, seizure-free rates have not dramatically improved, particularly for patients without focal lesions. This is in part because it is often unclear where to intervene in these cases. To address this clinical need, several research groups have published methods to map epileptic networks but applying them to improve patient care remains a challenge. In this study we advance clinical translation of these methods by: (i) presenting and sharing a robust pipeline to rigorously quantify the boundaries of the resection zone and determining which intracranial EEG electrodes lie within it; (ii) validating a brain network model on a retrospective cohort of 28 patients with drug-resistant epilepsy implanted with intracranial electrodes prior to surgical resection; and (iii) sharing all neuroimaging, annotated electrophysiology, and clinical metadata to facilitate future collaboration. Our network methods accurately forecast whether patients are likely to benefit from surgical intervention based on synchronizability of intracranial EEG (area under the receiver operating characteristic curve of 0.89) and provide novel information that traditional electrographic features do not. We further report that removing synchronizing brain regions is associated with improved clinical outcome, and postulate that sparing desynchronizing regions may further be beneficial. Our findings suggest that data-driven network-based methods can identify patients likely to benefit from resective or ablative therapy, and perhaps prevent invasive interventions in those unlikely to do so.

A practical guide to methodological considerations in the controllability of structural brain networks

OBJECTIVE

Predicting how the brain can be driven to specific states by means of internal or external control requires a fundamental understanding of the relationship between neural connectivity and activity. Network control theory is a powerful tool from the physical and engineering sciences that can provide insights regarding that relationship; it formalizes the study of how the dynamics of a complex system can arise from its underlying structure of interconnected units.

APPROACH

Given the recent use of network control theory in neuroscience, it is now timely to offer a practical guide to methodological considerations in the controllability of structural brain networks. Here we provide a systematic overview of the framework, examine the impact of modeling choices on frequently studied control metrics, and suggest potentially useful theoretical extensions. We ground our discussions, numerical demonstrations, and theoretical advances in a dataset of high-resolution diffusion imaging with 730 diffusion directions acquired over approximately 1 h of scanning from ten healthy young adults.

MAIN RESULTS

Following a didactic introduction of the theory, we probe how a selection of modeling choices affects four common statistics: average controllability, modal controllability, minimum control energy, and optimal control energy. Next, we extend the current state-of-the-art in two ways: first, by developing an alternative measure of structural connectivity that accounts for radial propagation of activity through abutting tissue, and second, by defining a complementary metric quantifying the complexity of the energy landscape of a system. We close with specific modeling recommendations and a discussion of methodological constraints.

SIGNIFICANCE

Our hope is that this accessible account will inspire the neuroimaging community to more fully exploit the potential of network control theory in tackling pressing questions in cognitive, developmental, and clinical neuroscience.

Data-Driven Approaches to Neuroimaging Analysis to Enhance Psychiatric Diagnosis and Therapy

Combining advanced neuroimaging with novel computational methods in network science and machine learning has led to increasingly meaningful descriptions of structure and function in both the normal and the abnormal brain, thereby contributing significantly to our understanding of psychiatric disorders as circuit dysfunctions. Despite its marked potential for psychiatric care, this approach has not yet extended beyond the research setting to any clinically useful applications. Here we review current developments in the study of neuroimaging data using network models and machine learning methods, with a focus on their promise in offering a framework for clinical translation. We discuss 3 potential contributions of these methods to psychiatric care: 1) a better understanding of psychopathology beyond current diagnostic boundaries; 2) individualized prediction of treatment response and prognosis; and 3) formal theories to guide the development of novel interventions. Finally, we highlight current obstacles and sketch a forward-looking perspective of how the application of machine learning and network modeling methods should proceed to accelerate their potential transformation of clinically useful tools.

Model-based robust suppression of epileptic seizures without sensory measurements

Uncontrolled seizures may lead to irreversible damages in the brain and various limitations in the patient's life. There exist experimental studies to stabilize the patient seizures. However, the experimental setups have many sensory devices to measure the dynamics of the brain cortex. These equipments prevent to produce small portable stabilizers for patients in everyday life. Recently, a comprehensive cortex model is introduced to apply model-based observers and controllers. However, this cortex model can be uncertain and have time-varying parameters. Therefore, in this paper, a robust Takagi-Sugeno (TS) controller and observer are designed to suppress the epileptic seizures without sensory measurements. The unavailable sensory measurements are provided by the designed nonlinear observer. The exponential convergence of the observer and controller is satisfied by the feedback parameter design using linear matrix inequalities. In addition, TS fuzzy observer-controller design has been compared with the conventional PID method in terms of control performance and design problem. The numerical computations show that the epileptic seizures are more effectively suppressed by the TS fuzzy observer-based controller under uncertain membrane potential dynamics.

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.

Functional control of electrophysiological network architecture using direct neurostimulation in humans

Chronically implantable neurostimulation devices are becoming a clinically viable option for treating patients with neurological disease and psychiatric disorders. Neurostimulation offers the ability to probe and manipulate distributed networks of interacting brain areas in dysfunctional circuits. Here, we use tools from network control theory to examine the dynamic reconfiguration of functionally interacting neuronal ensembles during targeted neurostimulation of cortical and subcortical brain structures. By integrating multimodal intracranial recordings and diffusion-weighted imaging from patients with drug-resistant epilepsy, we test hypothesized structural and functional rules that predict altered patterns of synchronized local field potentials. We demonstrate the ability to predictably reconfigure functional interactions depending on stimulation strength and location. Stimulation of areas with structurally weak connections largely modulates the functional hubness of downstream areas and concurrently propels the brain towards more difficult-to-reach dynamical states. By using focal perturbations to bridge large-scale structure, function, and markers of behavior, our findings suggest that stimulation may be tuned to influence different scales of network interactions driving cognition.

Brain stimulation devices capable of perturbing the physiological state of neural systems are rapidly gaining popularity for their potential to treat neurological and psychiatric disease. A root problem is that underlying dysfunction spans a large-scale network of brain regions, requiring the ability to control the complex interactions between multiple brain areas. Here, we use tools from network control theory to examine the dynamic reconfiguration of functionally interacting neuronal ensembles during targeted neurostimulation of cortical and subcortical brain structures. We demonstrate the ability to predictably reconfigure patterns of interactions between functional brain areas by modulating the strength and location of stimulation. Our findings have high significance for designing stimulation protocols capable of modulating distributed neural circuits in the human brain.

Robust closed-loop control of spike-and-wave discharges in a thalamocortical computational model of absence epilepsy

In this paper, we investigate the abatement of spike-and-wave discharges in a thalamocortical model using a closed-loop brain stimulation method. We first explore the complex states and various transitions in the thalamocortical computational model of absence epilepsy by using bifurcation analysis. We demonstrate that the Hopf and double cycle bifurcations are the key dynamical mechanisms of the experimental observed bidirectional communications during absence seizures through top-down cortical excitation and thalamic feedforward inhibition. Then, we formulate the abatement of epileptic seizures to a closed-loop tracking control problem. Finally, we propose a neural network based sliding mode feedback control system to drive the dynamics of pathological cortical area to track the desired normal background activities. The control system is robust to uncertainties and disturbances, and its stability is guaranteed by Lyapunov stability theorem. Our results suggest that the seizure abatement can be modeled as a tracking control problem and solved by a robust closed-loop control method, which provides a promising brain stimulation strategy.

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.

Network neuroscience for optimizing brain-computer interfaces

Human-machine interactions are being increasingly explored to create alternative ways of communication and to improve our daily life. Based on a classification of the user's intention from the user's underlying neural activity, brain-computer interfaces (BCIs) allow direct interactions with the external environment while bypassing the traditional effector of the musculoskeletal system. Despite the enormous potential of BCIs, there are still a number of challenges that limit their societal impact, ranging from the correct decoding of a human's thoughts, to the application of effective learning strategies. Despite several important engineering advances, the basic neuroscience behind these challenges remains poorly explored. Indeed, BCIs involve complex dynamic changes related to neural plasticity at a diverse range of spatiotemporal scales. One promising antidote to this complexity lies in network science, which provides a natural language in which to model the organizational principles of brain architecture and function as manifest in its interconnectivity. Here, we briefly review the main limitations currently affecting BCIs, and we offer our perspective on how they can be addressed by means of network theoretic approaches. We posit that the emerging field of network neuroscience will prove to be an effective tool to unlock human-machine interactions.

Seizure pathways: A model-based investigation

We present the results of a model inversion algorithm for electrocorticography (ECoG) data recorded during epileptic seizures. The states and parameters of neural mass models were tracked during a total of over 3000 seizures from twelve patients with focal epilepsy. These models provide an estimate of the effective connectivity within intracortical circuits over the time course of seizures. Observing the dynamics of effective connectivity provides insight into mechanisms of seizures. Estimation of patients seizure dynamics revealed: 1) a highly stereotyped pattern of evolution for each patient, 2) distinct sub-groups of onset mechanisms amongst patients, and 3) different offset mechanisms for long and short seizures. Stereotypical dynamics suggest that, once initiated, seizures follow a deterministic path through the parameter space of a neural model. Furthermore, distinct sub-populations of patients were identified based on characteristic motifs in the dynamics at seizure onset. There were also distinct patterns between long and short duration seizures that were related to seizure offset. Understanding how these different patterns of seizure evolution arise may provide new insights into brain function and guide treatment for epilepsy, since specific therapies may have preferential effects on the various parameters that could potentially be individualized. Methods that unite computational models with data provide a powerful means to generate testable hypotheses for further experimental research. This work provides a demonstration that the hidden connectivity parameters of a neural mass model can be dynamically inferred from data. Our results underscore the power of theoretical models to inform epilepsy management. It is our hope that this work guides further efforts to apply computational models to clinical data.

On the nature and use of models in network neuroscience

Network theory provides an intuitively appealing framework for studying relationships among interconnected brain mechanisms and their relevance to behaviour. As the space of its applications grows, so does the diversity of meanings of the term network model. This diversity can cause confusion, complicate efforts to assess model validity and efficacy, and hamper interdisciplinary collaboration. In this Review, we examine the field of network neuroscience, focusing on organizing principles that can help overcome these challenges. First, we describe the fundamental goals in constructing network models. Second, we review the most common forms of network models, which can be described parsimoniously along the following three primary dimensions: from data representations to first-principles theory; from biophysical realism to functional phenomenology; and from elementary descriptions to coarse-grained approximations. Third, we draw on biology, philosophy and other disciplines to establish validation principles for these models. We close with a discussion of opportunities to bridge model types and point to exciting frontiers for future pursuits.

Design of optimal nonlinear network controllers for Alzheimer's disease

Brain stimulation can modulate the activity of neural circuits impaired by Alzheimer’s disease (AD), having promising clinical benefit. However, all individuals with the same condition currently receive identical brain stimulation, with limited theoretical basis for this generic approach. In this study, we introduce a control theory framework for obtaining exogenous signals that revert pathological electroencephalographic activity in AD at a minimal energetic cost, while reflecting patients’ biological variability. We used anatomical networks obtained from diffusion magnetic resonance images acquired by the Alzheimer’s Disease Neuroimaging Initiative (ADNI) as mediators for the interaction between Duffing oscillators. The nonlinear nature of the brain dynamics is preserved, given that we extend the so-called state-dependent Riccati equation control to reflect the stimulation objective in the high-dimensional neural system. By considering nonlinearities in our model, we identified regions for which control inputs fail to correct abnormal activity. There are changes to the way stimulated regions are ranked in terms of the energetic cost of controlling the entire network, from a linear to a nonlinear approach. We also found that limbic system and basal ganglia structures constitute the top target locations for stimulation in AD. Patients with highly integrated anatomical networks–namely, networks having low average shortest path length, high global efficiency–are the most suitable candidates for the propagation of stimuli and consequent success on the control task. Other diseases associated with alterations in brain dynamics and the self-control mechanisms of the brain can be addressed through our framework.

This work aims to close the knowledge gap between theory and experiment in brain stimulation. Previous modeling approaches for stimulation have overlooked the nonlinear dynamical nature of the brain and failed to shed light on efficient mechanisms for the exogenous control of the brain. Amid the current efforts for developing personalized medicine, we introduce a framework for producing tailored stimulation signals, based on individual neuroimaging data and innovative modeling. This is the first time, to our knowledge, that brain stimulation for the most common cause of dementia, Alzheimer’s disease, is theoretically addressed. Our approach leads to the identification of potential target regions and subjects to successfully respond to brain stimulation therapies and yields various disease-reverting signals. Although focused on Alzheimer’s in this study, our methodology could be applied to other clinical conditions characterized by abnormalities in brain dynamics, like epilepsy and Parkinson’s, the treatment of which can benefit from the use of optimal control strategies.

Designing Patient-Specific Optimal Neurostimulation Patterns for Seizure Suppression

Neurostimulation is a promising therapy for abating epileptic seizures. However, it is extremely difficult to identify optimal stimulation patterns experimentally. In this study, human recordings are used to develop a functional 24 neuron network statistical model of hippocampal connectivity and dynamics. Spontaneous seizure-like activity is induced in silico in this reconstructed neuronal network. The network is then used as a testbed to design and validate a wide range of neurostimulation patterns. Commonly used periodic trains were not able to permanently abate seizures at any frequency. A simulated annealing global optimization algorithm was then used to identify an optimal stimulation pattern, which successfully abated 92% of seizures. Finally, in a fully responsive, or closed-loop, neurostimulation paradigm, the optimal stimulation successfully prevented the network from entering the seizure state. We propose that the framework presented here for algorithmically identifying patient-specific neurostimulation patterns can greatly increase the efficacy of neurostimulation devices for seizures.

From Maps to Multi-dimensional Network Mechanisms of Mental Disorders

The development of advanced neuroimaging techniques and their deployment in large cohorts has enabled an assessment of functional and structural brain network architecture at an unprecedented level of detail. Across many temporal and spatial scales, network neuroscience has emerged as a central focus of intellectual efforts, seeking meaningful descriptions of brain networks and explanatory sets of network features that underlie circuit function in health and dysfunction in disease. However, the tools of network science commonly deployed provide insight into brain function at a fundamentally descriptive level, often failing to identify (patho-)physiological mechanisms that link system-level phenomena to the multiple hierarchies of brain function. Here we describe recently developed techniques stemming from advances in complex systems and network science that have the potential to overcome this limitation, thereby contributing mechanistic insights into neuroanatomy, functional dynamics, and pathology. Finally, we build on the Research Domain Criteria framework, highlighting the notion that mental illnesses can be conceptualized as dysfunctions of neural circuitry present across conventional diagnostic boundaries, to sketch how network-based methods can be combined with pharmacological, intermediate phenotype, genetic, and magnetic stimulation studies to probe mechanisms of psychopathology.

Controlling mechanism of absence seizures by deep brain stimulus applied on subthalamic nucleus

Based on a classical model of the basal ganglia thalamocortical network, in this paper, we employed a type of the deep brain stimulus voltage on the subthalamic nucleus to study the control mechanism of absence epilepsy seizures. We found that the seizure can be well controlled by turning the period and the duration of current stimulation into suitable ranges. It is the very interesting bidirectional periodic adjustment phenomenon. These parameters are easily regulated in clinical practice, therefore, the results obtained in this paper may further help us to understand the treatment mechanism of the epilepsy seizure.

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.

Control of Absence Seizures by the Thalamic Feed-Forward Inhibition

As a subtype of idiopathic generalized epilepsies, absence epilepsy is believed to be caused by pathological interactions within the corticothalamic (CT) system. Using a biophysical mean-field model of the CT system, we demonstrate here that the feed-forward inhibition (FFI) in thalamus, i.e., the pathway from the cerebral cortex (Ctx) to the thalamic reticular nucleus (TRN) and then to the specific relay nuclei (SRN) of thalamus that are also directly driven by the Ctx, may participate in controlling absence seizures. In particular, we show that increasing the excitatory Ctx-TRN coupling strength can significantly suppress typical electrical activities during absence seizures. Further, investigation demonstrates that the GABA_A- and GABA_B-mediated inhibitions in the TRN-SRN pathway perform combination roles in the regulation of absence seizures. Overall, these results may provide an insightful mechanistic understanding of how the thalamic FFI serves as an intrinsic regulator contributing to the control of absence seizures.

Predictive modeling of EEG time series for evaluating surgery targets in epilepsy patients

During the last 20 years, predictive modeling in epilepsy research has largely been concerned with the prediction of seizure events, whereas the inference of effective brain targets for resective surgery has received surprisingly little attention. In this exploratory pilot study, we describe a distributional clustering framework for the modeling of multivariate time series and use it to predict the effects of brain surgery in epilepsy patients. By analyzing the intracranial EEG, we demonstrate how patients who became seizure free after surgery are clearly distinguished from those who did not. More specifically, for 5 out of 7 patients who obtained seizure freedom (= Engel class I) our method predicts the specific collection of brain areas that got actually resected during surgery to yield a markedly lower posterior probability for the seizure related clusters, when compared to the resection of random or empty collections. Conversely, for 4 out of 5 Engel class III/IV patients who still suffer from postsurgical seizures, performance of the actually resected collection is not significantly better than performances displayed by random or empty collections. As the number of possible collections ranges into billions and more, this is a substantial contribution to a problem that today is still solved by visual EEG inspection. Apart from epilepsy research, our clustering methodology is also of general interest for the analysis of multivariate time series and as a generative model for temporally evolving functional networks in the neurosciences and beyond. Hum Brain Mapp 38:2509-2531, 2017. © 2017 Wiley Periodicals, Inc.

Predicting neurosurgical outcomes in focal epilepsy patients using computational modelling

See Eissa and Schevon (doi:10.1093/aww332) for a scientific commentary on this article.

Surgery can be a last resort for patients with intractable, medically refractory epilepsy. For many of these patients, however, there is substantial risk that the surgery will be ineffective. The prediction of who is likely to benefit from a surgical approach is crucial for being able to inform patients better, conduct principled prospective clinical trials, and ultimately tailor therapeutic approaches to these patients more effectively. Dynamical computational models, informed with patient data, can be used to make predictions and give mechanistic insight. In this study, we develop patient-specific dynamical network models of epileptogenic cortex. We infer the network connectivity matrix from non-seizure electrographic recordings of patients and use these connectivity matrices as the network structure in our model. The model simulates the dynamics of a bi-stable switch at every node in this network, meaning that every node starts in a background state, but has the ability to transit to a co-existing seizure state. Whether a transition happens in a node is partly determined by the stochastic nature of the input to the node, but also by the input the node receives from other connected nodes in the network. By conducting simulations with such a model, we can detect the average transition time for nodes in a given network, and therefore define nodes with a short transition time as highly epileptogenic. In a retrospective study, we found that in some patients the regions with high epileptogenicity in the model overlap with those identified clinically as the seizure onset zone. Moreover, it was found that the resection of these regions in the model reduces the overall likelihood of a seizure. Following removal of these regions in the model, we predicted surgical outcomes and compared these to actual patient outcomes. Our predictions were found to be 81.3% accurate on a dataset of 16 patients with intractable epilepsy. Intriguingly, in patients with unsuccessful outcomes, the proposed computational approach is able to suggest alternative resection sites. The model presented here gives mechanistic insight as to why surgery may be unsuccessful in some patients. This may aid clinicians in presurgical evaluation by providing a tool to explore various surgical options, offering complementary information to existing clinical techniques.

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.

Virtual Cortical Resection Reveals Push-Pull Network Control Preceding Seizure Evolution

In ∼20 million people with drug-resistant epilepsy, focal seizures originating in dysfunctional brain networks will often evolve and spread to surrounding tissue, disrupting function in otherwise normal brain regions. To identify network control mechanisms that regulate seizure spread, we developed a novel tool for pinpointing brain regions that facilitate synchronization in the epileptic network. Our method measures the impact of virtually resecting putative control regions on synchronization in a validated model of the human epileptic network. By applying our technique to time-varying functional networks, we identified brain regions whose topological role is to synchronize or desynchronize the epileptic network. Our results suggest that greater antagonistic push-pull interaction between synchronizing and desynchronizing brain regions better constrains seizure spread. These methods, while applied here to epilepsy, are generalizable to other brain networks and have wide applicability in isolating and mapping functional drivers of brain dynamics in health and disease.

Optimally controlling the human connectome: the role of network topology

To meet ongoing cognitive demands, the human brain must seamlessly transition from one brain state to another, in the process drawing on different cognitive systems. How does the brain's network of anatomical connections help facilitate such transitions? Which features of this network contribute to making one transition easy and another transition difficult? Here, we address these questions using network control theory. We calculate the optimal input signals to drive the brain to and from states dominated by different cognitive systems. The input signals allow us to assess the contributions made by different brain regions. We show that such contributions, which we measure as energy, are correlated with regions' weighted degrees. We also show that the network communicability, a measure of direct and indirect connectedness between brain regions, predicts the extent to which brain regions compensate when input to another region is suppressed. Finally, we identify optimal states in which the brain should start (and finish) in order to minimize transition energy. We show that the optimal target states display high activity in hub regions, implicating the brain's rich club. Furthermore, when rich club organization is destroyed, the energy cost associated with state transitions increases significantly, demonstrating that it is the richness of brain regions that makes them ideal targets.

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.

Critical Roles of the Direct GABAergic Pallido-cortical Pathway in Controlling Absence Seizures

The basal ganglia (BG), serving as an intermediate bridge between the cerebral cortex and thalamus, are believed to play crucial roles in controlling absence seizure activities generated by the pathological corticothalamic system. Inspired by recent experiments, here we systematically investigate the contribution of a novel identified GABAergic pallido-cortical pathway, projecting from the globus pallidus externa (GPe) in the BG to the cerebral cortex, to the control of absence seizures. By computational modelling, we find that both increasing the activation of GPe neurons and enhancing the coupling strength of the inhibitory pallido-cortical pathway can suppress the bilaterally synchronous 2-4 Hz spike and wave discharges (SWDs) during absence seizures. Appropriate tuning of several GPe-related pathways may also trigger the SWD suppression, through modulating the activation level of GPe neurons. Furthermore, we show that the previously discovered bidirectional control of absence seizures due to the competition between other two BG output pathways also exists in our established model. Importantly, such bidirectional control is shaped by the coupling strength of this direct GABAergic pallido-cortical pathway. Our work suggests that the novel identified pallido-cortical pathway has a functional role in controlling absence seizures and the presented results might provide testable hypotheses for future experimental studies.

Hippocampal closed-loop modeling and implications for seizure stimulation design

OBJECTIVE

Traditional hippocampal modeling has focused on the series of feedforward synapses known as the trisynaptic pathway. However, feedback connections from CA1 back to the hippocampus through the entorhinal cortex (EC) actually make the hippocampus a closed-loop system. By constructing a functional closed-loop model of the hippocampus, one may learn how both physiological and epileptic oscillations emerge and design efficient neurostimulation patterns to abate such oscillations.

APPROACH

Point process input-output models where estimated from recorded rodent hippocampal data to describe the nonlinear dynamical transformation from CA3 → CA1, via the schaffer-collateral synapse, and CA1 → CA3 via the EC. Each Volterra-like subsystem was composed of linear dynamics (principal dynamic modes) followed by static nonlinearities. The two subsystems were then wired together to produce the full closed-loop model of the hippocampus.

MAIN RESULTS

Closed-loop connectivity was found to be necessary for the emergence of theta resonances as seen in recorded data, thus validating the model. The model was then used to identify frequency parameters for the design of neurostimulation patterns to abate seizures.

SIGNIFICANCE

Deep-brain stimulation (DBS) is a new and promising therapy for intractable seizures. Currently, there is no efficient way to determine optimal frequency parameters for DBS, or even whether periodic or broadband stimuli are optimal. Data-based computational models have the potential to be used as a testbed for designing optimal DBS patterns for individual patients. However, in order for these models to be successful they must incorporate the complex closed-loop structure of the seizure focus. This study serves as a proof-of-concept of using such models to design efficient personalized DBS patterns for epilepsy.