A simulation of the SA node by a phase response curve-based model of a two-dimensional pacemaker cells array

This paper presents a simulation of the sino-atrial (SA) node by a two-dime nsional pacemaker cells array model, based on phase response curve (PRC) in teraction. This simple model of the cardiac pacemaker cells, involves only the most basic functional properties, which play a direct role in the deter mination of the SA node rhythm, The two most relevant functional properties of the pacemaker cells are: The intrinsic cycle length, an "internal" feat ure of each pacemaker cell, and the PRC, an "overall collective" function. The PRC contains the "information" about the type of interactions of each p acemaker cell with the outside world (i.e., interaction with neighboring ce lls, external stimulus, etc.), and "strength" of the interaction (strong, w eak, etc.), We studied the spatial interaction among a large number of pacemaker cells (15 x 15), as a function of the regional variation of cells properties, the "electrical" coupling between cells (the PRC), and the appearance of regio ns with abnormal cycle lengths. We investigated the influence of those para meters on the mutual interaction between the pacemaker cells, on the activa tion pattern and conduction time of the array, and on a pseudo-electrocardi oigram (ECG) signal, This study demonstrates that by representing the pacemaker cells in the SA node by only two fundamental features, and by applying a simple physical-ma thematical model, we can create a global picture of the SA node system. Thi s enables us to explore physiological phenomena related to the genesis and maintenance of the SA node activity, and to gain insight into the condition s which predispose the SA node instability, and conduction disturbances.