scholarly journals Extracranial estimation of neural mass model parameters using the Unscented Kalman Filter

2017 ◽  
Author(s):  
Lara Escuain-Poole ◽  
Jordi Garcia-Ojalvo ◽  
Antonio J. Pons

AbstractData assimilation, defined as the fusion of data with preexisting knowledge, is particularly suited to elucidating underlying phenomena from noisy/insufficient observations. Although this approach has been widely used in diverse fields, only recently have efforts been directed to problems in neuroscience, using mainly intracranial data and thus limiting its applicability to invasive measurements involving electrode implants. Here we intend to apply data assimilation to non-invasive electroencephalography (EEG) measurements to infer brain states and their characteristics. For this purpose, we use Kalman filtering to combine synthetic EEG data with a coupled neural-mass model together with Ary’s model of the head, which projects intracranial signals onto the scalp. Our results show that using several extracranial electrodes allows to successfully estimate the state and parameters of the neural masses and their interactions, whereas one single electrode provides only a very partial and insufficient view of the system. The superiority of using multiple extracranial electrodes over using only one, be it intra- or extracranial, is shown over a wide variety of dynamical behaviours. Our results show potential towards future clinical applications of the method.Author SummaryTo completely understand brain function, we will need to integrate experimental information into a consistent theoretical framework. Invasive techniques as EcoG recordings, together with models that describe the brain at the mesoscale, provide valuable information about the brain state and its dynamical evolution when combined with techniques coming from control theory, such as the Kalman filter. This method, which is specifically designed to deal with systems with noisy or imperfect data, combines experimental data with theoretical models assuming Bayesian inference. So far, implementations of the Kalman filter have not been suited for non-invasive measures like EEG. Here we attempt to overcome this situation by introducing a model of the head that allows to transfer the intracranial signals produced by a mesoscopic model to the scalp in the form of EEG recordings. Our results show the advantages of using multichannel EEG recordings, which are extended in space and allow to discriminate signals produced by the interaction of coupled columns. The extension of the Kalman method presented here can be expected to expand the applicability of the technique to all situations where EEG recordings are used, including the routine monitoring of illnesses or rehabilitation tasks, brain-computer interface protocols, and transcranial stimulation.

Author(s):  
Sheikh Md. Rabiul Islam ◽  
◽  
Md. Shakibul Islam ◽  

The electroencephalogram (EEG) is an electrophysiological monitoring strategy that records the spontaneous electrical movement of the brain coming about from ionic current inside the neurons of the brain. The importance of the EEG signal is mainly the diagnosis of different mental and brain neurodegenerative diseases and different abnormalities like seizure disorder, encephalopathy, dementia, memory problem, sleep disorder, stroke, etc. The EEG signal is very useful for someone in case of a coma to determine the level of brain activity. So, it is very important to study EEG generation and analysis. To reduce the complexity of understanding the pathophysiological mechanism of EEG signal generation and their changes, different simulation-based EEG modeling has been developed which are based on anatomical equivalent data. In this paper, Instead of a detailed model a neural mass model has been used to implement different simulation-based EEG models for EEG signal generation which refers to the simplified and straightforward method. This paper aims to introduce obtained EEG signals of own implementation of the Lopes da Silva model, Jansen-Rit model, and Wendling model in Simulink and to compare characteristic features with real EEG signals and better understanding the EEG abnormalities especially the seizure-like signal pattern.


2017 ◽  
Vol 28 (01) ◽  
pp. 1750027 ◽  
Author(s):  
Zhen Ma

