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.

[1]  D. Luenberger Optimization by Vector Space Methods , 1968 .

[2]  F. Takens Detecting strange attractors in turbulence , 1981 .

[3]  J. Penry,et al.  Some Clinical and EEG Aspects of Benign Juvenile Myoclonic Epilepsy , 1984, Epilepsia.

[4]  A. Babloyantz,et al.  Low-dimensional chaos in an instance of epilepsy. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[5]  M. Sillanpää,et al.  Distribution of Seizure Types in an Epileptic Population , 1988, Epilepsia.

[6]  P. Rajna,et al.  Sensory Stimulation for Inhibition of Epileptic Seizures , 1989, Epilepsia.

[7]  Bengt Fornberg,et al.  A practical guide to pseudospectral methods: Introduction , 1996 .

[8]  A. Destexhe Spike-and-Wave Oscillations Based on the Properties of GABAB Receptors , 1998, The Journal of Neuroscience.

[9]  Olivia Freeman,et al.  Talking points personal outcomes approach: practical guide. , 2012 .

[10]  P. Robinson,et al.  Dynamics of large-scale brain activity in normal arousal states and epileptic seizures. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[11]  Hal Blumenfeld,et al.  Why do Seizures Cause Loss of Consciousness? , 2003, The Neuroscientist : a review journal bringing neurobiology, neurology and psychiatry.

[12]  THE IMPORTANCE OF MODELLING , 2004 .

[13]  S. Amari Dynamics of pattern formation in lateral-inhibition type neural fields , 1977, Biological Cybernetics.

[14]  F. L. D. Silva,et al.  Dynamics of non-convulsive epileptic phenomena modeled by a bistable neuronal network , 2004, Neuroscience.

[15]  F. H. Lopes da Silva,et al.  Evolving concepts on the pathophysiology of absence seizures: the cortical focus theory. , 2005, Archives of neurology.

[16]  Olaf Sporns,et al.  The Human Connectome: A Structural Description of the Human Brain , 2005, PLoS Comput. Biol..

[17]  H. Blumenfeld Cellular and Network Mechanisms of Spike‐Wave Seizures , 2005, Epilepsia.

[18]  V. Fuster,et al.  Evolving Concepts , 2005 .

[19]  Terence J O'Brien,et al.  Cellular and network mechanisms of genetically-determined absence seizures. , 2005, Thalamus & related systems.

[20]  John R. Terry,et al.  A unifying explanation of primary generalized seizures through nonlinear brain modeling and bifurcation analysis. , 2006, Cerebral cortex.

[21]  I. Scheffer,et al.  Electroclinical features of absence seizures in childhood absence epilepsy , 2006, Neurology.

[22]  Anders M. Dale,et al.  An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest , 2006, NeuroImage.

[23]  Eric T. Shea-Brown,et al.  Optimal Inputs for Phase Models of Spiking Neurons , 2006 .

[24]  Mark A Kramer,et al.  Bifurcation control of a seizing human cortex. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[25]  Alfredo Berardelli,et al.  Effects of repetitive transcranial magnetic stimulation on spike-and-wave discharges , 2007, Neuroscience Research.

[26]  H. Siebner,et al.  Simultaneous EEG‐fMRI in drug‐naive children with newly diagnosed absence epilepsy , 2008, Epilepsia.

[27]  P. Ossenblok,et al.  Onset and propagation of spike and slow wave discharges in human absence epilepsy: A MEG study , 2009, Epilepsia.

[28]  O Sporns,et al.  Predicting human resting-state functional connectivity from structural connectivity , 2009, Proceedings of the National Academy of Sciences.

[29]  John R. Terry,et al.  Derivation and analysis of an ordinary differential equation mean-field model for studying clinically recorded epilepsy dynamics. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  Olaf Sporns,et al.  What Is the Human Connectome , 2009 .

[31]  Fang-Cheng Yeh,et al.  Generalized ${ q}$-Sampling Imaging , 2010, IEEE Transactions on Medical Imaging.

[32]  V. Wedeen,et al.  Generalized -Sampling Imaging , 2010 .

[33]  T. A. Zang,et al.  Spectral Methods: Fundamentals in Single Domains , 2010 .

[34]  Sheng-Fu Liang,et al.  A closed-loop brain computer interface for real-time seizure detection and control , 2010, 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology.

[35]  Hamid Soltanian-Zadeh,et al.  Multi-area neural mass modeling of EEG and MEG signals , 2010, NeuroImage.

