Distributed control in a mean-field cortical network model: implications for seizure suppression

Finding influential nodes in networks using pinning control: Centrality measures confirmed with electrochemical oscillators

The spatiotemporal organization of networks of dynamical units can break down resulting in diseases (e.g., in the brain) or large-scale malfunctions (e.g., power grid blackouts). Re-establishment of function then requires identification of the optimal intervention site from which the network behavior is most efficiently re-stabilized. Here, we consider one such scenario with a network of units with oscillatory dynamics, which can be suppressed by sufficiently strong coupling and stabilizing a single unit, i.e., pinning control. We analyze the stability of the network with hyperbolas in the control gain vs coupling strength state space and identify the most influential node (MIN) as the node that requires the weakest coupling to stabilize the network in the limit of very strong control gain. A computationally efficient method, based on the Moore-Penrose pseudoinverse of the network Laplacian matrix, was found to be efficient in identifying the MIN. In addition, we have found that in some networks, the MIN relocates when the control gain is changed, and thus, different nodes are the most influential ones for weakly and strongly coupled networks. A control theoretic measure is proposed to identify networks with unique or relocating MINs. We have identified real-world networks with relocating MINs, such as social and power grid networks. The results were confirmed in experiments with networks of chemical reactions, where oscillations in the networks were effectively suppressed through the pinning of a single reaction site determined by the computational method.

Developing control-theoretic objectives for large-scale brain dynamics and cognitive enhancement

The development of technologies for brain stimulation provides a means for scientists and clinicians to directly actuate the brain and nervous system. Brain stimulation has shown intriguing potential in terms of modifying particular symptom clusters in patients and behavioral characteristics of subjects. The stage is thus set for optimization of these techniques and the pursuit of more nuanced stimulation objectives, including the modification of complex cognitive functions such as memory and attention. Control theory and engineering will play a key role in the development of these methods, guiding computational and algorithmic strategies for stimulation. In particular, realizing this goal will require new development of frameworks that allow for controlling not only brain activity, but also latent dynamics that underlie neural computation and information processing. In the current opinion, we review recent progress in brain stimulation and outline challenges and potential research pathways associated with exogenous control of cognitive function.

Modeling seizures: From single neurons to networks

Dynamical system tools offer a complementary approach to detailed biophysical seizure modeling, with a high potential for clinical applications. This review describes the theoretical framework that provides a basis for theorizing certain properties of seizures and for their classification according to their dynamical properties at onset and offset. We describe various modeling approaches spanning different scales, from single neurons to large-scale networks. This narrative review provides an accessible overview of this field, including non-exhaustive examples of key recent works.

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.

Real-time Inference and Detection of Disruptive EEG Networks for Epileptic Seizures

Recent studies in brain science and neurological medicine paid a particular attention to develop machine learning-based techniques for the detection and prediction of epileptic seizures with electroencephalogram (EEG). As a noninvasive monitoring method to record brain electrical activities, EEG has been widely used for capturing the underlying dynamics of disruptive neuronal responses across the brain in real-time to provide clinical guidance in support of epileptic seizure treatments in practice. In this study, we introduce a novel dynamic learning method that first infers a time-varying network constituted by multivariate EEG signals, which represents the overall dynamics of the brain network, and subsequently quantifies its topological property using graph theory. We demonstrate the efficacy of our learning method to detect relatively strong synchronization (characterized by the algebraic connectivity metric) caused by abnormal neuronal firing during a seizure onset. The computational results for a realistic scalp EEG database show a detection rate of 93.6% and a false positive rate of 0.16 per hour (FP/h); furthermore, our method observes potential pre-seizure phenomena in some cases.

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.

Emergence of spatially periodic diffusive waves in small-world neuronal networks

