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

Circadian regulation of sinoatrial nodal cell pacemaking function: Dissecting the roles of autonomic control, body temperature, and local circadian rhythmicity

Strong circadian (~24h) rhythms in heart rate (HR) are critical for flexible regulation of cardiac pacemaking function throughout the day. While this circadian flexibility in HR is sustained in diverse conditions, it declines with age, accompanied by reduced maximal HR performance. The intricate regulation of circadian HR involves the orchestration of the autonomic nervous system (ANS), circadian rhythms of body temperature (CRBT), and local circadian rhythmicity (LCR), which has not been fully understood. Here, we developed a mathematical model describing ANS, CRBT, and LCR in sinoatrial nodal cells (SANC) that accurately captures distinct circadian patterns in adult and aged mice. Our model underscores how the alliance among ANS, CRBT, and LCR achieves circadian flexibility to cover a wide range of firing rates in SANC, performance to achieve maximal firing rates, while preserving robustness to generate rhythmic firing patterns irrespective of external conditions. Specifically, while ANS dominates in promoting SANC flexibility and performance, CRBT and LCR act as primary and secondary boosters, respectively, to further enhance SANC flexibility and performance. Disruption of this alliance with age results in impaired SANC flexibility and performance, but not robustness. This unexpected outcome is primarily attributed to the age-related reduction in parasympathetic activities, which maintains SANC robustness while compromising flexibility. Our work sheds light on the critical alliance of ANS, CRBT, and LCR in regulating time-of-day cardiac pacemaking function and dysfunction, offering insights into novel therapeutic targets for the prevention and treatment of cardiac arrhythmias.

Age-dependent contribution of intrinsic mechanisms to sinoatrial node function in humans

Average beat interval (BI) and beat interval variability (BIV) are primarily determined by mutual entrainment between the autonomic-nervous system (ANS) and intrinsic mechanisms that govern sinoatrial node (SAN) cell function. While basal heart rate is not affected by age in humans, age-dependent reductions in intrinsic heart rate have been documented even in so-called healthy individuals. The relative contributions of the ANS and intrinsic mechanisms to age-dependent deterioration of SAN function in humans are not clear. We recorded ECG on patients (n = 16 < 21 years and n = 23 41-78 years) in the basal state and after ANS blockade (propranolol and atropine) in the presence of propofol and dexmedetomidine anesthesia. Average BI and BIV were analyzed. A set of BIV features were tested to designated the "signatures" of the ANS and intrinsic mechanisms and also the anesthesia "signature". In young patients, the intrinsic mechanisms and ANS mainly contributed to long- and short-term BIV, respectively. In adults, both ANS and intrinsic mechanisms contributed to short-term BIV, while the latter also contributed to long-term BIV. Furthermore, anesthesia affected ANS function in young patients and both mechanisms in adult. The work also showed that intrinsic mechanism features can be calculated from BIs, without intervention.

Human Sinoatrial Node Pacemaker Activity: Role of the Slow Component of the Delayed Rectifier K+ Current, IKs

