Relative importance of funny current in human versus rabbit sinoatrial node

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

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

The 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.

HCN4 current during human sinoatrial node-like action potentials

BACKGROUND

Despite the many studies carried out over the past 40 years, the contribution of the HCN4 encoded hyperpolarization-activated 'funny' current (I_f) to pacemaker activity in the mammalian sinoatrial node (SAN), and the human SAN in particular, is still controversial and not fully established.

OBJECTIVE

To study the contribution of I_f to diastolic depolarization of human SAN cells and its dependence on heart rate, cAMP levels, and atrial load.

METHODS

HCN4 channels were expressed in human cardiac myocyte progenitor cells (CMPCs) and HCN4 currents assessed using perforated patch-clamp in traditional voltage clamp mode and during action potential clamp with human SAN-like action potential waveforms with 500-1500 ms cycle length, in absence or presence of forskolin to mimic β-adrenergic stimulation and a -15 mV command potential offset to mimic atrial load.

RESULTS

Forskolin significantly increased the fully-activated HCN4 current density at -140 mV by 14% and shifted the steady-state activation curve by +7.4 mV without affecting its slope. In addition, forskolin significantly accelerated current activation but slowed deactivation. The HCN4 current did not completely deactivate before the subsequent diastolic depolarization during action potential clamp. The amplitude of HCN4 current increased with increasing cycle length, was significantly larger in the presence of forskolin at all cycle lengths, and was significantly increased upon the negative offset to the command potential.

CONCLUSIONS

I_f is active during a human SAN action potential waveform and its amplitude is modulated by heart rate, β-adrenergic stimulation, and diastolic voltage range, such that I_f is under delicate control.

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.

Pacemaker Created in Human Ventricle by Depressing Inward-Rectifier K⁺ Current: A Simulation Study

Cardiac conduction disorders are common diseases which cause slow heart rate and syncope. The best way to treat these diseases by now is to implant electronic pacemakers, which, yet, have many disadvantages, such as the limited battery life and infection. Biopacemaker has been expected to replace the electronic devices. Automatic ventricular myocytes (VMs) could show pacemaker activity, which was induced by depressing inward-rectifier K⁺ current (I_K1). In this study, a 2D model of human biopacemaker was created from the ventricular endocardial myocytes. We examined the stability of the created biopacemaker and investigated its driving capability by finding the suitable size and spatial distribution of the pacemaker for robust pacing and driving the surrounding quiescent cardiomyocytes. Our results suggest that the rhythm of the pacemaker is similar to that of the single cell at final stable state. The driving force of the biopacemaker is closely related to the pattern of spatial distribution of the pacemaker.

Ca(2+) cycling properties are conserved despite bradycardic effects of heart failure in sinoatrial node cells

BACKGROUND

In animal models of heart failure (HF), heart rate decreases due to an increase in intrinsic cycle length (CL) of the sinoatrial node (SAN). Pacemaker activity of SAN cells is complex and modulated by the membrane clock, i.e., the ensemble of voltage gated ion channels and electrogenic pumps and exchangers, and the Ca(2+) clock, i.e., the ensemble of intracellular Ca(2+) ([Ca(2+)]i) dependent processes. HF in SAN cells results in remodeling of the membrane clock, but few studies have examined its effects on [Ca(2+)]i homeostasis.

METHODS

SAN cells were isolated from control rabbits and rabbits with volume and pressure overload-induced HF. [Ca(2+)]i concentrations, and action potentials (APs) and Na(+)-Ca(2+) exchange current (INCX) were measured using indo-1 and patch-clamp methodology, respectively.

RESULTS

The frequency of spontaneous [Ca(2+)]i transients was significantly lower in HF SAN cells (3.0 ± 0.1 (n = 40) vs. 3.4 ± 0.1 Hz (n = 45); mean ± SEM), indicating that intrinsic CL was prolonged. HF slowed the [Ca(2+)]i transient decay, which could be explained by the slower frequency and reduced sarcoplasmic reticulum (SR) dependent rate of Ca(2+) uptake. Other [Ca(2+)]i transient parameters, SR Ca(2+) content, INCX density, and INCX-[Ca(2+)]i relationship were all unaffected by HF. Combined AP and [Ca(2+)]i recordings demonstrated that the slower [Ca(2+)]i transient decay in HF SAN cells may result in increased INCX during the diastolic depolarization, but that this effect is likely counteracted by the HF-induced increase in intracellular Na(+). β-adrenergic and muscarinic stimulation were not changed in HF SAN cells, except that late diastolic [Ca(2+)]i rise, a prominent feature of the Ca(2+) clock, is lower during β-adrenergic stimulation.

CONCLUSIONS

HF SAN cells have a slower [Ca(2+)]i transient decay with limited effects on pacemaker activity. Reduced late diastolic [Ca(2+)]i rise during β-adrenergic stimulation may contribute to an impaired increase in intrinsic frequency in HF SAN cells.

Pacemaker activity of the human sinoatrial node: an update on the effects of mutations in HCN4 on the hyperpolarization-activated current

Since 2003, several loss-of-function mutations in the HCN4 gene, which encodes the HCN4 protein, have been associated with sinus node dysfunction. In human sinoatrial node (SAN), HCN4 is the most abundant of the four isoforms of the HCN family. Tetramers of HCN subunits constitute the ion channels that conduct the hyperpolarization-activated "funny" current (If), which plays an important modulating role in SAN pacemaker activity. Voltage-clamp experiments on HCN4 channels expressed in COS-7, CHO and HEK-293 cells, as well as in Xenopus oocytes have revealed changes in the expression and kinetics of mutant channels, but the extent to which especially the kinetic changes would affect If flowing during a human SAN action potential often remains unresolved. In our contribution to the Topical Collection on Human Single Nucleotide Polymorphisms and Disease Diagnostics, we provide an updated review of the mutation-induced changes in the expression and kinetics of HCN4 channels and provide an overview of their effects on If during the time course of a human SAN action potential, as assessed in simulated action potential clamp experiments. Future research may solve apparent inconsistencies between data from clinical studies and data from in vitro and in silico experiments.

Pacemaker activity of the human sinoatrial node: effects of HCN4 mutations on the hyperpolarization-activated current

The hyperpolarization-activated 'funny' current, If, plays an important modulating role in the pacemaker activity of the human sinoatrial node (SAN). If is carried by hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, which are tetramers built of four HCN subunits. In human SAN, HCN4 is the most abundant of the four isoforms of the HCN family. Since 2003, several loss-of-function mutations in the HCN4 gene, which encodes the HCN4 protein, or in the KCNE2 gene, which encodes the MiRP1 accessory β-subunit, have been associated with sinus node dysfunction. Voltage-clamp experiments on HCN4 channels expressed in COS-7 cells, Xenopus oocytes, or HEK-293 cells have revealed changes in the expression and kinetics of mutant channels, but the extent to which these changes would affect If flowing during a human SAN action potential is unresolved. Here, we review the changes in expression and kinetics of HCN4 mutant channels and provide an overview of their effects on If during the time course of a human SAN action potential, both under resting conditions and upon adrenergic stimulation. These effects are assessed in simulated action potential clamp experiments, with action potentials recorded from isolated human SAN pacemaker cells as command potential and kinetics of If based on voltage-clamp data from these cells. Results from in vitro and in silico experiments show several inconsistencies with clinical observations, pointing to challenges for future research.

Hyperpolarization-activated current, If, in mathematical models of rabbit sinoatrial node pacemaker cells

A typical feature of sinoatrial (SA) node pacemaker cells is the presence of an ionic current that activates upon hyperpolarization. The role of this hyperpolarization-activated current, I_f, which is also known as the “funny current” or “pacemaker current,” in the spontaneous pacemaker activity of SA nodal cells remains a matter of intense debate. Whereas some conclude that I_f plays a fundamental role in the generation of pacemaker activity and its rate control, others conclude that the role of I_f is limited to a modest contribution to rate control. The ongoing debate is often accompanied with arguments from computer simulations, either to support one's personal view or to invalidate that of the antagonist. In the present paper, we review the various mathematical descriptions of I_f that have been used in computer simulations and compare their strikingly different characteristics with our experimental data. We identify caveats and propose a novel model for I_f based on our experimental data.

Calcium transient and sodium-calcium exchange current in human versus rabbit sinoatrial node pacemaker cells

There is an ongoing debate on the mechanism underlying the pacemaker activity of sinoatrial node (SAN) cells, focusing on the relative importance of the “membrane clock” and the “Ca²⁺ clock” in the generation of the small net membrane current that depolarizes the cell towards the action potential threshold. Specifically, the debate centers around the question whether the membrane clock-driven hyperpolarization-activated current, I_f, which is also known as the “funny current” or “pacemaker current,” or the Ca²⁺ clock-driven sodium-calcium exchange current, I_NaCa, is the main contributor to diastolic depolarization. In our contribution to this journal's “Special Issue on Cardiac Electrophysiology,” we present a numerical reconstruction of I_f and I_NaCa in isolated rabbit and human SAN pacemaker cells based on experimental data on action potentials, I_f, and intracellular calcium concentration ([Ca²⁺]_i) that we have acquired from these cells. The human SAN pacemaker cells have a smaller I_f, a weaker [Ca²⁺]_i transient, and a smaller I_NaCa than the rabbit cells. However, when compared to the diastolic net membrane current, I_NaCa is of similar size in human and rabbit SAN pacemaker cells, whereas I_f is smaller in human than in rabbit cells.

An updated computational model of rabbit sinoatrial action potential to investigate the mechanisms of heart rate modulation

