The "funny" current (I(f)) inhibition by ivabradine at membrane potentials encompassing spontaneous depolarization in pacemaker cells

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.

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.

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.

Improved Generation of Human Induced Pluripotent Stem Cell-Derived Cardiac Pacemaker Cells Using Novel Differentiation Protocols

Current protocols for the differentiation of human-induced pluripotent stem cells (hiPSC) into cardiomyocytes only generate a small amount of cardiac pacemaker cells. In previous work, we reported the generation of high amounts of cardiac pacemaker cells by co-culturing hiPSC with mouse visceral endoderm-like (END2) cells. However, potential medical applications of cardiac pacemaker cells generated according to this protocol, comprise an incalculable xenogeneic risk. We thus aimed to establish novel protocols maintaining the differentiation efficiency of the END2 cell-based protocol, yet eliminating the use of END2 cells. Three protocols were based on the activation and inhibition of the Wingless/Integrated (Wnt) signaling pathway, supplemented either with retinoic acid and the Wnt activator CHIR99021 (protocol B) or with the NODAL inhibitor SB431542 (protocol C) or with a combination of all three components (protocol D). An additional fourth protocol (protocol E) was used, which was originally developed by the manufacturer STEMCELL Technologies for the differentiation of hiPSC or hESC into atrial cardiomyocytes. All protocols (B, C, D, E) were compared to the END2 cell-based protocol A, serving as reference, in terms of their ability to differentiate hiPSC into cardiac pacemaker cells. Our analysis revealed that protocol E induced upregulation of 12 out of 15 cardiac pacemaker-specific genes. For comparison, reference protocol A upregulated 11, while protocols B, C and D upregulated 9, 10 and 8 cardiac pacemaker-specific genes, respectively. Cells differentiated according to protocol E displayed intense fluorescence signals of cardiac pacemaker-specific markers and showed excellent rate responsiveness to adrenergic and cholinergic stimulation. In conclusion, we characterized four novel and END2 cell-independent protocols for the differentiation of hiPSC into cardiac pacemaker cells, of which protocol E was the most efficient.

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.

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.

No evidence for pericardial restraint in the snapping turtle (Chelydra serpentina) following pharmacologically-induced bradycardia at rest or during exercise

Most animals elevate cardiac output during exercise through a rise in heart rate (f_H), whilst stroke volume (V_S) remains relatively unchanged. Cardiac pacing reveals that elevating f_H alone does not alter cardiac output, which is instead largely regulated by the peripheral vasculature. In terms of myocardial oxygen demand, an increase in f_H is more costly than that which would incur if V_S instead were to increase. We hypothesized that f_H must increase because any substantial rise in V_S would be constrained by the pericardium. To investigate this hypothesis, we explored the effects of pharmacologically-induced bradycardia, with ivabradine treatment, on V_S at rest and during exercise in the common snapping turtle (Chelydra serpentina) with intact or opened pericardium. We first showed that, in isolated myocardial preparations, ivabradine exerted a pronounced positive inotropic effect on atrial tissue, but only minor effects on ventricle. Ivabradine reduced f_H in vivo, such that exercise tachycardia was attenuated. Pulmonary and systemic V_S rose in response to ivabradine. The rise in pulmonary V_S largely compensated for the bradycardia at rest, leaving total pulmonary flow unchanged by ivabradine, although ivabradine reduced pulmonary blood flow during swimming (exercise x ivabradine interaction, P<0.05). Although systemic V_S increased, systemic blood flow was reduced by ivabradine both at rest and during exercise, in spite of ivabradine's potential to increase cardiac contractility. Opening the pericardium had no effect on f_H, V_S or blood flows before or after ivabradine, indicating that the pericardium does not constrain V_S in turtles, even during pharmacologically-induced bradycardia.

Dexmedetomidine Exerts a Negative Chronotropic Action on Sinoatrial Node Cells Through the Activation of Imidazoline Receptors

ABSTRACT

