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

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.

Deciphering cellular signals in adult mouse sinoatrial node cells

Sinoatrial node (SAN) cells are the pacemakers of the heart. This study describes a method for culturing and infection of adult mouse SAN cells with FRET-based biosensors that can be exploited to examine signaling events. SAN cells cultured in media with blebbistatin or (S)-nitro-blebbistatin retain their morphology, protein distribution, action potential (AP) waveform, and cAMP dynamics for at least 40 h. SAN cells expressing targeted cAMP sensors show distinct β-adrenergic-mediated cAMP pools. Cyclic GMP, protein kinase A, Ca2+/CaM kinase II, and protein kinase D in SAN cells also show unique dynamics to different stimuli. Heart failure SAN cells show a decrease in cAMP and cGMP levels. In summary, a reliable method for maintaining adult mouse SAN cells in culture is presented, which facilitates studies of signaling networks and regulatory mechanisms during physiological and pathological conditions.

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

The effect of forskolin on membrane clock and calcium clock in the hypoxic/reoxygenation of sinoatrial node cells and its mechanism

BACKGROUND

In this study, we investigated the effect of forskolin (FSK, a selective adenylate cyclase agonist) on the automatic diastolic depolarization of sinus node cells (SNC) with hypoxia/reoxygenation (H/R) injury.

METHODS

The SNC of the newborn rat was randomly assigned into the control group, the H/R (H/R injury) group, or the H/R + FSK (H/R injury + FSK treatment) group. Patch-clamp was performed to record the action potential and electrophysiological changes. The cellular distribution of intracellular calcium concentration was analyzed by fluorescence staining.

RESULTS

Compared with the control cells, spontaneous pulsation frequency (SPF) and diastolic depolarization rate (DDR) of H/R cells were reduced from 244.3 ± 10.6 times/min and 108.7 ± 7.8 mV/s to 130.5 ± 7.6 times/min and 53.4 ± 6.5 mV/s, respectively. FSK significantly increased SPF and DDR of H/R cells to 208.3 ± 8.3 times/min and 93.2 ± 8.9 mV/s (n = 15, both p < 0.01), respectively. H/R reduced the current densities of I_f, I_Ca,T and inward I_NCX, which were significantly increased by 10 μM FSK treatment (n = 15, p < 0.01). Furthermore, reduced expression of HCN4 and NCX1.1 channel protein were significantly increased by FSK. Inhibitor studies showed that both SQ22536 (a selective adenylate cyclase inhibitor) and H89 (a selective protein kinases A [PKA] inhibitor) blocked the effects of FSK on SPF and DDR.

CONCLUSIONS

H/R causes pacemaker dysfunction in newborn rat sinoatrial node cells leading to divergence of the DD and the slow of spontaneous APs, which change can be dramatically reversed by FSK through increasing I_NCX and I_f current in H/R injury.

Distinct mechanisms mediate pacemaker dysfunction associated with catecholaminergic polymorphic ventricular tachycardia mutations: Insights from computational modeling

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is a stress-induced ventricular arrhythmia associated with rhythm disturbance and impaired sinoatrial node cell (SANC) automaticity (pauses). Mutations associated with dysfunction of Ca²⁺-related mechanisms have been shown to be present in CPVT. These dysfunctions include impaired Ca²⁺ release from the ryanodine receptor (i.e., RyR2^R4496C mutation) or binding to calsequestrin 2 (CASQ2). In SANC, Ca²⁺ signaling directly and indirectly mediates pacemaker function. We address here the following research questions: (i) what coupled-clock mechanisms and pathways mediate pacemaker mutations associated with CPVT in basal and in response to β-adrenergic stimulation? (ii) Can different mechanisms lead to the same CPVT-related pacemaker pauses? (iii) Can the mutation-induced deteriorations in SANC function be reversed by drug intervention or gene manipulation? We used a numerical model of mice SANC that includes membrane and intracellular mechanisms and their interconnected signaling pathways. In the basal state of RyR2^R4496C SANC, the model predicted that the Na⁺-Ca²⁺ exchanger current (I_NCX) and T-type Ca²⁺ current (I_CaT) mediate between changes in Ca²⁺ signaling and SANC dysfunction. Under β-adrenergic stimulation, changes in cAMP-PKA signaling and the sodium currents (I_Na), in addition to I_NCX and I_CaT, mediate between changes in Ca²⁺ signaling and SANC automaticity pauses. Under basal conditions in Casq2^-/-, the same mechanisms drove changes in Ca²⁺ signaling and subsequent pacemaker dysfunction. However, SANC automaticity pauses in response to β-AR stimulation were mediated by I_CaT and I_Na. Taken together, distinct mechanisms can lead to CPVT-associated SANC automaticity pauses. In addition, we predict that specifically increasing SANC cAMP-PKA activity by either a pharmacological agent (IBMX, a phosphodiesterase (PDE) inhibitor), gene manipulation (overexpression of adenylyl cyclase 1/8) or direct manipulation of the SERCA phosphorylation target through changes in gene expression, compensate for the impairment in SANC automaticity. These findings suggest new insights for understanding CPVT and its therapeutic approach.

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.

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.

Understanding the Phosphorylation Mechanism by Using Quantum Chemical Calculations and Molecular Dynamics Simulations

Phosphorylation is one of the most frequent post-translational modifications on proteins. It regulates many cellular processes by modulation of phosphorylation on protein structure and dynamics. However, the mechanism of phosphorylation-induced conformational changes of proteins is still poorly understood. Here, we report a computational study of three representative groups of tyrosine in ADP-ribosylhydrolase 1, serine in BTG2, and serine in Sp100C by using six molecular dynamics (MD) simulations and quantum chemical calculations. Added phosphorylation was found to disrupt hydrogen bond, and increase new weak interactions (hydrogen bond and hydrophobic interaction) during MD simulations, leading to conformational changes. Quantum chemical calculations further indicate that the phosphorylation on tyrosine, threonine, and serine could decrease the optical band gap energy (Egap), which can trigger electronic transitions to form or disrupt interactions easily. Our results provide an atomic and electronic description of how phosphorylation facilitates conformational and dynamic changes in proteins, which may be useful for studying protein function and protein design.

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.