This a model from the article: Ion currents underlying sinoatrial node pacemaker activity: a new single cell mathematical model. Dokos S, Celler B, Lovell N. J Theor Biol 1996 Aug 7;181(3):245-72 8869126 , Abstract: The ionic currents underlying autorhythmicity of the mammalian sinoatrial node and their wider contribution to each phase of the action potential have been investigated in this study using a new single cell mathematical model. The new model provides a review and update of existing formulations of sinoatrial node membrane currents, derived from a wide range of electrophysiological data available in the literature. Simulations of spontaneous activity suggest that the dominant mechanism underlying pacemaker depolarisation is the inward background Na+ current, ib,Na. In contrast to previous models, the decay of the delayed rectifying K+ current, iK, was insignificant during this phase. Despite the presence of a pseudo-outward background current throughout the pacemaker range of potentials (Na-K pump+leak currents), the hyperpolarisation-activated current i(f) was not essential to pacemaker activity. A closer inspection of the current-voltage characteristics of the model revealed that the "instantaneous" time-independent current was inward for holding potentials in the pacemaker range, which rapidly became outward within 2 ms due to the inactivation of the L-type Ca2+ current, iCa,L. This suggests that reports in the literature in which the net background current is outward throughout the pacemaker range of potentials may be exaggerated. The magnitudes of the action potential overshoot and the maximum diastolic potential were determined largely by the reversal potentials of iCa,L and iK respectively. The action potential was sustained by the incomplete deactivation of iCa,L and the Na-Ca exchanger, iNaCa. Despite the incorporation of "square-root" activation by [K]o of all K+ currents, the model was unable to correctly simulate the response to elevated [K]o. This model was taken from the CellML repository and automatically converted to SBML. The original model was: Dokos S, Celler B, Lovell N. (1996) - version06 The original CellML model was created by: Lloyd, Catherine, May c.lloyd@aukland.ac.nz The University of Auckland The Bioengineering Institute This model originates from BioModels Database: A Database of Annotated Published Models (http://www.ebi.ac.uk/biomodels/). It is copyright (c) 2005-2011 The BioModels.net Team. To the extent possible under law, all copyright and related or neighbouring rights to this encoded model have been dedicated to the public domain worldwide. Please refer to CC0 Public Domain Dedication for more information. In summary, you are entitled to use this encoded model in absolutely any manner you deem suitable, verbatim, or with modification, alone or embedded it in a larger context, redistribute it, commercially or not, in a restricted way or not.. To cite BioModels Database, please use: Li C, Donizelli M, Rodriguez N, Dharuri H, Endler L, Chelliah V, Li L, He E, Henry A, Stefan MI, Snoep JL, Hucka M, Le Novère N, Laibe C (2010) BioModels Database: An enhanced, curated and annotated resource for published quantitative kinetic models. BMC Syst Biol., 4:92.