Biological Modeling of Neural Networks
|
|
- Melvyn Lawson
- 5 years ago
- Views:
Transcription
1 Week 3 part 1 : Rection of the Hodgkin-Huxley Model 3.1 From Hodgkin-Huxley to 2D Biological Modeling of Neural Netorks - Overvie: From 4 to 2 dimensions - MathDetour 1: Exploiting similarities - MathDetour 2: Separation of time scales 3.2 Phase Plane Analysis - Role of nullclines Week 3 Recing detail: To-dimensional neuron models Wulfram Gerstner EPFL, Lausanne, Sitzerland 3.3 Analysis of a 2D Neuron Model - constant input vs pulse input - MathDetour 3: Stability of fixed points 3.4 TypeI and II Neuron Models next eek!
2 Neuronal Dynamics 3.1. Revie :Hodgkin-Huxley Model cortical neuron It () u d u... Hodgkin-Huxley model Compartmental models
3 Neuronal Dynamics 3.1 Revie :Hodgkin-Huxley Model Week 2: Cell membrane contains - ion channels - ion pumps -70mV Dendrites (eek 4): Active processes? assumption: passive dendrite point neuron Na + spike generation action potential K + Ca 2+ Ions/proteins
4 Neuronal Dynamics 3.1. Revie :Hodgkin-Huxley Model 100 inside K Ka mv 0 Ion channels Ion pump Na outside ln nu ( ) nu ( ) u u u Reversal potential kt ion pumps concentration difference voltage difference q 2
5 Neuronal Dynamics 3.1. Revie: Hodgkin-Huxley Model 100 I Hodgkin and Huxley, 1952 inside Ka C mv g K 0 g Na g l Ion channels Ion pump Na outside stimulus I Na I K I leak 4 equations C g Na m 3 h( u E Na ) g K n 4 ( u E K ) g l ( u E l ) I ( t) = 4D system dh dm dn hn m h n0( m u0( ) u) ( u () u ) hn m n 0 (u) n (u) u u
6 Neuronal Dynamics 3.1. Overvie and aims Can e understand the dynamics of the HH model? - mathematical principle of Action Potential generation? - constant input current vs pulse input? - Types of neuron model (type I and II)? (next eek) - threshold behavior? (next eek) Rece from 4 to 2 equations Type I and type II models ramp input/ constant input f f-i curve f f-i curve I0 I0 I0
7 Neuronal Dynamics 3.1. Overvie and aims Can e understand the dynamics of the HH model? Rece from 4 to 2 equations
8 Neuronal Dynamics Quiz 3.1. A - A biophysical point neuron model ith 3 ion channels, each ith activation and inactivation, has a total number of equations equal to [ ] 3 or [ ] 4 or [ ] 6 or [ ] 7 or [ ] 8 or more
9 Neuronal Dynamics 3.1. Overvie and aims Toard a to-dimensional neuron model -Rection of Hodgkin-Huxley to 2 dimension -step 1: separation of time scales -step 2: exploit similarities/correlations
10 Neuronal Dynamics 3.1. Rection of Hodgkin-Huxley model stimulus I I I Na K leak C g m 3 h( u E ) g n 4 ( u E ) g ( u E ) I ( t) Na Na K K l l dm m m0( u) m 0 (u) (u) h m ( u) dh h h0( u) h 0 (u) m (u) h ( u) u u dn n n0( u) n ( u) MathDetour 1) dynamics of m are fast m( t) m ( u( t)) 0
11 Neuronal Dynamics MathDetour 3.1: Separation of time scales dx x c() t 1 To coupled differential equations 1 dx x h( y) dy Exercise 1 (eek 3) 2 f ( y) g( x) (later!) Separation of time scales x h( y) 1 2 Reced 1-dimensional system dy f ( y) g( h( y)) 2
12 Neuronal Dynamics 3.1. Rection of Hodgkin-Huxley model stimulus I I I Na K leak C g m 3 h( u E ) g n 4 ( u E ) g ( u E ) I ( t) Na Na K K l l dm m m0( u) m ( u) n 0 (u) n h (u) dh h h0( u) h ( u) h 0 (u) u m (u) u dn n n0( u) n ( u) 1) dynamics of m are fast m( t) m ( u( t)) 0 2) dynamics of h and n are similar
13 Neuronal Dynamics 3.1. Rection of Hodgkin-Huxley model Rection of Hodgkin-Huxley Model to 2 Dimension -step 1: separation of time scales -step 2: exploit similarities/correlations No!
14 Neuronal Dynamics 3.1. Rection of Hodgkin-Huxley model stimulus C g m 3 h( u E ) g n 4 ( u E ) g ( u E ) I ( t) Na Na K K l l 2) dynamics of h and n are similar 1 h( t) a n ( t) u MathDetour t 1 hrest nrest 0 h(t) n(t) t
15 Neuronal Dynamics Detour 3.1. Exploit similarities/correlations dynamics of h and n are similar h 1 h( t) a n( t ) 1-h Math. argument n 1 hrest nrest 0 h(t) n(t) t
16 Neuronal Dynamics Detour 3.1. Exploit similarities/correlations dynamics of h and n are similar 1 h( t) a n( t ) at rest n dh h h0( u) h ( u) 1 hrest nrest h(t) n(t) dn n n0( u) n ( u) 0 t
17 Neuronal Dynamics Detour 3.1. Exploit similarities/correlations dynamics of h and n are similar (i) Rotate coordinate system (ii) Suppress one coordinate (iii) Express dynamics in ne variable 1 h( t) a n( t) ( t) 1 hrest nrest n h(t) n(t) dh dn h h0( u) h n n0( u) n ( u) ( u) d 0 ( u ) eff ( u) 0 t
18 Neuronal Dynamics 3.1. Rection of Hodgkin-Huxley model I I Na K C g [ m( t)] 3 h( t)( u( t) E ) g [ n( t)] 4 ( u( t) E ) g ( u( t) E ) I( t) Na Na K K l l Ileak C g m ( u) 3 (1 )( u E ) g [ ] 4 ( u E ) g ( u E ) I( t) Na 0 Na K K l l a 1) dynamics of m are fast m( t) m ( u( t)) 0 2) dynamics of h and n are similar 1 h( t) a n( t) (t) (t) dh h h0( u) h ( u) d 0 ( u ) dn n n0( u) eff ( u) n ( u)
19 Neuronal Dynamics 3.1. Rection of Hodgkin-Huxley model I Na I K Ileak C g m ( u) 3 (1 )( u E ) g ( ) 4 ( u E ) g ( u E ) I( t) Na 0 Na K K l l a d 0 ( u ) eff ( u) F( u( t), ( t)) R I( t) d G ( u ( t ), ( t ))
20 NOW Exercise : separation of time scales C g m 3 h( u E ) g n 4 ( u E ) g ( u E ) I ( t) Na Na K K l l dm Exercises: m m0( u) m ( u) A: dx x c(t) - calculate x(t)! x(t) c(t) 0 1 t no! - hat if is small?? 1.5 homeork Exerc. 9h50-10h00 Next lecture: 10h15 dm m c( u) m f ( u) m B: -calculate m(t) if is small! - rece to 1 eq. t
21 Week 3 part 1 : Rection of the Hodgkin-Huxley Model 3.1 From Hodgkin-Huxley to 2D Biological Modeling of Neural Netorks - Overvie: From 4 to 2 dimensions - MathDetour 1: Exploiting similarities - MathDetour 2: Separation of time scales 3.2 Phase Plane Analysis - Role of nullclines Week 3 Recing detail: To-dimensional neuron models Wulfram Gerstner EPFL, Lausanne, Sitzerland 3.3 Analysis of a 2D Neuron Model - constant input vs pulse input - MathDetour 3: Stability of fixed points 3.4 TypeI and II Neuron Models next eek!
22 Neuronal Dynamics MathDetour 3.1: Separation of time scales To coupled differential equations dx Ex. 1-A x c() t 1 1 dx x c() t dy Exercise 1 (eek 3) 2 f ( y) g( x) Separation of time scales 1 2 Reced 1-dimensional system dy f ( y) g( c( t)) 2
23 Neuronal Dynamics MathDetour 3.2: Separation of time scales Linear differential equation dx x c() t 1 step
24 Neuronal Dynamics MathDetour 3.2 Separation of time scales To coupled differential equations 1 dx x c() t x 2 dc c f ( x) I( t) 1 2 c slo drive I
25 Neuronal Dynamics Rection of Hodgkin-Huxley model stimulus I I I Na K leak C g m 3 h( u E ) g n 4 ( u E ) g ( u E ) I ( t) Na Na K K l l dm m m0( u) m ( u) n 0 (u) n h (u) dh h h0( u) h ( u) h 0 (u) u m (u) u dn n n0( u) n ( u) dynamics of m is fast m( t) m ( u( t)) 0 Fast compared to hat?
26 Neuronal Dynamics MathDetour 3.2: Separation of time scales To coupled differential equations dx 1 x h( y) dy Exercise 1 (eek 3) 2 f ( y) g( x) Separation of time scales even more general x h( y) 1 2 Reced 1-dimensional system dy f ( y) g( h( y)) 2
27 Neuronal Dynamics Quiz 3.2. A- Separation of time scales: We start ith to equations dx x y I() t 1 2 dy [ ] If then the system can be reded to 2 y x A 1 2 dy y [ y I( t)] 2 A 2 Attention I(t) can move rapidly, therefore choice [1] not correct 2 1 [ ] If then the system can be reded to dx x x 2 A I() t 1 [ ] None of the above is correct.
28 Neuronal Dynamics Quiz 3.2-similar dynamics Exploiting similarities: A sufficient condition to replace to gating variables r,s by a single gating variable is [ ] Both r and s have the same time constant (as a function of u) [ ] Both r and s have the same activation function [ ] Both r and s have the same time constant (as a function of u) AND the same activation function [ ] Both r and s have the same time constant (as a function of u) AND activation functions that are identical after some additive rescaling [ ] Both r and s have the same time constant (as a function of u) AND activation functions that are identical after some multiplicative rescaling
29 Neuronal Dynamics 3.1. Rection to 2 dimensions 2-dimensional equation C f ( u( t), ( t)) I( t) d g ( u ( t ), ( t )) Enables graphical analysis! Phase plane analysis -Discussion of threshold - Constant input current vs pulse input -Type I and II - Repetitive firing
30 Week 3 part 1 : Rection of the Hodgkin-Huxley Model 3.1 From Hodgkin-Huxley to 2D - Overvie: From 4 to 2 dimensions Biological Modeling of Neural Netorks - MathDetour 1: Exploiting similarities - MathDetour 2: Separation of time scales 3.2 Phase Plane Analysis - Role of nullclines 3.3 Analysis of a 2D Neuron Model Week 3 Recing detail: To-dimensional neuron models Wulfram Gerstner EPFL, Lausanne, Sitzerland - constant input vs pulse input - MathDetour 3: Stability of fixed points 3.4 TypeI and II Neuron Models next eek!
31 Neuronal Dynamics 3.2. Reced Hodgkin-Huxley model I Na I K Ileak C g m ( u) 3 (1 )( u E ) g ( ) 4 ( u E ) g ( u E ) I( t) Na 0 Na K K l l a d 0 ( ( u) u) stimulus F( u, ) RI( t) d G( u, )
32 Neuronal Dynamics 3.2. Phase Plane Analysis/nullclines 2-dimensional equation stimulus F( u, ) RI( t) d G( u, ) u-nullcline Enables graphical analysis! -Discussion of threshold -Type I and II -nullcline
33 Neuronal Dynamics 3.2. FitzHugh-Nagumo Model F( u, ) RI ( t) 1 u u 3 RI () t 3 MathAnalysis, blackboard d G ( u, ) b b u 0 1 u-nullcline -nullcline
34 Neuronal Dynamics 3.2. flo arros F( u, ) RI( t) Stimulus I=0 d 0 d G( u, ) Consider change in small time step Flo on nullcline I(t)=0 u Flo in regions beteen nullclines Stable fixed point 0
35 Neuronal Dynamics Quiz 3.3 A. u-nullclines [ ] On the u-nullcline, arros are alays vertical [ ] On the u-nullcline, arros point alays vertically upard [ ] On the u-nullcline, arros are alays horizontal [ ] On the u-nullcline, arros point alays to the left [ ] On the u-nullcline, arros point alays to the right Take 1 minute, continue at 10:55 B. -Nullclines [ ] On the -nullcline, arros are alays vertical [ ] On the -nullcline, arros point alays vertically upard [ ] On the -nullcline, arros are alays horizontal [ ] On the -nullcline, arros point alays to the left [ ] On the -nullcline, arros point alays to the right [ ] On the -nullcline, arros can point in an arbitrary direction
36 Neuronal Dynamics 4.2. FitzHugh-Nagumo Model F( u, ) RI ( t) 1 u u 3 RI () t 3 d G ( u, ) b b u 0 1 change b1
37 Neuronal Dynamics 3.2. Nullclines of reced HH model d 0 stimulus F( u, ) RI( t) d u-nullcline G( u, ) 0 -nullcline Stable fixed point
38 Neuronal Dynamics 3.2. Phase Plane Analysis 2-dimensional equation stimulus F( u, ) RI( t) d G( u, ) Enables graphical analysis! Important role of Application to neuron models - nullclines - flo arros
39 Week 3 part 3: Analysis of a 2D neuron model 3.1 From Hodgkin-Huxley to 2D 3.2 Phase Plane Analysis - Role of nullcline 3.3 Analysis of a 2D Neuron Model - pulse input - constant input -MathDetour 3: Stability of fixed points 3.4 Type I and II Neuron Models (next eek)
40 Neuronal Dynamics 3.3. Analysis of a 2D neuron model 2-dimensional equation stimulus F( u, ) RI( t) 2 important input scenarios d G( u, ) - Pulse input Enables graphical analysis! - Constant input
41 Neuronal Dynamics D neuron model : Pulse input d u F ( u, ) RI d G ( u, ) pulse input
42 Neuronal Dynamics 3.3. FitzHugh-Nagumo Model : Pulse input 1 F( u, ) RI( t) u u 3 RI( t) 3 d 0 I(t)=0 d G ( u, ) b b u 0 1 u pulse input I(t) Pulse input: jump of voltage 0
43 Neuronal Dynamics 3.3. FitzHugh-Nagumo Model : Pulse input FN model ith b 0.9; b Pulse input: jump of voltage/initial condition Image: Neuronal Dynamics, Gerstner et al., Cambridge Univ. Press (2014)
44 Neuronal Dynamics 3.3. FitzHugh-Nagumo Model Pulse input: DONE! d 0 - jump of voltage - ne initial condition I(t)=0 - spike generation for large input pulses 2 important input scenarios constant input: u - graphics? - spikes? No - repetitive firing? 0
45 Neuronal Dynamics 3.3. FitzHugh-Nagumo Model: Constant input F( u, ) RI 0 d 0 -nullcline 1 u u 3 RI 3 0 d G ( u, ) b b u 0 1 Intersection point (fixed point) I(t)=I 0 u -moves -changes Stability 0 u-nullcline
46 NOW Exercise 2.1: Stability of Fixed Point in 2D u d 0 d u - calculate stability Exercises: - compare 2.1no! dx x I(t)=I 0 u 2.2 homeork Next lecture: 11:42 0
47 Week 3 part 3: Analysis of a 2D neuron model 3.1 From Hodgkin-Huxley to 2D 3.2 Phase Plane Analysis - Role of nullcline 3.3 Analysis of a 2D Neuron Model - pulse input - constant input -MathDetour 3: Stability of fixed points 3.4 Type I and II Neuron Models (next eek)
48 Neuronal Dynamics Detour 3.3 : Stability of fixed points. d b b u 0 1 d 0 -nullcline stable? I(t)=I 0 u F( u, ) RI 0 0 u-nullcline
49 Neuronal Dynamics 3.3 Detour. Stability of fixed points 2-dimensional equation stimulus F( u, ) RI 0 d G( u, ) Ho to determine stability of fixed point?
50 Neuronal Dynamics 3.3 Detour. Stability of fixed points stimulus d 0 au I 0 d cu I(t)=I 0 u 0 unstable saddle stable
51 Neuronal Dynamics 3.3 Detour. Stability of fixed points F( u, ) RI 0 d 0 -nullcline d G ( u, ) stable? zoom in: I(t)=I 0 u 0 unstable stable saddle Math derivation u-nullcline no
52 Neuronal Dynamics 3.3 Detour. Stability of fixed points F( u, ) RI 0 Fixed point at ( u, ) 0 0 At fixed point d G ( u, ) 0 F( u, ) RI G( u, ) 0 0 zoom in: x u u 0 y 0
53 Neuronal Dynamics 3.3 Detour. Stability of fixed points F( u, ) RI 0 Fixed point at ( u, ) 0 0 At fixed point d G ( u, ) 0 F( u, ) RI G( u, ) 0 0 zoom in: x u u 0 y 0 dx F x u F y dy G x G y u
54 Neuronal Dynamics 3.3 Detour. Stability of fixed points Linear matrix equation Search for solution To solution ith Eigenvalues, F u G F G F G u u
55 Neuronal Dynamics 3.3 Detour. Stability of fixed points Linear matrix equation Search for solution Stability requires: 0 and 0 To solution ith Eigenvalues, F u G F u and G 0 F G F G u u F G F G u u 0
56 Neuronal Dynamics 3.3 Detour. Stability of fixed points stimulus d 0 au I 0 d cu / I(t)=I 0 u 0 unstable saddle stable
57 Neuronal Dynamics 3.3 Detour. Stability of fixed points 2-dimensional equation F( u, ) RI 0 No Back: d G( u, ) Application to our neuron model Stability characterized by Eigenvalues of linearized equations
58 Neuronal Dynamics 3.3. FitzHugh-Nagumo Model: Constant input F( u, ) RI 0 d 0 -nullcline 1 u u 3 RI 3 0 d G ( u, ) b b u 0 1 Intersection point (fixed point) I(t)=I 0 u -moves -changes Stability 0 u-nullcline
59 Neuronal Dynamics 3.3. FitzHugh-Nagumo Model: Constant input F( u, ) RI 0 d 0 -nullcline 1 u u 3 RI 3 0 d G ( u, ) b b u 0 1 Intersection point (fixed point) I(t)=I 0 u -moves -changes Stability 0 u-nullcline
60 Neuronal Dynamics 3.3. FitzHugh-Nagumo Model : Constant input b 0.9; b 1.0; RI 2 FN model ith constant input: u-nullcline moves limit cycle f f-i curve Image: Neuronal Dynamics, Gerstner et al., CUP (2014) I0
61 Neuronal Dynamics Quiz 3.4. A. Short current pulses. In a 2-dimensional neuron model, the effect of a delta current pulse can be analyzed [ ] By moving the u-nullcline vertically upard [ ] By moving the -nullcline vertically upard [ ] As a potential change in the stability or number of the fixed point(s) [ ] As a ne initial condition [ ] By folloing the flo of arros in the appropriate phase plane diagram B. Constant current. In a 2-dimensional neuron model, the effect of a constant current can be analyzed [ ] By moving the u-nullcline vertically upard [ ] By moving the -nullcline vertically upard [ ] As a potential change in the stability or number of the fixed point(s) [ ] By folloing the flo of arros in the appropriate phase plane diagram
62 Computer exercise no Can e understand the dynamics of the 2D model? The END for today No: computer exercises Type I and type II models ramp input/ constant input f f-i curve f f-i curve I0 I0 I0
Neuronal Dynamics: Computational Neuroscience of Single Neurons
Week 4 part 1 : Redction of the Hodgkin-Hxley Model Neronal Dynamics: Comptational Neroscience of Single Nerons Week 4 Redcing detail: Two-dimensional neron models Wlfram Gerstner EPFL, Lasanne, Switzerland
More informationNeuronal Dynamics: Computational Neuroscience of Single Neurons
Week 4 part 5: Nonlinear Integrate-and-Fire Model 4.1 From Hodgkin-Huxley to 2D Neuronal Dynamics: Computational Neuroscience of Single Neurons Week 4 Recing detail: Two-dimensional neuron models Wulfram
More informationBiological Modeling of Neural Networks
Week 4 part 2: More Detail compartmental models Biological Modeling of Neural Networks Week 4 Reducing detail - Adding detail 4.2. Adding detail - apse -cable equat Wulfram Gerstner EPFL, Lausanne, Switzerland
More informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best Koch, Chapters 6, 8, 9 See also Hodgkin and Huxley, J. Physiol. 117:500-544 (1952. (the source Clay, J. Neurophysiol. 80:903-913
More informationSingle-Compartment Neural Models
Single-Compartment Neural Models BENG/BGGN 260 Neurodynamics University of California, San Diego Week 2 BENG/BGGN 260 Neurodynamics (UCSD) Single-Compartment Neural Models Week 2 1 / 18 Reading Materials
More informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best) Koch, Chapters 6, 8, 9 See also Clay, J. Neurophysiol. 80:903-913 (1998) (for a recent version of the HH squid axon model) Rothman
More informationNeuronal Dynamics: Computational Neuroscience of Single Neurons
Week 5 part 3a :Three definitions of rate code Neuronal Dynamics: Computational Neuroscience of Single Neurons Week 5 Variability and Noise: The question of the neural code Wulfram Gerstner EPFL, Lausanne,
More informationModeling Action Potentials in Cell Processes
Modeling Action Potentials in Cell Processes Chelsi Pinkett, Jackie Chism, Kenneth Anderson, Paul Klockenkemper, Christopher Smith, Quarail Hale Tennessee State University Action Potential Models Chelsi
More informationDecoding. How well can we learn what the stimulus is by looking at the neural responses?
Decoding How well can we learn what the stimulus is by looking at the neural responses? Two approaches: devise explicit algorithms for extracting a stimulus estimate directly quantify the relationship
More information3 Action Potentials - Brutal Approximations
Physics 172/278 - David Kleinfeld - Fall 2004; Revised Winter 2015 3 Action Potentials - Brutal Approximations The Hodgkin-Huxley equations for the behavior of the action potential in squid, and similar
More informationNeural Modeling and Computational Neuroscience. Claudio Gallicchio
Neural Modeling and Computational Neuroscience Claudio Gallicchio 1 Neuroscience modeling 2 Introduction to basic aspects of brain computation Introduction to neurophysiology Neural modeling: Elements
More information1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model
1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model Alan Hodgkin and Andrew Huxley developed the first quantitative model of the propagation of an electrical signal (the action potential)
More informationFrom neuronal oscillations to complexity
1/39 The Fourth International Workshop on Advanced Computation for Engineering Applications (ACEA 2008) MACIS 2 Al-Balqa Applied University, Salt, Jordan Corson Nathalie, Aziz Alaoui M.A. University of
More informationDynamical Systems in Neuroscience: Elementary Bifurcations
Dynamical Systems in Neuroscience: Elementary Bifurcations Foris Kuang May 2017 1 Contents 1 Introduction 3 2 Definitions 3 3 Hodgkin-Huxley Model 3 4 Morris-Lecar Model 4 5 Stability 5 5.1 Linear ODE..............................................
More information9 Generation of Action Potential Hodgkin-Huxley Model
9 Generation of Action Potential Hodgkin-Huxley Model (based on chapter 12, W.W. Lytton, Hodgkin-Huxley Model) 9.1 Passive and active membrane models In the previous lecture we have considered a passive
More informationNonlinear Dynamics of Neural Firing
Nonlinear Dynamics of Neural Firing BENG/BGGN 260 Neurodynamics University of California, San Diego Week 3 BENG/BGGN 260 Neurodynamics (UCSD) Nonlinear Dynamics of Neural Firing Week 3 1 / 16 Reading Materials
More informationChapter 2 The Hodgkin Huxley Theory of Neuronal Excitation
Chapter 2 The Hodgkin Huxley Theory of Neuronal Excitation Hodgkin and Huxley (1952) proposed the famous Hodgkin Huxley (hereinafter referred to as HH) equations which quantitatively describe the generation
More informationMATH 723 Mathematical Neuroscience Spring 2008 Instructor: Georgi Medvedev
MATH 723 Mathematical Neuroscience Spring 28 Instructor: Georgi Medvedev 2 Lecture 2. Approximate systems. 2. Reduction of the HH model to a 2D system. The original HH system consists of 4 differential
More informationTopics in Neurophysics
Topics in Neurophysics Alex Loebel, Martin Stemmler and Anderas Herz Exercise 2 Solution (1) The Hodgkin Huxley Model The goal of this exercise is to simulate the action potential according to the model
More informationNeural Spike Train Analysis 1: Introduction to Point Processes
SAMSI Summer 2015: CCNS Computational Neuroscience Summer School Neural Spike Train Analysis 1: Introduction to Point Processes Uri Eden BU Department of Mathematics and Statistics July 27, 2015 Spikes
More informationCompartmental Modelling
Modelling Neurons Computing and the Brain Compartmental Modelling Spring 2010 2 1 Equivalent Electrical Circuit A patch of membrane is equivalent to an electrical circuit This circuit can be described
More information/639 Final Solutions, Part a) Equating the electrochemical potentials of H + and X on outside and inside: = RT ln H in
580.439/639 Final Solutions, 2014 Question 1 Part a) Equating the electrochemical potentials of H + and X on outside and inside: RT ln H out + zf 0 + RT ln X out = RT ln H in F 60 + RT ln X in 60 mv =
More informationDynamical Systems in Neuroscience: The Geometry of Excitability and Bursting
Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting Eugene M. Izhikevich The MIT Press Cambridge, Massachusetts London, England Contents Preface xv 1 Introduction 1 1.1 Neurons
More informationSingle neuron models. L. Pezard Aix-Marseille University
Single neuron models L. Pezard Aix-Marseille University Biophysics Biological neuron Biophysics Ionic currents Passive properties Active properties Typology of models Compartmental models Differential
More informationBME 5742 Biosystems Modeling and Control
BME 5742 Biosystems Modeling and Control Hodgkin-Huxley Model for Nerve Cell Action Potential Part 1 Dr. Zvi Roth (FAU) 1 References Hoppensteadt-Peskin Ch. 3 for all the mathematics. Cooper s The Cell
More informationSingle-Cell and Mean Field Neural Models
1 Single-Cell and Mean Field Neural Models Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 The neuron
More informationLecture 11 : Simple Neuron Models. Dr Eileen Nugent
Lecture 11 : Simple Neuron Models Dr Eileen Nugent Reading List Nelson, Biological Physics, Chapter 12 Phillips, PBoC, Chapter 17 Gerstner, Neuronal Dynamics: from single neurons to networks and models
More informationMATH 3104: THE HODGKIN-HUXLEY EQUATIONS
MATH 3104: THE HODGKIN-HUXLEY EQUATIONS Parallel conductance model A/Prof Geoffrey Goodhill, Semester 1, 2009 So far we have modelled neuronal membranes by just one resistance (conductance) variable. We
More informationDynamical systems in neuroscience. Pacific Northwest Computational Neuroscience Connection October 1-2, 2010
Dynamical systems in neuroscience Pacific Northwest Computational Neuroscience Connection October 1-2, 2010 What do I mean by a dynamical system? Set of state variables Law that governs evolution of state
More informationLecture 10 : Neuronal Dynamics. Eileen Nugent
Lecture 10 : Neuronal Dynamics Eileen Nugent Origin of the Cells Resting Membrane Potential: Nernst Equation, Donnan Equilbrium Action Potentials in the Nervous System Equivalent Electrical Circuits and
More informationComputational Neuroscience. Session 4-2
Computational Neuroscience. Session 4-2 Dr. Marco A Roque Sol 06/21/2018 Two-Dimensional Two-Dimensional System In this section we will introduce methods of phase plane analysis of two-dimensional systems.
More informationMath 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction
Math 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction Junping Shi College of William and Mary November 8, 2018 Neuron Neurons Neurons are cells in the brain and other subsystems
More information2013 NSF-CMACS Workshop on Atrial Fibrillation
2013 NSF-CMACS Workshop on A Atrial Fibrillation Flavio H. Fenton School of Physics Georgia Institute of Technology, Atlanta, GA and Max Planck Institute for Dynamics and Self-organization, Goettingen,
More informationElectrophysiology of the neuron
School of Mathematical Sciences G4TNS Theoretical Neuroscience Electrophysiology of the neuron Electrophysiology is the study of ionic currents and electrical activity in cells and tissues. The work of
More informationComputational Neuroscience Summer School Neural Spike Train Analysis. An introduction to biophysical models (Part 2)
Computational Neuroscience Summer School Neural Spike Train Analysis Instructor: Mark Kramer Boston University An introduction to biophysical models (Part 2 SAMSI Short Course 2015 1 Goal: Model this,
More informationChapter 24 BIFURCATIONS
Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point Saddle-Node Bifurcation Diagram Codimension-1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses
More informationBasic elements of neuroelectronics -- membranes -- ion channels -- wiring. Elementary neuron models -- conductance based -- modelers alternatives
Computing in carbon Basic elements of neuroelectronics -- membranes -- ion channels -- wiring Elementary neuron models -- conductance based -- modelers alternatives Wiring neurons together -- synapses
More informationSignal processing in nervous system - Hodgkin-Huxley model
Signal processing in nervous system - Hodgkin-Huxley model Ulrike Haase 19.06.2007 Seminar "Gute Ideen in der theoretischen Biologie / Systembiologie" Signal processing in nervous system Nerve cell and
More informationDynamical Systems for Biology - Excitability
Dynamical Systems for Biology - Excitability J. P. Keener Mathematics Department Dynamical Systems for Biology p.1/25 Examples of Excitable Media B-Z reagent CICR (Calcium Induced Calcium Release) Nerve
More informationIntroduction and the Hodgkin-Huxley Model
1 Introduction and the Hodgkin-Huxley Model Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 Reference:
More informationBiological Modeling of Neural Networks:
Week 14 Dynamics and Plasticity 14.1 Reservoir computing - Review:Random Networks - Computing with rich dynamics Biological Modeling of Neural Networks: 14.2 Random Networks - stationary state - chaos
More informationBIOELECTRIC PHENOMENA
Chapter 11 BIOELECTRIC PHENOMENA 11.3 NEURONS 11.3.1 Membrane Potentials Resting Potential by separation of charge due to the selective permeability of the membrane to ions From C v= Q, where v=60mv and
More informationFRTF01 L8 Electrophysiology
FRTF01 L8 Electrophysiology Lecture Electrophysiology in general Recap: Linear Time Invariant systems (LTI) Examples of 1 and 2-dimensional systems Stability analysis The need for non-linear descriptions
More informationMemory and hypoellipticity in neuronal models
Memory and hypoellipticity in neuronal models S. Ditlevsen R. Höpfner E. Löcherbach M. Thieullen Banff, 2017 What this talk is about : What is the effect of memory in probabilistic models for neurons?
More informationMathematical Foundations of Neuroscience - Lecture 3. Electrophysiology of neurons - continued
Mathematical Foundations of Neuroscience - Lecture 3. Electrophysiology of neurons - continued Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland
More informationDeconstructing Actual Neurons
1 Deconstructing Actual Neurons Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 Reference: The many ionic
More informationBasic elements of neuroelectronics -- membranes -- ion channels -- wiring
Computing in carbon Basic elements of neuroelectronics -- membranes -- ion channels -- wiring Elementary neuron models -- conductance based -- modelers alternatives Wires -- signal propagation -- processing
More information6.3.4 Action potential
I ion C m C m dφ dt Figure 6.8: Electrical circuit model of the cell membrane. Normally, cells are net negative inside the cell which results in a non-zero resting membrane potential. The membrane potential
More informationBursting Oscillations of Neurons and Synchronization
Bursting Oscillations of Neurons and Synchronization Milan Stork Applied Electronics and Telecommunications, Faculty of Electrical Engineering/RICE University of West Bohemia, CZ Univerzitni 8, 3064 Plzen
More informationElectronics 101 Solving a differential equation Incorporating space. Numerical Methods. Accuracy, stability, speed. Robert A.
Numerical Methods Accuracy, stability, speed Robert A. McDougal Yale School of Medicine 21 June 2016 Hodgkin and Huxley: squid giant axon experiments Top: Alan Lloyd Hodgkin; Bottom: Andrew Fielding Huxley.
More informationSTUDENT PAPER. Santiago Santana University of Illinois, Urbana-Champaign Blue Waters Education Program 736 S. Lombard Oak Park IL, 60304
STUDENT PAPER Differences between Stochastic and Deterministic Modeling in Real World Systems using the Action Potential of Nerves. Santiago Santana University of Illinois, Urbana-Champaign Blue Waters
More informationActivity Driven Adaptive Stochastic. Resonance. Gregor Wenning and Klaus Obermayer. Technical University of Berlin.
Activity Driven Adaptive Stochastic Resonance Gregor Wenning and Klaus Obermayer Department of Electrical Engineering and Computer Science Technical University of Berlin Franklinstr. 8/9, 187 Berlin fgrewe,obyg@cs.tu-berlin.de
More informationThe Spike Response Model: A Framework to Predict Neuronal Spike Trains
The Spike Response Model: A Framework to Predict Neuronal Spike Trains Renaud Jolivet, Timothy J. Lewis 2, and Wulfram Gerstner Laboratory of Computational Neuroscience, Swiss Federal Institute of Technology
More informationIntegral Equation and PDE Neuronal Models:
Integral Equation and PDE Neuronal Models: 1952-2015 W. C. Troy Integral Equation and PDE Neuronal Models: 1952-2015 p.1/51 Part 1. 1D Single Neuron Models. (i) Hodgkin-Huxley (1952), Fitzhugh-Nagumo (1961)
More informationProblem Set Number 02, j/2.036j MIT (Fall 2018)
Problem Set Number 0, 18.385j/.036j MIT (Fall 018) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139) September 6, 018 Due October 4, 018. Turn it in (by 3PM) at the Math. Problem Set
More informationAction Potential Propagation
Action Potential Propagation 2 Action Potential is a transient alteration of transmembrane voltage (or membrane potential) across an excitable membrane generated by the activity of voltage-gated ion channels.
More information1 Introduction and neurophysiology
Dynamics of Continuous, Discrete and Impulsive Systems Series B: Algorithms and Applications 16 (2009) 535-549 Copyright c 2009 Watam Press http://www.watam.org ASYMPTOTIC DYNAMICS OF THE SLOW-FAST HINDMARSH-ROSE
More informationPresented by Sarah Hedayat. Supervised by Pr.Cappy and Dr.Hoel
1 Presented by Sarah Hedayat Supervised by Pr.Cappy and Dr.Hoel Outline 2 Project objectives Key elements Membrane models As simple as possible Phase plane analysis Review of important Concepts Conclusion
More informationNeural Coding: Integrate-and-Fire Models of Single and Multi-Neuron Responses
Neural Coding: Integrate-and-Fire Models of Single and Multi-Neuron Responses Jonathan Pillow HHMI and NYU http://www.cns.nyu.edu/~pillow Oct 5, Course lecture: Computational Modeling of Neuronal Systems
More informationBalance of Electric and Diffusion Forces
Balance of Electric and Diffusion Forces Ions flow into and out of the neuron under the forces of electricity and concentration gradients (diffusion). The net result is a electric potential difference
More informationCSE/NB 528 Final Lecture: All Good Things Must. CSE/NB 528: Final Lecture
CSE/NB 528 Final Lecture: All Good Things Must 1 Course Summary Where have we been? Course Highlights Where do we go from here? Challenges and Open Problems Further Reading 2 What is the neural code? What
More informationResting Distribution of Ions in Mammalian Neurons. Outside Inside (mm) E ion Permab. K Na Cl
Resting Distribution of Ions in Mammalian Neurons Outside Inside (mm) E ion Permab. K + 5 100-81 1.0 150 15 +62 0.04 Cl - 100 10-62 0.045 V m = -60 mv V m approaches the Equilibrium Potential of the most
More informationSample HHSim Exercises
I. Equilibrium Potential II. Membrane Potential III. The Action Potential IV. The Fast Sodium Channel V. The Delayed Rectifier VI. Voltage-Gated Channel Parameters Part I: Equilibrium Potential Sample
More informationNeuronal Dynamics: Computational Neuroscience of Single Neurons
Nonlinear Integrate-and-Fire Model Neronal Dynamics: Comptational Neroscience of Single Nerons Nonlinear Integrate-and-fire (NLIF) - Definition - qadratic and expon. IF - Extracting NLIF model from data
More informationModeling cells in tissue by connecting. electrodiffusion and Hodgkin Huxley-theory.
Modeling cells in tissue by connecting electrodiffusion and Hodgkin Huxley-theory. Brigham Young University Department of Mathematics Viktoria R.T. Hsu University of Utah Department of Mathematics Outline
More informationCellular Electrophysiology. Cardiac Electrophysiology
Part 1: Resting and Action Potentials Cardiac Electrophysiology Theory Simulation Experiment Scale The membrane: structure, channels and gates The cell: resting potential, whole cell currents, cardiac
More informationDynamics and complexity of Hindmarsh-Rose neuronal systems
Dynamics and complexity of Hindmarsh-Rose neuronal systems Nathalie Corson and M.A. Aziz-Alaoui Laboratoire de Mathématiques Appliquées du Havre, 5 rue Philippe Lebon, 766 Le Havre, FRANCE nathalie.corson@univ-lehavre.fr
More informationIntegrating mathematics and modelling into life science programs
MATRIC MODELING SEMINAR August 9, 2016 Integrating mathematics and modelling into life science programs André Heck Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents 1. What is
More informationSimple models of neurons!! Lecture 4!
Simple models of neurons!! Lecture 4! 1! Recap: Phase Plane Analysis! 2! FitzHugh-Nagumo Model! Membrane potential K activation variable Notes from Zillmer, INFN 3! FitzHugh Nagumo Model Demo! 4! Phase
More informationACTION POTENTIAL. Dr. Ayisha Qureshi Professor MBBS, MPhil
ACTION POTENTIAL Dr. Ayisha Qureshi Professor MBBS, MPhil DEFINITIONS: Stimulus: A stimulus is an external force or event which when applied to an excitable tissue produces a characteristic response. Subthreshold
More informationIntroduction to Neural Networks U. Minn. Psy 5038 Spring, 1999 Daniel Kersten. Lecture 2a. The Neuron - overview of structure. From Anderson (1995)
Introduction to Neural Networks U. Minn. Psy 5038 Spring, 1999 Daniel Kersten Lecture 2a The Neuron - overview of structure From Anderson (1995) 2 Lect_2a_Mathematica.nb Basic Structure Information flow:
More informationLIMIT CYCLE OSCILLATORS
MCB 137 EXCITABLE & OSCILLATORY SYSTEMS WINTER 2008 LIMIT CYCLE OSCILLATORS The Fitzhugh-Nagumo Equations The best example of an excitable phenomenon is the firing of a nerve: according to the Hodgkin
More informationMODEL NEURONS: FROM HODGKIN-HUXLEY TO HOPFIELD. L.F. Abbott and Thomas B. Kepler. Physics and Biology Departments. and. Center for Complex Systems
MODEL NEURONS: FROM HODGKIN-HUXLEY TO HOPFIELD L.F. Abbott and Thomas B. Kepler Physics and Biology Departments and Center for Complex Systems Brandeis University Waltham, MA 02254 Abstract A new technique
More information9 Generation of Action Potential Hodgkin-Huxley Model
9 Generation of Action Potential Hodgkin-Huxley Model (based on chapter 2, W.W. Lytton, Hodgkin-Huxley Model) 9. Passive and active membrane models In the previous lecture we have considered a passive
More informationWave Pinning, Actin Waves, and LPA
Wave Pinning, Actin Waves, and LPA MCB 2012 William R. Holmes Intercellular Waves Weiner et al., 2007, PLoS Biology Dynamic Hem1 waves in neutrophils Questions? How can such waves / pulses form? What molecular
More informationNeurons. The Molecular Basis of their Electrical Excitability
Neurons The Molecular Basis of their Electrical Excitability Viva La Complexity! Consider, The human brain contains >10 11 neurons! Each neuron makes 10 3 (average) synaptic contacts on up to 10 3 other
More informationDynamical phase transitions in periodically driven model neurons
Dynamical phase transitions in periodically driven model neurons Jan R. Engelbrecht 1 and Renato Mirollo 2 1 Department of Physics, Boston College, Chestnut Hill, Massachusetts 02467, USA 2 Department
More informationEffects of synaptic conductance on the voltage distribution and firing rate of spiking neurons
PHYSICAL REVIEW E 69, 051918 (2004) Effects of synaptic conductance on the voltage distribution and firing rate of spiking neurons Magnus J. E. Richardson* Laboratory of Computational Neuroscience, Brain
More informationBo Deng University of Nebraska-Lincoln UNL Math Biology Seminar
Mathematical Model of Neuron Bo Deng University of Nebraska-Lincoln UNL Math Biology Seminar 09-10-2015 Review -- One Basic Circuit By Kirchhoff's Current Law 0 = I C + I R + I L I ext By Kirchhoff s Voltage
More informationNeuron. Detector Model. Understanding Neural Components in Detector Model. Detector vs. Computer. Detector. Neuron. output. axon
Neuron Detector Model 1 The detector model. 2 Biological properties of the neuron. 3 The computational unit. Each neuron is detecting some set of conditions (e.g., smoke detector). Representation is what
More informationTransmission of Nerve Impulses (see Fig , p. 403)
How a nerve impulse works Transmission of Nerve Impulses (see Fig. 12.13, p. 403) 1. At Rest (Polarization) outside of neuron is positively charged compared to inside (sodium ions outside, chloride and
More information2182. Hopf bifurcation analysis and control of three-dimensional Prescott neuron model
2182. Hopf bifurcation analysis and control of three-dimensional Prescott neuron model Chunhua Yuan 1, Jiang Wang 2 School of Electrical Engineering and Automation, Tianjin University, Tianjin, China 2
More informationMathematical analysis of a 3D model of cellular electrophysiology
Mathematical analysis of a 3D model of cellular electrophysiology Hiroshi Matano (Univ. of Tokyo) joint work with Yoichiro Mori (Univ. of Minnesota) Seoul-Tokyo Conference on Elliptic and Parabolic PDEs
More informationOn Parameter Estimation for Neuron Models
On Parameter Estimation for Neuron Models Abhijit Biswas Department of Mathematics Temple University November 30th, 2017 Abhijit Biswas (Temple University) On Parameter Estimation for Neuron Models November
More information3 Detector vs. Computer
1 Neurons 1. The detector model. Also keep in mind this material gets elaborated w/the simulations, and the earliest material is often hardest for those w/primarily psych background. 2. Biological properties
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 203 Outline
More information3.3 Simulating action potentials
6 THE HODGKIN HUXLEY MODEL OF THE ACTION POTENTIAL Fig. 3.1 Voltage dependence of rate coefficients and limiting values and time constants for the Hodgkin Huxley gating variables. (a) Graphs of forward
More informationarxiv: v1 [q-bio.nc] 13 Feb 2018
Gain control with A-type potassium current: I A as a switch between divisive and subtractive inhibition Joshua H Goldwyn 1*, Bradley R Slabe 2, Joseph B Travers 3, David Terman 2 arxiv:182.4794v1 [q-bio.nc]
More informationNonlinear Observer Design and Synchronization Analysis for Classical Models of Neural Oscillators
Nonlinear Observer Design and Synchronization Analysis for Classical Models of Neural Oscillators Ranjeetha Bharath and Jean-Jacques Slotine Massachusetts Institute of Technology ABSTRACT This work explores
More informationΝευροφυσιολογία και Αισθήσεις
Biomedical Imaging & Applied Optics University of Cyprus Νευροφυσιολογία και Αισθήσεις Διάλεξη 5 Μοντέλο Hodgkin-Huxley (Hodgkin-Huxley Model) Response to Current Injection 2 Hodgin & Huxley Sir Alan Lloyd
More informationLokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt
Periodic solutions A periodic solution is a solution (x(t), y(t)) of dx = f(x, y) dt dy = g(x, y) dt such that x(t + T ) = x(t) and y(t + T ) = y(t) for any t, where T is a fixed number which is a period
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 2013 Outline
More informationChimera states in networks of biological neurons and coupled damped pendulums
in neural models in networks of pendulum-like elements in networks of biological neurons and coupled damped pendulums J. Hizanidis 1, V. Kanas 2, A. Bezerianos 3, and T. Bountis 4 1 National Center for
More informationNeuroscience applications: isochrons and isostables. Alexandre Mauroy (joint work with I. Mezic)
Neuroscience applications: isochrons and isostables Alexandre Mauroy (joint work with I. Mezic) Outline Isochrons and phase reduction of neurons Koopman operator and isochrons Isostables of excitable systems
More informationA FINITE STATE AUTOMATON MODEL FOR MULTI-NEURON SIMULATIONS
A FINITE STATE AUTOMATON MODEL FOR MULTI-NEURON SIMULATIONS Maria Schilstra, Alistair Rust, Rod Adams and Hamid Bolouri Science and Technology Research Centre, University of Hertfordshire, UK Department
More informationNerve Signal Conduction. Resting Potential Action Potential Conduction of Action Potentials
Nerve Signal Conduction Resting Potential Action Potential Conduction of Action Potentials Resting Potential Resting neurons are always prepared to send a nerve signal. Neuron possesses potential energy
More informationAn Introductory Course in Computational Neuroscience
An Introductory Course in Computational Neuroscience Contents Series Foreword Acknowledgments Preface 1 Preliminary Material 1.1. Introduction 1.1.1 The Cell, the Circuit, and the Brain 1.1.2 Physics of
More informationResting membrane potential,
Resting membrane potential Inside of each cell is negative as compared with outer surface: negative resting membrane potential (between -30 and -90 mv) Examination with microelectrode (Filled with KCl
More informationPhys 102 Lecture 9 RC circuits
Phys 102 Lecture 9 RC circuits 1 Recall from last time... We solved various circuits with resistors and batteries (also capacitors and batteries) ε R 1 R 2 R 3 R 1 ε 1 ε 2 R 3 What about circuits that
More informationLecture Notes 8C120 Inleiding Meten en Modelleren. Cellular electrophysiology: modeling and simulation. Nico Kuijpers
Lecture Notes 8C2 Inleiding Meten en Modelleren Cellular electrophysiology: modeling and simulation Nico Kuijpers nico.kuijpers@bf.unimaas.nl February 9, 2 2 8C2 Inleiding Meten en Modelleren Extracellular
More information