Electroencephalography (EEG) is an important method to investigate the neurophysiological mechanism underlying epileptogenesis to identify new therapies for the treatment of epilepsy. The neurophysiologically based neural mass model (NMM) can build a bridge between signal processing and neurophysiology, which can be used as a platform to explore the neurophysiological mechanism of epileptogenesis. Most EEG signals cannot be regarded as the outputs of a single NMM with identical model parameters. The outputs of NMM are simple because the diversity of neural signals in the same NMM is ignored. To improve the simulation of EEG signals, a multiple NMM is proposed, the output of which is the linear combination of the outputs of all NMMs. The NMM number is not fixed and is minimized under the premise of guaranteeing the fitting effect. Orthogonal matching pursuit is used to solve a constrained [Formula: see text] norm minimization problem for NMM number and the strength of every NMM. The results showed that the NMM number was significantly lower during the ictal period than during the interictal period, and the strength of major NMMs increased. This indicates that neural masses fuse into fewer larger neural masses with greater strength. The distribution of excitatory and inhibitory strength during the ictal and interictal periods was similar, whereas the excitation/inhibition ratio was higher during the ictal period than during the interictal period.


2017 ◽  
Vol 29 (2) ◽  
pp. 485-501 ◽  
Author(s):  
Xian Liu ◽  
Jing Gao ◽  
Guan Wang ◽  
Zhi-Wang Chen

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


2021 ◽  
Author(s):  
Viktor Sip ◽  
Spase Petkoski ◽  
Meysam Hashemi ◽  
Viktor K Jirsa

Model-based data analysis of whole-brain dynamics links the observed data to model parameters in a network of neural masses. In recent years a special focus was placed on the role of regional variance of model parameters for the emergent activity. Such analyses however depend on the properties of the employed neural mass model, which is often obtained through a series of major simplifications or analogies. Here we propose a data-driven approach where the neural mass model needs not to be specified. Building on the recent progresses in identification of dynamical systems with neural networks, we propose a method to infer from the functional data both the neural mass model representing the regional dynamics as well as the region- and subject-specific parameters, while respecting the known network structure. We demonstrate on two synthetic data sets that our method is able to recover the original model parameters, and that the trained generative model produces dynamics resembling the training data both on the regional level and on the whole-brain level. The present approach opens a novel way to the analysis of resting-state fMRI with possible applications in understanding the changes of whole-brain dynamics during aging or in neurodegenerative diseases.


2021 ◽  
Author(s):  
Edmundo Lopez-Sola ◽  
Roser Sanchez-Todo ◽  
Èlia Lleal ◽  
Elif Köksal-Ersöz ◽  
Maxime Yochum ◽  
...  

The prospect of personalized computational modeling in neurological disorders, and in particular in epilepsy, is poised to revolutionize the field. Work in the last two decades has demonstrated that neural mass models (NMM) can realistically reproduce and explain epileptic seizure transitions as recorded by electrophysiological methods (EEG, SEEG). In previous work, advances were achieved by i) increasing excitation in NMM and ii) heuristically varying network inhibitory coupling parameters or, equivalently, inhibitory synaptic gains. Based on those studies, we provide here a laminar neural mass model capable of realistically reproducing the electrical activity recorded by SEEG in the epileptogenic zone during interictal to ictal states. With the exception of the external noise input onto the pyramidal cell population, the model dynamics are autonomous --- all model parameters are static. By setting the system at a point close to bifurcation, seizure-like transitions are generated, including pre-ictal spikes, low voltage fast activity, and ictal rhythmic activity. A novel element in the model is a physiologically plausible algorithm for chloride accumulation dynamics: the gain of GABAergic post-synaptic potentials is modulated by the pathological accumulation of Cl$^-$ in pyramidal cells, due to high inhibitory input and/or dysfunctional chloride transport. In addition, in order to simulate SEEG signals to compare with real recordings performed in epileptic patients, the NMM is embedded first in a layered model of the neocortex and then in a realistic physical model. We compare modeling results with data from four epilepsy patient cases. By including key pathophysiological mechanisms, the proposed framework captures succinctly the electrophysiological phenomenology observed in ictal states, paving the way for robust personalization methods using brain network models based on NMMs.


2021 ◽  
Author(s):  
Saba Tabatabaee ◽  
Fariba Bahrami ◽  
Mahyar Janahmadi

