Evidence-Based Complementary and Alternative Medicine

Evidence-Based Complementary and Alternative Medicine / 2020 / Article
Special Issue

Complementary and Alternative Medicine for Ameliorating Bone and Joint Disorders

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 7151634 | https://doi.org/10.1155/2020/7151634

Wei Liu, Yihua Fan, Chunying Tian, Yue Jin, Shaopeng Du, Ping Zeng, Aihua Wang, "Deciphering the Molecular Targets and Mechanisms of HGWD in the Treatment of Rheumatoid Arthritis via Network Pharmacology and Molecular Docking", Evidence-Based Complementary and Alternative Medicine, vol. 2020, Article ID 7151634, 13 pages, 2020. https://doi.org/10.1155/2020/7151634

Deciphering the Molecular Targets and Mechanisms of HGWD in the Treatment of Rheumatoid Arthritis via Network Pharmacology and Molecular Docking

Academic Editor: Arham Shabbir
Received07 Jul 2020
Accepted18 Aug 2020
Published26 Aug 2020

Abstract

Background. Huangqi Guizhi Wuwu Decoction (HGWD) has been applied in the treatment of joint pain for more than 1000 years in China. Currently, most physicians use HGWD to treat rheumatoid arthritis (RA), and it has proved to have high efficacy. Therefore, it is necessary to explore the potential mechanism of action of HGWD in RA treatment based on network pharmacology and molecular docking methods. Methods. The active compounds of HGWD were collected, and their targets were identified from the Traditional Chinese Medicine Systems Pharmacology Database (TCMSP) and DrugBank database, respectively. The RA-related targets were retrieved by analyzing the differentially expressed genes between RA patients and healthy individuals. Subsequently, the compound-target network of HGWD was constructed and visualized through Cytoscape 3.8.0 software. Protein-protein interaction (PPI) network was constructed to explore the potential mechanisms of HGWD on RA using the plugin BisoGenet of Cytoscape 3.8.0 software. Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed in R software (Bioconductor, clusterProfiler). Afterward, molecular docking was used to analyze the binding force of the top 10 active compounds with target proteins of VCAM1, CTNNB1, and JUN. Results. Cumulatively, 790 active compounds and 1006 targets of HGWD were identified. A total of 4570 differentially expressed genes of RA with a value <0.05 and  > 0.5 were collected. Moreover, 739 GO entries of HGWD on RA were identified, and 79 pathways were screened based on GO and KEGG analysis. The core target gene of HGWD in RA treatment was JUN. Other key target genes included FOS, CCND1, IL6, E2F2, and ICAM1. It was confirmed that the TNF signaling pathway and IL-17 signaling pathway are important pathways of HGWD in the treatment of RA. The molecular docking results revealed that the top 10 active compounds of HGWD had a strong binding to the target proteins of VCAM1, CTNNB1, and JUN. Conclusion. HGWD has important active compounds such as quercetin, kaempferol, and beta-sitosterol, which exert its therapeutic effect on multiple targets and multiple pathways.

1. Introduction

Rheumatoid arthritis (RA), a common chronic systemic autoimmune disease, is characterized by synovial hyperproliferation and inflammatory/immune cell infiltration [1, 2]. Typical clinical symptoms of RA include tender and swelling of the joints, accompanied by morning stiffness of the affected joints. In severe cases, large joints can be injured leading to joint deformity and loss of function. The current global incidence of RA is 0.24%, and it is expected to increase [36]. Besides, RA is ranked 74th as a social burden and 42nd as a disability [36]. Current international guidelines for the management of RA recommend the use of disease-modifying antirheumatic drugs (DMARDs), with methotrexate (MTX) as the first-line drug [7, 8]. However, MTX is not an ideal therapeutic agent. It has damaging side effects on the neuronal, gastrointestinal, and immune systems [9]. Therefore, it is imperative to explore safe and effective clinical treatment for RA.

The traditional Chinese medicine (TCM) is based on the theory of syndrome differentiation and has long been established as an effective treatment of RA. Huangqi Guizhi Wuwu Decoction (HGWD), a classical prescription described in Jingui Yaolue, has been used over 1000 years in China to treat joint pain. HGWD was formulated by Zhang Zhongjing during the Han dynasty. It is a mixture of five Chinese medicines including Astragalus membranaceus (Huangqi), Cassia Twig (Guizhi), Radix Paeoniae Alba (Baishao), ginger (Shengjiang), and Chinese date (Dazao). HGWD is considered to have the ability to reinforce Qi and nourish blood, warm and smoothen meridian, and dredge arthritic pain traditionally [10]. In a study of rats with RA, HGWD administration reduced joint inflammation, synovial hyperplasia, and cartilage damage [11]. Previous studies have revealed that HGWD improves the clinical symptoms and signs, as well as laboratory indices, in RA patients [12]. Moreover, a systematic review and meta-analysis showed that when HGWD was combined with Western medicine therapy to treat RA, better efficacy, improved morning stiffness, reduced C-reactive protein, and rheumatoid factor content were achieved compared to Western medicine only [13]. However, the mechanisms of HGWD in RA treatment are not clear, a major limiting factor for its extensive application.

Network pharmacology is partly bioinformatics and was first proposed by Hopkins [14]. It has been successfully used to study complex TCM formulations. This is because it not only combines system network analysis and pharmacology but is also based on the connotation of holistic theory, multicomponents, multitargets, and multipathways of Chinese medicine [15, 16]. Network pharmacology can elucidate the mechanisms of HGWD in the RA treatment at the molecular level via compound-compound network, compound-target network, and target-disease network.

In the present study, network pharmacology was used to establish a compound-target-disease network for exploring the potential HGWD mechanism of action in RA treatment. This study provides a reference for future pharmacological studies and clinical applications. The flow diagram of the network is shown in Figure 1.

2. Methods

2.1. Collection of Active HGWD Compound Information

