Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 03 February 2021
Sec. Vascular Physiology
This article is part of the Research Topic Diabetes, Hypertension and Cardiovascular Diseases View all 14 articles

Downregulation of Soluble Guanylate Cyclase and Protein Kinase G With Upregulated ROCK2 in the Pulmonary Artery Leads to Thromboxane A2 Sensitization in Monocrotaline-Induced Pulmonary Hypertensive Rats

\r\nSuhan ChoSuhan Cho1Hyun NamgoongHyun Namgoong1Hae Jin Kim,Hae Jin Kim1,2Rany VornRany Vorn3Hae Young YooHae Young Yoo3Sung Joon Kim,*Sung Joon Kim1,2*
  • 1Department of Physiology, College of Medicine, Seoul National University, Seoul, South Korea
  • 2Ischemic/Hypoxic Disease Institute, College of Medicine, Seoul National University, Seoul, South Korea
  • 3Department of Nursing, Chung-Ang University, Seoul, South Korea

Thromboxane A2 (TXA2) promotes various physiological responses including pulmonary artery (PA) contraction, and pathophysiological implications have been suggested in cardiovascular diseases including pulmonary hypertension. Here, we investigated the role of TXA2 receptor (TP)-mediated signaling in the pathophysiology of pulmonary arterial hypertension (PAH). The sensitivity of PA to the contractile agonist could be set by relaxing signals such as the nitric oxide (NO), soluble guanylate cyclase (sGC), and cGMP-dependent kinase (PKG) pathways. Changes in the TP agonist (U46619)-induced PA contraction and its modulation by NO/cGMP signaling were analyzed in a monocrotaline-induced PAH rat model (PAH-MCT). In the myograph study, PA from PAH-MCT showed higher responsiveness to U46619, that is decreased EC50. Immunoblot analysis revealed a lower expression of eNOS, sGC, and PKG, while there was a higher expression of RhoA-dependent kinase 2 (ROCK2) in the PA from PAH-MCT than in the control. In PAH-MCT, the higher sensitivity to U46619 was reversed by 8-Br-cGMP, a membrane-permeable cGMP analog, but not by the NO donor, sodium nitroprusside (SNP 30 μM). In contrast, in the control PA, inhibition of sGC by its inhibitor (1H− [1,2,4] oxadiazolo [4,3−a] quinoxalin-1-one (ODQ), 10 μM) lowered the threshold of U46619-induced contraction. In the presence of ODQ, SNP treatment had no effect whereas the addition of 8-Br-cGMP lowered the sensitivity to U46619. The inhibition of ROCK by Y-27632 attenuated the sensitivity to U46619 in both control and PAH-MCT. The study suggests that the attenuation of NO/cGMP signaling and the upregulation of ROCK2 increase the sensitivity to TXA2 in the PAH animal, which might have pathophysiological implications in patients with PAH.

Introduction

Thromboxane A2 (TXA2), a metabolic product of arachidonic acid, is involved in various physiological activities such as platelet aggregation, airway narrowing, and contraction of various types of arteries, including the pulmonary artery (PA). TXA2 is synthesized by platelets and parenchymal cells of the intestine, kidney, and lung. In the pulmonary circulation system, TXA2 appears to act as a physiological modulator of blood flow distribution. For instance, pretreatment with TXA2 greatly facilitates hypoxic pulmonary vasoconstriction in vitro (Park et al., 2012; Yoo et al., 2012) and in vivo (Kylhammar and Rådegran, 2012).

Due to the short half-life (<30 s) of TXA2, a stable synthetic analog U46619 is widely used to investigate its physiological and pathophysiological roles (Hamberg et al., 1975; Nakahata, 2008). TXA2 and U46619 act through the thromboxane-prostanoid (TP) receptor which is a G-protein coupled receptor. The plasma TXA2 levels correlated with the prevalence of cardiovascular diseases, including hypertension (Francois et al., 2004), atherosclerosis (Mehta et al., 1988), ischemic heart disease (Serneri et al., 1981), and stroke (Stier et al., 1988). In addition, increased TXA2 activity was observed in pulmonary hypertension (Zamora et al., 1993).

Against the contractile signals in vascular smooth muscle cells (VSMCs) via phospholipase C (PLC)-coupled G-protein coupled receptors (GPCR) and voltage-operated L-type Ca2+ channels (VOCCL), the tone of arteries is regulated by endothelium-derived relaxing factors (EDRFs) and endothelium-derived hyperpolarizing (EDH) mechanisms. Among EDRFs, the release of NO potently regulates vascular tone through the soluble guanylate cyclase (sGC)/cGMP/protein kinase G (PKG) pathway in VSMCs (Buvinic and Huidobro-Toro, 2001). The production of NO via an endothelial type of NO synthase (eNOS) occurs primarily through an increase in cytosolic Ca2+ concentration ([Ca2+]c) in the endothelium and Akt-dependent phosphorylation of eNOS (Zhao et al., 2015).

The stimulation of the TP receptor in vascular smooth muscle increases [Ca2+]c via both stored Ca2+ release via PLC/IP3 signaling and VOCCL-dependent Ca2+ influx mechanisms including the activation of non-selective cation channels (Dorn and Becker, 1993; Nakahata, 2008; Yoo et al., 2012). In addition to the direct contractile signals to VSMCs, the stimulation of TP receptors in endothelial cells inhibited K+ channels (e.g., SKCa) through an unknown mechanism, and thus reduced the EDH signals (Ellinsworth et al., 2014). Gα12/13-dependent RhoA-mediated activation of kinase (ROCK) is activated by TP receptor-mediated signals that produce vascular contraction (Kozasa et al., 1998). The activation of ROCK increases the Ca2+-sensitivity of smooth muscle via inhibition of myosin light chain phosphatase (MLCP) (Butler et al., 2013). The U46619-induced arterial contraction was significantly attenuated by Y-27632, an inhibitor of ROCK (Fu et al., 1998). In addition to contractile signaling, our recent study suggested that TP receptor-mediated stimulation of eNOS in the PA smooth muscle layer, counterbalances the potent vasoconstrictive effect of TXA2 via the NO/cGMP-dependent signaling pathway (Kim et al., 2016).