The pacemaker activity of the sinoatrial node (SAN) has been studied extensively in animal species but is virtually unexplored in humans. Here we assess the role of the slowly activating component of the delayed rectifier K+ current (IKs) in human SAN pacemaker activity and its dependence on heart rate and β-adrenergic stimulation. HEK-293 cells were transiently transfected with wild-type KCNQ1 and KCNE1 cDNA, encoding the α- and β-subunits of the IKs channel, respectively. KCNQ1/KCNE1 currents were recorded both during a traditional voltage clamp and during an action potential (AP) clamp with human SAN-like APs. Forskolin (10 µmol/L) was used to increase the intracellular cAMP level, thus mimicking β-adrenergic stimulation. The experimentally observed effects were evaluated in the Fabbri-Severi computer model of an isolated human SAN cell. Transfected HEK-293 cells displayed large IKs-like outward currents in response to depolarizing voltage clamp steps. Forskolin significantly increased the current density and significantly shifted the half-maximal activation voltage towards more negative potentials. Furthermore, forskolin significantly accelerated activation without affecting the rate of deactivation. During an AP clamp, the KCNQ1/KCNE1 current was substantial during the AP phase, but relatively small during diastolic depolarization. In the presence of forskolin, the KCNQ1/KCNE1 current during both the AP phase and diastolic depolarization increased, resulting in a clearly active KCNQ1/KCNE1 current during diastolic depolarization, particularly at shorter cycle lengths. Computer simulations demonstrated that IKs reduces the intrinsic beating rate through its slowing effect on diastolic depolarization at all levels of autonomic tone and that gain-of-function mutations in KCNQ1 may exert a marked bradycardic effect during vagal tone. In conclusion, IKs is active during human SAN pacemaker activity and has a strong dependence on heart rate and cAMP level, with a prominent role at all levels of autonomic tone.

Measuring and modeling the effects of vagus nerve stimulation on heart rate and laryngeal muscles

BACKGROUND

Reduced heart rate (HR) during vagus nerve stimulation (VNS) is associated with therapy for heart failure, but stimulation frequency and amplitude are limited by patient tolerance. An understanding of physiological responses to parameter adjustments would allow differential control of therapeutic and side effects. To investigate selective modulation of the physiological responses to VNS, we quantified the effects and interactions of parameter selection on two physiological outcomes: one related to therapy (reduced HR) and one related to side effects (laryngeal muscle EMG).

METHODS

We applied a broad range of stimulation parameters (mean pulse rates (MPR), intra-burst frequencies, and amplitudes) to the vagus nerve of anesthetized mice. We leveraged the in vivo recordings to parameterize and validate computational models of HR and laryngeal muscle activity across amplitudes and temporal patterns of VNS. We constructed a finite element model of excitation of fibers within the mouse cervical vagus nerve.

RESULTS

HR decreased with increased amplitude, increased MPR, and decreased intra-burst frequency. EMG increased with increased MPR. Preferential HR effects over laryngeal EMG effects required combined adjustments of amplitude and MPR. The model of HR responses highlighted contributions of ganglionic filtering to VNS-evoked changes in HR at high stimulation frequencies. Overlap in activation thresholds between small and large modeled fibers was consistent with the overlap in dynamic ranges of related physiological measures (HR and EMG).

CONCLUSION

The present study provides insights into physiological responses to VNS required for informed parameter adjustment to modulate selectively therapeutic effects and side effects.

Progress on role of ion channels of cardiac fibroblasts in fibrosis

Cardiac fibrosis is defined as excessive deposition of extracellular matrix (ECM) in pathological conditions. Cardiac fibroblasts (CFs) activated by injury or inflammation differentiate into myofibroblasts (MFs) with secretory and contractile functions. In the fibrotic heart, MFs produce ECM which is composed mainly of collagen and is initially involved in maintaining tissue integrity. However, persistent fibrosis disrupts the coordination of excitatory contractile coupling, leading to systolic and diastolic dysfunction, and ultimately heart failure. Numerous studies have demonstrated that both voltage- and non-voltage-gated ion channels alter intracellular ion levels and cellular activity, contributing to myofibroblast proliferation, contraction, and secretory function. However, an effective treatment strategy for myocardial fibrosis has not been established. Therefore, this review describes the progress made in research related to transient receptor potential (TRP) channels, Piezo1, Ca²⁺ release-activated Ca²⁺ (CRAC) channels, voltage-gated Ca²⁺ channels (VGCCs), sodium channels, and potassium channels in myocardial fibroblasts with the aim of providing new ideas for treating myocardial fibrosis.

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.

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

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

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

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

Inherited and Acquired Rhythm Disturbances in Sick Sinus Syndrome, Brugada Syndrome, and Atrial Fibrillation: Lessons from Preclinical Modeling

