Electrophysiological heterogeneity of pacemaker cells in the rabbit intercaval region, including the SA node: insights from recording multiple ion currents in each cell

Analysis of changes in the action potential morphology of the mouse sinoatrial node true pacemaker cells during ontogenetic development in vitro and in silico

BACKGROUND

Maturation of the mouse is accompanied by the increase in heart rate. However, the mechanisms underlying this process remain unclear. We performed an action potentials (APs) recordings in mouse sinoatrial node (SAN) true pacemaker cells and in silico analysis to clarify the mechanisms underlying pre-postnatal period heart rate changes.

RESULTS

The APs of true pacemaker cells at different stages had similar configurations and dV/dtmax values. The cycle length, action potential duration (APD90), maximal diastolic potential (MDP), and AP amplitude decreased, meanwhile the velocity of diastolic depolarization (DDR) increased from E12.5 stage to adult. Using a pharmacological approach we found that in SAN true pacemaker cells ivabradine reduces the DDR and the cycle length significantly stronger in E12.5 than in newborn and adult mice, whereas the effects of Ni2+ and nifedipine were significantly stronger in adult mice. Computer simulations further suggested that the density of the hyperpolarization-activated pacemaker сurrent (If) decreased during development, whereas transmembrane and intracellular Ca2+ flows increased.

CONCLUSIONS

The ontogenetic decrease in IK1 density from E12.5 to adult leads to depolarization of MDP to the voltage range in which calcium currents are activated, thereby shifting the balance from the "membrane-clock" to the "calcium-clock."

A novel conceptual model of heart rate autonomic modulation based on a small-world modular structure of the sinoatrial node

The present view on heartbeat initiation is that a primary pacemaker cell or a group of cells in the sinoatrial node (SAN) center paces the rest of the SAN and the atria. However, recent high-resolution imaging studies show a more complex paradigm of SAN function that emerges from heterogeneous signaling, mimicking brain cytoarchitecture and function. Here, we developed and tested a new conceptual numerical model of SAN organized similarly to brain networks featuring a modular structure with small-world topology. In our model, a lower rate module leads action potential (AP) firing in the basal state and during parasympathetic stimulation, whereas a higher rate module leads during β-adrenergic stimulation. Such a system reproduces the respective shift of the leading pacemaker site observed experimentally and a wide range of rate modulation and robust function while conserving energy. Since experimental studies found functional modules at different scales, from a few cells up to the highest scale of the superior and inferior SAN, the SAN appears to feature hierarchical modularity, i.e., within each module, there is a set of sub-modules, like in the brain, exhibiting greater robustness, adaptivity, and evolvability of network function. In this perspective, our model offers a new mainframe for interpreting new data on heterogeneous signaling in the SAN at different scales, providing new insights into cardiac pacemaker function and SAN-related cardiac arrhythmias in aging and disease.

Emergent activity, heterogeneity, and robustness in a calcium feedback model of the sinoatrial node

The sinoatrial node (SAN) is the primary pacemaker of the heart. SAN activity emerges at an early point in life and maintains a steady rhythm for the lifetime of the organism. The ion channel composition and currents of SAN cells can be influenced by a variety of factors. Therefore, the emergent activity and long-term stability imply some form of dynamical feedback control of SAN activity. We adapt a recent feedback model-previously utilized to describe control of ion conductances in neurons-to a model of SAN cells and tissue. The model describes a minimal regulatory mechanism of ion channel conductances via feedback between intracellular calcium and an intrinsic target calcium level. By coupling a SAN cell to the calcium feedback model, we show that spontaneous electrical activity emerges from quiescence and is maintained at steady state. In a 2D SAN tissue model, spatial variability in intracellular calcium targets lead to significant, self-organized heterogeneous ion channel expression and calcium transients throughout the tissue. Furthermore, multiple pacemaking regions appear, which interact and lead to time-varying cycle length, demonstrating that variability in heart rate is an emergent property of the feedback model. Finally, we demonstrate that the SAN tissue is robust to the silencing of leading cells or ion channel knockouts. Thus, the calcium feedback model can reproduce and explain many fundamental emergent properties of activity in the SAN that have been observed experimentally based on a minimal description of intracellular calcium and ion channel regulatory networks.

Synergy between Membrane Currents Prevents Severe Bradycardia in Mouse Sinoatrial Node Tissue

Bradycardia is initiated by the sinoatrial node (SAN), which is regulated by a coupled-clock system. Due to the clock coupling, reduction in the 'funny' current (I_f), which affects SAN automaticity, can be compensated, thus preventing severe bradycardia. We hypothesize that this fail-safe system is an inherent feature of SAN pacemaker cells and is driven by synergy between I_f and other ion channels. This work aimed to characterize the connection between membrane currents and their underlying mechanisms in SAN cells. SAN tissues were isolated from C57BL mice and Ca²⁺ signaling was measured in pacemaker cells within them. A computational model of SAN cells was used to understand the interactions between cell components. Beat interval (BI) was prolonged by 54 ± 18% (N = 16) and 30 ± 9% (N = 21) in response to I_f blockade, by ivabradine, or sodium current (I_Na) blockade, by tetrodotoxin, respectively. Combined drug application had a synergistic effect, manifested by a BI prolonged by 143 ± 25% (N = 18). A prolongation in the local Ca²⁺ release period, which reports on the level of crosstalk within the coupled-clock system, was measured and correlated with the prolongation in BI. The computational model predicted that I_Na increases in response to I_f blockade and that this connection is mediated by changes in T and L-type Ca²⁺ channels.

The virtual sinoatrial node: What did computational models tell us about cardiac pacemaking?

Since its discovery, the sinoatrial node (SAN) has represented a fascinating and complex matter of research. Despite over a century of discoveries, a full comprehension of pacemaking has still to be achieved. Experiments often produced conflicting evidence that was used either in support or against alternative theories, originating intense debates. In this context, mathematical descriptions of the phenomena underlying the heartbeat have grown in importance in the last decades since they helped in gaining insights where experimental evaluation could not reach. This review presents the most updated SAN computational models and discusses their contribution to our understanding of cardiac pacemaking. Electrophysiological, structural and pathological aspects - as well as the autonomic control over the SAN - are taken into consideration to reach a holistic view of SAN activity.

