Cardiac beta-adrenergic signaling: from subcellular microdomains to heart failure

Molecular and cellular neurocardiology in heart disease

This paper updates and builds on a previous White Paper in this journal that some of us contributed to concerning the molecular and cellular basis of cardiac neurobiology of heart disease. Here we focus on recent findings that underpin cardiac autonomic development, novel intracellular pathways and neuroplasticity. Throughout we highlight unanswered questions and areas of controversy. Whilst some neurochemical pathways are already demonstrating prognostic viability in patients with heart failure, we also discuss the opportunity to better understand sympathetic impairment by using patient specific stem cells that provides pathophysiological contextualization to study 'disease in a dish'. Novel imaging techniques and spatial transcriptomics are also facilitating a road map for target discovery of molecular pathways that may form a therapeutic opportunity to treat cardiac dysautonomia.

Action potential responses to changes in stimulation frequency and isoproterenol in rat ventricular myocytes

PURPOSE

Current understanding of ventricular action potential adaptation to physiological stress is generally based on protocols using non-physiological rates and conditions isolating rate effects from escalating adrenergic stimulation. To permit refined understanding, ventricular action potentials were assessed across physiological pacing frequencies in the presence and absence of adrenergic stimuli. Isolated and combined effects were analyzed to assess their ability to replicate in-vivo responses.

METHODS

Steady-state action potentials from ventricular myocytes isolated from male Wistar rats (3 months; N = 8 animals) were recorded at 37°C with steady-state pacing at 1, 2, 4, 6, 8 and 10 Hz using whole-cell patch-clamp. Action potential repolarization to 25, 50, 75, 90 and 100% of full repolarization (APD25-100 ) was compared before and after 5 nM, 100 nM and 1 µM isoproterenol doses.

RESULTS

A Repeated measures ANOVA found APD50-90  shortened with 5 nM isoproterenol infusion by 6-25% (but comparable across doses) (p ≤ 0.03). Pacing frequencies emulating a normal rat heart rate (6 Hz) prolonged APD50 23% compared with 1 Hz pacing. Frequencies emulating exercise or stress (10 Hz) shortened APD90 (29%).

CONCLUSION

These results demonstrate modest action potential shortening in response to adrenergic stimulation and elevations in pacing beyond physiological resting rates. Our findings indicate changes in action potential plateau and late repolarization predominantly underlie simulated exercise responses in the rat heart. This work provides novel action potential reference data and will help model cardiac responses to physiological stimuli in the rat heart via computational techniques.

Electrophysiological Remodeling: Cardiac T-Tubules and ß-Adrenoceptors

Beta-adrenoceptors (βAR) are often viewed as archetypal G-protein coupled receptors. Over the past fifteen years, investigations in cardiovascular biology have provided remarkable insights into this receptor family. These studies have shifted pharmacological dogma, from one which centralized the receptor to a new focus on structural micro-domains such as caveolae and t-tubules. Important studies have examined, separately, the structural compartmentation of ion channels and βAR. Despite links being assumed, relatively few studies have specifically examined the direct link between structural remodeling and electrical remodeling with a focus on βAR. In this review, we will examine the nature of receptor and ion channel dysfunction on a substrate of cardiomyocyte microdomain remodeling, as well as the likely ramifications for cardiac electrophysiology. We will then discuss the advances in methodologies in this area with a specific focus on super-resolution microscopy, fluorescent imaging, and new approaches involving microdomain specific, polymer-based agonists. The advent of powerful computational modelling approaches has allowed the science to shift from purely empirical work, and may allow future investigations based on prediction. Issues such as the cross-reactivity of receptors and cellular heterogeneity will also be discussed. Finally, we will speculate as to the potential developments within this field over the next ten years.

Reduced cardiac response to the adrenergic system is a key limiting factor for physical capacity in old age

Ageing is associated with a progressive reduction in physical capacity reducing quality of life. One key physiological limitation of physical capacity that deteriorates in a progressive age-dependent manner is cardiac reserve. Peak cardiac output falls progressively with advancing age such that in extreme old age there is limited ability to enhance cardiac output beyond basal function as is required to support the increased metabolic needs of physical activity. This loss of dynamic range in cardiac output associates with a progressive reduction in the heart's response to adrenergic stimulation. A combination of decreases in the expression and functioning of beta1 adrenergic receptors partially underlies this change. Changes in end effector proteins also have a role to play in this decline. Alterations in the efficiency of excitation-contraction coupling contribute to the reduced chronotropic, inotropic and lusitropic responses of the aged heart. Moderate to vigorous endurance exercise training however has some potential to counter elements of these changes. Further studies are required to fully elucidate the key pivotal mechanisms involved in the age-related loss of response to adrenergic signalling to allow targeted therapeutic strategies to be developed with the aim of preserving physical capacity in advanced old age.

Cardiac sensory afferents modulate susceptibility to anxio-depressive behaviour in a mouse model of chronic heart failure

AIM

Impairments in cerebral structure and cognitive performance in chronic heart failure (CHF) are critical components of its comorbidity spectrum. Autonomic afferents that arise from cardiac sensory fibres show enhanced activity with CHF. Desensitization of these fibres by local application of resiniferatoxin (RTX) during myocardial infarction (MI) is known to prevent cardiac hypertrophy, sympathetic hyperactivity and CHF. Whether these afferents mediate cerebral allostasis is unknown.

METHODS

CHF was induced by myocardial infarction. To evaluate if cardiac afferents contribute to cerebral allostasis, RTX was acutely applied to the pericardial space in controls (RTX) and in MI treated animals (MI/RTX). Subjects were then evaluated in a series of behavioural tests recapitulating different symptoms of depressive disorders. Proteomics of the frontal cortices (FC) was performed to identify contributing proteins and pathways responsible for behavioural allostasis.

RESULTS

Desensitization of cardiac afferents relieves hallmarks of an anxio/depressive-like state in mice. Unique protein signatures and regulatory pathways in FCs isolated from each treatment reveal the degree of complexity inherent in the FC response to stresses originating in the heart. While cortices from the combined treatment (MI/RTX) did not retain protein signatures from the individual treatment groups, all three groups suffer dysregulation in circadian entrainment.

CONCLUSION

CHF is comorbid with an anxio/depressive-like state and ablation of cardiac afferents relieves the despair phenotype. The strikingly different proteomic profiles observed in FCs suggest that MI and RTX lead to unique brain-signalling patterns and that the combined treatment, potentially through destructive interference mechanisms, most closely resembles controls.

Computational modeling approaches to cAMP/PKA signaling in cardiomyocytes

The cAMP/PKA pathway is a fundamental regulator of excitation-contraction coupling in cardiomyocytes. Activation of cAMP has a variety of downstream effects on cardiac function including enhanced contraction, accelerated relaxation, adaptive stress response, mitochondrial regulation, and gene transcription. Experimental advances have shed light on the compartmentation of cAMP and PKA, which allow for control over the varied targets of these second messengers and is disrupted in heart failure conditions. Computational modeling is an important tool for understanding the spatial and temporal complexities of this system. In this review article, we outline the advances in computational modeling that have allowed for deeper understanding of cAMP/PKA dynamics in the cardiomyocyte in health and disease, and explore new modeling frameworks that may bring us closer to a more complete understanding of this system. We outline various compartmental and spatial signaling models that have been used to understand how β-adrenergic signaling pathways function in a variety of simulation conditions. We also discuss newer subcellular models of cardiovascular function that may be used as templates for the next phase of computational study of cAMP and PKA in the heart, and outline open challenges which are important to consider in future models.

