Voltage-clamp and Hodgkin-Huxley models

Similar documents
Voltage-clamp and Hodgkin-Huxley models

Single-Compartment Neural Models

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

9 Generation of Action Potential Hodgkin-Huxley Model

BME 5742 Biosystems Modeling and Control

3.3 Simulating action potentials

MATH 3104: THE HODGKIN-HUXLEY EQUATIONS

Simulation of Cardiac Action Potentials Background Information

ACTION POTENTIAL. Dr. Ayisha Qureshi Professor MBBS, MPhil

5.4 Modelling ensembles of voltage-gated ion channels

Electrophysiology of the neuron

Topics in Neurophysics

Channels can be activated by ligand-binding (chemical), voltage change, or mechanical changes such as stretch.

Deconstructing Actual Neurons

9 Generation of Action Potential Hodgkin-Huxley Model

SUPPLEMENTARY INFORMATION

Physiology Unit 2. MEMBRANE POTENTIALS and SYNAPSES

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

Lecture 11 : Simple Neuron Models. Dr Eileen Nugent

Decoding. How well can we learn what the stimulus is by looking at the neural responses?

Νευροφυσιολογία και Αισθήσεις

3 Action Potentials - Brutal Approximations

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

Introduction and the Hodgkin-Huxley Model

Effects of Betaxolol on Hodgkin-Huxley Model of Tiger Salamander Retinal Ganglion Cell

Particles with opposite charges (positives and negatives) attract each other, while particles with the same charge repel each other.

On Parameter Estimation for Neuron Models

Neural Modeling and Computational Neuroscience. Claudio Gallicchio

Signal processing in nervous system - Hodgkin-Huxley model

Bo Deng University of Nebraska-Lincoln UNL Math Biology Seminar

From neuronal oscillations to complexity

tre of Mark Louie D. Lop

Physiology Unit 2. MEMBRANE POTENTIALS and SYNAPSES

Rahaf Nasser mohammad khatatbeh

they give no information about the rate at which repolarization restores the

Lecture goals: Learning Objectives

Nervous System: Part II How A Neuron Works

Hodgkin-Huxley model simulator

BIOELECTRIC PHENOMENA

Chapter 2 The Hodgkin Huxley Theory of Neuronal Excitation

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

6.3.4 Action potential

Resting Membrane Potential

STUDIES OF CHANGING MEMBRANE POTENTIAL * : 1. BASIC ELECTRICAL THEORY, 2. GRADED AND ACTION POTENTIALS 3. THE VOLTAGE CLAMP AND MEMBRANE POTENTIALS

Basic elements of neuroelectronics -- membranes -- ion channels -- wiring. Elementary neuron models -- conductance based -- modelers alternatives

UC San Diego UC San Diego Electronic Theses and Dissertations

Supplementary Information

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

Quantitative Electrophysiology

Action Potential (AP) NEUROEXCITABILITY II-III. Na + and K + Voltage-Gated Channels. Voltage-Gated Channels. Voltage-Gated Channels

Compartmental Modelling

Neural Conduction. biologyaspoetry.com

Parameters for Minimal Model of Cardiac Cell from Two Different Methods: Voltage-Clamp and MSE Method

1 Introduction and neurophysiology

Basic elements of neuroelectronics -- membranes -- ion channels -- wiring

LESSON 2.2 WORKBOOK How do our axons transmit electrical signals?

Limitations of the Hodgkin-Huxley Formalism: Effects of Single Channel Kinetics on Transmembrane Voltage Dynamics

Chapter 24 BIFURCATIONS

Quantitative Electrophysiology

Membrane Protein Channels

Resting membrane potential,

Dynamics and complexity of Hindmarsh-Rose neuronal systems

Biological Modeling of Neural Networks

BRIEF COMMUNICATION OF ASYMMETRY CURRENT SQUID AXON MEMBRANE FREQUENCY DOMAIN ANALYSIS

Cellular Electrophysiology. Cardiac Electrophysiology

Supplemental Information. C. elegans AWA Olfactory Neurons Fire. Calcium-Mediated All-or-None Action Potentials

4. Active Behavior of the Cell Membrane 4.1 INTRODUCTION