The paradigm shift: Heartbeat initiation without "the pacemaker cell"

The current dogma about the heartbeat origin is based on "the pacemaker cell," a specialized cell residing in the sinoatrial node (SAN) that exhibits spontaneous diastolic depolarization triggering rhythmic action potentials (APs). Recent high-resolution imaging, however, demonstrated that Ca signals and APs in the SAN are heterogeneous, with many cells generating APs of different rates and rhythms or even remaining non-firing (dormant cells), i.e., generating only subthreshold signals. Here we numerically tested a hypothesis that a community of dormant cells can generate normal automaticity, i.e., "the pacemaker cell" is not required to initiate rhythmic cardiac impulses. Our model includes 1) non-excitable cells generating oscillatory local Ca releases and 2) an excitable cell lacking automaticity. While each cell in isolation was not "the pacemaker cell", the cell system generated rhythmic APs: The subthreshold signals of non-excitable cells were transformed into respective membrane potential oscillations via electrogenic Na/Ca exchange and further transferred and integrated (computed) by the excitable cells to reach its AP threshold, generating rhythmic pacemaking. Cardiac impulse is an emergent property of the SAN cellular network and can be initiated by cells lacking intrinsic automaticity. Cell heterogeneity, weak coupling, subthreshold signals, and their summation are critical properties of the new pacemaker mechanism, i.e., cardiac pacemaker can operate via a signaling process basically similar to that of "temporal summation" happening in a neuron with input from multiple presynaptic cells. The new mechanism, however, does not refute the classical pacemaker cell-based mechanism: both mechanisms can co-exist and interact within SAN tissue.

Adenosine reduces sinoatrial node cell action potential firing rate by uncoupling its membrane and calcium clocks

The spontaneous action potential (AP) firing rate of sinoatrial nodal cells (SANC) is regulated by a system of intracellular Ca²⁺ and membrane ion current clocks driven by Ca²⁺-calmodulin-activated adenylyl cyclase-protein kinase-A signaling. The mean AP-cycle length (APCL) and APCL variability inform on the effectiveness of clock coupling. Endogenous ATP metabolite adenosine binds to adenosine receptors (A₁, A₃) that couple to G_i protein-coupled receptors, reducing spontaneous AP firing rate via G_βγ signaling that activates I_KAch,Ado. Adenosine also inhibits adenylyl cyclase activity via G_αi signaling, impacting cAMP-mediated protein kinase-A-dependent protein phosphorylation. We hypothesize that in addition to I_KAch,Ado activation, adenosine impacts also Ca²⁺ via G_αi signaling and that both effects reduce AP firing rate by reducing the effectiveness of the Ca²⁺ and membrane clock coupling. To this end, we measured Ca²⁺ and membrane potential characteristics in enzymatically isolated single rabbit SANC. 10 µM adenosine substantially increased both the mean APCL (on average by 43%, n = 10) and AP beat-to-beat variability from 5.1 ± 1.7% to 7.2 ± 2.0% (n = 10) measured via membrane potential and 5.0 ± 2.2% to 10.6 ± 5.9% (n = 40) measured via Ca²⁺ (assessed as the coefficient of variability = SD/mean). These effects were mediated by hyperpolarization of the maximum diastolic membrane potential (membrane clock effect) and suppression of diastolic local Ca²⁺releases (LCRs) (Ca²⁺-clock effect): as LCR size distributions shifted to smaller values, the time of LCR occurrence during diastolic depolarization (LCR period) became prolonged, and the ensemble LCR signal became reduced. The tight linear relationship of coupling between LCR period to the APCL in the presence of adenosine "drifted" upward and leftward, i.e. for a given LCR period, APCL was prolonged, becoming non-linear indicating clock uncoupling. An extreme case of uncoupling occurred at higher adenosine concentrations (>100 µM): small stochastic LCRs failed to self-organize and synchronize to the membrane clock, thus creating a failed attempt to generate an AP resulting in arrhythmia and cessation of AP firing. Thus, the effects of adenosine to activate G_βγ and I_KACh,Ado and to activate G_αi, suppressing adenylyl cyclase activity, both contribute to the adenosine-induced increase in the mean APCL and APCL variability by reducing the fidelity of clock coupling and AP firing rate.

Coupling and heterogeneity modulate pacemaking capability in healthy and diseased two-dimensional sinoatrial node tissue models

Both experimental and modeling studies have attempted to determine mechanisms by which a small anatomical region, such as the sinoatrial node (SAN), can robustly drive electrical activity in the human heart. However, despite many advances from prior research, important questions remain unanswered. This study aimed to investigate, through mathematical modeling, the roles of intercellular coupling and cellular heterogeneity in synchronization and pacemaking within the healthy and diseased SAN. In a multicellular computational model of a monolayer of either human or rabbit SAN cells, simulations revealed that heterogenous cells synchronize their discharge frequency into a unique beating rhythm across a wide range of heterogeneity and intercellular coupling values. However, an unanticipated behavior appeared under pathological conditions where perturbation of ionic currents led to reduced excitability. Under these conditions, an intermediate range of intercellular coupling (900-4000 MΩ) was beneficial to SAN automaticity, enabling a very small portion of tissue (3.4%) to drive propagation, with propagation failure occurring at both lower and higher resistances. This protective effect of intercellular coupling and heterogeneity, seen in both human and rabbit tissues, highlights the remarkable resilience of the SAN. Overall, the model presented in this work allowed insight into how spontaneous beating of the SAN tissue may be preserved in the face of perturbations that can cause individual cells to lose automaticity. The simulations suggest that certain degrees of gap junctional coupling protect the SAN from ionic perturbations that can be caused by drugs or mutations.

