New evidence for coupled clock regulation of the normal automaticity of sinoatrial nodal pacemaker cells: bradycardic effects of ivabradine are linked to suppression of intracellular Ca²⁺ cycling

The role of P21-activated kinase (Pak1) in sinus node function

Sinoatrial node (SAN) dysfunction (SND) and atrial arrhythmia frequently occur simultaneously with a hazard ratio of 4.2 for new onset atrial fibrillation (AF) in SND patients. In the atrial muscle attenuated activity of p21-activated kinase 1 (Pak1) increases the risk for AF by enhancing NADPH oxidase 2 dependent production of reactive oxygen species (ROS). However, the role of Pak1 dependent ROS regulation in SAN function has not yet been determined. We hypothesize that Pak1 activity maintains SAN activity by regulating the expression of the hyperpolarization activated cyclic nucleotide gated cation channel (HCN). To determine Pak1 dependent changes in heart rate (HR) regulation we quantified the intrinsic sinus rhythm in wild type (WT) and Pak1 deficient (Pak1-/-) mice of both sexes in vivo and in isolated Langendorff perfused hearts. Pak1-/- hearts displayed an attenuated HR in vivo after autonomic blockage and in isolated hearts. The contribution of the Ca2+ clock to pacemaker activity remained unchanged, but Ivabradine (3 μM), a blocker of HCN channels that are a membrane clock component, eliminated the differences in SAN activity between WT and Pak1-/- hearts. Reduced HCN4 expression was confirmed in Pak1-/- right atria. The reduced HCN activity in Pak1-/- could be rescued by class II HDAC inhibition (LMK235), ROS scavenging (TEMPOL) or attenuation of Extracellular Signal-Regulated Kinase (ERK) 1/2 activity (SCH772984). No sex specific differences in Pak1 dependent SAN regulation were determined. Our results establish Pak1 as a class II HDAC regulator and a potential therapeutic target to attenuate SAN bradycardia and AF susceptibility.

STIM1 interacts with HCN4 channels to coordinate diastolic depolarization in the mouse Sinoatrial node

Cardiomyocytes in the sinoatrial node (SAN) are specialized to undergo spontaneous diastolic depolarization (DD) to create action potentials (AP) that serve as the origin of the heartbeat. Two cellular clocks govern DD: the membrane clock where ion channels contribute ionic conductance to create DD and the Ca ²⁺ clock where rhythmic Ca ²⁺ release from sarcoplasmic reticulum (SR) during diastole contributes pacemaking. How the membrane and Ca ²⁺ clocks interact to synchronize and drive DD is not well understood. Here, we identified stromal interaction molecule 1 (STIM1), the activator of store operated Ca ²⁺ entry (SOCE), in the P-cell cardiomyocytes of the SAN. Functional studies from STIM1 KO mice reveal dramatic changes in properties of AP and DD. Mechanistically, we show that STIM1 regulates the funny currents and HCN4 channels that are required to initiate DD and maintain sinus rhythm in mice. Taken together, our studies suggest that STIM1 acts as a sensor for both the Ca ²⁺ and membrane clocks for mouse SAN for cardiac pacemaking.

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.

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.

cAMP-PKA signaling modulates the automaticity of human iPSC-derived cardiomyocytes

Human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) have been used to screen and characterize drugs and to reveal mechanisms underlying cardiac diseases. However, before hiPSC-CMs can be used as a reliable experimental model, the physiological mechanisms underlying their normal function should be further explored. Accordingly, a major feature of hiPSC-CMs is automaticity, which is regulated by both Ca2+ and membrane clocks. To investigate the mechanisms coupling these clocks, we tested three hypotheses: (1) normal automaticity of spontaneously beating hiPSC-CMs is regulated by local Ca2+ releases (LCRs) and cAMP/PKA-dependent coupling of Ca2+ clock to M clock; (2) the LCR period indicates the level of crosstalk within the coupled-clock system; and (3) perturbing the activity of even one clock can lead to hiPSC-CM-altered automaticity due to diminished crosstalk within the coupled-clock system. By measuring the local and global Ca2+ transients, we found that the LCRs properties are correlated with the spontaneous beat interval. Changes in cAMP-dependent coupling of the Ca2+ and M clocks, caused by a pharmacological intervention that either activates the β-adrenergic or cholinergic receptor or upregulates/downregulates PKA signaling, affected LCR properties, which in turn altered hiPSC-CMs automaticity. Clocks' uncoupling by attenuating the pacemaker current If or the sarcoplasmic reticulum Ca2+ kinetics, decreased hiPSC-CMs beating rate, and prolonged the LCR period. Finally, LCR characteristics of spontaneously beating (at comparable rates) hiPSC-CMs and rabbit SAN are similar. In conclusion, hiPSC-CM automaticity is controlled by the coupled-clock system whose function is mediated by Ca2+-cAMP-PKA signaling.

Heart failure in mice induces a dysfunction of the sinus node associated with reduced CaMKII signaling

Dysfunction of the sinoatrial node (SAN), the natural heart pacemaker, is common in heart failure (HF) patients. SAN spontaneous activity relies on various ion currents in the plasma membrane (voltage clock), but intracellular Ca2+ ([Ca2+]i) release via ryanodine receptor 2 (RYR2; Ca2+ clock) plays an important synergetic role. Whereas remodeling of voltage-clock components has been revealed in HF, less is known about possible alterations to the Ca2+ clock. Here, we analyzed [Ca2+]i handling in SAN from a mouse HF model after transverse aortic constriction (TAC) and compared it with sham-operated animals. ECG data from awake animals showed slower heart rate in HF mice upon autonomic nervous system blockade, indicating intrinsic sinus node dysfunction. Confocal microscopy analyses of SAN cells within whole tissue showed slower and less frequent [Ca2+]i transients in HF. This correlated with fewer and smaller spontaneous Ca2+ sparks in HF SAN cells, which associated with lower RYR2 protein expression level and reduced phosphorylation at the CaMKII site. Moreover, PLB phosphorylation at the CaMKII site was also decreased in HF, which could lead to reduced sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) function and lower sarcoplasmic reticulum Ca2+ content, further depressing the Ca2+ clock. The inhibition of CaMKII with KN93 slowed [Ca2+]i transient rate in both groups, but this effect was smaller in HF SAN, consistent with less CaMKII activation. In conclusion, our data uncover that the mechanism of intrinsic pacemaker dysfunction in HF involves reduced CaMKII activation.

Increase in Ca2+-Activated cAMP/PKA Signaling Prevents Hydroxychloroquine-Induced Bradycardia of the Cardiac Pacemaker

Bradycardia or tachycardia are known side effects of drugs that limit their clinical use. The heart pacemaker function which control the heart rate under normal conditions is determined by coupled clock system. Thus, interfering with specific clock mechanism will affect other clock mechanisms through changes in interconnected signaling and can lead to rhythm disturbance. However, upregulation of a different clock components can compensate for this change. We focus here on hydroxychloroquine (HCQ), which has been shown effective in treating COVID-19 patients, however its bradycardic side effect limits its clinical use. We aim to decipher the mechanisms underlying the effect of HCQ on pacemaker automaticity, to identify a potential drug that will eliminate the bradycardia. We used isolated rabbit sinoatrial node (SAN) cells, human-induced pluripotent stem cell–derived cardiomyocytes (hiPSC-CMs) and mouse SAN cells residing in SAN tissue. Further, we employed SAN cell computational model to suggest mechanistic insights of the effect of HCQ on pacemaker function. HCQ increased mean spontaneous beat interval and variability in all three models in parallel to slower intracellular kinetics. The computational model suggested that HCQ affects the pacemaker (funny) current (If), L-type Ca2+ current (ICa,L), transient outward potassium (Ito) and due to changes in Ca2+ kinetics, the sodium-calcium exchanger current (INCX). Co-application of 3’-isobutylmethylxanthine (IBMX) and HCQ prevented the increase in beat interval and variability in all three experimental models. The HCQ-induced increase in rabbit and mice SAN cell and hiPSC-CM spontaneous beat interval, can be prevented by a phosphodiester inhibitor that restores automaticity due to slower intracellular Ca2+ kinetics.

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.

SARS-CoV-2 Infection Induces Ferroptosis of Sinoatrial Node Pacemaker Cells

BACKGROUND

Increasing evidence suggests that cardiac arrhythmias are frequent clinical features of coronavirus disease 2019 (COVID-19). Sinus node damage may lead to bradycardia. However, it is challenging to explore human sinoatrial node (SAN) pathophysiology due to difficulty in isolating and culturing human SAN cells. Embryonic stem cells (ESCs) can be a source to derive human SAN-like pacemaker cells for disease modeling.

METHODS

We used both a hamster model and human ESC (hESC)-derived SAN-like pacemaker cells to explore the impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection on the pacemaker cells of the heart. In the hamster model, quantitative real-time polymerase chain reaction and immunostaining were used to detect viral RNA and protein, respectively. We then created a dual knock-in SHOX2:GFP;MYH6:mCherry hESC reporter line to establish a highly efficient strategy to derive functional human SAN-like pacemaker cells, which was further characterized by single-cell RNA sequencing. Following exposure to SARS-CoV-2, quantitative real-time polymerase chain reaction, immunostaining, and RNA sequencing were used to confirm infection and determine the host response of hESC-SAN-like pacemaker cells. Finally, a high content chemical screen was performed to identify drugs that can inhibit SARS-CoV-2 infection, and block SARS-CoV-2-induced ferroptosis.

