Sinoatrial node heterogeneity and fibroblasts increase atrial driving capability in a two-dimensional human computational model.

Cardiac pacemaking remains an unsolved matter from many perspectives. Extensive experimental and computational studies have been performed to describe the sinoatrial physiology across different scales, from the molecular to clinical levels. Nevertheless, the mechanism by which a heartbeat is generated inside the sinoatrial node and propagated to the working myocardium is not fully understood at present. This work aims to provide quantitative information about this fascinating phenomenon, especially regarding the contributions of cellular heterogeneity and fibroblasts to sinoatrial node automaticity and atrial driving. We developed a bidimensional computational model of the human right atrial tissue, including the sinoatrial node. State-of-the-art knowledge of the anatomical and physiological aspects was adopted during the design of the baseline tissue model. The novelty of this study is the consideration of cellular heterogeneity and fibroblasts inside the sinoatrial node for investigating the manner by which they tune the robustness of stimulus formation and conduction under different conditions (baseline, ionic current blocks, autonomic modulation, and external high-frequency pacing). The simulations show that both heterogeneity and fibroblasts significantly increase the safety factor for conduction by more than 10% in almost all the conditions tested and shorten the sinus node recovery time after overdrive suppression by up to 60%. In the human model, especially under challenging conditions, the fibroblasts help the heterogeneous myocytes to synchronise their rate (e.g. -82% in under 25 nM of acetylcholine administration) and capture the atrium (with 25% L-type calcium current block). However, the anatomical and gap junctional coupling aspects remain the most important model parameters that allow effective atrial excitations. Despite the limitations to the proposed model, this work suggests a quantitative explanation to the astonishing overall heterogeneity shown by the sinoatrial node.