Dexmedetomidine (DEX), an α2-adrenoreceptor (α2-AR) and imidazoline receptor agonist, is most often used for the sedation of patients in the intensive care unit. Its administration is associated with an increased incidence of bradycardia; however, the precise mechanism of DEX-induced bradycardia has yet to be fully elucidated. This study was undertaken to examine whether DEX modifies pacemaker activity and the underlying ionic channel function through α2-AR and imidazoline receptors. The whole-cell patch-clamp techniques were used to record action potentials and related ionic currents of sinoatrial node cells in guinea pigs. DEX (≥10 nM) reduced sinoatrial node automaticity and the diastolic depolarization rate. DEX reduced the amplitude of hyperpolarization-activated cation current (If or Ih) the pacemaker current, even within the physiological pacemaker potential range. DEX slowed the If current activation kinetics and caused a significant shift in the voltage dependence of channel activation to negative potentials. In addition, efaroxan, an α2-AR and imidazoline I1 receptor antagonist, attenuated the inhibitory effects of DEX on sinoatrial node automaticity and If current activity, whereas yohimbine, an α2-AR-selective antagonist, did not. DEX did not affect the current activities of other channels, including rapidly and slowly activating delayed rectifier K+ currents (IKr and IKs), L-type Ca2+ current (ICa,L), Na+/Ca2+ exchange current (INCX), and muscarinic K+ current (IK,ACh). Our results indicate that DEX, at clinically relevant concentrations, induced a negative chronotropic effect on the sinoatrial node function through the downregulation of If current through an imidazoline I1 receptor other than the α2-AR in the clinical setting.

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.

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.

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.

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.

Automatic Optimization of an in Silico Model of Human iPSC Derived Cardiomyocytes Recapitulating Calcium Handling Abnormalities

The growing importance of human induced pluripotent stem cell-derived cardiomyoyctes (hiPSC-CMs), as patient-specific and disease-specific models for studying cellular cardiac electrophysiology or for preliminary cardiotoxicity tests, generated better understanding of hiPSC-CM biophysical mechanisms and great amount of action potential and calcium transient data. In this paper, we propose a new hiPSC-CM in silico model, with particular attention to Ca²⁺ handling. We used (i) the hiPSC-CM Paci2013 model as starting point, (ii) a new dataset of Ca²⁺ transient measurements to tune the parameters of the inward and outward Ca²⁺ fluxes of sarcoplasmic reticulum, and (iii) an automatic parameter optimization to fit action potentials and Ca²⁺ transients. The Paci2018 model simulates, together with the typical hiPSC-CM spontaneous action potentials, more refined Ca²⁺ transients and delayed afterdepolarizations-like abnormalities, which the old Paci2013 was not able to predict due to its mathematical formulation. The Paci2018 model was validated against (i) the same current blocking experiments used to validate the Paci2013 model, and (ii) recently published data about effects of different extracellular ionic concentrations. In conclusion, we present a new and more versatile in silico model, which will provide a platform for modeling the effects of drugs or mutations that affect Ca²⁺ handling in hiPSC-CMs.

Pharmacological analysis of the transmembrane action potential configuration in myoepithelial cells of the spontaneously beating heart of the ascidian Styela rustica in vitro

The mechanisms of action potential (AP) generation in the myoepithelial cells of the tunicate heart are not yet well understood. Here, an attempt was made to elucidate these mechanisms by analyzing the effects of specific blockers of K⁺, Na⁺ and Ca²⁺ currents on the configuration of transmembrane APs and their frequency in the spontaneously beating ascidian heart. In addition, an immunocytochemical analysis of heart myoepithelial cells was performed. Staining with anti-FMRF-amide and anti-tubulin antibodies did not reveal any nerve elements within the heart tube. Treatment with 1 mmol l^-1 TEA (I_K blocker) resulted in depolarization of heart cell sarcolemma by 10 mV, and inhibition of APs generation was recorded after 3 min of exposure. Prior to this moment, the frequency of AP generation in a burst decreased from 16-18 to 2 beats min^-1 owing to prolongation of the diastole. After application of ivabradine (3 or 10 µmol l^-1), the spontaneous APs generation frequency decreased by 24%. Based on these results and published data, it is concluded that the key role in the automaticity of the ascidian heart is played by the outward K⁺ currents, Na⁺ currents, activated hyperpolarization current I_f and a current of unknown nature I_X.

Subtype-specific differentiation of cardiac pacemaker cell clusters from human induced pluripotent stem cells

Background

Human induced pluripotent stem cells (hiPSC) harbor the potential to differentiate into diverse cardiac cell types. Previous experimental efforts were primarily directed at the generation of hiPSC-derived cells with ventricular cardiomyocyte characteristics. Aiming at a straightforward approach for pacemaker cell modeling and replacement, we sought to selectively differentiate cells with nodal-type properties.

Methods

hiPSC were differentiated into spontaneously beating clusters by co-culturing with visceral endoderm-like cells in a serum-free medium. Subsequent culturing in a specified fetal bovine serum (FBS)-enriched cell medium produced a pacemaker-type phenotype that was studied in detail using quantitative real-time polymerase chain reaction (qRT-PCR), immunocytochemistry, and patch-clamp electrophysiology. Further investigations comprised pharmacological stimulations and co-culturing with neonatal cardiomyocytes.

Results

hiPSC co-cultured in a serum-free medium with the visceral endoderm-like cell line END-2 produced spontaneously beating clusters after 10–12 days of culture. The pacemaker-specific genes HCN4, TBX3, and TBX18 were abundantly expressed at this early developmental stage, while levels of sarcomeric gene products remained low. We observed that working-type cardiomyogenic differentiation can be suppressed by transfer of early clusters into a FBS-enriched cell medium immediately after beating onset. After 6 weeks under these conditions, sinoatrial node (SAN) hallmark genes remained at high levels, while working-type myocardial transcripts (NKX2.5, TBX5) were low. Clusters were characterized by regular activity and robust beating rates (70–90 beats/min) and were triggered by spontaneous Ca²⁺ transients recapitulating calcium clock properties of genuine pacemaker cells. They were responsive to adrenergic/cholinergic stimulation and able to pace neonatal rat ventricular myocytes in co-culture experiments. Action potential (AP) measurements of cells individualized from clusters exhibited nodal-type (63.4%) and atrial-type (36.6%) AP morphologies, while ventricular AP configurations were not observed.

Conclusion

We provide a novel culture media-based, transgene-free approach for targeted generation of hiPSC-derived pacemaker-type cells that grow in clusters and offer the potential for disease modeling, drug testing, and individualized cell-based replacement therapy of the SAN.

Electronic supplementary material

The online version of this article (doi:10.1186/s13287-017-0681-4) contains supplementary material, which is available to authorized users.

Small functional If current in sinoatrial pacemaker cells of the brown trout (Salmo trutta fario) heart despite strong expression of HCN channel transcripts

Funny current (I_f), formed by hyperpolarization-activated cyclic nucleotide-gated channels (HCN channels), is supposed to be crucial for the membrane clock regulating the cardiac pacemaker mechanism. We examined the presence and activity of HCN channels in the brown trout (Salmo trutta fario) sinoatrial (SA) pacemaker cells and their putative role in heart rate (f_H) regulation. Six HCN transcripts (HCN1, HCN2a, HCN2ba, HCN2bb, HCN3, and HCN4) were expressed in the brown trout heart. The total HCN transcript abundance was 4.0 and 4.9 times higher in SA pacemaker tissue than in atrium and ventricle, respectively. In the SA pacemaker, HCN3 and HCN4 were the main isoforms representing 35.8 ± 2.7 and 25.0 ± 1.5%, respectively, of the total HCN transcripts. Only a small I_f with a mean current density of -1.2 ± 0.37 pA/pF at -140 mV was found in 4 pacemaker cells out of 16 spontaneously beating cells examined, despite the optimization of recording conditions for I_f activity. I_f was not found in any of the 24 atrial myocytes and 21 ventricular myocytes examined. HCN4 coexpressed with the MinK-related peptide 1 (MiRP1) β-subunit in CHO cells generated large I_f currents. In contrast, HCN3 (+MiRP1) failed to produce I_f in the same expression system. Cs⁺ (2 mM), which blocked 84 ± 12% of the native I_f, reversibly reduced f_H 19.2 ± 3.6% of the excised multicellular pacemaker tissue from 53 ± 5 to 44 ± 5 beats/min (P < 0.05). However, this effect was probably due to the reduction of I_Kr, which was also inhibited (63.5 ± 4.6%) by Cs⁺ These results strongly suggest that f_H regulation in the brown trout heart is largely independent on I_f.

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.

Hydroxychloroquine reduces heart rate by modulating the hyperpolarization-activated current If: Novel electrophysiological insights and therapeutic potential

BACKGROUND

Bradycardic agents are of interest for the treatment of ischemic heart disease and heart failure, as heart rate is an important determinant of myocardial oxygen consumption.

OBJECTIVES

The purpose of this study was to investigate the propensity of hydroxychloroquine (HCQ) to cause bradycardia.

METHODS

We assessed the effects of HCQ on (1) cardiac beating rate in vitro (mice); (2) the "funny" current (If) in isolated guinea pig sinoatrial node (SAN) myocytes (1, 3, 10 µM); (3) heart rate and blood pressure in vivo by acute bolus injection (rat, dose range 1-30 mg/kg), (4) blood pressure and ventricular function during feeding (mouse, 100 mg/kg/d for 2 wk, tail cuff plethysmography, anesthetized echocardiography).

RESULTS