A phenotype-based forward genetic screen identifies Dnajb6 as a sick sinus syndrome gene

Previously we showed the generation of a protein trap library made with the gene-break transposon (GBT) in zebrafish (Danio rerio) that could be used to facilitate novel functional genome annotation towards understanding molecular underpinnings of human diseases (Ichino et al, 2020). Here, we report a significant application of this library for discovering essential genes for heart rhythm disorders such as sick sinus syndrome (SSS). SSS is a group of heart rhythm disorders caused by malfunction of the sinus node, the heart's primary pacemaker. Partially owing to its aging-associated phenotypic manifestation and low expressivity, molecular mechanisms of SSS remain difficult to decipher. From 609 GBT lines screened, we generated a collection of 35 zebrafish insertional cardiac (ZIC) mutants in which each mutant traps a gene with cardiac expression. We further employed electrocardiographic measurements to screen these 35 ZIC lines and identified three GBT mutants with SSS-like phenotypes. More detailed functional studies on one of the arrhythmogenic mutants, GBT411, in both zebrafish and mouse models unveiled Dnajb6 as a novel SSS causative gene with a unique expression pattern within the subpopulation of sinus node pacemaker cells that partially overlaps with the expression of hyperpolarization activated cyclic nucleotide gated channel 4 (HCN4), supporting heterogeneity of the cardiac pacemaker cells.

Biological noise is a key determinant of the reproducibility and adaptability of cardiac pacemaking and EC coupling

Each heartbeat begins with the generation of an action potential in pacemaking cells in the sinoatrial node. This signal triggers contraction of cardiac muscle through a process termed excitation-contraction (EC) coupling. EC coupling is initiated in dyadic structures of cardiac myocytes, where ryanodine receptors in the junctional sarcoplasmic reticulum come into close apposition with clusters of CaV1.2 channels in invaginations of the sarcolemma. Cooperative activation of CaV1.2 channels within these clusters causes a local increase in intracellular Ca2+ that activates the juxtaposed ryanodine receptors. A salient feature of healthy cardiac function is the reliable and precise beat-to-beat pacemaking and amplitude of Ca2+ transients during EC coupling. In this review, we discuss recent discoveries suggesting that the exquisite reproducibility of this system emerges, paradoxically, from high variability at subcellular, cellular, and network levels. This variability is attributable to stochastic fluctuations in ion channel trafficking, clustering, and gating, as well as dyadic structure, which increase intracellular Ca2+ variance during EC coupling. Although the effects of these large, local fluctuations in function and organization are sometimes negligible at the macroscopic level owing to spatial-temporal summation within and across cells in the tissue, recent work suggests that the "noisiness" of these intracellular Ca2+ events may either enhance or counterintuitively reduce variability in a context-dependent manner. Indeed, these noisy events may represent distinct regulatory features in the tuning of cardiac contractility. Collectively, these observations support the importance of incorporating experimentally determined values of Ca2+ variance in all EC coupling models. The high reproducibility of cardiac contraction is a paradoxical outcome of high Ca2+ signaling variability at subcellular, cellular, and network levels caused by stochastic fluctuations in multiple processes in time and space. This underlying stochasticity, which counterintuitively manifests as reliable, consistent Ca2+ transients during EC coupling, also allows for rapid changes in cardiac rhythmicity and contractility in health and disease.

Disorder in Ca2+ release unit locations confers robustness but cuts flexibility of heart pacemaking

Excitation-contraction coupling kinetics is dictated by the action potential rate of sinoatrial-nodal cells. These cells generate local Ca releases (LCRs) that activate Na/Ca exchanger current, which accelerates diastolic depolarization and determines the pace. LCRs are generated by clusters of ryanodine receptors, Ca release units (CRUs), residing in the sarcoplasmic reticulum. While CRU distribution exhibits substantial heterogeneity, its functional importance remains unknown. Using numerical modeling, here we show that with a square lattice distribution of CRUs, Ca-induced-Ca-release propagation during diastolic depolarization is insufficient for pacemaking within a broad range of realistic ICaL densities. Allowing each CRU to deviate randomly from its lattice position allows sparks to propagate, as observed experimentally. As disorder increases, the CRU distribution exhibits larger empty spaces and simultaneously CRU clusters, as in Poisson clumping. Propagating within the clusters, Ca release becomes synchronized, increasing action potential rate and reviving pacemaker function of dormant/nonfiring cells. However, cells with fully disordered CRU positions could not reach low firing rates and their β-adrenergic-receptor stimulation effect was substantially decreased. Inclusion of Cav1.3, a low-voltage activation L-type Ca channel isoform into ICaL, strongly increases recruitment of CRUs to fire during diastolic depolarization, increasing robustness of pacemaking and complementing effects of CRU distribution. Thus, order/disorder in CRU locations along with Cav1.3 expression regulates pacemaker function via synchronization of CRU firing. Excessive CRU disorder and/or overexpression of Cav1.3 boosts pacemaker function in the basal state, but limits the rate range, which may contribute to heart rate range decline with age and disease.

Functional Heterogeneity of Cell Populations Increases Robustness of Pacemaker Function in a Numerical Model of the Sinoatrial Node Tissue