Genetic background influences expression and function of the cation channel TRPM4 in the mouse heart

Transient receptor potential melastatin 4 (TRPM4) cation channels act in cardiomyocytes as a negative modulator of the L-type Ca²⁺ current. Ubiquitous Trpm4 deletion in mice leads to an increased β-adrenergic inotropy in healthy mice as well as after myocardial infarction. In this study, we set out to investigate cardiac inotropy in mice with cardiomyocyte-specific Trpm4 deletion. The results guided us to investigate the relevance of TRPM4 for catecholamine-evoked Ca²⁺ signaling in cardiomyocytes and inotropy in vivo in TRPM4-deficient mouse models of different genetic background. Cardiac hemodynamics were investigated using pressure-volume analysis. Surprisingly, an increased β-adrenergic inotropy was observed in global TRPM4-deficient mice on a 129SvJ genetic background, but the inotropic response was unaltered in mice with global and cardiomyocyte-specific TRPM4 deletion on the C57Bl/6N background. We found that the expression of TRPM4 proteins is about 78 ± 10% higher in wild-type mice on the 129SvJ versus C57Bl/6N background. In accordance with contractility measurements, our analysis of the intracellular Ca²⁺ transients revealed an increase in ISO-evoked Ca²⁺ rise in Trpm4-deficient cardiomyocytes of the 129SvJ strain, but not of the C57Bl/6N strain. No significant differences were observed between the two mouse strains in the expression of other regulators of cardiomyocyte Ca²⁺ homeostasis. We conclude that the relevance of TRPM4 for cardiac contractility depends on homeostatic TRPM4 expression levels or the genetic endowment in different mouse strains as well as on the health/disease status. Therefore, the concept of inhibiting TRPM4 channels to improve cardiac contractility needs to be carefully explored in specific strains and species and prospectively in different genetically diverse populations of patients.

Therapeutic Implications for PDE2 and cGMP/cAMP Mediated Crosstalk in Cardiovascular Diseases

Phosphodiesterases (PDEs) are the principal superfamily of enzymes responsible for degrading the secondary messengers 3',5'-cyclic nucleotides cAMP and cGMP. Their refined subcellular localization and substrate specificity contribute to finely regulate cAMP/cGMP gradients in various cellular microdomains. Redistribution of multiple signal compartmentalization components is often perceived under pathological conditions. Thereby PDEs have long been pursued as therapeutic targets in diverse disease conditions including neurological, metabolic, cancer and autoimmune disorders in addition to numerous cardiovascular diseases (CVDs). PDE2 is a unique member of the broad family of PDEs. In addition to its capability to hydrolyze both cAMP and cGMP, PDE2 is the sole isoform that may be allosterically activated by cGMP increasing its cAMP hydrolyzing activity. Within the cardiovascular system, PDE2 serves as an integral regulator for the crosstalk between cAMP/cGMP pathways and thereby may couple chronically adverse augmented cAMP signaling with cardioprotective cGMP signaling. This review provides a comprehensive overview of PDE2 regulatory functions in multiple cellular components within the cardiovascular system and also within various subcellular microdomains. Implications for PDE2- mediated crosstalk mechanisms in diverse cardiovascular pathologies are discussed highlighting the prospective use of PDE2 as a potential therapeutic target in cardiovascular disorders.

Slow Delayed Rectifier Current Protects Ventricular Myocytes From Arrhythmic Dynamics Across Multiple Species: A Computational Study

BACKGROUND

The slow and rapid delayed rectifier K⁺ currents (I_Ks and I_Kr, respectively) are responsible for repolarizing the ventricular action potential (AP) and preventing abnormally long APs that may lead to arrhythmias. Although differences in biophysical properties of the 2 currents have been carefully documented, the respective physiological roles of I_Kr and I_Ks are less established. In this study, we sought to understand the individual roles of these currents and quantify how effectively each stabilizes the AP and protects cells against arrhythmias across multiple species.

METHODS

We compared 10 mathematical models describing ventricular myocytes from human, rabbit, dog, and guinea pig. We examined variability within heterogeneous cell populations, tested the susceptibility of cells to proarrhythmic behavior, and studied how I_Ks and I_Kr responded to changes in the AP.

RESULTS

We found that (1) models with higher baseline I_Ks exhibited less cell-to-cell variability in AP duration; (2) models with higher baseline I_Ks were less susceptible to early afterdepolarizations induced by depolarizing perturbations; (3) as AP duration is lengthened, I_Ks increases more profoundly than I_Kr, thereby providing negative feedback that resists excessive AP prolongation; and (4) the increase in I_Ks that occurs during β-adrenergic stimulation is critical for protecting cardiac myocytes from early afterdepolarizations under these conditions.

CONCLUSIONS

Slow delayed rectifier current is uniformly protective across a variety of cell types. These results suggest that I_Ks enhancement could potentially be an effective antiarrhythmic strategy.

A Microwell Cell Capture Device Reveals Variable Response to Dobutamine in Isolated Cardiomyocytes

Isolated ventricular cardiomyocytes exhibit substantial cell-to-cell variability, even when obtained from the same small volume of myocardium. In this study, we investigated the possibility that cardiomyocyte responses to β-adrenergic stimulus are also highly heterogeneous. To achieve the throughput and measurement duration desired for these experiments, we designed and validated a novel microwell system that immobilizes and uniformly orients isolated adult cardiomyocytes. In this configuration, detailed drug responses of dozens of cells can be followed for intervals exceeding 1 h. At the conclusion of an experiment, specific cells can also be harvested via a precision aspirator for single-cell gene expression profiling. Using this system, we followed changes in Ca²⁺ signaling and contractility of individual cells under sustained application of either dobutamine or omecamtiv mecarbil. Both compounds increased average cardiomyocyte contractility over the course of an hour, but responses of individual cells to dobutamine were significantly more variable. Surprisingly, some dobutamine-treated cardiomyocytes augmented Ca²⁺ release without increasing contractility. Other cells responded with increased contractility despite unchanged Ca²⁺ release. Single-cell gene expression analysis revealed significant co-expression of β-adrenergic pathway genes PKA regulatory subunit type I, PKA regulatory subunit type II, and Ca²⁺/calmodulin-dependent protein kinase II across cardiomyocytes. Other data supported a connection between the effects of dobutamine on relaxation rate and the expression of protein phosphatase 2. These findings suggest that variable drug responses among cells are not merely experimental artifacts. By enabling direct comparison of the functional behavior of an individual cell and the genes it expresses, this new system constitutes a unique tool for interrogating cardiomyocyte drug responses and discovering the genes that modulate them.

Studying β1 and β2 adrenergic receptor signals in cardiac cells using FRET-based sensors

