This a model from the article: Parasympathetic modulation of sinoatrial node pacemaker activity in rabbit heart: a unifying model. Demir SS, Clark JW, Giles WR. Am J Physiol 1999 Jun;276(6 Pt 2):H2221-44 10362707 , Abstract: We have extended our compartmental model [Am. J. Physiol. 266 (Cell Physiol. 35): C832-C852, 1994] of the single rabbit sinoatrial node (SAN) cell so that it can simulate cellular responses to bath applications of ACh and isoprenaline as well as the effects of neuronally released ACh. The model employs three different types of muscarinic receptors to explain the variety of responses observed in mammalian cardiac pacemaking cells subjected to vagal stimulation. The response of greatest interest is the ACh-sensitive change in cycle length that is not accompanied by a change in action potential duration or repolarization or hyperpolarization of the maximum diastolic potential. In this case, an ACh-sensitive K+ current is not involved. Membrane hyperpolarization occurs in response to much higher levels of vagal stimulation, and this response is also mimicked by the model. Here, an ACh-sensitive K+ current is involved. The well-known phase-resetting response of the SAN cell to single and periodically applied vagal bursts of impulses is also simulated in the presence and absence of the beta-agonist isoprenaline. Finally, the responses of the SAN cell to longer continuous trains of periodic vagal stimulation are simulated, and this can result in the complete cessation of pacemaking. Therefore, this model is 1) applicable over the full range of intensity and pattern of vagal input and 2) can offer biophysically based explanations for many of the phenomena associated with the autonomic control of cardiac pacemaking. This model was taken from the CellML repository and automatically converted to SBML. The original model was: Demir SS, Clark JW, Giles WR. (1999) - version03 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.