Information on the HGWD compounds was retrieved from the Traditional Chinese Medicine Systems Pharmacology Database (TCMSP, http://tcmspw.com/tcmsp.php) [17]. TCMSP is a unique system pharmacology database of Chinese herbal medicines with data on absorption, distribution, metabolism, and excretion (ADME)-related parameters of herbal ingredients as well as the relationships between diseases, targets, ingredients, and drugs [18]. The active compounds of HGWD were primarily screened based on oral bioavailability (OB) and drug-like (DL) properties, the two important indicators for bioinformatics evaluation of ADME characteristics [19]. The OB is a major pharmacokinetic parameter of oral drugs and is the proportion of oral drug dose in the systemic circulation [20]. DL properties are the physical and chemical properties that qualitatively evaluate whether a compound is similar to existing drugs [21]. Subsequently, the compounds with OB ≥ 30% and DL ≥ 0.18 were chosen as the candidate compounds for further analysis.

2.2. Identification of HGWD Potential Targets

Target identification is an important aspect of drug exploration [22]. The target candidate compounds were imported into the DrugBank database (https://www.drugbank.ca/) to identify the corresponding potential targets of HGWD [23].

2.3. Known RA-Related Targets

The differentially expressed genes in RA were retrieved from the GEO database (https://www.ncbi.nlm.nih.gov; series: GSE21959; samples: there were 36 samples, 18 for healthy individuals and 18 for those with rheumatoid arthritis). The genes with a value <0.05 and  > 0.5 were regarded as differentially expressed genes and RA-related targets.

2.4. Network Construction

The compound-target network of HGWD was constructed and visualized via Cytoscape 3.8.0 software [24]. Protein-protein interaction (PPI) networks were constructed to explore the potential mechanisms of HGWD on RA using the plugin BisoGenet of Cytoscape 3.8.0 software [25]. Afterward, the topological importance of each node by calculating degree centrality (DC), betweenness centrality (BC), closeness centrality (CC), eigenvector centrality (EC), local average connectivity-based method (LAC), and network centrality (NC) was evaluated with a Cytoscape plugin CytoNCA. Their definitions and computational formulas have been reported and are the topological importance representative of each node [26].

2.5. Gene Ontology and Pathway Enrichment Analysis

Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) are important methods that describe the features of candidate targets. The two were performed in R software (Bioconductor, clusterProfiler) with the standard value cutoff of 0.05 and the q value of 0.05 [27, 28]. The GO analysis was applied for target protein analysis, and the top 20 functional categories in biological process (BP), cellular component (CC), and molecular function (MF) were chosen. Based on the targets of HGWD in RA treatment, KEGG enrichment analysis was carried out. The top 20 KEGG pathways were selected for plotting a histogram. Meanwhile, the network diagram of the gene pathway was drawn.

2.6. Molecular Docking of the Main Active Constituents of HGWD and Core Proteins

The 3D protein structure of the three core proteins corresponding to the core targets, VCAM1, CTNNB1, and JUN, was downloaded from the UniProt database. Subsequently, the structure of the active ingredient in HGWD (the top 10 places in the number of targets) was downloaded from the PubChem database and saved in the PDB format. Using PyMOL software, the three proteins were virtually dehydrated and hydrogenated, the original ligands were extracted in each protein, and then were stored separately. AutoDockTools 1.5.6 was utilized to convert compounds, ligands, and proteins into the “pdbqt” format and to define if the location of each protein or its ligands was the active pocket of the protein. Finally, Vina 1.5.6 was run to assess molecular docking. At a binding energy value <0, the molecular proteins were considered to spontaneously bind and interact with each other. Accordingly, the lower the energy is, the more stable the molecular conformation is.

3. Results

3.1. Target Screening of HGWD and RA

A total of 790 compounds of the five herb medicines in HGWD were retrieved from the TCMSP database. This included 87 compounds in Huangqi, 220 in Guizhi, 85 in Baishao, 265 in Shengjiang, and 133 in Dazao. Among them, 74 compounds passed OB ≥ 30% and DL ≥ 0.18 filtering. Specifically, the numbers of candidate compounds in Huangqi, Guizhi, Baishao, Shengjiang, and Dazao were 20, 7, 13, 5, and 29, respectively. The candidate compounds in HGWD used for further analysis are shown in Table 1. The DrugBank database retrieval predicted a total of 1006 potential targets. The potential targets linked to Huangqi, Guizhi, Baishao, Shengjiang, and Dazao were 405, 57, 104, 60, and 380, respectively. Moreover, a total of 4570 differentially expressed genes in RA were collected from the GEO database. As shown in Figure 2, a volcano plot was drawn to show the distribution of the differentially expressed genes. The genes are represented by the red and green dots in the plot. We compared the target genes regulated by the active compounds in HGWD, and different genes in RA were compared, obtaining 49 common target genes. These 49 target genes were found to be regulated by 28 active compounds.


Molecule IDNameOBDLSource

MOL000211Mairin55.380.78Huangqi, Baishao, Dazao
MOL000239Jaranol50.830.29Huangqi
MOL000296Hederagenin36.910.75Huangqi
MOL000033(3S,8S,9S,10R,13R,14S,17R)-10,13-Dimethyl-17-[(2R,5S)-5-propan-2-yloctan-2-yl]-2,3,4,7,8,9,11,12,14,15,16,17-dodecahydro-1H-cyclopenta[a]phenanthren-3-ol36.230.78Huangqi
MOL000354Isorhamnetin49.60.31Huangqi
MOL0003713,9-Di-O-methylnissolin53.740.48Huangqi
MOL0003745′-Hydroxyiso-muronulatol-2′,5′-di-O-glucoside41.720.69Huangqi
MOL0003787-O-Methylisomucronulatol74.690.3Huangqi
MOL0003799,10-Dimethoxypterocarpan-3-O-β-D-glucoside36.740.92Huangqi
MOL000380(6aR,11aR)-9,10-Dimethoxy-6a,11a-dihydro-6H-benzofuro[3,2-c]chromen-3-ol64.260.42Huangqi
MOL000387Bifendate31.10.67Huangqi
MOL000392Formononetin69.670.21Huangqi
MOL000398Isoflavanone109.990.3Huangqi
MOL000417Calycosin47.750.24Huangqi
MOL000422Kaempferol41.880.24Huangqi, Baishao
MOL000433FA68.960.71Huangqi
MOL000438(3R)-3-(2-Hydroxy-3,4-dimethoxyphenyl)chroman-7-ol67.670.26Huangqi
MOL000439Isomucronulatol-7,2′-di-O-glucosiole49.280.62Huangqi
MOL0004421,7-Dihydroxy-3,9-dimethoxy pterocarpene39.050.48Huangqi
MOL000098Quercetin46.430.28Huangqi, Dazao
MOL001736(-)-Taxifolin60.510.27Guizhi
MOL000358Beta-sitosterol36.910.75Guizhi、Baishao、Shengjiang、Dazao
MOL000359Sitosterol36.910.75Guizhi, Baishao
MOL000492(+)-Catechin54.830.24Guizhi, Baishao, Dazao
MOL000073ent-Epicatechin48.960.24Guizhi
MOL004576Taxifolin57.840.27Guizhi
MOL011169Peroxyergosterol44.390.82Guizhi
MOL001928Albiflorin_qt66.640.33Baishao
MOL001918Paeoniflorgenone87.590.37Baishao
MOL00191011alpha,12alpha-Epoxy-3beta-23-dihydroxy-30-norolean-20-en-28,12beta-olide64.770.38Baishao
MOL001925Paeoniflorin_qt68.180.4Baishao
MOL001919(3S,5R,8R,9R,10S,14S)-3,17-Dihydroxy-4,4,8,10,14-pentamethyl-2,3,5,6,7,9-hexahydro-1H-cyclopenta[a]phenanthrene-15,16-dione43.560.53Baishao
MOL001930Benzoyl paeoniflorin31.270.75Baishao
MOL001924Paeoniflorin53.870.79Baishao
MOL001921Lactiflorin49.120.8Baishao
MOL0061296-Methylgingediacetate248.730.32Shengjiang
MOL000449Stigmasterol43.830.76Shengjiang, Dazao
MOL001771Poriferast-5-en-3beta-ol36.910.75Shengjiang
MOL008698Dihydrocapsaicin47.070.19Shengjiang
MOL012921Stepharine31.550.33Dazao
MOL012940Spiradine A113.520.61Dazao
MOL012946Zizyphus saponin I_qt32.690.62Dazao
MOL012961Jujuboside A_qt36.670.62Dazao
MOL012976Coumestrol32.490.34Dazao
MOL012980Daechuine S646.480.79Dazao
MOL012981Daechuine S744.820.83Dazao
MOL012986Jujubasaponin V_qt36.990.63Dazao
MOL012989Jujuboside C_qt40.260.62Dazao
MOL012992Mauritine D89.130.45Dazao
MOL001454Berberine36.860.78Dazao
MOL001522(S)-Coclaurine42.350.24Dazao
MOL003410Ziziphin_qt66.950.62Dazao
MOL004350Ruvoside_qt36.120.76Dazao
MOL005360Malkangunin57.710.63Dazao
MOL000627Stepholidine33.110.54Dazao
MOL007213Nuciferin34.430.4Dazao
MOL000783Protoporphyrin30.860.56Dazao
MOL000787Fumarine59.260.83Dazao
MOL00803421302-79-473.520.77Dazao
MOL008647Moupinamide86.710.26Dazao
MOL002773Beta-carotene37.180.58Dazao
MOL000096(-)-Catechin49.680.24Dazao
MOL013357(3S,6R,8S,9S,10R,13R,14S,17R)-17-[(1R,4R)-4-Ethyl-1,5-dimethylhexyl]-10,13-dimethyl-2,3,6,7,8,9,11,12,14,15,16,17-dodecahydro-1H-cyclopenta[a]phenanthrene-3,6-diol34.370.78Dazao

Note. OB, oral bioavailability; DL, drug-likeness.
3.2. Compound-Target Network Analysis

The compound-target network of HGWD was established with the collected compounds and their targets as shown in Figure 3. The network contains 77 nodes (28 compounds in HGWD and 49 compound targets) and 130 edges elucidating the compound-target interactions. Quercetin, kaempferol, and beta-sitosterol acted on 33, 14, and 8 targets, respectively. The OB of quercetin, kaempferol, and beta-sitosterol was 43.43, 41.88, and 36.91%, respectively. Therefore, they are potential key active compounds due to their relative positioning in the network. Besides, it has been revealed that Guizhi, Baishao, and Dazao share a common component (+)-catechin (ID: MOL000492); Shengjiang and Dazao have the same ingredient stigmasterol (ID: MOL000449); Guizhi, Baishao, Shengjiang, and Dazao share a common ingredient beta-sitosterol (ID: MOL000358); Huangqi and Baishao share the same constituent kaempferol (ID: MOL000422); and Huangqi and Dazao share the same component quercetin (ID: MOL000098).

3.3. The Candidate Targets for HGWD against RA

To elucidate the mechanism by which HGWD ameliorates RA, the potential target network was merged with the RA-related target network to form a core PPI network of 2284 nodes and 53,119 edges (Figure 4(a)). The previous research showed that the median degree of all nodes was more than 74 degrees which were identified as significant targets [29]. Accordingly, a network of the significant targets for HGWD against RA with 426 nodes and 17,112 edges with a DC ≥ 74 was constructed (Figure 4(b)).

During the second screening, since the number of genes was limited, only BC average value was used. The candidate targets were further identified, and 111 targets with BC ≥ 348.07 (BC average value) were chosen (Figure 4(c)). Eventually, 111 HGWD target genes against RA were identified. VCAM1 (degree: 452), CTNNB1 (degree: 390), and JUN (degree: 274) might be the most important among 111 target genes for HGWD against RA.

3.4. GO and KEGG Pathway Enrichment Analysis

To further confirm the biological responses from RA treatment with HGWD, GO analysis of the 49 RA-related potential therapeutic target genes was performed based on BP, CC, and MF. The analysis results revealed a total of 739 entries. In BP enrichment analysis, 666 entries were obtained including response to antibiotic, response to alcohol, and response to steroid hormone. In CC enrichment analysis, 34 entries involved in membrane raft, membrane microdomain, membrane region, etc. were obtained. The MF enrichment analysis revealed 39 entries, including kinase regulator activity, protein kinase regulator activity, and serine hydrolase activity. The top 20 terms are shown in Figure 5.

To further show the biological processes of these targets, the KEGG pathway analysis was performed. The analysis results showed that these targets were enriched in 79 pathways with an adjusted value <0.05. The top 20 KEGG pathways’ enrichment analysis is shown in Figure 6. Among these pathways, the top three were fluid shear stress and atherosclerosis, TNF signaling pathway, and IL-17 signaling pathway based on the number of the pathway target genes.

3.5. The Gene-Pathway Network

The gene-pathway network analysis was constructed according to the significantly enriched pathways and genes which regulate these pathways as shown in Figure 7. The diagram shows the network relationship between the top 20 pathways and their regulated target genes. According to the network analysis, JUN had the highest volume, hence the core target gene. Additionally, other genes had a relatively high volume including FOS, CCND1, IL6, E2F2, and ICAM1. These genes are potential key target genes involved in the HGWD treatment of RA.

3.6. The Molecular Docking of the Main Active Constituents of HGWD and Core Proteins

Normally, high connectivity compounds are associated with more targets. According to the number of related targets, the top 10 active compounds in HGWD with high connectivity in the compound-target network were selected for molecular docking (Table 2, binding energy of the main potential active ingredients in HGWD). Based on the binding energy value, the lower the binding energy value, the stronger the binding to the target protein. Among them, stigmasterol and stepholidine had the strongest binding force with VCAM1 (PDB: 1ij9), beta-carotene and quercetin had the strongest binding force with CTNNB1 (PDB: 1jdh), and beta-carotene and stigmasterol had the strongest binding force with JUN (PDB: 2g01). The corresponding molecular docking diagram of stepholidine, stigmasterol, and VCAM1, beta-carotene, quercetin and CTNNB1, and beta-carotene, stigmasterol, and JUN is illustrated in Figure 8.


CompoundMolecular formulaBinding energy (VCAM1)Binding energy (CTNNB1)Binding energy (JUN)

QuercetinC15H10O7−7.1−7.3−7.3
KaempferolC15H10O6−7.1−6.5−7.4
Beta-sitosterolC29H50O−7.3−7.1−7.8
FormononetinC16H12O4−7−6.3−6.9
7-O-MethylisomucronulatolC18H20O5−6.7−6−6.9
IsorhamnetinC16H12O7−6.9−6.5−7.5
Beta-caroteneC40H56−6.5−7.9−8.1
StepholidineC19H21NO4−7.5−6.7−6.9
(S)-CoclaurineC17H19NO3−7.1−6.3−7.1
StigmasterolC29H48O−7.6−7−8.1

4. Discussion

TCM, characterized by multicompound and multitarget medicines, cures diseases via multiple targets, multiple pathways, and multiple links. Due to the complex chemical ingredients of TCM, conventional research methods such as clinical and pharmacological research are not capable of fully elucidating the mechanism of action of TCM. Fortunately, network pharmacology provides a solution to this challenge since it is suitable for multicompound and multitarget research. In the present study, the mechanism of action of HGWD on RA was explored via network pharmacology methods. This provided a clear direction for further basic and clinical research.

Modern pharmacology has shown that HGWD has a specific therapeutic effect on RA, where it elicits obvious anti-inflammatory and analgesic effects [30]. Specifically, Huangqi has an anti-inflammatory effect and improves the immunity [31]. The volatile oil of Guizhi has good anti-inflammatory, immunological, and chondrocyte proliferation effects [32]. The total glycosides, the main components of Baishao, have anti-inflammatory, analgesic, and autoimmunological effects [33]. Shengjiang resists oxidation, scavenges free radicals, and relieves pain [34]. The polysaccharide of Dazao can significantly inhibit proinflammatory cytokines, such as IL-6 and TNF, with anti-inflammatory activity [35]. It is suggested that this prescription has anti-inflammatory, antioxidation, analgesic, and autoimmune reaction pharmacological effects, thus providing specific pharmacological basis for its clinical function in the treatment of RA.

In this study, the compound-target network of HGWD was constructed using 28 compounds and their 49 targets. The network diagram demonstrated that most of the HGWD compounds affected multiple targets. For instance, quercetin, kaempferol, and beta-sitosterol acted on 33, 14, and 8 targets, respectively. Therefore, they were the likely key active compounds for HGWD. Quercetin is a common active component of Huangqi and Dazao. Kaempferol is a common active component of Huangqi and Baishao. Beta-sitosterol is a common active component of Guizhi, Baishao, Shengjiang, and Dazao. These drugs have been mainly used together in Chinese medical history. It indicates that the compatibility between these drugs has a synergistic effect and increases their efficacy.

Quercetin is a flavonol that has unique therapeutic biological properties including anticarcinogenic, anti-inflammatory, antioxidant, antiviral, and psychostimulant activities [36, 37]. Kaempferol, a representative flavonoid, has been known to exert a range of pharmacological actions including the mediation of antioxidant, antimicrobial, and anti-inflammatory effects [38, 39]. Kaempferol is a potent immunosuppressant that reduces the harmful immune responses including autoimmunity and chronic inflammation [40]. Beta-sitosterol, a main constituent of plants and vegetables, has versatile activities that impact cell activities including anti-inflammatory effect [41, 42]. In this study, quercetin, kaempferol, and beta-sitosterol regulated most RA-related targets and exhibited anti-inflammatory effects. Besides, they have a high oral bioavailability. Therefore, they were identified as the representative compounds of HGWD.

The PPI networks of HGWD targets and RA-related targets were constructed and merged to identify the candidate targets of HGWD against RA. To accurately obtain the targets, two parameters including DC and BC were set to construct a new network. Through PPI analysis, VCAM1, CTNNB1, and JUN were established as the important targets in RA treatment. VCAM-1 is a glycoprotein expressed in vascular endothelial cells, whose serum is positively correlated with RA [43]. Studies have found that [44] TNF-α upregulates the expression of VCAM-1 in the endothelial cell membrane at the small joint synovium. VCAM-1 binds to α4β1 on the surface of free immune cells in the blood vessels, inducing the migration of immune cells to inflammatory joints thus expanding the inflammatory response cascade and aggravating small joint injury. Proliferation and invasion of fibroblast synovial cells are important mechanisms that result in the thickened synovial lining, increased secretion of inflammatory cytokines, and cartilage and bone injury [45]. Studies have shown that inhibition of CTNNB1 transcription reduces the proliferation of fibroblast synovial cells and the levels of proinflammatory cytokines such as interleukin-6 (IL-6), IL-10, and TNF-α [46]. Elsewhere, it has been shown that CD44 can activate the transcription factor AP-1 (the whole is a protein) expressed by the gene JUN, thus promoting the activation of T cells and aggravating RA synovitis [47].

The targets of HGWD against RA were enriched in BP, CC, and MF through GO analysis. Results suggested that HGWD mainly regulated response to antibiotics, response to steroid hormones, and response to radiation. Furthermore, it was shown to affect some cellular components and molecular functions including membrane raft, membrane microdomain, membrane region, and kinase regulator activity. In the present study, 79 KEGG pathways including TNF signaling pathway and IL-17 signaling pathway were significantly enriched. IL-17 and TNF-α are classic proinflammatory cytokines. IL-17 is mainly secreted by helper T-cells type 17 and has strong proinflammatory effects. Analysis of synovial fluid in patients with RA found that IL-17 content was significantly increased and positively correlated with Disease Activity Score-28 (DAS28) [48]. This was a confirmation that IL-17 is involved in the occurrence of RA. Fischer et al. [49] found that, after IL-17 treatment of fibroblast synovial cells, the expression of proinflammatory cytokines such as IL-6, IL-8, and granulocyte-macrophage colony-stimulating factor (GM-CSF) increases, aggravating the inflammatory damage of fibroblast synovial cells. Besides, IL-17 can also induce joint bone and cartilage damages [50]. It has been revealed that IL-17 upregulates the levels of nuclear factor-kappa B (NF-κB) receptor activator ligand RANKL. Consequently, this increases the content of matrix metalloproteinases, possibly promoting cartilage degradation, inhibiting chondroprogenitor cell formation, and stimulating osteoclast bone damaging [5153].

TNF-α plays an important role in the development of RA. On the one hand, it can bind to the TNFR1 receptor on fibroblast synovial cells, promoting the release of inflammatory cytokines such as IL-1, IL-6, and IL-8 and aggravating the damage of the articular cartilage and bone [54]. On the other hand, TNF-α can stimulate immune cells in blood to enter the joint through VCAM-1, aggravating the joint injury [44]. Besides, the other pathways in the first 20 pathways have not been reported to be related to RA, providing new ideas and clues for future research.

The gene-pathway network was constructed to explore the main target genes of HGWD against RA. The results revealed that JUN, FOS, CCND1, IL6, E2F2, and ICAM1 are important target genes. A study found that the ERK-JNK-P38 signaling pathways in autoimmune disease are activated, leading to high levels of downstream JUN and Fos protein. Subsequently, the two combine to form dimer transcription factor-activating protein 1 (AP-1), which is involved in the occurrence and development of RA. It regulates the transformation of initial T cells into effector T cells, hence regulating the immune response and inflammatory process [5558]. Inflammatory response runs throughout RA. The proinflammatory cytokines such as TNF-α, IL-6, and IL-17 play a crucial role in synovial fibroblast inflammation. Studies have demonstrated that TNF-α upregulates the E2F2 expression by activating the nuclear transcription factor NF-κB pathway. This enhances synovial cell proliferation and synovial tissue thickening leading to joint injury [59]. The serum level of intercellular adhesion molecule-1 (ICAM-1) is high in RA patients. It has been found that ICAM-1 induces inflammatory cell aggregation and promotes synovial cell inflammation once it combines with ligand lymphocyte function antigen-1 (LFA-1) [60]. Through the molecular docking analysis of the effects of the chief active components in HGWD and the core target proteins of RA, the molecular mechanism of the treatment of RA by the active components of TCM can be predicted. This is of great reference significance for subsequent research and development of targeted drugs. The results showed that the main compounds in HGWD have a strong binding force to the core proteins VCAM1, CTNNB1, and JUN. Among them, stigmasterol and stepholidine closely link to VCAM1 through hydrogen bonds and hydrophobic forces. Beta-carotene and quercetin closely link to CTNNB1 through hydrogen bonds and hydrophobic interactions. Stigmasterol closely links to JUN through hydrogen bonds and hydrophobic interactions. Beta-carotene closely links to JUN only through hydrophobic interactions. Based on this, it can be concluded that stigmasterol, stepholidine, beta-carotene, and quercetin are the key compounds in the treatment of RA.

5. Conclusion

Our study systematically elucidated the mechanisms of action and molecular targets of HGWD against RA via the network pharmacology approach. Quercetin, kaempferol, and beta-sitosterol regulated most of the targets related to RA. The HGWD can regulate gene function through their related pathways including TNF signaling pathway and IL-17 signaling pathway. The key target genes in the gene-pathway network of HGWD against RA were JUN, FOS, CCND1, IL6, E2F2, and ICAM1. Furthermore, according to molecular docking analysis, important compounds such as stepholidine, stigmasterol, beta-carotene, quercetin, and the core protein CTNNB1, VCAM1, and JUN all have good binding ability. The network pharmacology is a promising suitable approach for the study of TCM formulations.

Abbreviations

RA:Rheumatoid arthritis
DMARDs:Disease-modifying antirheumatic drugs
MTX:Methotrexate
TCM:Traditional Chinese medicine
HGWD:Huangqi Guizhi Wuwu Decoction
TCMSP:Traditional Chinese Medicine Systems Pharmacology Database
ADME:Absorption, distribution, metabolism, and excretion
OB:Oral bioavailability
DL:Drug-like properties
PPI:Protein-protein interaction
DC:Degree centrality
BC:Betweenness centrality
CC:Closeness centrality
EC:Eigenvector centrality
LAC:Local average connectivity-based method
NC:Network centrality
GO:Gene ontology
KEGG:Kyoto Encyclopedia of Genes and Genomes
BP:Biological process
CC:Cellular component
MF:Molecular function
IL-6:Interleukin-6
DAS28:Disease Activity Score-28
GM-CSF:Granulocyte-macrophage colony-stimulating factor
NF-κB:Nuclear factor-kappa B
AP-1:Activating protein-1
ICAM-1:Intercellular adhesion molecule-1
LFA-1:Lymphocyte function antigen-1.

Data Availability

The data supporting the findings of this study are available from the corresponding author upon request.

Disclosure

Wei Liu and Yihua Fan are the co-first authors.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Wei Liu and Yihua Fan contributed equally to this work.

Acknowledgments

The authors would like to thank FreeScience (http://www.home-for-researchers.com) for the help with the English language. This research was financially supported by the National Natural Science Foundation of China (no. 81673927), the National Key Research and Development Project (no. 2018YFC1705201), the Twelfth Five-Year Plan Key Disciplines of Traditional Chinese Medicine, organized by the State Administration of Traditional Chinese Medicine of China (no. 2018ZDXK001), the Tianjin Science and Technology Commission Program (2019ZD12), and Tianjin Key Specialty Program (20200602–1).

References

  1. E. Zampeli, P. G. Vlachoyiannopoulos, and A. G. Tzioufas, “Treatment of rheumatoid arthritis: unraveling the conundrum,” Journal of Autoimmunity, vol. 65, pp. 1–18, 2015. View at: Publisher Site | Google Scholar
  2. Q. Guo, Y. Wang, D. Xu, J. Nossent, N. J. Pavlos, and J. Xu, “Rheumatoid arthritis: pathological mechanisms and modern pharmacologic therapies,” Bone Research, vol. 6, no. 1, p. 15, 2018. View at: Publisher Site | Google Scholar
  3. M. Cross, E. Smith, D. Hoy et al., “The global burden of rheumatoid arthritis: estimates from the global burden of disease 2010 study,” Annals of the Rheumatic Diseases, vol. 73, no. 7, pp. 1316–1322, 2014. View at: Publisher Site | Google Scholar
  4. J. S. Smolen, D. Aletaha, M. Koeller, M. H. Weisman, and P. Emery, “New therapies for treatment of rheumatoid arthritis,” The Lancet, vol. 370, no. 9602, pp. 1861–1874, 2007. View at: Publisher Site | Google Scholar
  5. D. Symmons, G. Turner, R. Webb et al., “The prevalence of rheumatoid arthritis in the United Kingdom: new estimates for a new century,” Rheumatology, vol. 41, no. 7, p. 793, 2002. View at: Publisher Site | Google Scholar
  6. S. E. Gabriel, C. S. Crowson, and M. O’Fallon, “The epidemiology of rheumatoid arthritis in rochester, minnesota, 1955–1985,” Arthritis & Rheumatism, vol. 42, no. 3, pp. 415–420, 1999. View at: Publisher Site | Google Scholar
  7. J. S. Smolen, R. Landewe, J. Bijlsma et al., “Eular recommendations for the management of rheumatoid arthritis with synthetic and biological disease-modifying antirheumatic drugs: 2016 update,” Annals of the Rheumatic Diseases, vol. 76, no. 6, pp. 960–977, 2017. View at: Publisher Site | Google Scholar
  8. A. Mian, F. Ibrahim, and D. L. Scott, “A systematic review of guidelines for managing rheumatoid arthritis,” BMC Rheumatology, vol. 3, p. 42, 2019. View at: Publisher Site | Google Scholar
  9. V. C. Romão, A. Lima, M. Bernardes, H. Canhão, and J. E. Fonseca, “Three decades of low-dose methotrexate in rheumatoid arthritis: can we predict toxicity?” Immunologic Research, vol. 60, no. 2-3, pp. 289–310, 2014. View at: Publisher Site | Google Scholar
  10. B. Chen, P. W. Yuan, W. L. Kang et al., “Application of HGWD in orthopedics,” Chinese Journal of Traditional Medical Traumatology, vol. 23, pp. 71–74, 2015. View at: Google Scholar
  11. J. W. Liu, Y. H. Wang, Y. Y. Li et al., “Effect of HGWD on JAK-STAT signal pathway in rheumatoid arthritis CIA model rats,” Shizhen Guoyi Guoyao, vol. 4, pp. 811–814, 2019. View at: Google Scholar
  12. Q. W. Yang, X. L. Li, Q. C. Huang et al., “Therapeutic effect of HGWD on rheumatoid arthritis,” World Traditional Chinese Medicine, vol. 4, pp. 861–864, 2018. View at: Google Scholar
  13. Z. Jiang, F. Tang, W. Ma et al., “Meta analysis of the clinical effect of HGWD on rheumatoid arthritis,” Rheumatism and Arthritis, vol. 8, no. 5, pp. 33–38, 2019, in Chinese. View at: Google Scholar
  14. A. L. Hopkins, “Network pharmacology: the next paradigm in drug discovery,” Nature Chemical Biology, vol. 4, no. 11, pp. 682–690, 2008. View at: Publisher Site | Google Scholar
  15. J.-S. You, C.-Y. Li, W. Chen et al., “A network pharmacology-based study on alzheimer disease prevention and treatment of qiong yu gao,” BioData Mining, vol. 13, no. 1, 2020. View at: Publisher Site | Google Scholar
  16. Y. Jiang, N. Liu, S. Zhu, X. Hu, D. Chang, and J. Liu, “Elucidation of the mechanisms and molecular targets of yiqi shexue formula for treatment of primary immune thrombocytopenia based on network pharmacology,” Frontiers in Pharmacology, vol. 10, p. 1136, 2019. View at: Publisher Site | Google Scholar
  17. 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, p. 13, 2014. View at: Publisher Site | Google Scholar
  18. W.-Y. Lee, C.-Y. Lee, Y.-S. Kim, and C.-E. Kim, “The methodological trends of traditional herbal medicine employing network pharmacology,” Biomolecules, vol. 9, no. 8, p. 362, 2019. View at: Publisher Site | Google Scholar
  19. B. Gong, Y. Kao, C. Zhang, F. Sun, and H. Zhao, “Systematic investigation of scutellariae barbatae herba for treating hepatocellular carcinoma based on network pharmacology,” Evidence-Based Complementary and Alternative Medicine, vol. 2018, Article ID 4365739, 12 pages, 2018. View at: Publisher Site | Google Scholar
  20. X. Xu, W. Zhang, C. Huang et al., “A novel chemometric method for the prediction of human oral bioavailability,” International Journal of Molecular Sciences, vol. 13, no. 6, pp. 6964–6982, 2012. View at: Publisher Site | Google Scholar
  21. W. Wang, T. Liu, L. Yang et al., “Study on the multi-targets mechanism of triphala on cardio-cerebral vascular diseases based on network pharmacology,” Biomedicine & Pharmacotherapy, vol. 116, Article ID 108994, 2019. View at: Publisher Site | Google Scholar
  22. M. Xu, J. Shi, Z. Min, H. Zhu, and W. Sun, “A network pharmacology approach to uncover the molecular mechanisms of herbal formula kang-Bai-ling for treatment of vitiligo,” Evidence-Based Complementary and Alternative Medicine, vol. 2019, Article ID 3053458, 11 pages, 2019. View at: Publisher Site | Google Scholar
  23. V. Law, C. Knox, Y. Djoumbou et al., “Drugbank 4.0: shedding new light on drug metabolism,” Nucleic Acids Research, vol. 42, no. 1, pp. D1091–D1097, 2014. View at: Publisher Site | Google Scholar
  24. P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003. View at: Publisher Site | Google Scholar
  25. A. Martin, M. E. Ochagavia, L. C. Rabasa, J. Miranda, J. Fernandez-de-Cossio, and R. Bringas, “Bisogenet: a new tool for gene network building, visualization and analysis,” BMC Bioinformatics, vol. 11, no. 1, p. 91, 2010. View at: Publisher Site | Google Scholar
  26. Y. Tang, M. Li, J. Wang, Y. Pan, and F.-X. Wu, “CytoNCA: a cytoscape plugin for centrality analysis and evaluation of protein interaction networks,” Biosystems, vol. 127, pp. 67–72, 2015. View at: Publisher Site | Google Scholar
  27. R. Wang and J. Lin, “Analysis of the mechanism of zhichuanling oral liquid in treating bronchial asthma based on network pharmacology,” Evidence-Based Complementary and Alternative Medicine, vol. 2020, Article ID 1875980, 10 pages, 2020. View at: Publisher Site | Google Scholar
  28. W. D. Lu, L. Li, Y. J. Shen et al., “Network pharmacology study of xiaoxuming decoction based on vasodilatory and vasoconstrictory related GPCR targets,” Zhongguo Zhong Yao Za Zhi, vol. 43, no. 23, pp. 4698–4708, 2018, in Chinese. View at: Google Scholar
  29. Y. Zhang, Z. Li, M. Yang et al., “Identification of GRB2 and GAB1 coexpression as an unfavorable prognostic factor for hepatocellular carcinoma by a combination of expression profile and network analysis,” PLoS One, vol. 8, no. 12, Article ID e85170, 2013. View at: Publisher Site | Google Scholar
  30. Z. Huang, X. Shi, W. Zhu, Z. Zhao, and Z. Wang, “A comparative study on HGWD and its compatibility with anti inflammation and analgesia,” New Chinese Medicine and Clinical Pharmacology, vol. 2, pp. 93–96, 2005, in Chinese. View at: Google Scholar
  31. Y. Xiang, L. Xiao, J. Zhang et al., “Research progress on the effect and mechanism of Astragalus polysaccharides in the treatment of nervous system diseases,” Chinese Journal of Hospital Pharmacy, vol. 8, pp. 687–691, 2016, in Chinese. View at: Google Scholar
  32. L. Zhang, D. Ma, J. Cheng, and L. Ni, “Comparison of cell anti-inflammatory, immune and bone cell repair activities of volatile oil/water extract of traditional Chinese medicine,” New Chinese Medicine and Clinical Pharmacology, vol. 1, pp. 34–39, 2015, in Chinese. View at: Google Scholar
  33. L. Li, D. Zhang, and F. Kunmiao, “Research progress of pharmacological action and mechanism of TGP in the treatment of rheumatoid arthritis,” Shandong Medicine, vol. 59, no. 26, pp. 111–113, 2019, in Chinese. View at: Google Scholar
  34. L. Zhou, J. Zhao, H. Xu, M. Shi, D. Chen, and Y. Luo, “Comparative study on antioxidant effect of turmeric extract of different commodity specifications and grades,” Pharmacology and Clinical of Traditional Chinese Medicine, vol. 1, pp. 110–112, 2016, in Chinese. View at: Google Scholar
  35. X. Ji, Q. Peng, H. Li et al., “Chemical characterization and anti-inflammatory activity of polysaccharides from Zizyphus jujube cv. muzao,” International Journal of Food Engineering, vol. 13, no. 7, 2017. View at: Publisher Site | Google Scholar
  36. J. M. Davis, E. A. Murphy, and M. D. Carmichael, “Effects of the dietary flavonoid quercetin upon performance and health,” Current Sports Medicine Reports, vol. 8, no. 4, pp. 206–213, 2009. View at: Publisher Site | Google Scholar
  37. Y. Li, J. Yao, C. Han et al., “Quercetin, inflammation and immunity,” Nutrients, vol. 8, no. 3, p. 167, 2016. View at: Publisher Site | Google Scholar
  38. M.-S. Tsai, Y.-H. Wang, Y.-Y. Lai et al., “Kaempferol protects against propacetamol-induced acute liver injury through CYP2E1 inactivation, UGT1A1 activation, and attenuation of oxidative stress, inflammation and apoptosis in mice,” Toxicology Letters, vol. 290, pp. 97–109, 2018. View at: Publisher Site | Google Scholar
  39. J. M. Calderón-Montaño, E. Burgos-Morón, C. Pérez-Guerrero, and M. López-Lázaro, “A review on the dietary flavonoid kaempferol,” Mini-Reviews in Medicinal Chemistry, vol. 11, pp. 298–344, 2011. View at: Publisher Site | Google Scholar
  40. M.-K. Lin, Y.-L. Yu, K.-C. Chen et al., “Kaempferol from semen cuscutae attenuates the immune function of dendritic cells,” Immunobiology, vol. 216, no. 10, pp. 1103–1109, 2011. View at: Publisher Site | Google Scholar
  41. M. Valerio and A. B. Awad, “β-Sitosterol down-regulates some pro-inflammatory signal transduction pathways by increasing the activity of tyrosine phosphatase SHP-1 in J774A. 1 murine macrophages,” International Immunopharmacology, vol. 11, no. 8, pp. 1012–1017, 2011. View at: Publisher Site | Google Scholar
  42. M. S. Bin Sayeed and S. S. Ameen, “Beta-sitosterol: a promising but orphan nutraceutical to fight against cancer,” Nutrition and Cancer, vol. 67, no. 8, pp. 1214–1220, 2015. View at: Publisher Site | Google Scholar
  43. L. Wang, Y. Ding, X. Guo, and Q. Zhao, “Role and mechanism of vascular cell adhesion molecule-1 in the development of rheumatoid arthritis,” Experimental and Therapeutic Medicine, vol. 10, no. 3, p. 1229, 2015. View at: Publisher Site | Google Scholar
  44. D. H. Kong, Y. K. Kim, M. R. Kim, J. H. Jang, and S. Lee, “Emerging roles of vascular cell adhesion molecule-1 (VCAM-1) in immunological disorders and cancer,” International Journal of Molecular Sciences, vol. 19, no. 4, p. 1057, 2018. View at: Publisher Site | Google Scholar
  45. B. Bartok and G. S. Firestein, “Fibroblast-like synoviocytes: key effector cells in rheumatoid arthritis,” Immunological Reviews, vol. 233, no. 1, p. 233, 2010. View at: Publisher Site | Google Scholar
  46. G.-Q. Li, Y.-X. Fang, Y. Liu et al., “MALAT1-driven inhibition of wnt signal impedes proliferation and inflammation in fibroblast-like synoviocytes through CTNNB1 promoter methylation in rheumatoid arthritis,” Human Gene Therapy, vol. 30, no. 8, p. 1008, 2019. View at: Publisher Site | Google Scholar
  47. K. Fujii, Y. Tanaka, S. Hubscher et al., “Cross-linking of CD44 on rheumatoid synovial cells up-regulates VCAM-1,” The Journal of Immunology, vol. 167, no. 3, pp. 2391–2398, 1999. View at: Publisher Site | Google Scholar
  48. G. Liu and L. Deng, “Correlation between IL-17 content in synovial fluid and inflammatory injury in patients with rheumatoid arthritis,” Chinese Journal of Experimental Diagnostics, vol. 20, no. 4, pp. 669-670, 2016, in Chinese. View at: Google Scholar
  49. J. A. A. Fischer, A. J. Hueber, S. Wilson et al., “Combined inhibition of tumor necrosis factor α and interleukin-17 as a therapeutic opportunity in rheumatoid arthritis: development and characterization of a novel bispecific antibody,” Arthritis & Rheumatology, vol. 67, no. 1, pp. 51–62, 2015. View at: Publisher Site | Google Scholar
  50. X. Wang, X. Lu, Y. Lu et al., “Correlation between IL-17 level and bone and joint damage in early rheumatoid arthritis patients,” Shandong Medical Journal, vol. 51, no. 1, pp. 102-103, 2011, in Chinese. View at: Google Scholar
  51. J. Yao, G. Han, and M. Hu, “The effect of IL-17 on the expression of MMP-9 in osteoclasts and its significance,” Journal of Jilin University, vol. 42, no. 3, pp. 462–466, 2016, in Chinese. View at: Google Scholar
  52. B. Schminke, S. Trautmann, B. Mai et al., “Interleukin 17 inhabits progenitor cells in rheumatoid arthritis cartilage,” European Journal of Immunology, vol. 46, no. 2, pp. 400–405, 2016. View at: Publisher Site | Google Scholar
  53. M. I. Koenders, R. J. Marijnissen, I. Devesa et al., “Tumor necrosis factor-interleukin-17 interplay induces S100A8, interleukin-1β, and matrix metalloproteinases, and drives irreversible cartilage destruction in murine arthritis: rationale for combination treatment during arthritis,” Arthritis & Rheumatism, vol. 63, no. 8, pp. 2329–2339, 2011. View at: Publisher Site | Google Scholar
  54. H. Zhang and W. Xiao, “TNFR1 and TNFR2 differentially mediate TNF-α-induced inflammatory responses in rheumatoid arthritis fibroblast-like synoviocytes,” Cell Biology International, vol. 41, no. 4, p. 415, 2017. View at: Publisher Site | Google Scholar
  55. R. Eferl and E. F. Wagner, “AP-1: a double-edged sword in tumorigenesis,” Nature Reviews Cancer, vol. 3, no. 11, p. 859, 2003. View at: Publisher Site | Google Scholar
  56. E. F. Wagner and Á. R. Nebreda, “Signal integration by JNK and p38 MAPK pathways in cancer development,” Nature Reviews Cancer, vol. 9, no. 8, p. 537, 2009. View at: Publisher Site | Google Scholar
  57. H. Bazzazi, M. Aghaei, A. Memarian, H. Asgarian-Omran, N. Behnampour, and Y. Yazdani, “Th1–Th17 ratio as a new insight in rheumatoid arthritis disease,” Iranian Journal of Allergy, Asthma and Immunology, vol. 17, no. 1, pp. 68–77, 2018. View at: Google Scholar
  58. A. Tanel, S. G. Fonseca, B. Yassine-Diab et al., “Cellular and molecular mechanisms of memory T-cell survival,” Expert Review of Vaccines, vol. 8, no. 3, p. 299, 2009. View at: Publisher Site | Google Scholar
  59. R. Zhang, L. Wang, and J. Pan, “Expression of E2F2 in rheumatoid arthritis synovial fibroblast like cells,” Journal of Medical Research, vol. 3, pp. 73–77, 2017. View at: Google Scholar
  60. L. Xie, D. Ou, W. Wang et al., “Study on the change characteristics and clinical significance of serum sICAM-1 in patients with rheumatoid arthritis related interstitial lung disease,” Modern Medicine and Health, vol. 20, pp. 3117–3119, 2019. View at: Google Scholar

Copyright © 2020 Wei Liu 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.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views489
Downloads322
Citations

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.