RESULTS

Viral RNA and spike protein were detected in SAN cells in the hearts of infected hamsters. We established an efficient strategy to derive from hESCs functional human SAN-like pacemaker cells, which express pacemaker markers and display SAN-like action potentials. Furthermore, SARS-CoV-2 infection causes dysfunction of human SAN-like pacemaker cells and induces ferroptosis. Two drug candidates, deferoxamine and imatinib, were identified from the high content screen, able to block SARS-CoV-2 infection and infection-associated ferroptosis.

CONCLUSIONS

Using a hamster model, we showed that primary pacemaker cells in the heart can be infected by SARS-CoV-2. Infection of hESC-derived functional SAN-like pacemaker cells demonstrates ferroptosis as a potential mechanism for causing cardiac arrhythmias in patients with COVID-19. Finally, we identified candidate drugs that can protect the SAN cells from SARS-CoV-2 infection.

The Inhibition of the Small-Conductance Ca2+-Activated Potassium Channels Decreases the Sinus Node Pacemaking during Beta-Adrenergic Activation

Sinus pacemaking is based on tight cooperation of intracellular Ca²⁺ handling and surface membrane ion channels. An important player of this synergistic crosstalk could be the small-conductance Ca²⁺-activated K⁺-channel (I_SK) that could contribute to the sinoatrial node (SAN) pacemaking driven by the intracellular Ca²⁺ changes under normal conditions and beta-adrenergic activation, however, the exact role is not fully clarified. SK2 channel expression was verified by immunoblot technique in rabbit SAN cells. Ionic currents and action potentials were measured by patch-clamp technique. The ECG R-R intervals were obtained by Langendorff-perfusion method on a rabbit heart. Apamin, a selective inhibitor of SK channels, was used during the experiments. Patch-clamp experiments revealed an apamin-sensitive current. When 100 nM apamin was applied, we found no change in the action potential nor in the ECG R-R interval. In experiments where isoproterenol was employed, apamin increased the cycle length of the SAN action potentials and enhanced the ECG R-R interval. Apamin did not amplify the cycle length variability or ECG R-R interval variability. Our data indicate that I_SK has no role under normal condition, however, it moderately contributes to the SAN automaticity under beta-adrenergic activation.

Drivers of Sinoatrial Node Automaticity in Zebrafish: Comparison With Mechanisms of Mammalian Pacemaker Function

Cardiac excitation originates in the sinoatrial node (SAN), due to the automaticity of this distinct region of the heart. SAN automaticity is the result of a gradual depolarisation of the membrane potential in diastole, driven by a coupled system of transarcolemmal ion currents and intracellular Ca²⁺ cycling. The frequency of SAN excitation determines heart rate and is under the control of extra- and intracardiac (extrinsic and intrinsic) factors, including neural inputs and responses to tissue stretch. While the structure, function, and control of the SAN have been extensively studied in mammals, and some critical aspects have been shown to be similar in zebrafish, the specific drivers of zebrafish SAN automaticity and the response of its excitation to vagal nerve stimulation and mechanical preload remain incompletely understood. As the zebrafish represents an important alternative experimental model for the study of cardiac (patho-) physiology, we sought to determine its drivers of SAN automaticity and the response to nerve stimulation and baseline stretch. Using a pharmacological approach mirroring classic mammalian experiments, along with electrical stimulation of intact cardiac vagal nerves and the application of mechanical preload to the SAN, we demonstrate that the principal components of the coupled membrane- Ca²⁺ pacemaker system that drives automaticity in mammals are also active in the zebrafish, and that the effects of extra- and intracardiac control of heart rate seen in mammals are also present. Overall, these results, combined with previously published work, support the utility of the zebrafish as a novel experimental model for studies of SAN (patho-) physiological function.

The Functional Role of Hyperpolarization Activated Current (I f) on Cardiac Pacemaking in Human vs. in the Rabbit Sinoatrial Node: A Simulation and Theoretical Study

The cardiac hyperpolarization-activated “funny” current (I_f), which contributes to sinoatrial node (SAN) pacemaking, has a more negative half-maximal activation voltage and smaller fully-activated macroscopic conductance in human than in rabbit SAN cells. The consequences of these differences for the relative roles of I_f in the two species, and for their responses to the specific bradycardic agent ivabradine at clinical doses have not been systematically explored. This study aims to address these issues, through incorporating rabbit and human I_f formulations developed by Fabbri et al. into the Severi et al. model of rabbit SAN cells. A theory was developed to correlate the effect of I_f reduction with the total inward depolarising current (I_total) during diastolic depolarization. Replacing the rabbit I_f formulation with the human one increased the pacemaking cycle length (CL) from 355 to 1,139 ms. With up to 20% I_f reduction (a level close to the inhibition of I_f by ivabradine at clinical concentrations), a modest increase (~5%) in the pacemaking CL was observed with the rabbit I_f formulation; however, the effect was doubled (~12.4%) for the human I_f formulation, even though the latter has smaller I_f density. When the action of acetylcholine (ACh, 0.1 nM) was considered, a 20% I_f reduction markedly increased the pacemaking CL by 37.5% (~27.3% reduction in the pacing rate), which is similar to the ivabradine effect at clinical concentrations. Theoretical analysis showed that the resultant increase of the pacemaking CL is inversely proportional to the magnitude of I_total during diastolic depolarization phase: a smaller I_f in the model resulted in a smaller I_total amplitude, resulting in a slower pacemaking rate; and the same reduction in I_f resulted in a more significant change of CL in the cell model with a smaller I_total. This explained the mechanism by which a low dose of ivabradine slows pacemaking rate more in humans than in the rabbit. Similar results were seen in the Fabbri et al. model of human SAN cells, suggesting our observations are model-independent. Collectively, the results of study explain why low dose ivabradine at clinically relevant concentrations acts as an effective bradycardic agent in modulating human SAN pacemaking.

Dual Activation of Phosphodiesterase 3 and 4 Regulates Basal Cardiac Pacemaker Function and Beyond

The sinoatrial (SA) node is the physiological pacemaker of the heart, and resting heart rate in humans is a well-known risk factor for cardiovascular disease and mortality. Consequently, the mechanisms of initiating and regulating the normal spontaneous SA node beating rate are of vital importance. Spontaneous firing of the SA node is generated within sinoatrial nodal cells (SANC), which is regulated by the coupled-clock pacemaker system. Normal spontaneous beating of SANC is driven by a high level of cAMP-mediated PKA-dependent protein phosphorylation, which rely on the balance between high basal cAMP production by adenylyl cyclases and high basal cAMP degradation by cyclic nucleotide phosphodiesterases (PDEs). This diverse class of enzymes includes 11 families and PDE3 and PDE4 families dominate in both the SA node and cardiac myocardium, degrading cAMP and, consequently, regulating basal cardiac pacemaker function and excitation-contraction coupling. In this review, we will demonstrate similarities between expression, distribution, and colocalization of various PDE subtypes in SANC and cardiac myocytes of different species, including humans, focusing on PDE3 and PDE4. Here, we will describe specific targets of the coupled-clock pacemaker system modulated by dual PDE3 + PDE4 activation and provide evidence that concurrent activation of PDE3 + PDE4, operating in a synergistic manner, regulates the basal cardiac pacemaker function and provides control over normal spontaneous beating of SANCs through (PDE3 + PDE4)-dependent modulation of local subsarcolemmal Ca²⁺ releases (LCRs).

Dual Activation of Phosphodiesterase 3 and 4 Regulates Basal Cardiac Pacemaker Function and Beyond

The sinoatrial (SA) node is the physiological pacemaker of the heart, and resting heart rate in humans is a well-known risk factor for cardiovascular disease and mortality. Consequently, the mechanisms of initiating and regulating the normal spontaneous SA node beating rate are of vital importance. Spontaneous firing of the SA node is generated within sinoatrial nodal cells (SANC), which is regulated by the coupled-clock pacemaker system. Normal spontaneous beating of SANC is driven by a high level of cAMP-mediated PKA-dependent protein phosphorylation, which rely on the balance between high basal cAMP production by adenylyl cyclases and high basal cAMP degradation by cyclic nucleotide phosphodiesterases (PDEs). This diverse class of enzymes includes 11 families and PDE3 and PDE4 families dominate in both the SA node and cardiac myocardium, degrading cAMP and, consequently, regulating basal cardiac pacemaker function and excitation-contraction coupling. In this review, we will demonstrate similarities between expression, distribution, and colocalization of various PDE subtypes in SANC and cardiac myocytes of different species, including humans, focusing on PDE3 and PDE4. Here, we will describe specific targets of the coupled-clock pacemaker system modulated by dual PDE3 + PDE4 activation and provide evidence that concurrent activation of PDE3 + PDE4, operating in a synergistic manner, regulates the basal cardiac pacemaker function and provides control over normal spontaneous beating of SANCs through (PDE3 + PDE4)-dependent modulation of local subsarcolemmal Ca²⁺ releases (LCRs).

