Simulations of the Cardiac Action Potential Based on the Hodgkin-Huxley Kinetics with the Use of Microsoft Excel Spreadsheets

Sheng-Nan Wu

simulation ； cardiac action potential ； Excel ； mathematical model

The Chinese Journal of Physiology

47卷1期（2004 / 03 / 01）

15 - 22

英文

The purpose of this study was to develop a method to simulate the cardiac action potential using a Microsoft Excel spreadsheet. The mathematical model contained voltage-gated ionic currents that were modeled using either Beeler-Reuter (B-R) or Luo-Rudy (L-R) phase 1 kinetics. The simulation protocol involves the use of in-cell formulas directly typed into a spreadsheet. The capability of spreadsheet iteration was used in these simulations. It does not require any prior knowledge of computer programming, although the use of the macro language can speed up the calculation. The normal configuration of the cardiac ventricular action potential can be well simulated in the B-R model that is defined by four individual ionic currents, each representing the diffusion of ions through channels in the membrane. The contribution of Na(superscript +) inward current to the rate of depolarization is reproduced in this model. After removal of Na(superscript +) current from the model, a constant current stimulus elicits an oscillatory change in membrane potential. In the L-R phase 1 model where six types of ionic currents were defined, the effect of extracellular K(superscript +) concentration on changes both in the time course of repolarization and in the time-independent K(superscript +) current can be demonstrated, when the solutions are implemented in Excel. Using the simulation protocols described here, the users can readily study and graphically display the underlying properties of ionic currents to see how changes in these properties determine the behavior of the heart cell. The method employed in these simulation protocols may also be extended or modified to other biological simulation programs.