In mouse atria, spontaneous beating rate was significantly (P < .05) reduced (by 9% ± 3% and 15% ± 2% at 3 and 10 µM HCQ, n = 7). In guinea pig isolated SAN cells, HCQ conferred a significant reduction in spontaneous action potential firing rate (17% ± 6%, 1 μM dose) and a dose-dependent reduction in If (13% ± 3% at 1 µM; 19% ± 2% at 3 µM). Effects were also observed on L-type calcium ion current (ICaL) (12% ± 4% reduction) and rapid delayed rectifier potassium current (IKr) (35% ± 4%) at 3 µM. Intravenous HCQ decreased heart rate in anesthetized rats (14.3% ± 1.1% at 15mg/kg; n = 6) without significantly reducing mean arterial blood pressure. In vivo feeding studies in mice showed no significant change in systolic blood pressure nor left ventricular function.

CONCLUSIONS

We have shown that HCQ acts as a bradycardic agent in SAN cells, in atrial preparations, and in vivo. HCQ slows the rate of spontaneous action potential firing in the SAN through multichannel inhibition, including that of If.

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.

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.

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.

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

Beneficial clinical bradycardic effects of ivabradine (IVA) have been interpreted solely on the basis of If inhibition, because IVA specifically inhibits If in sinoatrial nodal pacemaker cells (SANC). However, it has been recently hypothesized that SANC normal automaticity is regulated by crosstalk between an "M clock," the ensemble of surface membrane ion channels, and a "Ca(2+) clock," the sarcoplasmic reticulum (SR). We tested the hypothesis that crosstalk between the two clocks regulates SANC automaticity, and that indirect suppression of the Ca(2+) clock further contributes to IVA-induced bradycardia. IVA (3 μM) not only reduced If amplitude by 45 ± 6% in isolated rabbit SANC, but the IVA-induced slowing of the action potential (AP) firing rate was accompanied by reduced SR Ca(2+) load, slowed intracellular Ca(2+) cycling kinetics, and prolonged the period of spontaneous local Ca(2+) releases (LCRs) occurring during diastolic depolarization. Direct and specific inhibition of SERCA2 by cyclopiazonic acid (CPA) had effects similar to IVA on LCR period and AP cycle length. Specifically, the LCR period and AP cycle length shift toward longer times almost equally by either direct perturbations of the M clock (IVA) or the Ca(2+) clock (CPA), indicating that the LCR period reports the crosstalk between the clocks. Our numerical model simulations predict that entrainment between the two clocks that involves a reduction in INCX during diastolic depolarization is required to explain the experimentally AP firing rate reduction by IVA. In summary, our study provides new evidence that a coupled-clock system regulates normal cardiac pacemaker cell automaticity. Thus, IVA-induced bradycardia includes a suppression of both clocks within this system.

Genetic inhibition of Na+-Ca2+ exchanger current disables fight or flight sinoatrial node activity without affecting resting heart rate

RATIONALE

The sodium-calcium exchanger 1 (NCX1) is predominantly expressed in the heart and is implicated in controlling automaticity in isolated sinoatrial node (SAN) pacemaker cells, but the potential role of NCX1 in determining heart rate in vivo is unknown.

OBJECTIVE

To determine the role of Ncx1 in heart rate.

METHODS AND RESULTS

We used global myocardial and SAN-targeted conditional Ncx1 knockout (Ncx1(-/-)) mice to measure the effect of the NCX current on pacemaking activity in vivo, ex vivo, and in isolated SAN cells. We induced conditional Ncx1(-/-) using a Cre/loxP system. Unexpectedly, in vivo and ex vivo hearts and isolated SAN cells showed that basal rates in Ncx1(-/-) (retaining ≈20% of control level NCX current) and control mice were similar, suggesting that physiological NCX1 expression is not required for determining resting heart rate. However, increases in heart rate and SAN cell automaticity in response to isoproterenol or the dihydropyridine Ca(2+) channel agonist BayK8644 were significantly blunted or eliminated in Ncx1(-/-) mice, indicating that NCX1 is important for fight or flight heart rate responses. In contrast, the pacemaker current and L-type Ca(2+) currents were equivalent in control and Ncx1(-/-) SAN cells under resting and isoproterenol-stimulated conditions. Ivabradine, a pacemaker current antagonist with clinical efficacy, reduced basal SAN cell automaticity similarly in control and Ncx1(-/-) mice. However, ivabradine decreased automaticity in SAN cells isolated from Ncx1(-/-) mice more effectively than in control SAN cells after isoproterenol, suggesting that the importance of NCX current in fight or flight rate increases is enhanced after pacemaker current inhibition.

CONCLUSIONS

Physiological Ncx1 expression is required for increasing sinus rates in vivo, ex vivo, and in isolated SAN cells, but not for maintaining resting heart rate.