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

Pacemaker Channels and the Chronotropic Response in Health and Disease

Loss or dysregulation of the normally precise control of heart rate via the autonomic nervous system plays a critical role during the development and progression of cardiovascular disease-including ischemic heart disease, heart failure, and arrhythmias. While the clinical significance of regulating changes in heart rate, known as the chronotropic effect, is undeniable, the mechanisms controlling these changes remain not fully understood. Heart rate acceleration and deceleration are mediated by increasing or decreasing the spontaneous firing rate of pacemaker cells in the sinoatrial node. During the transition from rest to activity, sympathetic neurons stimulate these cells by activating β-adrenergic receptors and increasing intracellular cyclic adenosine monophosphate. The same signal transduction pathway is targeted by positive chronotropic drugs such as norepinephrine and dobutamine, which are used in the treatment of cardiogenic shock and severe heart failure. The cyclic adenosine monophosphate-sensitive hyperpolarization-activated current (If) in pacemaker cells is passed by hyperpolarization-activated cyclic nucleotide-gated cation channels and is critical for generating the autonomous heartbeat. In addition, this current has been suggested to play a central role in the chronotropic effect. Recent studies demonstrate that cyclic adenosine monophosphate-dependent regulation of HCN4 (hyperpolarization-activated cyclic nucleotide-gated cation channel isoform 4) acts to stabilize the heart rate, particularly during rapid rate transitions induced by the autonomic nervous system. The mechanism is based on creating a balance between firing and recently discovered nonfiring pacemaker cells in the sinoatrial node. In this way, hyperpolarization-activated cyclic nucleotide-gated cation channels may protect the heart from sinoatrial node dysfunction, secondary arrhythmia of the atria, and potentially fatal tachyarrhythmia of the ventricles. Here, we review the latest findings on sinoatrial node automaticity and discuss the physiological and pathophysiological role of HCN pacemaker channels in the chronotropic response and beyond.

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.

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.

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.

Paradigm shift: new concepts for HCN4 function in cardiac pacemaking

Hyperpolarization-activated cyclic nucleotide–gated (HCN) channels are the molecular correlate of the I_f current and are critically involved in controlling neuronal excitability and the autonomous rhythm of the heart. The HCN4 isoform is the main HCN channel subtype expressed in the sinoatrial node (SAN), a tissue composed of specialized pacemaker cells responsible for generating the intrinsic heartbeat. More than 40 years ago, the I_f current was first discovered in rabbit SAN tissue. Along with this discovery, a theory was proposed that cyclic adenosine monophosphate–dependent modulation of I_f mediates heart rate regulation by the autonomic nervous system—a process called chronotropic effect. However, up to the present day, this classical theory could not be reliably validated. Recently, new concepts emerged confirming that HCN4 channels indeed play an important role in heart rate regulation. However, the cellular mechanism by which HCN4 controls heart rate turned out to be completely different than originally postulated. Here, we review the latest findings regarding the physiological role of HCN4 in the SAN. We describe a newly discovered mechanism underlying heart rate regulation by HCN4 at the tissue and single cell levels, and we discuss these observations in the context of results from previously studied HCN4 mouse models.

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.

Frequency-Dependent Properties of the Hyperpolarization-Activated Cation Current, If, in Adult Mouse Heart Primary Pacemaker Myocytes