Cyclic 3'-5' adenosine monophosphate (cAMP) is a key modulator of cardiac function. Thanks to the sophisticated organization of its pathway in distinct functional units called microdomains, cAMP is involved in the regulation of both inotropy and chronotropy as well as transcription and cardiac death. While visualization of cAMP microdomains can be achieved thanks to cAMP-sensitive FRET-based sensors, the molecular mechanisms through which cAMP-generating stimuli are coupled to distinct functional outcomes are not well understood. One possibility is that each stimulus activates multiple microdomains in order to generate a spatiotemporal code that translates into function. To test this hypothesis here we propose a series of experimental protocols that allow to simultaneously follow cAMP or Protein Kinase A (PKA)-dependent phosphorylation in different subcellular compartments of living cells. We investigate the responses of β Adrenergic receptors (β₁AR and β₂AR) challenged with selective drugs that enabled us to measure the actions of each receptor independently. At the whole cell level, we used a combination of co-culture with selective βAR stimulation and were able to molecularly separate cardiac fibroblasts from neonatal rat ventricular myocytes based on their cAMP responses. On the other hand, at the subcellular level, these experimental protocols allowed us to dissect the relative weight of β₁ and β₂ adrenergic receptors on cAMP signalling at the cytosol and outer mitochondrial membrane of NRVMs. We propose that experimental procedures that allow the collection of multiparametric data are necessary in order to understand the molecular mechanisms underlying the coupling between extracellular signals and cellular responses.

Genetically Encoded Fluorescent Biosensors Illuminate the Spatiotemporal Regulation of Signaling Networks

Cellular signaling networks are the foundation which determines the fate and function of cells as they respond to various cues and stimuli. The discovery of fluorescent proteins over 25 years ago enabled the development of a diverse array of genetically encodable fluorescent biosensors that are capable of measuring the spatiotemporal dynamics of signal transduction pathways in live cells. In an effort to encapsulate the breadth over which fluorescent biosensors have expanded, we endeavored to assemble a comprehensive list of published engineered biosensors, and we discuss many of the molecular designs utilized in their development. Then, we review how the high temporal and spatial resolution afforded by fluorescent biosensors has aided our understanding of the spatiotemporal regulation of signaling networks at the cellular and subcellular level. Finally, we highlight some emerging areas of research in both biosensor design and applications that are on the forefront of biosensor development.

Arrhythmia initiation in catecholaminergic polymorphic ventricular tachycardia type 1 depends on both heart rate and sympathetic stimulation

Aims

Catecholaminergic polymorphic ventricular tachycardia type 1 (CPVT1) predisposes to ventricular tachyarrhythmias (VTs) during high heart rates due to physical or psychological stress. The essential role of catecholaminergic effects on ventricular cardiomyocytes in this situation is well documented, but the importance of heart rate per se for arrhythmia initiation in CPVT1 is largely unexplored.

Methods and results

Sixteen CPVT1 patients performed a bicycle stress-test. Occurrence of VT triggers, i.e. premature ventricular complexes (PVC), depended on high heart rate, with individual thresholds. Atrial pacing above the individual PVC threshold in three patients did not induce PVCs.

The underlying mechanism for the clinical observation was explored using cardiomyocytes from mice with the RyR2-R2474S (RyR2-RS) mutation, which exhibit exercise-induced VTs. While rapid pacing increased the number of Ca²⁺ waves in both RyR2-RS and wild-type (p<0.05), β-adrenoceptor (βAR) stimulation induced more Ca²⁺ waves in RyR2-RS (p<0.05). Notably, Ca²⁺ waves occurred despite decreased sarcoplasmic reticulum (SR) Ca²⁺ content in RyR2-RS (p<0.05), suggesting increased cytosolic RyR2 Ca²⁺ sensitivity.

A computational model of mouse ventricular cardiomyocyte electrophysiology reproduced the cellular CPVT1 phenotype when RyR2 Ca²⁺ sensitivity was increased. Importantly, diastolic fluctuations in phosphorylation of RyR2 and SR Ca²⁺ content determined Ca²⁺ wave initiation. These factors were modulated towards increased propensity for arrhythmia initiation by increased pacing rates, but even more by βAR stimulation.

Conclusion

In CPVT1, VT propensity depends on individual heart rate thresholds for PVCs. Through converging data from clinical exercise stress-testing, cellular studies and computational modelling, we confirm the heart rate-independent pro-arrhythmic effects of βAR stimulation in CPVT1, but also identify an independent and synergistic contribution from effects of high heart rate.

Investigating β-adrenergic-induced cardiac hypertrophy through computational approach: classical and non-classical pathways

The chronic stimulation of β-adrenergic receptors plays a crucial role in cardiac hypertrophy and its progression to heart failure. In β-adrenergic signaling, in addition to the well-established classical pathway, Gs/AC/cAMP/PKA, activation of non-classical pathways such as Gi/PI3K/Akt/GSK3β and Gi/Ras/Raf/MEK/ERK contribute in cardiac hypertrophy. The signaling network of β-adrenergic-induced hypertrophy is very complex and not fully understood. So, we use a computational approach to investigate the dynamic response and contribution of β-adrenergic mediators in cardiac hypertrophy. The proposed computational model provides insights into the effects of β-adrenergic classical and non-classical pathways on the activity of hypertrophic transcription factors CREB and GATA4. The results illustrate that the model captures the dynamics of the main signaling mediators and reproduces the experimental observations well. The results also show that despite the low portion of β2 receptors out of total cardiac β-adrenergic receptors, their contribution in the activation of hypertrophic mediators and regulation of β-adrenergic-induced hypertrophy is noticeable and variations in β1/β2 receptors ratio greatly affect the ISO-induced hypertrophic response. The model results illustrate that GSK3β deactivation after β-adrenergic receptor stimulation has a major influence on CREB and GATA4 activation and consequent cardiac hypertrophy. Also, it is found through sensitivity analysis that PKB (Akt) activation has both pro-hypertrophic and anti-hypertrophic effects in β-adrenergic signaling.

Saturated high-fat diet-induced obesity increases adenylate cyclase of myocardial β-adrenergic system and does not compromise cardiac function

Obesity is a worldwide pandemic associated with high incidence of cardiovascular disease. The mechanisms by which the obesity leads cardiac dysfunction are not fully elucidated and few studies have evaluated the relationship between obesity and proteins involved in myocardial β‐adrenergic (βA) system. The purpose of this study was to evaluate the cardiac function and βA pathway components in myocardium of obese rats. Male Wistar rats were distributed into two groups: control (n = 17; standard diet) and obese (n = 17; saturated high‐fat diet) fed for 33 weeks. Nutritional profile and comorbidities were assessed. Cardiac structure and function was evaluated by macroscopic postmortem, echocardiographic and isolated papillary muscle analyzes. Myocardial protein expression of β₁‐ and β₂‐adrenergic receptors, Gαs protein, adenylate cyclase (AC) and protein kinase A (PKA) was performed by Western blot. Cardiac cyclic adenosine monophosphate (cAMP) levels and PKA activity were assessed by ELISA. Obese rats showed increased adiposity index (P < 0.001) and several comorbidities as hypertension, glucose intolerance, insulin resistance, and dyslipidemia compared with control rats. Echocardiographic assessment revealed increased left atrium diameter (C: 4.98 ± 0.38 vs. Ob: 5.47 ± 0.53, P = 0.024) and posterior wall shortening velocity (C: 37.1 ± 3.6 vs. Ob: 41.8 ± 3.8, P = 0.007) in obese group. Papillary muscle evaluation indicated that baseline data and myocardial responsiveness to isoproterenol stimulation were similar between the groups. Protein expression of myocardial AC was higher in obese group than in the control (C: 1.00 ± 0.21 vs. Ob: 1.25 ± 0.10, P = 0.025), whereas the other components were unchanged. These results suggest that saturated high‐fat diet‐induced obesity was not effective in triggering cardiac dysfunction and impair the beta‐adrenergic signaling.