[36]  Demetrios N. Velis,et al.  Stimulation-based anticipation and control of state transitions in the epileptic brain , 2010, Epilepsy & Behavior.

[37]  Matthew N. DeSalvo,et al.  Dynamic Time Course of Typical Childhood Absence Seizures: EEG, Behavior, and Functional Magnetic Resonance Imaging , 2010, The Journal of Neuroscience.

[38]  Gerold Baier,et al.  A spatially extended model for macroscopic spike-wave discharges , 2011, Journal of Computational Neuroscience.

[39]  Albert-László Barabási,et al.  Controllability of complex networks , 2011, Nature.

[40]  Kaspar Anton Schindler,et al.  Forbidden ordinal patterns of periictal intracranial EEG indicate deterministic dynamics in human epileptic seizures , 2011, Epilepsia.

[41]  Mark A Kramer,et al.  Distributed control in a mean-field cortical network model: implications for seizure suppression. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[42]  Sheng-Fu Liang,et al.  Live demonstration: A portable closed-loop seizure controller in freely moving rats , 2012, 2012 IEEE Biomedical Circuits and Systems Conference (BioCAS).

[43]  M. Belluscio,et al.  Closed-Loop Control of Epilepsy by Transcranial Electrical Stimulation , 2012, Science.

[44]  M. Kramer,et al.  Epilepsy as a Disorder of Cortical Network Organization , 2012, The Neuroscientist : a review journal bringing neurobiology, neurology and psychiatry.

[45]  P. Taylor,et al.  The importance of modeling epileptic seizure dynamics as spatio-temporal patterns , 2012, Front. Physio..

[46]  Gerold Baier,et al.  Towards a large-scale model of patient-specific epileptic spike-wave discharges , 2012, Biological Cybernetics.

[47]  P. Taylor,et al.  Phase space approach for modeling of epileptic dynamics. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[48]  Oscar Benjamin,et al.  Seizure generation: The role of nodes and networks , 2012, Epilepsia.

[49]  Jr-Shin Li,et al.  Optimal Control of Inhomogeneous Ensembles , 2011, IEEE Transactions on Automatic Control.

[50]  Philippe Kahane,et al.  Probabilistic functional tractography of the human cortex , 2013, NeuroImage.

[51]  Walter J. Freeman,et al.  Open loop optogenetic control of simulated cortical epileptiform activity , 2013, Journal of Computational Neuroscience.

[52]  Jr-Shin Li,et al.  Control and Synchronization of Neuron Ensembles , 2011, IEEE Transactions on Automatic Control.

[53]  Justin Dauwels,et al.  A model of stimulus induced epileptic spike-wave discharges , 2013, 2013 IEEE Symposium on Computational Intelligence, Cognitive Algorithms, Mind, and Brain (CCMB).

[54]  Peng Li,et al.  The emergence of abnormal hypersynchronization in the anatomical structural network of human brain , 2013, NeuroImage.

[55]  Maurizio Corbetta,et al.  Resting-State Functional Connectivity Emerges from Structurally and Dynamically Shaped Slow Linear Fluctuations , 2013, The Journal of Neuroscience.

[56]  Dante R Chialvo,et al.  Brain organization into resting state networks emerges at criticality on a model of the human connectome. , 2012, Physical review letters.

[57]  Guillaume Charvet,et al.  Neural Adaptation to Responsive Stimulation: A Comparison of Auditory and Deep Brain Stimulation in a Rat Model of Absence Epilepsy , 2013, Brain Stimulation.

[58]  U. Stephani,et al.  A Computational Study of Stimulus Driven Epileptic Seizure Abatement , 2014, PloS one.

[59]  W. Stacey,et al.  On the nature of seizure dynamics. , 2014, Brain : a journal of neurology.

[60]  Derek Ruths,et al.  Control Profiles of Complex Networks , 2014, Science.

[61]  Marcus Kaiser,et al.  Structural connectivity based whole brain modelling in epilepsy , 2014, Journal of Neuroscience Methods.

[62]  Habib Benali,et al.  Relating Structure and Function in the Human Brain: Relative Contributions of Anatomy, Stationary Dynamics, and Non-stationarities , 2014, PLoS Comput. Biol..

[63]  Renaud Lopes,et al.  Intra-subject reliability of the high-resolution whole-brain structural connectome , 2014, NeuroImage.

[64]  Justin Dauwels,et al.  Optimal Control of an Epileptic Neural Population Model , 2014 .

[65]  Spatiotemporal Dynamics of Large-Scale Brain Activity , 2015 .