Tracking neuronal dynamics during seizures
© Ullah and Schiff; licensee BioMed Central Ltd. 2010
Published: 20 July 2010
We demonstrate that a meaningful estimation of the neuronal network dynamics from single measured variable is now possible by incorporating biophysical neuronal models into a model-based control framework. Specifically, we assimilate noisy membrane potential (V) measurements from individual hippocampal neurons to reconstruct the dynamics of networks of these cells, their extracellular microenvironment, and the activities of different neuronal types during seizures.
We used two-compartmental models for the pyramidal cells (PCs) and interneurons (INs): a cellular compartment and the surrounding extracellular microenvironment. The membrane potentials of both cells were modeled by Hodgkin-Huxley equations containing sodium, potassium, calcium-gated potassium, and leak currents. The current equations were augmented with dynamic variables representing the intra- and extracellular ion concentrations (K+, Na+, and Ca2+). To estimate and track the dynamics of the neuronal networks, we applied a nonlinear ensemble version of the Kalman filter, the Unscented Kalman Filter (UKF) . Details of the model and UKF implementation can be found in [2, 3].
In Fig. 1 we show an example where we use V measurements to estimate the rest of the dynamics of CA1 PC. Fig. 1A shows an intracellular recording from a PC during seizures, and plot the estimated [K]o in Fig. 1B. As is clear from the figure the [K]o oscillates as the cell goes into and out of seizures. The [K]o begins to rise as the cell enters seizures and peaks with the maximal firing frequency, followed by decreasing [K]o as the firing rate decreases and the seizure terminates. Higher [K]o makes the PC more excitable by raising the reversal potential for K+ currents. The increased K+ reversal potential causes the cell to burst fire spontaneously. Changes in the [Na]i are closely coupled with the changes of [K]o (Fig. 1C). As shown in panels (1D-F) we reconstructed the parameters controlling the microenvironment of the cell. These parameters included the diffusion constant of K+ in the extracellular space, K+ buffering strength of glia, and K+ concentration in the reservoir of the perfusing solution in vitro (or in the vasculature in vivo) during seizures.
- Kalman RE: A new approach to linear filtering and prediction problems. Trans ASME J Basic Eng. 1960, 82: 35-45.View ArticleGoogle Scholar
- Cressman JR, Ullah G, Ziburkus J, Schiff SJ, Barreto E: The influence of sodium and potassium dynamics on excitability, seizures, and the stability of persistent states: I. Single neuron dynamics. J Comput Neurosci. 2009, 26: 159-170. 10.1007/s10827-008-0132-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Ullah G, Schiff SJ: Tracking and control of neuronal Hodgkin-Huxley dynamics. Phys Rev E Stat Nonlin Soft Matter Phys. 2009, 79: 040901-10.1103/PhysRevE.79.040901.PubMed CentralView ArticlePubMedGoogle Scholar
- Ziburkus J, Cressman JR, Barreto E, Schiff SJ: Interneuron and pyramidal cell interplay during in vitro seizure-like events. J Neurophysiol. 2006, 95: 3948-3954. 10.1152/jn.01378.2005.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd.