From neuronal oscillations to complexity

Size: px
Start display at page:

Download "From neuronal oscillations to complexity"

Transcription

1 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 Le Havre, LMAH July 24th, 2008

2 2/39 Biological introduction Table of contents Biological introduction Mathematical models Numerical analysis Conclusion and perspectives

3 3/39 Biological introduction Neuron a cell body, extended by : dendrites an axon an axon terminal

4 4/39 Biological introduction Ionic channels Membrane is composed of : a double lipidic layer ionic channels Ionic channels : voltage-dependant open / close / inactive conductance

5 5/39 Biological introduction Some important channels Sodium channels : 4 gates : 3 to control the opening 1 to control the inactivity particularity : inactivity period Potassium channels : 4 gates all control the opening particularity : delay Leakage channels : particularity : always opened

6 6/39 Biological introduction Resting potential Potential difference between both sides of the membrane ; more Na + outside more K + inside Experimentally, 90mV < ddp < 50mV intra-cellular space more negative extra-cellular space more positive

7 7/39 Biological introduction Action potential

8 8/39 Mathematical models Table of contents Biological introduction Mathematical models Numerical analysis Conclusion and perspectives

9 9/39 Mathematical models Hodgkin-Huxley model (1) Alan Lloyd Hodgkin and Andrew Fielding Huxley : neurophysiologists experiments on squid giant neuron ; discovery of ionic channels ; description of electric activity of a neuron by mathematical equations ; Nobel price of physiology and medecine in 1952 ;

10 10/39 Mathematical models Hodgkin-Huxley model (2) C dv = I K + I Na + I L I dt dm = α m (1 m) β m m dt dn = α n (1 n) β n n dt dh = α h (1 h) β h h dt

11 11/39 Mathematical models Hodgkin-Huxley model (3) I = C dv dt + I K + I Na + I L C dv dt = I K + I Na + I L I I K = n 4 g K (E E K ) ; n = one of the 4 gates of a potassium channel is open ; g K = conductance of a potassium channel number of opened channels ; E K = equilibrium potential of potassium (= Nernst potential) ;

12 12/39 Mathematical models Hodgkin-Huxley model (4) Similarly, we have : I Na = m 3 hg Na (E E Na ) ; I L = g L (E E L ) ; α i = opening coefficient, i = m, n, h ; β i = closing coefficient, i = m, n, h ;

13 13/39 Mathematical models Hodgkin-Huxley model (5) C dv = I Na + I K + I L I dt dm = α m (1 m) β m m dt dn = α n (1 n) β n n dt dh = α h (1 h) β h h dt

14 14/39 Mathematical models Hindmarsh-Rose 1982 model (1) From Hodgkin-Huxley model to Hindmarsh-Rose 1982 model ; Mathematical simplification thanks to biological observations : m is substituted by a constant because of its fast activation ; h + n = 0.8 (experimentally) ; C dv = I K + I Na + I L I dt dm = α m (1 m) β m m dt dn = α n (1 n) β n n dt dh = α h (1 h) β h h dt

15 15/39 Mathematical models Hindmarsh-Rose 1982 model (2) A Fitzhugh-Nagumo type model : dx dt dy dt = y x 3 + ax 2 + I = 1 dx 2 y a, d : parameters determined experimentally. I : applied current.

16 16/39 Mathematical models Hindmarsh-Rose 1984 model A third equation to be closer to reality : an adaptation equation dx = y x 3 + ax 2 z + I dt dy = 1 dx 2 y dt dz = ɛ(b(x x c ) z), ɛ << 1 dt slow-fast system

17 17/39 Numerical analysis Table of contents Biological introduction Mathematical models Numerical analysis Conclusion and perspectives

18 18/39 Numerical analysis Numerical analysis (1) Time series (t, x) : bursting, and projections of the phase portrait on the planes (x, y), (x, z) et (y, z)

19 19/39 Numerical analysis Numerical analysis (2) Bifurcation Qualitative change of solutions due to parameters changes I = 1.5: 3-cycle I = 2 : 5-cycle I = 3: 10-cycle 3.25 I : chaos I : stop generating bursts

20 20/39 Numerical analysis Numerical analysis (3) I = 1.5, I = 2, I = 2.5, I = 3 I = 3.25, I = 3.33 I = 3.35 I = 3.50

21 21/39 Table of contents Biological introduction Mathematical models Numerical analysis Conclusion and perspectives

22 22/39 -Generalities characteristic of many processes in natural systems and non linear science syn : common, chronos : time to share the same motion at the same time there exist many kind of synchronization (identical, generalized, phase...) neurons : optimal transmission of information

23 Coupling systems (1) Coupling Making different entities be dependant one from the others. Linear coupling function Mutual coupling of 2 neurons by a linear function : electrical connection x 1 = ax1 2 x 1 3 y 1 z 1 k(x 2 x 1 ) y 1 = (a + α)x1 2 y 1 z 1 = µ(bx 1 + c z 1 ) x 2 = ax2 2 x 2 3 y 2 z 2 k(x 1 x 2 ) y 2 = (a + α)x2 2 y 2 z 2 = µ(bx 2 + c z 2 ) 23/39 k(x 1 x 2 ) : coupling function, k : coupling strength

24 24/39 Coupling systems (2) Sigmoid coupling function Mutual coupling of 2 neurons by a sigmoid function : chemical connection x 1 = ax1 2 x 1 3 y 1 1 z 1 (x 1 V s )k exp( λ(x 2 Θ s )) y 1 = (a + α)x1 2 y 1 z 1 = µ(bx 1 + c z 1 ) x 2 = ax2 2 x 2 3 y 1 2 z 2 (x 2 V s )k exp( λ(x 1 Θ s )) y 2 = (a + α)x2 2 y 2 z 2 = µ(bx 2 + c z 2 ) Θ s : threshold, k 2 : coupling strength.

25 25/39 Coupling 2 neurons Two neurons identically synchronize if : x 2 x 1 y 2 y 1 z 2 z 1 x = x 1 x 2 0 y = y 1 y 2 0 z = z 1 z 2 0 Stability study of the transverse system : x = x 1 x 2 y = y 1 y 2 z = z 1 Example of theoretical tool to study this system : Lyapunov function. z 2

26 26/39 HR neurons mutually coupled by a sigmoid function x 1 = ax1 2 x 1 3 y 1 1 z 1 (x 1 V s )k exp( λ(x 2 Θ s )) y 1 = (a + α)x1 2 y 1 z 1 = µ(bx 1 + c z 1 ) x 2 = ax2 2 x 2 3 y 1 2 z 2 (x 2 V s )k exp( λ(x 1 Θ s )) y 2 = (a + α)x2 2 y 2 z 2 = µ(bx 2 + c z 2 )

27 27/39 Numerical results Remark Diagonal x 1 = x 2 synchronization x 1 vs x 2 for k 2 = 0.8,k 2 = 1,k 2 = 1.22,k 2 = Same minimal coupling strength to obtain synchronization between y 1 and y 2 and between z 1 and z 2.

28 28/39 Coupling 3 neurons x 1 vs x 2 for k 3 = 0.2,k 3 = 0.4,k 3 = 0.6,k 3 = same resuts for x i vs x j, y i vs y j, z i vs z j ( i, j = 1, 2, 3)

29 29/39 Coupling 4 neurons x 1 vs x 2 for k 4 = 0.1,k 4 = 0.3,k 4 = 0.42,k 4 = 0.5. same results for x i vs x j, y i vs y j, z i vs z j ( i, j = 1,..., 4)

30 30/39 Coupling 5 neurons x 1 vs x 2 for k 5 = 0.1,k 5 = 0.2,k 5 = 0.29,k 5 = 0.3. same results for x i vs x j, y i vs y j, z i vs z j ( i, j = 1,..., 5)

31 Coupling 6 neurons x 1 vs x 2 for k 6 = 0.1,k 6 = 0.2,k 6 = 0.23,k 6 = /39 same results for x i vs x j, y i vs y j, z i vs z j ( i, j = 1,..., 6)

32 Coupling 7 neurons x 1 vs x 2 for k 7 = 0.1,k 7 = 0.17,k 7 = 0.2,k 7 = /39 same results for x i vs x j, y i vs y j, z i vs z j ( i, j = 1,..., 7)

33 Coupling 8 neurons x 1 vs x 2 for k 8 = 0.1,k 8 = 0.15,k 8 = 0.16,k 8 = /39 same results for x i vs x j, y i vs y j, z i vs z j ( i, j = 1,..., 8)

34 34/39 Coupling n neurons each neuron is connected to all the others in a network of n neurons each neuron is connected to n 1 neurons (n 1 inputs) k n = k 2 n 1 k2 : 2 neurons mutually coupled synchronization threshold n 1 : number of inputs for each neuron

35 35/39 classical law which is found in many self-organized complex systems

36 36/39 Examples earthquakes : 1000 of magnitude 4 for 100 of magnitude 5 and 10 of magnitude 6 linguistic : for 1000 occurrences of the in text, 100 occurrences of I and 10 of say urban systems : big cities are rare and small ones are frequent in an exponential way.

37 37/39 Conclusion and perspectives Table of contents Biological introduction Mathematical models Numerical analysis Conclusion and perspectives

38 38/39 Conclusion and perspectives Conclusion and perspectives intersection between synchronization in coupled nonlinear systems and complex networks study conditions under which a complex network of dynamical systems synchronizes study of the different topological structures of networks and their consequences on the synchronization phenomenon synchronization and complex networks theoretical study

39 39/39 Conclusion and perspectives THANK YOU!

Dynamics and complexity of Hindmarsh-Rose neuronal systems

Dynamics 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 information

1 Introduction and neurophysiology

1 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 information

Nonlinear Observer Design and Synchronization Analysis for Classical Models of Neural Oscillators

Nonlinear 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

Lecture 10 : Neuronal Dynamics. Eileen Nugent

Lecture 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 information

Voltage-clamp and Hodgkin-Huxley models

Voltage-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 information

Math 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 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 information

Electrophysiology of the neuron

Electrophysiology 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 information

Voltage-clamp and Hodgkin-Huxley models

Voltage-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 information

Single neuron models. L. Pezard Aix-Marseille University

Single 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 information

Signal processing in nervous system - Hodgkin-Huxley model

Signal 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 information

On Parameter Estimation for Neuron Models

On 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 information

Mathematical Foundations of Neuroscience - Lecture 3. Electrophysiology of neurons - continued

Mathematical 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 information

Modeling Action Potentials in Cell Processes

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 information

Lecture Notes 8C120 Inleiding Meten en Modelleren. Cellular electrophysiology: modeling and simulation. Nico Kuijpers

Lecture 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

Topics in Neurophysics

Topics 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 information

STUDENT PAPER. Santiago Santana University of Illinois, Urbana-Champaign Blue Waters Education Program 736 S. Lombard Oak Park IL, 60304

STUDENT 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 information

Dynamical Systems in Neuroscience: Elementary Bifurcations

Dynamical 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 information

9 Generation of Action Potential Hodgkin-Huxley Model

9 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 information

Chimera states in networks of biological neurons and coupled damped pendulums

Chimera 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 information

BME 5742 Biosystems Modeling and Control

BME 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 information

Bursting Oscillations of Neurons and Synchronization

Bursting 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 information

Biomedical Instrumentation

Biomedical Instrumentation ELEC ENG 4BD4: Biomedical Instrumentation Lecture 5 Bioelectricity 1. INTRODUCTION TO BIOELECTRICITY AND EXCITABLE CELLS Historical perspective: Bioelectricity first discovered by Luigi Galvani in 1780s

More information

1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model

1 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

Dynamical 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 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 information

Electronics 101 Solving a differential equation Incorporating space. Numerical Methods. Accuracy, stability, speed. Robert A.

Electronics 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 information

6.3.4 Action potential

6.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 information

BIOELECTRIC PHENOMENA

BIOELECTRIC 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 information

ACTION POTENTIAL. Dr. Ayisha Qureshi Professor MBBS, MPhil

ACTION 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 information

Introduction and the Hodgkin-Huxley Model

Introduction 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 information

Action Potential Propagation

Action 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 information

Simulation of Cardiac Action Potentials Background Information

Simulation 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 information

Neural Modeling and Computational Neuroscience. Claudio Gallicchio

Neural 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 information

Single-Compartment Neural Models

Single-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 information

Compartmental Modelling

Compartmental 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

3.3 Simulating action potentials

3.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 information

FRTF01 L8 Electrophysiology

FRTF01 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 information

MATH 3104: THE HODGKIN-HUXLEY EQUATIONS

MATH 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 information

Chapter 2 The Hodgkin Huxley Theory of Neuronal Excitation

Chapter 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 information

Bo Deng University of Nebraska-Lincoln UNL Math Biology Seminar

Bo 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 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 information

Neuroscience applications: isochrons and isostables. Alexandre Mauroy (joint work with I. Mezic)

Neuroscience 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 information

Presented by Sarah Hedayat. Supervised by Pr.Cappy and Dr.Hoel

Presented 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 information

Slow Manifold of a Neuronal Bursting Model

Slow Manifold of a Neuronal Bursting Model Slow Manifold of a Neuronal Bursting Model Jean-Marc Ginoux 1 and Bruno Rossetto 2 1 PROTEE Laboratory, Université du Sud, B.P. 2132, 83957, La Garde Cedex, France, ginoux@univ-tln.fr 2 PROTEE Laboratory,

More information

9 Generation of Action Potential Hodgkin-Huxley Model

9 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 information

Biological Modeling of Neural Networks

Biological 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 information

All-or-None Principle and Weakness of Hodgkin-Huxley Mathematical Model

All-or-None Principle and Weakness of Hodgkin-Huxley Mathematical Model All-or-None Principle and Weakness of Hodgkin-Huxley Mathematical Model S. A. Sadegh Zadeh, C. Kambhampati International Science Index, Mathematical and Computational Sciences waset.org/publication/10008281

More information

Mathematical analysis of a 3D model of cellular electrophysiology

Mathematical 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 information

Lokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt

Lokta-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 information

Chapter 24 BIFURCATIONS

Chapter 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 information

Resting membrane potential,

Resting 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 information

Probabilistic Modeling of Action Potential Generation by Neurons

Probabilistic Modeling of Action Potential Generation by Neurons Probabilistic Modeling of Action Potential Generation by Neurons Student: Bobby Sena under the direction of Professors Helen Wearing and Eric Toolson University of New Mexico bsena@unm.edu May 6, 21 Abstract

More information

Quantitative Electrophysiology

Quantitative Electrophysiology ECE 795: Quantitative Electrophysiology Notes for Lecture #1 Wednesday, September 13, 2006 1. INTRODUCTION TO EXCITABLE CELLS Historical perspective: Bioelectricity first discovered by Luigi Galvani in

More information

LIMIT CYCLE OSCILLATORS

LIMIT 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 information

Dynamical Systems for Biology - Excitability

Dynamical 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 information

Deconstructing Actual Neurons

Deconstructing 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 information

Nernst Equilibrium Potential. p. 1

Nernst Equilibrium Potential. p. 1 Nernst Equilibrium Potential p. 1 Diffusion The conservation law for a compound with concentration c: rate change of c = local production + accumulation due to transport. Model: d c dv = p dv J n da dt

More information

Single-Cell and Mean Field Neural Models

Single-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 information

Computational Neuroscience Summer School Neural Spike Train Analysis. An introduction to biophysical models (Part 2)

Computational 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 information

Simulating Hodgkin-Huxley-like Excitation using Comsol Multiphysics

Simulating Hodgkin-Huxley-like Excitation using Comsol Multiphysics Presented at the COMSOL Conference 2008 Hannover Simulating Hodgkin-Huxley-like Excitation using Comsol Multiphysics Martinek 1,2, Stickler 2, Reichel 1 and Rattay 2 1 Department of Biomedical Engineering

More information

Nerve Excitation. L. David Roper This is web page

Nerve Excitation. L. David Roper  This is web page Nerve Excitation L. David Roper http://arts.bev.net/roperldavid This is web page http://www.roperld.com/science/nerveexcitation.pdf Chapter 1. The Action Potential A typical action potential (inside relative

More information

Introduction to Physiology II: Control of Cell Volume and Membrane Potential

Introduction to Physiology II: Control of Cell Volume and Membrane Potential Introduction to Physiology II: Control of Cell Volume and Membrane Potential J. P. Keener Mathematics Department Math Physiology p.1/23 Basic Problem The cell is full of stuff: Proteins, ions, fats, etc.

More information

Quantitative Electrophysiology

Quantitative Electrophysiology ECE 795: Quantitative Electrophysiology Notes for Lecture #1 Tuesday, September 18, 2012 1. INTRODUCTION TO EXCITABLE CELLS Historical perspective: Bioelectricity first discovered by Luigi Galvani in 1780s

More information

Cellular Electrophysiology. Cardiac Electrophysiology

Cellular 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 information

Final 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 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 information

5.4 Modelling ensembles of voltage-gated ion channels

5.4 Modelling ensembles of voltage-gated ion channels 5.4 MODELLING ENSEMBLES 05 to as I g (Hille, 200). Gating currents tend to be much smaller than the ionic currents flowing through the membrane. In order to measure gating current, the ionic current is

More information

Integral Equation and PDE Neuronal Models:

Integral 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 information

Chaos in the Hodgkin-Huxley Equations: The Takens-Bodganov Cusp Bifurcation. Rami Mohieddine

Chaos in the Hodgkin-Huxley Equations: The Takens-Bodganov Cusp Bifurcation. Rami Mohieddine CORNELL UNIVERSITY MATHEMATICS DEPARTMENT SENIOR THESIS Chaos in the Hodgkin-Huxley Equations: The Takens-Bodganov Cusp Bifurcation A THESIS PRESENTED IN PARTIAL FULFILLMENT OF CRITERIA FOR HONORS IN MATHEMATICS

More information

Bifurcation Analysis of Neuronal Bursters Models

Bifurcation Analysis of Neuronal Bursters Models Bifurcation Analysis of Neuronal Bursters Models B. Knowlton, W. McClure, N. Vu REU Final Presentation August 3, 2017 Overview Introduction Stability Bifurcations The Hindmarsh-Rose Model Our Contributions

More information

Memory and hypoellipticity in neuronal models

Memory 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 information

2013 NSF-CMACS Workshop on Atrial Fibrillation

2013 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 information

Strange Nonchaotic Spiking in the Quasiperiodically-forced Hodgkin-Huxley Neuron

Strange Nonchaotic Spiking in the Quasiperiodically-forced Hodgkin-Huxley Neuron Journal of the Korean Physical Society, Vol. 57, o. 1, July 2010, pp. 23 29 Strange onchaotic Spiking in the Quasiperiodically-forced Hodgkin-Huxley euron Woochang Lim and Sang-Yoon Kim Department of Physics,

More information

1 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 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 information

Modelling biological oscillations

Modelling biological oscillations Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van

More information

Supratim Ray

Supratim Ray Supratim Ray sray@cns.iisc.ernet.in Biophysics of Action Potentials Passive Properties neuron as an electrical circuit Passive Signaling cable theory Active properties generation of action potential Techniques

More information

General Details Acknowledgements

General Details Acknowledgements Prof Waters Mathematical Physiology Notes. 1 General Details Acknowledgements The following notes are based on material generously passed to me by Professors Gaffney, Chapman and Fowler and Dr Hinch, which

More information

Dynamic Systems: Ordinary Differential Equations. Ordinary Differential Equations

Dynamic Systems: Ordinary Differential Equations. Ordinary Differential Equations Dynamic Systems: Ordinary Differential Equations Adapted From: Numerical Methods in Biomedical Engineering Stanley M. Dunn, Alkis Constantinides, Prabhas V. Moghe Chapter 7 Kim Ferlin and John Fisher Ordinary

More information

Learning Cycle Linear Hybrid Automata for Excitable Cells

Learning Cycle Linear Hybrid Automata for Excitable Cells Learning Cycle Linear Hybrid Automata for Excitable Cells Sayan Mitra Joint work with Radu Grosu, Pei Ye, Emilia Entcheva, I V Ramakrishnan, and Scott Smolka HSCC 2007 Pisa, Italy Excitable Cells Outline

More information

Lecture 11 : Simple Neuron Models. Dr Eileen Nugent

Lecture 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 information

1 Single neuron computation and the separation of processing and signaling

1 Single neuron computation and the separation of processing and signaling Physics 178/278 - David Kleinfeld - Winter 2017 - Week 1 1 Single neuron computation and the separation of processing and signaling This overview serves to connect simple notions of statistical physics

More information

Membrane Potentials, Action Potentials, and Synaptic Transmission. Membrane Potential

Membrane Potentials, Action Potentials, and Synaptic Transmission. Membrane Potential Cl Cl - - + K + K+ K + K Cl - 2/2/15 Membrane Potentials, Action Potentials, and Synaptic Transmission Core Curriculum II Spring 2015 Membrane Potential Example 1: K +, Cl - equally permeant no charge

More information

An Investigation of the Coefficient of Variation Using the Dissipative Stochastic Mechanics Based Neuron Model

An Investigation of the Coefficient of Variation Using the Dissipative Stochastic Mechanics Based Neuron Model Available Online at www.ijcsmc.com International Journal of Computer Science and Mobile Computing A Monthly Journal of Computer Science and Information Technology IJCSMC, Vol. 2, Issue. 10, October 2013,

More information

Transmission of Nerve Impulses (see Fig , p. 403)

Transmission 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 information

Ch. 5. Membrane Potentials and Action Potentials

Ch. 5. Membrane Potentials and Action Potentials Ch. 5. Membrane Potentials and Action Potentials Basic Physics of Membrane Potentials Nerve and muscle cells: Excitable Capable of generating rapidly changing electrochemical impulses at their membranes

More information

Tracking the State of the Hindmarsh-Rose Neuron by Using the Coullet Chaotic System Based on a Single Input

Tracking the State of the Hindmarsh-Rose Neuron by Using the Coullet Chaotic System Based on a Single Input ISSN 1746-7659, England, UK Journal of Information and Computing Science Vol. 11, No., 016, pp.083-09 Tracking the State of the Hindmarsh-Rose Neuron by Using the Coullet Chaotic System Based on a Single

More information

Parameter Reduction of the Hodgkin-Huxley Model. of Action Potential Propagation. Tyler A. Bahr

Parameter Reduction of the Hodgkin-Huxley Model. of Action Potential Propagation. Tyler A. Bahr Parameter Reduction of the Hodgkin-Huxley Model of Action Potential Propagation Tyler A. Bahr A senior thesis submitted to the faculty of Brigham Young University in partial fulfillment of the requirements

More information

Electricity in biology

Electricity in biology Electricity in biology Maintaining an osmotic balance Different strategies Animal: pumping ions out Plants/some bacteria: using a rigid shell Amoeba: extruding water bubble periodically Animal Plant/Bacteria

More information

UC San Diego UC San Diego Electronic Theses and Dissertations

UC San Diego UC San Diego Electronic Theses and Dissertations UC San Diego UC San Diego Electronic Theses and Dissertations Title Single neuron modeling and data assimilation in BNST neurons Permalink https://escholarship.org/uc/item/5vr6560f Author Farsian, Reza

More information

Action Potentials and Synaptic Transmission Physics 171/271

Action Potentials and Synaptic Transmission Physics 171/271 Action Potentials and Synaptic Transmission Physics 171/271 Flavio Fröhlich (flavio@salk.edu) September 27, 2006 In this section, we consider two important aspects concerning the communication between

More information

MEMBRANE POTENTIALS AND ACTION POTENTIALS:

MEMBRANE POTENTIALS AND ACTION POTENTIALS: University of Jordan Faculty of Medicine Department of Physiology & Biochemistry Medical students, 2017/2018 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Review: Membrane physiology

More information

Action Potentials & Nervous System. Bio 219 Napa Valley College Dr. Adam Ross

Action Potentials & Nervous System. Bio 219 Napa Valley College Dr. Adam Ross Action Potentials & Nervous System Bio 219 Napa Valley College Dr. Adam Ross Review: Membrane potentials exist due to unequal distribution of charge across the membrane Concentration gradients drive ion

More information

/639 Final Solutions, Part a) Equating the electrochemical potentials of H + and X on outside and inside: = RT ln H in

/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 information

How and why neurons fire

How and why neurons fire How and why neurons fire 1 Neurophysiological Background The Neuron Contents: Structure Electrical Mebrane Properties Ion Channels Actionpotential Signal Propagation Synaptic Transission 2 Structure of

More information

Nervous System: Part II How A Neuron Works

Nervous System: Part II How A Neuron Works Nervous System: Part II How A Neuron Works Essential Knowledge Statement 3.E.2 Continued Animals have nervous systems that detect external and internal signals, transmit and integrate information, and

More information

Modeling cells in tissue by connecting. electrodiffusion and Hodgkin Huxley-theory.

Modeling 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 information

9.01 Introduction to Neuroscience Fall 2007

9.01 Introduction to Neuroscience Fall 2007 MIT OpenCourseWare http://ocw.mit.edu 9.01 Introduction to Neuroscience Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 9.01 Recitation (R02)

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY 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 information

Ch 7. The Nervous System 7.1 & 7.2

Ch 7. The Nervous System 7.1 & 7.2 Ch 7 The Nervous System 7.1 & 7.2 SLOs Describe the different types of neurons and supporting cells, and identify their functions. Identify the myelin sheath and describe how it is formed in the CNS and

More information

Passive Membrane Properties

Passive Membrane Properties Passive Membrane Properties Communicating through a leaky garden hose... Topics I Introduction & Electrochemical Gradients Passive Membrane Properties Action Potentials Voltage-Gated Ion Channels Topics

More information

Neuronal Dynamics: Computational Neuroscience of Single Neurons

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 information

Bursting and Chaotic Activities in the Nonlinear Dynamics of FitzHugh-Rinzel Neuron Model

Bursting and Chaotic Activities in the Nonlinear Dynamics of FitzHugh-Rinzel Neuron Model Bursting and Chaotic Activities in the Nonlinear Dynamics of FitzHugh-Rinzel Neuron Model Abhishek Yadav *#, Anurag Kumar Swami *, Ajay Srivastava * * Department of Electrical Engineering, College of Technology,

More information