It has been observed in experiment that the anatomical structure of neuronal networks in the brain possesses the feature of small-world networks. Yet how the small-world structure affects network dynamics remains to be fully clarified. Here we study the dynamics of a class of small-world networks consisting of pulse-coupled integrate-and-fire (I&F) neurons. Under stochastic Poisson drive, we find that the activity of the entire network resembles diffusive waves. To understand its underlying mechanism, we analyze the simplified regular-lattice network consisting of firing-rate-based neurons as an approximation to the original I&F small-world network. We demonstrate both analytically and numerically that, with strongly coupled connections, in the absence of noise, the activity of the firing-rate-based regular-lattice network spatially forms a static grating pattern that corresponds to the spatial distribution of the firing rate observed in the I&F small-world neuronal network. We further show that the spatial grating pattern with different phases comprise the continuous attractor of both the I&F small-world and firing-rate-based regular-lattice network dynamics. In the presence of input noise, the activity of both networks is perturbed along the continuous attractor, which gives rise to the diffusive waves. Our numerical simulations and theoretical analysis may potentially provide insights into the understanding of the generation of wave patterns observed in cortical networks.

The role of network structure and time delay in a metapopulation Wilson--Cowan model

We study the dynamics of a network Wilson--Cowan model (a system of connected Wilson--Cowan oscillators) for interacting excitatory and inhibitory neuron populations with time delays. Each node in this model corresponds to a population of neurons, including excitatory and inhibitory subpopulations, and hence it can be viewed as a metapopulation model. It is known that information transfer within each cortical area is not instantaneous, and therefore we consider a system of delay differential equations with two different kinds of discrete delays. We account for the time delay in information propagation between individual excitatory and inhibitory subpopulations at each node via intra-node time delays, and we account for time delay in information propagation between neuron populations at different nodes with inter-node time delays. The biologically relevant resting state solutions are oscillatory (stable limit cycles). After determining the influence of the coupling parameters between nodes, the intra-node delays, and the inter-node delays on the dynamics of the two coupled Wilson--Cowan oscillators, we then explore a variety of larger networks of 16 and 100 nodes, in order to determine how the network topology will influence time delayed Wilson--Cowan dynamics. We find that network structure can regularize or deregularize the dynamics, with networks of higher mean degree permitting stable limit cycles and networks with smaller mean degree yielding less regular dynamics (which may range from chaotic solutions, to solutions for which limit cycles collapse into steady states, which are biologically undesirable compared with the preferred stable limit cycles). Furthermore, heterogeneity in the degree distribution of the network (resulting from networks with nodes of varying degree) can result in asynchronous dynamics, even if at each node the local dynamics are that of a limit cycle, in contrast to the synchronization of dynamics between nodes seen when the degree of all nodes is equal. This suggests that homogeneous and well-connected networks permit robust limit cycles under time-delayed Wilson--Cowan dynamics, whereas heterogeneous or poorly connected networks may fail to provide such desirable dynamics, a phenomena akin to structural loss of neuron connections in neurodegenerative diseases.

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.

Stability of entrainment of a continuum of coupled oscillators

Complex natural and engineered systems are ubiquitous, and their behavior is challenging to characterize and control. We examine the design of the entrainment process for an uncountably infinite collection of coupled phase oscillators that are all subject to the same periodic driving signal. In the absence of coupling, an appropriately designed input can result in each oscillator attaining the frequency of the driving signal, with a phase offset determined by its natural frequency. We consider a special case of interacting oscillators in which the coupling tends to destabilize the phase configuration to which the driving signal would send the collection in the absence of coupling. In this setting, we derive stability results that characterize the trade-off between the effects of driving and coupling, and compare these results to the well-known Kuramoto model of a collection of free-running coupled oscillators.

Role of Graph Architecture in Controlling Dynamical Networks with Applications to Neural Systems

Networked systems display complex patterns of interactions between components. In physical networks, these interactions often occur along structural connections that link components in a hard-wired connection topology, supporting a variety of system-wide dynamical behaviors such as synchronization. While descriptions of these behaviors are important, they are only a first step towards understanding and harnessing the relationship between network topology and system behavior. Here, we use linear network control theory to derive accurate closed-form expressions that relate the connectivity of a subset of structural connections (those linking driver nodes to non-driver nodes) to the minimum energy required to control networked systems. To illustrate the utility of the mathematics, we apply this approach to high-resolution connectomes recently reconstructed from Drosophila, mouse, and human brains. We use these principles to suggest an advantage of the human brain in supporting diverse network dynamics with small energetic costs while remaining robust to perturbations, and to perform clinically accessible targeted manipulation of the brain's control performance by removing single edges in the network. Generally, our results ground the expectation of a control system's behavior in its network architecture, and directly inspire new directions in network analysis and design via distributed control.