Mechanisms of the cyclic nucleotide cross-talk signaling network in cardiac L-type calcium channel regulation

Regulation of L-type Calcium (Ca2+) Channel (LCC) gating is critical to shaping the cardiac action potential (AP) and triggering the initiation of excitation-contraction (EC) coupling in cardiac myocytes. The cyclic nucleotide (cN) cross-talk signaling network, which encompasses the β-adrenergic and the Nitric Oxide (NO)/cGMP/Protein Kinase G (PKG) pathways and their interaction (cross-talk) through distinctively-regulated phosphodiesterase isoenzymes (PDEs), regulates LCC current via Protein Kinase A- (PKA) and PKG-mediated phosphorylation. Due to the tightly-coupled and intertwined biochemical reactions involved, it remains to be clarified how LCC gating is regulated by the signaling network from receptor to end target. In addition, the large number of EC coupling-related phosphorylation targets of PKA and PKG makes it difficult to quantify and isolate changes in L-type Ca2+ current (ICaL) responses regulated by the signaling network. We have developed a multi-scale, biophysically-detailed computational model of LCC regulation by the cN signaling network that is supported by experimental data. LCCs are modeled with functionally distinct PKA- and PKG-phosphorylation dependent gating modes. The model exhibits experimentally observed single channel characteristics, as well as whole-cell LCC currents upon activation of the cross-talk signaling network. Simulations show 1) redistribution of LCC gating modes explains changes in whole-cell current under various stimulation scenarios of the cN cross-talk network; 2) NO regulation occurs via potentiation of a gating mode characterized by prolonged closed times; and 3) due to compensatory actions of cross-talk and antagonizing functions of PKA- and PKG-mediated phosphorylation of LCCs, the effects of individual inhibitions of PDEs 2, 3, and 4 on ICaL are most pronounced at low levels of β-adrenergic stimulation. Simulations also delineate the contribution of the following two mechanisms to overall LCC regulation, which have otherwise been challenging to distinguish: 1) regulation of PKA and PKG activation via cN cross-talk (Mechanism 1); and 2) LCC interaction with activated PKA and PKG (Mechanism 2). These results provide insights into how cN signals transduced via the cN cross-talk signaling network are integrated via LCC regulation in the heart.

Over-expression of microRNA-1 causes arrhythmia by disturbing intracellular trafficking system

Dysregulation of intracellular trafficking system plays a fundamental role in the progression of cardiovascular disease. Up-regulation of miR-1 contributes to arrhythmia, we sought to elucidate whether intracellular trafficking contributes to miR-1-driven arrhythmia. By performing microarray analyses of the transcriptome in the cardiomyocytes-specific over-expression of microRNA-1 (miR-1 Tg) mice and the WT mice, we found that these differentially expressed genes in miR-1 Tg mice were significantly enrichment with the trafficking-related biological processes, such as regulation of calcium ion transport. Also, the qRT-PCR and western blot results validated that Stx6, Braf, Ube3a, Mapk8ip3, Ap1s1, Ccz1 and Gja1, which are the trafficking-related genes, were significantly down-regulated in the miR-1 Tg mice. Moreover, we found that Stx6 was decreased in the heart of mice after myocardial infarction and in the hypoxic cardiomyocytes, and further confirmed that Stx6 is a target of miR-1. Meanwhile, knockdown of Stx6 in cardiomyocytes resulted in the impairments of PLM and L-type calcium channel, which leads to the increased resting ([Ca²⁺]_i). On the contrary, overexpression of Stx6 attenuated the impairments of miR-1 or hypoxia on PLM and L-type calcium channel. Thus, our studies reveals that trafficking-related gene Stx6 may regulate intracellular calcium and is involved in the occurrence of cardiac arrhythmia, which provides new insights in that miR-1 participates in arrhythmia by regulating the trafficking-related genes and pathway.

Role of protein phosphatase inhibitor-1 in cardiac beta adrenergic pathway

Phosphoproteomic studies have shown that about one third of all cardiac proteins are reversibly phosphorylated, affecting virtually every cellular signaling pathway. The reversibility of this process is orchestrated by the opposing enzymatic activity of kinases and phosphatases. Conversely, imbalances in subcellular protein phosphorylation patterns are a hallmark of many cardiovascular diseases including heart failure and cardiac arrhythmias. While numerous studies have revealed excessive beta-adrenergic signaling followed by deregulated kinase expression or activity as a major driver of the latter cardiac pathologies, far less is known about the beta-adrenergic regulation of their phosphatase counterparts. In fact, most of the limited knowledge stems from the detailed analysis of the endogenous inhibitor of the protein phosphatase 1 (I-1) in cellular and animal models. I-1 acts as a nodal point between adrenergic and putatively non-adrenergic cardiac signaling pathways and is able to influence widespread cellular functions of protein phosphatase 1 which are contributing to cardiac health and disease, e.g. Ca²⁺ handling, sarcomere contractility and glucose metabolism. Finally, nearly all of these studies agree that I-1 is a promising drug target on the one hand but the outcome of its pharmacological regulation maybe extremely context-dependent on the other hand, thus warranting for careful interpretation of past and future experimental results. In this respect we will: 1) comprehensively review the current knowledge about structural, functional and regulatory properties of I-1 within the heart 2) highlight current working hypothesis and potential I-1 mediated disease mechanisms 3) discuss state-of-the-art knowledge and future prospects of a potential therapeutic strategy targeting I-1 by restoring the balance of cardiac protein phosphorylation.

Roles of phosphodiesterases in the regulation of the cardiac cyclic nucleotide cross-talk signaling network

The balanced signaling between the two cyclic nucleotides (cNs) cAMP and cGMP plays a critical role in regulating cardiac contractility. Their degradation is controlled by distinctly regulated phosphodiesterase isoenzymes (PDEs), which in turn are also regulated by these cNs. As a result, PDEs facilitate communication between the β-adrenergic and Nitric Oxide (NO)/cGMP/Protein Kinase G (PKG) signaling pathways, which regulate the synthesis of cAMP and cGMP respectively. The phenomena in which the cAMP and cGMP pathways influence the dynamics of each other are collectively referred to as cN cross-talk. However, the cross-talk response and the individual roles of each PDE isoenzyme in shaping this response remain to be fully characterized. We have developed a computational model of the cN cross-talk network that mechanistically integrates the β-adrenergic and NO/cGMP/PKG pathways via regulation of PDEs by both cNs. The individual model components and the integrated network model replicate experimentally observed activation-response relationships and temporal dynamics. The model predicts that, due to compensatory interactions between PDEs, NO stimulation in the presence of sub-maximal β-adrenergic stimulation results in an increase in cytosolic cAMP accumulation and corresponding increases in PKA-I and PKA-II activation; however, the potentiation is small in magnitude compared to that of NO activation of the NO/cGMP/PKG pathway. In a reciprocal manner, β-adrenergic stimulation in the presence of sub-maximal NO stimulation results in modest cGMP elevation and corresponding increase in PKG activation. In addition, we demonstrate that PDE2 hydrolyzes increasing amounts of cAMP with increasing levels of β-adrenergic stimulation, and hydrolyzes increasing amounts of cGMP with decreasing levels of NO stimulation. Finally, we show that PDE2 compensates for inhibition of PDE5 both in terms of cGMP and cAMP dynamics, leading to cGMP elevation and increased PKG activation, while maintaining whole-cell β-adrenergic responses similar to that prior to PDE5 inhibition. By defining and quantifying reactions comprising cN cross-talk, the model characterizes the cross-talk response and reveals the underlying mechanisms of PDEs in this non-linear, tightly-coupled reaction system.

