AJP - Heart and Circulatory Physiology, Vol 259, Issue 2 352-H369, Copyright © 1990 by American Physiological Society

ARTICLES

A mathematical model of a bullfrog cardiac pacemaker cell

R. L. Rasmusson, J. W. Clark, W. R. Giles, E. F. Shibata and D. L. Campbell
Department of Electrical and Computer Engineering, Rice University, Houston, Texas 77251-1892.

Previous models of cardiac cellular electrophysiology have been based largely on voltage-clamp measurements obtained from multicellular preparations and often combined data from different regions of the heart and a variety of species. We have developed a model of cardiac pacemaking based on a comprehensive set of voltage-clamp measurements obtained from single cells isolated from one specific tissue type, the bullfrog sinus venosus (SV). Consequently, sarcolemmal current densities and kinetics are not influenced by secondary phenomena associated with multicellular preparations, allowing us to realistically simulate processes thought to be important in pacemaking, including the Na(+)-K+ pump and Na(+)-Ca2+ exchanger. The membrane is surrounded extracellularly by a diffusion-limited space and intracellularly by a limited myoplasmic volume containing Ca2(+)-binding proteins (calmodulin, troponin). The model makes several predictions regarding mechanisms involved in pacing. 1) Primary pacemaking cannot be attributed to any single current but arises from both the lack of a background K+ current and a complex interaction between Ca2+, delayed-rectifier K+, and background leak currents. 2) Ca2+ current displays complex behavior and is important during repolarization. 3) Because of Ca2+ buffering by myoplasmic proteins, the Na(+)-Ca2+ exchanger current is small and has little influence on action potential repolarization but may modulate the maximum diastolic potential. 4) The Na(+)-K+ pump current does not play an active role in repolarization but is of sufficient size to modulate the rate of diastolic depolarization. 5) K+ accumulation and Ca2+ depletion may occur in the extracellular spaces but play no role in either the diastolic depolarization or repolarization of a single action potential. This model illustrates the importance of basing simulations on quantitative measurements of ionic currents in myocytes and of including both electrogenic transporter mechanisms and Ca2+ buffering by myoplasmic Ca2(+)-binding proteins.

This article has been cited by other articles:

				J. J. Torres, L. N. Cornelisse, E. G. A. Harks, W. P. M. van Meerwijk, A. P. R. Theuvenet, and D. L. Ypey Modeling action potential generation and propagation in NRK fibroblasts Am J Physiol Cell Physiol, October 1, 2004; 287(4): C851 - C865. [Abstract] [Full Text] [PDF]

				V. E. Bondarenko, G. P. Szigeti, G. C. L. Bett, S.-J. Kim, and R. L. Rasmusson Computer model of action potential of mouse ventricular myocytes Am J Physiol Heart Circ Physiol, September 1, 2004; 287(3): H1378 - H1403. [Abstract] [Full Text] [PDF]

				S. Miyoshi, H. Mitamura, K. Fujikura, Y. Fukuda, K. Tanimoto, Y. Hagiwara, M. Ita, and S. Ogawa A mathematical model of phase 2 reentry: role of L-type Ca current Am J Physiol Heart Circ Physiol, April 1, 2003; 284(4): H1285 - H1294. [Abstract] [Full Text] [PDF]

				T. Mitsuiye, Y. Shinagawa, and A. Noma Sustained Inward Current During Pacemaker Depolarization in Mammalian Sinoatrial Node Cells Circ. Res., July 21, 2000; 87(2): 88 - 91. [Abstract] [Full Text] [PDF]

				M. Courtemanche, R. J. Ramirez, and S. Nattel Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model Am J Physiol Heart Circ Physiol, July 1, 1998; 275(1): H301 - H321. [Abstract] [Full Text] [PDF]

				A. Nygren, C. Fiset, L. Firek, J. W. Clark, D. S. Lindblad, R. B. Clark, and W. R. Giles Mathematical Model of an Adult Human Atrial Cell : The Role of K+ Currents in Repolarization Circ. Res., January 23, 1998; 82(1): 63 - 81. [Abstract] [Full Text] [PDF]