Self-Similar Synchronization of Calcium and Membrane Potential Transitions During Action Potential Cycles Predict Heart Rate Across Species

OBJECTIVES

The purpose of this study was to discover regulatory universal mechanisms of normal automaticity in sinoatrial nodal (SAN) pacemaker cells that are self-similar across species.

BACKGROUND

Translation of knowledge of SAN automaticity gleaned from animal studies to human dysrhythmias (e.g., "sick sinus" syndrome [SSS]) requiring electronic pacemaker insertion has been suboptimal, largely because heart rate varies widely across species.

METHODS

Subcellular Ca²⁺ releases, whole cell action potential (AP)-induced Ca²⁺ transients, and APs were recorded in isolated mouse, guinea pig, rabbit, and human SAN cells. Ca²⁺-Vm kinetic parameters during phases of AP cycles from their ignition to recovery were quantified.

RESULTS

Although both action potential cycle lengths (APCLs) and Ca²⁺-Vm kinetic parameters during AP cycles differed across species by 10-fold, trans-species scaling of these during AP cycles and scaling of these to APCL in cells in vitro, electrocardiogram RR intervals in vivo, and body mass (BM) were self-similar (obeyed power laws) across species. Thus, APCL in vitro, heart rate in vivo, and BM of any species can be predicted by Ca²⁺-Vm kinetics during AP cycles in SAN cells measured in any single species in vitro.

CONCLUSIONS

In designing optimal heart rate to match widely different BM and energy requirements from mice to humans, nature did not "reinvent pacemaker cell wheels," but differentially scaled kinetics of gears that regulate the rates at which the "wheels spin." This discovery will facilitate the development of novel pharmacological and biological pacemakers featuring a normal, wide-range rate regulation in animal models and the translation of these to humans to target recalcitrant human SSS.

Genetic Ablation of G Protein-Gated Inwardly Rectifying K+ Channels Prevents Training-Induced Sinus Bradycardia

Background: Endurance athletes are prone to bradyarrhythmias, which in the long-term may underscore the increased incidence of pacemaker implantation reported in this population. Our previous work in rodent models has shown training-induced sinus bradycardia to be due to microRNA (miR)-mediated transcriptional remodeling of the HCN4 channel, leading to a reduction of the "funny" (I _f) current in the sinoatrial node (SAN). Objective: To test if genetic ablation of G-protein-gated inwardly rectifying potassium channel, also known as I _KACh channels prevents sinus bradycardia induced by intensive exercise training in mice. Methods: Control wild-type (WT) and mice lacking GIRK4 (Girk4 ^-/-), an integral subunit of I _KACh were assigned to trained or sedentary groups. Mice in the trained group underwent 1-h exercise swimming twice a day for 28 days, 7 days per week. We performed electrocardiogram recordings and echocardiography in both groups at baseline, during and after the training period. At training cessation, mice were euthanized and SAN tissues were isolated for patch clamp recordings in isolated SAN cells and molecular profiling by quantitative PCR (qPCR) and western blotting. Results: At swimming cessation trained WT mice presented with a significantly lower resting HR that was reversible by acute I _KACh block whereas Girk4 ^-/- mice failed to develop a training-induced sinus bradycardia. In line with HR reduction, action potential rate, density of I _f, as well as of T- and L-type Ca²⁺ currents (I _CaT and I _CaL ) were significantly reduced only in SAN cells obtained from WT-trained mice. I _f reduction in WT mice was concomitant with downregulation of HCN4 transcript and protein, attributable to increased expression of corresponding repressor microRNAs (miRs) whereas reduced I _CaL in WT mice was associated with reduced Ca_v1.3 protein levels. Strikingly, I _KACh ablation suppressed all training-induced molecular remodeling observed in WT mice. Conclusion: Genetic ablation of cardiac I _KACh in mice prevents exercise-induced sinus bradycardia by suppressing training induced remodeling of inward currents I _f, I _CaT and I _CaL due in part to the prevention of miR-mediated transcriptional remodeling of HCN4 and likely post transcriptional remodeling of Ca_v1.3. Strategies targeting cardiac I _KACh may therefore represent an alternative to pacemaker implantation for bradyarrhythmias seen in some veteran athletes.

Concomitant genetic ablation of L-type Cav1.3 (α1D) and T-type Cav3.1 (α1G) Ca2+ channels disrupts heart automaticity

Cardiac automaticity is set by pacemaker activity of the sinus node (SAN). In addition to the ubiquitously expressed cardiac voltage-gated L-type Ca_v1.2 Ca²⁺ channel isoform, pacemaker cells within the SAN and the atrioventricular node co-express voltage-gated L-type Ca_v1.3 and T-type Ca_v3.1 Ca²⁺ channels (SAN-VGCCs). The role of SAN-VGCCs in automaticity is incompletely understood. We used knockout mice carrying individual genetic ablation of Ca_v1.3 (Ca_v1.3^−/−) or Ca_v3.1 (Ca_v3.1^−/−) channels and double mutant Ca_v1.3^−/−/Ca_v3.1^−/− mice expressing only Ca_v1.2 channels. We show that concomitant loss of SAN-VGCCs prevents physiological SAN automaticity, blocks impulse conduction and compromises ventricular rhythmicity. Coexpression of SAN-VGCCs is necessary for impulse formation in the central SAN. In mice lacking SAN-VGCCs, residual pacemaker activity is predominantly generated in peripheral nodal and extranodal sites by f-channels and TTX-sensitive Na⁺ channels. In beating SAN cells, ablation of SAN-VGCCs disrupted late diastolic local intracellular Ca²⁺ release, which demonstrates an important role for these channels in supporting the sarcoplasmic reticulum based “Ca²⁺clock” mechanism during normal pacemaking. These data implicate an underappreciated role for co-expression of SAN-VGCCs in heart automaticity and define an integral role for these channels in mechanisms that control the heartbeat.

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.

Neurohumoral Control of Sinoatrial Node Activity and Heart Rate: Insight From Experimental Models and Findings From Humans

The sinoatrial node is perhaps one of the most important tissues in the entire body: it is the natural pacemaker of the heart, making it responsible for initiating each-and-every normal heartbeat. As such, its activity is heavily controlled, allowing heart rate to rapidly adapt to changes in physiological demand. Control of sinoatrial node activity, however, is complex, occurring through the autonomic nervous system and various circulating and locally released factors. In this review we discuss the coupled-clock pacemaker system and how its manipulation by neurohumoral signaling alters heart rate, considering the multitude of canonical and non-canonical agents that are known to modulate sinoatrial node activity. For each, we discuss the principal receptors involved and known intracellular signaling and protein targets, highlighting gaps in our knowledge and understanding from experimental models and human studies that represent areas for future research.

Novel Na+/Ca2+ Exchanger Inhibitor ORM-10962 Supports Coupled Function of Funny-Current and Na+/Ca2+ Exchanger in Pacemaking of Rabbit Sinus Node Tissue

Background and Purpose

The exact mechanism of spontaneous pacemaking is not fully understood. Recent results suggest tight cooperation between intracellular Ca²⁺ handling and sarcolemmal ion channels. An important player of this crosstalk is the Na⁺/Ca²⁺ exchanger (NCX), however, direct pharmacological evidence was unavailable so far because of the lack of a selective inhibitor. We investigated the role of the NCX current in pacemaking and analyzed the functional consequences of the I_f-NCX coupling by applying the novel selective NCX inhibitor ORM-10962 on the sinus node (SAN).

Experimental Approach

Currents were measured by patch-clamp, Ca²⁺-transients were monitored by fluorescent optical method in rabbit SAN cells. Action potentials (AP) were recorded from rabbit SAN tissue preparations. Mechanistic computational data were obtained using the Yaniv et al. SAN model.

Key Results

ORM-10962 (ORM) marginally reduced the SAN pacemaking cycle length with a marked increase in the diastolic Ca²⁺ level as well as the transient amplitude. The bradycardic effect of NCX inhibition was augmented when the funny-current (I_f) was previously inhibited and vice versa, the effect of I_f was augmented when the Ca²⁺ handling was suppressed.

Conclusion and Implications

We confirmed the contribution of the NCX current to cardiac pacemaking using a novel NCX inhibitor. Our experimental and modeling data support a close cooperation between I_f and NCX providing an important functional consequence: these currents together establish a strong depolarization capacity providing important safety factor for stable pacemaking. Thus, after individual inhibition of I_f or NCX, excessive bradycardia or instability cannot be expected because each of these currents may compensate for the reduction of the other providing safe and rhythmic SAN pacemaking.

Inhibitory Effects of Cyclopiazonic Acid on the Pacemaker Current in Sinoatrial Nodal Cells

OBJECTIVE