Each heartbeat is initiated by specialized pacemaker cells operating within the sinoatrial node (SAN). While individual cells within SAN tissue exhibit substantial heterogeneity of their electrophysiological parameters and Ca cycling, the role of this heterogeneity for cardiac pacemaker function remains mainly unknown. Here we investigated the problem numerically in a 25 × 25 square grid of connected coupled-clock Maltsev-Lakatta cell models. The tissue models were populated by cells with different degree of heterogeneity of the two key model parameters, maximum L-type Ca current conductance (g_CaL) and sarcoplasmic reticulum Ca pumping rate (P_up). Our simulations showed that in the areas of P_up-g_CaL parametric space at the edge of the system stability, where action potential (AP) firing is absent or dysrhythmic in SAN tissue models populated with identical cells, rhythmic AP firing can be rescued by populating the tissues with heterogeneous cells. This robust SAN function is synergistic with respect to heterogeneity in g_CaL and P_up and can be further strengthened by clustering of cells with similar properties. The effect of cell heterogeneity is not due to a simple summation of activity of intrinsically firing cells naturally present in heterogeneous SAN; rather AP firing cells locally and critically interact with non-firing/dormant cells. When firing cells prevail, they recruit many dormant cells to fire, strongly enhancing overall SAN function; and vice versa, prevailing dormant cells suppress AP firing in cells with intrinsic automaticity and halt SAN function. The transitions between firing and non-firing states of the system are sharp, resembling phase transitions in statistical physics. Furthermore, robust function of heterogeneous SAN tissue requires weak cell coupling, a known property of the central area of SAN where cardiac impulse emerges; stronger cell coupling reduces AP firing rate and ultimately halts SAN automaticity at the edge of stability.

Principles, mechanisms and functions of entrainment in biological oscillators

Entrainment is a phenomenon in which two oscillators interact with each other, typically through physical or chemical means, to synchronize their oscillations. This phenomenon occurs in biology to coordinate processes from the molecular to organismal scale. Biological oscillators can be entrained within a single cell, between cells or to an external input. Using six illustrative examples of entrainable biological oscillators, we discuss the distinctions between entrainment and synchrony and explore features that contribute to a system's propensity to entrain. Entrainment can either enhance or reduce the heterogeneity of oscillations within a cell population, and we provide examples and mechanisms of each case. Finally, we discuss the known functions of entrainment and discuss potential functions from an evolutionary perspective.

Implementing Biological Pacemakers: Design Criteria for Successful

Each heartbeat that pumps blood throughout the body is initiated by an electrical impulse generated in the sinoatrial node (SAN). However, a number of disease conditions can hamper the ability of the SAN's pacemaker cells to generate consistent action potentials and maintain an orderly conduction path, leading to arrhythmias. For symptomatic patients, current treatments rely on implantation of an electronic pacing device. However, complications inherent to the indwelling hardware give pause to categorical use of device therapy for a subset of populations, including pediatric patients or those with temporary pacing needs. Cellular-based biological pacemakers, derived in vitro or in situ, could function as a therapeutic alternative to current electronic pacemakers. Understanding how biological pacemakers measure up to the SAN would facilitate defining and demonstrating its advantages over current treatments. In this review, we discuss recent approaches to creating biological pacemakers and delineate design criteria to guide future progress based on insights from basic biology of the SAN. We emphasize the need for long-term efficacy in vivo via maintenance of relevant proteins, source-sink balance, a niche reflective of the native SAN microenvironment, and chronotropic competence. With a focus on such criteria, combined with delivery methods tailored for disease indications, clinical implementation will be attainable.

Continuous cardiorespiratory monitoring is a dominant source of predictive signal in machine learning for risk stratification and clinical decision support. Physiol Meas 2021;42. [PMID: 34580243 DOI: 10.1088/1361-6579/ac2130]

Beaulieu-Jones and coworkers propose a litmus test for the field of predictive analytics-performance improvements must be demonstrated to be the result of non-clinician-initiated data, otherwise, there should be caution in assuming that predictive models could improve clinical decision-making (Beaulieu-Joneset al2021). They demonstrate substantial prognostic information in unsorted physician orders made before the first midnight of hospital admission, and we are persuaded that it is fair to ask-if the physician thought of it first, what exactly is machine learning for in-patient risk stratification learning about? While we want predictive analytics to represent the leading indicators of a patient's illness, does it instead merely reflect the lagging indicators of clinicians' actions? We propose that continuous cardiorespiratory monitoring-'routine telemetry data,' in Beaulieu-Jones' terms-represents the most valuable non-clinician-initiated predictive signal present in patient data, and the value added to patient care justifies the efforts and expense required. Here, we present a clinical and a physiological point of view to support our contention.

β-Adrenergic Stimulation Synchronizes a Broad Spectrum of Action Potential Firing Rates of Cardiac Pacemaker Cells toward a Higher Population Average

The heartbeat is initiated by pacemaker cells residing in the sinoatrial node (SAN). SAN cells generate spontaneous action potentials (APs), i.e., normal automaticity. The sympathetic nervous system increases the heart rate commensurate with the cardiac output demand via stimulation of SAN β-adrenergic receptors (βAR). While SAN cells reportedly represent a highly heterogeneous cell population, the current dogma is that, in response to βAR stimulation, all cells increase their spontaneous AP firing rate in a similar fashion. The aim of the present study was to investigate the cell-to-cell variability in the responses of a large population of SAN cells. We measured the βAR responses among 166 single SAN cells isolated from 33 guinea pig hearts. In contrast to the current dogma, the SAN cell responses to βAR stimulation substantially varied. In each cell, changes in the AP cycle length were highly correlated (R2 = 0.97) with the AP cycle length before βAR stimulation. While, as expected, on average, the cells increased their pacemaker rate, greater responses were observed in cells with slower basal rates, and vice versa: cells with higher basal rates showed smaller responses, no responses, or even decreased their rate. Thus, βAR stimulation synchronized the operation of the SAN cell population toward a higher average rate, rather than uniformly shifting the rate in each cell, creating a new paradigm of βAR-driven fight-or-flight responses among individual pacemaker cells.