A number of distinct electrophysiological mechanisms that modulate the myogenic spontaneous pacemaker activity in the sinoatrial node (SAN) of the mammalian heart have been investigated extensively. There is agreement that several (3 or 4) different transmembrane ionic current changes (referred to as the voltage clock) are involved; and that the resulting net current interacts with direct and indirect effects of changes in intracellular Ca²⁺ (the calcium clock). However, significant uncertainties, and important knowledge gaps, remain concerning the functional roles in SAN spontaneous pacing of many of the individual ion channel- or exchanger-mediated transmembrane current changes. We report results from patch clamp studies and mathematical modeling of the hyperpolarization-activated current, I_f, in the generation/modulation of the diastolic depolarization, or pacemaker potential, produced by individual myocytes that were enzymatically isolated from the adult mouse sinoatrial node (SAN). Amphotericin-mediated patch microelectrode recordings at 35 °C were made under control conditions and in the presence of 5 or 10 nM isoproterenol (ISO). These sets of results were complemented and integrated with mathematical modeling of the current changes that take place in the range of membrane potentials (−70 to −50 mV), which corresponds to the ‘pacemaker depolarization’ in the adult mouse SAN. Our results reveal a very small, but functionally important, approximately steady-state or time-independent current generated by residual activation of I_f channels that are expressed in these pacemaker myocytes. Recordings of the pacemaker depolarization and action potential, combined with measurements of changes in I_f, and the well-known increases in the L-type Ca²⁺ current, I_CaL, demonstrated that I_CaL activation, is essential for myogenic pacing. Moreover, after being enhanced (approximately 3-fold) by 5 or 10 nM ISO, I_CaL contributes significantly to the positive chronotropic effect. Our mathematical model has been developed in an attempt to better understand the underlying mechanisms for the pacemaker depolarization and action potential in adult mouse SAN myocytes. After being updated with our new experimental data describing I_f, our simulations reveal a novel functional component of I_f in adult mouse SAN. Computational work carried out with this model also confirms that in the presence of ISO the residual activation of I_f and opening of I_CaL channels combine to generate a net current change during the slow diastolic depolarization phase that is essential for the observed accelerated pacemaking rate of these SAN myocytes.

L-Type Cav1.3 Calcium Channels Are Required for Beta-Adrenergic Triggered Automaticity in Dormant Mouse Sinoatrial Pacemaker Cells

Background: Sinoatrial node cells (SANC) automaticity is generated by functional association between the activity of plasmalemmal ion channels and local diastolic intracellular Ca2+ release (LCR) from ryanodine receptors. Strikingly, most isolated SANC exhibit a “dormant” state, whereas only a fraction shows regular firing as observed in intact SAN. Recent studies showed that β-adrenergic stimulation can initiate spontaneous firing in dormant SANC, though this mechanism is not entirely understood. Methods: To investigate the role of L-type Cav1.3 Ca2+ channels in the adrenergic regulation of automaticity in dormant SANC, we used a knock-in mouse strain in which the sensitivity of L-type Cav1.2 α1 subunits to dihydropyridines (DHPs) was inactivated (Cav1.2DHP−/−), enabling the selective pharmacological inhibition of Cav1.3 by DHPs. Results: In dormant SANC, β-adrenergic stimulation with isoproterenol (ISO) induced spontaneous action potentials (AP) and Ca2+ transients, which were completely arrested with concomitant perfusion of the DHP nifedipine. In spontaneously firing SANC at baseline, Cav1.3 inhibition completely reversed the effect of β-adrenergic stimulation on AP and the frequency of Ca2+ transients. Confocal calcium imaging of SANC showed that the β-adrenergic-induced synchronization of LCRs is regulated by the activity of Cav1.3 channels. Conclusions: Our study shows a novel role of Cav1.3 channels in initiating and maintaining automaticity in dormant SANC upon β-adrenergic stimulation.

Mechanisms of Sinoatrial Node Dysfunction in Heart Failure With Preserved Ejection Fraction

BACKGROUND

The ability to increase heart rate during exercise and other stressors is a key homeostatic feature of the sinoatrial node (SAN). When the physiological heart rate response is blunted, chronotropic incompetence limits exercise capacity, a common problem in patients with heart failure with preserved ejection fraction (HFpEF). Despite its clinical relevance, the mechanisms of chronotropic incompetence remain unknown.

METHODS

Dahl salt-sensitive rats fed a high-salt diet and C57Bl6 mice fed a high-fat diet and an inhibitor of constitutive nitric oxide synthase (Nω-nitro-L-arginine methyl ester [L-NAME]; 2-hit) were used as models of HFpEF. Myocardial infarction was created to induce HF with reduced ejection fraction. Rats and mice fed with a normal diet or those that had a sham surgery served as respective controls. A comprehensive characterization of SAN function and chronotropic response was conducted by in vivo, ex vivo, and single-cell electrophysiologic studies. RNA sequencing of SAN was performed to identify transcriptomic changes. Computational modeling of biophysically-detailed human HFpEF SAN was created.

RESULTS