The spontaneous action potential of isolated sinoatrial node (SAN) cells is regulated by a coupled-clock system of two clocks: the calcium clock and membrane clock. However, it remains unclear whether calcium clock inhibitors have a direct effect on the membrane clock. The purpose of this study was to investigate the direct effect of cyclopiazonic acid (CPA), a selective calcium clock inhibitor, on the function of the membrane clock of SAN cells.

METHODS

at SAN cells were isolated by trypsinization and identified based on morphology and electrophysiology. I_f and HCN currents were recorded via patch clamp technique. The expression of the HCN channel protein was determined by Western blotting analysis.

RESULTS

The diastolic depolarization rate of spontaneous action potentials and the current densities of I_f were reduced by exposure to 10 μM CPA. The inhibitory effect of CPA was concentration-dependent with an IC₅₀ value of 16.3 μM and a Hill coefficient of 0.98. The effect of CPA on I_f current was also time-dependent, and the I_f current amplitude was partially restored after washout. Furthermore, the steady-state activation curve of the I_f current was shifted to a negative potential, indicating that channel activation slowed down. Finally, the protein expression of HCN4 in HEK293 cells was markedly downregulated by CPA.

CONCLUSIONS

These results indicate that the direct inhibition effect of CPA on the I_f current in SAN cells is both concentration- and time-dependent. The underlying mechanisms may involve slowing down steady-state activation and the downregulation of pacemaker channel protein expression.

The Hyperpolarization-Activated Cyclic-Nucleotide-Gated Channel Blocker Ivabradine Does Not Prevent Arrhythmias in Catecholaminergic Polymorphic Ventricular Tachycardia

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is an inherited, stressed-provoked ventricular arrhythmia. CPVT is treated by β-adrenergic receptor blockers, Na⁺ channel inhibitors, sympathetic denervation, or by implanting a defibrillator. We showed recently that blockers of SK4 Ca²⁺-activated K⁺ channels depolarize the maximal diastolic potential, reduce the heart rate, and attenuate ventricular arrhythmias in CPVT. The aim of the present study was to examine whether the pacemaker channel inhibitor, ivabradine could demonstrate anti-arrhythmic properties in CPVT like other bradycardic agents used in this disease and to compare them with those of the SK4 channel blocker, TRAM-34. The effects of ivabradine were examined on the arrhythmic beating of human induced pluripotent stem cells derived cardiomyocytes (hiPSC-CMs) from CPVT patients, on sinoatrial node (SAN) calcium transients, and on ECG measurements obtained from transgenic mice model of CPVT. Ivabradine did neither prevent the arrhythmic pacing of hiPSC-CMs derived from CPVT patients, nor preclude the aberrant SAN calcium transients. In contrast to TRAM-34, ivabradine was unable to reduce in vivo the ventricular premature complexes and ventricular tachyarrhythmias in transgenic CPVT mice. In conclusion, ivabradine does not exhibit anti-arrhythmic properties in CPVT, which indicates that this blocker cannot be used as a plausible treatment for CPVT ventricular arrhythmias.

Sinus node-like pacemaker mechanisms regulate ectopic pacemaker activity in the adult rat atrioventricular ring

In adult mammalian hearts, atrioventricular rings (AVRs) surround the atrial orifices of atrioventricular valves and are hotbed of ectopic activity in patients with focal atrial tachycardia. Experimental data offering mechanistic insights into initiation and maintenance of ectopic foci is lacking. We aimed to characterise AVRs in structurally normal rat hearts, identify arrhythmia predisposition and investigate mechanisms underlying arrhythmogenicity. Extracellular potential mapping and intracellular action potential recording techniques were used for electrophysiology, qPCR for gene and, Western blot and immunohistochemistry for protein expression. Conditions favouring ectopic foci were assessed by simulations. In right atrial preparations, sinus node (SN) was dominant and AVRs displayed 1:1 impulse conduction. Detaching SN unmasked ectopic pacemaking in AVRs and pacemaker action potentials were SN-like. Blocking pacemaker current I_f, and disrupting intracellular Ca²⁺ release, prolonged spontaneous cycle length in AVRs, indicating a role for SN-like pacemaker mechanisms. AVRs labelled positive for HCN4, and SERCA2a was comparable to SN. Pacemaking was potentiated by isoproterenol and abolished with carbachol and AVRs had abundant sympathetic nerve endings. β₂-adrenergic and M₂-muscarinic receptor mRNA and β₂-receptor protein were comparable to SN. In computer simulations of a sick SN, ectopic foci in AVR were unmasked, causing transient suppression of SN pacemaking.

Overexpression of a Neuronal Type Adenylyl Cyclase (Type 8) in Sinoatrial Node Markedly Impacts Heart Rate and Rhythm

Heart rate (HR) and HR variability (HRV), predictors of over-all organism health, are widely believed to be driven by autonomic input to the sinoatrial node (SAN), with sympathetic input increasing HR and reducing HRV. However, variability in spontaneous beating intervals in isolated SAN tissue and single SAN cells, devoid of autonomic neural input, suggests that clocks intrinsic to SAN cells may also contribute to HR and HRV in vivo. We assessed contributions of both intrinsic and autonomic neuronal input mechanisms of SAN cell function on HR and HRV via in vivo, telemetric EKG recordings. This was done in both wild type (WT) mice, and those in which adenylyl cyclase type 8 (ADCY8), a main driver of intrinsic cAMP-PKA-Ca²⁺ mediated pacemaker function, was overexpressed exclusively in the heart (TG^AC8). We hypothesized that TG^AC8 mice would: (1) manifest a more coherent pattern of HRV in vivo, i.e., a reduced HRV driven by mechanisms intrinsic to SAN cells, and less so to modulation by autonomic input and (2) utilize unique adaptations to limit sympathetic input to a heart with high levels of intrinsic cAMP-Ca²⁺ signaling. Increased adenylyl cyclase (AC) activity in TG^AC8 SAN tissue was accompanied by a marked increase in HR and a concurrent marked reduction in HRV, both in the absence or presence of dual autonomic blockade. The marked increase in intrinsic HR and coherence of HRV in TG^AC8 mice occurred in the context of: (1) reduced HR and HRV responses to β-adrenergic receptor (β-AR) stimulation; (2) increased transcription of genes and expression of proteins [β-Arrestin, G Protein-Coupled Receptor Kinase 5 (GRK5) and Clathrin Adaptor Protein (Dab2)] that desensitize β-AR signaling within SAN tissue, (3) reduced transcripts or protein levels of enzymes [dopamine beta-hydorxylase (DBH) and phenylethanolamine N-methyltransferase (PNMT)] required for catecholamine production in intrinsic cardiac adrenergic cells, and (4) substantially reduced plasma catecholamine levels. Thus, mechanisms driven by cAMP-PKA-Ca²⁺ signaling intrinsic to SAN cells underlie the marked coherence of TG^AC8 mice HRV. Adaptations to limit additional activation of AC signaling, via decreased neuronal sympathetic input, are utilized to ensure the hearts survival and prevent Ca²⁺ overload.

Positive Feedback Mechanisms among Local Ca Releases, NCX, and ICaL Ignite Pacemaker Action Potentials

Recent data suggest that cardiac pacemaker cell function is determined by numerous time-, voltage-, and Ca-dependent interactions of cell membrane electrogenic proteins (M-clock) and intracellular Ca cycling proteins (Ca-clock), forming a coupled-clock system. Many aspects of the coupled-clock system, however, remain underexplored. The key players of the system are Ca release channels (ryanodine receptors), generating local Ca releases (LCRs) from sarcoplasmic reticulum, electrogenic Na/Ca exchanger (NCX) current, and L-type Ca current (I_CaL). We combined numerical model simulations with experimental simultaneous recordings of action potentials (APs) and Ca to gain further insight into the complex interactions within the system. Our simulations revealed a positive feedback mechanism, dubbed AP ignition, which accelerates the diastolic depolarization (DD) to reach AP threshold. The ignition phase begins when LCRs begin to occur and the magnitude of inward NCX current begins to increase. The NCX current, together with funny current and T-type Ca current accelerates DD, bringing the membrane potential to I_CaL activation threshold. During the ignition phase, I_CaL-mediated Ca influx generates more LCRs via Ca-induced Ca release that further activates inward NCX current, creating a positive feedback. Simultaneous recordings of membrane potential and confocal Ca images support the model prediction of the positive feedback among LCRs and I_CaL, as diastolic LCRs begin to occur below and continue within the voltage range of I_CaL activation. The ignition phase onset (identified within the fine DD structure) begins when DD starts to notably accelerate (∼0.15 V/s) above the recording noise. Moreover, the timing of the ignition onset closely predicted the duration of each AP cycle in the basal state, in the presence of autonomic receptor stimulation, and in response to specific inhibition of either the M-clock or Ca-clock, thus indicating general importance of the new coupling mechanism for regulation of the pacemaker cell cycle duration, and ultimately the heart rate.

SOCE and STIM1 signaling in the heart: Timing and location matter