Structural Identification of the Pacemaker Cells and Expression of Hyperpolarization-Activated Cyclic Nucleotide-Gated (HCN) Channels in the Heart of the Wild Atlantic Cod, Gadus morhua (Linnaeus, 1758)

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are proteins that contain highly conserved functional domains and sequence motifs that are correlated with their unique biophysical activities, to regulate cardiac pacemaker activity and synaptic transmission. These pacemaker proteins have been studied in mammalian species, but little is known now about their heart distribution in lower vertebrates and c-AMP modulation. Here, we characterized the pacemaker system in the heart of the wild Atlantic cod (Gadus morhua), with respect to primary pacemaker molecular markers. Special focus is given to the structural, ultrastructural and molecular characterization of the pacemaker domain, through the expression of HCN channel genes and the immunohistochemistry of HCN isoforms, including the location of intracardiac neurons that are adjacent to the sinoatrial region of the heart. Similarly to zebrafish and mammals, these neurons are immunoreactive to ChAT, VAChT and nNOS. It has been shown that cardiac pacemaking can be modulated by sympathetic and parasympathetic pathways, and the existence of intracardiac neurons projecting back to the central nervous system provide a plausible link between them.

Adherens junction engagement regulates functional patterning of the cardiac pacemaker cell lineage

Cardiac pacemaker cells (CPCs) rhythmically initiate the electrical impulses that drive heart contraction. CPCs display the highest rate of spontaneous depolarization in the heart despite being subjected to inhibitory electrochemical conditions that should theoretically suppress their activity. While several models have been proposed to explain this apparent paradox, the actual molecular mechanisms that allow CPCs to overcome electrogenic barriers to their function remain poorly understood. Here, we have traced CPC development at single-cell resolution and uncovered a series of cytoarchitectural patterning events that are critical for proper pacemaking. Specifically, our data reveal that CPCs dynamically modulate adherens junction (AJ) engagement to control characteristics including surface area, volume, and gap junctional coupling. This allows CPCs to adopt a structural configuration that supports their overall excitability. Thus, our data have identified a direct role for local cellular mechanics in patterning critical morphological features that are necessary for CPC electrical activity.

Cellular and Molecular Mechanisms of Functional Hierarchy of Pacemaker Clusters in the Sinoatrial Node: New Insights into Sick Sinus Syndrome

The sinoatrial node (SAN), the primary pacemaker of the heart, consists of a heterogeneous population of specialized cardiac myocytes that can spontaneously produce action potentials, generating the rhythm of the heart and coordinating heart contractions. Spontaneous beating can be observed from very early embryonic stage and under a series of genetic programing, the complex heterogeneous SAN cells are formed with specific biomarker proteins and generate robust automaticity. The SAN is capable to adjust its pacemaking rate in response to environmental and autonomic changes to regulate the heart's performance and maintain physiological needs of the body. Importantly, the origin of the action potential in the SAN is not static, but rather dynamically changes according to the prevailing conditions. Changes in the heart rate are associated with a shift of the leading pacemaker location within the SAN and accompanied by alterations in P wave morphology and PQ interval on ECG. Pacemaker shift occurs in response to different interventions: neurohormonal modulation, cardiac glycosides, pharmacological agents, mechanical stretch, a change in temperature, and a change in extracellular electrolyte concentrations. It was linked with the presence of distinct anatomically and functionally defined intranodal pacemaker clusters that are responsible for the generation of the heart rhythm at different rates. Recent studies indicate that on the cellular level, different pacemaker clusters rely on a complex interplay between the calcium (referred to local subsarcolemmal Ca²⁺ releases generated by the sarcoplasmic reticulum via ryanodine receptors) and voltage (referred to sarcolemmal electrogenic proteins) components of so-called "coupled clock pacemaker system" that is used to describe a complex mechanism of SAN pacemaking. In this review, we examine the structural, functional, and molecular evidence for hierarchical pacemaker clustering within the SAN. We also demonstrate the unique molecular signatures of intranodal pacemaker clusters, highlighting their importance for physiological rhythm regulation as well as their role in the development of SAN dysfunction, also known as sick sinus syndrome.

cAMP-Dependent Signaling Restores AP Firing in Dormant SA Node Cells via Enhancement of Surface Membrane Currents and Calcium Coupling

Action potential (AP) firing rate and rhythm of sinoatrial nodal cells (SANC) are controlled by synergy between intracellular rhythmic local Ca2+ releases (LCRs) ("Ca2+ clock") and sarcolemmal electrogenic mechanisms ("membrane clock"). However, some SANC do not fire APs (dormant SANC). Prior studies have shown that β-adrenoceptor stimulation can restore AP firing in these cells. Here we tested whether this relates to improvement of synchronization of clock coupling. We characterized membrane potential, ion currents, Ca2+ dynamics, and phospholamban (PLB) phosphorylation, regulating Ca2+ pump in enzymatically isolated single guinea pig SANC prior to, during, and following β-adrenoceptor stimulation (isoproterenol) or application of cell-permeant cAMP (CPT-cAMP). Phosphorylation of PLB (Serine 16) was quantified in the same cells following Ca2+ measurement. In dormant SANC LCRs were small and disorganized at baseline, membrane potential was depolarized (-38 ± 1 mV, n = 46), and ICaL, If, and IK densities were smaller vs SANC firing APs. β-adrenoceptor stimulation or application of CPT-cAMP led to de novo spontaneous AP generation in 44 and 46% of dormant SANC, respectively. The initial response was an increase in size, rhythmicity and synchronization of LCRs, paralleled with membrane hyperpolarization and small amplitude APs (rate ∼1 Hz). During the transition to steady-state AP firing, LCR size further increased, while LCR period shortened. LCRs became more synchronized resulting in the growth of an ensemble LCR signal peaked in late diastole, culminating in AP ignition; the rate of diastolic depolarization, AP amplitude, and AP firing rate increased. ICaL, IK, and If amplitudes in dormant SANC increased in response to β-adrenoceptor stimulation. During washout, all changes reversed in order. Total PLB was higher, but the ratio of phosphorylated PLB (Serine 16) to total PLB was lower in dormant SANC. β-adrenoceptor stimulation increased this ratio in AP-firing cells. Thus, transition of dormant SANC to AP firing is linked to the increased functional coupling of membrane and Ca2+ clock proteins. The transition occurs via (i) an increase in cAMP-mediated phosphorylation of PLB accelerating Ca2+ pumping, (ii) increased spatiotemporal LCR synchronization, yielding a larger diastolic LCR ensemble signal resulting in an earlier increase in diastolic INCX; and (iii) increased current densities of If, ICaL, and IK.