Emerging Frontiers of Neuroengineering: A Network Science of Brain Connectivity

Neuroengineering is faced with unique challenges in repairing or replacing complex neural systems that are composed of many interacting parts. These interactions form intricate patterns over large spatiotemporal scales and produce emergent behaviors that are difficult to predict from individual elements. Network science provides a particularly appropriate framework in which to study and intervene in such systems by treating neural elements (cells, volumes) as nodes in a graph and neural interactions (synapses, white matter tracts) as edges in that graph. Here, we review the emerging discipline of network neuroscience, which uses and develops tools from graph theory to better understand and manipulate neural systems from micro- to macroscales. We present examples of how human brain imaging data are being modeled with network analysis and underscore potential pitfalls. We then highlight current computational and theoretical frontiers and emphasize their utility in informing diagnosis and monitoring, brain-machine interfaces, and brain stimulation. A flexible and rapidly evolving enterprise, network neuroscience provides a set of powerful approaches and fundamental insights that are critical for the neuroengineer's tool kit.

Emergence of coupling-induced oscillations and broken symmetries in heterogeneously driven nonlinear reaction networks

Many natural systems including the brain comprise coupled elements that are stimulated non-uniformly. In this paper we show that heterogeneously driven networks of excitatory-inhibitory units exhibit a diverse range of collective phenomena, including the appearance of spontaneous oscillations upon coupling quiescent elements. On varying the coupling strength a previously unreported transition is seen wherein the symmetries of the synchronization patterns in the stimulated and unstimulated groups undergo mutual exchange. The system also exhibits coexisting chaotic and non-chaotic attractors - a result that may be of interest in connection to earlier reports of varying degrees of chaoticity in the brain.

A network engineering perspective on probing and perturbing cognition with neurofeedback

Network science and engineering provide a flexible and generalizable tool set to describe and manipulate complex systems characterized by heterogeneous interaction patterns among component parts. While classically applied to social systems, these tools have recently proven to be particularly useful in the study of the brain. In this review, we describe the nascent use of these tools to understand human cognition, and we discuss their utility in informing the meaningful and predictable perturbation of cognition in combination with the emerging capabilities of neurofeedback. To blend these disparate strands of research, we build on emerging conceptualizations of how the brain functions (as a complex network) and how we can develop and target interventions or modulations (as a form of network control). We close with an outline of current frontiers that bridge neurofeedback, connectomics, and network control theory to better understand human cognition.

Human seizures couple across spatial scales through travelling wave dynamics

Epilepsy—the propensity toward recurrent, unprovoked seizures—is a devastating disease affecting 65 million people worldwide. Understanding and treating this disease remains a challenge, as seizures manifest through mechanisms and features that span spatial and temporal scales. Here we address this challenge through the analysis and modelling of human brain voltage activity recorded simultaneously across microscopic and macroscopic spatial scales. We show that during seizure large-scale neural populations spanning centimetres of cortex coordinate with small neural groups spanning cortical columns, and provide evidence that rapidly propagating waves of activity underlie this increased inter-scale coupling. We develop a corresponding computational model to propose specific mechanisms—namely, the effects of an increased extracellular potassium concentration diffusing in space—that support the observed spatiotemporal dynamics. Understanding the multi-scale, spatiotemporal dynamics of human seizures—and connecting these dynamics to specific biological mechanisms—promises new insights to treat this devastating disease.

The authors record both local and long-range neural activity during human epileptic seizures to study the underlying multi-scale dynamics. They find that coupling of activity across spatial scales increases during seizures through propagating waves that are fit by a model that combines neural activity and potassium concentration dynamics.

Network neuroscience