Store operated Ca²⁺ entry (SOCE) is an ancient and ubiquitous Ca²⁺ signaling pathway discovered decades ago, but the function of SOCE in human physiology is only now being revealed. The relevance of this pathway to striated muscle was solidified with the description of skeletal myopathies that result from mutations in STIM1 and Orai1, the two SOCE components. Here, we consider the evidence for STIM1 and SOCE in cardiac muscle and the sinoatrial node. We highlight recent studies revealing a role for STIM1 in cardiac growth in response to developmental and pathologic cues. We also review the role of STIM1 in the regulation of SOCE and Ca²⁺ store refilling in a non-Orai dependent manner. Finally, we discuss the importance of this pathway in ventricular cardiomyocytes where SOCE contribute to developmental growth and in pacemaker cells where SOCE likely has a fundamental to generating the cardiac rhythm.

Modernized Classification of Cardiac Antiarrhythmic Drugs

BACKGROUND

Among his major cardiac electrophysiological contributions, Miles Vaughan Williams (1918-2016) provided a classification of antiarrhythmic drugs that remains central to their clinical use.

METHODS

We survey implications of subsequent discoveries concerning sarcolemmal, sarcoplasmic reticular, and cytosolic biomolecules, developing an expanded but pragmatic classification that encompasses approved and potential antiarrhythmic drugs on this centenary of his birth.

RESULTS

We first consider the range of pharmacological targets, tracking these through to cellular electrophysiological effects. We retain the original Vaughan Williams Classes I through IV but subcategorize these divisions in light of more recent developments, including the existence of Na+ current components (for Class I), advances in autonomic (often G protein-mediated) signaling (for Class II), K+ channel subspecies (for Class III), and novel molecular targets related to Ca2+ homeostasis (for Class IV). We introduce new classes based on additional targets, including channels involved in automaticity, mechanically sensitive ion channels, connexins controlling electrotonic cell coupling, and molecules underlying longer-term signaling processes affecting structural remodeling. Inclusion of this widened range of targets and their physiological sequelae provides a framework for a modernized classification of established antiarrhythmic drugs based on their pharmacological targets. The revised classification allows for the existence of multiple drug targets/actions and for adverse, sometimes actually proarrhythmic, effects. The new scheme also aids classification of novel drugs under investigation.

CONCLUSIONS

We emerge with a modernized classification preserving the simplicity of the original Vaughan Williams framework while aiding our understanding and clinical management of cardiac arrhythmic events and facilitating future developments in this area.

Unique Ca2+-Cycling Protein Abundance and Regulation Sustains Local Ca2+ Releases and Spontaneous Firing of Rabbit Sinoatrial Node Cells

Spontaneous beating of the heart pacemaker, the sinoatrial node, is generated by sinoatrial node cells (SANC) and caused by gradual change of the membrane potential called diastolic depolarization (DD). Submembrane local Ca²⁺ releases (LCR) from sarcoplasmic reticulum (SR) occur during late DD and activate an inward Na⁺/Ca²⁺ exchange current, which accelerates the DD rate leading to earlier occurrence of an action potential. A comparison of intrinsic SR Ca²⁺ cycling revealed that, at similar physiological Ca²⁺ concentrations, LCRs are large and rhythmic in permeabilized SANC, but small and random in permeabilized ventricular myocytes (VM). Permeabilized SANC spontaneously released more Ca²⁺ from SR than VM, despite comparable SR Ca²⁺ content in both cell types. In this review we discuss specific patterns of expression and distribution of SR Ca²⁺ cycling proteins (SR Ca²⁺ ATPase (SERCA2), phospholamban (PLB) and ryanodine receptors (RyR)) in SANC and ventricular myocytes. We link ability of SANC to generate larger and rhythmic LCRs with increased abundance of SERCA2, reduced abundance of the SERCA inhibitor PLB. In addition, an increase in intracellular [Ca²⁺] increases phosphorylation of both PLB and RyR exclusively in SANC. The differences in SR Ca²⁺ cycling protein expression between SANC and VM provide insights into diverse regulation of intrinsic SR Ca²⁺ cycling that drives automaticity of SANC.

Spontaneous, local diastolic subsarcolemmal calcium releases in single, isolated guinea-pig sinoatrial nodal cells

Uptake and release calcium from the sarcoplasmic reticulum (SR) (dubbed “calcium clock”), in the form of spontaneous, rhythmic, local diastolic calcium releases (LCRs), together with voltage-sensitive ion channels (membrane clock) form a coupled system that regulates the action potential (AP) firing rate. LCRs activate Sodium/Calcium exchanger (NCX) that accelerates diastolic depolarization and thus participating in regulation of the time at which the next AP will occur. Previous studies in rabbit SA node cells (SANC) demonstrated that the basal AP cycle length (APCL) is tightly coupled to the basal LCR period (time from the prior AP-induced Ca²⁺ transient to the diastolic LCR occurrence), and that this coupling is further modulated by autonomic receptor stimulation. Although spontaneous LCRs during diastolic depolarization have been reported in SANC of various species (rabbit, cat, mouse, toad), prior studies have failed to detect LCRs in spontaneously beating SANC of guinea-pig, a species that has been traditionally used in studies of cardiac pacemaker cell function. We performed a detailed investigation of whether guinea-pig SANC generate LCRs and whether they play a similar key role in regulation of the AP firing rate. We used two different approaches, 2D high-speed camera and classical line-scan confocal imaging. Positioning the scan-line beneath sarcolemma, parallel to the long axis of the cell, we found that rhythmically beating guinea-pig SANC do, indeed, generate spontaneous, diastolic LCRs beneath the surface membrane. The average key LCR characteristics measured in confocal images in guinea-pig SANC were comparable to rabbit SANC, both in the basal state and in the presence of β-adrenergic receptor stimulation. Moreover, the relationship between the LCR period and APCL was subtended by the same linear function. Thus, LCRs in guinea-pig SANC contribute to the diastolic depolarization and APCL regulation. Our findings indicate that coupled-clock system regulation of APCL is a general, species-independent, mechanism of pacemaker cell normal automaticity. Lack of LCRs in prior studies is likely explained by technical issues, as individual LCRs are small stochastic events occurring mainly near the cell border.

Age-related pacemaker deterioration is due to impaired intracellular and membrane mechanisms: Insights from numerical modeling

Pacemaker function deteriorates in advanced age. Behar and Yaniv show that both intracellular and membrane mechanisms are responsible for age-associated pacemaker function deterioration and explain why the maximal beating rate is restored as a result of changes in sensitivity of HCN₄ to cAMP and phospholamban to PKA.

Age-related deterioration of pacemaker function has been documented in mammals, including humans. In aged isolated sinoatrial node tissues and cells, reduction in the spontaneous action potential (AP) firing rate was associated with deterioration of intracellular and membrane mechanisms; however, their relative contribution to age-associated deficient pacemaker function is not known. Interestingly, pharmacological interventions that increase posttranslation modification signaling activities can restore the basal and maximal AP firing rate, but the identities of the protein targets responsible for AP firing rate restoration are not known. Here, we developed a numerical model that simulates the function of a single mouse pacemaker cell. In addition to describing membrane and intracellular mechanisms, the model includes descriptions of autonomic receptor activation pathways and posttranslation modification signaling cascades. The numerical model shows that age-related deterioration of pacemaker function is related to impaired intracellular and membrane mechanisms: HCN₄, T-type channels, and phospholamban functions, as well as the node connecting these mechanisms, i.e., intracellular Ca²⁺ and posttranslation modification signaling. To explain the restored maximal beating rate in response to maximal phosphodiesterase (PDE) inhibition, autonomic receptor stimulation, or infused cyclic adenosine monophosphate (cAMP), the model predicts that phospholamban phosphorylation by protein kinase A (PKA) and HCN₄ sensitivity to cAMP are altered in advanced age. Moreover, alteration in PKA and cAMP sensitivity can also explain age-reduced sensitivity to PDE inhibition and autonomic receptor stimulation. Finally, the numerical model suggests two pharmacological approaches and one gene manipulation method to restore the basal beating rate of aged pacemaker cells to that of normal adult cells. In conclusion, our numerical model shows that impaired membrane and intracellular mechanisms and the nodes that couple them can lead to deteriorated pacemaker function. By increasing posttranslation modification signaling, the deteriorated basal and maximal age-associated beating rate can be restored to adult levels.

Computational analysis of the human sinus node action potential: model development and effects of mutations

KEY POINTS

We constructed a comprehensive mathematical model of the spontaneous electrical activity of a human sinoatrial node (SAN) pacemaker cell, starting from the recent Severi-DiFrancesco model of rabbit SAN cells. Our model is based on electrophysiological data from isolated human SAN pacemaker cells and closely matches the action potentials and calcium transient that were recorded experimentally. Simulated ion channelopathies explain the clinically observed changes in heart rate in corresponding mutation carriers, providing an independent qualitative validation of the model. The model shows that the modulatory role of the 'funny current' (If ) in the pacing rate of human SAN pacemaker cells is highly similar to that of rabbit SAN cells, despite its considerably lower amplitude. The model may prove useful in the design of experiments and the development of heart-rate modulating drugs.

ABSTRACT

