Suppression of seizures based on the multi-coupled neural mass model

Closed-loop seizure modulation via extreme learning machine based extended state observer

Neuromodulation is a promising way in clinical treatment of epilepsy, but the existing methods cannot adjust stimulations according to patients' real-time reactions. Therefore, it is necessary to acquire a systematic and a scientific regulation method based on patients' real-time reactions. The linear active disturbance rejection control can adapt to complex epileptic dynamics and improve the epilepsy regulation, even if little model information is available, and various uncertainties and external disturbances exist. However, a linear extended state observer estimates the time-varying total disturbance with a steady-state error. To improve regulation, it is crucial to estimate the total disturbance in a more accurate manner. An extreme learning machine is capable of approximating any nonlinear function. Its initial parameter generation is more convenient, adjustable parameters are fewer, and learning speed is faster. Thus, a nonlinear time-varying function can be estimated more timely and accurately. Then, an extreme learning machine based extended state observer is proposed to get a more satisfactory total disturbance estimation and more desired closed-loop regulation. The convergence of the extreme learning machine based extended state observer is verified and the stability of the closed-loop system is analyzed. Numerical results show that the proposed extended state observer is much better than a linear extended state observer in estimating the total disturbance. It guarantees a more satisfied closed-loop neuromodulation.

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

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

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.

Model Predictive Control for Seizure Suppression Based on Nonlinear Auto-Regressive Moving-Average Volterra Model

This article investigates a closed-loop brain stimulation method based on model predictive control strategy to suppress epileptic seizures. A neural mass model (NMM), exhibiting the normal and various epileptic seizures by changing physiologically meaningful parameters, is used as a black-box model of the brain. Based on system identification, an auto-regressive moving-average Volterra model is established to reveal the relationship between stimulation and neuronal responses. Then, the model predictive control strategy is implemented based the Volterra model, which can generate an optimal stimulation waveform to eliminate epileptiform waves. The computational simulation results indicate the proposed closed-loop control strategy can optimize the stimulation waveform without particular knowledge of the physiological properties of the system. The robustness of the proposed control strategy to system disturbances makes it more appropriate for future clinical application.

Principal dynamic mode analysis of neural mass model for the identification of epileptic states

The detection of epileptic seizures in Electroencephalography (EEG) signals is significant for the diagnosis and treatment of epilepsy. In this paper, in order to obtain characteristics of various epileptiform EEGs that may differentiate different states of epilepsy, the concept of Principal Dynamic Modes (PDMs) was incorporated to an autoregressive model framework. First, the neural mass model was used to simulate the required intracerebral EEG signals of various epileptiform activities. Then, the PDMs estimated from the nonlinear autoregressive Volterra models, as well as the corresponding Associated Nonlinear Functions (ANFs), were used for the modeling of epileptic EEGs. The efficient PDM modeling approach provided physiological interpretation of the system. Results revealed that the ANFs of the 1st and 2nd PDMs for the auto-regressive input exhibited evident differences among different states of epilepsy, where the ANFs of the sustained spikes' activity encountered at seizure onset or during a seizure were the most differentiable from that of the normal state. Therefore, the ANFs may be characteristics for the classification of normal and seizure states in the clinical detection of seizures and thus provide assistance for the diagnosis of epilepsy.

Particle swarm optimization algorithm based parameters estimation and control of epileptiform spikes in a neural mass model

This paper proposes an epilepsy detection and closed-loop control strategy based on Particle Swarm Optimization (PSO) algorithm. The proposed strategy can effectively suppress the epileptic spikes in neural mass models, where the epileptiform spikes are recognized as the biomarkers of transitions from the normal (interictal) activity to the seizure (ictal) activity. In addition, the PSO algorithm shows capabilities of accurate estimation for the time evolution of key model parameters and practical detection for all the epileptic spikes. The estimation effects of unmeasurable parameters are improved significantly compared with unscented Kalman filter. When the estimated excitatory-inhibitory ratio exceeds a threshold value, the epileptiform spikes can be inhibited immediately by adopting the proportion-integration controller. Besides, numerical simulations are carried out to illustrate the effectiveness of the proposed method as well as the potential value for the model-based early seizure detection and closed-loop control treatment design.