Rats with phenotypically-verified HFpEF exhibited limited chronotropic response associated with intrinsic SAN dysfunction, including impaired β-adrenergic responsiveness and an alternating leading pacemaker within the SAN. Prolonged SAN recovery time and reduced SAN sensitivity to isoproterenol were confirmed in the 2-hit mouse model. Adenosine challenge unmasked conduction blocks within the SAN, which were associated with structural remodeling. Chronotropic incompetence and SAN dysfunction were also found in rats with HF with reduced ejection fraction. Single-cell studies and transcriptomic profiling revealed HFpEF-related alterations in both the "membrane clock" (ion channels) and the "Ca2+ clock" (spontaneous Ca2+ release events). The physiologic impairments were reproduced in silico by empirically-constrained quantitative modeling of human SAN function.

CONCLUSIONS

Chronotropic incompetence and SAN dysfunction were seen in both models of HF. We identified that intrinsic abnormalities of SAN structure and function underlie the chronotropic response in HFpEF.

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.

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.

Synchronized Cardiac Impulses Emerge From Heterogeneous Local Calcium Signals Within and Among Cells of Pacemaker Tissue

OBJECTIVES

This study sought to identify subcellular Ca²⁺ signals within and among cells comprising the sinoatrial node (SAN) tissue.

BACKGROUND

The current paradigm of SAN impulse generation: 1) is that full-scale action potentials (APs) of a common frequency are initiated at 1 site and are conducted within the SAN along smooth isochrones; and 2) does not feature fine details of Ca²⁺ signaling present in isolated SAN cells, in which small subcellular, subthreshold local Ca²⁺ releases (LCRs) self-organize to generate cell-wide APs.

METHODS

Immunolabeling was combined with a novel technique to detect the occurrence of LCRs and AP-induced Ca²⁺ transients (APCTs) in individual pixels (chronopix) across the entire mouse SAN images.

RESULTS

At high magnification, Ca²⁺ signals appeared markedly heterogeneous in space, amplitude, frequency, and phase among cells comprising an HCN4⁺/CX43^- cell meshwork. The signaling exhibited several distinguishable patterns of LCR/APCT interactions within and among cells. Rhythmic APCTs that were apparently conducted within the meshwork were transferred to a truly conducting HCN4^-/CX43⁺ network of striated cells via narrow functional interfaces where different cell types intertwine, that is, the SAN anatomic/functional unit. At low magnification, the earliest APCT of each cycle occurred within a small area of the HCN4 meshwork, and subsequent APCT appearance throughout SAN pixels was discontinuous and asynchronous.

CONCLUSIONS

The study has discovered a novel, microscopic Ca²⁺ signaling paradigm of SAN operation that has escaped detection using low-resolution, macroscopic tissue isochrones employed in prior studies: synchronized APs emerge from heterogeneous subcellular subthreshold Ca²⁺ signals, resembling multiscale complex processes of impulse generation within clusters of neurons in neuronal networks.

Intracellular Na+ Modulates Pacemaking Activity in Murine Sinoatrial Node Myocytes: An In Silico Analysis

Background: The mechanisms underlying dysfunction in the sinoatrial node (SAN), the heart's primary pacemaker, are incompletely understood. Electrical and Ca2+-handling remodeling have been implicated in SAN dysfunction associated with heart failure, aging, and diabetes. Cardiomyocyte [Na+]i is also elevated in these diseases, where it contributes to arrhythmogenesis. Here, we sought to investigate the largely unexplored role of Na+ homeostasis in SAN pacemaking and test whether [Na+]i dysregulation may contribute to SAN dysfunction. Methods: We developed a dataset-specific computational model of the murine SAN myocyte and simulated alterations in the major processes of Na+ entry (Na+/Ca2+ exchanger, NCX) and removal (Na+/K+ ATPase, NKA). Results: We found that changes in intracellular Na+ homeostatic processes dynamically regulate SAN electrophysiology. Mild reductions in NKA and NCX function increase myocyte firing rate, whereas a stronger reduction causes bursting activity and loss of automaticity. These pathologic phenotypes mimic those observed experimentally in NCX- and ankyrin-B-deficient mice due to altered feedback between the Ca2+ and membrane potential clocks underlying SAN firing. Conclusions: Our study generates new testable predictions and insight linking Na+ homeostasis to Ca2+ handling and membrane potential dynamics in SAN myocytes that may advance our understanding of SAN (dys)function.

Rad-GTPase contributes to heart rate via L-type calcium channel regulation