Despite substantial recent progress, our understanding of the principles and mechanisms underlying complex brain function and cognition remains incomplete. Network neuroscience proposes to tackle these enduring challenges. Approaching brain structure and function from an explicitly integrative perspective, network neuroscience pursues new ways to map, record, analyze and model the elements and interactions of neurobiological systems. Two parallel trends drive the approach: the availability of new empirical tools to create comprehensive maps and record dynamic patterns among molecules, neurons, brain areas and social systems; and the theoretical framework and computational tools of modern network science. The convergence of empirical and computational advances opens new frontiers of scientific inquiry, including network dynamics, manipulation and control of brain networks, and integration of network processes across spatiotemporal domains. We review emerging trends in network neuroscience and attempt to chart a path toward a better understanding of the brain as a multiscale networked system.

Optimal trajectories of brain state transitions

The complexity of neural dynamics stems in part from the complexity of the underlying anatomy. Yet how white matter structure constrains how the brain transitions from one cognitive state to another remains unknown. Here we address this question by drawing on recent advances in network control theory to model the underlying mechanisms of brain state transitions as elicited by the collective control of region sets. We find that previously identified attention and executive control systems are poised to affect a broad array of state transitions that cannot easily be classified by traditional engineering-based notions of control. This theoretical versatility comes with a vulnerability to injury. In patients with mild traumatic brain injury, we observe a loss of specificity in putative control processes, suggesting greater susceptibility to neurophysiological noise. These results offer fundamental insights into the mechanisms driving brain state transitions in healthy cognition and their alteration following injury.

Small-World Brain Networks Revisited

It is nearly 20 years since the concept of a small-world network was first quantitatively defined, by a combination of high clustering and short path length; and about 10 years since this metric of complex network topology began to be widely applied to analysis of neuroimaging and other neuroscience data as part of the rapid growth of the new field of connectomics. Here, we review briefly the foundational concepts of graph theoretical estimation and generation of small-world networks. We take stock of some of the key developments in the field in the past decade and we consider in some detail the implications of recent studies using high-resolution tract-tracing methods to map the anatomical networks of the macaque and the mouse. In doing so, we draw attention to the important methodological distinction between topological analysis of binary or unweighted graphs, which have provided a popular but simple approach to brain network analysis in the past, and the topology of weighted graphs, which retain more biologically relevant information and are more appropriate to the increasingly sophisticated data on brain connectivity emerging from contemporary tract-tracing and other imaging studies. We conclude by highlighting some possible future trends in the further development of weighted small-worldness as part of a deeper and broader understanding of the topology and the functional value of the strong and weak links between areas of mammalian cortex.

The Geometry of Plasticity-Induced Sensitization in Isoinhibitory Rate Motifs

A well-known phenomenon in sensory perception is desensitization, wherein behavioral responses to persistent stimuli become attenuated over time. In this letter, our focus is on studying mechanisms through which desensitization may be mediated at the network level and, specifically, how sensitivity changes arise as a function of long-term plasticity. Our principal object of study is a generic isoinhibitory motif: a small excitatory-inhibitory network with recurrent inhibition. Such a motif is of interest due to its overrepresentation in laminar sensory network architectures. Here, we introduce a sensitivity analysis derived from control theory in which we characterize the fixed-energy reachable set of the motif. This set describes the regions of the phase-space that are more easily (in terms of stimulus energy) accessed, thus providing a holistic assessment of sensitivity. We specifically focus on how the geometry of this set changes due to repetitive application of a persistent stimulus. We find that for certain motif dynamics, this geometry contracts along the stimulus orientation while expanding in orthogonal directions. In other words, the motif not only desensitizes to the persistent input, but heightens its responsiveness (sensitizes) to those that are orthogonal. We develop a perturbation analysis that links this sensitization to both plasticity-induced changes in synaptic weights and the intrinsic dynamics of the network, highlighting that the effect is not purely due to weight-dependent disinhibition. Instead, this effect depends on the relative neuronal time constants and the consequent stimulus-induced drift that arises in the motif phase-space. For tightly distributed (but random) parameter ranges, sensitization is quite generic and manifests in larger recurrent E-I networks within which the motif is embedded.

Computational modeling of neurostimulation in brain diseases

