Modelling Effects of Electrical Stimulation on Seizure BENG/BGGN - - PowerPoint PPT Presentation

Modelling Effects of Electrical Stimulation on Seizure

BENG/BGGN 260 final project By: Carissa Gunawan

Introduction

• Seizure is caused by abnormal

excessive synchronous neural activity in the brain

• Therapy through drugs, VNS,

and DBS

• Advantage:
• Reversible
• More localized area

Fisher, Robert S., and Ana Luisa Velasco. "Electrical brain stimulation for epilepsy." Nature Reviews Neurology 10.5 (2014): 261-270. Fisher, Robert S., et al. "Epileptic seizures and epilepsy: definitions proposed by the International League Against Epilepsy (ILAE) and the International Bureau for Epilepsy (IBE)." Epilepsia 46.4 (2005): 470-472.

Purpose

• The goal of this project is

to find out the effects of electrical stimulation on various neural network during seizure

• Expected result:
• Stop over excitation during

seizure

• alter activity of normal

network

DeGiorgio, Christopher M., and Scott E. Krahl. "Neurostimulation for drug-resistant epilepsy." Continuum: Lifelong Learning in Neurology 19.3 Epilepsy (2013): 743.

Stimulation model

• Extracellular medium assumptions:
• Homogeneous electrical property
• Homogeneous density
• Point source
• Biphasic pulse sequence
• f=120-180Hz, pw= 0.06-0.2ms, V=1-5V
• Relationship of current with distance
• 𝐽 𝑢, 𝑆 = ' (

)*+,

R=distance and ρ=resistance/mm

Monfared, Omid, et al. "Electrical stimulation of neural tissue modeled as a cellular composite: Point Source electrode in an isotropic tissue." Engineering in Medicine and Biology Society (EMBC), 2014 36th Annual International Conference of the IEEE. IEEE, 2014.

mm mm (b) 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 −1 −0.8 −0.6 −0.4 −0.2 0.2 0.4 0.6 0.8 1 0.5 1 1.5 2 2.5 3 3.5 4 4.5 2 4 6 8 10 12 14 16 18 20 −50 −40 −30 −20 −10 10 20 30 40 50 Ie(t) (mA) t (ms) (a)

Neuron model

• Model: Hodgkin and Huxley
• Sodium channel, potassium channel, chlorine channel, Na-K

pump, leak currents, glia cells, inhibitory synapse, and excitatory synapse

• Assumptions:
• m is fast as compared to the voltage change
• Total amount of sodium ion is conserved
• Spherical cell body

Cressman, John R., et al. "The influence of sodium and potassium dynamics on excitability, seizures, and the stability of persistent states:

• I. Single neuron dynamics." Journal of computational neuroscience 26.2 (2009): 159-170.

Neuron model: Ionic Currents

Leak Current: 𝐽-. = 𝑕-. 𝑛1 2ℎ 𝑊 − 𝐹-. + 𝑕-.8 𝑊 − 𝐹-. 𝐽9 = 𝑕9𝑜) 𝑊 − 𝐹9 + 𝑕98 𝑊 − 𝐹9 𝐽;< = 𝑕;<8 𝑊 − 𝐹;< Current by Na-K pump, glia cells, and diffusion: 𝐽=>?= = 𝜍 1 + exp 25 − 𝑂𝑏 I 3 × 1 1 + exp 5.5 − 𝐿 N 𝐽O<I. = 𝐻O<I. 1 + exp 18 − 𝐿 N 2.5 𝐽RISS = 𝜗 𝐿 N − 𝑙1

Barreto, Ernest, and John R. Cressman. "Ion concentration dynamics as a mechanism for neuronal bursting." Journal of biological physics 37.3 (2011): 361-373.

Neuron model: Change in [K]o and [Na]i