Assembly of the Cardiac Pacemaking Complex: Electrogenic Principles of Sinoatrial Node Morphogenesis

Cardiac pacemaker cells located in the sinoatrial node initiate the electrical impulses that drive rhythmic contraction of the heart. The sinoatrial node accounts for only a small proportion of the total mass of the heart yet must produce a stimulus of sufficient strength to stimulate the entire volume of downstream cardiac tissue. This requires balancing a delicate set of electrical interactions both within the sinoatrial node and with the downstream working myocardium. Understanding the fundamental features of these interactions is critical for defining vulnerabilities that arise in human arrhythmic disease and may provide insight towards the design and implementation of the next generation of potential cellular-based cardiac therapeutics. Here, we discuss physiological conditions that influence electrical impulse generation and propagation in the sinoatrial node and describe developmental events that construct the tissue-level architecture that appears necessary for sinoatrial node function.

Heart rate fragmentation gives novel insights into non-autonomic mechanisms governing beat-to-beat control of the heart's rhythm

To demonstrate how heart rate fragmentation gives novel insights into non-autonomic mechanisms of beat-to-beat variability in cycle length, and predicts survival of cardiology clinic patients, over and above traditional clinical risk factors and measures of heart rate variability.

Approach: We studied 2893 patients seen by cardiologists with clinical data including 24-hour Holter monitoring. Novel measures of heart rate fragmentation alongside canonical time and frequency domain measures of heart rate variability, as well as an existing local dynamics score were calculated. A proportional hazards model was utilized to relate the results to survival.

Main results: The novel heart rate fragmentation measures were validated and characterized with respect to the effects of age, ectopy and atrial fibrillation. Correlations between parameters were determined. Critically, heart rate fragmentation results could not be accounted for by undersampling respiratory sinus arrhythmia. Increased heart rate fragmentation was associated with poorer survival (p ≪ 0.01 in the univariate model). In multivariable analyses, increased heart rate fragmentation and more abnormal local dynamics (p 0.045), along with increased clinical risk factors (age (p ≪ 0.01), tobacco use (p ≪ 0.01) and history of heart failure (p 0.019)) and lower low- to high-frequency ratio (p 0.022) were all independent predictors of 2-year mortality.

Significance: Analysis of continuous ECG data with heart rate fragmentation indices yields information regarding non-autonomic control of beat-to-beat variability in cycle length that is independent of and additive to established parameters for investigating heart rate variability, and predicts mortality in concert with measures of local dynamics, frequency content of heart rate, and clinical risk factors.

Sinoatrial Node Structure, Mechanics, Electrophysiology and the Chronotropic Response to Stretch in Rabbit and Mouse

The rhythmic electrical activity of the heart's natural pacemaker, the sinoatrial node (SAN), determines cardiac beating rate (BR). SAN electrical activity is tightly controlled by multiple factors, including tissue stretch, which may contribute to adaptation of BR to changes in venous return. In most animals, including human, there is a robust increase in BR when the SAN is stretched. However, the chronotropic response to sustained stretch differs in mouse SAN, where it causes variable responses, including decreased BR. The reasons for this species difference are unclear. They are thought to relate to dissimilarities in SAN electrophysiology (particularly action potential morphology) between mouse and other species and to how these interact with subcellular stretch-activated mechanisms. Furthermore, species-related differences in structural and mechanical properties of the SAN may influence the chronotropic response to SAN stretch. Here we assess (i) how the BR response to sustained stretch of rabbit and mouse isolated SAN relates to tissue stiffness, (ii) whether structural differences could account for observed differences in BR responsiveness to stretch, and (iii) whether pharmacological modification of mouse SAN electrophysiology alters stretch-induced chronotropy. We found disparities in the relationship between SAN stiffness and the magnitude of the chronotropic response to stretch between rabbit and mouse along with differences in SAN collagen structure, alignment, and changes with stretch. We further observed that pharmacological modification to prolong mouse SAN action potential plateau duration rectified the direction of BR changes during sustained stretch, resulting in a positive chronotropic response akin to that of other species. Overall, our results suggest that structural, mechanical, and background electrophysiological properties of the SAN influence the chronotropic response to stretch. Improved insight into the biophysical determinants of stretch effects on SAN pacemaking is essential for a comprehensive understanding of SAN regulation with important implications for studies of SAN physiology and its dysfunction, such as in the aging and fibrotic heart.

Desmin is essential for the structure and function of the sinoatrial node: implications for increased arrhythmogenesis

Our objective was to investigate the effect of desmin depletion on the structure and function of the sinoatrial pacemaker complex (SANcl) and its implication in arrhythmogenesis. Analysis of mice and humans (SANcl) indicated that the sinoatrial node exhibits high amounts of desmin, desmoplakin, N-cadherin, and β-catenin in structures we call "lateral intercalated disks" connecting myocytes side by side. Examination of the SANcl from an arrhythmogenic cardiomyopathy model, desmin-deficient (Des^-/-) mouse, by immunofluorescence, ultrastructural, and Western blot analysis showed that the number of these lateral intercalated disks was diminished. Also, electrophysiological recordings of the isolated compact sinoatrial node revealed increased pacemaker systolic potential and higher diastolic depolarization rate compared with wild-type mice. Prolonged interatrial conduction expressed as a longer P wave duration was also observed in Des^-/- mice. Upregulation of mRNA levels of both T-type Ca²⁺ current channels, Cav3.1 and Cav3.2, in the Des^-/- myocardium (1.8- and 2.3-fold, respectively) and a 1.9-fold reduction of funny hyperpolarization-activated cyclic nucleotide-gated K⁺ channel 1 could underlie these functional differences. To investigate arrhythmogenicity, electrocardiographic analysis of Des-deficient mice revealed a major increase in supraventricular and ventricular ectopic beats compared with wild-type mice. Heart rate variability analysis indicated a sympathetic predominance in Des^-/- mice, which may further contribute to arrhythmogenicity. In conclusion, our results indicate that desmin elimination leads to structural and functional abnormalities of the SANcl. These alterations may be enhanced by the sympathetic component of the cardiac autonomic nervous system, which is predominant in the desmin-deficient heart, thus leading to increased arrhythmogenesis.NEW & NOTEWORTHY The sinoatrial node exhibits high amounts of desmin and desmoplakin in structures we call "lateral intercalated disks," connecting side-by-side adjacent cardiomyocytes. These structures are diminished in desmin-deficient mouse models. Misregulation of T-type Ca²⁺ current and hyperpolarization-activated cyclic nucleotide-gated K⁺ channel 1 was proved along with prolonged interatrial conduction and cardiac autonomic nervous system dysfunction.

