A single pacemaker cell model based on the phase response curve

A single pacemaker cell model and its response to repetitive external depolarization stimulations is described in this paper. This model is a simple model based on the two most important functional properties of the cardiac pacemaker cells. The first property is the intrinsic pacemaker cycle length, which is an `internal' parameter of the cell, describing the most important feature of a pacemaker cell. The second functional property is the phase response curve (PRC), which is an `overall collective' function: it contains all the `information' about the possible interactions of the pacemaker cell with the outside world (external stimulus, interaction with surrounding cells, etc.). This study demonstrates that by representing the pacemaker cell only by two fundamental features, and by applying a simple physical-mathematical model, a global picture of the system can be achieved, allowing us to explore qualitatively various physiological phenomena related to the pacemaker function. For example, we demonstrated that the PRC is a crucial parameter in the prediction of the entrainment phenomena of a single pacemaker cell in response to a periodic train of depolarization pulses. Specifically, the PRC permits a quantitative determination of the 1:1 synchronization range for a single pacemaker cell and an external depolarization pulse. Moreover, we show that the PRC can be used to represent the type of external stimulus applied to the pacemaker (e.g. depolarization pulse) and its intensity. Therefore, the PRC emerges as an important determinant and a useful `tool' for the understanding of the dynamic interaction of pacemaker cells with the outside world. As a result of our simulations, we unveil a new important parameter: the `degree of influence', which determines the range of 1:1 synchronization between an external depolarization pulse and a pacemaker cell. This interaction parameter is a direct function of the PRC parameters. It appears to be a helpful `tool' for the understanding of synchronization and mutual entrainment mechanisms between the pacemaker cell and an external stimulus, and therefore it supports the basic importance of the PRC in the description and determination of these mechanisms.

[1]  Sarit Abramovich-Sivan,et al.  A pacemaker cell pair model based on the phase response curve , 1998, Biological Cybernetics.

[2]  J Jalife,et al.  Effects of current flow on pacemaker activity of the isolated kitten sinoatrial node. , 1980, The American journal of physiology.

[3]  Sarit Abramovich-Sivan,et al.  A simulation of the SA node by a phase response curve-based model of a two-dimensional pacemaker cells array , 2000, IEEE Transactions on Biomedical Engineering.

[4]  L. Glass,et al.  Global bifurcations of a periodically forced biological oscillator , 1984 .

[5]  T Sato,et al.  Difference equation model of ventricular parasystole as an interaction between cardiac pacemakers based on the phase response curve. , 1983, Journal of theoretical biology.

[6]  L. Glass,et al.  Phase locking, period-doubling bifurcations, and irregular dynamics in periodically stimulated cardiac cells. , 1981, Science.

[7]  C Antzelevitch,et al.  Electrotonic Modulation of Pacemaker Activity Further Biological and Mathematical Observations on the Behavior of Modulated Parasystole , 1982, Circulation.

[8]  T Sano,et al.  Mechanism of rhythm determination among pacemaker cells of the mammalian sinus node. , 1978, The American journal of physiology.

[9]  C Antzelevitch,et al.  Phase resetting and annihilation in a mathematical model of sinus node. , 1985, The American journal of physiology.

[10]  L Glass,et al.  Phase-locked rhythms in periodically stimulated heart cell aggregates. , 1988, The American journal of physiology.

[11]  S Akselrod,et al.  A PRC based model of a pacemaker cell: effect of vagal activity and investigation of the respiratory sinus arrhythmia. , 1998, Journal of theoretical biology.

[12]  Sarit Abramovich-Sivan,et al.  Simulation of atrial activity by a phase response curve based model of a two-dimensional pacemaker cells array: the transition from a normal activation pattern to atrial fibrillation , 1999, Biological Cybernetics.

[13]  L. Glass,et al.  From Clocks to Chaos: The Rhythms of Life , 1988 .

[14]  A E Becker,et al.  Functional and Morphological Organization of the Rabbit Sinus Node , 1980, Circulation research.

[15]  J Jalife,et al.  Effect of Electrotonic Potentials on Pacemaker Activity of Canine Purkinje Fibers in Relation to Parasystole , 1976, Circulation research.

[16]  J Jalife,et al.  Dynamic Interactions and Mutual Synchronization of Sinoatrial Node Pacemaker Cells A Mathematical Model , 1986, Circulation research.

[17]  J Jalife,et al.  A Mathematical Model of Parasystole and its Application to Clinical Arrhythmias , 1977, Circulation.

[18]  J Jalife,et al.  Mutual entrainment and electrical coupling as mechanisms for synchronous firing of rabbit sino‐atrial pace‐maker cells. , 1984, The Journal of physiology.

[19]  A. Winfree The geometry of biological time , 1991 .