Interaction between phosphodiesterases in the regulation of the cardiac β-adrenergic pathway

In cardiac myocytes, the second messenger cAMP is synthesized within the β-adrenergic signaling pathway upon sympathetic activation. It activates Protein Kinase A (PKA) mediated phosphorylation of multiple target proteins that are functionally critical to cardiac contractility. The dynamics of cAMP are also controlled indirectly by cGMP-mediated regulation of phosphodiesterase isoenzymes (PDEs). The nature of the interactions between cGMP and the PDEs, as well as between PDE isoforms, and how these ultimately transduce the cGMP signal to regulate cAMP remains unclear. To better understand this, we have developed mechanistically detailed models of PDEs 1-4, the primary cAMP-hydrolyzing PDEs in cardiac myocytes, and integrated them into a model of the β-adrenergic signaling pathway. The PDE models are based on experimental studies performed on purified PDEs which have demonstrated that cAMP and cGMP bind competitively to the cyclic nucleotide (cN)-binding domains of PDEs 1, 2, and 3, while PDE4 regulation occurs via PKA-mediated phosphorylation. Individual PDE models reproduce experimentally measured cAMP hydrolysis rates with dose-dependent cGMP regulation. The fully integrated model replicates experimentally observed whole-cell cAMP activation-response relationships and temporal dynamics upon varying degrees of β-adrenergic stimulation in cardiac myocytes. Simulations reveal that as a result of network interactions, reduction in the level of one PDE is partially compensated for by increased activation of others. PDE2 and PDE4 exert the strongest compensatory roles among all PDEs. In addition, PDE2 competes with other PDEs to bind and hydrolyze cAMP and is a strong regulator of PDE interactions. Finally, an increasing level of cGMP gradually out-competes cAMP for the catalytic sites of PDEs 1, 2, and 3, suppresses their cAMP hydrolysis rates, and results in amplified cAMP signaling. These results provide insights into how PDEs transduce cGMP signals to regulate cAMP and how PDE interactions affect cardiac β-adrenergic response.

Cardiac Dysfunction Induced by Obesity Is Not Related to β-Adrenergic System Impairment at the Receptor-Signalling Pathway

Obesity has been shown to impair myocardial performance. Some factors have been suggested as responsible for possible cardiac abnormalities in models of obesity, among them beta-adrenergic (βA) system, an important mechanism of regulation of myocardial contraction and relaxation. The objective of present study was to evaluate the involvement of βA system components in myocardial dysfunction induced by obesity. Thirty-day-old male Wistar rats were distributed in control (C, n = 25) and obese (Ob, n = 25) groups. The C group was fed a standard diet and Ob group was fed four unsaturated high-fat diets for 15 weeks. Cardiac function was evaluated by isolated papillary muscle preparation and βA system evaluated by using cumulative concentrations of isoproterenol and Western blot. After 15 weeks, the Ob rats developed higher adiposity index than C rats and several comorbidities; however, were not associated with changes in systolic blood pressure. Obesity caused structural changes and the myocardial responsiveness to post-rest contraction stimulus and increased extracellular calcium (Ca²⁺) was compromised. There were no changes in cardiac function between groups after βA stimulation. The obesity was not accompanied by changes in protein expression of G protein subunit alpha (Gsα) and βA receptors (β₁AR and β₂AR). In conclusion, the myocardial dysfunction caused by unsaturated high-fat diet-induced obesity, after 15 weeks, is not related to βAR system impairment at the receptor-signalling pathway.

Modeling the effects of β1-adrenergic receptor blockers and polymorphisms on cardiac myocyte Ca2+ handling

β-Adrenergic receptor blockers (β-blockers) are commonly used to treat heart failure, but the biologic mechanisms governing their efficacy are still poorly understood. The complexity of β-adrenergic signaling coupled with the influence of receptor polymorphisms makes it difficult to intuit the effect of β-blockers on cardiac physiology. While some studies indicate that β-blockers are efficacious by inhibiting β-adrenergic signaling, other studies suggest that they work by maintaining β-adrenergic responsiveness. Here, we use a systems pharmacology approach to test the hypothesis that in ventricular myocytes, these two apparently conflicting mechanisms for β-blocker efficacy can occur concurrently. We extended a computational model of the β(1)-adrenergic pathway and excitation-contraction coupling to include detailed receptor interactions for 19 ligands. Model predictions, validated with Ca(2+) and Förster resonance energy transfer imaging of adult rat ventricular myocytes, surprisingly suggest that β-blockers can both inhibit and maintain signaling depending on the magnitude of receptor stimulation. The balance of inhibition and maintenance of β(1)-adrenergic signaling is predicted to depend on the specific β-blocker (with greater responsiveness for metoprolol than carvedilol) and β(1)-adrenergic receptor Arg389Gly polymorphisms.

Fibroblast growth factor receptor 1 signaling in adult cardiomyocytes increases contractility and results in a hypertrophic cardiomyopathy

Fibroblast growth factors (FGFs) and their receptors are highly conserved signaling molecules that have been implicated in postnatal cardiac remodeling. However, it is not known whether cardiomyocyte-expressed FGF receptors are necessary or sufficient for ventricular remodeling in the adult heart. To determine whether cardiomyocytes were competent to respond to an activated FGF receptor, and to determine if this signal would result in the development of hypertrophy, we engineered a doxycycline (DOX)-inducible, cardiomyocyte-specific, constitutively active FGF receptor mouse model (αMHC-rtTA, TRE-caFgfr1-myc). Echocardiographic and hemodynamic analysis indicated that acute expression of caFGFR1 rapidly and directly increased cardiac contractility, while chronic expression resulted in significant hypertrophy with preservation of systolic function. Subsequent histologic analysis showed increased cardiomyocyte cross-sectional area and regions of myocyte disarray and fibrosis, classic features of hypertrophic cardiomyopathy (HCM). Analysis of downstream pathways revealed a lack of clear activation of classical FGF-mediated signaling pathways, but did demonstrate a reduction in Serca2 expression and troponin I phosphorylation. Isolated ventricular myocytes showed enhanced contractility and reduced relaxation, an effect that was partially reversed by inhibition of actin-myosin interactions. We conclude that adult cardiomyocytes are competent to transduce FGF signaling and that FGF signaling is sufficient to promote increased cardiomyocyte contractility in vitro and in vivo through enhanced intrinsic actin-myosin interactions. Long-term, FGFR overexpression results in HCM with a dynamic outflow tract obstruction, and may serve as a unique model of HCM.

PKA catalytic subunit compartmentation regulates contractile and hypertrophic responses to β-adrenergic signaling