Complexities in cardiovascular rhythmicity: perspectives on circadian normality, ageing and disease

Biological rhythms exist in organisms at all levels of complexity, in most organs and at myriad time scales. Our own biological rhythms are driven by energy emitted by the sun, interacting via our retinas with brain stem centres, which then send out complex messages designed to synchronize the behaviour of peripheral non-light sensing organs, to ensure optimal physiological responsiveness and performance of the organism based on the time of day. Peripheral organs themselves have autonomous rhythmic behaviours that can act independently from central nervous system control but is entrainable. Dysregulation of biological rhythms either through environment or disease has far-reaching consequences on health that we are only now beginning to appreciate. In this review, we focus on cardiovascular rhythms in health, with ageing and under disease conditions.

ArrhythmoGenoPharmacoTherapy

This review is focusing on the understanding of various factors and components governing and controlling the occurrence of ventricular arrhythmias including (i) the role of various ion channel-related changes in the action potential (AP), (ii) electrocardiograms (ECGs), (iii) some important arrhythmogenic mediators of reperfusion, and pharmacological approaches to their attenuation. The transmembrane potential in myocardial cells is depending on the cellular concentrations of several ions including sodium, calcium, and potassium on both sides of the cell membrane and active or inactive stages of ion channels. The movements of Na+, K+, and Ca2+ via cell membranes produce various currents that provoke AP, determining the cardiac cycle and heart function. A specific channel has its own type of gate, and it is opening and closing under specific transmembrane voltage, ionic, or metabolic conditions. APs of sinoatrial (SA) node, atrioventricular (AV) node, and Purkinje cells determine the pacemaker activity (depolarization phase 4) of the heart, leading to the surface manifestation, registration, and evaluation of ECG waves in both animal models and humans. AP and ECG changes are key factors in arrhythmogenesis, and the analysis of these changes serve for the clarification of the mechanisms of antiarrhythmic drugs. The classification of antiarrhythmic drugs may be based on their electrophysiological properties emphasizing the connection between basic electrophysiological activities and antiarrhythmic properties. The review also summarizes some important mechanisms of ventricular arrhythmias in the ischemic/reperfused myocardium and permits an assessment of antiarrhythmic potential of drugs used for pharmacotherapy under experimental and clinical conditions.

The Cardiac Pacemaker Story-Fundamental Role of the Na+/Ca2+ Exchanger in Spontaneous Automaticity

The electrophysiological mechanism of the sinus node automaticity was previously considered exclusively regulated by the so-called "funny current". However, parallel investigations increasingly emphasized the importance of the Ca²⁺-homeostasis and Na⁺/Ca²⁺ exchanger (NCX). Recently, increasing experimental evidence, as well as insight through mechanistic in silico modeling demonstrates the crucial role of the exchanger in sinus node pacemaking. NCX had a key role in the exciting story of discovery of sinus node pacemaking mechanisms, which recently settled with a consensus on the coupled-clock mechanism after decades of debate. This review focuses on the role of the Na⁺/Ca²⁺ exchanger from the early results and concepts to recent advances and attempts to give a balanced summary of the characteristics of the local, spontaneous, and rhythmic Ca²⁺ releases, the molecular control of the NCX and its role in the fight-or-flight response. Transgenic animal models and pharmacological manipulation of intracellular Ca²⁺ concentration and/or NCX demonstrate the pivotal function of the exchanger in sinus node automaticity. We also highlight where specific hypotheses regarding NCX function have been derived from computational modeling and require experimental validation. Nonselectivity of NCX inhibitors and the complex interplay of processes involved in Ca²⁺ handling render the design and interpretation of these experiments challenging.

Ultrastructure of primary pacemaking cells in rabbit sino-atrial node cells indicates limited sarcoplasmic reticulum content

The main mammalian heart pacemakers are spindle-shaped cells compressed into tangles within protective layers of collagen in the sino-atrial node (SAN). Two cell types, "dark" and "light," differ on their high or low content of intermediate filaments, but share scarcity of myofibrils and a high content of glycogen. Sarcoplasmic reticulum (SR) is scarce. The free SR (fSR) occupies 0.04% of the cell volume within ~0.4 µm wide peripheral band. The junctional SR (jSR), constituting peripheral couplings (PCs), occupies 0.03% of the cell volume. Total fSR + jSR volume is 0.07% of cell volume, lower than the SR content of ventricular myocytes. The average distance between PCs is 7.6 µm along the periphery. On the average, 30% of the SAN cells surfaces is in close proximity to others. Identifiable gap junctions are extremely rare, but small sites of close membrane-to-membrane contacts are observed. Possibly communication occurs via these very small sites of contact if conducting channels (connexons) are located within them. There is no obvious anatomical detail that might support ephaptic coupling. These observations have implications for understanding of SAN cell physiology, and require incorporation into biophysically detailed models of SAN cell behavior that currently do not include such features.

Shift of leading pacemaker site during reflex vagal stimulation and altered electrical source-to-sink balance

KEY POINTS

