A correction has been applied to this article in:
Corrigendum: Impact of a TAK-1 inhibitor as a single or as an add-on therapy to riociguat on the metabolic reprograming and pulmonary hypertension in the SUGEN5416/hypoxia rat model
Daniel Morales-Cano1,2,3,4,5†
Jose Luis Izquierdo-García2,6,7†
Bianca Barreira1,2,3
Sergio Esquivel-Ruiz1,2,3Maria Callejo1,2,3
Rachele Pandolfi1,2,3
Palmira Villa-Valverde2,8Ignacio Rodríguez2,6
Angel Cogolludo1,2,3Jesus Ruiz-Cabello2,6,9‡
Francisco Perez-Vizcaino1,2,3‡
Laura Moreno1,2,3*‡- 1Department of Pharmacology and Toxicology, School of Medicine, Universidad Complutense de Madrid, Madrid, Spain
- 2Ciber Enfermedades Respiratorias (Ciberes), Madrid, Spain
- 3Instituto de Investigación Sanitaria Gregorio Marañón (IISGM), Madrid, Spain
- 4Centro Nacional de Investigaciones Cardiovasculares (CNIC), Madrid, Spain
- 5Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
- 6Department of Chemistry in Pharmaceutical Sciences, School of Pharmacy, Universidad Complutense de Madrid, Madrid, Spain
- 7Instituto Pluridisciplinar, Universidad Complutense de Madrid, Madrid, Spain
- 8ICTS Bioimagen Complutense, Universidad Complutense de Madrid, Madrid, Spain
- 9Center for Cooperative Research in Biomaterials (CIC biomaGUNE), Basque Research and Technology Alliance (BRTA), Donostia San Sebastián, Spain
Background: Despite increasing evidence suggesting that pulmonary arterial hypertension (PAH) is a complex disease involving vasoconstriction, thrombosis, inflammation, metabolic dysregulation and vascular proliferation, all the drugs approved for PAH mainly act as vasodilating agents. Since excessive TGF-β signaling is believed to be a critical factor in pulmonary vascular remodeling, we hypothesized that blocking TGFβ-activated kinase 1 (TAK-1), alone or in combination with a vasodilator therapy (i.e., riociguat) could achieve a greater therapeutic benefit.
Methods: PAH was induced in male Wistar rats by a single injection of the VEGF receptor antagonist SU5416 (20 mg/kg) followed by exposure to hypoxia (10%O2) for 21 days. Two weeks after SU5416 administration, vehicle, riociguat (3 mg/kg/day), the TAK-1 inhibitor 5Z-7-oxozeaenol (OXO, 3 mg/kg/day), or both drugs combined were administered for 7 days. Metabolic profiling of right ventricle (RV), lung tissues and PA smooth muscle cells (PASMCs) extracts were performed by magnetic resonance spectroscopy, and the differences between groups analyzed by multivariate statistical methods.
Results: In vitro, riociguat induced potent vasodilator effects in isolated pulmonary arteries (PA) with negligible antiproliferative effects and metabolic changes in PASMCs. In contrast, 5Z-7-oxozeaenol effectively inhibited the proliferation of PASMCs characterized by a broad metabolic reprogramming but had no acute vasodilator effects. In vivo, treatment with riociguat partially reduced the increase in pulmonary arterial pressure (PAP), RV hypertrophy (RVH), and pulmonary vascular remodeling, attenuated the dysregulation of inosine, glucose, creatine and phosphocholine (PC) in RV and fully abolished the increase in lung IL-1β expression. By contrast, 5Z-7-oxozeaenol significantly reduced pulmonary vascular remodeling and attenuated the metabolic shifts of glucose and PC in RV but had no effects on PAP or RVH. Importantly, combined therapy had an additive effect on pulmonary vascular remodeling and induced a significant metabolic effect over taurine, amino acids, glycolysis, and TCA cycle metabolism via glycine-serine-threonine metabolism. However, it did not improve the effects induced by riociguat alone on pulmonary pressure or RV remodeling. None of the treatments attenuated pulmonary endothelial dysfunction and hyperresponsiveness to serotonin in isolated PA.
Conclusion: Our results suggest that inhibition of TAK-1 induces antiproliferative effects and its addition to short-term vasodilator therapy enhances the beneficial effects on pulmonary vascular remodeling and RV metabolic reprogramming in experimental PAH.
Introduction
Pulmonary arterial hypertension (PAH) is a chronic and progressive disease characterized by a sustained elevation of pulmonary arterial pressure (PAP) and associated with high morbidity and mortality (Simonneau et al., 2013; Nicolls and Voelkel, 2017; Humbert et al., 2019). As disease severity progresses, right ventricular (RV) function worsens due to sustained pressure overload, leading to RV failure. Indeed, RV failure remains the primary cause of death in patients with PAH (Vonk Noordegraaf et al., 2017).
The pathophysiology of PAH is complex, involving vasoconstriction, thrombosis, inflammation, metabolic dysregulation and vascular proliferation. While much progress has been made in the treatment of PAH, currently approved drugs for PAH act as primary vasodilating agents targeting endothelin, nitric oxide, and prostacyclin pathways (Rabinovitch, 2012; Nicolls and Voelkel, 2017; Prins et al., 2019a; Humbert et al., 2019; Koudstaal et al., 2020). Although these therapies have led to significant improvements in overall medial survival and quality of life of patients, none of them are curative, and PAH remains a fatal disease (Montani et al., 2014; Scott et al., 2021). Furthermore, despite RV dysfunction determines the prognosis of PAH patients, therapies targeting the RV are still needed (Vonk Noordegraaf et al., 2017). Therefore, novel therapies must go beyond vasodilatation and start targeting maladaptive vascular and RV remodeling if we aim to improve the outcomes in PAH patients and progress toward a cure (Prins et al., 2019a; Prins et al., 2019b; Scott et al., 2021).
Extensive research during the last decades has identified new cellular and molecular pathways that have yet to be exploited therapeutically (Prins et al., 2019b; Andre et al., 2021; Scott et al., 2021; Sharmin et al., 2021). For example, dysregulation of the TGF-β superfamily pathways has been proven critical in the pathogenesis of PAH (Rol et al., 2018; Andre et al., 2021; Sharmin et al., 2021). Thus, in PAH, deficient bone morphogenetic protein (BMP) type II receptor (BMPR-II) impairs the anti-proliferative SMAD 1/5/8 signaling pathways whilst the overactive TGFβ pathway leads to activation of the SMAD2/3 pathway signaling and results in excessive proliferation, attenuation of apoptosis, production of inflammatory cytokines, mitochondrial dysfunction (Marsboom et al., 2012; Freund-Michel et al., 2014) and a metabolic reprogramming mainly characterized by a shift in glucose metabolism from oxidative phosphorylation to glycolysis (also known as the Warburg effect (Rol et al., 2018; Hernandez-Saavedra et al., 2020; Scott et al., 2021; Sharmin et al., 2021). These pieces of evidence suggest that therapeutic interventions eliciting anti-TGF-β and pro-BMP effects have the potential to be disease modifying. Indeed, sotatercept, a ligand trap for activin (another member of the TGF-β superfamily) which is able to restore the balance between BMP/TGF-β is currently undergoing a Phase III clinical trial following the encouraging results from the Phase II PULSAR trial (Humbert et al., 2021).
In this regard, TGF-β-associated kinase 1 (TAK1) seems crucial in the cross-regulation between the BMP and the TGFβ arms. Thus, with fully functional BMPR-II-mediated signaling, TAK1 is retained by the BMPR-II receptor complex. By contrast, through BMPR2 dysregulation, the TAK1 intracellular pool becomes available for the TGFβ receptor complex, leading to abnormal proliferation and apoptosis (Nasim et al., 2012). In line with this, TAK1 has been shown to mediate myocardial and vascular remodeling in systemic hypertension (Song et al., 2014; Garfield et al., 2019; Meijles et al., 2020).
Based on these evidences, we hypothesized that sustained inhibition of TAK1 with the low molecular weight inhibitor 5Z-7-oxozeaenol in combination with a vasodilator therapy such as the soluble guanylate cyclase (sGC) stimulator riociguat could ameliorate or even reverse pulmonary vascular remodeling and RV dysfunction in PAH by their ability to limit inflammation and metabolic disruption and to mitigate vascular and cardiac remodeling. Pharmacological treatment was evaluated in pulmonary arterial smooth muscle cells (PASMCs) and a relevant PAH animal model to validate our hypothesis. Finally, cell extract and intact lung and heart tissues were analyzed by Nuclear Magnetic Resonance spectroscopy (NMR) to characterize the metabolic state after treatment.
Materials and methods
All experimental procedures involving animals were approved by the Institutional and Regional Ethical Committees of Madrid (Spain; PROEX 251/15). Male Wistar rats (8–10 weeks of age; body weight ≈220 g) obtained from Envigo (Barcelona, Spain) were kept in the Animal Core Facility of University Complutense of Madrid under standard conditions of temperature (24°C ± 1°C) and 12 h dark/light cycle with ad libitum access to food and water. Animal studies are reported in compliance with the ARRIVE guidelines (Kilkenny et al., 2010).
Experimental procedure for the rat SU5416/Hypoxia model
Adult male Wistar rats were injected subcutaneously with a single dose of the VEGFR type 2 inhibitor SUGEN 5416 (SU5416; 20 mg/kg) or the equivalent volume of the vehicle CMC (0.5% [w/v] carboxymethyl cellulose sodium, 0.9% [w/v] sodium chloride, 0.4% [v/v] Tween 80, 0.9% [v/v] benzyl alcohol in deionized water) before being exposed to normobaric hypoxia (10% oxygen) for 3 weeks (SuHyp model) (Moral-Sanz et al., 2012; Mondejar-Parreno et al., 2019; Callejo et al., 2020). 2 weeks after the administration of SU5416, rats were further randomly assigned to receive the treatment with vehicle, the TAK-1 inhibitor 5Z-7-oxozeaenol (3 mg/kg/day), riociguat (3 mg/kg/day) or both drugs combined for an additional period of 7 days. Drugs were administered by daily intraperitoneal injection following the study flowchart shown in Supplementary Figure S1A.
Hemodynamic measurements
Twenty-four hours after the last administration of each treatment, the rats were anesthetized (80 mg/kg ketamine and 8 mg/kg xylazine i. p.), intubated through a tracheostomy and ventilated with room air (tidal volume 9 ml/kg, 60 breaths/min, and a positive end-expiratory pressure of 2 cm H2O, Nemi Scientific Inc., Medway, United States of America (Cogolludo et al., 2009; Moral-Sanz et al., 2012; Morales-Cano et al., 2014). Using an open-chest approach, a pressure transducer was placed, via a catheter, into the RV and then advanced into the pulmonary artery (PA). Systolic, end-diastolic and mean RV pressure (RVSP, RVEDP and RVMP) and systolic, diastolic and mean PA pressures (SPAP, DPAP and MPAP) were measured over three consecutive stable recordings. It should be noted that the pressure and heart rate values obtained in anaesthetized open-chest rats may represent an underestimation of the values in closed-chest conscious animals (Cogolludo et al., 2009; Moral-Sanz et al., 2012; Morales-Cano et al., 2014).
Assessment of RV hypertrophy
At the end of the recordings, animals were sacrificed by overdose of anesthesia. Hearts were excised and the RV and the left ventricle plus septum (LV + S) were carefully dissected and weighed before being snap frozen for subsequent analysis. The ratios RV to body weight and RV to LV + S (Fulton Index) were calculated to assess the RV hypertrophy.
Lung histology
The left lung lobes were immediately collected for vascular reactivity or snap frozen for metabolomic profiling. The right lung was inflated in situ with formol saline through the right bronchus and embedded in paraffin. Lung sections were stained with eosin and hematoxylin, and elastin was visualized by its green autofluorescence. Small arteries (50–200 μm outer diameter) were analyzed in a blinded fashion and categorized as muscular, partly muscular, or non-muscular as previously described (Moral-Sanz et al., 2012; Morales-Cano et al., 2014). The medial wall thickness was calculated as external elastic lamina diameter minus the internal lamina diameter using ImageJ (Vs. 1.41, NIH, Bethesda, MD, United States of America).
Vascular reactivity
Intrapulmonary arteries were carefully dissected free of surrounding tissue from the inferior right lung lobe and cut into rings (1.8–2 mm length). PA rings were mounted on a wire myograph in Krebs physiological solution continuously bubbled with 21% O2, 5% CO2, and 74% N2 (Normoxia) and stretched to a transmural pressure equivalent to 30 mmHg, as previously described (Cogolludo et al., 2009; Morales-Cano et al., 2014). To confirm smooth muscle viability, arteries were first stimulated by raising the potassium concentration of the buffer to 80 mM and then allowed to recover. After washout, cumulative concentration responses curves to serotonin (5-HT, 30 nM–30 µM) were performed. Thereafter, concentration-response curves to acetylcholine (ACh; 1 nM−10 µM) and sodium nitroprusside (SNP, 10 pM–10 µM) were performed by cumulative addition to assess the endothelium-dependent and endothelium-independent vasodilatation, respectively.
In another set of experiments, rat PA rings derived from healthy animals were submaximally precontracted with the thromboxane A2 mimetic U46619 (0.1 µM), and cumulative concentration-response curves (10 nM–10 µM) to 5Z-7-oxozeaenol or riociguat in the absence or presence of 1 µM 5Z-7-oxozeaenol, were determined.
Western blotting analysis
Lung homogenates were applied to a sulfate-polyacrylamide gel, and proteins were transferred to polyvinylidene difluoride membranes and incubated with primary antibodies directed against the following proteins: TAK-1, ERK1/2, p38 and JNK or the phosphorylated forms of ERK1/2, and JNK (all: Cell Signaling Technology, MA), p-TAK-1 (Abcam, United Kingdom), p-p38 (Santa Cruz, United States of America) or β-actin (Sigma Aldrich, Spain). Membranes were then incubated with secondary peroxidase-conjugated antibodies. Antibody binding was detected by an ECL system (Amersham Pharmacia Biotech, Amersham, United Kingdom), and images were acquired using Odyssey Fc System (Li-COR Biosciences, United States) and densitometry analysis performed using Quantity One software. Results were normalized by the relative expression of smooth muscle β-actin.
Analysis of IL-1β production
Levels of IL-1β in whole lung homogenates were quantified by specific Rat DuoSet ELISA Development Systems (R&D System, United States) (Pandolfi et al., 2017).
Cell proliferation and metabolic extraction
Rat PA smooth muscle cells (PASMCs) primary cell cultures were generated as previously described from SuHyp rats (Pandolfi et al., 2017). PASMC, seeded at 30.000 cell/ml (approximately 50% confluence) in 96-well plates, were growth arrested for 24 h in 0.1% fetal calf serum (FCS) DMEM medium and then exposed to 5Z-7-oxozeaenol, riociguat (1nM–1 µM) or vehicle for 48 h in 10% FCS or fresh 0.1% FCS. Cell proliferation was assessed using the BrdU incorporation assay. In another set of experiments, cells were washed with saline to remove residual medium and fixed with 10 ml each of ice-cold methanol, chloroform, and deionized water. Following phase separation, the aqueous phase was lyophilized and resuspended in 400 µl deuterium oxide (Cambridge Isotope Laboratories) for metabolic profiling.
NMR data acquisition
Intact lung and RV tissue samples were examined by HR-MAS NMR using a Bruker AMX500 spectrometer (11.7 T). Samples were placed into a 50-μl zirconium oxide rotor using a rinsed cylindrical insert with 15 µl of a 0.1 mM solution of TSP in deuterium water (D2O) and spun at 4,000 Hz to remove the effects of spinning sidebands from the acquired spectra. Cell extracts were analyzed using a 500 MHz Bruker Ultrashield Plus spectrometer.
Shimming and NMR preparation times were reduced to a minimum, while the sample was chilled to 4°C to minimize metabolic changes. Several two-dimensional homonuclear and heteronuclear experiments, such as standard gradient-enhanced correlation spectroscopy (COSY), 1H–1H total correlated spectroscopy (TOCSY) and gradient-selected heteronuclear single quantum correlation (HSQC), were performed to carry out metabolomic assignments. A control 1H NMR spectrum was measured between consecutive two-dimensional (2D) spectra. No gross degradation was noted in the signals of multiple spectra acquired under the same conditions.
Standard solvent-suppressed spectra were grouped into 32,000 data points, averaged over 256 acquisitions. The data acquisition lasted a total of 15 min using a sequence based on the first increment of the nuclear Overhauser effect spectroscopy (NOESY) pulse sequence to suppress the effects of water. Sample acquisition was performed using a spectral width of 6,024 Hz prior to Fourier transformation, and the free induction decay (FID) signals were multiplied by an exponential weight function corresponding to a line broadening of 0.3 Hz. The spectra were referenced to the TSP singlet at a chemical shift of 0 ppm.
The NMR spectra were processed as described previously (Izquierdo-Garcia et al., 2011). Briefly, 1H-NMR spectra were binned to equal-length integral segments (δ = 0.01 ppm), and normalized to the total sum of the spectral regions.
Metabolomic analysis
Principal components analysis (PCA) and Partial Least Square Discriminant Analysis (PLS-DA) (Hotelling, 1933) were applied to 0.01 binned spectra to extract the most discriminative spectral subset from the total pool of metabolites and to remove outliers. The data obtained from this analysis were centered and scaled. The potential biomarkers selected from the PLS-DA loading matrix were confirmed by Hoteling’s T2 test (Hotelling, 1933). Multivariate statistical analysis was performed with ‘stats’ R package version 4.3.0 (Team, 2016). In addition, the resonances highlighted by the PCA were deconvoluted for metabolic quantification using the Global Spectral Deconvolution algorithm of MestRenova v. 8.1 (Mestrelab Research S.L., Santiago de Compostela, Spain) (Xia et al., 2009; Chong et al., 2018). The resonances were identified according to Chenomx NMR suit metabolic database (software version 9.0; Chenomx Inc. Edmonton, Canada)), the Human Metabolome Database (Wishart et al., 2012), and characteristic cross-peaks from 2D spectra to help in unequivocal assignation of these metabolites. Supplementary Table S1 display the complete list of metabolites assigned to detected NMR signals. For metabolic quantification, statistical significance was determined using a Bonferroni-corrected Student’s t-test (Vinaixa et al., 2012) assuming unequal variance, and variations were considered statistically significant when the adjusted p-value was less than 0.05.
Metabolic pathways analyses
Metabolic pathway analysis was performed using the Pathway Analysis module (Xia and MetPA, 2010) of Metaboanalyst v.5.0 (Chong et al., 2018), which combines results from robust pathway enrichment analysis (Kankainen et al., 2011) with pathway topology analysis (Aittokallio and Schwikowski, 2006) to help researchers identify the most relevant pathways involved in the conditions under study. Briefly, pathway enrichment analysis examines whether metabolites in predefined pathways are at the top or bottom of a ranked list. In contrast, pathway topology analysis applies graph theory to measure the importance of an experimentally identified metabolite in a predefined metabolic pathway. KEGG metabolic pathways were used as the backend knowledgebase. The selected pathway enrichment analysis method was GlobalAncova (Hummel et al., 2008), node importance measure for topological analysis was out-degree centrality. Centrality is a standard metric used in graph theory to estimate the relative importance of individual nodes to the overall network (Liu et al., 2019). Out-degree is the number of outgoing links or successor nodes (Xia et al., 2009; Chong et al., 2018).
Statistical analysis
Results are expressed as mean ± SEM. Technical replicates were averaged to provide a single data point before further analysis. Statistical analysis was performed using GraphPad Prism 5 as detailed in each figure legend. Individual cumulative concentration-response curves were fitted to a logistic equation and maximal responses (Emax) and the negative log of IC50 (concentration of the drug that induced half-maximal response) were calculated. For multiple comparisons, normally distributed data were analyzed by one-way ANOVA followed by Dunnet’s post hoc test (Dunnett, 1955). For metabolic quantification, statistical significance was determined using a Bonferroni-corrected Student’s t-test (Vinaixa et al., 2012), assuming unequal variance. Adjusted p-value less than 0.05 was considered statistically significant.
Results
In vitro characterization of the vasodilator and antiproliferative effects of riociguat and 5Z-7-oxozeaenol in pulmonary vascular tissues.
We first analyzed the vasodilator and antiproliferative effects in vitro of the two drugs and the potential interactions between them. Riociguat induced a concentration-dependent relaxation in U46619-stimulated rat PA (Figure 1A; Emax = 92 ± 8%; pIC50 = 6.8 ± 0.4) but had negligible effects on the proliferation of rat PASMCs induced by 10% of FCS (Figure 1B). In contrast, the TAK-1 inhibitor 5Z-7-oxozeaenol only exerted a minor relaxant response at the highest concentration tested (10 μM; Emax = 21 ± 5%; Figure 1A) whereas it inhibited PASMCs proliferation in a concentration-dependent manner (Emax = 46 ± 3%; pIC30 = 6.9 ± 0.2; Figure 1B). Notably, neither 5Z-7-oxozeaenol (1 µM) affected the relaxant responses induced by riociguat (Figure 1A), nor activation of sGC modulated the antiproliferative effects exerted by 5Z-7-oxozeaenol (Figure 1B).
FIGURE 1
FIGURE 1. Vasodilator and antiproliferative effects of the TAK-1 inhibitor 5Z-7-oxozeaenol and the soluble guanylate cyclase activator riociguat in healthy rat pulmonary arteries. (A) Concentration-dependent relaxant effects of 5Z-7-oxozeaenol (OXO) and riociguat (RIO) in the absence or presence of OXO (1 µM) in PA stimulated with 0.1 µM of U46619. (B) Concentration-dependent effects of RIO and OXO in the absence or presence of RIO (1 µM) in the proliferation induced by 10% of FCS in PASMCs. n = 6-7 vascular rings (A) or experiments run in triplicate (B) for each group *p < 0,05 versus control (two-way ANOVA followed by Dunnett’s post hoc test).
These in vitro studies suggested that riociguat acts mainly as a vasodilator agent whereas 5Z-7-oxozeaenol is an effective antiproliferative agent, targeting two central pathological features in PH, but there were no interactions between them. Therefore, we next investigated the potential therapeutic effects of 5Z-7-oxozeaenol as a single or as an add-on therapy to riociguat in the SuHyp rat in vivo model of PAH.
In vivo study of riociguat and 5Z-7-oxozeaenol in the SuHyp rat model of PH
PH was induced in rats by a single injection of SU5416 (20 mg/kg, sc) and chronic exposure to normobaric hypoxia (10% O2) for 3 weeks (Supplementary Figure S1A). 2 weeks after SU5416 injection, rats started receiving daily intraperitoneal injections of vehicle, 5Z-7-oxozeaenol, riociguat or combined therapy for 1 week. Body weight in the SuHyp group was significantly smaller than in the normoxic control group (Supplementary Figure S1B). Treatment with 5Z-7-oxozeaenol, riociguat or combined therapy did not affect body weight compared to vehicle-treated SuHyp rats.
Effects of riociguat and 5Z-7-oxozeaenol on pulmonary and cardiac hemodynamics
Three weeks following SU5416 administration and chronic exposure to hypoxia, SuHyp rats developed PH with PAP mean values of 43 ± 2 mmHg (Figure 2). Some animals were sacrificed 2 weeks after administration of SU5416 to confirm the onset of development of PH prior to initiating interventional treatments (Mean PAP values 12 ± 0.7 mmHg vs. 23 ± 0.7 in the normoxic and SuHyp groups; n = 3; p < 0.005). Treatment with 5Z-7-oxozeaenol for 1 week had no effect on PAP (Figures 2A–D) but attenuated RV dysfunction, as evidenced by a significant reduction in RVSP and RVEDP (Figures 2E,F). In contrast, riociguat significantly reduced PAP and RVSP but did not modify RVEDP. Combined treatment with riociguat and 5Z-7-oxozeaenol resulted in a significant reduction in all these parameters but they remained above control values and no additive effects were observed. None of the treatments induced changes in heart rate (data not shown).
FIGURE 2
FIGURE 2. Treatment with riociguat alone or combined with 5-z-7oxozeaenol attenuates pulmonary hypertension in the Sugen 5426/hypoxia (SuHyp) rat model. (A) Mean pulmonary arterial pressure (mPAP), (B) diastolic PAP (dPAP) and (C) systolic PAP (sPAP) in control (NMX-Vh) and SuHyp rats treated with vehicle, 5Z-7-oxozeaenol (SuHyp-OXO; 3 mg·kg–1·day–1), riociguat (SuHyp-RIO; 3 mg·kg–1·day–1) or both drugs combined (SuHyp-OXO-RIO). (D) Representative PAP recordings in each treatment group. (E) Right systolic and (F) end-diastolic ventricular pressure (RVSP and RVEDP). (G) Fulton index (ratio between RV and left ventricle plus septum weight) and (H) right ventricular weight relative to total body weight in the different treatment groups. Each bar shows the mean +SEM (n = 6–10 animals in each group). *p < 0.05 versus NMX-Vh and #p < 0,05 versus SuHyp vehicle (One-way ANOVA followed by Bonferroni’s post hoc test).
Effects of riociguat and 5Z-7-oxozeaenol on cardiac hypertrophy and cardiac metabolic reprogramming
In line with the hemodynamic alterations described above, SuHyp rats developed a marked RV hypertrophy (Figures 2G,H). Riociguat induced a minor but significant reduction in the Fulton index, whereas treatment with 5Z-7-oxozeaenol neither induced effects on its own nor modified the effects induced by riociguat (Figures 2G,H).
Metabolomic fingerprinting detected 10 metabolites significantly altered in the RV of SuHyp rats compared to control rats (Figure 3A and Supplementary Table S2). We found a full metabolic reprograming in SuHyp rats (Figure 3A) characterized by a shift in glucose metabolism, with higher metabolic concentrations of lactate and lower glucose; an increase of glutaminolysis with higher concentrations of glutamate and glutathione and lower concentration of glutamine; alterations in choline metabolism with a higher concentration of glycerophosphocholine (GPC) and lower of phosphocholine (PC). We also found an increase in taurine and GPC concentrations, whereas a decrease in inosine and creatine concentrations in SuHyp rats. PCA of RV NMR spectra provided perfect separation between SuHyp and control rats along the first two principal components (Figure 3B). The metabolites that changed between groups are highlighted in the PCA and PLS-DA loading plots shown in Supplementary Figures S2, 3. Metabolic pathway analysis was performed to identify the most relevant pathways involved in characteristic RV metabolic reprograming (Figure 3F). This pathway analysis identified alterations in energy metabolism, including glycolysis, Warburg effect and anaplerotic metabolism, amino acids biosynthesis and metabolism, glycerolipid metabolism and fatty acid degradation, and taurine metabolism.
FIGURE 3
FIGURE 3. Treatment with riociguat and 5-Z-7oxozeaenol reduces metabolic reprogramming in the right ventricle (RV) from SuHyp rats. (A) Metabolomic profiling of RV samples showing the quantification of metabolites significantly altered between the five intervention groups. Statistical significance was determined using a Bonferroni-corrected Student’s t-test (*p < 0.05 versus normoxia + vehicle and #p < 0.05 versus SuHyp + vehicle). (B) Metabolic profiling of RV samples discriminated between SuHyp and control groups as shown by principal components analysis (PCA) performed on 1H-MRS data. Treatment with (C) 5Z-7-oxozeaenol, (D) riociguat or (E) both drugs combined induced a significant metabolic shift in the SuHyp RV. (F) Metabolic pathway analysis of the set of metabolites found to be significantly dysregulated between control and SuHyp rats (G) Incidence of the combined treatment with 5Z-7-oxozeaenol and riociguat on the SuHyp group over the specific metabolic pathways. Y-axis represents the statistical p values from enrichment analysis, and the X-axis represents the pathway impact value calculated from pathway topology analysis. The node colors represent the p-values (lowest: light yellow; highest: dark red) and the node radius indicate the pathway impact values. Dysregulated metabolic pathways are labeled. Control group: NMX-Vh, Sugen 5,416/hypoxia: SuHyp, SuHyp + 5Z-7-oxozeaenol treatment: SuHyp-OXO; SuHyp + riociguat treatment: SuHyp-RIO; SuHyp+ 5Z-7-oxozeaenol and riociguat treatment: SuHyp-OXO-RIO, PC: Principal Component, GPC: Glycero-phosphocholine and PC: Phosphocholine.
Treatment with 5z-7-oxozeaenol, riociguat or both drugs combined, induced an apparent metabolic normalization in RV tissue (Figure 3A). Thus, PCA of RV NMR spectra provided clear discrimination between SuHyp and treated rats (Figures 3C–E). Treatment with 5Z-7-oxozeaenol induced the normalization of the metabolic levels of glucose, GPC, and PC whereas riociguat treatment induced a trend to normalize inosine, glucose, and GPC metabolic levels (Figure 3A). The combined treatment produced a further metabolic regularization characterized by the normalization of inosine, glucose, creatine, glutamine, GPC and PC metabolic levels. Finally, treatments with riociguat or with the combination of riociguat and 5Z-7-oxozeaenol induced an increment in taurine levels (Figure 3A). Metabolic pathway analysis highlighted that the combination of 5Z-7-oxozeaenol with riociguat induced a significant metabolic effect over taurine metabolism, amino acids biosynthesis and metabolism, glycolysis and TCA cycle metabolism via glycine, serine-threonine metabolism (Figure 3G).
Effects of riociguat and 5Z-7-oxozeaenol on pulmonary vascular remodeling and lung metabolic reprogramming
Lung histological evaluation confirmed the development of pulmonary vascular remodeling in SuHyp rats, with an increase in the percentage of muscular PA and the corresponding decrease in partially and non-muscular PA (Figures 4A,B) and a significant increase in the medial wall thickness of distal PA (Figure 4C). In contrast to our hemodynamic findings, combined therapy had an additive protective effect on pulmonary vascular remodeling. Thus, treatment with riociguat or 5Z-7-oxozeaenol significantly reduced the percentage of muscular PA but only combined treatment with riociguat plus 5Z-7-oxozeaenol significantly increased the percentage of non-muscular PA and reduced the medial wall thickness in SuHyp rats (Figure 4).
FIGURE 4
FIGURE 4. Riociguat and 5-z-7oxozeaenol exert additive effects on preventing pulmonary vascular remodeling. (A) Percentage of muscular, partially muscular and non-muscular arteries in the different treatment groups. (B) Medial wall thickness of pulmonary arteries (50–200 μm external diameter). (C) Representative images of lung histology. Each bar shows the mean +SEM (n = 5-6 animals in each group). *p < 0.05 versus normoxia + vehicle and #p < 0.05 versus hypoxia + vehicle (One-way ANOVA followed by Bonferroni’s post hoc test).
Metabolomic fingerprinting of lung samples detected eight metabolites significantly altered (Figure 5A and Supplementary Table S2), enabling the discrimination between SuHyp and control groups (Figure 5B). The PCA and PLS-DA loading plots highlighted the most significant metabolites that distinguish between groups (Supplementary Figures S4, 5). Specifically, we observed a reduction in the concentration of glucose, glycine and free fatty acids, and an increase in alanine, glutamate, glutamine, creatine, threonine, myo-inositol, taurine, carnitine, valine, 2-hydroxybutyrate and leucine (Figure 5A). However, the metabolic impact of the intervention protocols in the lung (Figures 5C–E) was less than that detected in the RV tissue (Figure 3). 5Z-7-oxozeaenol treatment induced a trend to the normalization of glutamine levels and 2-hydroxybutyrate, whereas riociguat treatment induced the normalization of myo-inositol and 2-hydroxybutyrate. 5Z-7-oxozeaenol and riociguat combined treatment induced a significant trend to the normalization of glucose and 2-hydroxybutyrate concentrations. In addition, a significant decrease in glycine and free fatty acids, and an increase in lactate were detected after the treatment with riociguat, either alone or in combination with 5Z-7-oxozeaenol (Figure 5A).
FIGURE 5
FIGURE 5. Treatment with riociguat and 5-z-7oxozeaenol has a modest effect on the lung metabolic reprogramming in the SuHyp rats. (A) Metabolomic profiling of lung samples showing the quantification of metabolites significantly altered between the five intervention groups. Statistical significance was determined using a Bonferroni-corrected Student’s t-test (*p < 0.05 versus normoxia vehicle and #p < 0.05 versus SuHyp vehicle). (B) Metabolic profiling of lung samples discriminated between SuHyp and control groups as shown by principal components analysis (PCA) performed on 1H-MRS data. Treatment with (C) 5Z-7-oxozeaenol, (D) riociguat or (E) both drugs induced a modest effect on the metabolic changes observed in the lungs from SuHPX rats. Control group: NMX-Vh, Sugen 5,416/hypoxia: SuHyp, SuHyp + 5Z-7-oxozeaenol treatment: SuHyp-OXO; SuHyp + riociguat treatment: SuHyp-RIO; SuHyp + 5Z-7-oxozeaenol and riociguat treatment: SuHyp-OXO-RIO, PC: Principal Component, FA: Fatty Acids.
Effects of riociguat and 5Z-7-oxozeaenol on pulmonary vascular dysfunction
PA from SuHyp rats exhibited a significant increase in the contractile responses induced by serotonin and a marked reduction in the relaxant responses induced by acetylcholine (Figure 6). Administration of riociguat, 5Z-7-oxozeaenol or the combination of both drugs did not reduce the hyperresponsiveness to serotonin (Figure 6A). In contrast, the combination of riociguat and 5Z-7-oxozeaenol attenuated the endothelial dysfunction induced by SuHyp, as evidenced by a significant increase in the relaxation induced by acetylcholine (Figure 6B).
FIGURE 6
FIGURE 6. Riociguat plus 5-z-7oxozeaenol combination attenuates pulmonary endothelial dysfunction in rats with PH. (A) Contractile responses induced by serotonin is isolated intrapulmonary arteries from control (normoxic) rats and rats with PAH (SuHyp) treated with vehicle (Vh), 5-z-7-oxozeaenol (OXO), riociguat (RIO) or a combination of both drugs (OXO + RIO). (B) Concentration-dependent relaxant responses induced by the endothelium-dependent vasodilator acetylcholine in pulmonary arteries isolated from the different experimental groups. N = 12-15, *p < 0.05 versus normoxia and # vs. SuHyp + vehicle (two-way ANOVA followed by Bonferroni’s post hoc test).
Effects of riociguat and 5Z-7-oxozeaenol on molecular markers of pulmonary hypertension
To gain insight into the underlying pathways mediating the changes induced by our intervention protocol, we assessed the expression of molecular markers involved in the pathophysiology of PH, including those associated with the TGFβ signaling pathway, inflammation, or ionic remodeling (Rabinovitch, 2012; Humbert et al., 2019). As previously reported in the monocrotaline and chronic hypoxia PH rat models (Nasim et al., 2012), lungs from SuHyp rats also displayed reduced protein levels of BMPR2 and increased activated TAK-1 (Figures 7A,B). As expected, treatment with 5Z-7-oxozeaenol but not riociguat significantly reduced the levels of activated TAK1 (Figures 7B,C). In contrast, administration of riociguat but not 5Z-7-oxozeaenol resulted in an increase in lung BMPR2 levels, reaching values above those found in control animals (Figure 7A). Downregulation of voltage-gated K+ channels (KV) channels (notably KV 1.5), is also considered to be an early contributor to the pathophysiology of PAH leading to vasoconstriction of PA, hypertrophy and remodeling (Antigny et al., 2016; Thenappan et al., 2018). In line with previous findings (Wang et al., 2005; Moral-Sanz et al., 2012; Morales-Cano et al., 2014; Mondejar-Parreno et al., 2019), KV 1.5 protein expression was significantly reduced in the lungs from SuHyp rats (Figure 7D). However, none of the pharmacological treatments was able to restore the expression of these channels. Recently, the relevance of an immunoinflammatory component (characterized by infiltration and activation of macrophages and various T-lymphocyte subpopulations) in the development of PAH has also been demonstrated (Rabinovitch, 2012; Nicolls and Voelkel, 2017; Koudstaal et al., 2020). Thus, increased proinflammatory cytokines (IL-1β, IL-6, IFNγ or IL-17) have been associated with reduced quality of life and increased mortality (Chaouat et al., 2009; Pellicelli et al., 2010; Soon et al., 2010). Accordingly, lungs from SuHyp rats exhibited significantly increased levels of IL-1β (Figure 7E). Treatment with 5Z-7-oxozeaenol induced a modest reduction of IL-1β whereas treatment with riociguat, either alone or combined with the TAK1 inhibitor, fully abolished the increase in pulmonary IL-1β. Next, we examined the activation of two common inflammatory signaling pathways, nuclear factor-kappa B (NF-κB) and mitogen-activated protein kinases (MAPK) pathways (Shi and Sun, 2018). To evaluate the activation of MAPK and NF-κB signaling pathways, we measured p38, extracellular signal-regulated kinases (ERK), c-Jun N-terminal kinases (JNK) and NF-κB p65 protein phosphorylation by Western blot analysis. As shown in Figure 7F, only JNK was found to be activated in the lungs from SuHyp rats but treatment with 5Z-7-oxozeaenol or riociguat further increased the phosphorylation rate of this kinase. A rising trend in ERK1 and NF-κB p65 subunit phosphorylation was also observed in the lungs from SuHyp rats but not in those treated with 5Z-7-oxozeaenol or riociguat.
FIGURE 7
FIGURE 7. Treatment with riociguat and 5-z-7oxozeaenol restores the expression of BMPR2 and inhibits IL-β in the lungs from SuHyp rats. Densitometry analysis of (A) BMPR2/β-actin and (B) pTAK1/TAK1 protein expression in lung homogenates from control (normoxic) rats and rats with PAH (SuHyp) treated with vehicle (Vh), 5-z-7-oxozeaenol (OXO), riociguat (RIO) or a combination of both drugs (OXO + RIO). (C) Representative image of Western blot showing p-TAK1, TAK1 and β-actin. (D) Densitometry analysis of Kv1.5/β-actin in lung homogenates. (E) Quantification of IL-1β levels in lung homogenates in the five experimental groups. (F) Heatmap showing the densitometry analysis of Western blot performed to assess the activation of the NF-κB and MAPK signaling pathway. N = 6-10 animals in each group, *p < 0.05 versus normoxia and # vs. SuHyp + vehicle (One sample t-test).
Effects of 5Z-7-oxozeaenol and riociguat on the metabolic reprogramming in PASMCs
To further characterize the metabolic effects of 5Z-7-oxozeaenol and riociguat, primary cultures of PASMC were established from SuHyp rats and treated in vitro with both drugs. Treatment with 5Z-7-oxozeaenol, significantly attenuated the metabolic shift induced by SuHyp. Specifically, alanine, aspartate, creatine, creatinine, glutamate, glutamine, glycine, inosine monophosphate, 2-methylglutarate, isoleucine, lactate, PC, threonine and valine metabolic pools were found significantly decreased after the exposition to 5Z-7-oxozeaenol or the combination of 5Z-7-oxozeaenol and riociguat (Figures 8A,B). 5Z-7-oxozeaenol and riociguat combined treatment also induced a decrease in myo-inositol and proline concentrations. In contrast, treatment with riociguat only induced a decrease in lactate concentration (Figure 8A). Metabolic pathway analysis was performed to identify the most relevant pathways affected by 5Z-7-oxozeaenol and riociguat (Figure 8C). This pathway analysis identified alterations in energy metabolism including glycolysis and pyruvate metabolism, amino acids biosynthesis, TCA cycle metabolism via glycine-serine-threonine metabolism, and anaplerotic metabolism, including glutaminolysis, and ketone bodies metabolism.
FIGURE 8
FIGURE 8. Ex vivo treatment with 5-z-7oxozeaenol attenuates the metabolic reprogramming in cultured pulmonary artery smooth muscle cells (PASMCs) from SuHyp rats. (A,B) Metabolomic profiling of cell extracts showing the quantification of metabolites significantly altered in PASMCs isolated from SuHyp rats and treated with vehicle (DMSO), 5Z-7-oxozeaenol (OXO; 1 µM), riociguat (RIO; 1 µM) or both drugs combined (OXO + RIO). Statistical significance was determined using a Bonferroni-corrected Student’s t-test (*p < 0.05 versus normoxia + vehicle and #p < 0.05 versus SuHyp + vehicle). (C) Metabolic pathway analysis of the set of metabolites found to be significantly dysregulated following treatment with the combination of OXO and RIO on PASMCs isolated from SuHyp rats. Y-axis represents the statistical p values from enrichment analysis, and the X-axis represents the pathway impact value calculated from pathway topology analysis. The node colors represent the p-values (lowest: light yellow; highest: dark red) and the node radius indicate the pathway impact values. Dysregulated metabolic pathways are labeled.
Discussion
Although several new therapies have become available over the last 2 decades, there is no cure for PAH, and the prognosis remains poor. In this study, we investigated the therapeutic efficacy of the combination of the TAK-1 inhibitor 5Z-7-oxozeaenol and the sGC stimulator riociguat in the experimental model of PH induced by the combined exposure to the VEGF receptor antagonist SU5416 and chronic hypoxia. We first demonstrated that in vitro, riociguat acts mainly as a vasodilator agent whereas 5Z-7-oxozeaenol is an effective antiproliferative drug. In vivo, treatment with 3 mg/kg riociguat for 1 week induced a modest reduction in PAP, pulmonary vascular remodeling and RV hypertrophy accompanied by an increase in BMPR2 expression and a reduction in lung inflammation. Riociguat also induced a partial reduction in metabolic reprogramming of RV and PASMCs. By contrast, the TAK-1 inhibitor 5Z-7-oxozeaenol significantly reduced pulmonary vascular remodeling and limited lung inflammation, induced a modest improvement in vascular remodeling and attenuated the metabolic reprogramming in RV and PASMCs in the rat SuHyp model of PAH. However, 5Z-7-oxozeaenol failed to reduce PAP or RV hypertrophy, did not increase BMPR2 expression and only had a modest metabolic effect in lung tissue. Importantly, combination therapy significantly improved hemodynamic parameters and cardiac hypertrophy, but no additive effects were observed. In contrast, the combination of riociguat and 5Z-7-oxozeaenol was associated with a greater reduction in pulmonary vascular remodeling and endothelial dysfunction than those observed with either drug alone. Moreover, the combined treatment was also associated with more pronounced metabolic effects, leading to the normalization of the metabolic dysregulation found in the RV and PASMCs of PH rats.
Riociguat isa sGC stimulator and is the only drug approved for both PAH and chronic thromboembolic pulmonary hypertension (CTEPH) (Ghofrani et al., 2013; Simonneau et al., 2015). sGC stimulators can directly stimulate the sGC enzyme in a NO-independent fashion and produce their sensitization to low endogenous NO levels, increasing cGMP biosynthesis and leading to vascular relaxation and antiproliferative effects (Cogolludo et al., 2007). In the present study, we confirmed the ability of riociguat to induce almost full relaxation in isolated rat PA. These results are in line with our previous report showing the efficacy of riociguat as a vasodilator in isolated human and rat PA under various oxygen conditions (Morales-Cano et al., 2021). Surprisingly, riociguat had no significant effects in the proliferative ability of PASMCs from neither control nor SuHyp rats nor human PASMCs (unpublished observations). In vitro antiproliferative effects by other direct stimulators of sGC, such as YC-1 or BAY 41-2,272, have been previously demonstrated (Huh et al., 2011; Joshi et al., 2011). However, studies assessing the in vitro antiproliferative capacity of riociguat are scarce. Thus, Patel et al., reported the ability of riociguat to inhibit ET-1-induced growth of PASMCs derived from PAH patients (Patel et al., 2014). In contrast, other groups, including ours, did not find a significant attenuation of cell proliferation at concentration at which riociguat induces strong vasodilator effects (Morales-Cano et al., 2015; Nawa et al., 2016). In contrast, we found that low-dose of riociguat (3 mg/kg/day) reduced the increase in PAP, RV hypertrophy, and pulmonary vascular remodeling in the rat SuHyp model of PAH. These findings are consistent with those reported by other groups in different experimental models of PH (Schermuly et al., 2008; Lang et al., 2012; Rai et al., 2018) but suggest that the ability of riociguat to reduce the progression of pulmonary vascular remodeling and RV hypertrophy are strongly associated with its ability to reduce PAP rather than with direct antiproliferative effects.
5Z-7-oxozeaenol, a resorcylic acid lactone derived from fungus, is a potent and selective inhibitor of TAK1 (Ninomiya-Tsuji et al., 2003). TAK1, also known as MAP3K7, is a key member of the TGFβ non-canonical signaling pathway (Sato et al., 2005; Mihaly et al., 2014). A role for TAK1 in the pathogenesis of PAH has been proposed (Nasim et al., 2012; Dalvi et al., 2017). Thus, Nassim et al., demonstrated that BMPR-II dysfunction promotes the activation of TAK1 leading to aberrant activation of the TGFβ-MAPK axis in the lung and resulting in a pro-proliferative and anti-apoptotic response (Nasim et al., 2012). Notably, TAK1 expression and activity were shown to be increased in murine PASMCs carrying a PAH-associated bmpr2 non-sense mutation but also in hypoxia and monocrotaline-induced PAH rat models, which display reduced levels of bmpr2 transcripts. Prof. Schneider’s group showed that an activating mutation of TAK1 is sufficient to induce myocardial hypertrophy and fulminant heart failure (Zhang et al., 2000). In the present study, we found that phosphorylated TAK1 is also increased in the lungs from SuHyp rats and that treatment with the TAK1 inhibitor 5Z-7-Oxozeanol significantly reduced PASMCs proliferation in vitro and limited pulmonary vascular remodeling in vivo. However, 5Z-7-Oxozeanol lacked acute vasodilator effects and it failed to reduce PAP or to limit RV hypertrophy in the SuHyp model. These results align with those found with sotatercept analogues (Yung et al., 2020; Joshi et al., 2022). Thus, despite neutralizing antibodies against activins induced modest effects on PAP, treatment with these analogues caused regression of vascular remodeling either as monotherapy or combined with a vasodilator.
Our results are also consistent with findings from other groups demonstrating the antiproliferative effects of TAK-1 inhibitors in several cell types (Nasim et al., 2012; Zippel et al., 2013; Song et al., 2014; Chen et al., 2016; Pandolfi et al., 2017) in the absence of significant hemodynamic effects (Garfield et al., 2019). Thus, Garfield et al. recently reported that 5Z-7-oxozeaenol had no effect on RVSP or Fulton index but prevented muscle loss in monocrotaline-treated rats (Garfield et al., 2019). Since the interventional treatments were initiated after the onset of pulmonary hypertension and were carried out only for 1 week, it remains to be determined whether treatment over a longer period of time would induce more robust hemodynamic effects. However, our study identifies novel protective effects induced by TAK-1 inhibition in the RV metabolic changes which could also contribute to alleviate PH-related RV dysfunction (Figures 2, 3).
Metabolic reprogramming in PAH is a major factor in the pathogenesis of the pulmonary vascular disease (Assad and Hemnes, 2015; Xu et al., 2021). The pulmonary vasculature in PAH shows normoxic activation of hypoxia-inducible factor 1-alpha (HIF-1α), which creates a ¨pseudo-hypoxic¨ environment despite normal oxygen supply T (Ryan and Archer, 2015). HIF-1α activation induces a metabolic shift toward aerobic glycolysis and a reduction of glucose flux into the mitochondria. To maintain mitochondria activity, metabolism through anaplerotic pathways increases, including fatty acid oxidation and choline-conjugated phospholipid metabolism (Sutendra et al., 2010; Kudryashova et al., 2015) and glutaminolysis (Egnatchik et al., 2017). In addition, RV hypertrophy induced by the sustained PH is characterized by a similar metabolic reprograming (Graham et al., 2018), although these changes are likely induced by ischemia rather than by impaired oxygen sensing (Ryan and Archer, 2015). Our findings in the rat SuHyp model reproduce the main metabolic alterations previously reported in lung and RV tissue in several PAH animal models (Zungu-Edmondson et al., 2017; Graham et al., 2018; Izquierdo-Garcia et al., 2018; Sakao et al., 2021). In the lung tissue, we confirmed: I) the enhancement of the glycolytic pathway, including lower glucose and higher lactate concentrations; II) alteration of one-carbon metabolism including higher glycine metabolism; III) upregulation of anaplerotic pathways, including glutaminolysis and upregulation of fatty acid oxidation. In RV tissue, we reproduced: I) the Warburg effect characterized by lower glucose and higher lactate concentrations; II) increase of anaplerotic metabolism including amino acids metabolism, glycerolipid metabolism and fatty acid degradation; III) and taurine metabolism.
Our data revealed that the impact of treatment with riociguat on lung tissue metabolism is limited inducing only a significant reduction in the concentration of myo-inositol, 2-hydroxybutyrate and leucine in the lungs from SuHyp rats. Myo-inositol has been proposed to modulate reduced interleukin-6 levels (Bizzarri et al., 2020), a cytokine involved in the development and progression of (Steiner et al., 2009), 2-hydroxybutyrate is a marker of the ketone metabolism that characterizes PAH (Mey et al., 2020), whereas leucine reduction may be related to glutamine metabolism (Bertero et al., 2019). The effect of riociguat in cultured rat PASMC is also scarce since we only found a reduction in lactate concentration but no other metabolites of the glycolytic pathway were normalized. In contrast, treatment with riociguat induced protective effects in RV metabolic disruption, inducing the normalization of choline-conjugated phospholipid metabolism, inosine, creatine, and glucose levels. Since heart failure is characterized by reduced myocardial inosine and creatine (Archer et al., 1995; Luptak et al., 2018), the normalization of these metabolic concentrations may be markers of an improvement of heart function following treatment with riociguat. To the best of our knowledge, this is the first report demonstrating a protective effect of riociguat on PAH metabolic alterations.
Treatment with 5Z-7-oxozeaenol in rat PASMC induced a significant effect in several metabolic pathways, including in glycolysis, citrate cycle, glutaminolysis, or branched-chain amino acids biosynthesis (Figure 8C), suggesting a protective effect on the metabolic dysregulation characteristic found in PH patients. For example, the reduction of lactate, and alanine concentrations is related to the normalization of the glycolysis pathway. Lower glutamine, glutamate and PC concentrations results from a lesser requirement of metabolic supplies by anaplerotic pathways. However, treatment with 5Z-7-oxozeaenol only normalized glutamine levels in lung tissue suggesting that its effects are more moderate in vivo. By contrast, 5Z-7-oxozeaenol induced the normalization of glucose and choline-conjugated phospholipid metabolism in RV tissue.
The combined treatment of 5Z-7-oxozeaenol and riociguat reinforces the metabolic effects of both treatments separately. In rat PASMC, combined treatment reproduces the metabolic alteration caused by 5Z-7-oxozeaenol treatment, but we also observed a statically significant reduction in glutamate, myo-inositol and proline concentrations. In lung tissue, combined treatment did not improve the metabolic effect of riociguat treatment alone. However, the combined treatment potentiates the metabolic effect induced by each drug alone in the RV tissue. Indeed, the levels of creatine and inosine from animals treated with 5Z-7-oxozeaenol and riociguat are more similar to the levels detected in normoxic animals than from animals treated with riociguat. These metabolic markers may reflect a better cardiac function in combined treatment animals.
In summary, the addition of the TAK1 inhibitor 5Z-7-oxozeaenol to riociguat significantly improved hemodynamic parameters and cardiac hypertrophy in the SuHyp rat model. Greater protective effects on pulmonary vascular remodeling and endothelial dysfunction were observed when both drugs were combined. Moreover, only the combination of riociguat with 5Z-7-oxozeaenol was able to normalize the metabolic dysregulation found in the RV rats with PAH suggesting that the inhibition of TAK1 might represent a cardioprotective target for the treatment of RV failure associated to PAH.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was reviewed and approved by Institutional (UCM) and Regional Ethical Committees of Madrid (Spain).
Author contributions
Involved in the conception, hypothesis delineation, and design of the study: DM-C, JI, JR-C, FP-V, LM; Performed the experiments or analysed the data: DM-C, JI, BB, MC, SE-R, RP, PV-V, IR. Interpretation of the data: DM-C, JI, AC, JR-C, FP-V, LM. Wrote the paper or made substantial contribution prior to submission: DM-C, JI, JR-C, FP-V, LM.
Funding
This work was supported by the Instituto de Salud Carlos III-ISCIII (Grant numbers: PI15/01100 and PI19/01616 to LM), the Spanish Ministry of Science and Innovation MCIN (Grant numbers: PID 2019-107363RB-I00 to FP-V, PID 2020-117939RB-I00 to AC and PID 2021-123238OB-I00, PDC 2021-121696-I00 to JR-C and PID2019-106564RJ-I00 to JI-G), the Comunidad de Madrid-CAM (CM S2017/BMD-3727 to AC and LM and B2017/BMD3875 to JI-G) and, as appropriate, by “ERDF A way of making Europe”, co-funded by the “European Union”. FP-V received funding from Fundación Contra la Hipertensión Pulmonar (Empathy grant) and JR-C from La Caixa Foundation (Health Research Call 2020: HR20-00075). This work was performed under the Maria de Maeztu Units of Excellence Programme–Grant MDM-2017-0720 funded by MCIN/AEI/10.13039/501100011033.
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.
Publisher’s note
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2023.1021535/full#supplementary-material
References
Aittokallio, T., and Schwikowski, B. (2006). Graph-based methods for analysing networks in cell biology. Brief. Bioinform 7, 243–255. doi:10.1093/bib/bbl022
Andre, P., Joshi, S. R., Briscoe, S. D., Alexander, M. J., Li, G., and Kumar, R. (2021). Therapeutic approaches for treating pulmonary arterial hypertension by correcting imbalanced TGF-beta superfamily signaling. Front. Med. (Lausanne) 8, 814222. doi:10.3389/fmed.2021.814222
Antigny, F., Hautefort, A., Meloche, J., Belacel-Ouari, M., Manoury, B., Rucker-Martin, C., et al. (2016). Potassium Channel subfamily K member 3 (KCNK3) contributes to the development of pulmonary arterial hypertension. Circulation 133, 1371–1385. doi:10.1161/CIRCULATIONAHA.115.020951
Archer, D. F., Fahy, G. E., Viniegra-Sibal, A., Anderson, F. D., Snipes, W., and Foldesy, R. G. (1995). Initial and steady-state pharmacokinetics of a vaginally administered formulation of progesterone. Am. J. Obstet. Gynecol. 173, 471–477; discussion 477-478. ; discussion 477-478. doi:10.1016/0002-9378(95)90268-6
Assad, T. R., and Hemnes, A. R. (2015). Metabolic dysfunction in pulmonary arterial hypertension. Curr. Hypertens. Rep. 17, 20. doi:10.1007/s11906-014-0524-y
Bertero, T., Perk, D., and Chan, S. Y. (2019). The molecular rationale for therapeutic targeting of glutamine metabolism in pulmonary hypertension. Expert Opin. Ther. Targets 23, 511–524. doi:10.1080/14728222.2019.1615438
Bizzarri, M., Lagana, A. S., Aragona, D., and Unfer, V. (2020). Inositol and pulmonary function. Could myo-inositol treatment downregulate inflammation and cytokine release syndrome in SARS-CoV-2? Eur. Rev. Med. Pharmacol. Sci. 24, 3426–3432. doi:10.26355/eurrev_202003_20715
Callejo, M., Mondejar-Parreno, G., Morales-Cano, D., Barreira, B., Esquivel-Ruiz, S., Olivencia, M. A., et al. (2020). Vitamin D deficiency downregulates TASK-1 channels and induces pulmonary vascular dysfunction. Am. J. Physiol. Lung Cell Mol. Physiol. 319, L627–L640. doi:10.1152/ajplung.00475.2019
Chaouat, A., Savale, L., Chouaid, C., Tu, L., Sztrymf, B., Canuet, M., et al. (2009). Role for interleukin-6 in COPD-related pulmonary hypertension. Chest 136, 678–687. doi:10.1378/chest.08-2420
Chen, Z., Mei, Y., Lei, H., Tian, R., Ni, N., Han, F., et al. (2016). LYTAK1, a TAK1 inhibitor, suppresses proliferation and epithelialmesenchymal transition in retinal pigment epithelium cells. Mol. Med. Rep. 14, 145–150. doi:10.3892/mmr.2016.5275
Chong, J., Soufan, O., Li, C., Caraus, I., Li, S., Bourque, G., et al. (2018). MetaboAnalyst 4.0: Towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 46, W486–W494. doi:10.1093/nar/gky310
Cogolludo, A., Moreno, L., Frazziano, G., Moral-Sanz, J., Menendez, C., Castaneda, J., et al. (2009). Activation of neutral sphingomyelinase is involved in acute hypoxic pulmonary vasoconstriction. Cardiovasc Res. 82, 296–302. doi:10.1093/cvr/cvn349
Cogolludo, A., Moreno, L., and Villamor, E. (2007). Mechanisms controlling vascular tone in pulmonary arterial hypertension: Implications for vasodilator therapy. Pharmacology 79, 65–75. doi:10.1159/000097754
Dalvi, P., Sharma, H., Konstantinova, T., Sanderson, M., Brien-Ladner, A. O., and Dhillon, N. K. (2017). Hyperactive TGF-beta signaling in smooth muscle cells exposed to HIV-protein(s) and cocaine: Role in pulmonary vasculopathy. Sci. Rep. 7, 10433. doi:10.1038/s41598-017-10438-3
Dunnett, C. W. (1955). A multiple comparison procedure for comparing several treatments with a control. J. Am. Stat. Assoc. 50, 1096–1121. doi:10.1080/01621459.1955.10501294
Egnatchik, R. A., Brittain, E. L., Shah, A. T., Fares, W. H., Ford, H. J., Monahan, K., et al. (2017). Dysfunctional BMPR2 signaling drives an abnormal endothelial requirement for glutamine in pulmonary arterial hypertension. Pulm. Circ. 7, 186–199. doi:10.1086/690236
Freund-Michel, V., Khoyrattee, N., Savineau, J. P., Muller, B., and Guibert, C. (2014). Mitochondria: Roles in pulmonary hypertension. Int. J. Biochem. Cell Biol. 55, 93–97. doi:10.1016/j.biocel.2014.08.012
Garfield, B. E., Crosby, A., Shao, D., Yang, P., Read, C., Sawiak, S., et al. (2019). Growth/differentiation factor 15 causes TGFβ-activated kinase 1-dependent muscle atrophy in pulmonary arterial hypertension. Thorax 74, 164–176. doi:10.1136/thoraxjnl-2017-211440
Ghofrani, H. A., Galie, N., Grimminger, F., Grunig, E., Humbert, M., Jing, Z. C., et al. (2013). Riociguat for the treatment of pulmonary arterial hypertension. N. Engl. J. Med. 369, 330–340. doi:10.1056/NEJMoa1209655
Graham, B. B., Kumar, R., Mickael, C., Kassa, B., Koyanagi, D., Sanders, L., et al. (2018). Vascular adaptation of the right ventricle in experimental pulmonary hypertension. Am. J. Respir. Cell Mol. Biol. 59, 479–489. doi:10.1165/rcmb.2018-0095OC
Hernandez-Saavedra, D., Sanders, L., Freeman, S., Reisz, J. A., Lee, M. H., Mickael, C., et al. (2020). Stable isotope metabolomics of pulmonary artery smooth muscle and endothelial cells in pulmonary hypertension and with TGF-beta treatment. Sci. Rep. 10, 413. doi:10.1038/s41598-019-57200-5
Hotelling, H. (1933). Analysis of a complex of statistical variables into principal components. J. Educ. Psychol. 24, 417–441. doi:10.1037/h0071325
Huh, J. W., Kim, S. Y., Lee, J. H., and Lee, Y. S. (2011). YC-1 attenuates hypoxia-induced pulmonary arterial hypertension in mice. Pulm. Pharmacol. Ther. 24, 638–646. doi:10.1016/j.pupt.2011.09.003
Humbert, M., Guignabert, C., Bonnet, S., Dorfmuller, P., Klinger, J. R., Nicolls, M. R., et al. (2019). Pathology and pathobiology of pulmonary hypertension: State of the art and research perspectives. Eur. Respir. J. 53, 1801887. doi:10.1183/13993003.01887-2018
Humbert, M., McLaughlin, V., Gibbs, J. S. R., Gomberg-Maitland, M., Hoeper, M. M., Preston, I. R., et al. (2021). Sotatercept for the treatment of pulmonary arterial hypertension. N. Engl. J. Med. 384, 1204–1215. doi:10.1056/NEJMoa2024277
Hummel, M., Meister, R., and Mansmann, U. (2008). GlobalANCOVA: Exploration and assessment of gene group effects. Bioinformatics 24, 78–85. doi:10.1093/bioinformatics/btm531
Izquierdo-Garcia, J. L., Arias, T., Rojas, Y., Garcia-Ruiz, V., Santos, A., Martin-Puig, S., et al. (2018). Metabolic reprogramming in the heart and lung in a murine model of pulmonary arterial hypertension. Front. Cardiovasc Med. 5, 110. doi:10.3389/fcvm.2018.00110
Izquierdo-Garcia, J. L., Nin, N., Ruiz-Cabello, J., Rojas, Y., de Paula, M., Lopez-Cuenca, S., et al. (2011). A metabolomic approach for diagnosis of experimental sepsis. Intensive Care Med. 37, 2023–2032. doi:10.1007/s00134-011-2359-1
Joshi, C. N., Martin, D. N., Fox, J. C., Mendelev, N. N., Brown, T. A., and Tulis, D. A. (2011). The soluble guanylate cyclase stimulator BAY 41-2272 inhibits vascular smooth muscle growth through the cAMP-dependent protein kinase and cGMP-dependent protein kinase pathways. J. Pharmacol. Exp. Ther. 339, 394–402. doi:10.1124/jpet.111.183400
Joshi, S. R., Liu, J., Bloom, T., Karaca Atabay, E., Kuo, T. H., Lee, M., et al. (2022). Sotatercept analog suppresses inflammation to reverse experimental pulmonary arterial hypertension. Sci. Rep. 12, 7803. doi:10.1038/s41598-022-11435-x
Kankainen, M., Gopalacharyulu, P., Holm, L., and Oresic, M. (2011). MPEA-metabolite pathway enrichment analysis. Bioinformatics 27, 1878–1879. doi:10.1093/bioinformatics/btr278
Kilkenny, C., Browne, W. J., Cuthill, I. C., Emerson, M., and Altman, D. G. (2010). Improving bioscience research reporting: The ARRIVE guidelines for reporting animal research. PLoS Biol. 8, e1000412. doi:10.1371/journal.pbio.1000412
Koudstaal, T., Boomars, K. A., and Kool, M. (2020). Pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension: An immunological perspective. J. Clin. Med. 9, 561. doi:10.3390/jcm9020561
Kudryashova, T. V., Goncharov, D. A., Pena, A., Ihida-Stansbury, K., DeLisser, H., Kawut, S. M., et al. (2015). Profiling the role of mammalian target of rapamycin in the vascular smooth muscle metabolome in pulmonary arterial hypertension. Pulm. Circ. 5, 667–680. doi:10.1086/683810
Lang, M., Kojonazarov, B., Tian, X., Kalymbetov, A., Weissmann, N., Grimminger, F., et al. (2012). The soluble guanylate cyclase stimulator riociguat ameliorates pulmonary hypertension induced by hypoxia and SU5416 in rats. PLoS One 7, e43433. doi:10.1371/journal.pone.0043433
Liu, G., Lee, D. P., Schmidt, E., and Prasad, G. L. (2019). Pathway analysis of global metabolomic profiles identified enrichment of caffeine, energy, and arginine metabolism in smokers but not moist snuff consumers. Bioinform Biol. Insights 13, 1177932219882961. doi:10.1177/1177932219882961
Luptak, I., Sverdlov, A. L., Panagia, M., Qin, F., Pimentel, D. R., Croteau, D., et al. (2018). Decreased ATP production and myocardial contractile reserve in metabolic heart disease. J. Mol. Cell Cardiol. 116, 106–114. doi:10.1016/j.yjmcc.2018.01.017
Marsboom, G., Toth, P. T., Ryan, J. J., Hong, Z., Wu, X., Fang, Y. H., et al. (2012). Dynamin-related protein 1-mediated mitochondrial mitotic fission permits hyperproliferation of vascular smooth muscle cells and offers a novel therapeutic target in pulmonary hypertension. Circ. Res. 110, 1484–1497. doi:10.1161/CIRCRESAHA.111.263848
Meijles, D. N., Cull, J. J., Markou, T., Cooper, S. T. E., Haines, Z. H. R., Fuller, S. J., et al. (2020). Redox regulation of cardiac ASK1 (apoptosis signal-regulating kinase 1) controls p38-MAPK (Mitogen-Activated protein kinase) and orchestrates cardiac remodeling to hypertension. Hypertension 76, 1208–1218. doi:10.1161/HYPERTENSIONAHA.119.14556
Mey, J. T., Hari, A., Axelrod, C. L., Fealy, C. E., Erickson, M. L., Kirwan, J. P., et al. (2020). Lipids and ketones dominate metabolism at the expense of glucose control in pulmonary arterial hypertension: A hyperglycaemic clamp and metabolomics study. Eur. Respir. J. 55, 1901700. doi:10.1183/13993003.01700-2019
Mihaly, S. R., Ninomiya-Tsuji, J., and Morioka, S. (2014). TAK1 control of cell death. Cell Death Differ. 21, 1667–1676. doi:10.1038/cdd.2014.123
Mondejar-Parreno, G., Callejo, M., Barreira, B., Morales-Cano, D., Esquivel-Ruiz, S., Moreno, L., et al. (2019). miR-1 is increased in pulmonary hypertension and downregulates Kv1.5 channels in rat pulmonary arteries. J. Physiol. 597, 1185–1197. doi:10.1113/JP276054
Montani, D., Chaumais, M. C., Guignabert, C., Gunther, S., Girerd, B., Jais, X., et al. (2014). Targeted therapies in pulmonary arterial hypertension. Pharmacol. Ther. 141, 172–191. doi:10.1016/j.pharmthera.2013.10.002
Moral-Sanz, J., Lopez-Lopez, J. G., Menendez, C., Moreno, E., Barreira, B., Morales-Cano, D., et al. (2012). Different patterns of pulmonary vascular disease induced by type 1 diabetes and moderate hypoxia in rats. Exp. Physiol. 97, 676–686. doi:10.1113/expphysiol.2011.062257
Morales-Cano, D., Barreira, B., Chamorro, V., Pandolfi, R., De Olaiz Navarro, B., Esteban, A., et al. (2015). Antiproliferative effects of drugs affecting different signalling pathways on rat and human pulmonary artery smooth muscle cells. Eur. Respir. J. 46, PA584.
Morales-Cano, D., Barreira, B., Navarro, B. O., Callejo, M., Mondejar-Parreno, G., Esquivel-Ruiz, S., et al. (2021). Oxygen-sensitivity and pulmonary selectivity of vasodilators as potential drugs for pulmonary hypertension. Antioxidants (Basel) 10, 155. doi:10.3390/antiox10020155
Morales-Cano, D., Menendez, C., Moreno, E., Moral-Sanz, J., Barreira, B., Galindo, P., et al. (2014). The flavonoid quercetin reverses pulmonary hypertension in rats. PLoS One 9, e114492. doi:10.1371/journal.pone.0114492
Nasim, M. T., Ogo, T., Chowdhury, H. M., Zhao, L., Chen, C. N., Rhodes, C., et al. (2012). BMPR-II deficiency elicits pro-proliferative and anti-apoptotic responses through the activation of TGFβ-TAK1-MAPK pathways in PAH. Hum. Mol. Genet. 21, 2548–2558. doi:10.1093/hmg/dds073
Nawa, N., Ishida, H., Katsuragi, S., Baden, H., Takahashi, K., Higeno, R., et al. (2016). Constitutively active form of natriuretic peptide receptor 2 ameliorates experimental pulmonary arterial hypertension. Mol. Ther. Methods Clin. Dev. 3, 16044. doi:10.1038/mtm.2016.44
Nicolls, M. R., and Voelkel, N. F. (2017). The roles of immunity in the prevention and evolution of pulmonary arterial hypertension. Am. J. Respir. Crit. Care Med. 195, 1292–1299. doi:10.1164/rccm.201608-1630PP
Ninomiya-Tsuji, J., Kajino, T., Ono, K., Ohtomo, T., Matsumoto, M., Shiina, M., et al. (2003). A resorcylic acid lactone, 5Z-7-oxozeaenol, prevents inflammation by inhibiting the catalytic activity of TAK1 MAPK kinase kinase. J. Biol. Chem. 278, 18485–18490. doi:10.1074/jbc.M207453200
Pandolfi, R., Barreira, B., Moreno, E., Lara-Acedo, V., Morales-Cano, D., Martinez-Ramas, A., et al. (2017). Role of acid sphingomyelinase and IL-6 as mediators of endotoxin-induced pulmonary vascular dysfunction. Thorax 72, 460–471. doi:10.1136/thoraxjnl-2015-208067
Patel, J., Abraham, D., Nelson, A., Silverstein, A., and Clapp, L. (2014). Comparison of current therapies against endothelin-induced growth in pulmonary smooth arterial muscle cells (PASMCs) derived from patients with pulmonary arterial hypertension. Eur. Respir. J. 44, P2355.
Pellicelli, A. M., Barbaro, G., Puoti, C., Guarascio, P., Lusi, E. A., Bellis, L., et al. (2010). Plasma cytokines and portopulmonary hypertension in patients with cirrhosis waiting for orthotopic liver transplantation. Angiology 61, 802–806. doi:10.1177/0003319710369101
Prins, K. W., Rose, L., Archer, S. L., Pritzker, M., Weir, E. K., Olson, M. D., et al. (2019). Clinical determinants and prognostic implications of right ventricular dysfunction in pulmonary hypertension caused by chronic lung disease. J. Am. Heart Assoc. 8, e011464. doi:10.1161/JAHA.118.011464
Prins, K. W., Thenappan, T., Weir, E. K., Kalra, R., Pritzker, M., and Archer, S. L. (2019). Repurposing medications for treatment of pulmonary arterial hypertension: What's old is new again. J. Am. Heart Assoc. 8, e011343. doi:10.1161/JAHA.118.011343
Rabinovitch, M. (2012). Molecular pathogenesis of pulmonary arterial hypertension. J. Clin. Invest. 122, 4306–4313. doi:10.1172/JCI60658
Rai, N., Veeroju, S., Schymura, Y., Janssen, W., Wietelmann, A., Kojonazarov, B., et al. (2018). Effect of riociguat and sildenafil on right heart remodeling and function in pressure overload induced model of pulmonary arterial banding. Biomed. Res. Int. 2018, 3293584. doi:10.1155/2018/3293584
Rol, N., Kurakula, K. B., Happe, C., Bogaard, H. J., and Goumans, M. J. (2018). TGF-Beta and BMPR2 signaling in PAH: Two black sheep in one family. Int. J. Mol. Sci. 19, 2585. doi:10.3390/ijms19092585
Ryan, J. J., and Archer, S. L. (2015). Emerging concepts in the molecular basis of pulmonary arterial hypertension: Part I: Metabolic plasticity and mitochondrial dynamics in the pulmonary circulation and right ventricle in pulmonary arterial hypertension. Circulation 131, 1691–1702. doi:10.1161/CIRCULATIONAHA.114.006979
Sakao, S., Kawakami, E., Shoji, H., Naito, A., Miwa, H., Suda, R., et al. (2021). Metabolic remodeling in the right ventricle of rats with severe pulmonary arterial hypertension. Mol. Med. Rep. 23, 227. doi:10.3892/mmr.2021.11866
Sato, S., Sanjo, H., Takeda, K., Ninomiya-Tsuji, J., Yamamoto, M., Kawai, T., et al. (2005). Essential function for the kinase TAK1 in innate and adaptive immune responses. Nat. Immunol. 6, 1087–1095. doi:10.1038/ni1255
Schermuly, R. T., Stasch, J. P., Pullamsetti, S. S., Middendorff, R., Muller, D., Schluter, K. D., et al. (2008). Expression and function of soluble guanylate cyclase in pulmonary arterial hypertension. Eur. Respir. J. 32, 881–891. doi:10.1183/09031936.00114407
Scott, T. E., Qin, C. X., Drummond, G. R., Hobbs, A. J., and Kemp-Harper, B. K. (2021). Innovative anti-inflammatory and pro-resolving strategies for pulmonary hypertension: High blood pressure research council of Australia award 2019. Hypertension 78, 1168–1184. doi:10.1161/HYPERTENSIONAHA.120.14525
Sharmin, N., Nganwuchu, C. C., and Nasim, M. T. (2021). Targeting the TGF-beta signaling pathway for resolution of pulmonary arterial hypertension. Trends Pharmacol. Sci. 42, 510–513. doi:10.1016/j.tips.2021.04.002
Shi, J. H., and Sun, S. C. (2018). Tumor necrosis factor receptor-associated factor regulation of nuclear factor κB and mitogen-activated protein kinase pathways. Front. Immunol. 9, 1849. doi:10.3389/fimmu.2018.01849
Simonneau, G., D'Armini, A. M., Ghofrani, H. A., Grimminger, F., Hoeper, M. M., Jansa, P., et al. (2015). Riociguat for the treatment of chronic thromboembolic pulmonary hypertension: A long-term extension study (CHEST-2). Eur. Respir. J. 45, 1293–1302. doi:10.1183/09031936.00087114
Simonneau, G., Gatzoulis, M. A., Adatia, I., Celermajer, D., Denton, C., Ghofrani, A., et al. (2013). Updated clinical classification of pulmonary hypertension. J. Am. Coll. Cardiol. 62, D34–D41. doi:10.1016/j.jacc.2013.10.029
Song, Z., Zhu, X., Jin, R., Wang, C., Yan, J., Zheng, Q., et al. (2014). Roles of the kinase TAK1 in CD40-mediated effects on vascular oxidative stress and neointima formation after vascular injury. PLoS One 9, e101671. doi:10.1371/journal.pone.0101671
Soon, E., Holmes, A. M., Treacy, C. M., Doughty, N. J., Southgate, L., Machado, R. D., et al. (2010). Elevated levels of inflammatory cytokines predict survival in idiopathic and familial pulmonary arterial hypertension. Circulation 122, 920–927. doi:10.1161/CIRCULATIONAHA.109.933762
Steiner, M. K., Syrkina, O. L., Kolliputi, N., Mark, E. J., Hales, C. A., and Waxman, A. B. (2009). Interleukin-6 overexpression induces pulmonary hypertension. Circ. Res. 104, 236–244, 28p following 244. 228p following 244. doi:10.1161/CIRCRESAHA.108.182014
Sutendra, G., Bonnet, S., Rochefort, G., Haromy, A., Folmes, K. D., Lopaschuk, G. D., et al. (2010). Fatty acid oxidation and malonyl-CoA decarboxylase in the vascular remodeling of pulmonary hypertension. Sci. Transl. Med. 2, 44ra58. doi:10.1126/scitranslmed.3001327
Team, R. S. (2016). RStudio: Integrated development environment. Boston, Massachusetts, USA: for R. RStudio.
Thenappan, T., Ormiston, M. L., Ryan, J. J., and Archer, S. L. (2018). Pulmonary arterial hypertension: Pathogenesis and clinical management. BMJ 360, j5492. doi:10.1136/bmj.j5492
Vinaixa, M., Samino, S., Saez, I., Duran, J., Guinovart, J. J., and Yanes, O. (2012). A guideline to univariate statistical analysis for LC/MS-Based untargeted metabolomics-derived data. Metabolites 2, 775–795. doi:10.3390/metabo2040775
Vonk Noordegraaf, A., Westerhof, B. E., and Westerhof, N. (2017). The relationship between the right ventricle and its load in pulmonary hypertension. J. Am. Coll. Cardiol. 69, 236–243. doi:10.1016/j.jacc.2016.10.047
Wang, J., Weigand, L., Wang, W., Sylvester, J. T., and Shimoda, L. A. (2005). Chronic hypoxia inhibits Kv channel gene expression in rat distal pulmonary artery. Am. J. Physiol. Lung Cell Mol. Physiol. 288, L1049–L1058. doi:10.1152/ajplung.00379.2004
Wishart, D. S., Jewison, T., Guo, A. C., Wilson, M., Knox, C., Liu, Y., et al. (2012). Hmdb 3.0—the human metabolome database in 2013. Nucleic Acids Res. 41, D801–D807. doi:10.1093/nar/gks1065
Xia, J., and MetPA, Wishart D. S. (2010). MetPA: A web-based metabolomics tool for pathway analysis and visualization. Bioinformatics 26, 2342–2344. doi:10.1093/bioinformatics/btq418
Xia, J., Psychogios, N., Young, N., and Wishart, D. S. (2009). MetaboAnalyst: A web server for metabolomic data analysis and interpretation. Nucleic Acids Res. 37, W652–W660. doi:10.1093/nar/gkp356
Xu, W., Janocha, A. J., and Erzurum, S. C. (2021). Metabolism in pulmonary hypertension. Annu. Rev. Physiol. 83, 551–576. doi:10.1146/annurev-physiol-031620-123956
Yung, L. M., Yang, P., Joshi, S., Augur, Z. M., Kim, S. S. J., Bocobo, G. A., et al. (2020). ACTRIIA-Fc rebalances activin/GDF versus BMP signaling in pulmonary hypertension. Sci. Transl. Med. 12, eaaz5660. doi:10.1126/scitranslmed.aaz5660
Zhang, D., Gaussin, V., Taffet, G. E., Belaguli, N. S., Yamada, M., Schwartz, R. J., et al. (2000). TAK1 is activated in the myocardium after pressure overload and is sufficient to provoke heart failure in transgenic mice. Nat. Med. 6, 556–563. doi:10.1038/75037
Zippel, N., Malik, R. A., Fromel, T., Popp, R., Bess, E., Strilic, B., et al. (2013). Transforming growth factor-β-activated kinase 1 regulates angiogenesis via AMP-activated protein kinase-α1 and redox balance in endothelial cells. Arterioscler. Thromb. Vasc. Biol. 33, 2792–2799. doi:10.1161/ATVBAHA.113.301848
Keywords: pulmonary hypertension, antiproliferative, metabolomics, combination therapy, right ventricle
Citation: Morales-Cano D, Izquierdo-García JL, Barreira B, Esquivel-Ruiz S, Callejo M, Pandolfi R, Villa-Valverde P, Rodríguez I, Cogolludo A, Ruiz-Cabello J, Perez-Vizcaino F and Moreno L (2023) Impact of a TAK-1 inhibitor as a single or as an add-on therapy to riociguat on the metabolic reprograming and pulmonary hypertension in the SUGEN5416/hypoxia rat model. Front. Pharmacol. 14:1021535. doi: 10.3389/fphar.2023.1021535
Received: 17 August 2022; Accepted: 13 March 2023;
Published: 29 March 2023.
Edited by:
Alessia Vignoli, University of Florence, ItalyReviewed by:
Valentina Biasin, Ludwig Boltzmann Institute for Lung Vascular Research, AustriaFrancesca Di Cesare, University of Florence, Italy
Copyright © 2023 Morales-Cano, Izquierdo-García, Barreira, Esquivel-Ruiz, Callejo, Pandolfi, Villa-Valverde, Rodríguez, Cogolludo, Ruiz-Cabello, Perez-Vizcaino and Moreno. 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: Laura Moreno, lmorenog@med.ucm.es
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share last authorship