β-Adrenergic signaling is spatiotemporally heterogeneous in the cardiac myocyte, conferring exquisite control to sympathetic stimulation. Such heterogeneity drives the formation of protein kinase A (PKA) signaling microdomains, which regulate Ca(2+) handling and contractility. Here, we test the hypothesis that the nucleus independently comprises a PKA signaling microdomain regulating myocyte hypertrophy. Spatially-targeted FRET reporters for PKA activity identified slower PKA activation and lower isoproterenol sensitivity in the nucleus (t50=10.6±0.7 min; EC50=89.0 nmol/L) than in the cytosol (t50=3.71±0.25 min; EC50=1.22 nmol/L). These differences were not explained by cAMP or AKAP-based compartmentation. A computational model of cytosolic and nuclear PKA activity was developed and predicted that differences in nuclear PKA dynamics and magnitude are regulated by slow PKA catalytic subunit diffusion, while differences in isoproterenol sensitivity are regulated by nuclear expression of protein kinase inhibitor (PKI). These were validated by FRET and immunofluorescence. The model also predicted differential phosphorylation of PKA substrates regulating cell contractility and hypertrophy. Ca(2+) and cell hypertrophy measurements validated these predictions and identified higher isoproterenol sensitivity for contractile enhancements (EC50=1.84 nmol/L) over cell hypertrophy (EC50=85.9 nmol/L). Over-expression of spatially targeted PKA catalytic subunit to the cytosol or nucleus enhanced contractile and hypertrophic responses, respectively. We conclude that restricted PKA catalytic subunit diffusion is an important PKA compartmentation mechanism and the nucleus comprises a novel PKA signaling microdomain, insulating hypertrophic from contractile β-adrenergic signaling responses.

Beta-adrenergic stimulation maintains cardiac function in Serca2 knockout mice

Previous studies on Serca2 knockout (KO) mice showed that cardiac function is sustained in vivo for several weeks after knockout, whereas SERCA protein levels decrease and calcium dynamics are significantly impaired. In this study, we reconcile observed cellular and organ level contractile function using a cardiac multiscale model. We identified and quantified the changes in cellular function that are both consistent with observations and able to compensate for the decrease in SERCA. Calcium transients were used as input for multiscale computational simulations to predict whole-organ response. Although this response matched experimental pressure-volume (PV) measurements in healthy mice, the reduced magnitude calcium transients observed in KO cells were insufficient to trigger ventricular ejection. To replicate the effects of elevated catecholamine levels observed in vivo, cells were treated with isoproterenol. Incorporation of the resulting measured β-adrenergically stimulated calcium transients into the model resulted in a close match with experimental PV loops. Changes in myofilament properties, when considered in isolation, were not able to increase tension development to levels consistent with measurements, further confirming the necessity of a high β-adrenergic state. Modeling additionally indicated that increased venous return observed in the KO mice helps maintain a high ejection fraction via the Frank-Starling effect. Our study shows that increased β-adrenergic stimulation is a potentially highly significant compensatory mechanism by which cardiac function is maintained in Serca2 KO mice, producing the increases in both systolic and diastolic calcium, consistent with the observed contractile function observed in experimental PV measurements.

Genome-wide association study of cardiac structure and systolic function in African Americans: the Candidate Gene Association Resource (CARe) study

BACKGROUND

Using data from 4 community-based cohorts of African Americans, we tested the association between genome-wide markers (single-nucleotide polymorphisms) and cardiac phenotypes in the Candidate-gene Association Resource study.

METHODS AND RESULTS

Among 6765 African Americans, we related age, sex, height, and weight-adjusted residuals for 9 cardiac phenotypes (assessed by echocardiogram or magnetic resonance imaging) to 2.5 million single-nucleotide polymorphisms genotyped using Genome-wide Affymetrix Human SNP Array 6.0 (Affy6.0) and the remainder imputed. Within the cohort, genome-wide association analysis was conducted, followed by meta-analysis across cohorts using inverse variance weights (genome-wide significance threshold=4.0 ×10(-7)). Supplementary pathway analysis was performed. We attempted replication in 3 smaller cohorts of African ancestry and tested lookups in 1 consortium of European ancestry (EchoGEN). Across the 9 phenotypes, variants in 4 genetic loci reached genome-wide significance: rs4552931 in UBE2V2 (P=1.43×10(-7)) for left ventricular mass, rs7213314 in WIPI1 (P=1.68×10(-7)) for left ventricular internal diastolic diameter, rs1571099 in PPAPDC1A (P=2.57×10(-8)) for interventricular septal wall thickness, and rs9530176 in KLF5 (P=4.02×10(-7)) for ejection fraction. Associated variants were enriched in 3 signaling pathways involved in cardiac remodeling. None of the 4 loci replicated in cohorts of African ancestry was confirmed in lookups in EchoGEN.

CONCLUSIONS

In the largest genome-wide association study of cardiac structure and function to date in African Americans, we identified 4 genetic loci related to left ventricular mass, interventricular septal wall thickness, left ventricular internal diastolic diameter, and ejection fraction, which reached genome-wide significance. Replication results suggest that these loci may be unique to individuals of African ancestry. Additional large-scale studies are warranted for these complex phenotypes.

Ca(2+)-dependent large conductance K(+) currents in thalamocortical relay neurons of different rat strains

Mutations in genes coding for Ca(2+) channels were found in patients with childhood absence epilepsy (CAE) indicating a contribution of Ca(2+)-dependent mechanisms to the generation of spike-wave discharges (SWD) in humans. Since the involvement of Ca(2+) signals remains unclear, the aim of the present study was to elucidate the function of a Ca(2+)-dependent K(+) channel (BKCa) under physiological conditions and in the pathophysiological state of CAE. The activation of BKCa channels is dependent on both voltage and intracellular Ca(2+) concentrations. Moreover, these channels exhibit an outstandingly high level of regulatory heterogeneity that builds the basis for the influence of BKCa channels on different aspects of neuronal activity. Here, we analyse the contribution of BKCa channels to firing of thalamocortical relay neurons, and we test the hypothesis that BKCa channel activity affects the phenotype of a genetic rat model of CAE. We found that the activation of the β2-adrenergic receptor/protein kinase A pathway resulted in BKCa channel inhibition. Furthermore, BKCa channels affect the number of action potentials fired in a burst and produced spike frequency adaptation during tonic activity. The latter result was confirmed by a computer modelling approach. We demonstrate that the β2-adrenergic inhibition of BKCa channels prevents spike frequency adaptation and, thus, might significantly support the tonic firing mode of thalamocortical relay neurons. In addition, we show that BKCa channel functioning differs in epileptic WAG/Rij and thereby likely contributes to highly synchronised, epileptic network activity.

β-Adrenergic cAMP signals are predominantly regulated by phosphodiesterase type 4 in cultured adult rat aortic smooth muscle cells

BACKGROUND

We investigated the role of cyclic nucleotide phosphodiesterases (PDEs) in the spatiotemporal control of intracellular cAMP concentrations in rat aortic smooth muscle cells (RASMCs).

METHODOLOGY/PRINCIPAL FINDINGS

The rank order of PDE families contributing to global cAMP-PDE activity was PDE4> PDE3  =  PDE1. PDE7 mRNA expression but not activity was confirmed. The Fluorescence Resonance Energy Transfer (FRET)-based cAMP sensor, Epac1-camps, was used to monitor the time course of cytosolic cAMP changes. A pulse application of the β-adrenoceptor (β-AR) agonist isoproterenol (Iso) induced a transient FRET signal. Both β(1)- and β(2)-AR antagonists decreased the signal amplitude without affecting its kinetics. The non-selective PDE inhibitor (IBMX) dramatically increased the amplitude and delayed the recovery phase of Iso response, in agreement with a role of PDEs in degrading cAMP produced by Iso. Whereas PDE1, PDE3 and PDE7 blockades [with MIMX, cilostamide (Cil) and BRL 50481 (BRL), respectively] had no or minor effect on Iso response, PDE4 inhibition [with Ro-20-1724 (Ro)] strongly increased its amplitude and delayed its recovery. When Ro was applied concomitantly with MIMX or Cil (but not with BRL), the Iso response was drastically further prolonged. PDE4 inhibition similarly prolonged both β(1)- and β(2)-AR-mediated responses. When a membrane-targeted FRET sensor was used, PDE3 and PDE4 acted in a synergistic manner to hydrolyze the submembrane cAMP produced either at baseline or after β-AR stimulation.