Neurostimulation as a therapeutic tool has been developed and used for a range of different diseases such as Parkinson's disease, epilepsy, and migraine. However, it is not known why the efficacy of the stimulation varies dramatically across patients or why some patients suffer from severe side effects. This is largely due to the lack of mechanistic understanding of neurostimulation. Hence, theoretical computational approaches to address this issue are in demand. This chapter provides a review of mechanistic computational modeling of brain stimulation. In particular, we will focus on brain diseases, where mechanistic models (e.g., neural population models or detailed neuronal models) have been used to bridge the gap between cellular-level processes of affected neural circuits and the symptomatic expression of disease dynamics. We show how such models have been, and can be, used to investigate the effects of neurostimulation in the diseased brain. We argue that these models are crucial for the mechanistic understanding of the effect of stimulation, allowing for a rational design of stimulation protocols. Based on mechanistic models, we argue that the development of closed-loop stimulation is essential in order to avoid inference with healthy ongoing brain activity. Furthermore, patient-specific data, such as neuroanatomic information and connectivity profiles obtainable from neuroimaging, can be readily incorporated to address the clinical issue of variability in efficacy between subjects. We conclude that mechanistic computational models can and should play a key role in the rational design of effective, fully integrated, patient-specific therapeutic brain stimulation.

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.

Cognitive network neuroscience

Network science provides theoretical, computational, and empirical tools that can be used to understand the structure and function of the human brain in novel ways using simple concepts and mathematical representations. Network neuroscience is a rapidly growing field that is providing considerable insight into human structural connectivity, functional connectivity while at rest, changes in functional networks over time (dynamics), and how these properties differ in clinical populations. In addition, a number of studies have begun to quantify network characteristics in a variety of cognitive processes and provide a context for understanding cognition from a network perspective. In this review, we outline the contributions of network science to cognitive neuroscience. We describe the methodology of network science as applied to the particular case of neuroimaging data and review its uses in investigating a range of cognitive functions including sensory processing, language, emotion, attention, cognitive control, learning, and memory. In conclusion, we discuss current frontiers and the specific challenges that must be overcome to integrate these complementary disciplines of network science and cognitive neuroscience. Increased communication between cognitive neuroscientists and network scientists could lead to significant discoveries under an emerging scientific intersection known as cognitive network neuroscience.

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.

Emergence from general anesthesia and the sleep-manifold

The electroencephalogram (EEG) during the re-establishment of consciousness after general anesthesia and surgery varies starkly between patients. Can the EEG during this emergence period provide a means of estimating the underlying biological processes underpinning the return of consciousness? Can we use a model to infer these biological processes from the EEG patterns? A frontal EEG was recorded from 84 patients. Ten patients were chosen for state-space analysis. Five showed archetypal emergences; which consisted of a progressive decrease in alpha power and increase peak alpha frequency before return of responsiveness. The five non-archetypal emergences showed almost no spectral EEG changes (even as the volatile general anesthetic decreased) and then an abrupt return of responsiveness. We used Bayesian methods to estimate the likelihood of an EEG pattern corresponding to the position of the patient on a 2-dimensional manifold in a state space of excitatory connection strength vs. change in intrinsic resting neuronal membrane conductivity. We could thus visualize the trajectory of each patient in the state-space during their emergence period. The patients who followed an archetypal emergence displayed a very consistent pattern; consisting of progressive increase in conductivity, and a temporary period of increased connection strength before return of responsiveness. The non-archetypal emergence trajectories remained fixed in a region of phase space characterized by a relatively high conductivity and low connection strength throughout emergence. This unexpected progressive increase in conductivity during archetypal emergence may be due to an abating of the surgical stimulus during this period. Periods of high connection strength could represent forays into dissociated consciousness, but the model suggests all patients reposition near the fold in the state space to take advantage of bi-stable cortical dynamics before transitioning to consciousness.

A Hamilton-Jacobi-Bellman approach for termination of seizure-like bursting

We use Hamilton-Jacobi-Bellman methods to find minimum-time and energy-optimal control strategies to terminate seizure-like bursting behavior in a conductance-based neural model. Averaging is used to eliminate fast variables from the model, and a target set is defined through bifurcation analysis of the slow variables of the model. This method is illustrated for a single neuron model and for a network model to illustrate its efficacy in terminating bursting once it begins. This work represents a numerical proof-of-concept that a new class of control strategies can be employed to mitigate bursting, and could ultimately be adapted to treat medically intractible epilepsy in patient-specific models.