|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1Department of Biomedical Engineering and Institute for Computational Medicine, Johns Hopkins University, Baltimore, Maryland; 2Center for Biomedical Computing, Simula Research Laboratory, Oslo, Norway; and 3Faculty of Kinesiology, University of Calgary, Calgary, Alberta, Canada
Submitted May 4, 2009 ; accepted in final form July 20, 2009
Ongoing investigation of the electrophysiology and pathophysiology of the human atria requires an accurate representation of the membrane dynamics of the human atrial myocyte. However, existing models of the human atrial myocyte action potential do not accurately reproduce experimental observations with respect to the kinetics of key repolarizing currents or rate dependence of the action potential and fail to properly enforce charge conservation, an essential characteristic in any model of the cardiac membrane. In addition, recent advances in experimental methods have resulted in new data regarding the kinetics of repolarizing currents in the human atria. The goal of this study was to develop a new model of the human atrial action potential, based on the Nygren et al. model of the human atrial myocyte and newly available experimental data, that ensures an accurate representation of repolarization processes and reproduction of action potential rate dependence and enforces charge conservation. Specifically, the transient outward K+ current (It) and ultrarapid rectifier K+ current (IKur) were newly formulated. The inwardly recitifying K+ current (IK1) was also reanalyzed and implemented appropriately. Simulations of the human atrial myocyte action potential with this new model demonstrated that early repolarization is dependent on the relative conductances of It and IKur, whereas densities of both IKur and IK1 underlie later repolarization. In addition, this model reproduces experimental measurements of rate dependence of It, IKur, and action potential duration. This new model constitutes an improved representation of excitability and repolarization reserve in the human atrial myocyte and, therefore, provides a useful computational tool for future studies involving the human atrium in both health and disease.
ionic model; repolarization; potassium current
Abbreviations: aur, Activation gating variable for IKur AP, Action potential APD, AP duration APD20, APD at 20% repolariation APD30, APD at 30% repolarization APD90, APD at 90% repolarization [Ca2+]c, Ca2+ concentration in the extracellular cleft space [Ca2+]d, Ca2+ concentration in the restricted subsarcolemmal space [Ca2+]i, Intracellular Ca2+ concentration [Ca2+]rel, Ca2+ concentration in the sarcoplasmic reticulum release compartment [Ca2+]up, Ca2+ concentration in the sarcoplasmic reticulum uptake compartment Cm, Cell capacitance dL, Activation gating variable for ICa,L dV/dtmax, Maximum upstroke velocity EK, Equilibrium (Nernst) potential for K+ fL1, Fast inactivation gating variable for ICa,L fL2, Slow inactivation gating variable for ICa,L F, Faraday's constant F1, Relative amount of "inactive precursor" in the Irel formulation F2, Relative amount of "activator" in the Irel formulation gKur, Maximum conductance for IKur gt, Maximum conductance for It gK1, Maximum conductance for IK1 h1, Fast inactivation variable for INa h2, Slow inactivation variable for INa hAMr model, Human atrial myocyte with new repolarization model iur, Inactivation gating variable for IKur IBCa, Background Ca2+ current IBNa, Background Na+ current ICa,L, L-type Ca2+ current ICaP, Ca2+ pump current, sarcolemmal IK1, Inwardly rectifying K+ current IK(ACh), ACh-sensitive K+ current IKr, Rapid delayed rectifier K+ current (human ether-a-go-go-related gene) IKs, Slow delayed rectifier K+ current IKur, Ultrarapid delayed rectifier K+ current INa, Na+ current INaCa, Na+/Ca2+ exchange current INaK, Na+-K+ pump current, sarcolemmal INaK,max, Maximum Na+-K+ pump current Istim, Stimulus current It, Transient outward K+ current kCa, Half-maximum Ca2+ binding concentration for fCa [K+]c, K+ concentration in the extracellular cleft space [K+]i, Intracellular K+ concentration m, Activation gating variable for INa n, Activation gating variable for IKs [Na+]c, Na+ concentration in the extracellular cleft space [Na+]i, Intracellular Na+ concentration OC, Fractional occupancy of the calmodulin-Ca2+ buffer by Ca2+ OCalse, Fractional occupancy of the calsequestrin-Ca2+ buffer by Ca2+ OTC, Fractional occupancy of the troponin-Ca2+ buffer by Ca2+ OTMgC, Fractional occupancy of the troponin-Mg2+ buffer by Ca2+ OTMgMg, Fractional occupancy of the troponin-Mg2+ buffer by Mg2+ pa, Activation gating variable for IKr PNa, Permeability for INa r, Activation gating variable for It R, Universal gas constant RMP, Resting membrane potential s, Inactivation gating variable for It T, Absolute temperature V, Membrane voltage Voli, Total myocyte volume
Na,en, Electroneutral Na+ influx
aur, Activation time constant for IKur
iur, Inactivation time constant for IKur
r, Activation time constant for It
s, Inactivation time constant for It
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |