Network Pharmacology and Molecular Docking Analysis of Shufeiya Recipe in the Treatment of Pulmonary Hypertension

Zhang, Zeyu; Wang, Xianliang; Wang, Shuai; Jia, Zhuangzhuang; Mao, Jingyuan

doi:https://doi.org/10.1155/2022/7864976

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 7864976 | https://doi.org/10.1155/2022/7864976

Network Pharmacology and Molecular Docking Analysis of Shufeiya Recipe in the Treatment of Pulmonary Hypertension

Zeyu Zhang,¹Xianliang Wang,¹Shuai Wang,¹Zhuangzhuang Jia,¹and Jingyuan Mao¹

Academic Editor: Nadeem Sheikh

Received25 Jul 2022

Revised16 Nov 2022

Accepted05 Dec 2022

Published28 Dec 2022

Abstract

Objective. This study is aimed at exploring the molecular mechanism of Shufeiya recipe in the treatment of pulmonary hypertension (PH) using network pharmacology and molecular docking analysis. Methods. Active components and their target proteins in the recipe were screened using the TCMSP database. PH-related core proteins were screened using GeneCards, STRING database, and Cytoscape-v3.8.2. Common proteins were obtained by intersection of the target proteins of these recipe active components and pH-related core proteins. Rx64 4.0.2 software was used to perform GO functional enrichment analysis and KEGG pathway enrichment analysis on the common proteins to obtain pathway-enriched proteins, and then core enriched proteins were further screened. We analyzed the relationship between the active components and pathway-enriched proteins using Cytoscape-v3.8.2. AutoDock Vina was used to dock their core proteins into the components. Results. Shufeiya recipe contained 67 active components. 61 common proteins of the target proteins of the active components and PH-related core proteins were obtained. The treatment involved both functional and pathway regulations. The core pathway-enriched proteins were prostaglandin G/H synthase 2 (PTGS2), epidermal growth factor receptor (EGFR), and RAC-alpha serine/threonine-protein kinase (AKT1), and their binding energies to the corresponding components were all less than -5 kJ•mol-1. Conclusion. It was found that the main mechanism might be the active components acting on the core pathway-enriched proteins to regulate related signaling pathways, thereby playing roles in anticoagulation, vasodilation, anti-PASMC proliferation, promotion of PAECs apoptosis, inhibition of oxidative stress, and anti-inflammatory effects.

1. Introduction

Pulmonary hypertension (PH) is a chronic progressive pulmonary vascular system disease caused by pulmonary arteriole obstruction, and it is characterized by pulmonary artery vasoconstriction and remodeling, increased pulmonary circulatory resistance, and increased pulmonary artery pressure ( mmHg). Its clinical manifestations include chest tightness, shortness of breath, fatigue, and dyspnea, which can eventually lead to right-sided heart failure and death [1]. Epidemiological surveys show that about 1% of people suffer from PH around the world, about 10% of people over 65 years old suffer from PH, the 1-year mortality rate of this disease is about 15%, and the 3-year mortality rate is about 30% [2, 3]. Therefore, it is important to seek effective methods for the prevention and treatment of PH.

The clinical practice of PH treatment in our country not only fully applies modern research results but also recognizes that traditional Chinese medicine (TCM) has the synergistic effects of improving clinical symptoms and quality of life, increasing exercise tolerance, and even improving prognosis. The combination of Chinese and western medicine has become a generally accepted treatment strategy for PH in our country. Based on his clinical practice and researches of many years, Professor Mao Jingyuan, as a “Qihuang Scholar,” and national leading talent of TCM, concludes that Shufeiya recipe can be used for the treatment of PH; this recipe is composed of platycodi radix (Jiegeng), Salvia miltiorrhiza (Danshen), safflower (Honghua), and Cornus officinalis (Shanzhuyu). Previous animal studies suggested that Shufeiya recipe could downregulate serum reactive oxygen species levels in rats with pulmonary hypertension induced by monocrotaline; upregulate serum NO, sGC, cGMP, PKG, PTGS1, PTGS2, and Mn-SOD activity; inhibit oxidative damage; promote the expression of SIRT3, FOXO3a, p-PI3K, p-AKT, and p-eNOS protein in the pulmonary artery; downregulate the protein expression of Ras, p-MEK1/2, p-ERK1/2, and c-fos in the pulmonary artery; and reduce the thickening of the pulmonary vascular wall and pulmonary artery pathological injury such as cavity stenosis, thereby reducing pulmonary artery pressure in the model group [4, 5]. However, the specific mechanism of action of the Shufeiya recipe has not been fully elucidated, and needs to be explored using more analytical methods and from various perspectives. Therefore, this study is aimed at exploring the molecular mechanism of Shufeiya recipe in the treatment of pulmonary hypertension (PH) using network pharmacology and molecular docking analysis.