The sinoatrial node (SAN) is the normal pacemaker of the mammalian heart.  Over several decades, a large amount of data on the ionic mechanisms underlying the spontaneous electrical activity of SAN pacemaker cells has been obtained, mostly in experiments on single cells isolated from rabbit SAN. This wealth of data has allowed the development of mathematical models of the electrical activity of rabbit SAN pacemaker cells. The present study aimed to construct a comprehensive model of the electrical activity of a human SAN pacemaker cell using recently obtained electrophysiological data from human SAN pacemaker cells.  We based our model on the recent Severi-DiFrancesco model of a rabbit SAN pacemaker cell. The action potential and calcium transient of the resulting model are close to the experimentally recorded values. The model has a much smaller 'funny current' (If ) than do rabbit cells, although its modulatory role is highly similar. Changes in pacing rate upon the implementation of mutations associated with sinus node dysfunction agree with the clinical observations. This agreement holds for both loss-of-function and gain-of-function mutations in the HCN4, SCN5A and KCNQ1 genes, underlying ion channelopathies in If , fast sodium current and slow delayed rectifier potassium current, respectively. We conclude that our human SAN cell model can be a useful tool in the design of experiments and the development of drugs that aim to modulate heart rate.

Semi-automated program for analysis of local Ca2+ spark release with application for classification of heart cell type

Local Ca²⁺ spark releases are essential to the Ca²⁺ cycling process. Thus, they play an important role in ventricular and atrial cell contraction, as well as in sinoatrial cell automaticity. Characterizing their properties in healthy cells from different regions in the heart can reveal the basic biophysical differences among these regions. We designed a semi-automatic Matlab Graphical User Interface (called Sparkalyzer) to characterize parameters of Ca²⁺ spark release from any major cardiac tissue, as recorded in line-scan mode with a confocal laser-scanning microscope. We validated the algorithm on experimental images from rabbit sinoatrial, atrial, and ventricular cells loaded with Fluo-4 AM. The program characterizes general image parameters of Ca²⁺ transients and sparks: spark duration, which indicates for how long the spark provides Ca²⁺ to the closed intracellular mechanisms (typical value: 25±1, 23±1, 26±1ms for sinoatrial, atrial, and ventricular cells, respectively); spark amplitude, which indicates the amount of Ca²⁺ released by a single spark (1.6±0.1, 1.6±0.2, 1.4±0.1F/F₀ for sinoatrial, atrial, and ventricular cells, respectively); spark length, which is the length of the Ca²⁺ wavelets fired out of a row of ryanodine receptors (5±0.1, 5±0.2, 3.4±0.3μm for sinoatrial, atrial, or ventricular cells, respectively) and number of sparks (0.14±0.02, 0.025±0.01, 0.02±0.01 for 1μm in 1s for sinoatrial, atrial, and ventricular cells, respectively). This method is reliable for Ca²⁺ spark analysis of sinoatrial, atrial, or ventricular cells. Moreover, by examining the average value of Ca²⁺ spark characteristics and their scattering around the mean, atrial, ventricular and sinoatrial cells can be differentiated.

Advances in the management of heart failure: the role of ivabradine

A high resting heart rate (≥70–75 b.p.m.) is a risk factor for patients with heart failure (HF) with reduced ejection fraction (EF), probably in the sense of accelerated atherosclerosis, with an increased morbidity and mortality. Beta-blockers not only reduce heart rate but also have negative inotropic and blood pressure-lowering effects, and therefore, in many patients, they cannot be given in the recommended dose. Ivabradine specifically inhibits the pacemaker current (funny current, I_f) of the sinoatrial node cells, resulting in therapeutic heart rate lowering without any negative inotropic and blood pressure-lowering effect. According to the European Society of Cardiology guidelines, ivabradine should be considered to reduce the risk of HF hospitalization and cardiovascular death in symptomatic patients with a reduced left ventricular EF ≤35% and sinus rhythm ≥70 b.p.m. despite treatment with an evidence-based dose of beta-blocker or a dose below the recommended dose (recommendation class “IIa” = weight of evidence/opinion is in favor of usefulness/efficacy: “should be considered”; level of evidence “B” = data derived from a single randomized clinical trial or large nonrandomized studies). Using a heart rate cutoff of ≥ 75 b.p.m., as licensed by the European Medicines Agency, treatment with ivabradine 5–7.5 mg b.i.d. reduces cardiovascular mortality by 17%, HF mortality by 39% and HF hospitalization rate by 30%. A high resting heart rate is not only a risk factor in HF with reduced EF but also at least a risk marker in HF with preserved EF, in acute HF and also in special forms of HF. In this review, we discuss the proven role of ivabradine in the validated indication “HF with reduced EF” together with interesting preliminary findings, and the potential role of ivabradine in further, specific forms of HF.

The end effector of circadian heart rate variation: the sinoatrial node pacemaker cell

Cardiovascular function is regulated by the rhythmicity of circadian, infradian and ultradian clocks. Specific time scales of different cell types drive their functions: circadian gene regulation at hours scale, activation-inactivation cycles of ion channels at millisecond scales, the heart's beating rate at hundreds of millisecond scales, and low frequency autonomic signaling at cycles of tens of seconds. Heart rate and rhythm are modulated by a hierarchical clock system: autonomic signaling from the brain releases neurotransmitters from the vagus and sympathetic nerves to the heart’s pacemaker cells and activate receptors on the cell. These receptors activating ultradian clock functions embedded within pacemaker cells include sarcoplasmic reticulum rhythmic spontaneous Ca²⁺ cycling, rhythmic ion channel current activation and inactivation, and rhythmic oscillatory mitochondria ATP production. Here we summarize the evidence that intrinsic pacemaker cell mechanisms are the end effector of the hierarchical brain-heart circadian clock system. [BMB Reports 2015; 48(12): 677-684]

The Autonomic Nervous System Regulates the Heart Rate through cAMP-PKA Dependent and Independent Coupled-Clock Pacemaker Cell Mechanisms

Sinoatrial nodal cells (SANCs) generate spontaneous action potentials (APs) that control the cardiac rate. The brain modulates SANC automaticity, via the autonomic nervous system, by stimulating membrane receptors that activate (adrenergic) or inactivate (cholinergic) adenylyl cyclase (AC). However, these opposing afferents are not simply additive. We showed that activation of adrenergic signaling increases AC-cAMP/PKA signaling, which mediates the increase in the SANC AP firing rate (i.e., positive chronotropic modulation). However, there is a limited understanding of the underlying internal pacemaker mechanisms involved in the crosstalk between cholinergic receptors and the decrease in the SANC AP firing rate (i.e., negative chronotropic modulation). We hypothesize that changes in AC-cAMP/PKA activity are crucial for mediating either decrease or increase in the AP firing rate and that the change in rate is due to both internal and membrane mechanisms. In cultured adult rabbit pacemaker cells infected with an adenovirus expressing the FRET sensor AKAR3, PKA activity and AP firing rate were tightly linked in response to either adrenergic receptor stimulation (by isoproterenol, ISO) or cholinergic stimulation (by carbachol, CCh). To identify the main molecular targets that mediate between PKA signaling and pacemaker function, we developed a mechanistic computational model. The model includes a description of autonomic-nervous receptors, post- translation signaling cascades, membrane molecules, and internal pacemaker mechanisms. Yielding results similar to those of the experiments, the model simulations faithfully reproduce the changes in AP firing rate in response to CCh or ISO or a combination of both (i.e., accentuated antagonism). Eliminating AC-cAMP-PKA signaling abolished the core effect of autonomic receptor stimulation on the AP firing rate. Specifically, disabling the phospholamban modulation of the SERCA activity resulted in a significantly reduced effect of CCh and a failure to increase the AP firing rate under ISO stimulation. Directly activating internal pacemaker mechanisms led to a similar extent of changes in the AP firing rate with respect to brain receptor stimulation. Thus, Ca²⁺ and cAMP/PKA-dependent phosphorylation limits the rate and magnitude of chronotropic changes in the spontaneous AP firing rate.

Ca(2+)/calmodulin-activated phosphodiesterase 1A is highly expressed in rabbit cardiac sinoatrial nodal cells and regulates pacemaker function