1. Beeler, G.W.,Reuter, H.(1977).Reconstruction of the action potential of ventricular myocardial fibres.J. Physiol,268,177-210.
2. Bower, J.M.,Beeman, D.(1997).The Book of Genesis, Springer.Berlin:
3. Brown, A.M.(1999).A methodology for simulating biological systems using Microsoft Excel.Comp. Methods Frog. Biomed,58,181-190.
4. Campbell, D.L.,Giles, W.R.,Robinson. K.,Shibata, E.F.(1988).Studies of the sodium-calcium exchanger in bull-frog atrial myocytes.J. Physiol,403,317-340.
5. Demir, S.S.,Clark, J.W.,Murphey, C.R.,Giles, W.R.(1994).A mathematical model of a rabbit sinoatrial node cells.Am. J. Physiol,266,C832-C852.
6. Ferrero, J.M. Jr.,Saiz, J.,Ferrero, J.M.,Thakor, N.V.(1996).Simulation of action potentials from metabolically impaired cardiac myocytes. Role of ATP-sensitive K+ current.Circ. Res,79,208-221.
7. Giles, W.,Nakajima, T.,Ono, K.,Shibata, E.F.(1989).Modulation of the delayed rectifier K+ current by isoprenaline in bull-frog atrial myocytes.J. Physiol,415,233-249.
8. Hines, M.L.,Carnevale, N.T.(1997).The NEURON simulation environment.Neural Comput,9,1179-1209.
9. Hodgkin, A.L.,Huxley, A.F.(1952).A quantitative description of membrane current and its application to current pulse.J. Physiol,117,500-544.
10. Irisawa, H.,Brown, H.F.,Giles, W.(1993).Cardiac pacemaking in the sinoatrial node.Physiol. Rev,73,197-227.
11. Isenberg, G.,Klockner, U.(1982).Calcium currents of isolated bovine ventricular myocytes are fast and of large amplitude.Pfiugers Arch,395,30-41.
12. Kiyosue, T.,Arita, M.,Muramatsu, H.,Soindler, A.J.,Noble, D.(1993).Ionic mechanisms of action potential prolongation at low temperatures in guinea-pig ventricular myocytes.J. Physiol,468,85-106.
13. Kootsey, J.M.,Kohn, M.C.,Feezor, M.D.,Mitchell, G.R.,Fletcher, P.R.(1986).SCoP: an interactive simulation control program for micro and minicomputers.Bull. Math. Biol,48,427-441.
14. Lindblad, D.S.,Murphey, C.R.,Clark, J.W.,Giles, W.R.(1996).A model of the action potential and underlying membrane currents in a rabbit atrial cell.Am. J. Physiol,271,H1666-H1696.
15. Luo, C.H.,Rudy, Y.(1991).A model of the ventricular cardiac action potential: depolarization, repolarization and their interaction.Circ. Res,68,1501-1526.
16. Luo, C.H.,Rudy, Y.(1994).A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes.Circ. Res,74,1071-1096.
17. McCullough, B.D.,Wilson, B.(2002).On the accuracy of statistical procedures in Microsoft Excel 2000 and Excel XP.Comput. Statist. Data Anal,40,713-721.
18. Nakajima, T.,Wu, S.,Irisawa, H.,Giles, W.(1990).Mechanism of acetylcholine-induced inhibition of Ca2+ current in bullfrog atrial myocytes.J. Gen. Physiol,96,865-885.
19. Noble, D.(2002).Modelling the heart: insights, failure and progress.Bioessays,24,1155-1163.
20. Nygren, A.,Giles, W.R.(2000).Mathematical simulation of slowing of cardiac conduction velocity by elevated extracellular potassium.Ann. Biomed. Eng,28,951-957.
21. Revest, P.(1995).Neurosim for Windows.Trends Neurosci,18,556.
22. Tseng, G.N.,Robinson, R.B.,Hoffmann, B.F.(1987).Passive proper ties and membrane currents of canine ventricular myocytes.J. Gen. Physiol,90,671-701.
23. Wu, S.N.,Liu, S.I.,Hwang, T.L.(1998).Activation of muscarinic K+ channels by extracellular ATP and UTP in rat atrial myocytes.J. Cardiovasc. Pharmacol,31,203-211.
24. Wu, S.N.,Lue, S.I.,Yang, S.L.,Hsu, H.K.,Liu. M.S.(1993).Electrophysiologic properties of isolated adult cardiomyocytcs from septic rats.Circ. Shock,41,239-247.
25. Wu, S.N.,Nakajima, T.,Yamashita, T.,Hamada, E.,Hazama. H.,Iwasawa, K.,Omata, M.,Kurachi, Y.(1994).Molecular niechanism of cibenzoline-induced anticholinergic action in single atrial myocytes: comparison with effect of disopyramide.J. Cardiovasc.Pharrnacol,23,618-623.
26. Yanagihara, K.,Noma, A.,Irisawa, H.(1980).Reconstruction of sinoatrial node pacemaker potential based on the voltage clamp experiments.Jpn. J. Physiol,30,841-885.
27. Zeng, J.,Laurita, K.R.,Rosenbaum, D.S.,Rudy, Y.(1995).Two components of the delayed rectifier K current in ventricular niyocytes of the guinea pig type. Theoretical formulation and their role in repolarization.Circ. Res,77,140-152.

1. Wu, Sheng-nan,Lo, Yi-ching,Liu, Yen-chin,Chen, Bing-shuo(2010).Effects of Transient Receptor Potential-Like Current on the Firing Pattern of Action Potentials in the Hodgkin-Huxley Neuron during Exposure to Sinusoidal External Voltage.中國生理學雜誌,53(6),423-429.
2. Wu, Sheng-Nan,Shen, Ai-Yu,Lo, Yi-Ching,Chen, Bing-Shuo(2011).Contribution of Non-Inactivating Na+ Current Induced by Oxidizing Agents to the Firing Behavior of Neuronal Action Potentials: Experimental and Theoretical Studies from NG108-15 Neuronal Cells.中國生理學雜誌,54(1),19-29.
3. Yang, Chia-Jung,Wu, Sheng-Nan,Lee, Yu-Chi,Chang, Tzu-Tung(2018).Stimulatory Action of Telmisartan, an Antagonist of Angiotensin II Receptor, on Voltage-Gated Na+ Current: Experimental and Theoretical Studies.中國生理學雜誌,61(1),1-13.