2. Materials and Methods

2.1. Active Components Screening and their Target Proteins Obtained

Based on the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) (http://tcmspw.com/tcmsp.php) [6] and the screening criteria of , , , , , and relative molecular weight (MW) <500, the active components and their target proteins of platycodi radix, Salvia miltiorrhiza, safflower, and Cornus officinalis were screened, and the target proteins were imported into the UniProt database (https://www.uniprot.org/) [7] to standardize the proteins.

2.2. The “Drug-Disease” Targets Obtained

Based on the GeneCards database (http://www.GeneCards.org/) [8], PH-related proteins were obtained, with PH as the search term and a relevance score of >6 as the screening criteria. The STRING database (http://string-db.org) [9] was used to predict protein-protein interactions with the criterion of . The Cytoscape-v3.8.2 software [10] was used to screen PH-related core proteins, with the condition that degree, betweenness, closeness, network, and local marginal connectivity (LAC) were all greater than the median value. Finally, the intersection of “PH-related core proteins” and “target proteins of the active components” was used to obtain the common targets, named “drug-disease” targets.

2.3. Gene Ontology (GO) Functional Enrichment Analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Enrichment Analysis

The “drug-disease” targets were subjected to GO functional enrichment analysis, and KEGG pathway enrichment analysis using Rx64 4.0.2 software [11], with as the screening criteria. GO functional enrichment analysis mainly included three aspects: biological processes (BP), molecular functions (MF), and cellular components (CC). The top ten functions obtained in each aspect by GO enrichment analysis were visualized using a bar chart. The top 30 signaling pathways obtained by KEGG pathway enrichment analysis were visualized using a bubble chart, where the ordinate represents the pathway name, the abscissa represents the gene ratio, the bubble size represents the number of enriched proteins, the color represents the degree of enrichment, and the darker the red, the stronger the enrichment. PH-related pathways were screened by referring to the literature, and enriched proteins were obtained from KEGG pathway enrichment analysis. The map of the pathways most significantly related to the treatment of PH by Shufeiya recipe was obtained from the KEGG database (https://www.kegg.jp/) [12]. The red nodes in this map represent the pathway-enriched proteins.

2.4. Screening of Core Enriched Proteins, and Construction of the Pathway-Target Topology Network

The enriched proteins of PH-related pathways were imported into the STRING database, and Cytoscape-v3.8.2 software was used to screen the core enriched proteins of PH-related pathways with the condition that the degree, betweenness, closeness, network, and LAC were all greater than the median value. Then, PH-related pathways and these enriched proteins were imported into Cytoscape-v3.8.2 software for topology analysis and visualization, and a pathway-target topology network was constructed.

2.5. Construction of the Component-Target Topology Network

The active components of Shufeiya recipe and pathway-enriched proteins were imported into Cytoscape-v3.8.2 software for topology analysis and visualization, and the component-target topology network was constructed, where nodes represented components and targets, and edges represented interactions between nodes.

2.6. Molecular Docking Verification and Screening

Molecular docking was performed with the core-enriched proteins as the receptor, and the corresponding drug components as the ligand. The 2D structure of the drug component was downloaded from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) [13] and imported into the Chem3D software to be converted into a 3D structure. Then, we selected the Minimize Energy of the MM2 item in the calculations to optimize the component structure, and saved it as lig.mol2. The protein entry was obtained from the UniProt database (https://www.uniprot.org/), and entered into the PDB database (http://www.rcsb.org/) [14] to obtain the protein crystal structure; then the protein crystal structure was imported into PyMOL 2.5 software [15], “remove solvent and remove organic” item was used to remove its H2O and small molecule ligands, and the file was saved as rep.pdb. We imported rep.pdb into AutoDockTools 1.5.6 software [16], added Hydrogens, and saved it in the PDB format. We selected Macromolecule of the Grid item to output the protein in rep.pdbqt format; we clicked ligand to import lig.mol2 and output it as lig.pdbqt format. We selected the secondary structure and undisplay lines of the display item to convert the protein and the component into the secondary structure, and then selected the Grid BOX to determine the binding site of protein and component using the number of points in , , , as conditions, and the output file was saved as grid.gpf. We saved rep.pdbqt, lig.pdbqt, and grid.gpf in the same folder, and then vina.exe, vina_split.exe, config.txt (, , center_x =, center_y =, center_z =, , , , , and ) were saved in this folder, and center_x, center_y, and center_z of config.txt were filled with the gridcenter data of the grid.gpf file. The path of the folder was entered into Perl, and then the command “vina.exe --config config.txt --log log.txt --out output.pdbqt” was entered to calculate affinity. Twenty conformations were generated for docking results. The conformation with the best affinity was selected as the final docking conformation and was visualized using PyMOL 2.5.

3. Results

3.1. Active Components of Shufeiya Recipe and their Target Proteins

A total of 67 active components of Shufeiya recipe were obtained from the TCMSP database, including 3 active components from platycodi radix, 52 active components from Salvia miltiorrhiza, 10 active components from safflower, and 4 active components from Cornus officinalis (Supplementary Table 1). After protein prediction and standardization, 214 target proteins were identified.

3.2. PH-Related Core Proteins Screening and the “Drug-Disease” Targets Obtained

The PH-related core protein screening process was visualized using GeneCards, STRING database and Cytoscape-v3.8.2 software (Figures 1(a)–1(c)). The red nodes in figure 1(a)–1(c) represents the proteins whose degree, betweenness, closeness, network, and LAC were all less than the median value, and the yellow nodes represent the proteins whose degree, betweenness, closeness, network, and LAC were all greater than the median value. Figure 1(a) shows the first screening result, Figure 1(b) shows the second screening result, and Figure 1(c) shows the PH-related core proteins. A total of 210 PH-related core proteins were obtained. After taking the intersection of 210 PH-related core proteins and 214 target proteins of drug active components, 61 common proteins (“drug-disease” targets) were obtained (Figure 2(a), Table 1).

(a)

(b)

(c)

(a)

(b)

(c)

3.3. GO and KEGG Enrichment Analyses

A total of 2469 functional terms were obtained by GO functional enrichment analysis, including 2337 terms for BP, 299 terms for CC, and 100 terms for MF. The top ten functional terms for BP, MF and CC are shown in Figure 2(b). BP mainly included response to oxidative stress, cellular response to chemical stress, and epithelial cell proliferation; CC mainly included membrane raft membrane rafts, membrane microdomain, and membrane region; MF mainly included cytokine receptor binding, DNA-binding transcription factor binding, and signaling receptor activator activity. A total of 161 signaling pathways were obtained by the KEGG pathway enrichment analysis. The top 30 signaling pathways are shown in Figure 2(c). Previous studies have reported that the PI3K-AKT signaling pathway [17], IL-17 signaling pathway [18], TNF pathway [19], MAPK signaling pathway [20], HIF-1 signaling pathway [21], and Toll-like receptor signaling pathway [22] were related to the occurrence of PH, among which the PI3K-AKT pathway was most closely related to the treatment of PH with Shufeiya recipe (Figure 3).

3.4. Core-Enriched Proteins Screening and Pathway-Target Topology Network Analysis

KEGG enrichment analysis showed that 46 proteins were enriched in the signaling pathways related to the treatment of pH by the Shufeiya recipe (Table 2). Based on the STRING database and Cytoscape-v3.8.2 software, these signaling pathways and enriched proteins were utilized to construct the pathway-target topology network for screening 19 core-enriched proteins (Figure 4(a)). In Figure 4(a), the blue arrows represent the signaling pathways related to the treatment of PH by Shufeiya recipe, the red nodes represent general enriched proteins, and the yellow nodes represent core enriched proteins, suggesting that prostaglandin G/H synthase 2 (PTGS2), epidermal growth factor receptor (EGFR), RAC-alpha serine/threonine-protein kinase (AKT1), vascular endothelial growth factor A (VEGFA), mitogen-activated protein kinase 1 (MAPK1), interleukin-6 (IL6), and receptor tyrosine-protein kinase erbB-2 (ERBB2) might play key roles in regulating these signaling pathways. Degree value refers to the number of pathway-enriched proteins. The PI3K-Akt signaling pathway was 22, the IL-17 signaling pathway was 20, the TNF signaling pathway was 20, the MAPK signaling pathway was 20, the HIF-1 signaling pathway was 16, and the Toll-like receptor signaling pathway was 15, indicating that the PI3K-Akt signaling pathway was the most regulated protein by Shufeiya recipe.

(a)

(b)

3.5. Component-Target Topology Network Analysis

The component-target topology network analysis showed that the network contained 110 nodes, including 64 nodes for drug components, 46 nodes for target proteins, and 167 edges. According to the data obtained by topological analysis, the components of Shufeiya recipe with are listed in Table 3, which indicates that the main components of this recipe for the treatment of PH were quercetin, luteolin, baicalein, kaempferol, tanshinone IIA, and cryptotanshinone. The component-target topology network is shown in Figure 4(b). The left circle in Figure 4(b) represents the active components of this recipe, where green nodes represented components from platycodi radix, purple nodes represent components from Salvia miltiorrhiza, blue nodes represent components from safflower, orange nodes represent components from Cornus officinalis, and nodes of multiple colors represent the common components of multiple drugs. The right circle in figure 4(b) represents proteins enriched in the signaling pathways related to the treatment of PH by Shufeiya recipe, where the red nodes represent the general enriched proteins, and the yellow nodes represent the core enriched proteins. The degree value refers to the number of active components of the regulatory proteins. Among the core enriched proteins, PTGS2 was 63, CASP3 was 6, TP53 was 5, AKT1 was 4, TNF was 4, JUN was 4, MMP9 was 4, FOS was 3, and ICAM1 was 3, VEGFA was 3, NFKBIA was 3, EGFR was 2, IL6 was 2, HIF1A was 2, MYC was 2, CCL2 was 1, IL1B was 1, STAT3 was 1, and EGF was 1, which indicated the core protein PTGS2 was most strongly regulated by the components of this recipe.

3.6. Molecular Docking of these Core Enriched Proteins and the Drug Active Components

When the binding energy between the receptor and ligand is lower than -5 kJ·mol^-1, it indicates that the active component can bind well to the target protein [23]. The docking accuracy of the core enriched protein and its corresponding components was verified through molecular docking. The results showed that the binding energies of the core-enriched proteins and their corresponding components were all less than -5 kJ•mol^-1 (Table 4). As shown in figure 5, a molecular docking model was constructed for the component and its target protein with the strongest affinity. Figure 5(a) shows the molecular docking model of quercetin and MMP9, Figure 5(b) shows the molecular docking model of luteolin and MMP9, Figure 5(c) shows the molecular docking model of baicalein and MMP9, Figure 5(d) shows the molecular docking model of kaoneferol and JUN, Figure 5(e) shows the molecular docking model of tanshinone IIA and PTGS2, and Figure 5(f) shows the molecular docking model of cryptotanshinone and PTGS2. This indicated that there was good binding ability between the drug components and their targets.

(a)

(b)

(c)

(d)

(e)

(f)

4. Discussion

The occurrence of PH is closely related to abnormalities in the pulmonary vascular function and/or structure. Long-term hypoxia or other factors damage the endothelial cell structure, resulting in an imbalance of vasoactive substances, such as NO, ion channel dysfunction, and abnormal endothelial cell barrier function, which eventually leads to abnormal contraction of pulmonary blood vessels, progressive stenosis, and occlusion of the lumen. Its pathological characteristics include pulmonary artery intima damage and proliferation, media hypertrophy, adventitial matrix deposition fibrosis and inflammatory infiltration, distal vasodilation, and in situ thrombus [1]. Traditional Chinese medicine (TCM) is widely used for the treatment of PH because of its unique efficacy and safety. Modern science has thoroughly studied the mechanism of action of the components of traditional Chinese medicine to promote their application, so that Chinese and Western medicines complement each other. After years of clinical practice, traditional Chinese medicine believes that “stasis and bloking of lung and heart” is the pathogenesis of PH, Shufeiya recipe has functions of dispersing lung and removing stasis, and nourishing yin and dredging pulse, and the clinical effects of this recipe are quite good. Network pharmacology can clarify the mechanism of action of drugs at multiple levels, and molecular docking can be used to calculate the affinity between receptors and ligands. Therefore, the present study used network pharmacology and molecular docking methods to explore the relationship between Shufeiya recipe and PH, and preliminarily explored the specific mechanisms of functions, signaling pathways, components, and targets involved in the treatment of PH using Shufeiya recipe.

GO functional enrichment analysis showed that Shufeiya recipe for treating PH involved the regulation of functions such as oxidative stress, chemical stress, epithelial cell proliferation, and signaling receptor agonist activity. Oxidative stress is one of the key causes of PH, and is a process of oxidative damage caused by the imbalance between the generation and metabolism of oxygen free radicals in tissue cells, resulting in the accumulation of reactive oxygen species (ROS) and an imbalance in the ratio of antioxidant enzymes to lipid peroxides [24, 25]. ROS can activate or inhibit a series of signaling pathways, resulting in pathological changes, such as pulmonary artery smooth muscle cells (PASMCs) proliferation and migration, pulmonary artery endothelial cells (PAECs) disorders, situ thrombosis, inflammation, pulmonary vascular remodeling, vasoconstriction response, and extracellular matrix deposition [24, 25]. Our previous study confirmed that Shufeiya recipe could downregulate the serum level of ROS, upregulate PTGS1 and PTGS2 concentrations, and Mn-SOD activity, and inhibit oxidative damage [5]. KEGG pathway enrichment analysis showed that Shufeiya recipe could inhibit the development of PH by regulating PI3K-AKT, IL-17, TNF, MAPK, HIF-1, and Toll-like receptor signaling pathways. Activation of the PI3K-AKT signaling pathway can inhibit hypoxia-induced PASMCs cycle progression, proliferation, migration, and autophagy; regulate PASMCs cytoskeletal rearrangement, and phenotypic transition; protect PAECs from hypoxia-induced apoptosis; inhibit hypoxia-induced pulmonary vascular remodeling; and promotes angiogenesis and survival, thereby attenuating pulmonary vascular remodeling, and the development of PH [17]. Our previous study confirmed that Shufeiya recipe can promote the protein expression of p-PI3K and p-AKT [5]. Activation of the IL-17 signaling pathway can induce vascular oxidase, endothelial nitric oxide synthase, vascular smooth muscle cells, and adventitial fibroblast dysfunction [18]. Overexpression of TNF-α can promote the inflammatory response in rats with PH and reduce the expression of vascular endothelial growth factor and flk-1, causing pulmonary vascular remodeling [19]. The MAPK signaling pathway mainly includes ERK1/2, p38MAPK, and JNK, and hypoxia-induced activation of Ca²⁺-dependent PKC/MAPK signaling pathway in PH promotes PASMCs proliferation [26]. PH caused by left heart disease can promote the expression of matrix metalloproteinase 9 and transforming growth factor β-1 proprotein by activating the MAPK signaling pathway, which is involved in pulmonary artery vascular remodeling [27]. p38 MAPKα-related protein inhibition can effectively reduce inflammation, increase NO production, reduce superoxide load, inhibit pulmonary artery fibroblast proliferation, restore PAECs function, and reverse pulmonary vascular remodeling [28]. ERK activation can promote the proliferation of pulmonary artery fibroblasts, while Shufeiya recipe can downregulate the expression levels of Ras, p-MEK1/2, p-ERK1/2, and c-fos proteins in the ERK signaling pathway, thereby inhibitting PAECs fibrosis [5]. Activation of the HIF-1 signaling pathway in PAECs can promote the proliferation of pulmonary artery medial cells, and the expression of factors such as endothelin 1 [21, 29]. The Toll-like receptor signaling pathway is closely related to the proliferation, apoptosis, migration, and interstitial transition of PAECs [22, 30]. For example, inhibition of Toll-like receptor 3 can promote caspase-dependent PAECs apoptosis, and inhibition of Toll-like receptor 4 can limit PASMCs proliferation and migration [22, 30]. Screening of core-enriched proteins, pathway-target topology analysis, component-target topology analysis, and molecular docking suggested that quercetin, luteolin, baicalein, kaempferol, tanshinone IIA, and cryptotanshinone, as the main components of Shufeiya recipe, might effectively intervene in core proteins, including PTGS2, EGFR, AKT1, VEGFA, and MAPK1, to regulate their signaling pathways. Modern pharmacology suggests that quercetin can improve pulmonary arterial pressure, right ventricular hypertrophy, and vascular remodeling through anticoagulation, vasodilator, antiproliferation, and induction of PASMCs apoptosis [31]; luteolin can downregulate platelet-derived growth factor-BB-induced AKT phosphorylation, thereby inhibiting the proliferation and migration of PASMCs [17]; baicalein can block the production of 12(S)-hydroxyeicosatetraenoic acid to resist ERK1/2 phosphorylation, and inhibit hypoxia-induced proliferation of PASMCs [32]; kaempferol and its derivatives have a wide range of pharmacological activities, including anti-inflammatory, antibacterial, antioxidant, and relaxation of pulmonary arteries in an endothelium-independent manner [33]; tanshinone IIA can activate the PI3K/Akt/mTOR signaling pathway, reduce oxidative stress, inhibit inflammatory response, and downregulate the expression of high mobility group box 1 protein, thereby inhibiting the proliferation of PAECs and inflammatory response; in addition, its water-soluble derivative of tanshinone IIA sodium sulfonate can relax blood vessels on this basis [34–36]; and cryptotanshinone has anti-inflammatory effects, and inhibits epithelial-mesenchymal transition [37, 38]. Therefore, the Shufeiya recipe for the treatment of PH might have multifunction, multichannel, multicomponent, and multitarget advantages.

However, there are some limitations in this study. Our results were derived only from the analysis and prediction of the public database and silicon computers, and predicting herb targets only used one platform (TCMSP). Therefore, we will explore more detailed information about the recipe by searching a few more platform, such as BATMAN-TCM, Herb, and TCMID, and in vivo and in vitro experiments need to be performed for further validation in the future.

In summary, the molecular mechanism of Shufeiya recipe in the treatment of PH was analyzed using network pharmacology. We screened the active components and their targets of Shufeiya recipe in the treatment of PH, constructed component-target networks, analyzed their biological process annotations, and signaling pathways, and systematically studied the potential mechanism of Shufeiya recipe in treating PH. It was found that the main mechanism of this recipe in the treatment of PH might be that the active components of this recipe, including quercetin, luteolin, baicalein, kaempferol, tanshinone IIA, and cryptotanshinone, acted on the core proteins, including PTGS2, EGFR, AKT1, VEGFA, and MAPK1, to regulate signaling pathways such as PI3K-AKT, IL-17, TNF, MAPK, HIF-1, and Toll-like receptors, thereby playing the roles of anticoagulation, vasodilator, anti-PASMC proliferation, promotion of PAEC apoptosis, inhibition of oxidative stress, and anti-inflammatory effects. Our present study identified the core targets of Shufeiya recipe in the treatment of PH at the molecular level, clarified its potential material basis and mechanism, and provided a general reference for the clinical application and experimental research of traditional Chinese medicine in the treatment of PH. However, these results need to be confirmed by more experiments and clinical practices in the future.

Data Availability

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Conflicts of Interest

All authors declare that they have no conflict of interests.

Authors’ Contributions

(1) Zeyu Zhang and Jingyuan Mao provided the substantial contributions to the conception and design of the work; (2) Zeyu Zhang, Xianliang Wang, Shuai Wang, and Zhuangzhuang Jia were responsible for the acquisition, analysis, and interpretation of data for the work; (3) Zeyu Zhang was assigned in drafting the work; (4) Jingyuan Mao was tasked in revising it critically for important intellectual content; (5) Zeyu Zhang, Xianliang Wang, Shuai Wang, Zhuangzhuang Jia, Jingyuan Mao were assigned in final approval of the version to be published; and (6) Zeyu Zhang, Xianliang Wang, Shuai Wang, Zhuangzhuang Jia, and Jingyuan Mao were responsible for the agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Acknowledgments

Innovation Team Development Program of the Ministry of Education, No. IRT-16R54; Research Project of Integrated Traditional Chinese and Western Medicine of Tianjin Health Commission and Tianjin Administration of Traditional Chinese Medicine, No. 2021048.

Supplementary Materials

Supplementary Table 1. Active components of Shufeiya recipe. (Supplementary Materials)

References

D. Poch and J. Mandel, “Pulmonary Hypertension,” Annals of Internal Medicine, vol. 174, no. 4, pp. ITC49–ITC64, 2021.
View at: Publisher Site | Google Scholar
M. M. Hoeper, M. Humbert, R. Souza et al., “A global view of pulmonary hypertension,” The Lancet Respiratory Medicine, vol. 4, no. 4, pp. 306–322, 2016.
View at: Publisher Site | Google Scholar
E. Spiekerkoetter, S. M. Kawut, and V. A. de Jesus Perez, “New and emerging therapies for pulmonary arterial hypertension,” Annual Review of Medicine, vol. 70, p. 45, 2019.
View at: Publisher Site | Google Scholar
H. Yan, L. Wang, X. Wang et al., “Effect of Shufeiya formula on pulmonary arterial pressure and oxidative stress in monocrotaline-induced pulmonary hypertension model rats,” Journal of Traditional Chinese Medicine, vol. 62, no. 20, pp. 1814–1819, 2021.
View at: Publisher Site | Google Scholar
Z. Jia, H. Yan, S. Wang et al., “Shufeiya recipe improves monocrotaline-induced pulmonary hypertension in rats by regulating SIRT3/FOXO3a and its downstream signaling pathways,” Disease Markers, vol. 2022, Article ID 3229888, 2022.
View at: Publisher Site | Google Scholar
J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, no. 1, 2014.
View at: Publisher Site | Google Scholar
“UniProt: the universal protein knowledgebase in 2021,” Nucleic Acids Research, vol. 49, no. D1, pp. D480–d489, 2021.
View at: Publisher Site | Google Scholar
R. Barshir, S. Fishilevich, T. Iny-Stein et al., “GeneCaRNA: a comprehensive gene-centric database of human non-coding RNAs in the GeneCards suite,” Journal of Molecular Biology, vol. 433, no. 11, article 166913, 2021.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, K. C. Nastou et al., “The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets,” Nucleic Acids Research, vol. 49, no. D1, pp. D605–d612, 2021.
View at: Publisher Site | Google Scholar
D. Otasek, J. H. Morris, J. Bouças, A. R. Pico, and B. Demchak, “Cytoscape automation: empowering workflow-based network analysis,” Genome Biology, vol. 20, no. 1, p. 185, 2019.
View at: Publisher Site | Google Scholar
B. K. C. Chan, “Data Analysis Using R Programming,” in Biostatistics for Human Genetic Epidemiology. Advances in Experimental Medicine and Biology, vol. 1082, Springer, Cham, 2018.
View at: Publisher Site | Google Scholar
M. Kanehisa, M. Furumichi, M. Tanabe, Y. Sato, and K. Morishima, “KEGG: new perspectives on genomes, pathways, diseases and drugs,” Nucleic Acids Research, vol. 45, no. D1, pp. D353–d361, 2017.
View at: Publisher Site | Google Scholar
S. Kim, J. Chen, T. Cheng et al., “PubChem in 2021: new data content and improved web interfaces,” Nucleic Acids Research, vol. 49, no. D1, pp. D1388–d1395, 2021.
View at: Publisher Site | Google Scholar
S. K. Burley, H. M. Berman, G. J. Kleywegt, J. L. Markley, H. Nakamura, and S. Velankar, “Protein Data Bank (PDB): The Single Global Macromolecular Structure Archive,” in Protein Crystallography. Methods in Molecular Biology, A. Wlodawer, Z. Dauter, and M. Jaskolski, Eds., vol. 1607, Humana Press, New York, NY, 2017.
View at: Publisher Site | Google Scholar
S. Yuan, H. C. S. Chan, S. Filipek, and H. Vogel, “PyMOL and Inkscape Bridge the Data and the Data Visualization,” Structure, vol. 24, no. 12, pp. 2041-2042, 2016.
View at: Publisher Site | Google Scholar
N. El-Hachem, B. Haibe-Kains, A. Khalil, F. H. Kobeissy, and G. Nemer, “AutoDock and AutoDockTools for Protein-Ligand Docking: Beta-Site Amyloid Precursor Protein Cleaving Enzyme 1(BACE1) as a Case Study,” in Neuroproteomics. Methods in Molecular Biology, F. Kobeissy and S. Stevens Jr., Eds., vol. 1598, Humana Press, New York, NY, 2017.
View at: Publisher Site | Google Scholar
W. Zuo, N. Liu, Y. Zeng et al., “Luteolin ameliorates experimental pulmonary arterial hypertension via suppressing hippo-YAP/PI3K/AKT signaling pathway,” Frontiers in Pharmacology, vol. 12, 2021.
View at: Publisher Site | Google Scholar
S. S. Batah, M. A. Alda, R. R. Figueira et al., “In situ evidence of collagen V and Interleukin-6/Interleukin-17 activation in vascular remodeling of experimental pulmonary hypertension,” Pathobiology, vol. 87, no. 6, pp. 356–366, 2020.
View at: Publisher Site | Google Scholar
M. Fujita, R. J. Mason, C. Cool, J. M. Shannon, N. Hara, and K. A. Fagan, “Pulmonary hypertension in TNF-alpha-overexpressing mice is associated with decreased VEGF gene expression,” Journal of Applied Physiology, vol. 93, no. 6, pp. 2162–2170, 2002.
View at: Publisher Site | Google Scholar
A. C. Church, D. H. Martin, R. Wadsworth et al., “The reversal of pulmonary vascular remodeling through inhibition of p38 MAPK-alpha: a potential novel anti-inflammatory strategy in pulmonary hypertension,” American Journal of Physiology - Lung Cellular and Molecular Physiology, vol. 309, no. 4, pp. L333–L347, 2015.
View at: Publisher Site | Google Scholar
D. Fu, A. Dai, R. Hu, Y. Chen, and L. Zhu, “Expression and role of factor inhibiting hypoxia-inducible factor-1 in pulmonary arteries of rat with hypoxia-induced hypertension,” Acta Biochimica et Biophysica Sinica, vol. 40, no. 10, pp. 883–892, 2008.
View at: Publisher Site | Google Scholar
D. Farkas, A. A. R. Thompson, A. R. Bhagwani et al., “Toll-like receptor 3 is a therapeutic target for pulmonary hypertension,” American Journal of Respiratory and Critical Care Medicine, vol. 199, no. 2, pp. 199–210, 2019.
View at: Publisher Site | Google Scholar
K. Y. Hsin, S. Ghosh, and H. Kitano, “Combining machine learning systems and multiple docking simulation packages to improve docking prediction reliability for network pharmacology,” PLoS One, vol. 8, no. 12, article e83922, 2013.
View at: Publisher Site | Google Scholar
A. C. Racanelli, S. A. Kikkers, A. M. K. Choi, and S. M. Cloonan, “Autophagy and inflammation in chronic respiratory disease,” Autophagy, vol. 14, no. 2, pp. 221–232, 2018.
View at: Publisher Site | Google Scholar
N. Ghasemzadeh, R. S. Patel, D. J. Eapen et al., “Oxidative stress is associated with increased pulmonary artery systolic pressure in humans,” Hypertension, vol. 63, no. 6, pp. 1270–1275, 2014.
View at: Publisher Site | Google Scholar
H. Tian, L. Liu, Y. Wu et al., “Resistin-like molecule β acts as a mitogenic factor in hypoxic pulmonary hypertension via the ca(2+)-dependent PI3K/Akt/mTOR and PKC/MAPK signaling pathways,” Respiratory Research, vol. 22, no. 1, p. 8, 2021.
View at: Publisher Site | Google Scholar
H. Zhang, W. Huang, H. Liu, Y. Zheng, and L. Liao, “Mechanical stretching of pulmonary vein stimulates matrix metalloproteinase-9 and transforming growth factor-β1 through stretch-activated channel/MAPK pathways in pulmonary hypertension due to left heart disease model rats,” PLoS One, vol. 15, no. 9, article e0235824, 2020.
View at: Publisher Site | Google Scholar
R. P. Weerackody, D. J. Welsh, R. M. Wadsworth, and A. J. Peacock, “Inhibition of p38 MAPK reverses hypoxia-induced pulmonary artery endothelial dysfunction,” American Journal of Physiology. Heart and Circulatory Physiology, vol. 296, no. 5, pp. H1312–H1320, 2009.
View at: Publisher Site | Google Scholar
G. B. Waypa and P. T. Schumacker, “Roles of HIF1 and HIF2 in pulmonary hypertension: it all depends on the context,” The European Respiratory Journal, vol. 54, no. 6, article 1901929, 2019.
View at: Publisher Site | Google Scholar
A. Bhagwani, A. A. R. Thompson, and L. Farkas, “When innate immunity meets angiogenesis-the role of toll-like receptors in endothelial cells and pulmonary hypertension,” Frontiers in Medicine, vol. 7, 2020.
View at: Publisher Site | Google Scholar
D. Morales-Cano, C. Menendez, E. Moreno et al., “The flavonoid quercetin reverses pulmonary hypertension in rats,” PLoS One, vol. 9, no. 12, article e114492, 2014.
View at: Publisher Site | Google Scholar
I. R. Preston, N. S. Hill, R. R. Warburton, and B. L. Fanburg, “Role of 12-lipoxygenase in hypoxia-induced rat pulmonary artery smooth muscle cell proliferation,” American Journal of Physiology - Lung Cellular and Molecular Physiology, vol. 290, no. 2, pp. L367–L374, 2006.
View at: Publisher Site | Google Scholar
A. Mahobiya, T. U. Singh, S. Rungsung et al., “Kaempferol-induces vasorelaxation _via_ endothelium-independent pathways in rat isolated pulmonary artery,” Pharmacological Reports, vol. 70, no. 5, pp. 863–874, 2018.
View at: Publisher Site | Google Scholar
J. Wang, Q. Jiang, L. Wan et al., “Sodium tanshinone IIA sulfonate inhibits canonical transient receptor potential expression in pulmonary arterial smooth muscle from pulmonary hypertensive rats,” American Journal of Respiratory Cell and Molecular Biology, vol. 48, no. 1, pp. 125–134, 2013.
View at: Publisher Site | Google Scholar
J. Wang, M. Q. Dong, M. L. Liu et al., “Tanshinone IIA modulates pulmonary vascular response to agonist and hypoxia primarily via inhibiting Ca²⁺ influx and release in normal and hypoxic pulmonary hypertension rats,” European Journal of Pharmacology, vol. 640, no. 1-3, pp. 129–138, 2010.
View at: Publisher Site | Google Scholar
Q. Jiang, W. Lu, K. Yang et al., “Sodium tanshinone IIA sulfonate inhibits hypoxia-induced enhancement of SOCE in pulmonary arterial smooth muscle cells via the PKG-PPAR-γ signaling axis,” American Journal of Physiology. Cell Physiology, vol. 311, no. 1, pp. C136–C149, 2016.
View at: Publisher Site | Google Scholar
X. X. Li, X. Zheng, Z. Liu et al., “Cryptotanshinone from salvia miltiorrhiza bunge (Danshen) inhibited inflammatory responses via TLR4/MyD88 signaling pathway,” Chinese Medicine, vol. 15, p. 20, 2020.
View at: Publisher Site | Google Scholar
L. Avila-Carrasco, P. Majano, J. A. Sánchez-Toméro et al., “Natural plants compounds as modulators of epithelial-to-mesenchymal transition,” Frontiers in Pharmacology, vol. 10, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Zeyu Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

435

Downloads

645

Citations