CONCLUSION/SIGNIFICANCE

Our study underlines the importance of cAMP-PDEs in the dynamic control of intracellular cAMP signals in RASMCs, and demonstrates the prominent role of PDE4 in limiting β-AR responses. PDE4 inhibition unmasks an effect of PDE1 and PDE3 on cytosolic cAMP hydrolyzis, and acts synergistically with PDE3 inhibition at the submembrane compartment. This suggests that mixed PDE4/PDE1 or PDE4/PDE3 inhibitors would be attractive to potentiate cAMP-related functions in vascular cells.

Phospholemman is a negative feed-forward regulator of Ca2+ in β-adrenergic signaling, accelerating β-adrenergic inotropy

Sympathetic stimulation enhances cardiac contractility by stimulating β-adrenergic signaling and protein kinase A (PKA). Recently, phospholemman (PLM) has emerged as an important PKA substrate capable of regulating cytosolic Ca(2+) transients. However, it remains unclear how PLM contributes to β-adrenergic inotropy. Here we developed a computational model to clarify PLM's role in the β-adrenergic signaling response. Simulating Na(+) and sarcoplasmic reticulum (SR) Ca(2+) clamps, we identify an effect of PLM phosphorylation on SR unloading as the key mechanism by which PLM confers cytosolic Ca(2+) adaptation to long-term β-adrenergic receptor (β-AR) stimulation. Moreover, we show that phospholamban (PLB) opposes and overtakes these actions on SR load, forming a negative feed-forward loop in the β-adrenergic signaling cascade. This network motif dominates the negative feedback conferred by β-AR desensitization and accelerates β-AR-induced inotropy. Model analysis therefore unmasks key actions of PLM phosphorylation during β-adrenergic signaling, indicating that PLM is a critical component of the fight-or-flight response.

Role of cardiac myofilament proteins titin and collagen in the pathogenesis of diastolic dysfunction in cirrhotic rats

BACKGROUND & AIMS

Significance of diastolic dysfunction in cirrhotic cardiomyopathy has been brought to the forefront with several reports of unexpected heart failure following liver transplantation and transjugular intrahepatic portosystemic stent-shunt, but the etiology remains unclear. The present study investigated the role of passive tension regulators - titin and collagen - in the pathogenesis of this condition.

METHODS

Cirrhosis was induced by bile duct ligation (BDL) in rats, while controls underwent bile duct inspection with no ligation. Four weeks after operation, cardiac mRNA and protein levels of titin, collagen, and protein kinase A (PKA) were determined. Diastolic function was examined in isolated right ventricular cardiomyocytes, while passive tension was examined in right ventricular trabeculae muscles.

RESULTS

In BDL animals, diastolic return velocity was significantly decreased, relaxation time increased and passive tension increased. However, no significant difference in mRNA and protein levels of titin was observed. PKA mRNA and protein levels were significantly decreased in BDL animals. Collagen levels were also significantly altered in the BDL group.

CONCLUSIONS

Therefore, diastolic dysfunction exists in cirrhosis with alterations in titin modulation, PKA levels, and collagen configuration contributing to the pathogenesis of this condition.

Robustness portraits of diverse biological networks conserved despite order-of-magnitude parameter uncertainty

MOTIVATION

Biological networks are robust to a wide variety of internal and external perturbations, yet fragile or sensitive to a small minority of perturbations. Due to this rare sensitivity of networks to certain perturbations, it is unclear how precisely biochemical parameters must be experimentally measured in order to accurately predict network function.

RESULTS

Here, we examined a model of cardiac β-adrenergic signaling and found that its robustness portrait, a global measure of steady-state network function, was well conserved even when all parameters were rounded to their nearest 1-2 orders of magnitude. In contrast, β-adrenergic network kinetics were more sensitive to parameter precision. This analysis was then extended to 10 additional networks, including Escherichia coli chemotaxis, stem cell differentiation and cytokine signaling, of which nine exhibited conserved robustness portraits despite the order-of-magnitude approximation of their biochemical parameters. Thus, both fragile and robust aspects of diverse biological networks are largely shaped by network topology and can be predicted despite order-of-magnitude uncertainty in biochemical parameters. These findings suggest an iterative strategy where order-of-magnitude models are used to prioritize experiments toward the fragile network elements that require precise measurements, efficiently driving model revision.

CONTACT

jsaucerman@virginia.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cardiac myocytes and local signaling in nano-domains

It is well known that calcium-induced calcium-release in cardiac myocytes takes place in spatially restricted regions known as dyads, where discrete patches of junctional sarcoplasmic reticulum tightly associate with the t-tubule membrane. The dimensions of a dyad are so small that it contains only a few Ca²⁺ ions at any given time. Ca²⁺ signaling in the dyad is therefore noisy, and dominated by the Brownian motion of Ca²⁺ ions in a potential field. Remarkably, from this complexity emerges the integrated behavior of the myocyte in which, under normal conditions, precise control of Ca²⁺ release and muscle contraction is maintained over the life of the cell. This is but one example of how signal processing within the cardiac myocyte and other cells often occurs in small "nano-domains" where proteins and protein complexes interact at spatial dimensions on the order of ∼1-10 nm and at time-scales on the order of nanoseconds to perform the functions of the cell. In this article, we will review several examples of local signaling in nano-domains, how it contributes to the integrative behavior of the cardiac myocyte, and present computational methods for modeling signal processing within these domains across differing spatio-temporal scales.

Computational models reduce complexity and accelerate insight into cardiac signaling networks

Cardiac signaling networks exhibit considerable complexity in size and connectivity. The intrinsic complexity of these networks complicates the interpretation of experimental findings. This motivates new methods for investigating the mechanisms regulating cardiac signaling networks and the consequences these networks have on cardiac physiology and disease. Next-generation experimental techniques are also generating a wealth of genomic and proteomic data that can be difficult to analyze or interpret. Computational models are poised to play a key role in addressing these challenges. Computational models have a long history in contributing to the understanding of cardiac physiology and are useful for identifying biological mechanisms, inferring multiscale consequences to cell signaling activities and reducing the complexity of large data sets. Models also integrate well with experimental studies to explain experimental observations and generate new hypotheses. Here, we review the contributions computational modeling approaches have made to the analysis of cardiac signaling networks and forecast opportunities for computational models to accelerate cardiac signaling research.

Cardiac and vascular gene profiles in an animal model of takotsubo cardiomyopathy

We investigated cardiac and vascular gene profiles in response to immobilization stress (IMO) in rats, an animal model of emotional stress-induced takotsubo cardiomyopathy using microarray analysis, followed by re-confirmation with real-time reverse transcription-polymerase chain reaction. Expression levels of the identified genes were further estimated by pretreatment with an α1-adrenoceptor blocker and/or a β1-adrenoceptor blocker. In response to IMO, expression of 46 genes was significantly altered in the heart and that of 49 genes was significantly altered in the aorta. Pathway analysis with DAVID Bioinformatics Resources indicated that regulation of transcription and response to endogenous stimulation were the top two scoring pathways. Altered expression of cardiac genes was blunted by pretreatment with a β1-adrenoceptor blocker or α1 + β1-adrenoceptor blockers. In contrast, that of aortic genes was blunted by pretreatment with an α1-adrenoceptor blocker or α1 + β1-adrenoceptor blockers. Activation of α1-adrenoceptor in the blood vessels or activation of β1-adrenoceptors in the heart were mainly responsible for emotional stress-induced alteration of cardiac and vascular gene profiles.