Electricity in biology

Nerve Signal Conduction. Resting Potential Action Potential Conduction of Action Potentials

Neurophysiology. Danil Hammoudi.MD

Math 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction

Nerve Excitation. L. David Roper This is web page

Dynamical Systems for Biology - Excitability

! Depolarization continued. AP Biology. " The final phase of a local action

9.01 Introduction to Neuroscience Fall 2007

Peripheral Nerve II. Amelyn Ramos Rafael, MD. Anatomical considerations

Neurons. The Molecular Basis of their Electrical Excitability

Conductance-Based Integrate-and-Fire Models

Neurons, Synapses, and Signaling

Numerical Simulation of Bistability between Regular Bursting and Chaotic Spiking in a Mathematical Model of Snail Neurons

Ch. 5. Membrane Potentials and Action Potentials

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

Theory of action potentials. Romain Brette

Lecture 10 : Neuronal Dynamics. Eileen Nugent

FRTF01 L8 Electrophysiology

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

Synaptic dynamics. John D. Murray. Synaptic currents. Simple model of the synaptic gating variable. First-order kinetics

Action Potentials and Synaptic Transmission Physics 171/271

Biomedical Instrumentation

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

Naseem Demeri. Mohammad Alfarra. Mohammad Khatatbeh

Ch 8: Neurons: Cellular and Network Properties, Part 1

Single-Cell and Mean Field Neural Models

Chapter 3. Voltage- Dependent Membrane Permeability. Overview. Ionic Currents Across Nerve Cell Membranes

Structure and Measurement of the brain lecture notes

لجنة الطب البشري رؤية تنير دروب تميزكم

Probabilistic Modeling of Action Potential Generation by Neurons

Neuronal Dynamics: Computational Neuroscience of Single Neurons

Modeling Action Potentials in Cell Processes

Membrane Currents in Mammalian Ventricular Heart Muscle Fibers Using a Voltage-Clamp Technique

