- 1Institute of Basic Research in Clinical Medicine, China Academy of Chinese Medical Sciences, Beijing, China
- 2Guang’anmen Hospital, China Academy of Chinese Medical Sciences, Beijing, China
- 3Research Department, Swiss University of Traditional Chinese Medicine, Bad Zurzach, Switzerland
Background: With the increasing global prevalence of hypertension, a condition that can severely affect multiple organs, there is a growing need for effective treatment options. Uncaria rhynchophylla-Alisma plantago-aquatica L. (UR-AP) is a traditional drug pair used for treating hypertension based on the liver-kidney synergy concept. However, the detailed molecular mechanisms underlying its efficacy remain unclear.
Methods: This study utilized an integrative approach combining network pharmacology, cluster analysis, and molecular docking to uncover the bioactive components and targets of UR-AP in the treatment of hypertension. Initially, we extracted data from public databases to identify these components and targets. A Protein-Protein Interaction (PPI) network was constructed, followed by enrichment analysis to pinpoint the bioactive components, core targets, and pivotal pathways. Cluster analysis helped in identifying key sub-networks and hypothesizing primary targets. Furthermore, molecular docking was conducted to validate the interaction between the core targets and major bioactive components, thus confirming their potential efficacy in hypertension treatment.
Results: Network pharmacological analysis identified 58 bioactive compounds in UR-AP, notably quercetin, kaempferol, beta-sitosterol (from Uncaria rhynchophylla), and Alisol B, alisol B 23-acetate (from Alisma plantago-aquatica L.), as pivotal bioactives. We pinpointed 143 targets common to both UR-AP and hypertension, highlighting MAPK1, IL6, AKT1, VEGFA, EGFR, and TP53 as central targets involved in key pathways like diastolic and endothelial function, anti-atherosclerosis, AGE-RAGE signaling, and calcium signaling. Cluster analysis emphasized IL6, TNF, AKT1, and VEGFA’s roles in atherosclerosis and inflammation. Molecular docking confirmed strong interactions between these targets and UR-AP’s main bioactives, underscoring their therapeutic potential.
Conclusion: This research delineates UR-AP’s pharmacological profile in hypertension treatment, linking traditional medicine with modern pharmacology. It highlights key bioactive components and their interactions with principal targets, suggesting UR-AP’s potential as a novel therapeutic option for hypertension. The evidence from molecular docking studies supports these interactions, indicating the relevance of these components in affecting hypertension pathways. However, the study acknowledges its limitations, including the reliance on in silico analyses and the need for in vivo validation. These findings pave the way for future clinical research, aiming to integrate traditional medicine insights with contemporary scientific approaches for developing innovative hypertension therapies.
1 Introduction
Globally, hypertension ranks as the most prevalent and burdensome chronic non-communicable disease (Grassi et al., 2024). The NCD Risk Factor Collaborative Group, analyzing 1,201 studies with 104 million participants from 1990 to 2019, reported a significant increase in hypertension prevalence (Collaboration et al., 2021). Furthermore, hypertension is known to cause damage to various target organs, leading to conditions such as stroke, kidney disease, and ischemic heart disease (Zhou et al., 2021). Consequently, enhancing the efficacy of hypertension treatment and actively controlling its risk factors are imperative to mitigate organ damage and improve patient quality of life.
Traditional Chinese Medicine (TCM) offers several advantages, including low drug resistance, therapeutic stability, protection of target organs, symptom relief, improved long-term survival chances, and fewer adverse effects (Li, 2022). Given these benefits, East Asian countries, particularly China, are increasingly focusing on medications derived from traditional medicine. In this context, herbal prescriptions for treating hypertension warrant attention. Data mining-based studies indicated that liver-kidney co-therapy is a common strategy in managing hypertensive disorders (Jie et al., 2021; Wenlong, 2022). Uncaria rhynchophylla (UR), targeting the liver meridian, and Alisma plantago-aquatica L. (AP), targeting the kidney meridian, from the UR-AP pair, emblematic of this approach Moreover, UR and AP are frequently utilized in prescriptions by renowned veteran Chinese medicine practitioners, as evidenced in data mining studies (Xiaotian et al., 2021; Yichu, 2021; Shudong, 2022). However, the complex chemical composition and obscure molecular mechanism of UR-AP limit its clinical application, necessitating a deeper understanding of its mechanism in treating hypertension.
Network pharmacology, analyzing drug-target-disease interactions (Jiashuo et al., 2022), is particularly suitable for Chinese medicines characterized by component complexity and multiple interaction pathways. Therefore, network pharmacology cannot be ignored in the mechanism research of TCM. This approach aligns well with the multi-target network philosophy of TCM. Moreover, molecular docking, which simulates ligand-receptor binding to assess interaction stability and optimal binding modes, effectively complements network pharmacology. Therefore, this study aims to identify the bioactive components, core targets, and key signaling pathways of UR-AP in combating hypertension by constructing a Protein-Protein Interaction (PPI) network, performing cluster and enrichment analysis, and validating interactions through molecular docking (Figure 1).
2 Methods
2.1 Identification of bioactive components in UR-AP
Bioactive components of UR-AP were identified using the TCMSP (http://tcmspw.com/tcmsp.php) (Ru et al., 2014), TCMID (http://119.3.41.228:8000/tcmid/) (Huang et al., 2018), and DCABM-TCM (http://bionet.ncpsb.org.cn/dcabm-tcm/) (Liu et al., 2023). The ADME model, alongside published literature (Bai et al., 2021a; Bai et al., 2021b; Zhai et al., 2021; Liu et al., 2022), guided the setting of thresholds for oral bioavailability (OB) ≥ 30% and drug-likeness (DL) ≥ 0.18 to screen for bioactive components.
2.2 Identification of core targets of UR-AP
Core targets related to UR-AP’s bioactive ingredients were gathered from TCMSP, Swiss Target Prediction (http://www.swisstargetprediction.ch/) (Gfeller et al., 2013), and PharmMapper (http://www.lilab-ecust.cn/pharmmapper/) (Wang et al., 2017). Canonical SMILES format was obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/) (Kim et al., 2023). Gene mapping for the corresponding target was performed using Uniprot (https://www.uniprot.org/) (Coudert et al., 2023; UniProt et al., 2023), where gene name, IDs, and symbols were standardized and duplicates removed.
2.3 HT-related targets
Targets related to hypertension (HT) were collated from GeneCards (https://www.genecards.org/) (Stelzer et al., 2016), the DisGeNet database (https://www.disgenet.org/) (Pinero et al., 2020), Online Mendelian Inheritance in Man (OMIM, https://www.omim.org/) (Rasmussen et al., 2020), and Drug Bank (https://go.drugbank.com/) (Wishart et al., 2006). “Hypertension” and “essential hypertension” served as keywords. Data from different databases were merged after removing duplicates, and standardization of HT-related target names was conducted using UniProt.
2.4 UR-AP-HT common-target network
Integration of data and construction of the common-target network were accomplished using Cytoscape 3.7.1. The Network Analyzer plugin computed the degree values. A Venn diagram, created using R 4.2.1, facilitated the identification common targets of UR-AP and HT. After that, import the common targets into Cytoscape to build the UR-AP-HT common target network. Analysis of the main targets and bioactive components of UR-AP against HT according to the degree.
2.5 PPI network and evaluation
To illuminate the role of targets at the phylum level, the String database (https://cn.string-db.org/) (Szklarczyk et al., 2021) was used to produce the protein-protein interaction (PPI) network of UR-AP-HT common targets. Set the minimum confidence score to 0.400 to obtain a highly reliable PPI network relationship. Afterward, degree (DC), betweenness (BC), closeness (CC), eigenvector (EC), information (IC), local average connectivity-based method (LAC), network (NC), subgragh (SC) in CytoNCA were used to screen core targets.
2.6 Cluster analysis
MCODE is a plug-in for applying the K-means clustering algorithm. Further analysis of the PPI network by applying the MCODE plug-in could not only obtain the clusters composed of similar targets but also infer the main targets. The conditions were set to Degree Cutoff = 2, Node Score Cutoff = 0.2, K-core = 2, and Max. Depth = 100.
2.7 GO and KEGG pathway enrichment analyses
Gene Set Enrichment Analysis (GSEA) aggregated functionally similar gene individuals, allowing for filtering and screening of relevant functionality and pathway information. Gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analyses were conducted using the GO knowledgebase and KEGG database (Liu et al., 2017). The GO knowledgebase (http://geneontology.org/) (Gene et al., 2021) covers cellular components, molecular functions, and biological processes. The KEGG database (https://www.genome.jp/kegg/) (Kanehisa et al., 2023) integrates information from genomics, biochemistry, and systemic functional omics, commonly used in the KEGG pathway. The BiocManager package of R 4.2.1 manages several other packages, among which the clusterProfiler package is used to interpret the functional characteristics of genomics data and provide up-to-date gene annotations (Wu et al., 2021). The org.Hs.e.g.,.db package was used to provide human genome-wide annotations. The screening was based on p < 0.05.
2.8 Molecular docking
Molecular docking assessed the binding energy and optimal binding modes between ligand (bioactive components) and receptor (core targets). The core targets and main bioactive components of UR-AP against HT were respectively selected as receptors and ligands for molecular docking. The RCSB Protein Data Bank (RCSB PDB) (http://www.rcsb.org/) (Burley et al., 2021) was used to download the 3D crystal structures and save them in pdb format, including MAPK1 (PDB ID:8AOJ), IL6 (PDB ID:1ALU), AKT1 (PDB ID:1UNQ), VEGFA (PDB ID:4KZN), EGFR (PDB ID:8A27), TP53 (PDB ID:3D06). Download mol2 structure files of bioactive components from TCMSP and convert them into pdb format files using OpenBabel 3.1.1. Remove water and add all hydrogens by using PyMol 2.2.0. Afterward, the pdb format file was imported into AutoDock Tools 4.2 for further preprocessing of receptors and ligands and exported as the qdbqt format file. Molecular docking using AutoDock Vina 1.2.0. The docking parameters were set to modes = 20 and exhaustiveness = 10. Binding energy less than -5kj/mol can be considered feasible for docking results (Li et al., 2019). Visualization of docking results with 2D and 3D diagrams was achieved by using LigPlot+ 2.2.7 and PyMol 2.2.0.
3 Results
3.1 UR-AP component-target network
The components of UR-AP were obtained from TCMSP, TCMID, DCABM-TCM databases (Supplementary Table S1). After that, a total of 58 bioactive components were obtained after de-weighting and screening (Supplementary Table S2). These included 43 components from UR, 16 from AP, and 1 common component (sitosterol). Action targets for these 58 bioactive components were predicted using TCMSP, Swiss Target Prediction, and PharmMapper. In total, 322 targets of UR-AP (except cyanide) were predicted after removing duplicates (Supplementary Table S3). The UR-AP common-target network comprised 380 nodes and 1062 edges (Figure 2), with quercetin (degree = 151), alisol B (degree = 95), alisol,b,23-acetate (degree = 81), and kaempferol (degree = 62) identified as the primary bioactive components (Supplementary Table S4).
FIGURE 2. UR-AP component-target network: The pink nodes represent the bioactive components of UR, the yellow nodes represent the bioactive components of AP, and the blue nodes represent the core targets of UR-AP.
3.2 HT-related targets
To ensure the comprehensiveness of HT-related targets obtained, 1141 HT-related targets were acquired from 4 databases. Of these, 945 targets were obtained from GeneCards, 416 targets from DisGeNET, 25 targets from OMIM, and 42 targets from Drug Bank (Supplementary Table S5).
3.3 UR-AP-HT common-target network
The intersection of UR-AP core targets and HT-related targets identified 143 common targets (Figure 3A). The resulting network consisted of 183 nodes and 477 edges (Figure 3B), involving 32 bioactive components from UR, 9 from AP, and 1 common component. The main targets and bioactive components were analyzed according to degree values (Supplementary Table S6). The results indicated that there were four critical bioactive components for UR-AP treatment of HT, quercetin (degree = 84) and kaempferol (degree = 35) from UR, and alisol,b,23-acetate (degree = 38) and alisol B (degree = 36) from AP. In addition, PTGS2, AR, PTGS1, KCNH2, ADRA1B, and F2R were the main targets of UR-AP for HT treatment (Table 1) (Figure 3C).
FIGURE 3. UR-AP-HT common-target network. (A) The intersection of the Venn diagram: 143 targets are common to HT and UR-AP. (B) UR-AP-HT common-target network. (C) Main bioactive components and potential targets. The purple nodes represent the bioactive components of UR, the green nodes represent the bioactive components of AP, and the blue nodes represent the bioactive components common to both UR and AP. The red nodes represent core targets. The node sizes were ordered by degree. The thickness and color shading of the edges were set according to the degree of interaction between the nodes.
3.4 UR-AP-HT PPI network
The 143 common targets were imported into String database, and then 1 free target was discarded. Afterward, the PPI network included 142 nodes and 2296 edges (Figure 4A). The topological parameters of the UR-AP-HT PPI network were calculated, and the significant targets of the UR-AP-HT PPI network were obtained based on DC, BC, CC, EC, IC, LAC, NC, and SC (Supplementary Table S7) (Figures 4B, C). Finally, six important targets were obtained, namely, Interleukin-6 (IL6), Tumor necrosis factor (TNF), RAC-alpha serine/threonine-protein kinase (AKT1), Vascular endothelial growth factor A (VEGFA), Interleukin-1 beta (IL1B), Cellular tumor antigen p53 (TP53) (Figure 4C).
FIGURE 4. UR-AP-HT PPI network. (A) Total UR-AP-HT PPI network. (B) PPI network by the screening criteria of DC ≥ 55, BC ≥ 48.71780123, CC ≥ 0.534090909, EC ≥ 0.059021739, IC ≥ 10.76018715, LAC ≥18.775, NC ≥ 20.805268, SC ≥ 18778800000000000000. (C) PPI network by the screening criteria of DC ≥ 80, BC ≥ 366.8913084, CC ≥ 0.684466019, EC ≥ 0.159114137, IC ≥ 14.1830368, LAC ≥36.4516129, NC ≥ 69.88953295, SC ≥ 136440000000000000000. The size and color shading of the nodes were ordered by degree. The thickness and color shading of the edges were set according to the degree of interaction between the nodes.
3.5 Cluster analysis
In order to further analyze the 143 common targets of UR-AP and HT, find the key sub-networks, and speculate the main targets, the UR-AP-HT PPI network was further computationally analyzed by MCODE. The PPI network was divided into 6 clusters (Supplementary Table S8) (Figures 5A–F). The score of cluster 1 (score = 37.023) was far higher than the other clusters. The degree of the targets in cluster 1 was higher, such as IL6 (degree = 98), TNF (degree = 93), AKT1 (degree = 88), VEGFA (degree = 88), IL1B (degree = 81), PTGS2 (degree = 78), PPARG (degree = 78).
FIGURE 5. Cluster analysis of UR-AP-HT PPI network. (A) Cluster 1, composed of 44 nodes and 796 edges (score = 37.023). (B) Cluster 2, composed of 23 nodes and 139 edges (score = 12.636). (C) Cluster 3, composed of 5 nodes and 7 edges (score = 3.5). (D) Cluster 4, composed of 3 nodes and 3 edges (score = 3). (E) Cluster 5, composed of 3 nodes and 3 edges (score = 3). (F) Cluster 6, composed of 3 nodes and 3 edges (score = 3). The size and color shading of the nodes were ordered by degree. The thickness and color shading of the edges were set according to the degree of interaction between the nodes. The arrow pointed to the target.
3.6 GO and KEGG pathway enrichment analyses of PPI network
Enrichment analysis of the UR-AP-HT PPI network was obtained by applying the BiocManager package of R 4.2.1 (Supplementary Table S9). In GO enrichment analysis, 143 UR-AP-HT common targets were highly enriched in 2569 biological processes (BP), 82 cellular fractions (CC), and 206 molecular functions (MF) by setting p-values less than 0.05. The top 20 significant results for each were taken as shown in the figure (Figures 6A–C). Specifically, BP had a tight correlation with response to decreased oxygen levels, response to oxygen levels, and response to hypoxia, which indicates a possible role of UR-AP on HT in antihypoxia. Besides, there was also a strong correlation with vascular process in circulatory system, regulation of tube diameter, blood vessel diameter maintenance, regulation of tube size, response to nutrient levels, and regulation of blood pressure, suggesting that UR-AP may regulate blood pressure by modulating tubular diameter. CC is closely correlated with membrane raft and membrane microdomain, reflecting the possibility that UR-AP regulates blood pressure by acting on the cell membrane. MF was associated with cytokine receptor binding, cytokine activity, signaling receptor activator activity, and receptor ligand activity, suggesting that UR-AP may affect cytokine, signaling receptor, and receptor ligand activities. In addition, a total of 156 pathways were enriched by KEGG enrichment analysis, and the top 20 important pathways were shown in the figure (Figure 6D). The results indicated that UR-AP may regulate atherosclerosis through Fluid shear stress and atherosclerosis, Lipid and atherosclerosis, as well as through AGE-RAGE signaling pathway in diabetic complications and calcium signaling pathway to regulate endocrine (Table 2). UR-AP also regulates signaling pathways such as the HIF-1 signaling pathway and TNF signaling pathway.
FIGURE 6. GO and KEGG pathway enrichment analyses of PPI network. (A) Bubble diagram of the top 20 biological processes. (B) Bubble diagram of the top 20 cellular components. (C) Bubble diagram of the top 20 molecular functions. The color was adjusted according to −1* log10 (p-value) and the size of the dots was adjusted according to the number of genes. (D) Bidirectional histogram of the top 20 KEGG pathways. Blue bars represent −1*log10 (p-value) and red bars represent the number of genes.
3.7 GO and KEGG pathway enrichment analyses of the cluster
GO and KEGG enrichment analysis was performed on cluster 1 (score = 37.023), which had the highest score in the cluster analysis, to further obtain the functional information and pathway information of cluster 1 (Supplementary Table S10). In the GO enrichment analysis, 44 targets in cluster 1 were highly enriched in 2026 biological processes, 27 cellular components, and 115 molecular functions by setting p-values less than 0.05. The top 20 significant results for each were taken as shown in the figure (Figure 7A). Specifically, BP had a strong correlation with response to lipopolysaccharide, cellular response to lipopolysaccharide, which suggests a possible effect on lipopolysaccharide. As well as a close correlation with response to reactive oxygen species, response to oxidative stress, and response to decreased oxygen levels, suggesting a possible antioxidant effect of UR-AP against HT. CC is associated with membrane raft, membrane microdomain, etc., suggesting that it may act on the cell membrane to exert the effect. MF is tightly related to cytokine receptor binding, cytokine activity, signaling receptor activator activity, and receptor ligand activity, indicating that UR-AP may influence cytokines and signaling in the pathogenesis of HT. Furthermore, a total of 132 pathways were enriched by KEGG enrichment analysis, and the top 20 important pathways were shown in the figure (Figure 7B), and the top 5 important pathways and gene distribution were displayed (Figure 7C) (Table 3).
FIGURE 7. GO and KEGG pathway enrichment analyses of cluster 1. (A) The top 20 GO enrichment analysis results. The red ones represent the results of biological processes and show the ID of -log10 (p-value) > 18. The blue ones represent the results of cellular components and show the ID of -log10 (p-value) > 6. The green ones represent the results of molecular functions and show the ID of -log10 (p-value) > 15. The size of the points was adjusted according to the number of genes. (B) Histogram of the top 20 KEGG pathways. Color shading was adjusted according to −1*log10 (p-value). (C) String diagram of KEGG enrichment analysis results for cluster 1.
3.8 Core component-target-pathway network
For further revealing the action mechanism of UR-AP for HT, 20 KEGG pathways, related targets, and bioactive components were constructed into a systematic and complex network, including 128 nodes and 731 edges (Figure 8). The core component-target-pathway network of UR-AP for HT was analyzed according to the degree values (Supplementary Table S11). The results indicated that there were mainly 3 key bioactive components in the core component-target-pathway network, including quercetin (degree = 40) and kaempferol (degree = 24) from UR, and Alisol B (degree = 21) from AP. Moreover, MAPK1, PIK3CA, AKT1, EGFR, VEGFA, MAPK8, BCL2, FOS, IL6, and TP53 were the top 10 significant targets.
FIGURE 8. Core component-target-pathway network. The yellow nodes represent HT. The orange nodes represent UR and AP. The red nodes represent the bioactive components of UR, and the green nodes represent the bioactive components of AP. The blue nodes represent relevant targets. The purple nodes represent the top 20 signaling pathways from the enrichment analysis. The node sizes were ordered by degree. The thickness and color shading of the edges were set according to the degree of interaction between the nodes.
The core component-target-pathway network identified relevant signaling pathways including Fluid shear stress and atherosclerosis (degree = 28), AGE-RAGE signaling pathway in diabetic complications (degree = 28), PI3K-Akt signaling pathway (degree = 27), Lipid and atherosclerosis (degree = 27). One of the most significant pathways considered for UR-AP in treating hypertension according to the degree value is Fluid shear stress and atherosclerosis, indicating that UR-AP may exert hypotensive and antiatherogenic effects by acting on vascular endothelial cells to exert stable laminar shear stress, thus participating in vasodilation, antioxidant, anti-inflammatory and antithrombotic Sclerosis. The action path is shown in the (Figure 9).
3.9 Molecular docking
Combining the top 10 most important targets from the results of the core component-target-pathway network (MAPK1, PIK3CA, AKT1, EGFR, VEGFA, MAPK8, BCL2, FOS, IL6, TP53) with the top 10 most significant targets from the results of the PPI network and the cluster analysis (IL6, TNF, AKT1, VEGFA, IL1B, TP53, PTGS2, PPARG, CASP3, EGFR), we found that the duplicated targets were IL6, AKT1, VEGFA, EGFR, TP53. Due to the fact that MAPK1 was the most prominent target in the core component-target-pathway network, we concluded that MAPK1, IL6, AKT1, VEGFA, EGFR, and TP53 were the core targets of UR-AP against HT. In addition, there were 8 duplications between each of the top 10 most important core bioactive components of the core component-target-pathway network and UR-AP-HT common-target network. Therefore, we selected that these 8 duplicated bioactive components (quercetin, alisol,b,23-acetate, alisol B, kaempferol, yohimbine, beta-sitosterol, methyl (E)-2-[(2S,3Z,12bS)-3-ethylidene ethylidene-2,4,6,7,12,12b-hexahydro-1H-indolo [3,2-h]quinolizin-2-yl]-3-methoxyprop-2-enoate, Tetrahydroalstonine) as small molecules (ligands), and MAPK1, IL6, AKT1, VEGFA, EGFR, TP53 as proteins (receptors) to perform molecular docking, with the aim of validating the interactions between core targets and bioactive components. Finally, 48 groups of receptor-ligand molecular docking results were obtained (Supplementary Table S12).
Binding energy less than −7.0 kcal/mol indicates superior binding activity, and the lower the binding energy, the better the docking effectiveness (Trott and Olson, 2010). Among the 48 groups of receptor-ligand docking results, there were 33 groups with the binding energy less than −7.0 kcal/mol, which indicated the strong binding interactions between the core targets and the bioactive components of UR-AP for HT treatment (Figure 10). The five groups with the lowest binding energy were selected to draw molecular docking models (Table 4) (Figure 11). The results indicated that beta-sitosterol was docked with 13 residues in EGFR (PRO914, SER912, TRP880, GLU906, ARG803, LEU799, LYS913, ARG841, CYS797, GLY719, NO31105, SER720, VAL726), and formed hydrogen bonds with ASP855 (2.91Å) and LYS745 (3.18Å) (Figures 11A–C). Tetrahydroalstonine was docked with seven residues forming hydrophobic interactions in EGFR (ASP761, ALA859, GLU762, GLU758, LEU858, KY91102, PHE723) (Figures 11D–F). Moreover, Tetrahydroalstonine was docked with six residues in TP53 (ARG158, ASP208, LEU206, MET160, ILE254, PRO98) and the hydrogen bonds were formed at PRO98 (3.22 Å) and THR256 (2.81 Å) (Figures 11G–I). Alisol,b,23-acetate was docked with 10 residues in EGFR (ASP800, GLY719, SER720, LEU718, EDO1104, VAL726, LEU844, CYS797, ASP855, LEU799), and formed hydrogen bonds with ARG803 (3.33Å, 3.00 Å, 2.84Å) and ARG841 (3.19Å) (Figures 11J–L). Alisol B was docked with 13 residues to form hydrophobic interactions in EGFR (PHE723, ALA722, GLY721, NO31105, ARG841, ASP855, SER720, GLY719, CYS797, LEU844, VAL726, EDO1104, LEU718) (Figures 11M–O). In addition, the docking energy of Alisol B with EGFR was the lowest (−8.9 kcal/mol) compared to other docking results, indicating that the Alisol B-EGFR complex is more stable than other targets.
FIGURE 11. Molecular docking models of the lowest binding energy. (A–C) Molecular docking models of beta-sitosterol and EGFR. (D–F) Molecular docking models of Tetrahydroalstonine and EGFR. (G–I) Molecular docking models of Tetrahydroalstonine and TP53. (J–L) Molecular docking models of alisol,b,23-acetate, and EGFR. (M–O) Molecular docking models of Alisol B and EGFR. (A, D, G, J, M) Diagrams of the overall effect of molecular docking. (A, D, G, J, M) 2D diagrams of molecular docking. (A, D, G, J, M) 3D diagrams of molecular docking.
4 Discussion
4.1 Hypertension: a global health challenge and the role of TCM
The rising incidence of hypertension, driven by factors such as dietary habits, lifestyle changes, and emotional stress, continues to be a global health concern (Rahimi et al., 2017). Symptoms like dizziness and headache significantly impact the quality of life of many patients. TCM, known for its efficacy in controlling blood pressure and alleviating symptoms, highlights the potential of developing new herbal treatments for HT. UR-AP, a common drug pair in TCM for HT treatment, is particularly notable. Although Uncaria rhynchophylla is widely recognized for its effectiveness in HT management, the complexities of the chemical composition of UR-AP and the specific components and targets involved in its action against hypertension are not fully understood (Dong et al., 2020; Tian et al., 2020; Deng et al., 2022). Therefore, this study is the first to use the network pharmacology approach to reveal the action mechanism of UR-AP in treating hypertension and to verify the interaction between core targets and bioactive components of UR-AP against hypertension with molecular docking.
4.2 Bioactive components and targets of UR-AP in hypertension management
In this study, 58 bioactive components and 322 targets of UR-AP were obtained by screening the online database. Additionally, 1,141 HT-related targets were gathered from disease target databases. The intersection of these targets using a Venn diagram led to the identification of 143 common targets, forming the basis for the UR-AP-HT common targets network. This network revealed that UR-AP’s therapeutic effects on hypertension are likely due to the synergistic action of various bioactive components. Based on the results of the UR-AP-HT common-target network and core component-target-pathway network, we first predicted that quercetin, kaempferol, beta-sitosterol from UR and Alisol B, alisol,b,23-acetate from AP were the main bioactive components of UR-AP against hypertension. The anti-hypertensive and cardioprotective effects of quercetin have been extensively studied. Quercetin can play a role in anti-hypertension by suppressing Ang II-stimulated VSMC proliferation (Wang et al., 2022). Quercetin may reduce oxidative stress, and minimize arterial endothelial cell dysfunction (Pereira et al., 2018). Quercetin was able to promote endothelial cell autophagy and reduce the blood pressure value of SHRs (Lin et al., 2020). Kaempferol (KFL) is a naturally occurring flavonoid with pleiotropic cardioprotective and antihypertensive effects. The pleiotropic cardioprotective effect of kaempferol can be attributed to the anti-inflammatory activity in vitro and animal studies (Dabeek and Marra, 2019; Zhao et al., 2023). Kaempferol can increase the phosphorylation of eNOS via the NO-cGMP-PKG pathway and induce endothelium-dependent vasorelaxation (Tettey et al., 2019). Beta-sitosterol has cholesterol-lowering and tissue-repairing effects and can be used to alleviate hyperlipidemia and atherosclerosis (Khan et al., 2022). Both alisol,b,23-acetate (AB23A) and alisol B, as alisol extracts, are effective ingredients in the treatment of atherosclerosis. AB23A and Alisol B may interfere with lipids and atherosclerosis by modulating PI3K-AKT and MAPK pathways, thereby treating Pregnancy-induced Hypertension and reducing liver and kidney injury (Liao et al., 2022). AB23A can regulate bile acid metabolism, lower cholesterol and alleviate atherosclerosis by increasing FXR-BSEP signaling (Fu et al., 2022).
To bridge the gap between these findings and clinical practice, it is essential to consider the feasibility and potential challenges of translating the identified gene targets into therapeutic interventions. The extensive research on the cardioprotective and antihypertensive effects of quercetin and kaempferol underscores their potential as foundational compounds in novel hypertension treatments. However, the clinical application of these bioactive components requires careful consideration of their bioavailability, pharmacokinetics, and safety profiles in human populations. Future studies should focus on the development of delivery mechanisms that enhance the efficacy of these compounds, rigorous clinical trials to establish therapeutic dosages, and comprehensive safety evaluations. This approach will ensure that the promising results of network pharmacology studies can be effectively translated into tangible benefits for patients with hypertension.
4.3 Mechanisms of action: core targets and pathways influenced by UR-AP
Synthesizing the results of the core component-target-pathway network, PPI network, and cluster analysis, we predicted that the core targets for UR-AP treatment of HT are MAPK1, IL6, AKT1, VEGFA, EGFR, and TP53. These have abundant interactions with other targets and are involved in additional signaling pathways. MAPK1 and AKT1 are able to facilitate vascular endothelial cell proliferation and angiogenesis, which is closely related to hypertension (Zhang et al., 2019). P53 is the protein product of TP53, which is implicated in atherosclerosis, VSMC growth, and cell death (Han et al., 2021). The p38 MAPK pathway is involved in blood pressure regulation and vascular injury (El-Sabbagh et al., 2023). The Akt pathway can be engaged in the differentiation of VSMC (Li et al., 2024). IL6 promotes smooth muscle cell proliferation and is implicated in hypertension through interactions with cytokines (Attiq et al., 2024). VEGFA is positively correlated with elevated blood pressure and is strongly associated with vascular damage and atherosclerosis (Braile et al., 2020). EGFR is related to cell proliferation, angiogenesis, and apoptosis inhibition (Luo et al., 2021).
4.4 Impact of UR-AP on vascular function and atherosclerosis
KEGG pathway enrichment analysis and core component-target-pathway network identified important signaling pathways involved in vasoconstriction, endothelial function, atherosclerosis, inflammatory response, and apoptosis (Figure 12). AGE-RAGE signaling pathway may mediate diabetic vascular complications (Burr et al., 2020). The binding of NAPDH, a product of RAGE, to AGEs enhances oxidative stress, stimulates the production of cytokines and growth factors, affects smooth muscle proliferation and apoptosis, causes atherosclerosis as well as vasoconstriction, and thus increases blood pressure (Eid et al., 2018). In addition, the combination of AGE and RAGE activates signaling pathways such as PI3K-Akt pathway and Calcium pathway, which can further activate NAPDH and lead to cell damage or even apoptosis (Sanajou et al., 2018). Calcium signaling pathway plays an important role in vasoconstriction and diastole. Elevated intracellular Ca2+ concentration induces vasoconstrictive effects (Gutierrez et al., 2023). Multiple signaling pathways can activate the PI3K-Akt signaling pathway, which can lead to increased transcription of NF-κB, and can further regulate the expression of cytokines, growth factors, and adhesion molecules (e.g., VEGF, IL6, VCAM-1), leading to cellular dysfunction (Li et al., 2023). The development of atherosclerosis is caused by impaired endothelial function, inflammation, and oxidative processes (Wolf and Ley, 2019). Fluid shear stress and atherosclerosis signaling pathway can activate PI3K-Akt signaling pathway, generate eNOS, and release NO by exerting stable laminar shear stress on vascular endothelial cells, thus participating in atherosclerosis, vasoconstriction, inflammatory response, cell proliferation, and angiogenesis, etc (Cheng et al., 2023). Lipid and atherosclerosis signaling pathways can be involved in the atherosclerotic process, allowing LDL-C to enter the vascular endothelium (Halasz et al., 2021). Vascular smooth muscle cells migrate and proliferate, which in turn form plaques. Damage to the vascular endothelium can cause an increase in blood pressure.
FIGURE 12. Illustration of the signaling pathways associated with the major bioactive components of UR-AP for the treatment of HT.
4.5 Cluster analysis and molecular docking: from prediction to validation
Cluster analysis was performed to analyze the common targets of UR-AP and HT, in order to find the key sub-networks and speculate the main targets. The results showed a high score for cluster 1 (score = 37.023). We predicted IL6, TNF, AKT1, VEGFA, IL1B, PTGS2, and PPARG as the main targets based on the degree. GO and KEGG pathway enrichment analysis of cluster 1 showed that cluster 1 may exert its efficacy by acting on the cell membrane and has an effect on lipopolysaccharides. Important signaling pathways including AGE-RAGE signaling pathway in diabetic complications, Fluid shear stress and atherosclerosis, TNF signaling pathway, IL-17 signaling pathway, Lipid and atherosclerosis were identified, which are mainly involved in atherosclerosis and inflammatory response.
Ultimately, we selected six core targets of UR-AP against HT and eight bioactive components for molecular docking to validate the interactions. The results showed strong binding interactions between the core targets and bioactive components with excellent stability.
4.6 Limitations and prospects for future research
This study presents an in-depth analysis of the mechanism of action of UR-AP in the context of hypertension, aiming to contribute significantly to TCM research, broaden clinical applications, and facilitate drug development. Despite its promise, challenges such as the intricate material basis of TCM and the gap between fundamental research and clinical application persist. Notably, the study faces limitations like the inability to perform experimental validations due to resource and financial constraints. Moving forward, our team is dedicated to keeping abreast of hypertension research advancements and plans to undertake animal or cellular experiments to further investigate UR-AP’s targets and mechanisms for hypertension. Additionally, pioneering efforts in network pharmacology will be made to elucidate the specific target modules of TCM’s effective population against hypertension, aiming to uncover TCM’s clinical efficacy and pharmacological mechanisms comprehensively.
5 Conclusion
This study systematically unraveled the intricate mechanisms underlying the therapeutic action of UR-AP in the management of HT, utilizing a novel integration of network pharmacology, cluster analysis, and molecular docking techniques. Key bioactive components from UR such as quercetin, kaempferol, and beta-sitosterol, along with Alisol B and alisol,b,23-acetate from AP, emerged as pivotal in mediating these beneficial effects. The study identified core targets including MAPK1, IL6, AKT1, VEGFA, EGFR, and TP53, highlighting their critical roles in the therapeutic efficacy of UR-AP.The investigation brought to light several vital signaling pathways, notably the AGE-RAGE signaling pathway, Fluid shear stress and atherosclerosis, Lipid and atherosclerosis, PI3K-Akt signaling pathway, and Calcium signaling pathway, which are instrumental in the action of UR-AP against HT. Our cluster analysis further pinpointed IL6, TNF, AKT1, and VEGFA as key targets, suggesting their crucial role in influencing cell membrane dynamics and responses to lipopolysaccharides. Predominantly, these pathways and targets are associated with the pathophysiology of atherosclerosis and inflammatory responses.Moreover, the molecular docking studies provided compelling evidence of strong binding interactions between the core targets and UR-AP’s major bioactive components, underscoring the potential of these compounds in HT therapy.
Building on this foundation, it is acknowledged that the study has its limitations, which have been discussed in a dedicated section within the document. Future research endeavors will seek to address these limitations through further experimental validations and explorations, particularly via animal models and cellular experiments, as well as employing cutting-edge techniques in network pharmacology to delve deeper into the action mechanisms of UR-AP and its potential applications in hypertension treatment.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
TY: Writing–original draft. HZ: Writing–original draft. XL: Writing–original draft. DW: Visualization, Writing–review and editing. CR: Data curation, Writing–review and editing. LC: Visualization, Writing–review and editing. YL: Data curation, Writing–review and editing. LnW: Visualization, Writing–review and editing. JW: Data curation, Writing–review and editing. ZZ: Data curation, Writing–review and editing. DL: Writing–review and editing. LyW: Writing–review and editing. XH: Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the science and technology innovation project of China Academy of Chinese Medical Sciences (grant number: CI2021A00105), independently selected project of China Academy of Chinese Medical Sciences (grant number: Z0814), and the Youth Program of National Natural Science Foundation of China (grant number: 82205320).
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/fchem.2024.1356458/full#supplementary-material
References
Attiq, A., Afzal, S., Ahmad, W., and Kandeel, M. (2024). Hegemony of inflammation in atherosclerosis and coronary artery disease. Eur. J. Pharmacol. 966, 176338. doi:10.1016/j.ejphar.2024.176338
Bai, X., Tang, Y., Li, Q., Chen, Y., Liu, D., Liu, G., et al. (2021a). Network pharmacology integrated molecular docking reveals the bioactive components and potential targets of Morinda officinalis-Lycium barbarum coupled-herbs against oligoasthenozoospermia. Sci. Rep. 11, 2220. doi:10.1038/s41598-020-80780-6
Bai, X., Tang, Y., Li, Q., Liu, D., Liu, G., Fan, X., et al. (2021b). An integrated analysis of network pharmacology, molecular docking, and experiment validation to explore the new candidate active component and mechanism of cuscutae semen-mori fructus coupled-herbs in treating oligoasthenozoospermia. Drug Des. Devel Ther. 15, 2059–2089. doi:10.2147/dddt.s307015
Braile, M., Marcella, S., Cristinziano, L., Galdiero, M. R., Modestino, L., Ferrara, A. L., et al. (2020). VEGF-A in cardiomyocytes and heart diseases. Int. J. Mol. Sci. 21, 5294. doi:10.3390/ijms21155294
Burley, S. K., Bhikadiya, C., Bi, C., Bittrich, S., Chen, L., Crichlow, G. V., et al. (2021). RCSB Protein Data Bank: powerful new tools for exploring 3D structures of biological macromolecules for basic and applied research and education in fundamental biology, biomedicine, biotechnology, bioengineering and energy sciences. Nucleic Acids Res. 49, D437–D451. doi:10.1093/nar/gkaa1038
Burr, S. D., Harmon, M. B., and Jr, J. A. S. (2020). The impact of diabetic conditions and AGE/RAGE signaling on cardiac fibroblast migration. Front. Cell Dev. Biol. 8, 112. doi:10.3389/fcell.2020.00112
Cheng, H., Zhong, W., Wang, L., Zhang, Q., Ma, X., Wang, Y., et al. (2023). Effects of shear stress on vascular endothelial functions in atherosclerosis and potential therapeutic approaches. Biomed. Pharmacother. 158, 114198. doi:10.1016/j.biopha.2022.114198
Collaboration, N. C. D. R. F., Carrillo-Larco, R. M., Danaei, G., Riley, L. M., Paciorek, C. J., Stevens, G. A., et al. (2021). Worldwide trends in hypertension prevalence and progress in treatment and control from 1990 to 2019: a pooled analysis of 1201 population-representative studies with 104 million participants. Lancet 398, 957–980. doi:10.1016/s0140-6736(21)01330-1
Coudert, E., Gehant, S., de Castro, E., Pozzato, M., Baratin, D., Neto, T., et al. (2023). Annotation of biologically relevant ligands in UniProtKB using ChEBI. Bioinformatics 39, btac793. doi:10.1093/bioinformatics/btac793
Dabeek, W. M., and Marra, M. V. (2019). Dietary quercetin and kaempferol: bioavailability and potential cardiovascular-related bioactivity in humans. Nutrients 11, 2288. doi:10.3390/nu11102288
Deng, L., Liu, W., Xu, Q., Guo, R., Zhang, D., Ni, J., et al. (2022). Tianma Gouteng Decoction regulates oxidative stress and inflammation in AngII-induced hypertensive mice via transcription factor EB to exert anti-hypertension effect. Biomed. Pharmacother. 145, 112383. doi:10.1016/j.biopha.2021.112383
Dong, H., Zhang, S., Du, W., Cong, H., and Zhang, L. (2020). Pharmacodynamics and metabonomics study of Tianma Gouteng Decoction for treatment of spontaneously hypertensive rats with liver-yang hyperactivity syndrome. J. Ethnopharmacol. 253, 112661. doi:10.1016/j.jep.2020.112661
Eid, B. G., Abu-Sharib, A. T., El-Bassossy, H. M., Balamash, K., and Smirnov, S. V. (2018). Enhanced calcium entry via activation of NOX/PKC underlies increased vasoconstriction induced by methylglyoxal. Biochem. Biophys. Res. Commun. 506, 1013–1018. doi:10.1016/j.bbrc.2018.10.171
El-Sabbagh, W. A., Fadel, N. A., El-Hazek, R. M., Osman, A. H., and Ramadan, L. A. (2023). Ubiquinol attenuates gamma-radiation induced coronary and aortic changes via PDGF/p38 MAPK/ICAM-1 related pathway. Sci. Rep. 13, 22959. doi:10.1038/s41598-023-50218-w
Fu, Y., Feng, H., Ding, X., Meng, Q. H., Zhang, S. R., Li, J., et al. (2022). Alisol B 23-acetate adjusts bile acid metabolisim via hepatic FXR-BSEP signaling activation to alleviate atherosclerosis. Phytomedicine 101, 154120. doi:10.1016/j.phymed.2022.154120
Gene, C., Douglass, E., Good, B. M., Unni, D. R., Harris, N. L., Mungall, C. J., et al. (2021). The Gene Ontology resource: enriching a GOld mine. Nucleic Acids Res. 49, D325–D334. doi:10.1093/nar/gkaa1113
Gfeller, D., Michielin, O., and Zoete, V. (2013). Shaping the interaction landscape of bioactive molecules. Bioinformatics 29, 3073–3079. doi:10.1093/bioinformatics/btt540
Grassi, G., Kreutz, R., and Mancia, G. (2024). The 2023 European Society of Hypertension guidelines for high blood pressure - what's new on diagnostic and therapeutic approaches. Pharmacol. Res. 199, 107003. doi:10.1016/j.phrs.2023.107003
Gutierrez, A., Gomez Del Val, A., Contreras, C., Olmos, L., Sanchez, A., and Prieto, D. (2023). Calcium handling coupled to the endothelin ET(A) and ET(B) receptor-mediated vasoconstriction in resistance arteries: differential regulation by PI3K, PKC and RhoK. Eur. J. Pharmacol. 956, 175948. doi:10.1016/j.ejphar.2023.175948
Halasz, G., Parati, G., and Piepoli, M. F. (2021). Editorial comments: focus on lipid and atherosclerosis. Eur. J. Prev. Cardiol. 28, 1975–1977. doi:10.1093/eurjpc/zwac010
Han, J. H., Heo, K. S., and Myung, C. S. (2021). Cytokine-induced apoptosis inhibitor 1 (CIAPIN1) accelerates vascular remodelling via p53 and JAK2-STAT3 regulation in vascular smooth muscle cells. Br. J. Pharmacol. 178, 4533–4551. doi:10.1111/bph.15631
Huang, L., Xie, D., Yu, Y., Liu, H., Shi, Y., Shi, T., et al. (2018). TCMID 2.0: a comprehensive resource for TCM. Nucleic Acids Res. 46, D1117–D1120. doi:10.1093/nar/gkx1028
Jiashuo, W. U., Fangqing, Z., Zhuangzhuang, L. I., Weiyi, J., and Yue, S. (2022). Integration strategy of network pharmacology in Traditional Chinese Medicine: a narrative review. J. Tradit. Chin. Med. 42, 479–486. doi:10.19852/j.cnki.jtcm.20220408.003
Jie, Z., Wende, T., Luxia, S., Jianqing, J., and Hao, X. (2021). Medication regularity of traditional Chinese patent medicine in treating hypertension based on data mining. Chin. J. Integr. Med. Cardio-cerebrovascular Dis. 19, 4216–4221. doi:10.12102/j.issn.1672-1349.2021.24.002
Kanehisa, M., Furumichi, M., Sato, Y., Kawashima, M., and Ishiguro-Watanabe, M. (2023). KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 51, D587–D592. doi:10.1093/nar/gkac963
Khan, Z., Nath, N., Rauf, A., Emran, T. B., Mitra, S., Islam, F., et al. (2022). Multifunctional roles and pharmacological potential of β-sitosterol: emerging evidence toward clinical applications. Chem. Biol. Interact. 365, 110117. doi:10.1016/j.cbi.2022.110117
Kim, S., Chen, J., Cheng, T., Gindulyte, A., He, J., He, S., et al. (2023). PubChem 2023 update. Nucleic Acids Res. 51, D1373–D1380. doi:10.1093/nar/gkac956
Li, J. (2022). Traditional Chinese medicine in treating hypertension. Circ. Cardiovasc Qual. Outcomes 15, e008723. doi:10.1161/circoutcomes.121.008723
Li, J., Dong, S., Quan, S., Ding, S., Zhou, X., Yu, Y., et al. (2023). Nuciferine reduces inflammation induced by cerebral ischemia-reperfusion injury through the PI3K/Akt/NF-κB pathway. Phytomedicine 125, 155312. doi:10.1016/j.phymed.2023.155312
Li, J., Fu, A., and Zhang, L. (2019). An overview of scoring functions used for protein-ligand interactions in molecular docking. Interdiscip. Sci. 11, 320–328. doi:10.1007/s12539-019-00327-w
Li, X. Y., Liu, J. Q., Wang, Y., Chen, Y., Hu, W. H., Lv, Y. X., et al. (2024). VNS improves VSMC metabolism and arteriogenesis in infarcted hearts through m/n-AChR-Akt-SDF-1α in adult male rats. J. Mol. Histol. 55, 51–67. doi:10.1007/s10735-023-10171-4
Liao, Y., Ding, Y., Yu, L., Xiang, C., and Yang, M. (2022). Exploring the mechanism of Alisma orientale for the treatment of pregnancy induced hypertension and potential hepato-nephrotoxicity by using network pharmacology, network toxicology, molecular docking and molecular dynamics simulation. Front. Pharmacol. 13, 1027112. doi:10.3389/fphar.2022.1027112
Lin, X., Han, T., Fan, Y., Wu, S., Wang, F., and Wang, C. (2020). Quercetin improves vascular endothelial function through promotion of autophagy in hypertensive rats. Life Sci. 258, 118106. doi:10.1016/j.lfs.2020.118106
Liu, B., Chen, L., Zhang, Y.-H., Wang, S., Zhang, Y., Huang, T., et al. (2017). Prediction and analysis of essential genes using the enrichments of gene ontology and KEGG pathways. Plos One 12, e0184129. doi:10.1371/journal.pone.0184129
Liu, H., Mohammed, S. A. D., Lu, F., Chen, P., Wang, Y., and Liu, S. (2022). Network pharmacology and molecular docking-based mechanism study to reveal antihypertensive effect of gedan jiangya decoction. Biomed. Res. Int. 2022, 1–17. doi:10.1155/2022/3353464
Liu, X., Liu, J., Fu, B., Chen, R., Jiang, J., Chen, H., et al. (2023). DCABM-TCM: a database of constituents absorbed into the blood and metabolites of traditional Chinese medicine. J. Chem. Inf. Model 63, 4948–4959. doi:10.1021/acs.jcim.3c00365
Luo, F., Wu, Y., Ding, Q., Yuan, Y., and Jia, W. (2021). Rs884225 polymorphism is associated with primary hypertension by compromising interaction between epithelial growth factor receptor (EGFR) and miR-214. J. Cell Mol. Med. 25, 3714–3723. doi:10.1111/jcmm.15976
Pereira, S. C., Parente, J. M., Belo, V. A., Mendes, A. S., Gonzaga, N. A., do Vale, G. T., et al. (2018). Quercetin decreases the activity of matrix metalloproteinase-2 and ameliorates vascular remodeling in renovascular hypertension. Atherosclerosis 270, 146–153. doi:10.1016/j.atherosclerosis.2018.01.031
Pinero, J., Ramirez-Anguita, J. M., Sauch-Pitarch, J., Ronzano, F., Centeno, E., Sanz, F., et al. (2020). The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 48, D845–D855. doi:10.1093/nar/gkz1021
Rahimi, K., Tucker, K. L., Sheppard, J. P., Stevens, R., Bosworth, H. B., Bove, A., et al. (2017). Self-monitoring of blood pressure in hypertension: a systematic review and individual patient data meta-analysis. PLOS Med. 14, e1002389. doi:10.1371/journal.pmed.1002389
Rasmussen, S. A., Hamosh, A., curators, O., Arnold, C., Bocchini, C., O‘Neill, M. J., et al. (2020). What's in a name? Issues to consider when naming Mendelian disorders. Genet. Med. 22, 1573–1575. doi:10.1038/s41436-020-0851-0
Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. (2014). TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform 6, 13. doi:10.1186/1758-2946-6-13
Sanajou, D., Ghorbani Haghjo, A., Argani, H., and Aslani, S. (2018). AGE-RAGE axis blockade in diabetic nephropathy: current status and future directions. Eur. J. Pharmacol. 833, 158–164. doi:10.1016/j.ejphar.2018.06.001
Shudong, M. (2022). Study on Traditional Chinese Medicine Master Lu Zhizheng’s experience in the treatment of hypertension, based on date mining. China: Chinese Academy of Traditional Chinese Medicine.
Stelzer, G., Rosen, N., Plaschkes, I., Zimmerman, S., Twik, M., Fishilevich, S., et al. (2016). The GeneCards suite: from gene data mining to disease genome sequence analyses. Curr. Protoc. Bioinforma. 54, 1. doi:10.1002/cpbi.5
Szklarczyk, D., Gable, A. L., Nastou, K. C., Lyon, D., Kirsch, R., Pyysalo, S., et al. (2021). The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 49, D605–D612. doi:10.1093/nar/gkaa1074
Tettey, C. O., Yang, I. J., and Shin, H. M. (2019). Vasodilatory effect of kaempferol-7-O-α-L-rhamnopyranoside via NO-cGMP-PKG signaling. Arch. Biochem. Biophys. 667, 1–5. doi:10.1016/j.abb.2019.04.001
Tian, Z., Zhang, S., Wang, H., Chen, Z., Sun, M., Sun, L., et al. (2020). Intervention of Uncaria and its components on liver lipid metabolism in spontaneously hypertensive rats. Front. Pharmacol. 11, 910. doi:10.3389/fphar.2020.00910
Trott, O., and Olson, A. J. (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 31, 455–461. doi:10.1002/jcc.21334
UniProt, C., Martin, M. J., Orchard, S., Magrane, M., Ahmad, S., Alpi, E., et al. (2023). UniProt: the universal protein knowledgebase in 2023. Nucleic Acids Res. 51, D523–D531. doi:10.1093/nar/gkac1052
Wang, D., Ali, F., Liu, H., Cheng, Y., Wu, M., Saleem, M. Z., et al. (2022). Quercetin inhibits angiotensin II-induced vascular smooth muscle cell proliferation and activation of JAK2/STAT3 pathway: a target based networking pharmacology approach. Front. Pharmacol. 13, 1002363. doi:10.3389/fphar.2022.1002363
Wang, X., Shen, Y., Wang, S., Li, S., Zhang, W., Liu, X., et al. (2017). PharmMapper 2017 update: a web server for potential drug target identification with a comprehensive target pharmacophore database. Nucleic Acids Res. 45, W356–W360. doi:10.1093/nar/gkx374
Wenlong, L. (2022). Analysis of medication patterns in the treatment of essential hypertension based on the theory of hepatic and renal homoeopathy. China: Heilongjiang University of Traditional Chinese Medicine.
Wishart, D. S., Knox, C., Guo, A. C., Shrivastava, S., Hassanali, M., Stothard, P., et al. (2006). DrugBank: a comprehensive resource for in silico drug discovery and exploration. Nucleic Acids Res. 34, D668–D672. doi:10.1093/nar/gkj067
Wolf, D., and Ley, K. (2019). Immunity and inflammation in atherosclerosis. Circ. Res. 124, 315–327. doi:10.1161/circresaha.118.313591
Wu, T., Hu, E., Xu, S., Chen, M., Guo, P., Dai, Z., et al. (2021). clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innov. (Camb) 2, 100141. doi:10.1016/j.xinn.2021.100141
Xiaotian, Y., Chun, W., Liping, D., and Zhiqing, C. (2021). CHENG Zhi-qing’s clinical experience in treating hypertension. Chin. J. Traditional Chin. Med. 36, 5942–5945.
Yichu, N. (2021). Data mining research on the prescription and medication pattern of Professor Shi Dazhuo in the treatment of hypertension. China: Beijing University of Chinese Medicine.
Zhai, Z., Tao, X., Alami, M. M., Shu, S., and Wang, X. (2021). Network pharmacology and molecular docking combined to analyze the molecular and pharmacological mechanism of pinellia ternata in the treatment of hypertension. Curr. Issues Mol. Biol. 43, 65–78. doi:10.3390/cimb43010006
Zhang, L., Zhang, Y., Wu, Y., Yu, J., Zhang, Y., Zeng, F., et al. (2019). Role of the balance of akt and MAPK pathways in the exercise-regulated phenotype switching in spontaneously hypertensive rats. Int. J. Mol. Sci. 20, 5690. doi:10.3390/ijms20225690
Zhao, J., Ling, L., Zhu, W., Ying, T., Yu, T., Sun, M., et al. (2023). M1/M2 re-polarization of kaempferol biomimetic NPs in anti-inflammatory therapy of atherosclerosis. J. Control Release 353, 1068–1083. doi:10.1016/j.jconrel.2022.12.041
Keywords: hypertension, traditional Chinese medicine, Uncaria rhynchophylla, Alisma plantago-aquatica L., network pharmacology, molecular docking, bioactive compounds, drug discovery
Citation: Yin T, Zhang H, Liu X, Wei D, Ren C, Cui L, Li Y, Wang L, Wang J, Zhao Z, Liu D, Wang L and Han X (2024) Elucidating the anti-hypertensive mechanisms of Uncaria rhynchophylla-Alisma plantago-aquatica L: an integrated network pharmacology, cluster analysis, and molecular docking approach. Front. Chem. 12:1356458. doi: 10.3389/fchem.2024.1356458
Received: 15 December 2023; Accepted: 19 February 2024;
Published: 29 February 2024.
Edited by:
Jorddy Neves Cruz, Federal University of Pará, BrazilReviewed by:
Tianlong Liu, The 940th Hospital of Joint Logistic Support Force of PLA, ChinaJian Chen, Shanghai University of Traditional Chinese Medicine, China
Copyright © 2024 Yin, Zhang, Liu, Wei, Ren, Cui, Li, Wang, Wang, Zhao, Liu, Wang and Han. 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: Xuejie Han, xuejieh@126.com; Liying Wang, coldmoon_ly@163.com; Dasheng Liu, zzdasheng@163.com
†These authors have contributed equally to this work