Modeling cardiac β-adrenergic signaling with normalized-Hill differential equations: comparison with a biochemical model

Background

New approaches are needed for large-scale predictive modeling of cellular signaling networks. While mass action and enzyme kinetic approaches require extensive biochemical data, current logic-based approaches are used primarily for qualitative predictions and have lacked direct quantitative comparison with biochemical models.

Results

We developed a logic-based differential equation modeling approach for cell signaling networks based on normalized Hill activation/inhibition functions controlled by logical AND and OR operators to characterize signaling crosstalk. Using this approach, we modeled the cardiac β₁-adrenergic signaling network, including 36 reactions and 25 species. Direct comparison of this model to an extensively characterized and validated biochemical model of the same network revealed that the new model gave reasonably accurate predictions of key network properties, even with default parameters. Normalized Hill functions improved quantitative predictions of global functional relationships compared with prior logic-based approaches. Comprehensive sensitivity analysis revealed the significant role of PKA negative feedback on upstream signaling and the importance of phosphodiesterases as key negative regulators of the network. The model was then extended to incorporate recently identified protein interaction data involving integrin-mediated mechanotransduction.

Conclusions

The normalized-Hill differential equation modeling approach allows quantitative prediction of network functional relationships and dynamics, even in systems with limited biochemical data.

Models of cardiac excitation-contraction coupling in ventricular myocytes

Mathematical and computational modeling of cardiac excitation-contraction coupling has produced considerable insights into how the heart muscle contracts. With the increase in biophysical and physiological data available, the modeling has become more sophisticated with investigations spanning in scale from molecular components to whole cells. These modeling efforts have provided insight into cardiac excitation-contraction coupling that advanced and complemented experimental studies. One goal is to extend these detailed cellular models to model the whole heart. While this has been done with mechanical and electrophysiological models, the complexity and fast time course of calcium dynamics have made inclusion of detailed calcium dynamics in whole heart models impractical. Novel methods such as the probability density approach and moment closure technique which increase computational efficiency might make this tractable.

Pharmacological PKA inhibition: all may not be what it seems

Signaling through the cyclic adenosine monophosphate-dependent protein kinase [protein kinase A (PKA)] is an important and widely studied area of signal transduction research. This signaling pathway is commonly investigated through the use of the pharmacological PKA inhibitors H89 and KT 5720. Both of these compounds are thought to block PKA actions through competitive inhibition of the adenosine triphosphate site on the PKA catalytic subunit. Recently, a number of studies have identified actions of H89 and KT 5720 that are independent of their effects on PKA. These nonspecific effects are widespread; they include actions on other protein kinases and signaling molecules and also on basic cellular functions, such as transcription. Here, I summarize the nonspecific effects of these two compounds and compare their actions with those of other PKA inhibitors.

Roles of PKA, PI3K, and cPLA2 in the NO-mediated negative inotropic effect of β2-adrenoceptor agonists in guinea pig right papillary muscles

Although β2-adrenoceptors represent 15–25% of β-adrenoceptors in the guinea pig heart, their functionality is controversial. We assessed the inotropic effects of β2-adrenoceptor partial agonists in right papillary muscles. Salbutamol induced a small but significant concentration-dependent negative inotropic effect (NIE, −5% at 60 nM) followed by a moderate positive inotropic effect (+36% at 6 μM) due to activation of β1-adrenoceptors. In the presence of 4 μM atenolol, the concentration-dependent NIE (−12% at 6 μM) was biphasic, best described by a double logistic equation with respective EC50 values of 3 and ∼420 nM, and was insensitive to SR59230A. In muscles from pertussis toxin-treated guinea pigs, the salbutamol-induced positive inotropic effect was sensitive to low concentrations of ICI-118551 in an unusual manner. Experiments in reserpinized animals revealed the importance of the phosphorylation-dephosphorylation processes. PKA inhibition reduced and suppressed the effects obtained at low and high concentrations, respectively, indicating that its activation was a prerequisite to the NIE. The effect occurring at nanomolar concentrations depended upon PKA/phosphatidylinositol 3-kinase/cytosolic phospholipase A2 (cPLA2) activations leading to nitric oxide (NO) release via the arachidonic acid/cyclooxygenase pathway. NO release via PKA-dependent phosphorylation of the receptor was responsible for the inotropic effect observed at submicromolar concentrations, which is negatively controlled by cPLA2. The possibility that these effects are due to an equilibrium between different affinity states of the receptor (Gs/Gi coupled and Gi independent with different signaling pathways) that can be displaced by ICI-118551 is discussed. We conclude that β2-adrenoceptors are functional in guinea pig heart and can modulate the inotropic state.

Epinephrine regulation of the endothelial nitric-oxide synthase: roles of RAC1 and beta3-adrenergic receptors in endothelial NO signaling

beta-Adrenergic receptors (betaAR) play an important role in vasodilation, but the mechanisms whereby adrenergic pathways regulate the endothelial isoform of nitric-oxide synthase (eNOS) are incompletely understood. We found that epinephrine significantly increases eNOS activity in cultured bovine aortic endothelial cells (BAEC). Epinephrine-dependent eNOS activation was accompanied by an increase in phosphorylation of eNOS at Ser(1179) and with decreased eNOS phosphorylation at the inhibitory phosphoresidues Ser(116) and Thr(497). Epinephrine promoted activation of the small G protein Rac1 and also led to the activation of protein kinase A. All of these responses to epinephrine in BAEC were blocked by the beta(3)AR blocker SR59230A. We transfected and validated duplex small interfering RNA (siRNA) constructs to selectively "knock down" specific signaling proteins in BAEC. siRNA-mediated knockdown of Rac1 completely blocked all beta(3)AR signaling to eNOS and also abrogated epinephrine-dependent cAMP-dependent protein kinase (PKA) and Akt activation. However, siRNA-mediated knockdown of PKA did not affect Rac1 activation by epinephrine but did attenuate Akt activation by epinephrine. These findings indicate that Rac1 is an upstream regulator of beta(3)AR signaling to PKA and to eNOS and identify a novel beta(3)AR --> Rac1 --> PKA --> Akt pathway in endothelium. We exploited the p21-activated kinase pulldown assay to identify proteins associated with activated Rac1 and found that epinephrine stimulated the association of eNOS with Rac1; epinephrine-stimulated eNOS-Rac1 interactions were blocked by the beta(3)AR antagonist SR59230A. Co-transfection of eNOS cDNA with constitutively active Rac1 enhanced beta(3)AR-promoted eNOS-Rac1 association; co-transfection of eNOS with dominant negative Rac1 completely blocked the eNOS-Rac1 association. We also found that epinephrine-induced Rac1 --> PKA --> Akt pathway mediates beta(3)AR-mediated endothelial cell migration. Taken together, our data establish that the small G protein Rac1 is a key regulator of beta(3)AR signaling in cultured aortic endothelial cells with potentially important implications for the pathways involved in adrenergic modulation of eNOS pathways in the vascular wall.