Constitutive Ca(2+)/calmodulin (CaM)-activation of adenylyl cyclases (ACs) types 1 and 8 in sinoatrial nodal cells (SANC) generates cAMP within lipid-raft-rich microdomains to initiate cAMP-protein kinase A (PKA) signaling, that regulates basal state rhythmic action potential firing of these cells. Mounting evidence in other cell types points to a balance between Ca(2+)-activated counteracting enzymes, ACs and phosphodiesterases (PDEs) within these cells. We hypothesized that the expression and activity of Ca(2+)/CaM-activated PDE Type 1A is higher in SANC than in other cardiac cell types. We found that PDE1A protein expression was 5-fold higher in sinoatrial nodal tissue than in left ventricle, and its mRNA expression was 12-fold greater in the corresponding isolated cells. PDE1 activity (nimodipine-sensitive) accounted for 39% of the total PDE activity in SANC lysates, compared to only 4% in left ventricular cardiomyocytes (LVC). Additionally, total PDE activity in SANC lysates was lowest (10%) in lipid-raft-rich and highest (76%) in lipid-raft-poor fractions (equilibrium sedimentation on a sucrose density gradient). In intact cells PDE1A immunolabeling was not localized to the cell surface membrane (structured illumination microscopy imaging), but located approximately within about 150nm inside of immunolabeling of hyperpolarization-activated cyclic nucleotide-gated potassium channels (HCN4), which reside within lipid-raft-rich microenvironments. In permeabilized SANC, in which surface membrane ion channels are not functional, nimodipine increased spontaneous SR Ca(2+) cycling. PDE1A mRNA silencing in HL-1 cells increased the spontaneous beating rate, reduced the cAMP, and increased cGMP levels in response to IBMX, a broad spectrum PDE inhibitor (detected via fluorescence resonance energy transfer microscopy). We conclude that signaling via cAMP generated by Ca(2+)/CaM-activated AC in SANC lipid raft domains is limited by cAMP degradation by Ca(2+)/CaM-activated PDE1A in non-lipid raft domains. This suggests that local gradients of [Ca(2+)]-CaM or different AC and PDE1A affinity regulate both cAMP production and its degradation, and this balance determines the intensity of Ca(2+)-AC-cAMP-PKA signaling that drives SANC pacemaker function.

Electrochemical Na+ and Ca2+ gradients drive coupled-clock regulation of automaticity of isolated rabbit sinoatrial nodal pacemaker cells

Coupling of an intracellular Ca(2+) clock to surface membrane ion channels, i.e., a "membrane clock, " via coupling of electrochemical Na(+) and Ca(2+) gradients (ENa and ECa, respectively) has been theorized to regulate sinoatrial nodal cell (SANC) normal automaticity. To test this hypothesis, we measured responses of [Na(+)]i, [Ca(2+)]i, membrane potential, action potential cycle length (APCL), and rhythm in rabbit SANCs to Na(+)/K(+) pump inhibition by the digitalis glycoside, digoxigenin (DG, 10-20 μmol/l). Initial small but significant increases in [Na(+)]i and [Ca(2+)]i and reductions in ENa and ECa in response to DG led to a small reduction in maximum diastolic potential (MDP), significantly enhanced local diastolic Ca(2+) releases (LCRs), and reduced the average APCL. As [Na(+)]i and [Ca(2+)]i continued to increase at longer times following DG exposure, further significant reductions in MDP, ENa, and ECa occurred; LCRs became significantly reduced, and APCL became progressively and significantly prolonged. This was accompanied by increased APCL variability. We also employed a coupled-clock numerical model to simulate changes in ENa and ECa simultaneously with ion currents not measured experimentally. Numerical modeling predicted that, as the ENa and ECa monotonically reduced over time in response to DG, ion currents (ICaL, ICaT, If, IKr, and IbNa) monotonically decreased. In parallel with the biphasic APCL, diastolic INCX manifested biphasic changes; initial INCX increase attributable to enhanced LCR ensemble Ca(2+) signal was followed by INCX reduction as ENCX (ENCX = 3ENa - 2ECa) decreased. Thus SANC automaticity is tightly regulated by ENa, ECa, and ENCX via a complex interplay of numerous key clock components that regulate SANC clock coupling.

Dynamics of PKA phosphorylation and gain of function in cardiac pacemaker cells: a computational model analysis

Cardiac pacemaker cell function is regulated by a coupled-clock system that integrates molecular cues on the cell-membrane surface (i.e., membrane clock) and on the sarcoplasmic reticulum (SR) (i.e., Ca(2+) clock). A recent study has shown that cotransfection of spontaneous beating cells (HEK293 cells and neonatal rat myocytes) with R524Q-mutant human hyperpolarization-activated cyclic nucleotide-gated molecules (the dominant component of funny channels) increases the funny channel's sensitivity to cAMP and leads to a decrease in spontaneous action potential (AP) cycle length (i.e., tachycardia). We hypothesize that in rabbit pacemaker cells, the same behavior is expected, and because of the coupled-clock system, the resultant steady-state decrease in AP cycle length will embody contributions from both clocks: the initial decrease in the spontaneous AP beating interval, arising from increased sensitivity of the f-channel to cAMP, will be accompanied by an increase in the adenylyl cyclase (AC)-cAMP-PKA-dependent phosphorylation activity, which will further decrease this interval. To test our hypothesis, we used the recently developed Yaniv-Lakatta pacemaker cell numerical model. This model predicts the cAMP signaling dynamics, as well as the kinetics and magnitude of protein phosphorylation in both normal and mutant pacemaker cells. We found that R524Q-mutant pacemaker cells have a shorter AP firing rate than that of wild-type cells and that gain in pacemaker function is the net effect of the R514Q mutation on the functioning of the coupled-clock system. Specifically, our results directly support the hypothesis that changes in Ca(2+)-activated AC-cAMP-PKA signaling are involved in the development of tachycardia in R524Q-mutant pacemaker cells.

Mechanisms underlying the cardiac pacemaker: the role of SK4 calcium-activated potassium channels

The proper expression and function of the cardiac pacemaker is a critical feature of heart physiology. The sinoatrial node (SAN) in human right atrium generates an electrical stimulation approximately 70 times per minute, which propagates from a conductive network to the myocardium leading to chamber contractions during the systoles. Although the SAN and other nodal conductive structures were identified more than a century ago, the mechanisms involved in the generation of cardiac automaticity remain highly debated. In this short review, we survey the current data related to the development of the human cardiac conduction system and the various mechanisms that have been proposed to underlie the pacemaker activity. We also present the human embryonic stem cell-derived cardiomyocyte system, which is used as a model for studying the pacemaker. Finally, we describe our latest characterization of the previously unrecognized role of the SK4 Ca(2+)-activated K(+) channel conductance in pacemaker cells. By exquisitely balancing the inward currents during the diastolic depolarization, the SK4 channels appear to play a crucial role in human cardiac automaticity.

Cell-specific Dynamic Clamp analysis of the role of funny If current in cardiac pacemaking

We used the Dynamic Clamp technique for i) comparative validation of conflicting computational models of the hyperpolarization-activated funny current, If, and ii) quantification of the role of If in mediating autonomic modulation of heart rate. Experimental protocols based on the injection of a real-time recalculated synthetic If current in sinoatrial rabbit cells were developed. Preliminary results of experiments mimicking the autonomic modulation of If demonstrated the need for a customization procedure to compensate for cellular heterogeneity. For this reason, we used a cell-specific approach, scaling the maximal conductance of the injected current based on the cell's spontaneous firing rate. The pacemaking rate, which was significantly reduced after application of Ivabradine, was restored by the injection of synthetic current based on the Severi-DiFrancesco formulation, while the injection of synthetic current based on the Maltsev-Lakatta formulation did not produce any significant variation. A positive virtual shift of the If activation curve, mimicking the Isoprenaline effects, led to a significant increase in pacemaking rate (+17.3 ± 6.7%, p < 0.01), although of lower magnitude than that induced by real Isoprenaline (+45.0 ± 26.1%). Similarly, a negative virtual shift of the activation curve significantly lowered the pacemaking rate (-11.8 ± 1.9%, p < 0.001), as did the application of real Acetylcholine (-20.5 ± 5.1%). The Dynamic Clamp approach, applied to the If study in cardiomyocytes for the first time and rate-adapted to manage intercellular variability, indicated that: i) the quantitative description of the If current in the Severi-DiFrancesco model accurately reproduces the effects of the real current on rabbit sinoatrial cell pacemaking rate and ii) a significant portion (50-60%) of the physiological autonomic rate modulation is due to the shift of the If activation curve.

Contributions of cardiac "funny" (f) channels and sarcoplasmic reticulum Ca2+ in regulating beating rate of mouse and guinea pig sinoatrial node

The aim of this study was to investigate the effects on spontaneous beating rate of mouse atrial preparations following selective block of cardiac "funny" (f) channels, I(f), and/or suppression of sarcoplasmic reticulum (SR) function in the absence and presence of β-adrenoceptor stimulation. ZD7288 [to block I(f)] caused a substantial reduction (222 ± 13 bpm) in beating rate from 431 ± 14 to 209 ± 14 bpm, ryanodine alone (to block SR Ca(2+) release) reduced beating rate by 105 ± 11 bpm, with subsequent addition of ZD7288 further reducing rate by 57 ± 9 bpm. Cyclopiazonic acid (CPA) alone (to inhibit Ca(2+) reuptake by the SR) reduced beating rate by 148 ± 13 bpm with subsequent addition of ZD7288 further reducing rate by 79 ± 12 bpm. In additional experiments measuring Ca(2+) transients in the SA node region using Rhod-2, effects of ivabradine and ZD7288 on rate were again substantially reduced after CPA. Effects of CPA alone on rate developed much more slowly than effects on Ca(2+) transient amplitude. ZD7288, ivabradine, and CPA reduced the slope and maximum response of the log(concentration)-response curves for effects of isoprenaline on beating rate. Very little response to isoprenaline remained after treatment with CPA followed by ZD7288. Similar changes in isoprenaline log(concentration)-response curves were seen in guinea pig preparations. These observations are consistent with a role for Ca(2+) released from the SR in regulating I(f) and therefore beating rate of SA node preparations; there appear to be additional contributions of SR-derived Ca(2+) to effects of β-adrenoceptor stimulation on beating rate that are independent of I(f).