Transcription:

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 (1998 (for a recent version of the HH squid axon model Rothman and Manis, J. Neurophysiol. 89:3070, 3083 and 3097 (2003 for examples of separation of currents by voltage clamp. Ion channel properties: Selectivity Rectification Saturation and block by toxins and other ions Gating voltage-gating (Hille chapts. 3-5 ligand-gating (Hille chapt. 6-7 sensory (Hille chapt. 8 today 1

The membrane model (parallel currents through all the channels present: I cap + I K + I Na + I leak = I ext C dv dt = I ext G K (V E K G Na (V E Na G leak (V E leak I ext inside I Cap I K I Na I leak G K G Na G leak V C + E K + E Na + E leak What are the properties of these conductances in nerve membrane? i.e. G na (V,t outside To separate the currents for study, first the capacitive current is eliminated using voltage clamp. A step of current I clamp is applied to force the membrane potential to adopt a constant value V clamp. V clamp V rest dv dt = 0 here C dv dt 0 = I clamp G K (V E K G Na (V E Na G leak (V E leak so I clamp = G K (V E K + G Na (V E Na + G leak (V E leak 2

Recording from single channels and ensembles of channels using voltage-clamp: single channels gate randomly, whole cell currents are the sum of the singlechannel currents and behave smoothly. I whole cell = j i jthchannel I whole cell N i single channel P open (V, t, V M (small capacitive transients at the time of the voltage step are subtracted out Hille, p. 90 Voltage clamp in the squid giant axon membrane from the original work of Hodgkin and Huxley. The trace 100% Na + is the recording in normal solutions. This is a mixture of inward Na currents and outward K currents. The currents are separated by setting the extracellular Na to 10% of its normal value, effectively eliminating the inward current and revealing the K current (I K. The difference of the 100% and 10% currents is the Na current (I Na. Hille, p. 39 3

Current-voltage relationships for the HH sodium and potassium channels. I K = G K (V,t (V E K and I Na = G Na (V,t (V E Na From these data, the conductances G K and G Na can be computed... V clamp E K Compare these currents to the expectation from the GHK constant-field equation. E Na Hille, p. 39... yielding the following dependence of conductance on time and V clamp V clamp Note that both the amplitude of the conductance change and its time constant change with V clamp Hille, p. 42 4

With the data shown in previous slides in hand, Hodgkin and Huxley modeled the conductances by assuming them to be proportional to one or two activation and inactivation variables, representing the fraction of gates that are open. In the case of the K conductance, only activation n is needed. G K = G K n 4 (V,t dn dt = n (V 1 n τ n (V 1 The behavior of G K is determined by functions n (V and t n (V. assuming that n(0 = n (V r [ ] 1 exp( t τ n (V 1 n(t= n (V r + n (V 1 n (V r V r = -60 mv n (V" V 1 = -30 mv n 0.5 τ n (V 1 n 4 n (V 1 n (V 1 r " resting value -10 0 10 20 Time re voltage step, ms The exponent 4 is chosen to make the rise of G K sigmoidal. An interpretation of the HH potassium channel model: Recall that the potassium channel has 4 subunits. If each subunit has an independent gate and the probability that any one gate is open is n then the probability that all 4 gates are open (i.e. the channel is open is n 4. 4α C 4 β C 3 O 3α C 2 O 2α 2 C 1 O α 3 O 4 2β 3β 4β Channel is open 5

For Na, an activation variable m and an inactivation variable h are needed. The former is a gate that opens with depolarization and the latter is a gate that closes with depolarization. G Na = G Na m 3 h dm dt = m (V m τ m (V and dh dt = h (V h τ h (V m tail current h Hille, p. 48 The HH model: C dv dt = I ext G Na (V E Na G K (V E K G leak (V E leak G Na = G Na m 3 h G K = G K n 4 dm dt = m (V m τ m (V dn dt = n (V n τ n (V dh dt = h (V h τ h (V 6

The functions n (V, m (V, and h (V determine whether gates serve to activate channels (conventionally, open the channel with depolarization or inactivate the channel (close the channel with depolarization. τ m, τ h, and τ n are the time constants. These shapes are expected from a barrier model of gating (see Homework. resting potential the m and n gates open with depolarization the h gate closes with depolarization Hille, p. 49 Reconstruction of the action potential by the HH model : 1 Depolarization of the cell (by an injected current in this case leads to 1 2 a self-sustaining increase in m (V, m, I Na, and V, which leads to 2 3 a decrease in h (V and an increase in n (V. The resulting decrease in h and increase in n terminate the action potential and repolarize the membrane. Note the difference in the response times of m (fast versus n and h (slow. (AP produced by a 1 ms, 9 ma current pulse at the heavy bar in the V plot 3 3 Milliseconds 7

An explanation for anode break excitation. During a hyperpolarization, membrane excitability is increased by decreases in K channel activation and Na channel inactivation. Note, this means that the model does not have a fixed threshold. The difference between instantaneous I-V curves and those dependent on gating. Consider the sodium channel: I Na = G Na (V,t (V E Na This equation applies instantaneously after a voltage clamp with a constant conductance, as well as over time as gating proceeds. E K E Na 8

In neuron cell bodies, there are often many more types of channels. The voltage-clamp data at right are from four types of cells in the cochlear nucleus. The Na + and Ca ++ currents were blocked using tetrodotoxin and Cd ++, respectively, so only K + currents are present. These records were decomposed into three different K + channels: I HT - like the HH K + channel, actually carried by two different channels. I LT - also like the HH K + channel, but with a low activation threshold. Specifically blocked by dendrotoxin I A - a K + channel with inactivation. I HT and I LT I HT I A and I HT Rothman and Manis, 2003 To see how the separation proceeds, consider the case of I A. No pharmacological blockers were available. This current is only activated when its inactivation is removed by hyperpolarizing. Hyperpolarizing prepulse, I A is present No hyperpolarizing prepulse, I A is not seen. The current is entirely I HT. I A derived from the difference between A and B. Rothman and Manis, 2003 9

The A-channel is fit by an HH model with activation a and two inactivations b and c. I A = G A a 4 bc(v E K Data Model Rothman and Manis, 2003 Action potentials are not always so simple. Many neurons fire complex spikes, often accompanied by simple spikes. This requires additional (Ca ++ channels, not included in HH. Manis et al. 1994 10