Sinoatrial node cardiomyocytes (SANcm) possess automatic, rhythmic electrical activity. SAN rate is influenced by autonomic nervous system input, including sympathetic nerve increases of heart rate (HR) via activation of β-adrenergic receptor signaling cascade (β-AR). L-type calcium channel (LTCC) activity contributes to membrane depolarization and is a central target of β-AR signaling. Recent studies revealed that the small G-protein Rad plays a central role in β-adrenergic receptor directed modulation of LTCC. These studies have identified a conserved mechanism in which β-AR stimulation results in PKA-dependent Rad phosphorylation: depletion of Rad from the LTCC complex, which is proposed to relieve the constitutive inhibition of CaV1.2 imposed by Rad association. Here, using a transgenic mouse model permitting conditional cardiomyocyte selective Rad ablation, we examine the contribution of Rad to the control of SANcm LTCC current (ICa,L) and sinus rhythm. Single cell analysis from a recent published database indicates that Rad is expressed in SANcm, and we show that SANcm ICa,L was significantly increased in dispersed SANcm following Rad silencing compared to those from CTRL hearts. Moreover, cRadKO SANcm ICa,L was not further increased with β-AR agonists. We also evaluated heart rhythm in vivo using radiotelemetered ECG recordings in ambulating mice. In vivo, intrinsic HR is significantly elevated in cRadKO. During the sleep phase cRadKO also show elevated HR, and during the active phase there is no significant difference. Rad-deletion had no significant effect on heart rate variability. These results are consistent with Rad governing LTCC function under relatively low sympathetic drive conditions to contribute to slower HR during the diurnal sleep phase HR. In the absence of Rad, the tonic modulated SANcm ICa,L promotes elevated sinus HR. Future novel therapeutics for bradycardia targeting Rad - LTCC can thus elevate HR while retaining βAR responsiveness.

Physiological Roles of the Rapidly Activated Delayed Rectifier K+ Current in Adult Mouse Heart Primary Pacemaker Activity

Robust, spontaneous pacemaker activity originating in the sinoatrial node (SAN) of the heart is essential for cardiovascular function. Anatomical, electrophysiological, and molecular methods as well as mathematical modeling approaches have quite thoroughly characterized the transmembrane fluxes of Na⁺, K⁺ and Ca²⁺ that produce SAN action potentials (AP) and ‘pacemaker depolarizations’ in a number of different in vitro adult mammalian heart preparations. Possible ionic mechanisms that are responsible for SAN primary pacemaker activity are described in terms of: (i) a Ca²⁺-regulated mechanism based on a requirement for phasic release of Ca²⁺ from intracellular stores and activation of an inward current-mediated by Na⁺/Ca²⁺ exchange; (ii) time- and voltage-dependent activation of Na⁺ or Ca²⁺ currents, as well as a cyclic nucleotide-activated current, I_f; and/or (iii) a combination of (i) and (ii). Electrophysiological studies of single spontaneously active SAN myocytes in both adult mouse and rabbit hearts consistently reveal significant expression of a rapidly activating time- and voltage-dependent K⁺ current, often denoted I_Kr, that is selectively expressed in the leading or primary pacemaker region of the adult mouse SAN. The main goal of the present study was to examine by combined experimental and simulation approaches the functional or physiological roles of this K⁺ current in the pacemaker activity. Our patch clamp data of mouse SAN myocytes on the effects of a pharmacological blocker, E4031, revealed that a rapidly activating K⁺ current is essential for action potential (AP) repolarization, and its deactivation during the pacemaker potential contributes a small but significant component to the pacemaker depolarization. Mathematical simulations using a murine SAN AP model confirm that well known biophysical properties of a delayed rectifier K⁺ current can contribute to its role in generating spontaneous myogenic 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.

Beyond pacemaking: HCN channels in sinoatrial node function

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are key proteins involved in the initiation and regulation of the heartbeat. Pacemaker cells within the sinoatrial node generate the electrical impulse that underlies the contraction of all atrial and ventricular cardiomyocytes. To generate a stable heart rhythm, it is necessary that the spontaneous activity of pacemaker cells is synchronized. Entrainment processes in the sinoatrial node create synchrony and also mediate heart rate regulation. In the past years it has become clear that the role of HCN channels goes beyond just pacemaking and that the channels play pivotal roles in these entrainment processes that coordinate and balance sinoatrial node network activity. Here, we review the role of HCN channels in the central pacemaker process and highlight new aspects of the contribution of HCN channels to stabilizing the electrical activity of the sinoatrial node network, especially during heart rate regulation by the autonomic nervous system.

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.

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.

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.