NO dysregulation appears to be important in vascular tone regulation and PA remodeling in pulmonary arterial hypertension (PAH). Decreased NO availability is a common phenomenon in patients with PAH with impairment in the biosynthesis of NO (Klinger et al., 2013). In PAH, PA remodeling occurs along with changes in the molecular properties of smooth muscle cells (Schermuly et al., 2011). In cultured adult rat pulmonary arterial smooth muscle cells (PASMCs), downregulation of sGC and PKG expression occurred in parallel to phenotypic changes ranging from contractile to synthetic type, which might explain the PA remodeling in PAH (Lincoln et al., 2006). Interestingly, in a Sugen5416/hypoxia-induced PAH rat model, inhibition of NOS by NG-Nitro-L-arginine-methyl ester (L-NAME), revealed biphasic changes in NO availability, and decreased and increased in the early and late phase PAH, respectively. The latter change might counterbalance the excessive contractile tone of the chronic hypoxia-induced PAH arteries (Tanaka et al., 2017).

As a classical animal model of PAH, monocrotaline-injected rats (MCT-PAH) show typical medial thickening of PA with right ventricular hypertrophy. In PAH-MCT rats, a molecular biological investigation of PA tissue revealed increased expression of ROCK2 and decreased sGC (Lee et al., 2016). These changes may decrease the efficiency of the NO/cGMP pathway to enhance responsiveness to contractile agonists such as TXA2. Taken together with the vasoactive signals of TXA2 including eNOS activation in rat PA, we hypothesized that the putative downregulation of sGC in PAH-MCT might increase the responsiveness of PA to TXA2. However, the functional changes, including the sensitivity to the relevant pharmacological agents, were not rigorously investigated in the remodeled PA. In the present study, the concentration-dependent contractile responses of rat PA to U46619 were analyzed using a Mulvany-type isometric tension recording system. In addition, an immunoblot assay was performed to confirm and further elucidate the changes in the levels of TP receptor, eNOS, sGC, and ROCK in PA from PAH-MCT.

Materials and Methods

Animal Model

All experimental procedures were conducted with the approval of the Institutional Animal Care and Use Committee of Seoul National University (approval number: SNU-190408-3). Monocrotaline (Sigma, St. Louis, MO, United States) was dissolved in 2 mL of 1 M HCl and adjusted to pH 7.4 using 2 M NaOH solution. This aqueous solution was diluted to 17 mL with distilled water. Seven-week-old male Sprague-Dawley rats were randomly assigned and treated with a single intraperitoneal injection of monocrotaline (60 mg/kg) to induce the PAH-MCT model or an appropriate amount of saline. After 21 days, both the PAH-MCT rats and the age-matched control were sacrificed for further analysis, after measuring their body weight (b.w.).

Arterial Tissue Preparation

Lung tissues were excised and stored in ice-cold normal Tyrode’s (NT) solution (140 mM of NaCl, 5.4 mM of KCl, 0.33 mM of NaH2PO4, 10 mM of HEPES, 10 mM of glucose, 1.8 mM of CaCl2, and 1 mM of MgCl2, pH 7.4 was adjusted with NaOH). The pulmonary artery bed was collected by removing the lung and bronchial tissues in NT solution. Excised pulmonary arteries were cleaned from perivascular adipose and other tissues under a microscope. Prepared arterial tissues were used for further analysis.

Isometric Tension Recordings

Isometric tension was measured using a dual-wire multi-channel myograph system (620 M; DMT, Aarhus, Denmark). Excised arteries were cut into 2 mm arterial ring segments and mounted with 25 μm tungsten wire on an NT solution-filled organ chamber for tension recording. For stabilization, mounted arteries were rested in physiological salt solution (PSS; 118 mM of NaCl, 4 mM of KCl, 24 mM of NaHCO3, 1 mM of MgSO4, 0.44 mM of NaH2PO4, 5.6 mM of glucose, and 1.8 mM of CaCl2) for at least 15 min with a gas mixture (21% O2, 5% CO2, N2 balance) after a basal tone was applied. The whole experiment was maintained at 37°C. For the experiment, 80 mM of KCl-PSS was used to induce contraction (80K) for the evaluation of vessel integrity, and the dose-dependent response of U46619 was evaluated with the application of differential dosages from 1 to 200 nM.

Histology

Pulmonary arterial ring segments were washed with phosphate-buffered saline (PBS) and fixed in 4% paraformaldehyde overnight. For histological analysis, paraffin-embedded tissues were cut and stained with hematoxylin and eosin (H&E) or Masson’s trichrome (MT). Digital images of the stained sections were obtained at 200 × magnification using the Aperio ImageScope 12.3 software.

Western Blotting

Whole pulmonary arteries from the left lung were excised and homogenized with RIPA buffer (Millipore, United States) and protease/phosphatase inhibitor cocktail (Roche Diagnostics, Germany) for 1 h at 4°C. The samples were centrifuged at 13,000 ×g for 30 min at 4°C and the supernatant was collected. Protein concentration was determined using the PierceTM BCA Protein Assay Kit (Thermo Fisher Scientific, United States). The loading samples were prepared with Laemmli sample buffer, resolved by 8% SDS-PAGE, and transferred to polyvinylidene difluoride membranes in 25 mM of Tris, 192 mM of glycine, 0.01% SDS, and 20% methanol. Membranes were blocked in 1 × TBS containing 1% Tween-20 and 5% skim milk (blocking solution) for 1 h at room temperature with gentle rocking and incubated overnight at 4°C with relevant primary antibodies for detecting protein expression. Mouse anti-eNOS monoclonal antibody 1:1000 (Abcam, ab76198,United Kingdom), rabbit anti-sGC-α polyclonal antibody 1:1000 (Cayman, 160895, United States), rabbit anti-sGC-β1 polyclonal antibody 1:1000 (Cayman, 160897, United States), rabbit anti-PKG monoclonal antibody 1:1000 (Cell signaling, 3248, United States), rabbit anti-ROCK1 monoclonal antibody 1:1000 (Abcam, ab45171, United Kingdom), rabbit anti-ROCK2 polyclonal antibody 1:1000 (Abcam, ab71598, United Kingdom), rabbit anti-TXA2R polyclonal antibody 1:1000 (Abcam, ab233288, United Kingdom), and mouse anti-β-actin monoclonal antibody 1:10000 (Sigma, A1978, United States) were used. The signals were determined using ECL Plus Western blotting detection reagents (Amersham Biosciences, United Kingdom) and detected images were obtained by the AmershamTM Imager 600 (Amersham Biosciences, United Kingdom). The intensity of each band was measured using the ImageJ analysis software.