The cellular basis of cardiac pacemaking is still debated. Reliable computational models of the sinoatrial node (SAN) action potential (AP) may help gain a deeper understanding of the phenomenon. Recently, novel models incorporating detailed Ca(2+)-handling dynamics have been proposed, but they fail to reproduce a number of experimental data, and more specifically effects of 'funny' (I(f)) current modifications. We therefore developed a SAN AP model, based on available experimental data, in an attempt to reproduce physiological and pharmacological heart rate modulation. Cell compartmentalization and intracellular Ca(2+)-handling mechanisms were formulated as in the Maltsev-Lakatta model, focusing on Ca(2+)-cycling processes. Membrane current equations were revised on the basis of published experimental data. Modifications of the formulation of currents/pumps/exchangers to simulate I(f) blockers, autonomic modulators and Ca(2+)-dependent mechanisms (ivabradine, caesium, acetylcholine, isoprenaline, BAPTA) were derived from experimental data. The model generates AP waveforms typical of rabbit SAN cells, whose parameters fall within the experimental ranges: 352 ms cycle length, 80 mV AP amplitude, -58 mV maximum diastolic potential (MDP), 108 ms APD(50), and 7.1 Vs(-1) maximum upstroke velocity. Rate modulation by I(f) -blocking drugs agrees with experimental findings: 20% and 22% caesium-induced (5mM) and ivabradine-induced (3 μM) rate reductions, respectively, due to changes in diastolic depolarization (DD) slope, with no changes in either MDP or take-off potential (TOP). The model consistently reproduces the effects of autonomic modulation: 20% rate decrease with 10 nM acetylcholine and 28%increase with 1 μM isoprenaline, again entirely due to increase in the DD slope,with no changes in either MDP or TOP. Model testing of BAPTA effects showed slowing of rate, -26%, without cessation of beating. Our up-to-date model describes satisfactorily experimental data concerning autonomic stimulation, funny-channel blockade and inhibition of the Ca(2+)-related system by BAPTA, making it a useful tool for further investigation. Simulation results suggest that a detailed description of the intracellular Ca(2+) fluxes is fully compatible with the observation that I(f) is a major component of pacemaking and rate modulation.

Rate control in atrial fibrillation by cooling: effect of temperature on dromotropy in perfused rabbit hearts

Background. Cooling has emerged as a therapeutic option in critically ill patients (especially after cardiac resuscitation) and might also have a negative dromotropic effect in atrial fibrillation. We sought to determine the impact of cooling on electrophysiologic properties of Langendorff-perfused rabbit hearts. Methods and Results. In 20 isolated Langendorff-perfused rabbit hearts, the temperature of the tissue bath was changed between 17 and 42°C. With decreasing temperature, significant increases of the spontaneous sinus cycle length, decreases of the mean ventricular heart rate during atrial fibrillation, and relevant increases of atrial and ventricular refractory periods were observed (ANOVA P < .01). Conclusions. Cardiac hypothermia leads to a significant drop of mean ventricular heart rate during atrial fibrillation. Negative chronotropy and dromotropy induced by moderate cardiac hypothermia might be a feasible therapeutic approach in patients with hemodynamically relevant tachyarrhythmias in a CCU/ICU setting.

I(f) inhibition in cardiovascular diseases

Heart rate (HR) is determined by the pacemaker activity of cells from the sinoatrial node (SAN), located in the right atria. Spontaneous electrical activity of SAN cells results from a diastolic depolarization (DD). Despite controversy in the exact contribution of funny current (I(f)) in pacemaking, it is a major contributor of DD. I(f) is an inward Na(+)/K(+) current, activated upon hyperpolarization and directly modulated by cyclic adenosine monophosphate. The f-proteins are hyperpolarization-activated cyclic nucleotide-gated channels, HCN4 being the main isoform of SAN. Ivabradine (IVA) decreases DD and inhibits I(f) in a use-dependent manner. Under normal conditions IVA selectively reduces HR and limits exercise-induced tachycardia, in animals and young volunteers. Reduction in HR with IVA both decreases myocardial oxygen consumption and increases its supply due to prolongation of diastolic perfusion time. In animal models and in human with coronary artery disease (CAD), IVA has anti-anginal and anti-ischemic efficacy, equipotent to classical treatments, β-blockers, or calcium channel blockers. As expected from its selectivity for I(f), the drug is safe and well tolerated with minor visual side effects. As a consequence, IVA is the first inhibitor of I(f) approved for the treatment of stable angina. Available clinical data indicate that IVA could improve the management of stable angina in all patients including those treated with β-blockers. As chronic elevation of resting HR is an independent predictor of mortality, pure HR reduction by inhibition of I(f) could, beyond the control of anti-anginal symptoms, improve the prognosis of CAD and heart failure; this therapeutic potential is currently under evaluation with IVA.