- Short Communication
- Open Access
Oscillation of cAMP and Ca2+ in cardiac myocytes: a systems biology approach
The Journal of Physiological Sciences volume 65, pages195–200(2015)
Cyclic adenosine monophosphate (cAMP) and Ca2+ levels may oscillate in harmony within excitable cells; a mathematical oscillation loop model, the Cooper model, of these oscillations was developed two decades ago. However, in that model all adenylyl cyclase (AC) isoforms were assumed to be inhibited by Ca2+, and it is now known that the heart expresses multiple AC isoforms, among which the type 5/6 isoforms are Ca2+-inhibitable whereas the other five (AC2, 3, 4, 7, and 9) are not. We used a computational systems biology approach with CellDesigner simulation software to develop a comprehensive graphical map and oscillation loop model for cAMP and Ca2+. This model indicated that Ca2+-mediated inhibition of AC is essential to create oscillations of Ca2+ and cAMP, and the oscillations were not altered by incorporation of phosphodiesterase-mediated cAMP hydrolysis or PKA-mediated inhibition of AC into the model. More importantly, they were created but faded out immediately in the co-presence of Ca2+-noninhibitable AC isoforms. Because the subcellular locations of AC isoforms are different, spontaneous cAMP and Ca2+ oscillations may occur within microdomains containing only Ca2+-inhibitable isoforms in cardiac myocytes, which might be necessary for fine tuning of excitation–contraction coupling.
β-Adrenergic receptor (β-AR) signaling is of crucial importance in regulating normal cardiac function, and abnormality of β-AR signaling contributes to the development of heart failure via altered cyclic AMP (cAMP) and calcium (Ca2+) signaling [1–3]. Cardiac excitation–contraction coupling (E–C coupling) is the process that links electrical excitation of cardiac myocytes to contraction of heart muscle. Ca2+ is essential for cardiac electrical activity and is a direct activator of myofilaments, causing both contraction and relaxation [4–6]. Dysregulation of cAMP and the subsequent Ca2+ oscillation are fundamental causes of both contractile and diastolic dysfunction, and arrhythmia among heart failure patients [1, 3, 5, 7–10].
Adenylyl cyclase (AC) is a membrane-bound enzyme that catalyzes the conversion of ATP to cAMP [1, 11]. cAMP, an intracellular second messenger, activates protein kinase A (PKA), leading to phosphorylation of multiple molecules involved in cardiac contraction, including the L-type Ca2+-channel . Phosphorylation of the L-type Ca2+-channel is known to increase the influx of Ca2+, resulting in increased intracellular levels of Ca2+. In the 1990s, AC isoforms directly inhibited by Ca2+ were identified, i.e., types 5 and 6 AC isoforms (AC5 and AC6) [11–14]. Because these isoforms are dominantly expressed in the heart, it was proposed that cAMP levels may oscillate in harmony with Ca2+ levels; an increase in cAMP, as generated by AC5/6, phosphorylates L-type Ca2+ channels, and induces influx of Ca2+ into the cytosol . An increase of cytosolic Ca2+ inhibits AC5/6 and reduces phosphorylation of the L-type Ca2+ channel. Thus, the activity of AC5/6 and the L-type Ca2+ channel may work synergistically to generate an oscillation loop of cAMP and Ca2+ in cardiac myocytes .
It is now well known that the heart expresses not only AC5 and AC6, but also many other AC isoforms (AC2, 3, 4, 5, 6, 7, and 9) . Although these isoforms are all expressed in the heart, recent studies have indicated they may have different subcellular locations [15, 16]. AC5 is a major cardiac isoform in adults, and AC6 is a fetal or neonatal cardiac AC isoform [1, 17–19]. AC5 and AC6 share most, if not all, of their biochemical properties, and are inhibited not only by Ca2+, but also by Gi and PKA [2, 11, 13, 20–22]. In contrast, the other AC isoforms (AC2, 3, 4, 7, and 9), which are ubiquitously expressed throughout the body, are not inhibited by Ca2+, Gi, or PKA , and their involvement in cAMP and Ca2+ oscillations is poorly understood.
We therefore examined the involvement of AC isoforms in cAMP and Ca2+ oscillations [1, 12]. Because it is difficult to examine these issues by means of traditional in vitro or in vivo biochemical approaches, we used a computational systems biology approach with CellDesigner software, a recently developed, structure diagram editor for drawing gene-regulatory and biochemical networks [23–25].
Materials and methods
CellDesigner version 4.2 (http://www.celldesigner.org/) enables users to describe molecular interactions by using well-defined and consistent graphical notions and to create a comprehensive model incorporating positive feedforward or negative feedback loops among AC, cAMP, Ca2+, phosphodiesterase (PDE), and PKA within the β-AR signaling pathway [23–25].
The CellDesigner notation used in this paper is briefly illustrated for a simple reaction scheme in Supplemental Fig. 1. Protein A is transformed to protein B and protein C promotes this transition (Supplemental Fig. 1a). Supplemental Fig. 1b shows the notation for degradation (upper) or production (lower) of protein A; their balance determines the concentration of protein A in cells under physiological conditions.
In this study, most of the variables were the same as in Cooper’s original model ; in future work, it would be desirable to optimize the variables used in the oscillation loop model to match physiological conditions.
The formulas and values used to generate the oscillation models shown in the figures can be found in the online Supplemental methods and Supplemental Tables 1–3, available on http://www.link.springer.com/journal/12576.
Cooper’s model mimicked by CellDesigner
We first mimicked Cooper’s original model by using CellDesigner to create a graphical comprehensive map (Fig. 1a) and oscillation loop model of cAMP and Ca2+ (Fig. 1b). We also incorporated the activity of PKA and AC. For AC, we used AC5/6 because they are the major cardiac isoforms and are directly inhibited by Ca2+ [2, 13]. We obtained stable and spontaneous oscillation curves for cAMP and Ca2+, as demonstrated in the original model . The activity of AC and PKA also oscillated (Fig. 1b).
The formulas and values used to generate this oscillation model are shown in online Supplemental methods and Supplemental Table 1.
Incorporation of the PDE-mediated cAMP hydrolysis
Intracellular cAMP concentration is determined by the balance between its production via AC and its hydrolysis via PDE in the heart under physiological and pathological conditions . Because PDE is activated by cAMP, we incorporated its activity into the model (Supplemental Fig. 2a). We found that oscillations of the molecules involved in the β-AR signaling pathway were maintained, even though their amplitudes were increased by approximately 1.6-fold, compared with those of Cooper’s original model (Supplemental Fig. 2b). Thus, PDE-mediated cAMP hydrolysis did not appear to change the behavior of the oscillation, but exaggerated its amplitude.
The formulas and values used to generate this oscillation model are shown in online Supplemental methods and Supplemental Table 2.
Effect of Ca2+-mediated inhibition of AC
The heart expresses multiple AC isoforms (AC2, 3, 4, 5, 6, 7, and 9) , of which AC5 and AC6 are directly inhibited by submicromolar Ca2+ [2, 13]. Thus, to examine the effect of Ca2+-mediated inhibition of AC, we modeled the situation in which all AC isoforms in the heart are not Ca2+-inhibitable in the heart (Fig. 2a). As shown in Fig. 2b, we found that no oscillation appeared. This result indicates that the presence of Ca2+-inhibitable AC isoforms is essential for stable and spontaneous cAMP and Ca2+ oscillations to occur.
The formulas and values used to generate this oscillation model are shown in online supplemental methods and supplemental Table 1.
Incorporation of PKA-mediated inhibition of cardiac AC isoforms into the model
Recent studies have indicated that AC5 and AC6 are inhibited not only by Ca2+, but also by PKA [13, 20–22]. Therefore, we next incorporated PKA-mediated inhibition into Cooper’s model (Fig. 3a). Oscillations of cAMP and Ca2+, as well as AC and PKA, were observed, but the amplitudes were reduced by approximately 13 %, compared with those in Cooper’s original model (Fig. 3b). Thus, PKA-mediated inhibition of cardiac AC isoforms did not seem to change the behavior of the oscillation, but reduced its amplitude.
The formulas and values used to generate this oscillation model are shown in online Supplemental methods and Supplemental Table 1.
Incorporation of Ca2+-mediated inhibition of type 5/6 AC isoforms into the model
We then examined the model incorporating both Ca2+-inhibitable (AC5/6) and non-inhibitable AC isoforms (AC2, 3, 4, 7, and 9) (Fig. 4a). Studies with transgenic mouse models in vivo have demonstrated that AC2, 3, 4, 7, and 9 contribute significantly to the total AC activity in the heart [2, 27], but, unlike AC5/6, are not subjected to PKA-mediated inhibition . Oscillations of cAMP, Ca2+, PKA, PDE, and AC5/6 were observed, and then faded. However, essentially, no oscillation of AC4/7 occurred (Fig. 4b). Thus, when both Ca2+-inhibitable and Ca2+-non-inhibitable AC isoforms coexist, continuous oscillation is not usually observed for cellular cAMP and Ca2+.
The formulas and values used to generate this oscillation model are shown in online Supplemental methods and Supplemental Table 3.
Control systems in vivo are dynamic and complex, and it is very difficult to predict systems behavior on the basis of biochemical studies of individual molecules. However, use of systems biology tools, for example CellDesigner, makes it feasible to simulate complex biochemical networks flexibly [28–30]. In this study, this software enabled us to study Ca2+ and cAMP oscillations under different conditions, e.g., in the presence or absence of Ca2+-inhibitable and Ca2+-non-inhibitable AC isoforms and other regulatory molecules, in silico, without the need for experimental assays .
First, we confirmed that the original oscillation model of cAMP and Ca2+ developed by Cooper in 1995  can be effectively simulated by use of CellDesigner [23–25]. When the signaling pathway contained only Ca2+-inhibitable AC isoforms (AC5/6), we found that stable and spontaneous oscillations occurred.
Inclusion of PDE-mediated cAMP hydrolysis or PKA-mediated inhibition of AC5/6 into Cooper’s model did not seem to change the oscillation behavior, but altered the amplitude to a greater or lesser extent [13, 20–22, 26].
Inclusion of PKA-mediated inhibition of AC5/6 induced a decrease of the amplitude by approximately 13 % (Fig. 3), whereas inclusion of PDE-mediated cAMP hydrolysis increased the amplitude by approximately 1.6-fold (Supplemental Fig. 2), compared with those of Cooper’s original model (Fig. 1).
Phosphorylation of the L-type Ca2+-channel increases the Ca2+ concentration, and might form an ascending loop. In turn, cardiac AC isoforms (AC5/6) are inhibited by Ca2+ and this might form a descending loop. Importantly, AC5/6 are inhibited by PKA, which might reduce the amplitude of the oscillation loop [13, 22]. Conversely, PDE, in association with Ca2+-mediated inhibition, forms the descending phase of cAMP and Ca2+ oscillation in cardiac myocytes [31, 32]. These data, together with our current findings, indicate that PKA may have both positive and negative regulatory effects on the amplitude of the cAMP and Ca2+ oscillation loop, whereas PDE may have a positive regulatory effect on the amplitude [31, 32].
In contrast, when AC isoforms were not Ca2+-inhibitable, no oscillation occurred. Interestingly, when both Ca2+-non-inhibitable and Ca2+-inhibitable AC isoforms were included in the model, oscillation occurred, but decayed very rapidly. Inclusion of PDE-mediated cAMP hydrolysis did not change this behavior. Therefore, our simulations indicate that for stable and spontaneous oscillation, the presence of Ca2+-inhibitable AC isoforms and the absence of Ca2+-non-inhibitable AC isoforms are both required.
Further studies will be required to incorporate the effects of newly identified AC5/6-associated proteins, including Snapin, a SNAP25-binding protein, and PAM, a protein associated with Myc, on the cAMP and Ca2+ oscillations, because the findings of this study show that Ca2+-inhibitable AC isoforms (AC5/6) are essential for the oscillations of cAMP and Ca2+ [2, 13, 15].
Because the heart expresses seven AC isoforms [1, 7], including both Ca2+-inhibitable and Ca2+-non-inhibitable isoforms, continuous cAMP and Ca2+ oscillation in cardiac myocytes may not always occur, on the basis of the above findings (Fig. 4). Studies using AC5-deficient mice from our laboratory have shown that nearly half of the AC activity within the heart may be because of Ca2+-non-inhibitable AC isoforms [1, 2, 7, 9, 33, 34]. However, microenvironments in which only Ca2+-inhibitable AC5/6 are accumulated, such as lipid rafts or caveolae, are believed to exist [15, 16]. Indeed, it has been reported that Ca2+-sensitive AC isoforms (AC1, 5, 6, and 8) and their associated proteins, such as PKA, A-kinase anchoring proteins (AKAPs), anchored PDEs, non-anchored PDE, and transient receptor potential (TRP) 1/3, are present in lipid rafts in many cell types, including cardiac myocytes, whereas the Ca2+-insensitive AC 2, 4, and 7 are excluded from the rafts . Because intracellular cAMP and Ca2+ mediate a diverse array of cellular functions, oscillation of cAMP and Ca2+ concentration might be involved in receptor-mediated signal transduction, not only in excitable cells, for example cardiac myocytes, but also in non-excitable cells . Further, the occurrence or disturbance of cAMP and Ca2+ oscillations might contribute to the development of cardiac dysfunction or arrhythmia.
Okumura S, Suzuki S, Ishikawa Y (2009) New aspects for the treatment of cardiac diseases based on the diversity of functional controls on cardiac muscles: effects of targeted disruption of the type 5 adenylyl cyclase gene. J Pharmacol Sci 109:354–359
Okumura S, Kawabe J, Yatani A, Takagi G, Lee MC, Hong C, Liu J, Takagi I, Sadoshima J, Vatner DE, Vatner SF, Ishikawa Y (2003) Type 5 adenylyl cyclase disruption alters not only sympathetic but also parasympathetic and calcium-mediated cardiac regulation. Circ Res 93:364–371
Okumura S, Fujita T, Cai W, Jin M, Namekata I, Mototani Y, Jin H, Ohnuki Y, Tsuneoka Y, Bai Y, Suzuki S, Hidaka Y, Umemura M, Ichikawa Y, Yokoyama U, Sato M, Ishikawa F, Izumi-Nakaseko H, Adachi-Akahane S, Tanaka H, Ishikawa Y (2014) EPAC1-dependent phospholamban phosphorylation mediates the cardiac response to stresses. J Clin Invest 124:2785–2801
Ohnuki Y, Nishimura S, Sugiura S, Saeki Y (2008) Phosphorylation status of regulatory proteins and functional characteristics in myocardium of dilated cardiomyopathy of Syrian hamsters. J Physiol Sci 58:15–20
Bers DM (2002) Calcium and cardiac rhythms: physiological and pathophysiological. Circ Res 90:14–17
Ohnuki Y, Yamada T, Mototani Y, Umeki D, Shiozawa K, Fujita T, Saeki Y, Okumura S (2013) Effects of protein kinase a on the phosphorylation status and transverse stiffness of cardiac myofibrils. J Pharmacol Sci 123:279–283
Okumura S, Takagi G, Kawabe J, Yang G, Lee MC, Hong C, Liu J, Vatner DE, Sadoshima J, Vatner SF, Ishikawa Y (2003) Disruption of type 5 adenylyl cyclase gene preserves cardiac function against pressure overload. Proc Natl Acad Sci USA 100:9986–9990
Marks AR (2013) Calcium cycling proteins and heart failure: mechanisms and therapeutics. J Clin Invest 123:46–52
Okumura S, Vatner DE, Kurotani R, Bai Y, Gao S, Yuan Z, Iwatsubo K, Ulucan C, Kawabe J, Gjosh K, Vatner SF, Ishikawa Y (2007) Disruption of type 5 adenylyl cyclase enhances desensitization of cyclic adenosine monophosphate signal and increases Akt signal with chronic catecholamine stress. Circulation 116:1776–1783
Ohnuki Y, Umeki D, Mototani Y, Jin H, Cai W, Shiozawa K, Suita K, Saeki Y, Fujita T, Ishikawa Y, Okumura S (2014) Role of cyclic AMP sensor epac1 in masseter muscle hypertrophy and myosin heavy chain transition induced by β2-adrenoceptor stimulation. J Physiol 592:5461–5475
Ishikawa Y, Homcy CJ (1997) The adenylyl cyclases as integrators of transmembrane signal transduction. Circ Res 80:297–304
Cooper DM, Mons N, Karpen JW (1995) Adenylyl cyclases and the interaction between calcium and cAMP signalling. Nature 374:421–424
Hanoune J, Defer N (2001) Regulation and role of adenylyl cyclase isoforms. Annu Rev Pharmacol Toxicol 41:145–174
Katsushika S, Chen L, Kawabe J, Nilakantan R, Halnon NJ, Homcy C, Ishikawa Y (1992) Cloning and characterization of a sixth adenylyl cyclase isoform: types V and VI constitute a subgroup within the mammalian adenylyl cyclase family. Proc Natl Acad Sci USA 89:8774–8778
Willoughby D, Cooper DM (2007) Organization and Ca2+ regulation of adenylyl cyclases in cAMP microdomains. Physiol Rev 87:965–1010
Timofeyev V, Myers RE, Kim HJ, Woltz RL, Sirish P, Heiserman JP, Li N, Singapuri A, Tang T, Yarov-Yarovoy V, Yamoah EN, Hammond HK, Chiamvimonvat N (2013) Adenylyl cyclase subtype-specific compartmentalization: differential regulation of L-type Ca2+ current in ventricular myocytes. Circ Res 112:1567–1576
Ishikawa Y, Katsushika S, Chen L, Halnon NJ, Kawabe J, Homcy CJ (1992) Isolation and characterization of a novel cardiac adenylylcyclase cDNA. J Biol Chem 267:13553–13557
Espinasse I, Iourgenko V, Defer N, Samson F, Hanoune J, Mercadier JJ (1995) Type V, but not type VI, adenylyl cyclase mRNA accumulates in the rat heart during ontogenic development. Correlation with increased global adenylyl cyclase activity. J Mol Cell Cardiol 27:1789–1795
Tobise K, Ishikawa Y, Holmer SR, Im MJ, Newell JB, Yoshie H, Fujita M, Susannie EE, Homcy CJ (1994) Changes in type VI adenylyl cyclase isoform expression correlate with a decreased capacity for cAMP generation in the aging ventricle. Circ Res 74:596–603
Premont RT, Jacobowitz O, Iyengar R (1992) Lowered responsiveness of the catalyst of adenylyl cyclase to stimulation by GS in heterologous desensitization: a role for adenosine 3′,5′-monophosphate-dependent phosphorylation. Endocrinology 131:2774–2784
Chen Y, Harry A, Li J, Smit MJ, Bai X, Magnusson R, Pieroni JP, Weng G, Iyengar R (1997) Adenylyl cyclase 6 is selectively regulated by protein kinase A phosphorylation in a region involved in Galphas stimulation. Proc Natl Acad Sci USA 94:14100–14104
Iwami G, Kawabe J, Ebina T, Cannon PJ, Homcy CJ (1995) Regulation of adenylyl cyclase by protein kinase A. J Biol Chem 270:12481–12484
Kitano H, Funahashi A, Matsuoka Y, Oda K (2005) Using process diagrams for the graphical representation of biological networks. Nat Biotechnol 23:961–966
Funahashi A, Jouraku A, Matsuoka Y, Kitano H (2007) Integration of Cell Designer and SABIO-RK. In Silico Biol 7:S81–S90
Matsuoka Y, Funahashi A, Ghosh S, Kitano H (2014) Modeling and simulation using Cell Designer. Methods Mol Biol 1164:121–145
Sato N, Asai K, Okumura S, Takagi G, Shannon RP, Fujita-Yamaguchi Y, Ishikawa Y, Vatner SF, Vatner DE et al (1999) Mechanisms of desensitization to a PDE inhibitor (milrinone) in conscious dogs with heart failure. Am J Physiol 276:H1699–H1705
Tang T, Gao MH, Lai NC, Firth AL, Takahashi T et al (2008) Adenylyl cyclase type 6 deletion decreases left ventricular function via impaired calcium handling. Circulation 117:61–69
Autiero I, Costantini S, Colonna G (2009) Modeling of the bacterial mechanism of methicillin-resistance by a systems biology approach. PLoS One 4:e6226
Grieco L, Calzone L, Bernard-Pierrot I, Radvanyi F, Kahn-Perles B, Thieffy D (2013) Integrative modelling of the influence of MAPK network on cancer cell fate decision. PLoS Comput Biol 9:e1003286
Wu G, Zhu L, Dent JE, Nardini C (2010) A comprehensive molecular interaction map for rheumatoid arthritis. PLoS One 5:e10137
Mika D, Richter W, Westenbroek RE, Catterall WA, Conti M (2014) PDE4B mediates local feedback regulation of β1-adrenergic cAMP signaling in a sarcolemmal compartment of cardiac myocytes. J Cell Sci 127:1033–1042
Masaki N, Fujimoto K, Honda-Kitahara M, Hada E, Sawai S (2013) Robustness of self-organizing chemoattractant field arising from precise pulse induction of its breakdown enzyme: a single-cell level analysis of PDE expression in Dictyostelium. Biophys J 104:1191–1202
Okumura S, Tsunematsu T, Bai Y, Jiao Q, Ono S, Suzuki S, Kurotani R, Sato M, Minamisawa S, Umemura S, Ishikawa Y (2008) Type 5 adenylyl cyclase plays a major role in stabilizing heart rate in response to microgravity induced by parabolic flight. J Appl Physiol 1985 105:173–179
Bai Y, Tsunematsu T, Jiao Q, Ohnuki Y, Mototani Y, Shiozawa K, Jin M, Cai W, Jin HL, Fujita T, Ichikawa Y, Suita K, Kurotani R, Yokoyama U, Sato M, Iwatsubo K, Isikawa Y (2012) Pharmacological stimulation of type 5 adenylyl cyclase stabilizes heart rate under both microgravity and hypergravity induced by parabolic flight. J Pharmacol Sci 119:381–389
Dolmetsch RE, Xu K, Lewis RS (1998) Calcium oscillations increase the efficiency and specificity of gene expression. Nature 392:933–936
This study was supported in part by grants from the Ministry of Health, Labor, and Welfare (Dr Ishikawa), a Grant-in-Aid for Scientific Research on Innovative Areas (22136009), and Grants from the Japanese Ministry of Education, Culture, Sports, Science, and Technology (Drs Ishikawa, Okumura, Fujita, Sato, Yokoyama, Mototani), the National Cerebral and Cardiovascular Center (Dr Ishikawa), the Takeda Science Foundation (Dr Okumura), the Yokohama Foundation for Advancement of Medical Science (Dr Okumura), the Yokohama Academic Foundation (Dr Ohnuki), the Research Foundation for Community Medicine (Dr Okumura), and the Suzuken Memorial Foundation (Dr Okumura).
Conflict of interest
The authors declare no potential conflicts of interest.
Electronic supplementary material
About this article
Cite this article
Kamide, T., Okumura, S., Ghosh, S. et al. Oscillation of cAMP and Ca2+ in cardiac myocytes: a systems biology approach. J Physiol Sci 65, 195–200 (2015). https://doi.org/10.1007/s12576-014-0354-3
- Adenlyl cyclase
- Computational analysis
- Cyclic AMP