AUTHOR=Yang Dongmei , Morrell Christopher H. , Lyashkov Alexey E. , Tagirova Sirenko Syevda , Zahanich Ihor , Yaniv Yael , Vinogradova Tatiana M. , Ziman Bruce D. , Maltsev Victor A. , Lakatta Edward G. TITLE=Ca2+ and Membrane Potential Transitions During Action Potentials Are Self-Similar to Each Other and to Variability of AP Firing Intervals Across the Broad Physiologic Range of AP Intervals During Autonomic Receptor Stimulation JOURNAL=Frontiers in Physiology VOLUME=Volume 12 - 2021 YEAR=2021 URL=https://www.frontiersin.org/journals/physiology/articles/10.3389/fphys.2021.612770 DOI=10.3389/fphys.2021.612770 ISSN=1664-042X ABSTRACT=Ca2+-Vm transitions occurring throughout AP cycles in sinoatrial nodal (SAN) cells are cues that: not only regulate activation states of molecules operating within criticality (Ca2+-domain) and limit-cycle (Vm-domain) mechanisms of a coupled-clock system that underlies SAN cell automaticity; but are also regulated by the activation states of the clock molecules they regulate. These cues are both causes and effects of clock molecular activation (recursion). Recently, we demonstrated that Ca2+-Vm transitions during AP cycles in single SAN cells from different species are self-similar (obey power-law) and self-similar to trans-species AP firing intervals in vitro, to heart rate in vivo, and to body mass. Neurotransmitter stimulation of β-adrenergic/cholinergic-receptor initiated signaling in SAN cells modulates their AP firing rate and rhythm by impacting on the degree to which SAN clocks couple to each other, creating the broad physiologic range of SAN cell mean AP firing intervals and firing interval variabilities. Here we show that Ca2+-Vm domain kinetic transitions (time-to-AP-ignition in diastole and 90% AP recovery) occurring within given AP, the mean AP firing intervals, and AP firing interval variabilities within time-series of APs in 230 individual SAN cells are self-similar (obey power-laws). These long-range correlations inform on self-similar distributions of order among SAN cells across the entire broad physiologic range of SAN AP firing intervals, regardless of whether autonomic receptors of these cells are stimulated or not, and regardless of the type (adrenergic or cholinergic) of autonomic receptor stimulation. These long-range correlations among distributions of Ca2+-Vm kinetic functions that regulate SAN cell clock coupling during each AP cycle in different individual, isolated SAN cells not in contact with each other. Our numerical model simulations further extended our perspectives to the molecular scale and demonstrated that many ion currents also behave self-similar across autonomic states. Thus, to ensure rapid flexibility of AP firing rates in response to different types and degrees of autonomic input, nature “did not reinvent molecular wheels within the coupled-clock system of pacemaker cells”, but differentially engaged or scaled the kinetics of gears that regulate the rate and rhythm at which the “wheels spin” in a given autonomic input context.