Rhythm disturbances are life-threatening cardiovascular diseases, accounting for many deaths annually worldwide. Abnormal electrical activity might arise in a structurally normal heart in response to specific triggers or as a consequence of cardiac tissue alterations, in both cases with catastrophic consequences on heart global functioning. Preclinical modeling by recapitulating human pathophysiology of rhythm disturbances is fundamental to increase the comprehension of these diseases and propose effective strategies for their prevention, diagnosis, and clinical management. In silico, in vivo, and in vitro models found variable application to dissect many congenital and acquired rhythm disturbances. In the copious list of rhythm disturbances, diseases of the conduction system, as sick sinus syndrome, Brugada syndrome, and atrial fibrillation, have found extensive preclinical modeling. In addition, the electrical remodeling as a result of other cardiovascular diseases has also been investigated in models of hypertrophic cardiomyopathy, cardiac fibrosis, as well as arrhythmias induced by other non-cardiac pathologies, stress, and drug cardiotoxicity. This review aims to offer a critical overview on the effective ability of in silico bioinformatic tools, in vivo animal studies, in vitro models to provide insights on human heart rhythm pathophysiology in case of sick sinus syndrome, Brugada syndrome, and atrial fibrillation and advance their safe and successful translation into the cardiology arena.

Phosphoprotein Phosphatase 1 but Not 2A Activity Modulates Coupled-Clock Mechanisms to Impact on Intrinsic Automaticity of Sinoatrial Nodal Pacemaker Cells

Spontaneous AP (action potential) firing of sinoatrial nodal cells (SANC) is critically dependent on protein kinase A (PKA) and Ca²⁺/calmodulin-dependent protein kinase II (CaMKII)-dependent protein phosphorylation, which are required for the generation of spontaneous, diastolic local Ca²⁺ releases (LCRs). Although phosphoprotein phosphatases (PP) regulate protein phosphorylation, the expression level of PPs and phosphatase inhibitors in SANC and the impact of phosphatase inhibition on the spontaneous LCRs and other players of the oscillatory coupled-clock system is unknown. Here, we show that rabbit SANC express both PP1, PP2A, and endogenous PP inhibitors I-1 (PPI-1), dopamine and cyclic adenosine 3′,5′-monophosphate (cAMP)-regulated phosphoprotein (DARPP-32), kinase C-enhanced PP1 inhibitor (KEPI). Application of Calyculin A, (CyA), a PPs inhibitor, to intact, freshly isolated single SANC: (1) significantly increased phospholamban (PLB) phosphorylation (by 2–3-fold) at both CaMKII-dependent Thr¹⁷ and PKA-dependent Ser¹⁶ sites, in a time and concentration dependent manner; (2) increased ryanodine receptor (RyR) phosphorylation at the Ser²⁸⁰⁹ site; (3) substantially increased sarcoplasmic reticulum (SR) Ca²⁺ load; (4) augmented L-type Ca²⁺ current amplitude; (5) augmented LCR’s characteristics and decreased LCR period in intact and permeabilized SANC, and (6) increased the spontaneous basal AP firing rate. In contrast, the selective PP2A inhibitor okadaic acid (100 nmol/L) had no significant effect on spontaneous AP firing, LCR parameters, or PLB phosphorylation. Application of purified PP1 to permeabilized SANC suppressed LCR, whereas purified PP2A had no effect on LCR characteristics. Our numerical model simulations demonstrated that PP inhibition increases AP firing rate via a coupled-clock mechanism, including respective increases in the SR Ca²⁺ pumping rate, L-type Ca²⁺ current, and Na⁺/Ca²⁺-exchanger current. Thus, PP1 and its endogenous inhibitors modulate the basal spontaneous firing rate of cardiac pacemaker cells by suppressing SR Ca²⁺ cycling protein phosphorylation, the SR Ca²⁺ load and LCRs, and L-type Ca²⁺ current.