Increasing evidence has shown that excitatory neurons in the brain play a significant role in seizure generation. However, spiny stellate cells are cortical excitatory non-pyramidal neurons in the brain which their 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 introduced originally by Taylor and colleagues in 2014. 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 due to the EI. Our results indicate that the transition is described by a supercritical Hopf bifurcation which shapes the preictal activity in the model and suggests why before seizure onset, the amplitude and frequency of neural activities increase gradually. 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 the ictal activity by decreasing the preictal period. Thereafter, we considered both sensory and cortical periodic inputs to drive the model responses to various harmonic stimulations. Our results from the bifurcation analysis of the model suggest that the initial stage of the brain 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 nonlinear resonance behavior result from the preictal stage of the brain. These results can be considered as a step forward to a deeper understanding of the mechanisms underlying the transition from brain normal activities to epileptic activities.


2020 ◽  
Vol 16 (11) ◽  
pp. e1008430
Author(s):  
Elif Köksal Ersöz ◽  
Julien Modolo ◽  
Fabrice Bartolomei ◽  
Fabrice Wendling

Epilepsy is a dynamic and complex neurological disease affecting about 1% of the worldwide population, among which 30% of the patients are drug-resistant. Epilepsy is characterized by recurrent episodes of paroxysmal neural discharges (the so-called seizures), which manifest themselves through a large-amplitude rhythmic activity observed in depth-EEG recordings, in particular in local field potentials (LFPs). The signature characterizing the transition to seizures involves complex oscillatory patterns, which could serve as a marker to prevent seizure initiation by triggering appropriate therapeutic neurostimulation methods. To investigate such protocols, neurophysiological lumped-parameter models at the mesoscopic scale, namely neural mass models, are powerful tools that not only mimic the LFP signals but also give insights on the neural mechanisms related to different stages of seizures. Here, we analyze the multiple time-scale dynamics of a neural mass model and explain the underlying structure of the complex oscillations observed before seizure initiation. We investigate population-specific effects of the stimulation and the dependence of stimulation parameters on synaptic timescales. In particular, we show that intermediate stimulation frequencies (>20 Hz) can abort seizures if the timescale difference is pronounced. Those results have the potential in the design of therapeutic brain stimulation protocols based on the neurophysiological properties of tissue.


2019 ◽  
Author(s):  
Andrea Ceni ◽  
Simona Olmi ◽  
Alessandro Torcini ◽  
David Angulo-Garcia

Coupling among neural rhythms is one of the most important mechanisms at the basis of cognitive processes in the brain. In this study we consider a neural mass model, rigorously obtained from the microscopic dynamics of an inhibitory spiking network with exponential synapses, able to autonomously generate collective oscillations (COs). These oscillations emerge via a super-critical Hopf bifurcation, and their frequencies are controlled by the synaptic time scale, the synaptic coupling and the excitability of the neural population. Furthermore, we show that two inhibitory populations in a master-slave configuration with different synaptic time scales can display various collective dynamical regimes: namely, damped oscillations towards a stable focus, periodic and quasi-periodic oscillations, and chaos. Finally, when bidirectionally coupled the two inhibitory populations can exhibit different types of θ-γ cross-frequency couplings (CFCs): namely, phase-phase and phase-amplitude CFC. The coupling between θ and γ COs is enhanced in presence of a external θ forcing, reminiscent of the type of modulation induced in Hippocampal and Cortex circuits via optogenetic drive.In healthy conditions, the brain’s activity reveals a series of intermingled oscillations, generated by large ensembles of neurons, which provide a functional substrate for information processing. How single neuron properties influence neuronal population dynamics is an unsolved question, whose solution could help in the understanding of the emergent collective behaviors arising during cognitive processes. Here we consider a neural mass model, which reproduces exactly the macroscopic activity of a network of spiking neurons. This mean-field model is employed to shade some light on an important and ubiquitous neural mechanism underlying information processing in the brain: the θ-γ cross-frequency coupling. In particular, we will explore in detail the conditions under which two coupled inhibitory neural populations can generate these functionally relevant coupled rhythms.


Sign in / Sign up

Export Citation Format

Share Document