MATH 723 Mathematical Neuroscience Spring 2008 Instructor: Georgi Medvedev
|
|
- Lawrence Morton
- 5 years ago
- Views:
Transcription
1 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 equations. In the present lecture, we introduce several approximations which allow to reduce the HH model to a system of 2 differential equations. The analysis of the approximate 2D model provides an insight into the mechanism for spike generation in the HH model. The following oservations lead to the desired reduction: Step. Note that τ m (V ) τ h (V ),τ n (V ) (see Figure.2),i.e., the dynamics of the activation of the sodium current has the fastest time scale. Therefore, to leading order, we approximate m(t) y its steady state value: m m (V ). The resulting system has the following form: V = g Na m 3 (V )h(v E Na ) g K n 4 (V E K ) g L (V E L ) + I, (2.) ṅ = n (V ) n, (2.2) τ n (V ) ḣ = h (V ) h. (2.3) τ h (V ) Step 2. The following oservation is due to Krinsky and Kokoz [KK]. The dynamics of h is similar to that of n. The numerics in Figure 2. show after a short period of time necessary for transients to settle down, h + n. More precisely,.n + h.89 (see Figure 2.). By approximating h l(n), l(n) =.89.n, we arrive at the following system of two differential equations: V = g Na m 3 (V )l(n)(v E Na ) g K n 4 (V E K ) g L (V E L ) + I, (2.4) ṅ = n (V ) n, (2.5) τ n (V ) Step 3. Finally, it is convenient to approximate steady-state functions m (v) and n (v) and time constant τ n (v) using hyperolic functions [BRS]: { ( )} V θs s (V ) = + exp, s = m,n, (2.6) σ s ( ) V θn τ n (V ) = τ n cosh, (2.7) 2σ n The original functions and their approximations (2.6) and (2.7) are shown in Figure 2.2a. The parameter values for the reduced model (2.4) and (2.5) are summarized in
2 V (mv) V(t) n(t) + h(t) h h =.89.n.5 h(t).2 a n(t) t (ms) n Figure 2.: Plots in a and show that variales n and h are approximately linearly dependent. Tale 2. g Na 2mS/cm 2 E Na 2mV θ m 25mV σ m 9.5 g K 36mS/cm 2 E K 2mV θ n.3mv σ n 9.5 τ n 2 g L.3mS/cm 2 E L.6mV C µf/cm The rescaled system. Variales v and n in the reduced HH model have different models of magnitude. During an AP, voltage varies from 2mV to aout mv, while the values of the gating variale range from.4 to. To analyze (2.4) and (2.5), it is convenient to rescale the dynamical variales so that they have the same orders of magnitude: V = k V v, k V = mv, τ = k t t, k t = ms. (2.8) We also rescale the maximal conductances and the reversal potentials: g s = g s mscm 2 and E s = Ẽs k v, s = {m,n}. Thus, we otain the following nondimensional system: ǫ v = g Na m 3 (v)h(v ẼNa) g K n 4 (v ẼK) g L (v ẼL) + Ĩ,, (2.9) ṅ = ñ (v) n. (2.) τ n (v) where s (v) = s (k V v), s = {m,n} and τ n (v) = τ n (k V v). The derivative of v is multiplied y a small nondimensional parameter ǫ = 2. This indicates that that the dynamics of v and n have different timescales. Equations (2.) and (2.) imply that typically: v = O ( ǫ ) while ṅ = O(), i.e., v changes 2
3 .8 m (v).8 n (v) τ n (v).2.2 a 5 V (mv) 5 V (mv) 5 5 V (mv) Figure 2.2: a. Steady-state functions m (v) and n (c), and and time constant τ n (v) (solid line) and their approximations y the hyperolic functions (dashed line).. Action potential generated y the 2D approximation of the HH model y applying a step current pulse. 3
4 a a Figure 2.3: The rescaled system: the phase plane and time series: a) excitale regime and ) oscillating regime. much faster than n. For this reason, v is called a fast variale. The v time series in Figure 2.3 show alternating intervals of fast and slow dynamics. Systems of ODEs with small parameters multiplying timederivatives in some of the equations are called singularly pertured or fast-slow systems. In the following susection, we descrie a geometric method for analyzing a 2D singularly pertured system. 2.3 The FitzHugh-Nagumo model: slow-decomposition decomposition Below we outline some ideas used for analyzing singularly pertured systems of ODEs. To make our discussion more explicit, we will use a specific example of a 2D singularly pertured system, the FitzHugh- Nagumo model: ǫ v = f(v) n + I, (2.) ṅ = v γn, (2.2) where v and n are the dynamical variales, and I,γ >, < ǫ are parameters. A smooth function f(v) is qualitatively cuic, i.e. it has two extrema and lim v ± f(v) = ±. A small parameter ǫ > is called a 4
5 5 4 N v N n L v 2 (v o, n o ) n (v * (n o ), n o ) 5 a c d n 5 a v 5 Figure 2.4: The FitzHugh-Nagumo model: a. The phase plane,. Solution plots. relaxation parameter. Along with (2.)-(2.2), we shall consider the following system of equations: dv dτ dn dτ = f(v) n + I, (2.3) = ǫ(v γn), (2.4) where τ = ǫ t is called the fast time. By formally taking the limit as ǫ in (2.3) and (2.4) we arrive at the fast susystem: dṽ dτ dñ dτ = f(ṽ) ñ + I, (2.5) =. (2.6) Let (v,n ) e the initial condition for (2.5), (2.6), as shown in Fig 2.3. For almost all initial conditions (v,n ), v elongs to the domain of attraction of one of the stale fixed points of the D ODE: dṽ dτ = f(ṽ) n + I. Denote this fixed point y v (n ). The trajectory of the fast system (2.5), (2.6) starting at (v,n ) tends to (v (n ),n ) as τ along the line n = n (see Figure 2.4). The vector field generated y (2.5),(2.6) is called the fast vector field. The theory for singularly pertured systems implies that for sufficiently small ǫ >, the solutions of the fast system (2.5) and (2.6) approximate those of the original system (2.3) and (2.4) with the (same initial condition) up to the time when the trajectory of the fast system enters some O(ǫ) neighorhood of (v (n ),n ). To study the dynamics of (2.3) and (2.4) after that moment of time, we consider original system (2.), and (2.2). Again, y formally taking the limit as ǫ, from (2.3) and (2.4) we derive the slow system: n = f(v) + I, (2.7) ṅ = v γn, (2.8) 5
6 The dynamics of the slow system is constrained y the two outer (stale) ranches of the slow curve n = f(v) + I. The set of points with coordinates satisfying (2.7) is also called the v nullcline: N = {(v,n) : n = f(v) + I} (see Figure 2.4). The motion along N is determined y the equation (2.8). By differentiating oth sides of (2.7) with respect to time and y comining the result with (2.8), we otain: By plugging in (2.7) into (2.9) we otain dv dt f (v) dv dt = v γn, (2.9) v γ(f(v) + I) = f, (2.2) (v) Equation (2.2) is valid outside of some neighorhoods of the points of extrema of f(v) where f (v) =. According to (2.2) on the slow manifold trajectory with initial condition, as shown in Figure 2.4, will move down along the left hand outer ranch of N v until it reaches an O(ǫ) neighorhood of a. After that the point will jump to the opposite stale ranch of N v under the action of the fast vector field. The description of the slow dynamics along this ranch is completely analogous to that for the left outer ranch. Therefore, after some finite time the phase trajectory remains in the neighorhood of the closed curve L composed of the segments of the slow manifold joined y horizontal segments (see Figure 2.4a). L is called a relaxation limit cycle. The alternation of the slow motions and fast jumps generates relaxation oscillations (see Figure 2.4). The aove discussion shows that for a 2D singularly pertured systems such as (2.3) and (2.4), the phase plane is to a large extent is determined y the geometries of the v and n nullclines: n = f(v) + I and n = γ v. In the FHN model, parameter I plays the same role as the applied current in the HH model. Let us consider how the variation of I affects the qualitative dynamics of the FHN model. Note that changing I results in the translation of the voltage nullcline in the vertical direction. In particular, y decreasing I we can achieve that the point of the intersection of the two nullclines ( v, n) lies on the left outer ranch of L (see Figure 2.4). Note that at (v,n) = ( v, n) the right-hand side of (2.) and (2.2) are zero. Therefore, ( v, n) is a steady state or a fixed point. For the FHN, if the fixed point is located on the stale ranch of the slow curve it attracts the trajectories following it. Therefore, after a finite time the dynamics of the system approaches the steady state ( v, n). The phase plane discussion for the FHN system gives us a first insight into the mechanism of the excitation in the HH system. By changing the applied current in the HH system, we change the character of the fixed (it ceases to e attracting) and create an attracting limit cycle similar to L in the present example. 2.4 Two types of neuronal excitaility We conclude this lecture with an introductory discussion of the two major mechanisms for neuronal excitaility (cf. [RE]). For this, we will use the phase plane techniques outlined in the present section. In 6
7 a c d Figure 2.5: Phase portraits plotted for the ML model with different parameter values, illustrating two mechanisms of excitaility: (a,) Type II, (c,d) Type I. the previous susection, we showed that y changing the value of the applied current, I, one can affect the position of the fixed point in the FH system. The fixed point moves from the left ranch to the middle ranch of the cuic nullcline under the increase in I. For low values of I, the fixed point lies on the left ranch and it locks the motion of the trajectories along the slow manifold. Shortly after the fixed point passes the lower fold of N v under the variation of I, it opens the way for the trajectories starting near the left ranch to jump toward the right ranch and gives rise to a limit cycle. This is an example of an excitale system, i.e., a system where a significant deviation from a stale steady state can e induced y a relatively small change of a control parameter. A similar mechanism is responsile for generation of the AP in the reduced HH system (Fig. 2.5a,). Next we discuss a different excitaility mechanism illustrated in Fig. 2.5c,d. As in the previous case, the system shown in Fig. 2.6c is in an excitale state: the fixed point is lying on the left ranch of the cuic nullcline. However, in contrast to the case discussed aove, under the increase of applied current the fixed point collides with another one and disappears (Fig. 2.5d). As a result the limit cycle is orn. The representative time series corresponding the two scenarios of transition to spiking are shown in Fig. 2.7 and 2.7. Note that in Fig. 2.6, the frequency of oscillations just near the transition is extremely low. This turns out to e a distinctive feature of a road class of models, Type I models. The FHN model and the 7
8 a c Figure 2.6: Phase portraits and the time series illustrating transition to spiking in Type I models. reduced HH model elong to the Type II class of models. The frequency of oscillations is ounded from zero near the transition to spiking in Type II models. Explanation of these effects e given later after a rief introduction to elementary ifurcation theory in Lecture 4. Homework project II. a. Plot the phase portrait in the ML model. Determine whether it is a type I or type II model. Change the value(s) of some parameter(s) in this model to make it of a different type.. For the ML model in the Type II regime, use the phase plane analysis (along the lines of Section 2.3) to descrie the mechanism of the AP generation under electrical stimulation. Derive the formula for the period of the relaxation limit cycle. Biliography BRS R.J. Butera, J. Rinzel, J.C. Smith, Models of Respiratory rhythm generation in the Pre-Botzinger 8
9 a c Figure 2.7: Phase portraits and the time series illustrating transition to spiking in Type II models. 9
10 Complex: I. Bursting pacemaker neurons, J. of Neurophys. 8: , 999) KK VI Krinsky and YM Kokoz, Analysis of equations of excitale memranes I. Reduction of the Hodgkin- Huxley equations to a second order system. Biofizika, 8, 56 5, 973. MKKR E.F. Mishchenko, Yu.S. Kolesov, A.Yu. Kolesov, and N.Kh. Rozov, Asymptotic Methods in Singularly pertured Systems, Consultants Bureau, New York, 994. RE J. Rinzel and G.B. Ermentrout, Analysis of neural excitaility and oscillations, in C. Koch and I. Segev, eds Methods in Neuronal Modeling, MIT Press, Camridge, MA, 989.
Modeling 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 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 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 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 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 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 informationarxiv:nlin/ v1 [nlin.cd] 13 May 2006
On transition to bursting via deterministic chaos arxiv:nlin/0605030v1 [nlin.cd] 13 May 2006 Georgi S. Medvedev September 23, 2018 Abstract We study statistical properties of the irregular bursting arising
More informationBiological Modeling of Neural Networks
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
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 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 informationEntrainment and Chaos in the Hodgkin-Huxley Oscillator
Entrainment and Chaos in the Hodgkin-Huxley Oscillator Kevin K. Lin http://www.cims.nyu.edu/ klin Courant Institute, New York University Mostly Biomath - 2005.4.5 p.1/42 Overview (1) Goal: Show that the
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 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 informationLinearization of F-I Curves by Adaptation
LETTER Communicated by Laurence Abbott Linearization of F-I Curves by Adaptation Bard Ermentrout Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260, U.S.A. We show that negative
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 informationCANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION
CANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION WARREN WECKESSER Department of Mathematics Colgate University Hamilton, NY 3346 E-mail: wweckesser@mail.colgate.edu Cartwright and Littlewood discovered
More informationCHAPTER 1. Bifurcations in the Fast Dynamics of Neurons: Implications for Bursting
CHAPTER 1 Bifurcations in the Fast Dynamics of Neurons: Implications for Bursting John Guckenheimer Mathematics Department, Cornell University Ithaca, NY 14853 Joseph H. Tien Center for Applied Mathematics,
More informationThe selection of mixed-mode oscillations in a Hodgkin-Huxley model with multiple timescales
CHAOS 18, 015105 2008 The selection of mixed-mode oscillations in a Hodgkin-Huxley model with multiple timescales Jonathan Rubin Department of Mathematics, University of Pittsburgh, Pittsburgh, Pennsylvania
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 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 informationHybrid Integrate-and-Fire Model of a Bursting Neuron
LETTER Communicated by John Rinzel Hybrid Integrate-and-Fire Model of a Bursting Neuron Barbara J. Breen bbreen@ece.gatech.edu William C. Gerken wgerken@ece.gatech.edu Robert J. Butera, Jr. rbutera@ece.gatech.edu
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 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 informationReduction of Conductance Based Models with Slow Synapses to Neural Nets
Reduction of Conductance Based Models with Slow Synapses to Neural Nets Bard Ermentrout October 1, 28 Abstract The method of averaging and a detailed bifurcation calculation are used to reduce a system
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 informationThe Phase Response Curve of Reciprocally Inhibitory Model Neurons Exhibiting Anti-Phase Rhythms
The Phase Response Curve of Reciprocally Inhibitory Model Neurons Exhibiting Anti-Phase Rhythms Jiawei Zhang Timothy J. Lewis Department of Mathematics, University of California, Davis Davis, CA 9566,
More information1 R.V k V k 1 / I.k/ here; we ll stimulate the action potential another way.) Note that this further simplifies to. m 3 k h k.
1. The goal of this problem is to simulate a propagating action potential for the Hodgkin-Huxley model and to determine the propagation speed. From the class notes, the discrete version (i.e., after breaking
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 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 informationDynamical modelling of systems of coupled oscillators
Dynamical modelling of systems of coupled oscillators Mathematical Neuroscience Network Training Workshop Edinburgh Peter Ashwin University of Exeter 22nd March 2009 Peter Ashwin (University of Exeter)
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 informationStochastic Oscillator Death in Globally Coupled Neural Systems
Journal of the Korean Physical Society, Vol. 52, No. 6, June 2008, pp. 19131917 Stochastic Oscillator Death in Globally Coupled Neural Systems Woochang Lim and Sang-Yoon Kim y Department of Physics, Kangwon
More informationDetermining the contributions of divisive and subtractive feedback in the Hodgkin-Huxley model
JComputNeurosci(2014)37:403 415 DOI 10.1007/s10827-014-0511-y Determining the contributions of divisive and subtractive feedback in the Hodgkin-Huxley model Sevgi Şengül & Robert Clewley & Richard Bertram
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 informationFinal Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger. Project I: Predator-Prey Equations
Final Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger Project I: Predator-Prey Equations The Lotka-Volterra Predator-Prey Model is given by: du dv = αu βuv = ρβuv
More informationFrequency Adaptation and Bursting
BioE332A Lab 3, 2010 1 Lab 3 January 5, 2010 Frequency Adaptation and Bursting In the last lab, we explored spiking due to sodium channels. In this lab, we explore adaptation and bursting due to potassium
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 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 informationNumerical Simulation of Bistability between Regular Bursting and Chaotic Spiking in a Mathematical Model of Snail Neurons
International Journal of Theoretical and Mathematical Physics 2015, 5(5): 145-150 DOI: 10.5923/j.ijtmp.20150505.08 Numerical Simulation of Bistability between Regular Bursting and Chaotic Spiking in a
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 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 information{ } is an asymptotic sequence.
AMS B Perturbation Methods Lecture 3 Copyright by Hongyun Wang, UCSC Recap Iterative method for finding asymptotic series requirement on the iteration formula to make it work Singular perturbation use
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 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 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 informationMathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells
Mathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University
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 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 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 informationEvaluating Bistability in a Mathematical Model of Circadian Pacemaker Neurons
International Journal of Theoretical and Mathematical Physics 2016, 6(3): 99-103 DOI: 10.5923/j.ijtmp.20160603.02 Evaluating Bistability in a Mathematical Model of Circadian Pacemaker Neurons Takaaki Shirahata
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 informationConductance-Based Integrate-and-Fire Models
NOTE Communicated by Michael Hines Conductance-Based Integrate-and-Fire Models Alain Destexhe Department of Physiology, Laval University School of Medicine, Québec, G1K 7P4, Canada A conductance-based
More informationCanards at Folded Nodes
Canards at Folded Nodes John Guckenheimer and Radu Haiduc Mathematics Department, Ithaca, NY 14853 For Yulij Il yashenko with admiration and affection on the occasion of his 60th birthday March 18, 2003
More informationMixed mode oscillations as a mechanism for pseudo-plateau bursting
J Comput Neurosci (2010) 28:443 458 DOI 10.1007/s10827-010-0226-7 Mixed mode oscillations as a mechanism for pseudo-plateau bursting Theodore Vo Richard Bertram Joel Tabak Martin Wechselberger Received:
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 informationSynchronization and Phase Oscillators
1 Synchronization and Phase Oscillators Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 Synchronization
More informationA Study of the Van der Pol Equation
A Study of the Van der Pol Equation Kai Zhe Tan, s1465711 September 16, 2016 Abstract The Van der Pol equation is famous for modelling biological systems as well as being a good model to study its multiple
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 informationA Mathematical Study of Electrical Bursting in Pituitary Cells
A Mathematical Study of Electrical Bursting in Pituitary Cells Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Collaborators on
More informationSimulation of Cardiac Action Potentials Background Information
Simulation of Cardiac Action Potentials Background Information Rob MacLeod and Quan Ni February 7, 2 Introduction The goal of assignments related to this document is to experiment with a numerical simulation
More information4 Insect outbreak model
4 Insect outbreak model In this lecture I will put to a good use all the mathematical machinery we discussed so far. Consider an insect population, which is subject to predation by birds. It is a very
More informationLocalized activity patterns in excitatory neuronal networks
Localized activity patterns in excitatory neuronal networks Jonathan Rubin Amitabha Bose February 3, 2004 Abstract. The existence of localized activity patterns, or bumps, has been investigated in a variety
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 informationChaos and Dynamical Systems
Chaos and Dynamical Systems y Megan Richards Astract: In this paper, we will discuss the notion of chaos. We will start y introducing certain mathematical concepts needed in the understanding of chaos,
More informationReduction of a model of an excitable cell to a one-dimensional map
Physica D 202 (2005) 37 59 Reduction of a model of an excitable cell to a one-dimensional map Georgi S. Medvedev Department of Mathematics, Drexel University, 3141 Chestnut Street, Philadelphia, PA 19104,
More informationTowards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University
Towards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University Dynamical systems with multiple time scales arise naturally in many domains. Models of neural systems
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 informationAPPM 2360 Project 3 Mathematical Investigation of Cardiac Dynamics
APPM 2360 Project 3 Mathematical Investigation of Cardiac Dynamics Due: Thursday, December 6, 2018 by 4:59 p.m. Submit as a PDF to Assignments on Canvas 1 Introduction Cardiac Arrhythmia, or irregular
More informationMultimodal regimes in a compartmental model of the dopamine neuron
Physica D 94 (4) 333 356 Multimodal regimes in a compartmental model of the dopamine neuron Georgi S. Medvedev a,, Jaime E. Cisternas b a Department of Mathematics, Drexel University, 34 Chestnut Street,
More informationThreshold Curve for the Excitability of Bidimensional Spiking Neurons
Threshold Curve for the Excitability of Bidimensional Spiking Neurons Arnaud Tonnelier INRIA-Grenoble 655 avenue de l Europe, Montbonnot 38334 Saint Ismier, France (Dated: July 29, 214) We shed light on
More informationOn the behavior of a neural oscillator electrically coupled to a bistable element
ELSEVIER Physica D 121 (1998) 367-395 On the behavior of a neural oscillator electrically coupled to a bistable element Nancy Kopell a,,, L.F. Abbott b, Cristina Soto-Trevifio a a Department of Mathematics,
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 informationJournal of Mathematical Neuroscience (2014) 4:11 DOI /
Journal of Mathematical Neuroscience (2014) 4:11 DOI 10.1186/2190-8567-4-11 RESEARCH OpenAccess Frequency Preference Response to Oscillatory Inputs in Two-dimensional Neural Models: A Geometric Approach
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 informationTravel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis
Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis DAVID KATOSHEVSKI Department of Biotechnology and Environmental Engineering Ben-Gurion niversity of the Negev
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 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 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 informationPhase Response Curves, Delays and Synchronization in Matlab
Phase Response Curves, Delays and Synchronization in Matlab W. Govaerts and B. Sautois Department of Applied Mathematics and Computer Science, Ghent University, Krijgslaan 281-S9, B-9000 Ghent, Belgium
More informationHow Adaptation Makes Low Firing Rates Robust
Journal of Mathematical Neuroscience (2017) 7:4 DOI 10.1186/s13408-017-0047-3 RESEARCH OpenAccess How Adaptation Makes Low Firing Rates Robust Arthur S. Sherman 1 Joon Ha 1 Received: 1 December 2016 /
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 informationarxiv: v3 [math.ds] 12 Jun 2013
Isostables, isochrons, and Koopman spectrum for the action-angle representation of stable fixed point dynamics A. Mauroy, I. Mezic, and J. Moehlis Department of Mechanical Engineering, University of California
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 informationEffects of Betaxolol on Hodgkin-Huxley Model of Tiger Salamander Retinal Ganglion Cell
Effects of Betaxolol on Hodgkin-Huxley Model of Tiger Salamander Retinal Ganglion Cell 1. Abstract Matthew Dunlevie Clement Lee Indrani Mikkilineni mdunlevi@ucsd.edu cll008@ucsd.edu imikkili@ucsd.edu Isolated
More informationSpike-adding canard explosion of bursting oscillations
Spike-adding canard explosion of bursting oscillations Paul Carter Mathematical Institute Leiden University Abstract This paper examines a spike-adding bifurcation phenomenon whereby small amplitude canard
More information7 Two-dimensional bifurcations
7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed
More informationA showcase of torus canards in neuronal bursters
Journal of Mathematical Neuroscience (2012) 2:3 DOI 10.1186/2190-8567-2-3 RESEARCH OpenAccess A showcase of torus canards in neuronal bursters John Burke Mathieu Desroches Anna M Barry Tasso J Kaper Mark
More informationA Quantitative Approximation Scheme for the Traveling Wave Solutions in the Hodgkin Huxley Model
Biophysical Journal Volume 79 December 2000 2893 2901 2893 A Quantitative Approximation Scheme for the Traveling Wave Solutions in the Hodgkin Huxley Model C. B. Muratov Department of Mathematical Sciences,
More informationarxiv: v3 [math.ds] 27 Nov 2013
Modeling the modulation of neuronal bursting: a singularity theory approach arxiv:1305.7364v3 [math.ds] 27 Nov 2013 Alessio Franci 1,, Guillaume Drion 2,3, & Rodolphe Sepulchre 2,4 1 INRIA Lille-Nord Europe,
More informationCHAPTER 4 RELAXATION OSCILLATORS WITH TIME DELAY COUPLING
CHAPTER 4 RELAXATION OSCILLATORS WITH TIME DELAY COUPLING 4.1 Introduction Thus far we have studied synchronization in networks of locally coupled neurobiologically based oscillator models. One additional
More informationNumerical Computation of Canards
Numerical Computation of Canards John Guckenheimer Kathleen Hoffman Warren Weckesser January 9, 23 Shortened title: Computing Canards Abstract Singularly perturbed systems of ordinary differential equations
More informationMath 216 Second Midterm 28 March, 2013
Math 26 Second Midterm 28 March, 23 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationNeuronal 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 informationLecture 8: Kinematics: Path and Trajectory Planning
Lecture 8: Kinematics: Path and Trajectory Planning Concept of Configuration Space c Anton Shiriaev. 5EL158: Lecture 8 p. 1/20 Lecture 8: Kinematics: Path and Trajectory Planning Concept of Configuration
More informationWhen Transitions Between Bursting Modes Induce Neural Synchrony
International Journal of Bifurcation and Chaos, Vol. 24, No. 8 (2014) 1440013 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414400136 When Transitions Between Bursting Modes Induce
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 information1. Introduction - Reproducibility of a Neuron. 3. Introduction Phase Response Curve. 2. Introduction - Stochastic synchronization. θ 1. c θ.
. Introduction - Reproducibility of a euron Science (995 Constant stimuli led to imprecise spike trains, whereas stimuli with fluctuations produced spike trains with timing reproducible to less than millisecond.
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Supplementary Figure S1. Pulses >3mJ reduce membrane resistance in HEK cells. Reversal potentials in a representative cell for IR-induced currents with laser pulses of 0.74 to
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 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