cGMP Concentration Measurement

cGMP concentration in the pulmonary arterial tissues was measured by a commercially available cGMP enzyme immunoassay (EIA) kit (Cayman). The excised pulmonary arteries were pulverized and homogenized with 5% trichloroacetic acid (TCA). The samples were centrifuged at 1,500 ×g for 30 min at 4°C and the supernatant was collected. The supernatants were analyzed according to the manufacturer’s instructions.

Statistical Analysis

All data are expressed as mean ± SEM and the number of tested arteries is indicated as n. For all comparisons, arteries were obtained from at least three rats per protocol. Statistical differences between PAH-MCT and normal rats were analyzed using a two-sample T test. ANOVA was used to investigate more than two groups, and Tukey correction was performed for post hoc testing for significant differences. Two-way ANOVA was used to show the statistical differences in the concentration-dependent contractile response of U46619 between groups. Statistical differences were analyzed using SPSS Version 23 (IBM, SPSS, Chicago) and differences were considered when the p-value was less than 0.05.

Results

Increased Contractile Sensitivity to U46619 in the PA of PAH-MCT

To characterize the PAH-MCT model, the b.w. changes, mass ratio of the right ventricle, and histological features of PA were analyzed. The b.w. on the 21st day after the MCT injection was lower in PAH-MCT than in CON (Figure 1A). The mass ratio of the right ventricles over the sum of the left ventricle and septum (RV/LV + S) was significantly higher in PAH-MCT than in CON (Figure 1B) (0.23 ± 0.006 vs. 0.57 ± 0.096; p = 0.0012). PA from PAH-MCT displayed a thicker medial layer than the PA from CON (Figure 1C, upper panel). MT staining revealed increased collagen deposition in the medial layer of PA from PAH-MCT (Figure 1C, lower panel).

FIGURE 1
www.frontiersin.org