Vagal reflexes slow heart rate and can change where the heartbeat originates within the sinoatrial node (SAN). The mechanisms responsible for this process - termed leading pacemaker (LP) shift - have not been investigated fully. We used optical mapping to measure the effects of baroreflex, chemoreflex and carbachol on pacemaker entrainment and electrical conduction across the SAN. All methods of stimulation triggered shifts in LP site from the central SAN to one or two caudal pacemaker regions. These shifts were associated with reduced current generation capacity centrally and increased electrical load caudally. Previous studies suggest LP shift is a rate-dependent phenomenon whereby acetylcholine slows central pacemaker rate disproportionately, enabling caudal cells that are less acetylcholine sensitive to assume control. However, our findings indicate the LP region is defined by both pacemaker rate and capacity to drive activation. Shifts in LP site provide an important homeostatic mechanism for rapid switches in heart rate.

ABSTRACT

Reflex vagal activity causes abrupt heart rate slowing with concomitant caudal shifts of the leading pacemaker (LP) site within the sinoatrial node (SAN). However, neither the mechanisms responsible nor their dynamics have been investigated fully. Therefore, the objective of this study was to elucidate the mechanisms driving cholinergic LP shift. Optical maps of right atrial activation were acquired in a rat working heart-brainstem preparation during baroreflex and chemoreflex stimulation or with carbachol. All methods of stimulation triggered shifts in LP site from the central SAN to caudal pacemaker regions, which were positive for HCN4 and received uniform cholinergic innervation. During baroreflex onset, the capacity of the central region to drive activation declined with a decrease in amplitude and gradient of optical action potentials (OAPs) in the surrounding myocardium. Accompanying this decline, there was altered entrainment in the caudal SAN as shown by decreased conduction velocity, OAP amplitude, gradient and activation time. Atropine abolished these responses. Chemoreflex stimulation produced similar effects but central capacity to drive activation was preserved before the LP shift. In contrast, carbachol produced a prolonged period of reduced capacity to drive and altered entrainment. Previous studies suggest LP shift is a rate-dependent phenomenon whereby acetylcholine slows central pacemaker rate disproportionately, enabling caudal cells that are less acetylcholine sensitive to assume control. Our findings indicate that cholinergic LP shifts are also determined by altered electrical source-to-sink balance in the SAN. We conclude that the LP region is defined by both rate and capacity to drive atrial activation.

Functional Microdomains in Heart's Pacemaker: A Step Beyond Classical Electrophysiology and Remodeling

Spontaneous beating of the sinoatrial node (SAN), the primary pacemaker of the heart, is initiated, sustained, and regulated by a complex system that integrates ion channels and transporters on the cell membrane surface (often referred to as "membrane clock") with subcellular calcium handling machinery (by parity of reasoning referred to as an intracellular "Ca²⁺ clock"). Stable, rhythmic beating of the SAN is ensured by a rigorous synchronization between these two clocks highlighted in the coupled-clock system concept of SAN timekeeping. The emerging results demonstrate that such synchronization of the complex pacemaking machinery at the cellular level depends on tightly regulated spatiotemporal signals which are restricted to precise sub-cellular microdomains and associated with discrete clusters of different ion channels, transporters, and regulatory receptors. It has recently become evident that within the microdomains, various proteins form an interacting network and work together as a part of a macromolecular signaling complex. These protein-protein interactions are tightly controlled and regulated by a variety of neurohormonal signaling pathways and the diversity of cellular responses achieved with a limited pool of second messengers is made possible through the organization of essential signal components in particular microdomains. In this review, we highlight the emerging understanding of the functionality of distinct subcellular microdomains in SAN myocytes and their functional role in the accumulation and neurohormonal regulation of proteins involved in cardiac pacemaking. We also demonstrate how changes in scaffolding proteins may lead to microdomain-targeted remodeling and regulation of pacemaker proteins contributing to SAN dysfunction.

Heterogeneity of calcium clock functions in dormant, dysrhythmically and rhythmically firing single pacemaker cells isolated from SA node

Current understanding of how cardiac pacemaker cells operate is based mainly on studies in isolated single sinoatrial node cells (SANC), specifically those that rhythmically fire action potentials similar to the in vivo behavior of the intact sinoatrial node. However, only a small fraction of SANC exhibit rhythmic firing after isolation. Other SANC behaviors have not been studied. Here, for the first time, we studied all single cells isolated from the sinoatrial node of the guinea pig, including traditionally studied rhythmically firing cells ('rhythmic SANC'), dysrhythmically firing cells ('dysrhythmic SANC') and cells without any apparent spontaneous firing activity ('dormant SANC'). Action potential-induced cytosolic Ca²⁺ transients and spontaneous local Ca²⁺ releases (LCRs) were measured with a 2D camera. LCRs were present not only in rhythmically firing SANC, but also in dormant and dysrhythmic SANC. While rhythmic SANC were characterized by large LCRs synchronized in space and time towards late diastole, dysrhythmic and dormant SANC exhibited smaller LCRs that appeared stochastically and were widely distributed in time. β-adrenergic receptor (βAR) stimulation increased LCR size and synchronized LCR occurrences in all dysrhythmic and a third of dormant cells (25 of 75 cells tested). In response to βAR stimulation, these dormant SANC developed automaticity, and LCRs became coupled to spontaneous action potential-induced cytosolic Ca²⁺ transients. Conversely, dormant SANC that did not develop automaticity showed no significant change in average LCR characteristics. The majority of dysrhythmic cells became rhythmic in response to βAR stimulation, with the rate of action potential-induced cytosolic Ca²⁺ transients substantially increasing. In summary, isolated SANC can be broadly categorized into three major populations: dormant, dysrhythmic, and rhythmic. We interpret our results based on simulations of a numerical model of SANC operating as a coupled-clock system. On this basis, the two previously unstudied dysrhythmic and dormant cell populations have intrinsically partially or completely uncoupled clocks. Such cells can be recruited to fire rhythmically in response to βAR stimulation via increased rhythmic LCR activity and ameliorated coupling between the Ca²⁺ and membrane clocks.