Change in external potassium and internal sodium concentration 𝜐 𝑒 𝐿 N 𝑒𝑢 = 𝛿𝛾𝐽9 − 2𝛾𝐽=>?= − 𝐽O<I. − 𝐽RISS>ZIN[ 𝜐 𝑒 𝑂𝑏 I 𝑒𝑢 = −𝛿𝐽-. − 3𝐽=>?= External Sodium and Internal Potassium concentration: 𝑂𝑏 N = 144 − 𝛾 𝑂𝑏 I − 18 𝐿 I = 140 + (18 − 𝑂𝑏 I)

Barreto, Ernest, and John R. Cressman. "Ion concentration dynamics as a mechanism for neuronal bursting." Journal of biological physics 37.3 (2011): 361-373.

Bifurcation factor

• Bifurcation occurs depending on the dynamics of potassium and sodium
• This could be because of:
• Increase in local K concentration due to drugs or physical injury
• Genetic disease that alter the capacity of glia cell

5 10 15 20 25 30 −100 −50 50 Vm (mV) 5 10 15 20 25 30 10 20 30 40 [K]o (mM) 5 10 15 20 25 30 26 28 30 32 34 [Na]i (mM) time (s) 5 10 15 20 25 30 −100 −50 50 100 Vm(mV) 5 10 15 20 25 30 2 4 6 8 10 [K]o(mM) 5 10 15 20 25 30 16 17 18 19 [Na]i(mM) time (s)

Result: Bifurcation of one neuron

• Bifurcation start when [K]o in normal condition increase
• Surrounding neurons in this environment also conduct action potential at the

same time, making synchronous excitation

• Stimulation reduce synchrony by changing its frequency

5 10 15 20 25 30 −100 −50 50 100 Vm(mV) 5 10 15 20 25 30 2 4 6 8 10 [K]o(mM) 5 10 15 20 25 30 16 17 18 19 [Na]i(mM) time (s)

22.9 23 23.1 23.2 23.3 23.4 23.5 23.6 −80 −60 −40 −20 20 40 60 Vm (mV) time(s) without stimulation with stimulation

Network formation

• Independent neural network
• Dependent neural network both inhibitory and excitatory with

respect to each other

• Random neural network

Beverlin II, Bryce, et al. "Dynamical changes in neurons during seizures determine tonic to clonic shift." Journal

• f computational neuroscience 33.1

(2012): 41-51.

Discussion: Effects on networks

• It takes time for electrical

stimulation to unsynchronized the excitation

• It is not guaranteed and all the

neuron will be not synchronize with each other

• Excitation pattern is less

synchronized when neurons in a network depends on each other

• In a totally random network, it is

hard to find condition that leads to seizure

5 10 15 20 25 30 −2000 −1000 1000 2000 Vsum(mV) t (s) (a) 5 10 15 20 25 30 −2000 −1000 1000 2000 Vsum(mV) t (s) (b) 5 10 15 20 25 30 −2000 −1000 1000 2000 Vsum(mV) t (s) (c)

(a) Independent network (b) dependent network (c) random network

EEG model

• Assuming all the affected neurons are the only target
• Attenuation and noise
• Amplitude of 10-100 µV
• Noise consist of noise from other neuron
• Random white noise
• Sampled at a sampling frequency 200-2000 Hz

Hämäläinen, Matti, et al. "Magnetoencephalography—theory, instrumentation, and applications to noninvasive studies of the working human brain." Reviews of modern Physics 65.2 (1993): 413. Aurlien, H., et al. "EEG background activity described by a large computerized database." Clinical Neurophysiology 115.3 (2004): 665-673.

Attenuation Amplification and Sampling + Noise from other neurons raw_eeg EEG

Result: EEG model

Raw EEG Constructed EEG (a) Independent network (b) dependent network (c) random network

5 10 15 20 25 30 −2000 −1000 1000 2000 Vsum(mV) t (s) (a) 5 10 15 20 25 30 −2000 −1000 1000 2000 Vsum(mV) t (s) (b) 5 10 15 20 25 30 −2000 −1000 1000 2000 Vsum(mV) t (s) (c) 5 10 15 20 25 30 −10 10 20 30 Vsum(µV) (a) 5 10 15 20 25 30 −20 20 40 Vsum(µV) (b) 5 10 15 20 25 30 −10 10 20 30 Vsum(µV) time(s) (c)

5 10 15 20 25 30 −10 10 20 30 Vsum(µV) (a) 5 10 15 20 25 30 −20 20 40 Vsum(µV) (b) 5 10 15 20 25 30 −10 10 20 30 Vsum(µV) time(s) (c)

Discussion: EEG comparison

Conclusion

• The simulated EEG model have a similar behavior as the EEG

recording in the study.

• Electrical stimulation reduce seizure
• Part where seizure did not happen have some changes when

electrical stimulation was applied

• Change in local concentration will most likely cause seizure in area

where neuron are less dependent on each other

• Synchrony can be reduced if they are dependent on each other

Future direction

• Test with more realistic neuron model
• The one with Ca2+ activated K+ channel
• Test with a more realistic neural network
• Variable inhibition and excitation strength
• Test with other conditions that cause seizure
• Example abnormal Ca2+ concentration
• Glia’s failure to regulate K+ ions

Thank You