Figure 1. Increased contractile responsiveness of PA to TXA2 in PAH-MCT. Characteristics of PAH-MCT model with lower b.w. (A), higher relative mass of the right ventricle (R/L + S mass ratio, (B), and medial thickening of PA revealed by the H&E stained light microscopy [(C), upper panel]. MT staining showed MCT-induced collagen deposition in PAH-MCT [(C), lower panel, magnification of 200×. Scale bar: 100 μm]. Isometric tension recording of PA from CON and PAH-MCT rats. Contractile force of Δ80K and U46619 were represented with g (D). The contractile responses to the stepwise increase in the concentration of U46619 (E). The contractile responsiveness was analyzed with concentration-dependent contractile response curve [C-R curve, (F)]. The bar graphs represent the means ± SEMs. Statistical difference between groups are indicated as ***P < 0.001; bar graph with two sample t-test, and C-R curve with two-way ANOVA. (G) Representative three pairs of immunoblotting analyses comparing the expression of TP receptors in PAs between CON and PAH-MCT.

In each vessel, the contractile response to 80 mM of KCl (80K-contraction) was initially measured. The PA of PAH-MCT and CON did not show a significant difference at 80K-contraction (Figure 1D). To analyze the sensitivity to U46619, the concentration-dependent increase in contractile tone was measured with the isometric tension recording system (Figure 1E). The maximum contraction induced by 200 nM of U46619 was higher than that induced by 80K-contraction in CON. The U46619-induced maximum contraction showed a tendency of smaller amplitude in PAH-MCT, and statistical significance was not found owing to the cases with sizable amplitude (Figure 1D). When normalized to the maximum U46619-contraction, the concentration-dependent contractile response curve (C-R curve) showed significantly higher sensitivity to U46619 in PAH-MCT than in CON. Notably, the contractile response at relatively low concentrations (5–20 nM) was more significant in PAH-MCT than in CON (Figure 1F). Despite the higher sensitivity to U46619 in PAH-MCT, the expression of TXA2 receptor (TP) was not different between CON and PAH-MCT (Figure 1G).

Decreased Expression of NO/cGMP-Regulating Proteins and Increased Expression of ROCK2 in the PA of PAH-MCT

The sGC and PKG are known to mediate the relaxing signal from eNOS and NO to downstream proteins such as myosin light chain phosphatase (MLCP) and ROCK (Figure 2A). The immunoblot study showed decreased eNOS expression in PA of PAH-MCT (Figure 2B). sGC is a heterodimer composed of one alpha (sGC-α) and one heme-binding beta subunit (sGC-β). Immunoblot analysis of sGC showed a tendency of lower expression in PAH-MCT whereas statistical significance was found in the level of sGC-β1 (Figures 2B,C). The expression of PKG was also lower in PAH-MCT than in CON (p = 0.030). In contrast, the expression of ROCK2 was increased in PAH-MCT (p = 0.043), while the level of ROCK1 was not changed (Figure 2C). cGMP level in the PA tissues was lower in PAH-MCT than in CON (Figure 2D).

FIGURE 2
www.frontiersin.org

Figure 2. Altered expression of NO/cGMP-regulating enzymes and ROCK in PAH-MCT. (A) Schematic drawing of NO/cGMP pathway molecules involved in the regulation of vascular smooth muscle contractility. (B) Representative western blot images of NO-sGC-cGMP pathway molecules with β-actin as a standard control. (C) Results of band intensity relative to β -actin (n = 4, respectively). (D) cGMP levels of PA tissues measured with EIA (n = 4). *P < 0.05 vs. PAH-MCT; Mann-Whitney U test. (E,F) Effect of SNP (30 μM on the C-R curves of CON (E) and PAH-MCT (F). ***P < 0.001; two-way ANOVA.

To investigate the contribution of relaxing signals distal to NO in the U46619-induced contraction, we analyzed the C-R curve after pretreatment with the NO donor sodium nitroprusside (SNP, 30 μM, 10 min). In the CON PA, the C-R curve was shifted to the right by SNP, implying lowered sensitivity to U46619 (Figure 2E). In contrast, the pretreatment with SNP did not affect the C-R curve in PAH-MCT (Figure 2F).

Pharmacological Regulation of eNOS, sGC, and ROCK in PAH-MCT and CON

To mimic the downregulation of eNOS and sGC proteins in the PAH-MCT, the PA of CON was pretreated with NOS inhibitor (L-NAME, 100 μM) or sGC inhibitor (ODQ, 10 μM). Both pharmacological inhibitors increased the sensitivity of PA to U46619 in CON, and the effect of L-NAME was relatively more prominent (Figures 3A,B). In the presence of L-NAME, the exogenous NO donor (SNP, 30 μM) shifted the C-R curve to the right, while SNP had no effect on the C-R curve in the presence of ODQ (Figures 3C,D). We compensated for sGC inhibition by using 8-bromoguanosine 3′, 5′-cyclic monophosphate (8-Br-cGMP), a membrane permeable analog of cGMP. The additional pretreatment with 5 μM of 8-Br-cGMP shifted the C-R curve to the right in the presence of ODQ (Figures 3E,F). Since the inhibition of ROCK2 was a plausible activity of PKG in the PA myocytes, we evaluated the effect of the ROCK inhibitor Y-27632. The additional pretreatment with 10 μM of Y-27632 also shifted the C-R curve to the right in the presence of ODQ (Figures 3G,H).

FIGURE 3
www.frontiersin.org

Figure 3. Mimicking the functional changes in PAH-MCT by the inhibitors of eNOS and sGC in the PA of CON. (A,B) Effect of preincubation with L-NAME (100 μM) and ODQ (10 μM) on the C-R curves of CON. ***P < 0.001 L-NAME vs CON and ##P < 0.01 ODQ vs CON; two-way ANOVA. (C,D) Effect of SNP (30 μM on the sensitized C-R curves by L-NAME (C) or by ODQ (D). (E,F) Effect of 8-Br-cGMP (30 μM) on the sensitized contraction and the C-R curve shifted by ODQ. (G,H) Effect of Y-27632 (10 μM) on the sensitized contraction and the C-R curve shifted by ODQ. Statistical difference between groups are indicated with **P < 0.01 or ***P < 0.001; two-way ANOVA.

Next, the effects of the above pharmacological agents were tested in the PA of PAH-MCT. In contrast to the PA of CON, neither the inhibition of eNOS (L-NAME) nor sGC (ODQ) had an effect on the C-R curve of U46619 (Figures 4A,B). Despite the insignificant effect of L-NAME and ODQ, pretreatment with 8-Br-cGMP shifted the C-R curve to the right in the PA of PAH-MCT (Figures 4C,D). In addition, pretreatment with Y-27632 shifted the C-R curve to the right (Figures 4E,F).

FIGURE 4
www.frontiersin.org

Figure 4. Effects of the inhibitors of eNOS, sGC, and ROCK on the C-R curves of U46619 in the PA of PAH-MCT. (A,B) No significant effect of pretreatment with L-NAME (100 μM) and ODQ (10 μM) on the C-R curves of PAH-MCT. (C,D) Rightward shift of the C-R curve by 8-Br-cGMP (30 μM). (E,F) Rightward shift of the C-R curve by Y-27632 (10 μM). Statistical difference between groups are indicated with ***P < 0.001; two-way ANOVA.

Since the C-R curves shown above are normalized values (% contraction induced by 200 nM of U46619), the pretreatment with vasodilators such as 8-Br-cGMP and Y-27632 affected the maximum tone in the tested vessels. Since the initial 80K contraction was measured in the absence of the pharmacological agents, we compared the normalized maximum contraction (% of 80K contraction) in each tested condition (Figure 5). To investigate this, a single dose of U46619 (200 nM) was applied in the absence or presence of 8-Br-cGMP and ODQ. Pretreatment with 8-Br-cGMP did not affect the maximum contraction induced by U46619. Furthermore, the combined treatment with ODQ (8-Br-cGMP + ODQ) did not have a significant effect on maximum contraction (Figure 5A).

FIGURE 5
www.frontiersin.org

Figure 5. Comparative analyzes of the signaling pathways determining the maximum tone of PA induced by 200 nM of U46619. (A) The pretreatment with 8-Br-cGMP alone or with ODQ did not affect the maximum contraction normalized to the 80K-contraction in each PA from CON and PAH-MCT. (B) Different effects of the pretreatment with Y27632 on the maximum tone of CON and PAH-MCT. Note that the inhibition by Y27632 was prevented by the co-treatment with ODQ in CON while not altered in PAH-MCT. (C) Effects of GF109203X (10 μM) and nifedipine (1 μM) on the maximum contraction induced by 200 nM of U46619. The bars represent the means ± SEMs. Statistical difference between groups are indicated as ***P < 0.001. Comparison of the effects of 8-Br-cGMP and Y27632.

In contrast to the effect of 8-Br-cGMP, the direct inhibition of ROCK (Y-27632) largely suppressed the maximum contraction by U46619. However, when pretreated with both ODQ and Y-27632, the maximum contraction was restored to the control level in the PA from CON (Figure 5B, left panel). Interestingly, in the PA from PAH-MCT, neither pretreatment with Y-27632 alone nor the co-treatment with ODQ affected the maximum contraction by U46619 (Figure 5B, right panel).

The resistance to Y-27632 in the maximum contraction of PAH-MCT might be due to a hidden compensatory upregulation of contractile signaling pathways other than ROCK, such as voltage-operated L-type Ca2+ channels (VOCCL) or protein kinase C (PKC). To test this possibility, we examined the effect of the VOCCL inhibitor (nifedipine, 1 μM) and an inhibitor of PKC (GF109203X, 10 μM). Both inhibitors decreased the maximum contraction by around 70% in PAH-MCT as well as in CON (Figure 5C).

Discussion

The present study showed that the PA segments from rat with PAH-MCT are contracted by nanomolar ranges of the TXA2 analog, which could be accounted for by the downregulation of NO-related signaling enzymes (eNOS, sGC, and PKG) and by the upregulation of ROCK2 (Figure 6). Although the membrane localization of TP receptors was not analyzed, the immunoblot analysis did not show changes in TP receptor expression (Figure 1G). Consistent with the schematic summary of the changes in PAH-MCT (Figure 6), the pharmacological inhibitors of eNOS and sGC shifted the C-R curves of CON PA in the rightward direction, but not in the PA from PAH-MCT (Figures 3, 4). In PAH-MCT, the inhibition of ROCK effectively restored the C-R curve similar to that in CON (Figure 4F). Although the PKG level in PA was downregulated in PAH-MCT, the application of the PKG activator 8-Br-cGMP reversed the increased responsiveness to U46619 in PAH-MCT (Figure 4D). Taken together with the results of the immunoblotting assay, the analysis of the C-R curves appears to be an appropriate approach to reveal the functional implication of the changes in the intracellular signals in PASMCs.

FIGURE 6
www.frontiersin.org

Figure 6. Schematic diagrams of the altered signaling pathways suggested for the increased TXA2 sensitivity in PAH-MCT. (A) Under the basal or stimulated signaling through NO-sGC-cGMP-PKG cascade, PKG activates MLCP but inhibits ROCK to prevent excessive contraction by TXA2. (B) Under the pathological condition of PAH, decreased signaling via NO-sGC-cGMP-PKG would increase the MLC phosphorylation, resulting in increased responsiveness to TXA2.

TP receptor signaling in the arterial contraction is mediated through Gα12/13 to RhoA/ROCK or through Gαq/11 to PLCβand InsP3/diacylglycerol (DAG)/PKC signaling pathways. It is generally known that activated ROCK induces MLC2 phosphorylation indirectly, which is mediated by the inhibitory phosphorylation of Thr residues (T696 and T853) of MYPT that also prevent the stimulatory phosphorylation of neighboring Ser residues (S695 and S852) (Butler et al., 2013). MYPT1 is a member of the myosin phosphatase targeting protein (MYPT) family consisting of five genes, MYPT1, MYPT2, MBS85, MYPT3, and TIMAP, which function as targeting and regulatory subunits to confer substrate specificity and subcellular localization on the catalytic subunit of type 1δ protein serine/threonine phosphatase (PP1cδ) (Grassie et al., 2011). Taken together, we propose that the upregulation of ROCK along with downregulated sGC/PKG in PAH-MCT are responsible for the sensitized contractile response to TXA2.

It was notable that the maximum contraction induced by 200 nM of U46619 was not affected by SNP or 8-Br-cGMP. In contrast, the inhibition of VOCCL and PKC significantly reduced the maximum contraction (Figures 5A,C). These results suggested that NO/cGMP-dependent signals have a modulatory role in MLC phosphorylation rather than a direct inhibition of Ca2+ signals and contractile mechanisms elicited by TP stimulation. The inhibition of ROCK showed differential effects on maximum contraction; Y-27632 significantly suppressed the maximum contraction in CON, but not in PAH-MCT. Interestingly, even in the PA of CON, Y-27632 did not decrease the maximum contraction in the presence of ODQ (Figure 5B). The results suggested a dominant relaxing signal from cGMP-PKG rather than a contractile signal via ROCK that converges on MYPT.

The phosphorylation states of Ser695/Ser852 and Thr696/Thr853 have mutually antagonistic effects on MYPT activity (Qiao et al., 2014). PKG phosphorylates and activates the Ser residues and induces translocation to MLC2, the dephosphorylation of which promotes vasorelaxation. In contrast, ROCK and other pro-contractile signals increase the phosphorylation of Thr residues in MYPT which induces membrane translocation and lowers the chances of interaction with MLC2 (Butler et al., 2013). In this scheme, the deficiency in cGMP-PKG signaling (e.g., ODQ treatment) would not allow for the control of MLC2 by MYPT, and the inhibitory effect of Y-27632 by preventing the interaction between MYPT and MLC2 could not be demonstrated. Since the downregulated eNOS/sGC/PKG in the PA of PAH-MCT would be equivalent to pharmacological inhibition, the effect of Y-27632 was not observed, especially in the conditions of maximum contraction.

Taken together, our present results could be interpreted based on the schematic model of the antagonistic influence on the MLCP activity between the sGC/cGMP/PKG and RhoA/ROCK pathways (Figure 6). The downregulation of sGC and PKG possibly minimized the role of MLCP, which is also the target of ROCK. Therefore, the pharmacological inhibition of ROCK might not affect the maximum level of contraction, while still shifting the C-R curve. Since the relative contribution from the other contractile signaling pathways was not altered in the PA from PAH-MCT (Figure 5C), we cautiously emphasized the critical counteractive regulatory effects of the eNOS-sGC-PKG and RhoA-ROCK pathways on MLCP in PA contraction by TP receptor stimulation. The altered expression of the above signaling molecules in the PASMCs would result in the dysregulation of vasodilation activity (Figure 6B).

The immunoblot assay revealed a significant decrease of eNOS expression in the PA of PAH-MCT. Downregulation of eNOS vascular endothelium is one of the major pathophysiological findings in cardiovascular diseases (Bauersachs et al., 1998; Chou et al., 1998; Kobayashi et al., 1999; Li et al., 2002; Lee et al., 2016). It is generally known that endothelial eNOS is a major source of NO which diffuses to neighboring tissues to activate signaling cascades, including sGC. In rat PA, we have previously reported that the medial layer, most likely the smooth muscle cells also express eNOS. Although the level of eNOS expression in the PA medial layer is minute compared to the endothelium, TP stimulation could induce the muscular eNOS phosphorylation to partly contribute to the relaxing signals intrinsic to the myocytes (Kim et al., 2016). However, in the present study, the C-R curves were obtained in PA with an intact endothelium since the process of endothelium removal by mechanical abrasion generally induces unstable responses during U46619-induced contraction experiments. Considering the relatively minute expression of eNOS in the medial layer of PA, the decreased eNOS expression in the PAH-MCT would mostly reflect the downregulation in the endothelium (Figure 2).

Along with eNOS, the downregulation of the sGC subunit, especially sGC-β1, was a prominent and consistent finding in PAH-MCT. The critical role of sGC in the regulation of vascular contractility could be demonstrated by the effect of ODQ. In this respect, the abolished influence of ODQ treatment on the C-R curve of PA from PAH-MCT (Figure 4B) suggested that the downregulation of sGC was functionally crucial and could be underestimated when simply interpreted from the decreased expression level (Figure 2C). The downregulation of sGC and the decreased cGMP (Figure 2D) observed in the present study could be consistent with the rationale for the clinical application of PDE5 inhibitors for the treatment of PAH. In addition, a direct activator of sGC, riociguat, has been recently introduced for clinical studies (Montani et al., 2014). The concomitant partial decrease in PKG expression should be considered to interpret the functional changes under pathological conditions. Decreased PKG expression has also been reported in primary aortic vascular smooth muscle cells treated with inflammatory cytokines (Browner et al., 2004; Choi et al., 2018).

ROCK2 expression was increased in PA from PAH-MCT and the C-R curve was restored with ROCK inhibition. The pathophysiological role of ROCK-dependent signaling has also been found in systemic arteries. In the mesenteric artery of spontaneous hypertensive rats, H2O2 induces c-Src-dependent TXA2 release and provokes vascular contractile responses through multiple signaling pathways, including ROCK (García-Redondo et al., 2015). Consistent with the pathophysiological implication of the TXA2 pathway, both genetic and pharmacological suppression of TXA2/TP signaling confers microvascular protection against oxidative injury in mesenteric arteries (Chiang, 2019).

In addition to the changes in PA contractility, increased synthesis of TXA2 and local availability might also underlie the higher tension under pathological conditions. Experimentally, it was reported that a direct injection of U46619 induced a severe prolonged increase in systolic blood pressure, and death of the tested mice, which was not observed in the TP knockout mice (Sparks et al., 2013). Clinically, there is evidence that increased TXA2 synthesis possibly contributes to increased vascular resistance in children with pulmonary hypertension and pregnancy-induced hypertension (Fitzgerald et al., 1990; Adatia et al., 1993). PAH patients showed higher TXA2 synthesis with decreased prostacyclin (PGI2) (Christman et al., 1992; Liel et al., 1993; Kazimierczyk and Kamiński, 2018) and platelet aggregation (Nakonechnicov et al., 1996; Varol et al., 2011), and increased thrombosis is one of the criteria for PAH (Farber and Loscalzo, 2004). Despite the repeatedly observed pathophysiological role of TXA2 in PAH, the pharmacological antagonist of the TP receptor has not been clinically investigated yet. Instead, the administration of PGI2 analog is one of the options for PAH treatment (Montani et al., 2014).

Conclusion

During the pathological condition of PAH, the combined increase in TXA2 availability and contractile sensitivity to TXA2 may worsen the burden of pulmonary circulation. Along with the functional changes, PA wall thickening dramatically increased peripheral resistance and right ventricular afterload in PAH. The markedly lowered threshold concentration of U46619 suggested vulnerability to collapse or uneven regional perfusion in the PAH lung. Since the C-R curves and their changes by pharmacological agents are consistent with the molecular biological findings than the simple comparison of maximum contraction levels, our present study emphasizes the implication of physiological parameters obtained from myography studies.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

The animal study was reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) of Seoul National University.

Author Contributions

SC and SK: conception and design of study. SC, HN, and RV: conducting experiments, acquisition of data, analysis, and interpretation of data. SC, HY, and HK: drafting the manuscript. SK: revising the manuscript critically for important intellectual content. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT of the Republic of Korea (grant nos. NRF-2018R1A5A2025964 and NRF-2018R1D1A1B07048998). In addition, this work was partly supported by the Education and Research Encouragement Fund of Seoul National University Hospital (2020).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Adatia, I., Barrow, S., Stratton, P., Miall-Allen, V., Ritter, J., and Haworth, S. (1993). Thromboxane A2 and prostacyclin biosynthesis in children and adolescents with pulmonary vascular disease. Circulation 88, 2117–2122. doi: 10.1161/01.cir.88.5.2117

CrossRef Full Text | Google Scholar

Bauersachs, J., Bouloumie, A., Mülsch, A., Wiemer, G., Fleming, I., and Busse, R. (1998). Vasodilator dysfunction in aged spontaneously hypertensive rats: changes in NO synthase III and soluble guanylyl cyclase expression, and in superoxide anion production. Cardiovasc. Res. 37, 772–779. doi: 10.1016/s0008-6363(97)00250-2

CrossRef Full Text | Google Scholar

Browner, N. C., Sellak, H., and Lincoln, T. M. (2004). Downregulation of cGMP-dependent protein kinase expression by inflammatory cytokines in vascular smooth muscle cells. Am. J. Physiol. Cell Physiol. 287, C88–C96.

Google Scholar

Butler, T., Paul, J., Europe-Finner, N., Smith, R., and Chan, E.-C. (2013). Role of serine-threonine phosphoprotein phosphatases in smooth muscle contractility. Am. J. Physiol. Cell Physiol. 304, C485–C504.

Google Scholar

Buvinic, S., and Huidobro-Toro, J. P. (2001). Basal tonic release of nitric oxide coupled to cGMP production regulates the vascular reactivity of the mesenteric bed. Eur. J. Pharmacol. 424, 221–227. doi: 10.1016/s0014-2999(01)01165-7

CrossRef Full Text | Google Scholar

Chiang, C.-Y. (2019). Genetic Depletion of Thromboxane A2/Thromboxane-Prostanoid Receptor Signalling Prevents Microvascular Dysfunction in Ischaemia/Reperfusion Injury. Taiwan: National Taiwan Normal University.

Google Scholar

Choi, S., Park, M., Kim, J., Park, W., Kim, S., Lee, D.-K., et al. (2018). TNF-α elicits phenotypic and functional alterations of vascular smooth muscle cells by miR-155-5p–dependent down-regulation of cGMP-dependent kinase 1. J. Biol. Chem. 293, 14812–14822. doi: 10.1074/jbc.ra118.004220

PubMed Abstract | CrossRef Full Text | Google Scholar

Chou, T.-C., Yen, M.-H., Li, C.-Y., and Ding, Y.-A. (1998). Alterations of nitric oxide synthase expression with aging and hypertension in rats. Hypertension 31, 643–648. doi: 10.1161/01.hyp.31.2.643

CrossRef Full Text | Google Scholar

Christman, B. W., McPherson, C. D., Newman, J. H., King, G. A., Bernard, G. R., Groves, B. M., et al. (1992). An imbalance between the excretion of thromboxane and prostacyclin metabolites in pulmonary hypertension. N. Engl. J. Med. 327, 70–75. doi: 10.1056/nejm199207093270202

PubMed Abstract | CrossRef Full Text | Google Scholar

Dorn, G.d., and Becker, M. W. (1993). Thromboxane A2 stimulated signal transduction in vascular smooth muscle. J. Pharmacol. Exp. Ther. 265, 447–456.

Google Scholar

Ellinsworth, D. C., Shukla, N., Fleming, I., and Jeremy, J. Y. (2014). Interactions between thromboxane A2, thromboxane/prostaglandin (TP) receptors, and endothelium-derived hyperpolarization. Cardiovasc. Res. 102, 9–16. doi: 10.1093/cvr/cvu015

PubMed Abstract | CrossRef Full Text | Google Scholar

Farber, H. W., and Loscalzo, J. (2004). Pulmonary arterial hypertension. N. Engl. J. Med. 351, 1655–1665.

Google Scholar

Fitzgerald, D., Rocki, W., Murray, R., Mayo, G., and FitzGerald, G. (1990). Thromboxane A2 synthesis in pregnancy-induced hypertension. Lancet 335, 751–754. doi: 10.1016/0140-6736(90)90869-7

CrossRef Full Text | Google Scholar

Francois, H., Athirakul, K., Mao, L., Rockman, H., and Coffman, T. M. (2004). Role for thromboxane receptors in angiotensin-II–induced hypertension. Hypertension 43, 364–369. doi: 10.1161/01.hyp.0000112225.27560.24

CrossRef Full Text | Google Scholar

Fu, X., Gong, M. C., Jia, T., Somlyo, A. V., and Somlyo, A. P. (1998). The effects of the Rho-kinase inhibitor Y-27632 on arachidonic acid-, GTPγ S-, and phorbol ester-induced Ca2 +-sensitization of smooth muscle. FEBS Lett. 440, 183–187. doi: 10.1016/s0014-5793(98)01455-0

CrossRef Full Text | Google Scholar

García-Redondo, A. B., Briones, A. M., Martínez-Revelles, S., Palao, T., Vila, L., Alonso, M. J., et al. (2015). c-Src, ERK1/2 and Rho kinase mediate hydrogen peroxide-induced vascular contraction in hypertension: role of TXA2, NAD(P)H oxidase and mitochondria. J. Hypertens. 33, 77–87. doi: 10.1097/hjh.0000000000000383

PubMed Abstract | CrossRef Full Text | Google Scholar

Grassie, M. E., Moffat, L. D., Walsh, M. P., and MacDonald, J. A. (2011). The myosin phosphatase targeting protein (MYPT) family: a regulated mechanism for achieving substrate specificity of the catalytic subunit of protein phosphatase type 1δ. Arch. Biochem. Biophys. 510, 147–159. doi: 10.1016/j.abb.2011.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamberg, M., Svensson, J., and Samuelsson, B. (1975). Thromboxanes: a new group of biologically active compounds derived from prostaglandin endoperoxides. Proc. Natl. Acad. Sci. U.S.A. 72, 2994–2998. doi: 10.1073/pnas.72.8.2994

PubMed Abstract | CrossRef Full Text | Google Scholar

Kazimierczyk, R., and Kamiński, K. (2018). The role of platelets in the development and progression of pulmonary arterial hypertension. Adv. Med. Sci. 63, 312–316. doi: 10.1016/j.advms.2018.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, H. J., Yoo, H. Y., Jang, J. H., Lin, H. Y., Seo, E. Y., Zhang, Y. H., et al. (2016). Wall stretch and thromboxane A2 activate NO synthase (eNOS) in pulmonary arterial smooth muscle cells via H2O2 and Akt-dependent phosphorylation. Pflügers Arch. Eur. J. Physiol. 468, 705–716. doi: 10.1007/s00424-015-1778-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Klinger, J. R., Abman, S. H., and Gladwin, M. T. (2013). Nitric oxide deficiency and endothelial dysfunction in pulmonary arterial hypertension. Am. J. Respirat. Crit. Care Med. 188, 639–646. doi: 10.1164/rccm.201304-0686pp

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobayashi, N., Kobayashi, K., Hara, K., Higashi, T., Yanaka, H., Yagi, S., et al. (1999). Benidipine stimulates nitric oxide synthase and improves coronary circulation in hypertensive rats. Am. J. Hypertens. 12, 483–491. doi: 10.1016/s0895-7061(98)00260-x

CrossRef Full Text | Google Scholar

Kozasa, T., Jiang, X., Hart, M. J., Sternweis, P. M., Singer, W. D., Gilman, A. G., et al. (1998). p115 RhoGEF, a GTPase activating protein for Gα12 and Gα13. Science 280, 2109–2111. doi: 10.1126/science.280.5372.2109

PubMed Abstract | CrossRef Full Text | Google Scholar

Kylhammar, D., and Rådegran, G. (2012). Cyclooxygenase-2 inhibition and thromboxane A2 receptor antagonism attenuate hypoxic pulmonary vasoconstriction in a porcine model. Acta Physiol. 205, 507–519. doi: 10.1111/j.1748-1716.2012.02437.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, M.-Y., Tsai, K.-B., Hsu, J.-H., Shin, S.-J., Wu, J.-R., and Yeh, J.-L. (2016). Liraglutide prevents and reverses monocrotaline-induced pulmonary arterial hypertension by suppressing ET-1 and enhancing eNOS/sGC/PKG pathways. Sci. Rep. 6:31788.

Google Scholar

Li, H., Wallerath, T., Münzel, T., and Förstermann, U. (2002). Regulation of endothelial-type NO synthase expression in pathophysiology and in response to drugs. Nitric Oxide 7, 149–164. doi: 10.1016/s1089-8603(02)00111-8

CrossRef Full Text | Google Scholar

Liel, N., Nathan, I., Yermiyahu, T., Zolotov, Z., Lieberman, J., Dvilansky, A., et al. (1993). Increased platelet thromboxane A2/prostaglandin H2 receptors in patients with pregnancy induced hypertension. Thrombosis Res. 70, 205–210. doi: 10.1016/0049-3848(93)90126-9

CrossRef Full Text | Google Scholar

Lincoln, T. M., Wu, X., Sellak, H., Dey, N., and Choi, C.-S. (2006). Regulation of vascular smooth muscle cell phenotype by cyclic GMP and cyclic GMP-dependent protein kinase. Front. Biosci. 11:356–367. doi: 10.2741/1803

PubMed Abstract | CrossRef Full Text | Google Scholar

Mehta, J. L., Lawson, D., Mehta, P., and Saldeen, T. (1988). Increased prostacyclin and thromboxane A2 biosynthesis in atherosclerosis. Proc. Natl. Acad. Sci. U.S.A. 85, 4511–4515. doi: 10.1073/pnas.85.12.4511

PubMed Abstract | CrossRef Full Text | Google Scholar

Montani, D., Chaumais, M.-C., Guignabert, C., Günther, S., Girerd, B., Jaïs, X., et al. (2014). Targeted therapies in pulmonary arterial hypertension. Pharmacol. Ther. 141, 172–191.

Google Scholar

Nakahata, N. (2008). Thromboxane A2: physiology/pathophysiology, cellular signal transduction and pharmacology. Pharmacol. Ther. 118, 18–35. doi: 10.1016/j.pharmthera.2008.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakonechnicov, S., Gabbasov, Z., Chazova, I., Popov, E., and Yu, B. (1996). Platelet aggregation in patients with primary pulmonary hypertension. Blood Coagul. Fibrinol. Int. J. Haemostasis Thrombosis 7, 225–227. doi: 10.1097/00001721-199603000-00029

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, S. J., Yoo, H. Y., Kim, H. J., Kim, J. K., Zhang, Y.-H., and Kim, S. J. (2012). Requirement of pretone by thromboxane a2 for hypoxic pulmonary vasoconstriction in precision-cut lung slices of rat. Kor. J. Physiol. Pharmacol. 16, 59–64. doi: 10.4196/kjpp.2012.16.1.59

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiao, Y. N., He, W. Q., Chen, C. P., Zhang, C. H., Zhao, W., Wang, P., et al. (2014). Myosin phosphatase target subunit 1 (MYPT1) regulates the contraction and relaxation of vascular smooth muscle and maintains blood pressure. J. Biol. Chem. 289, 22512–22523. doi: 10.1074/jbc.M113.525444

PubMed Abstract | CrossRef Full Text | Google Scholar

Schermuly, R. T., Ghofrani, H. A., Wilkins, M. R., and Grimminger, F. (2011). Mechanisms of disease: pulmonary arterial hypertension. Nat. Rev. Cardiol. 8:443.

Google Scholar

Serneri, G. G. N., Gensini, G. F., Abbate, R., Mugnaini, C., Favilla, S., Brunelli, C., et al. (1981). Increased fibrinopeptide A formation and thromboxane A2 production in patients with ischemic heart disease: relationships to coronary pathoanatomy, risk factors, and clinical manifestations. Am. Heart J. 101, 185–194. doi: 10.1016/0002-8703(81)90665-7

CrossRef Full Text | Google Scholar

Sparks, M. A., Makhanova, N. A., Griffiths, R. C., Snouwaert, J. N., Koller, B. H., and Coffman, T. M. (2013). Thromboxane receptors in smooth muscle promote hypertension, vascular remodeling, and sudden death. Hypertension 61, 166–173. doi: 10.1161/hypertensionaha.112.193250

PubMed Abstract | CrossRef Full Text | Google Scholar

Stier, C. Jr., Benter, I., and Levine, S. (1988). Thromboxane A2 in severe hypertension and stroke in stroke-prone spontaneously hypertensive rats. Stroke 19, 1145–1150. doi: 10.1161/01.str.19.9.1145

CrossRef Full Text | Google Scholar

Tanaka, M., Abe, K., Oka, M., Saku, K., Yoshida, K., Ishikawa, T., et al. (2017). Inhibition of nitric oxide synthase unmasks vigorous vasoconstriction in established pulmonary arterial hypertension. Physiol. Rep. 5:e13537. doi: 10.14814/phy2.13537

PubMed Abstract | CrossRef Full Text | Google Scholar

Varol, E., Uysal, B. A., and Ozaydin, M. (2011). Platelet indices in patients with pulmonary arterial hypertension. Clin. Appl. Thrombosis Hemostasis 17, E171–E174.

Google Scholar

Yoo, H. Y., Park, S. J., Seo, E.-Y., Park, K. S., Han, J.-A., Kim, K. S., et al. (2012). Role of thromboxane A2-activated nonselective cation channels in hypoxic pulmonary vasoconstriction of rat. Am. J. Physiol. Cell Physiol. 302, C307–C317.

Google Scholar

Zamora, C. A., Baron, D. A., and Heffner, J. E. (1993). Thromboxane contributes to pulmonary hypertension in ischemia-reperfusion lung injury. J. Appl. Physiol. 74, 224–229. doi: 10.1152/jappl.1993.74.1.224

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Vanhoutte, P. M., and Leung, S. W. (2015). Vascular nitric oxide: Beyond eNOS. J. Pharmacol. Sci. 129, 83–94. doi: 10.1016/j.jphs.2015.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: pulmonary hypertension, pulmonary artery smooth muscle, guanylate cyclase, protein kinase G, thromboxane A2 (TXA2), monocrotaline

Citation: Cho S, Namgoong H, Kim HJ, Vorn R, Yoo HY and Kim SJ (2021) Downregulation of Soluble Guanylate Cyclase and Protein Kinase G With Upregulated ROCK2 in the Pulmonary Artery Leads to Thromboxane A2 Sensitization in Monocrotaline-Induced Pulmonary Hypertensive Rats. Front. Physiol. 12:624967. doi: 10.3389/fphys.2021.624967

Received: 01 November 2020; Accepted: 13 January 2021;
Published: 03 February 2021.

Edited by:

Soo-Kyoung Choi, Yonsei University Health System, South Korea

Reviewed by:

María Galán, Sant Pau Institute for Biomedical Research, Spain
Young-Eun Cho, The University of Iowa, United States

Copyright © 2021 Cho, Namgoong, Kim, Vorn, Yoo and Kim. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Sung Joon Kim, physiolksj@gmail.com; sjoonkim@snu.ac.kr

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.