Page 9 - Cardiac Electrophysiology | A Modeling and Imaging Approach
P. 9

P. 9
                                                                      3
                                                                          3
                                             g g Na Na  Na  g •    • m m • h•   h                               (2.4)
                                                g = g g=
                                                          Na
                                                             Na Na
               Where m is an activation gate, h an inactivation gate and ḡ      Na is the maximum membrane
        conductance for the Na  ion species (when all gates are open). I  current activation (increase of
                                  +
                                                                             Na
        conductance) is simulated accurately with three identical m gates that transition from 0 (closed)
        to 1 (open) upon V  depolarization. The time dependence of m follows first-order kinetics:
                            m




                                          dm dm  = α   •  ( −1  m  )− β    •  m

                                           dt                                                                   (2.5)
                                           dt



               Where m and (1-m) are the gate open and closed probabilities, respectively; α and ß are
        opening and closing transition rates that depend on V . Because the transitions of each gate are
                                                                   m
        assumed independent of the other gates, the probability that all m gates are open is m . Upon
                                                                                                       3
        depolarization, these gates open very rapidly, generating a large inward I  that generates the rapid
                                                                                       Na
        upstroke of the action potential. The voltage clamp recordings of Hodgkin and Huxley indicated
        that I  decreases shortly after activation (a process termed inactivation). This process is
              Na
        represented by the parameter h, which is also governed by a first-order differential equation similar
        to equation (2.5). At hyperpolarized V , h is fully open; it closes when the membrane is depolarized
                                                m
        to cause a monoexponential decrease of I  during voltage-clamp protocols, as observed
                                                      Na
        experimentally. Based on the independent gating assumption, h transitions are independent of m
        and the open probability for I  is m h. For the axon, a similar formalism was used to compute a
                                               3
                                        Na
        potassium current I . The combination of inward I  and outward I  (together with a leakage
                             K                                Na                K
        current I ) was sufficient to reproduce the morphology of the action potential under a variety of
                 L
        conditions.


         2.2  Computing the Action Potential from Ion Channels Kinetic Properties



               In the Hodgkin Huxley scheme, I  in equation (2.1) is the macroscopic current through an
                                                   ion
        ensemble of many ion channels. The gating variables are model parameters that do not
        correspond to specific kinetic properties of the ion channel, nor to its molecular structure. In fact,

        at the time when Hodgkin and Huxley formulated the model, the mechanism of voltage and time
        dependent gating was not known and ion-channel proteins were not identified as the pathway
        for transmembrane flow of ions. Based on current knowledge of ion channels and their structure/
        function properties, it is apparent that the Hodgkin-Huxley formalism is limited in its ability to

        describe various aspects of channel behavior. The gating parameters do not represent specific
        kinetic states of the channel and the independent gating assumption precludes dependence of a
        given transition on occupancy of the different states of the channel (for example, sodium channel
   4   5   6   7   8   9   10   11   12   13   14