Excitation-contraction coupling kinetics are dictated by the rate and rhythm of the excitations generated by sinoatrial-nodal cells. These cells generate local Ca releases (LCRs) that activate Na/Ca exchanger current, which accelerates diastolic depolarization and determines the rate and rhythm of the excitations. The LCRs are generated by clusters of ryanodine receptors, Ca release units (CRUs), residing in the sarcoplasmic reticulum. While the spatial CRU distribution in pacemaker cells exhibits substantial heterogeneity, it remains unknown if it has any functional importance. Using numerical modeling, here we showed that with a square lattice distribution of CRUs, Ca-induced-Ca-release propagation during diastolic depolarization is insufficient for pacemaking within a broad lower range of realistic ICaL densities. Allowing each CRU to deviate from its original lattice position fundamentally changes the model behavior: during diastolic depolarization sparks propagate, forming LCRs observed experimentally. As disorder in the CRU positions increases, the CRU distribution exhibits larger empty spaces but simultaneously CRU clusters, as in Poisson clumping. Propagating within the clusters, Ca release becomes synchronized, increasing AP firing rate and reviving pacemaker function within lower ICaL densities at which cells with lattice CRU distribution were dormant/non-firing. However, cells with fully disordered CRU positions cannot reach low firing rates and their β-adrenergic receptor stimulation effect was substantially decreased. Thus, order/disorder in CRU locations regulates Ca release propagation and could be harnessed by pacemaker cells to regulate their function. Excessive disorder is expected to limit heart rate range that may contribute to heart rate range decline with age and in disease.SummaryThe present numerical modeling study shows that disorder in locations of Ca release units in cardiac pacemaker cells has substantial functional impact by creating release clusters, similar to Poisson clumping, and opportunity of Ca release to propagate within the clusters.