The role of irregularity in intercellular connections is studied in the first natural human pacemaker called the sinoatrial node by modeling with the Greenberg–Hastings cellular automata. Facts from modern physiology about the sinoatrial node drive modeling. Heterogeneity between cell connections is reproduced by a rewiring procedure applied to a square lattice. The Greenberg–Hastings rule, representing the intrinsic cellular dynamics, is modified to imitate self-excitation of each pacemaker cell. Moreover, interactions with nearest neighbors are changed to heterogeneous ones by enhancing horizontal connections. Stationary states of the modeled system emerge as self-organized robust oscillatory states. Since the sinoatrial node role relies on a single cell cyclic activity, properties of single cells are studied. It appears that the strength and diversity of cellular oscillations depend directly on properties of intrinsic cellular dynamics. But these oscillations also depend on the underlying topology. Moderate nonuniformity of intercellular connections are found vital for proper function of the sinoatrial node, namely, for producing robust oscillatory states that are able to respond effectively to the autonomic system control.