Real-time relationship between PKA biochemical signal network dynamics and increased action potential firing rate in heart pacemaker cells: Kinetics of PKA activation in heart pacemaker cells

cAMP-PKA protein kinase is a key nodal signaling pathway that regulates a wide range of heart pacemaker cell functions. These functions are predicted to be involved in regulation of spontaneous action potential (AP) generation of these cells. Here we investigate if the kinetics and stoichiometry of increase in PKA activity match the increase in AP firing rate in response to β-adrenergic receptor (β-AR) stimulation or phosphodiesterase (PDE) inhibition, that alters the AP firing rate of heart sinoatrial pacemaker cells. In cultured adult rabbit pacemaker cells infected with an adenovirus expressing the FRET sensor AKAR3, the EC50 in response to graded increases in the intensity of β-AR stimulation (by Isoproterenol) the magnitude of the increases in PKA activity and the spontaneous AP firing rate were similar (0.4±0.1nM vs. 0.6±0.15nM, respectively). Moreover, the kinetics (t1/2) of the increases in PKA activity and spontaneous AP firing rate in response to β-AR stimulation or PDE inhibition were tightly linked. We characterized the system rate-limiting biochemical reactions by integrating these experimentally derived data into a mechanistic-computational model. Model simulations predicted that phospholamban phosphorylation is a potent target of the increase in PKA activity that links to increase in spontaneous AP firing rate. In summary, the kinetics and stoichiometry of increases in PKA activity in response to a physiological (β-AR stimulation) or pharmacological (PDE inhibitor) stimuli match those of changes in the AP firing rate. Thus Ca(2+)-cAMP/PKA-dependent phosphorylation limits the rate and magnitude of increase in spontaneous AP firing rate.

Potential effects of intrinsic heart pacemaker cell mechanisms on dysrhythmic cardiac action potential firing

The heart's regular electrical activity is initiated by specialized cardiac pacemaker cells residing in the sinoatrial node. The rate and rhythm of spontaneous action potential firing of sinoatrial node cells are regulated by stochastic mechanisms that determine the level of coupling of chemical to electrical clocks within cardiac pacemaker cells. This coupled-clock system is modulated by autonomic signaling from the brain via neurotransmitter release from the vagus and sympathetic nerves. Abnormalities in brain-heart clock connections or in any molecular clock activity within pacemaker cells lead to abnormalities in the beating rate and rhythm of the pacemaker tissue that initiates the cardiac impulse. Dysfunction of pacemaker tissue can lead to tachy-brady heart rate alternation or exit block that leads to long atrial pauses and increases susceptibility to other cardiac arrhythmia. Here we review evidence for the idea that disturbances in the intrinsic components of pacemaker cells may be implemented in arrhythmia induction in the heart.

From two competing oscillators to one coupled-clock pacemaker cell system

At the beginning of this century, debates regarding “what are the main control mechanisms that ignite the action potential (AP) in heart pacemaker cells” dominated the electrophysiology field. The original theory which prevailed for over 50 years had advocated that the ensemble of surface membrane ion channels (i.e., “M-clock”) is sufficient to ignite rhythmic APs. However, more recent experimental evidence in a variety of mammals has shown that the sarcoplasmic reticulum (SR) acts as a “Ca²⁺-clock” rhythmically discharges diastolic local Ca²⁺ releases (LCRs) beneath the cell surface membrane. LCRs activate an inward current (likely that of the Na⁺/Ca²⁺ exchanger) that prompts the surface membrane “M-clock” to ignite an AP. Theoretical and experimental evidence has mounted to indicate that this clock “crosstalk” operates on a beat-to-beat basis and determines both the AP firing rate and rhythm. Our review is focused on the evolution of experimental definition and numerical modeling of the coupled-clock concept, on how mechanisms intrinsic to pacemaker cell determine both the heart rate and rhythm, and on future directions to develop further the coupled-clock pacemaker cell concept.

Stochasticity intrinsic to coupled-clock mechanisms underlies beat-to-beat variability of spontaneous action potential firing in sinoatrial node pacemaker cells

Recent evidence indicates that the spontaneous action potential (AP) of isolated sinoatrial node cells (SANCs) is regulated by a system of stochastic mechanisms embodied within two clocks: ryanodine receptors of the "Ca(2+) clock" within the sarcoplasmic reticulum, spontaneously activate during diastole and discharge local Ca(2+) releases (LCRs) beneath the cell surface membrane; clock crosstalk occurs as LCRs activate an inward Na(+)/Ca(2+) exchanger current (INCX), which together with If and decay of K(+) channels prompts the "M clock," the ensemble of sarcolemmal-electrogenic molecules, to generate APs. Prolongation of the average LCR period accompanies prolongation of the average AP beating interval (BI). Moreover, the prolongation of the average AP BI accompanies increased AP BI variability. We hypothesized that both the average AP BI and AP BI variability are dependent upon stochasticity of clock mechanisms reported by the variability of LCR period. We perturbed the coupled-clock system by directly inhibiting the M clock by ivabradine (IVA) or the Ca(2+) clock by cyclopiazonic acid (CPA). When either clock is perturbed by IVA (3, 10 and 30 μM), which has no direct effect on Ca(2+) cycling, or CPA (0.5 and 5 μM), which has no direct effect on the M clock ion channels, the clock system failed to achieve the basal AP BI and both AP BI and AP BI variability increased. The changes in average LCR period and its variability in response to perturbations of the coupled-clock system were correlated with changes in AP beating interval and AP beating interval variability. We conclude that the stochasticity within the coupled-clock system affects and is affected by the AP BI firing rate and rhythm via modulation of the effectiveness of clock coupling.

Mechanisms of beat-to-beat regulation of cardiac pacemaker cell function by Ca²⁺ cycling dynamics

Whether intracellular Ca(2+) cycling dynamics regulate cardiac pacemaker cell function on a beat-to-beat basis remains unknown. Here we show that under physiological conditions, application of low concentrations of caffeine (2-4 mM) to isolated single rabbit sinoatrial node cells acutely reduces their spontaneous action potential cycle length (CL) and increases Ca(2+) transient amplitude for several cycles. Numerical simulations, using a modified Maltsev-Lakatta coupled-clock model, faithfully reproduced these effects, and also the effects of CL prolongation and dysrhythmic spontaneous beating (produced by cytosolic Ca(2+) buffering) and an acute CL reduction (produced by flash-induced Ca(2+) release from a caged Ca(2+) buffer), which we had reported previously. Three contemporary numerical models (including the original Maltsev-Lakatta model) failed to reproduce the experimental results. In our proposed new model, Ca(2+) releases acutely change the CL via activation of the Na(+)/Ca(2+) exchanger current. Time-dependent CL reductions after flash-induced Ca(2+) releases (the memory effect) are linked to changes in Ca(2+) available for pumping into sarcoplasmic reticulum which, in turn, changes the sarcoplasmic reticulum Ca(2+) load, diastolic Ca(2+) releases, and Na(+)/Ca(2+) exchanger current. These results support the idea that Ca(2+) regulates CL in cardiac pacemaker cells on a beat-to-beat basis, and suggest a more realistic numerical mechanism of this regulation.

Modern perspectives on numerical modeling of cardiac pacemaker cell

Cardiac pacemaking is a complex phenomenon that is still not completely understood. Together with experimental studies, numerical modeling has been traditionally used to acquire mechanistic insights in this research area. This review summarizes the present state of numerical modeling of the cardiac pacemaker, including approaches to resolve present paradoxes and controversies. Specifically we discuss the requirement for realistic modeling to consider symmetrical importance of both intracellular and cell membrane processes (within a recent "coupled-clock" theory). Promising future developments of the complex pacemaker system models include the introduction of local calcium control, mitochondria function, and biochemical regulation of protein phosphorylation and cAMP production. Modern numerical and theoretical methods such as multi-parameter sensitivity analyses within extended populations of models and bifurcation analyses are also important for the definition of the most realistic parameters that describe a robust, yet simultaneously flexible operation of the coupled-clock pacemaker cell system. The systems approach to exploring cardiac pacemaker function will guide development of new therapies such as biological pacemakers for treating insufficient cardiac pacemaker function that becomes especially prevalent with advancing age.

Numerical Modeling Calcium and CaMKII Effects in the SA Node

Sinoatrial node (SAN) is the primary heart pacemaker which initiates each heartbeat under normal conditions. Numerous experimental data have demonstrated that Ca(2+-) and CaMKII-dependent processes are crucially important for regulation of SAN cells. However, specific mechanisms of this regulation and their relative contribution to pacemaker function remain mainly unknown. Our review summarizes available data and existing numerical modeling approaches to understand Ca(2+) and CaMKII effects on the SAN. Data interpretation and future directions to address the problem are given within the coupled-clock theory, i.e., a modern view on the cardiac pacemaker cell function generated by a system of sarcolemmal and intracellular proteins.