# Integrating Network Pharmacology with Molecular Docking to Unravel the Active Compounds and Potential Mechanism of Simiao Pill Treating Rheumatoid Arthritis

Authors: Mengshi Tang; Xi Xie; Pengji Yi; Jin Kang; Jiafen Liao; Wenqun Li; Fen Li
Journal: Evidence-Based Complementary and Alternative Medicine
 (2020)

---

## Abstract

Objective. To explore the main components and unravel the potential mechanism of simiao pill (SM) on rheumatoid arthritis (RA) based on network pharmacological analysis and molecular docking. Methods. Related compounds were obtained from TCMSP and BATMAN-TCM database. Oral bioavailability and drug-likeness were then screened by using absorption, distribution, metabolism, and excretion (ADME) criteria. Additionally, target genes related to RA were acquired from GeneCards and OMIM database. Correlations about SM-RA, compounds-targets, and pathways-targets-compounds were visualized through Cytoscape 3.7.1. The protein-protein interaction (PPI) network was constructed by STRING. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed via R packages. Molecular docking analysis was constructed by the Molecular Operating Environment (MOE). Results. A total of 72 potential compounds and 77 associated targets of SM were identified. The compounds-targets network analysis indicated that the 6 compounds, including quercetin, kaempferol, baicalein, wogonin, beta-sitosterol, and eugenol, were linked to ≥10 target genes, and the 10 target genes (PTGS1, ESR1, AR, PGR, CHRM3, PPARG, CHRM2, BCL2, CASP3, and RELA) were core target genes in the network. Enrichment analysis indicated that PI3K-Akt, TNF, and IL-17 signaling pathway may be a critical signaling pathway in the network pharmacology. Molecular docking showed that quercetin, kaempferol, baicalein, and wogonin have good binding activity with IL6, VEGFA, EGFR, and NFKBIA targets. Conclusion. The integrative investigation based on bioinformatics/network topology strategy may elaborate on the multicomponent synergy mechanisms of SM against RA and provide the way out to develop new combination medicines for RA.


---

## Body

## 1. Introduction

Rheumatoid arthritis (RA) is a chronic polyarticular symmetric disease. It is characterized by chronic inflammation of the synovial membrane, which can destroy articular cartilage and juxta-articular bone [1]. RA affects 0.3%–1% of the population worldwide [2]. If insufficiently treated, it usually leads to persistent joint inflammation, progressive joint destruction, continuing functional decline, extra-articular manifestations, disability, and increased mortality [3, 4]. Although current available therapeutic approaches against RA, including nonsteroidal anti-inflammatory drugs (NSAIDs), disease-modifying antirheumatic drugs (DMARDs), and corticosteroid, allow for excellent disease control, novel therapies are needed because RA remains incurable [5]. Furthermore, the long-term use of these drugs may cause multiple side effects and lead to limited therapeutic responses. Therefore, novel treatments are in urgent demand.Traditional Chinese medicine (TCM) has been extensively applied for the treatment of RA for centuries in Asia and has been gradually accepted for worldwide clinical applications [6, 7]. Numerous studies have indicated that TCM can be served as complementary and alternative RA drugs for therapeutic effects and with fewer side effects [8, 9].Simiao pill (SM), a traditional TCM formula, comprises four herbs, including Phellodendri Chinensis Cortex (Huang Bo),Atractylodes lancea (Thunb.) Dc. (Cang Zhu), achyranthis bidentatae radix (Niu Xi), and Coicis Semen (Yi Yi Ren). Previous studies have indicated the anti-inflammation pharmacological effect of SM [10] and that SM reduced proinflammatory cytokine production by suppressing nuclear factor kappaB (NF-κB)/pyrin domain containing 3 (NLRP3) inflammasome activation [11]. Recently, SM was demonstrated to exhibit anti-inflammatory and bone-protective effects by regulating autotaxin (ATX)-lysophosphatidic acid (LPA) and mitogen-activated protein kinase (MAPK) signaling pathways in collagen-induced arthritis (CIA) rats [12]. In addition, SM was recommended for the treatment of active RA (53.6%) in the expert consensus regarding the treatment of RA with various Chinese patent medicines (CPMs) [13]. However, because TCM formulas are characterized by multicomponents, multitargets, and multipathways [14], the therapeutic effect of SM against RA has not been fully elucidated. Therefore, it is necessary for further systematic investigation.Nowadays, network pharmacology integrates network biology and polypharmacology based on existing databases, providing a novel approach for exploring the mechanisms and synergistic effect of TCM formulas as disease treatments [14–16]. Combining network science with ancient TCM formulas to investigate multiple molecular mechanisms has achieved successful attempts in the previous researches [17–20].Therefore, in this study, a network pharmacology-based study was conducted to predict bioactive compounds and elucidate the comprehensive pharmacological mechanisms about the antirheumatic effect of SM. In addition, molecular docking analysis was performed to validate in silico to predict molecular interactions between compounds and targets.


## 2. Materials and Methods

Network pharmacology-based prediction of SM treating RA was constructed by the following (Figure1): (1) data collection and preparation, including retrieving the ingredients list of SM formula, screening for candidate compounds, identifying SM and RA targets, and intersecting the identified targets of compounds and disease; (2) topological analysis of network and protein-protein interaction (PPI) network construction; (3) enrichment analysis; and (4) molecular docking analysis.Figure 1
Workflow of network pharmacology analysis.


### 2.1. Data Collection and Preparation

#### 2.1.1. Composite Compounds of SM

The related composite compounds of SM were obtained from the Traditional Chinese Medicine Systems Pharmacology Database (TCMSP,http://lsp.nwu.edu.cn/tcmsp.php) and a Bioinformatics Analysis Tool for Molecular mechANism of Traditional Chinese Medicine (BATMAN-TCM, http://bionet.ncpsb.org/batman-tcm/).


#### 2.1.2. Pharmacokinetic ADME Evaluation

The in silico integrative ADME (absorption, distribution, metabolism, and excretion) model administrated by TCMSP is employed for pharmaceutical research. As an oral drug, two related-ADME models, oral bioavailability (OB), and drug-likeness (DL) are applied to identify the potential bioactive compounds in this study. Only the compounds with OB ≥ 30 and DL ≥ 0.18 that satisfied the criteria suggested by the TCMSP database (removed the duplicated) are retained as the candidate compounds for further study [21]. In addition, among the compounds with OB ＜30 or DL ＜0.18, which are searched with “compound (name)” and “rheumatoid arthritis” [all fields] in PubMed databases to find relevant researches, the compounds in purified form focused on anti-RA mechanisms are also considered to be bioactive compounds (removed the duplicated) and included for further study.


#### 2.1.3. Predictions of Target Genes Related to the Identified Compounds

All the potential compounds were input into TCMSP to capture the relationships between drugs and targets. Since the obtained targets include various biological species, all target names were also put into UniProt databases (http://www.uniprot.org/) to search for target gene names selected by human species.


#### 2.1.4. Potential Disease Target Genes

Information of known RA-related therapeutic target genes was collected by keywords “rheumatoid arthritis” as queries from The Human Gene Databases (GeneCards,https://www.genecards.org/, ver.4.9.0) and Online Mendelian Inheritance in Man (OMIM, http://www.omim.org/, updated June 6, 2019), and only “Homo sapiens” target genes linked to RA are selected.


#### 2.1.5. Venn Analysis

All target genes of identified compounds and RA are put into Bioinformatics and Evolutionary Genomics system (bioinformatics.psb.ugent.be/webtools/Venn/), respectively, to produce a Venn diagram, which indicates the intersection of identified targets of drug and disease.


### 2.2. Topological Analysis of Network and PPI Network Construction

#### 2.2.1. Topological Network Analysis

SM-RA mechanism network, compounds-targets network, and pathways-targets-compounds network were visualized through Cytoscape (https://cytoscape.org/, ver. 3.7.1) to systemically explore the molecular mechanisms of SM treating RA.


#### 2.2.2. PPI Network Construction

The above 77 target genes acquired from the Venn diagram intersection were imported into STRING (https://string-db.org/, version 11.0) to construct a PPI network for understanding protein interaction systematically. The PPI network is constructed by setting the organism as “human sapiens”, setting the minimum required interaction score to “medium confidence (0.40)”, and excluding the disconnected protein nodes. In addition, statistics of protein interactions are figured out according to the PPI network, and a related bar plot diagram is constructed with R 3.6.0 subsequently.


### 2.3. Enrichment Analysis

R 3.6.0 and relatedR packages (colorspace, stringi, DOSE, clusterProfiler, and pathview) are applied to carry out Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of intersection target genes of SM-RA. P values ＜0.05 and q values ＜0.05 are considered statistically significant based on Fisher’s test.


### 2.4. Molecular Docking Analysis

The 3D structures of candidate targets were obtained from the PDB database (http://www.rcsb.org/) in PDB format by setting the organism to “Homo sapiens only”. The 3D conformers of candidate compounds are acquired from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) with SDF format. Subsequently, they were imported to the Molecular Operating Environment (MOE) to get the docking score. The greater the absolute value of the docking score, the better.


## 2.1. Data Collection and Preparation

### 2.1.1. Composite Compounds of SM

The related composite compounds of SM were obtained from the Traditional Chinese Medicine Systems Pharmacology Database (TCMSP,http://lsp.nwu.edu.cn/tcmsp.php) and a Bioinformatics Analysis Tool for Molecular mechANism of Traditional Chinese Medicine (BATMAN-TCM, http://bionet.ncpsb.org/batman-tcm/).


### 2.1.2. Pharmacokinetic ADME Evaluation

The in silico integrative ADME (absorption, distribution, metabolism, and excretion) model administrated by TCMSP is employed for pharmaceutical research. As an oral drug, two related-ADME models, oral bioavailability (OB), and drug-likeness (DL) are applied to identify the potential bioactive compounds in this study. Only the compounds with OB ≥ 30 and DL ≥ 0.18 that satisfied the criteria suggested by the TCMSP database (removed the duplicated) are retained as the candidate compounds for further study [21]. In addition, among the compounds with OB ＜30 or DL ＜0.18, which are searched with “compound (name)” and “rheumatoid arthritis” [all fields] in PubMed databases to find relevant researches, the compounds in purified form focused on anti-RA mechanisms are also considered to be bioactive compounds (removed the duplicated) and included for further study.


### 2.1.3. Predictions of Target Genes Related to the Identified Compounds

All the potential compounds were input into TCMSP to capture the relationships between drugs and targets. Since the obtained targets include various biological species, all target names were also put into UniProt databases (http://www.uniprot.org/) to search for target gene names selected by human species.


### 2.1.4. Potential Disease Target Genes

Information of known RA-related therapeutic target genes was collected by keywords “rheumatoid arthritis” as queries from The Human Gene Databases (GeneCards,https://www.genecards.org/, ver.4.9.0) and Online Mendelian Inheritance in Man (OMIM, http://www.omim.org/, updated June 6, 2019), and only “Homo sapiens” target genes linked to RA are selected.


### 2.1.5. Venn Analysis

All target genes of identified compounds and RA are put into Bioinformatics and Evolutionary Genomics system (bioinformatics.psb.ugent.be/webtools/Venn/), respectively, to produce a Venn diagram, which indicates the intersection of identified targets of drug and disease.


## 2.1.1. Composite Compounds of SM

The related composite compounds of SM were obtained from the Traditional Chinese Medicine Systems Pharmacology Database (TCMSP,http://lsp.nwu.edu.cn/tcmsp.php) and a Bioinformatics Analysis Tool for Molecular mechANism of Traditional Chinese Medicine (BATMAN-TCM, http://bionet.ncpsb.org/batman-tcm/).


## 2.1.2. Pharmacokinetic ADME Evaluation

The in silico integrative ADME (absorption, distribution, metabolism, and excretion) model administrated by TCMSP is employed for pharmaceutical research. As an oral drug, two related-ADME models, oral bioavailability (OB), and drug-likeness (DL) are applied to identify the potential bioactive compounds in this study. Only the compounds with OB ≥ 30 and DL ≥ 0.18 that satisfied the criteria suggested by the TCMSP database (removed the duplicated) are retained as the candidate compounds for further study [21]. In addition, among the compounds with OB ＜30 or DL ＜0.18, which are searched with “compound (name)” and “rheumatoid arthritis” [all fields] in PubMed databases to find relevant researches, the compounds in purified form focused on anti-RA mechanisms are also considered to be bioactive compounds (removed the duplicated) and included for further study.


## 2.1.3. Predictions of Target Genes Related to the Identified Compounds

All the potential compounds were input into TCMSP to capture the relationships between drugs and targets. Since the obtained targets include various biological species, all target names were also put into UniProt databases (http://www.uniprot.org/) to search for target gene names selected by human species.


## 2.1.4. Potential Disease Target Genes

Information of known RA-related therapeutic target genes was collected by keywords “rheumatoid arthritis” as queries from The Human Gene Databases (GeneCards,https://www.genecards.org/, ver.4.9.0) and Online Mendelian Inheritance in Man (OMIM, http://www.omim.org/, updated June 6, 2019), and only “Homo sapiens” target genes linked to RA are selected.


## 2.1.5. Venn Analysis

All target genes of identified compounds and RA are put into Bioinformatics and Evolutionary Genomics system (bioinformatics.psb.ugent.be/webtools/Venn/), respectively, to produce a Venn diagram, which indicates the intersection of identified targets of drug and disease.


## 2.2. Topological Analysis of Network and PPI Network Construction

### 2.2.1. Topological Network Analysis

SM-RA mechanism network, compounds-targets network, and pathways-targets-compounds network were visualized through Cytoscape (https://cytoscape.org/, ver. 3.7.1) to systemically explore the molecular mechanisms of SM treating RA.


### 2.2.2. PPI Network Construction

The above 77 target genes acquired from the Venn diagram intersection were imported into STRING (https://string-db.org/, version 11.0) to construct a PPI network for understanding protein interaction systematically. The PPI network is constructed by setting the organism as “human sapiens”, setting the minimum required interaction score to “medium confidence (0.40)”, and excluding the disconnected protein nodes. In addition, statistics of protein interactions are figured out according to the PPI network, and a related bar plot diagram is constructed with R 3.6.0 subsequently.


## 2.2.1. Topological Network Analysis

SM-RA mechanism network, compounds-targets network, and pathways-targets-compounds network were visualized through Cytoscape (https://cytoscape.org/, ver. 3.7.1) to systemically explore the molecular mechanisms of SM treating RA.


## 2.2.2. PPI Network Construction

The above 77 target genes acquired from the Venn diagram intersection were imported into STRING (https://string-db.org/, version 11.0) to construct a PPI network for understanding protein interaction systematically. The PPI network is constructed by setting the organism as “human sapiens”, setting the minimum required interaction score to “medium confidence (0.40)”, and excluding the disconnected protein nodes. In addition, statistics of protein interactions are figured out according to the PPI network, and a related bar plot diagram is constructed with R 3.6.0 subsequently.


## 2.3. Enrichment Analysis

R 3.6.0 and relatedR packages (colorspace, stringi, DOSE, clusterProfiler, and pathview) are applied to carry out Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of intersection target genes of SM-RA. P values ＜0.05 and q values ＜0.05 are considered statistically significant based on Fisher’s test.


## 2.4. Molecular Docking Analysis

The 3D structures of candidate targets were obtained from the PDB database (http://www.rcsb.org/) in PDB format by setting the organism to “Homo sapiens only”. The 3D conformers of candidate compounds are acquired from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) with SDF format. Subsequently, they were imported to the Molecular Operating Environment (MOE) to get the docking score. The greater the absolute value of the docking score, the better.


## 3. Results

### 3.1. Data Collection and Preparation

#### 3.1.1. Identification of Compounds in SM

A total of 479 compounds were identified in SM, including 140 in Phellodendri Chinensis Cortex (Huang Bo), 49 inAtractylodes lancea (Thunb.) Dc. (Cang Zhu), 176 in achyranthis bidentatae radix (Niu Xi), and 38 in Coicis Semen (Yi Yi Ren), and in TCMSP, also including 37 in Phellodendri Chinensis Cortex (Huang Bo), 26 in Atractylodes lancea (Thunb.) Dc. (Cang Zhu), 10 in achyranthis bidentatae radix (Niu Xi), and 3 in Coicis Semen (Yi Yi Ren) in BATMAN-TCM.


#### 3.1.2. Selection of Compounds Using ADME Screening and Related Targets

All the identified compounds were selected through ADME screening, with 90 of 479 compounds satisfying the suggested criteria OB ≥ 30 and DL ≥ 0.18 [18–20]. Of the 90 compounds, 25 were duplicated and removed, and the remaining 65 compounds were included for further study. Moreover, of the excluded compounds that do not meet the suggested criteria, 7 compounds, including ferulic acid, beta-elemene, eugenol, and paeonol in Phellodendri Chinensis Cortex (Huang Bo) and geniposide, rutin, and astragalin in achyranthis bidentatae radix (Niu Xi), are considered bioactive compounds and included for further analysis, and the effects of ferulic acid [22], beta-elemene [23], eugenol [24–26], paeonol [27–29], geniposide [30–32], rutin [33, 34], and astragalin [35] on RA have been investigated. The final 72 compounds are selected from the four herbal medicines (Table 1). A total of 386 target genes related to the final identified compounds are obtained from the UniProt databases.Table 1
72 active compounds of SM.
Mol idMolecule nameOB (%)DLHerbMOL002636Kihadalactone A34.210.82Phellodendri Chinensis CortexMOL013352Obacunone43.290.77Phellodendri Chinensis CortexMOL002641Phellavin_qt35.860.44Phellodendri Chinensis CortexMOL002644Phellopterin40.190.28Phellodendri Chinensis CortexMOL002651Dehydrotanshinone II A43.760.40Phellodendri Chinensis CortexMOL002652delta7-dehydrosophoramine54.450.25Phellodendri Chinensis CortexMOL002656Dihydroniloticin36.430.81Phellodendri Chinensis CortexMOL002659Kihadanin A31.600.70Phellodendri Chinensis CortexMOL002660Niloticin41.410.82Phellodendri Chinensis CortexMOL002662Rutaecarpine40.300.60Phellodendri Chinensis CortexMOL002663Skimmianin40.140.20Phellodendri Chinensis CortexMOL002666Chelerythrine34.180.78Phellodendri Chinensis CortexMOL002668Worenine45.830.87Phellodendri Chinensis CortexMOL002670Cavidine35.640.81Phellodendri Chinensis CortexMOL002671Candletoxin A31.810.69Phellodendri Chinensis CortexMOL002672Hericenone H39.000.63Phellodendri Chinensis CortexMOL002673Hispidone36.180.83Phellodendri Chinensis CortexMOL000358Beta-sitosterol36.910.75Phellodendri Chinensis CortexMOL000622Magnograndiolide63.710.19Phellodendri Chinensis CortexMOL000762Palmidin A35.360.65Phellodendri Chinensis CortexMOL000787Fumarine59.260.83Phellodendri Chinensis CortexMOL000790Isocorypalmine35.770.59Phellodendri Chinensis CortexMOL001131phellamurin_qt56.600.39Phellodendri Chinensis CortexMOL001455(S)-canadine53.830.77Phellodendri Chinensis CortexMOL001771Poriferast-5-en-3beta-ol36.910.75Phellodendri Chinensis CortexMOL002894Berberrubine35.740.73Phellodendri Chinensis CortexMOL005438Campesterol37.580.71Phellodendri Chinensis CortexMOL006392Dihydroniloticin36.430.82Phellodendri Chinensis CortexMOL006401Melianone40.530.78Phellodendri Chinensis CortexMOL006413Phellochin35.410.82Phellodendri Chinensis CortexMOL006422Thalifendine44.410.73Phellodendri Chinensis CortexMOL002665Ferulic acid40.430.06Phellodendri Chinensis CortexMOL000908Beta-elemene25.630.06Phellodendri Chinensis CortexMOL000254Eugenol56.240.04Phellodendri Chinensis CortexMOL000874Paeonol28.790.04Phellodendri Chinensis CortexMOL0001792-Hydroxyisoxypropyl-3-hydroxy-7-isopentene-2,3-dihydrobenzofuran-5-carboxylic45.200.20Atractylodes lancea (Thunb.) Dc.MOL000184NSC6355139.250.76Atractylodes lancea (Thunb.) Dc.MOL000186Stigmasterol 3-O-beta-D-glucopyranoside_qt43.830.76Atractylodes lancea (Thunb.) Dc.MOL0001883β-acetoxyatractylone40.570.22Atractylodes lancea (Thunb.) Dc.MOL000088Beta-sitosterol 3-O-glucoside_qt36.910.75Atractylodes lancea (Thunb.) Dc.MOL000092daucosterin_qt36.910.76Atractylodes lancea (Thunb.) Dc.MOL000094daucosterol_qt36.910.76Atractylodes lancea (Thunb.) Dc.MOL001006Poriferasta-7,22E-dien-3beta-ol42.980.76Achyranthis bidentatae radixMOL01246128-Norolean-17-en-3-ol35.930.78Achyranthis bidentatae radixMOL012505bidentatoside, ii_qt31.760.59Achyranthis bidentatae radixMOL012537Spinoside A41.750.40Achyranthis bidentatae radixMOL012542β-ecdysterone44.230.82Achyranthis bidentatae radixMOL002714Baicalein33.520.21Achyranthis bidentatae radixMOL002776Baicalin40.120.75Achyranthis bidentatae radixMOL002897Epiberberine43.090.78Achyranthis bidentatae radixMOL003847Inophyllum E38.810.85Achyranthis bidentatae radixMOL000422Kaempferol41.880.24Achyranthis bidentatae radixMOL004355Spinasterol42.980.76Achyranthis bidentatae radixMOL012516Geniposide8.400.44Achyranthis bidentatae radixMOL000415Rutin3.200.68Achyranthis bidentatae radixMOL000561Astragalin14.030.74Achyranthis bidentatae radixMOL001323Sitosterol alpha143.280.78Coicis SemenMOL001494Mandenol42.000.19Coicis SemenMOL002372(6Z, 10 E, 14 E, 18 E)-2,6,10,15,19,23-hexamethyltetracosa-2,6,10,14,18,22-hexaene33.550.42Coicis SemenMOL002882[(2R)-2,3-dihydroxypropyl] (Z)-octadec-9-enoate34.130.30Coicis SemenMOL000359Sitosterol36.910.75Coicis SemenMOL008118Coixenolide32.400.43Coicis SemenMOL0081212-Monoolein34.230.29Coicis SemenMOL000953CLR37.870.68Coicis SemenMOL001454Berberine36.860.78Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL001458Coptisine30.670.86Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL002643Delta 7-stigmastenol37.420.75Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL000785Palmatine64.600.65Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL000098Quercetin46.430.28Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL000173Wogonin30.680.23Atractylodes lancea (Thunb.) Dc. Achyranthis bidentatae radixMOL000085Beta-daucosterol_qt36.910.75Atractylodes lancea (Thunb.) Dc. Achyranthis bidentatae radixMOL000449Stigmasterol43.830.76Phellodendri Chinensis Cortex, achyranthis bidentatae radix, Coicis Semen


#### 3.1.3. Identified Disease Target Genes

The target genes related to RA were searched in GeneCards and OMIM databases, which include 3768 genes in GeneCards and 1 gene in OMIM, with no overlapping target gene.


#### 3.1.4. Intersection of Identified Targets of Compounds and Disease

In the Venn diagram intersection of identified targets about identified compounds and of RA (Figure2), a total of 77 target genes are acquired.Figure 2
A Venn diagram showing intersection identified targets of identified compounds and RA.


### 3.2. Topological and PPI Network

#### 3.2.1. Topological Network Analysis

The SM-RA mechanism network (Figure3) consists of 77 target genes nodes (shared gene of SM and RA), 41 compound nodes, and 349 edges. Among the 14 compounds (Dehydrotanshinone II A, Stigmasterol, beta-sitosterol, Isocorypalmine, beta-elemene, quercetin, eugenol, paeonol, (S)-Canadine, wogonin, baicalein, Inophyllum E, rutin, and kaempferol) that connected to more than four genes, 55 target genes are associated with quercetin, 22 target genes are associated with kaempferol, 15 target genes are associated with baicalein, 14 target genes are associated with wogonin, and 10 target genes are associated with beta-sitosterol and eugenol, respectively (Table 2). In addition, 10 genes, including PTGS1, ESR1, AR, PGR, CHRM3, PPARG, CHRM2, BCL2, CASP3, and RELA, are related to more than five compounds, as shown in the compounds-targets network (Figure 4). These compounds and genes may be the key nodes in the network.Figure 3
The SM-RA mechanism network. The green octagon represents rheumatoid arthritis (RA), the purple hexagon represents the herbal medicine simiao pill (SM), while pink V’s represent compounds, and blue-purple ellipses represent genes.Table 2
Target genes interacting with compounds in the SM-RA network.
CompoundsTarget genesDelta 7-stigmastenolPGRPoriferast-5-en-3beta-olPGRCampesterolPGRNSC63551PGRBeta-daucosterol_qtPGRPoriferasta-7,22E-dien-3beta-olPGR28-Norolean-17-en-3-olPGRSpinasterolPGRSitosterol alpha1PGRSitosterolPGRCLRPGRChelerythrinePTGS1AstragalinPTGS1MandenolPTGS1Ferulic acidPTGS1, CHRM2phellamurin_qtESR1, NR3C1EpiberberineESR1, ARBerberinePTGS1, ESR1, ARCoptisinePTGS1, ESR1, ARWoreninePTGS1, ESR1, ARBerberrubinePTGS1, ESR1, ARThalifendinePTGS1, ESR1, ARFumarinePTGS1, CHRM3, HTR3ARutaecarpinePTGS1, AR, HTR3A, CYP3A4CavidinePTGS1, CHRM3, HTR3A, RXRBPalmatinePTGS1, ESR1, AR, ESR23β-acetoxyatractyloneCHRM3, AR, ACHE, CHRM2Dehydrotanshinone II ACHRM3, ESR1, AR, PPARG, ACHEInophyllum EPTGS1, ESR1, AR, ESR2, GSK3BStigmasterolPGR, PTGS1, AKR1B1, PLAU, CHRM3, CHRM2IsocorypalminePTGS1, CHRM3, HTR3A, CHRM2, DRD3, RXRB(S)-canadinePTGS1, CHRM3, HTR3A, CHRM2, DRD3, RXRBPaeonolPTGS1, CHRM2, RELA, BCL2, NFKBIA, ICAM1, TYRBeta-elemeneCHRM2, PTGS1, CHRM3, BCL2, RB1, TP63, CCNB1, RHOARutinRELA, IL6, CASP3, ALOX5, GSTP1, INS, FCER2, ITGB2Beta-sitosterolPGR, PTGS1, CHRM3, CHRM2, BCL2, CASP9, CASP3, CASP8, PRKCA, PON1EugenolPTGS1, CHRM3, CHRM2, PLAU, RELA, CYP1A1, ALOX5, AHR, ABCC2, MUC1WogoninPTGS1, ESR1, AR, PPARG, GSK3B, RELA, CCND1, BCL2, CASP9, IL6, CASP3, TP63, PTGER3, MCL1BaicaleinPTGS1, AR, RELA, VEGFA, BCL2, FOS, CASP3, TP63, HIF1A, FOSL1, CCNB1, AHR, IGF2, CYCS, NOX5KaempferolPTGS1, AR, PPARG, PGR, ACHE, CHRM2, RELA, IKBKB, BCL2, CASP3, MAPK8, PPARG, CYP3A4, CYP1A1, ICAM1, SELE, VCAM1, ALOX5, GSTP1, AHR, NR1I3, GSTM1QuercetinPTGS1, AR, PPARG, AKR1B1, ACHE, RELA, EGFR, VEGFA, CCND1, BCL2, FOS, CASP9, PLAU, RB1, IL6, CASP3, TP63, NFKBIA, CASP8,RAF1, PRKCA, HIF1A,ERBB2,PPARG,CYP3A4,CAV1,MYC,CYP1A1, ICAM1, SELE, VCAM1, PTGER3, BIRC5, NOS3, HSPB1, CCNB1, ALOX5, GSTP1, NFE2L2, NQO1, PARP1, AHR, COL3A1, DCAF5, NR1I3, HSF1, CRP, RUNX2, CTSD, IGFBP3, IGF2, IRF1, ERBB3, PON1, GSTM1Figure 4
The compounds-targets network. The pink V’s represent compounds and blue-purple ellipses represent target genes.


#### 3.2.2. PPI Network

The PPI network is established by setting the confidence level of more than 0.40 and hiding the independent target protein nodes. The PPI network nodes represent proteins and edges represent protein-protein interactions. The network has 75 nodes and 1604 edges (Figure5). In addition, we analyzed the importance prioritization (adjacent nodes count of each protein) of proteins according to the network, and the leading 30 genes with higher connection were visualized by constructing a bar plot diagram (Figure 6), which indicates the 30 genes or proteins that may play a bridge role in connecting other nodes in the PPI network. These 30 genes or proteins include inflammation-associated genes (IL6 [36], NFKBIA [37]), cell proliferation-, differentiation-, and transformation-related genes (FOS [38], EGFR [39], MAPK8 [40], NR3C1 [41], RHOA [42], and PARP1 [43]), cell apoptosis-related genes (CASP3, CASP8 [44], CASP9 [45], MYC [46], CYCS [47], HIF1A [48], MCL1 [49], and GSK3B [50]), cell cycle-related gene (CCND1 [51]), hormone-related genes (INS [52], ESR1 [53], AR [54], and PGR [55]), angiogenesis-related gene (VEGFA [56]), and transcription factor (RELA [57]).Figure 5
The PPI network of SM-RA. Each node represents the relevant gene, and the edges represent protein-protein associations.Figure 6
Hub top 30 genes of the PPI network. They-axis displays significant top 30 genes, and the x-axis shows line counts of these genes.


### 3.3. Enrichment Analysis

#### 3.3.1. GO Enrichment Analysis

GO analysis consisted of biological process (BP), cellular component (CC), and molecular function (MF). As showed in Figure7, the top 20 enrichment terms are visualized by the bar plot diagram. The results demonstrated that numerous targets are involved in various BPs associated with immune response and inflammation, such as the response to a steroid hormone, response to oxidative stress, and regulation of the apoptotic signaling pathway, which confirmed strongly the correlation with the pathogenesis in RA. The CC results showed that most of the targets are localized to the cellular membrane and nuclear chromatin part. The MF results indicated that many targets are associated with nuclear receptor activity and transcription factor activity.Figure 7
GO analysis of targets, the top 20 significant enrichment terms in BP (a), CC (b), and MF (c). They-axis shows significantly enriched biological process, cellular component, and molecular function categories of the target genes, respectively. The redder the color, the lower the P value. The x-axis displays the enrichment scores of these terms, and the length of the bar indicates the number of target genes in each pathway.
(a)(b)(c)


#### 3.3.2. KEGG Enrichment Analysis

The KEGG pathways are applied to examine the function and signaling pathways of the identified target genes, with the top 20 of the potential pathways (P<0.05 and q<0.05) shown by a bar plot diagram (Figure 8) and visualized with the pathways-targets-compounds network (Figure 9). The results showed that numerous targets are associated with certain virus infections (such as Epstein-Barr virus infection) and cancer, which are associated with the onset and prognosis of RA.Figure 8
KEGG enrichment pathways (top 20). They-axis displays the top 20 significantly enriched KEGG pathways of the target genes. The redder the color, the smaller the P value. The x-axis represents the target genes counts, and the length of the bar indicates the number of target genes in each pathway.Figure 9
The pathways-targets-compounds network. The green diamonds represent pathways, the blue-purple ellipses represent genes, and the pink V’s represent compounds.


### 3.4. Molecular Docking Analysis

The selected targets, including IL6, VEGFA, EGFR, and NFKBIA, play a significant role in the SM-RA network. The candidate compounds, including quercetin, kaempferol, baicalein, and wogonin, are the top 4 compounds (ranking by related target genes count) in the SM-RA network. These 4 target genes and 4 compounds are imported into MOE for molecular docking verification. The docking scores are shown in Table3. The action mode of NFKBIA and quercetin, kaempferol, baicalein, and wogonin and the action mode of wogonin and IL6, VEGFA, EGFR, and NFKBIA are shown in Figure 10.Table 3
Molecular docking scores.
IL6VEGFAEGFRNFKBIAQuercetin−4.7051−5.9131−6.0857−6.7291Kaempferol−4.9898−5.4844−5.7466−6.5524Baicalein−4.5010−5.5058−5.8372−6.2222Wogonin−4.6678−5.6466−6.1084−6.6169Figure 10
Molecular docking results. (a) The action mode of NFKBIA and quercetin, kaempferol, baicalein, and wogonin: (A) NFKBIA and quercetin; (B) NFKBIA and kaempferol; (C) NFKBIA and baicalein; (D) NFKBIA and wogonin. (b) The action mode of wogonin and IL6, VEGFA, EGFR, and NFKBIA: (A) wogonin and IL6; (B) wogonin and VEGFA; (C) wogonin and EGFR; (D) wogonin and NFKBIA.
(a)(b)


## 3.1. Data Collection and Preparation

### 3.1.1. Identification of Compounds in SM

A total of 479 compounds were identified in SM, including 140 in Phellodendri Chinensis Cortex (Huang Bo), 49 inAtractylodes lancea (Thunb.) Dc. (Cang Zhu), 176 in achyranthis bidentatae radix (Niu Xi), and 38 in Coicis Semen (Yi Yi Ren), and in TCMSP, also including 37 in Phellodendri Chinensis Cortex (Huang Bo), 26 in Atractylodes lancea (Thunb.) Dc. (Cang Zhu), 10 in achyranthis bidentatae radix (Niu Xi), and 3 in Coicis Semen (Yi Yi Ren) in BATMAN-TCM.


### 3.1.2. Selection of Compounds Using ADME Screening and Related Targets

All the identified compounds were selected through ADME screening, with 90 of 479 compounds satisfying the suggested criteria OB ≥ 30 and DL ≥ 0.18 [18–20]. Of the 90 compounds, 25 were duplicated and removed, and the remaining 65 compounds were included for further study. Moreover, of the excluded compounds that do not meet the suggested criteria, 7 compounds, including ferulic acid, beta-elemene, eugenol, and paeonol in Phellodendri Chinensis Cortex (Huang Bo) and geniposide, rutin, and astragalin in achyranthis bidentatae radix (Niu Xi), are considered bioactive compounds and included for further analysis, and the effects of ferulic acid [22], beta-elemene [23], eugenol [24–26], paeonol [27–29], geniposide [30–32], rutin [33, 34], and astragalin [35] on RA have been investigated. The final 72 compounds are selected from the four herbal medicines (Table 1). A total of 386 target genes related to the final identified compounds are obtained from the UniProt databases.Table 1
72 active compounds of SM.
Mol idMolecule nameOB (%)DLHerbMOL002636Kihadalactone A34.210.82Phellodendri Chinensis CortexMOL013352Obacunone43.290.77Phellodendri Chinensis CortexMOL002641Phellavin_qt35.860.44Phellodendri Chinensis CortexMOL002644Phellopterin40.190.28Phellodendri Chinensis CortexMOL002651Dehydrotanshinone II A43.760.40Phellodendri Chinensis CortexMOL002652delta7-dehydrosophoramine54.450.25Phellodendri Chinensis CortexMOL002656Dihydroniloticin36.430.81Phellodendri Chinensis CortexMOL002659Kihadanin A31.600.70Phellodendri Chinensis CortexMOL002660Niloticin41.410.82Phellodendri Chinensis CortexMOL002662Rutaecarpine40.300.60Phellodendri Chinensis CortexMOL002663Skimmianin40.140.20Phellodendri Chinensis CortexMOL002666Chelerythrine34.180.78Phellodendri Chinensis CortexMOL002668Worenine45.830.87Phellodendri Chinensis CortexMOL002670Cavidine35.640.81Phellodendri Chinensis CortexMOL002671Candletoxin A31.810.69Phellodendri Chinensis CortexMOL002672Hericenone H39.000.63Phellodendri Chinensis CortexMOL002673Hispidone36.180.83Phellodendri Chinensis CortexMOL000358Beta-sitosterol36.910.75Phellodendri Chinensis CortexMOL000622Magnograndiolide63.710.19Phellodendri Chinensis CortexMOL000762Palmidin A35.360.65Phellodendri Chinensis CortexMOL000787Fumarine59.260.83Phellodendri Chinensis CortexMOL000790Isocorypalmine35.770.59Phellodendri Chinensis CortexMOL001131phellamurin_qt56.600.39Phellodendri Chinensis CortexMOL001455(S)-canadine53.830.77Phellodendri Chinensis CortexMOL001771Poriferast-5-en-3beta-ol36.910.75Phellodendri Chinensis CortexMOL002894Berberrubine35.740.73Phellodendri Chinensis CortexMOL005438Campesterol37.580.71Phellodendri Chinensis CortexMOL006392Dihydroniloticin36.430.82Phellodendri Chinensis CortexMOL006401Melianone40.530.78Phellodendri Chinensis CortexMOL006413Phellochin35.410.82Phellodendri Chinensis CortexMOL006422Thalifendine44.410.73Phellodendri Chinensis CortexMOL002665Ferulic acid40.430.06Phellodendri Chinensis CortexMOL000908Beta-elemene25.630.06Phellodendri Chinensis CortexMOL000254Eugenol56.240.04Phellodendri Chinensis CortexMOL000874Paeonol28.790.04Phellodendri Chinensis CortexMOL0001792-Hydroxyisoxypropyl-3-hydroxy-7-isopentene-2,3-dihydrobenzofuran-5-carboxylic45.200.20Atractylodes lancea (Thunb.) Dc.MOL000184NSC6355139.250.76Atractylodes lancea (Thunb.) Dc.MOL000186Stigmasterol 3-O-beta-D-glucopyranoside_qt43.830.76Atractylodes lancea (Thunb.) Dc.MOL0001883β-acetoxyatractylone40.570.22Atractylodes lancea (Thunb.) Dc.MOL000088Beta-sitosterol 3-O-glucoside_qt36.910.75Atractylodes lancea (Thunb.) Dc.MOL000092daucosterin_qt36.910.76Atractylodes lancea (Thunb.) Dc.MOL000094daucosterol_qt36.910.76Atractylodes lancea (Thunb.) Dc.MOL001006Poriferasta-7,22E-dien-3beta-ol42.980.76Achyranthis bidentatae radixMOL01246128-Norolean-17-en-3-ol35.930.78Achyranthis bidentatae radixMOL012505bidentatoside, ii_qt31.760.59Achyranthis bidentatae radixMOL012537Spinoside A41.750.40Achyranthis bidentatae radixMOL012542β-ecdysterone44.230.82Achyranthis bidentatae radixMOL002714Baicalein33.520.21Achyranthis bidentatae radixMOL002776Baicalin40.120.75Achyranthis bidentatae radixMOL002897Epiberberine43.090.78Achyranthis bidentatae radixMOL003847Inophyllum E38.810.85Achyranthis bidentatae radixMOL000422Kaempferol41.880.24Achyranthis bidentatae radixMOL004355Spinasterol42.980.76Achyranthis bidentatae radixMOL012516Geniposide8.400.44Achyranthis bidentatae radixMOL000415Rutin3.200.68Achyranthis bidentatae radixMOL000561Astragalin14.030.74Achyranthis bidentatae radixMOL001323Sitosterol alpha143.280.78Coicis SemenMOL001494Mandenol42.000.19Coicis SemenMOL002372(6Z, 10 E, 14 E, 18 E)-2,6,10,15,19,23-hexamethyltetracosa-2,6,10,14,18,22-hexaene33.550.42Coicis SemenMOL002882[(2R)-2,3-dihydroxypropyl] (Z)-octadec-9-enoate34.130.30Coicis SemenMOL000359Sitosterol36.910.75Coicis SemenMOL008118Coixenolide32.400.43Coicis SemenMOL0081212-Monoolein34.230.29Coicis SemenMOL000953CLR37.870.68Coicis SemenMOL001454Berberine36.860.78Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL001458Coptisine30.670.86Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL002643Delta 7-stigmastenol37.420.75Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL000785Palmatine64.600.65Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL000098Quercetin46.430.28Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL000173Wogonin30.680.23Atractylodes lancea (Thunb.) Dc. Achyranthis bidentatae radixMOL000085Beta-daucosterol_qt36.910.75Atractylodes lancea (Thunb.) Dc. Achyranthis bidentatae radixMOL000449Stigmasterol43.830.76Phellodendri Chinensis Cortex, achyranthis bidentatae radix, Coicis Semen


### 3.1.3. Identified Disease Target Genes

The target genes related to RA were searched in GeneCards and OMIM databases, which include 3768 genes in GeneCards and 1 gene in OMIM, with no overlapping target gene.


### 3.1.4. Intersection of Identified Targets of Compounds and Disease

In the Venn diagram intersection of identified targets about identified compounds and of RA (Figure2), a total of 77 target genes are acquired.Figure 2
A Venn diagram showing intersection identified targets of identified compounds and RA.


## 3.1.1. Identification of Compounds in SM

A total of 479 compounds were identified in SM, including 140 in Phellodendri Chinensis Cortex (Huang Bo), 49 inAtractylodes lancea (Thunb.) Dc. (Cang Zhu), 176 in achyranthis bidentatae radix (Niu Xi), and 38 in Coicis Semen (Yi Yi Ren), and in TCMSP, also including 37 in Phellodendri Chinensis Cortex (Huang Bo), 26 in Atractylodes lancea (Thunb.) Dc. (Cang Zhu), 10 in achyranthis bidentatae radix (Niu Xi), and 3 in Coicis Semen (Yi Yi Ren) in BATMAN-TCM.


## 3.1.2. Selection of Compounds Using ADME Screening and Related Targets

All the identified compounds were selected through ADME screening, with 90 of 479 compounds satisfying the suggested criteria OB ≥ 30 and DL ≥ 0.18 [18–20]. Of the 90 compounds, 25 were duplicated and removed, and the remaining 65 compounds were included for further study. Moreover, of the excluded compounds that do not meet the suggested criteria, 7 compounds, including ferulic acid, beta-elemene, eugenol, and paeonol in Phellodendri Chinensis Cortex (Huang Bo) and geniposide, rutin, and astragalin in achyranthis bidentatae radix (Niu Xi), are considered bioactive compounds and included for further analysis, and the effects of ferulic acid [22], beta-elemene [23], eugenol [24–26], paeonol [27–29], geniposide [30–32], rutin [33, 34], and astragalin [35] on RA have been investigated. The final 72 compounds are selected from the four herbal medicines (Table 1). A total of 386 target genes related to the final identified compounds are obtained from the UniProt databases.Table 1
72 active compounds of SM.
Mol idMolecule nameOB (%)DLHerbMOL002636Kihadalactone A34.210.82Phellodendri Chinensis CortexMOL013352Obacunone43.290.77Phellodendri Chinensis CortexMOL002641Phellavin_qt35.860.44Phellodendri Chinensis CortexMOL002644Phellopterin40.190.28Phellodendri Chinensis CortexMOL002651Dehydrotanshinone II A43.760.40Phellodendri Chinensis CortexMOL002652delta7-dehydrosophoramine54.450.25Phellodendri Chinensis CortexMOL002656Dihydroniloticin36.430.81Phellodendri Chinensis CortexMOL002659Kihadanin A31.600.70Phellodendri Chinensis CortexMOL002660Niloticin41.410.82Phellodendri Chinensis CortexMOL002662Rutaecarpine40.300.60Phellodendri Chinensis CortexMOL002663Skimmianin40.140.20Phellodendri Chinensis CortexMOL002666Chelerythrine34.180.78Phellodendri Chinensis CortexMOL002668Worenine45.830.87Phellodendri Chinensis CortexMOL002670Cavidine35.640.81Phellodendri Chinensis CortexMOL002671Candletoxin A31.810.69Phellodendri Chinensis CortexMOL002672Hericenone H39.000.63Phellodendri Chinensis CortexMOL002673Hispidone36.180.83Phellodendri Chinensis CortexMOL000358Beta-sitosterol36.910.75Phellodendri Chinensis CortexMOL000622Magnograndiolide63.710.19Phellodendri Chinensis CortexMOL000762Palmidin A35.360.65Phellodendri Chinensis CortexMOL000787Fumarine59.260.83Phellodendri Chinensis CortexMOL000790Isocorypalmine35.770.59Phellodendri Chinensis CortexMOL001131phellamurin_qt56.600.39Phellodendri Chinensis CortexMOL001455(S)-canadine53.830.77Phellodendri Chinensis CortexMOL001771Poriferast-5-en-3beta-ol36.910.75Phellodendri Chinensis CortexMOL002894Berberrubine35.740.73Phellodendri Chinensis CortexMOL005438Campesterol37.580.71Phellodendri Chinensis CortexMOL006392Dihydroniloticin36.430.82Phellodendri Chinensis CortexMOL006401Melianone40.530.78Phellodendri Chinensis CortexMOL006413Phellochin35.410.82Phellodendri Chinensis CortexMOL006422Thalifendine44.410.73Phellodendri Chinensis CortexMOL002665Ferulic acid40.430.06Phellodendri Chinensis CortexMOL000908Beta-elemene25.630.06Phellodendri Chinensis CortexMOL000254Eugenol56.240.04Phellodendri Chinensis CortexMOL000874Paeonol28.790.04Phellodendri Chinensis CortexMOL0001792-Hydroxyisoxypropyl-3-hydroxy-7-isopentene-2,3-dihydrobenzofuran-5-carboxylic45.200.20Atractylodes lancea (Thunb.) Dc.MOL000184NSC6355139.250.76Atractylodes lancea (Thunb.) Dc.MOL000186Stigmasterol 3-O-beta-D-glucopyranoside_qt43.830.76Atractylodes lancea (Thunb.) Dc.MOL0001883β-acetoxyatractylone40.570.22Atractylodes lancea (Thunb.) Dc.MOL000088Beta-sitosterol 3-O-glucoside_qt36.910.75Atractylodes lancea (Thunb.) Dc.MOL000092daucosterin_qt36.910.76Atractylodes lancea (Thunb.) Dc.MOL000094daucosterol_qt36.910.76Atractylodes lancea (Thunb.) Dc.MOL001006Poriferasta-7,22E-dien-3beta-ol42.980.76Achyranthis bidentatae radixMOL01246128-Norolean-17-en-3-ol35.930.78Achyranthis bidentatae radixMOL012505bidentatoside, ii_qt31.760.59Achyranthis bidentatae radixMOL012537Spinoside A41.750.40Achyranthis bidentatae radixMOL012542β-ecdysterone44.230.82Achyranthis bidentatae radixMOL002714Baicalein33.520.21Achyranthis bidentatae radixMOL002776Baicalin40.120.75Achyranthis bidentatae radixMOL002897Epiberberine43.090.78Achyranthis bidentatae radixMOL003847Inophyllum E38.810.85Achyranthis bidentatae radixMOL000422Kaempferol41.880.24Achyranthis bidentatae radixMOL004355Spinasterol42.980.76Achyranthis bidentatae radixMOL012516Geniposide8.400.44Achyranthis bidentatae radixMOL000415Rutin3.200.68Achyranthis bidentatae radixMOL000561Astragalin14.030.74Achyranthis bidentatae radixMOL001323Sitosterol alpha143.280.78Coicis SemenMOL001494Mandenol42.000.19Coicis SemenMOL002372(6Z, 10 E, 14 E, 18 E)-2,6,10,15,19,23-hexamethyltetracosa-2,6,10,14,18,22-hexaene33.550.42Coicis SemenMOL002882[(2R)-2,3-dihydroxypropyl] (Z)-octadec-9-enoate34.130.30Coicis SemenMOL000359Sitosterol36.910.75Coicis SemenMOL008118Coixenolide32.400.43Coicis SemenMOL0081212-Monoolein34.230.29Coicis SemenMOL000953CLR37.870.68Coicis SemenMOL001454Berberine36.860.78Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL001458Coptisine30.670.86Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL002643Delta 7-stigmastenol37.420.75Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL000785Palmatine64.600.65Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL000098Quercetin46.430.28Phellodendri Chinensis Cortex, achyranthis bidentatae radixMOL000173Wogonin30.680.23Atractylodes lancea (Thunb.) Dc. Achyranthis bidentatae radixMOL000085Beta-daucosterol_qt36.910.75Atractylodes lancea (Thunb.) Dc. Achyranthis bidentatae radixMOL000449Stigmasterol43.830.76Phellodendri Chinensis Cortex, achyranthis bidentatae radix, Coicis Semen


## 3.1.3. Identified Disease Target Genes

The target genes related to RA were searched in GeneCards and OMIM databases, which include 3768 genes in GeneCards and 1 gene in OMIM, with no overlapping target gene.


## 3.1.4. Intersection of Identified Targets of Compounds and Disease

In the Venn diagram intersection of identified targets about identified compounds and of RA (Figure2), a total of 77 target genes are acquired.Figure 2
A Venn diagram showing intersection identified targets of identified compounds and RA.


## 3.2. Topological and PPI Network

### 3.2.1. Topological Network Analysis

The SM-RA mechanism network (Figure3) consists of 77 target genes nodes (shared gene of SM and RA), 41 compound nodes, and 349 edges. Among the 14 compounds (Dehydrotanshinone II A, Stigmasterol, beta-sitosterol, Isocorypalmine, beta-elemene, quercetin, eugenol, paeonol, (S)-Canadine, wogonin, baicalein, Inophyllum E, rutin, and kaempferol) that connected to more than four genes, 55 target genes are associated with quercetin, 22 target genes are associated with kaempferol, 15 target genes are associated with baicalein, 14 target genes are associated with wogonin, and 10 target genes are associated with beta-sitosterol and eugenol, respectively (Table 2). In addition, 10 genes, including PTGS1, ESR1, AR, PGR, CHRM3, PPARG, CHRM2, BCL2, CASP3, and RELA, are related to more than five compounds, as shown in the compounds-targets network (Figure 4). These compounds and genes may be the key nodes in the network.Figure 3
The SM-RA mechanism network. The green octagon represents rheumatoid arthritis (RA), the purple hexagon represents the herbal medicine simiao pill (SM), while pink V’s represent compounds, and blue-purple ellipses represent genes.Table 2
Target genes interacting with compounds in the SM-RA network.
CompoundsTarget genesDelta 7-stigmastenolPGRPoriferast-5-en-3beta-olPGRCampesterolPGRNSC63551PGRBeta-daucosterol_qtPGRPoriferasta-7,22E-dien-3beta-olPGR28-Norolean-17-en-3-olPGRSpinasterolPGRSitosterol alpha1PGRSitosterolPGRCLRPGRChelerythrinePTGS1AstragalinPTGS1MandenolPTGS1Ferulic acidPTGS1, CHRM2phellamurin_qtESR1, NR3C1EpiberberineESR1, ARBerberinePTGS1, ESR1, ARCoptisinePTGS1, ESR1, ARWoreninePTGS1, ESR1, ARBerberrubinePTGS1, ESR1, ARThalifendinePTGS1, ESR1, ARFumarinePTGS1, CHRM3, HTR3ARutaecarpinePTGS1, AR, HTR3A, CYP3A4CavidinePTGS1, CHRM3, HTR3A, RXRBPalmatinePTGS1, ESR1, AR, ESR23β-acetoxyatractyloneCHRM3, AR, ACHE, CHRM2Dehydrotanshinone II ACHRM3, ESR1, AR, PPARG, ACHEInophyllum EPTGS1, ESR1, AR, ESR2, GSK3BStigmasterolPGR, PTGS1, AKR1B1, PLAU, CHRM3, CHRM2IsocorypalminePTGS1, CHRM3, HTR3A, CHRM2, DRD3, RXRB(S)-canadinePTGS1, CHRM3, HTR3A, CHRM2, DRD3, RXRBPaeonolPTGS1, CHRM2, RELA, BCL2, NFKBIA, ICAM1, TYRBeta-elemeneCHRM2, PTGS1, CHRM3, BCL2, RB1, TP63, CCNB1, RHOARutinRELA, IL6, CASP3, ALOX5, GSTP1, INS, FCER2, ITGB2Beta-sitosterolPGR, PTGS1, CHRM3, CHRM2, BCL2, CASP9, CASP3, CASP8, PRKCA, PON1EugenolPTGS1, CHRM3, CHRM2, PLAU, RELA, CYP1A1, ALOX5, AHR, ABCC2, MUC1WogoninPTGS1, ESR1, AR, PPARG, GSK3B, RELA, CCND1, BCL2, CASP9, IL6, CASP3, TP63, PTGER3, MCL1BaicaleinPTGS1, AR, RELA, VEGFA, BCL2, FOS, CASP3, TP63, HIF1A, FOSL1, CCNB1, AHR, IGF2, CYCS, NOX5KaempferolPTGS1, AR, PPARG, PGR, ACHE, CHRM2, RELA, IKBKB, BCL2, CASP3, MAPK8, PPARG, CYP3A4, CYP1A1, ICAM1, SELE, VCAM1, ALOX5, GSTP1, AHR, NR1I3, GSTM1QuercetinPTGS1, AR, PPARG, AKR1B1, ACHE, RELA, EGFR, VEGFA, CCND1, BCL2, FOS, CASP9, PLAU, RB1, IL6, CASP3, TP63, NFKBIA, CASP8,RAF1, PRKCA, HIF1A,ERBB2,PPARG,CYP3A4,CAV1,MYC,CYP1A1, ICAM1, SELE, VCAM1, PTGER3, BIRC5, NOS3, HSPB1, CCNB1, ALOX5, GSTP1, NFE2L2, NQO1, PARP1, AHR, COL3A1, DCAF5, NR1I3, HSF1, CRP, RUNX2, CTSD, IGFBP3, IGF2, IRF1, ERBB3, PON1, GSTM1Figure 4
The compounds-targets network. The pink V’s represent compounds and blue-purple ellipses represent target genes.


### 3.2.2. PPI Network

The PPI network is established by setting the confidence level of more than 0.40 and hiding the independent target protein nodes. The PPI network nodes represent proteins and edges represent protein-protein interactions. The network has 75 nodes and 1604 edges (Figure5). In addition, we analyzed the importance prioritization (adjacent nodes count of each protein) of proteins according to the network, and the leading 30 genes with higher connection were visualized by constructing a bar plot diagram (Figure 6), which indicates the 30 genes or proteins that may play a bridge role in connecting other nodes in the PPI network. These 30 genes or proteins include inflammation-associated genes (IL6 [36], NFKBIA [37]), cell proliferation-, differentiation-, and transformation-related genes (FOS [38], EGFR [39], MAPK8 [40], NR3C1 [41], RHOA [42], and PARP1 [43]), cell apoptosis-related genes (CASP3, CASP8 [44], CASP9 [45], MYC [46], CYCS [47], HIF1A [48], MCL1 [49], and GSK3B [50]), cell cycle-related gene (CCND1 [51]), hormone-related genes (INS [52], ESR1 [53], AR [54], and PGR [55]), angiogenesis-related gene (VEGFA [56]), and transcription factor (RELA [57]).Figure 5
The PPI network of SM-RA. Each node represents the relevant gene, and the edges represent protein-protein associations.Figure 6
Hub top 30 genes of the PPI network. They-axis displays significant top 30 genes, and the x-axis shows line counts of these genes.


## 3.2.1. Topological Network Analysis

The SM-RA mechanism network (Figure3) consists of 77 target genes nodes (shared gene of SM and RA), 41 compound nodes, and 349 edges. Among the 14 compounds (Dehydrotanshinone II A, Stigmasterol, beta-sitosterol, Isocorypalmine, beta-elemene, quercetin, eugenol, paeonol, (S)-Canadine, wogonin, baicalein, Inophyllum E, rutin, and kaempferol) that connected to more than four genes, 55 target genes are associated with quercetin, 22 target genes are associated with kaempferol, 15 target genes are associated with baicalein, 14 target genes are associated with wogonin, and 10 target genes are associated with beta-sitosterol and eugenol, respectively (Table 2). In addition, 10 genes, including PTGS1, ESR1, AR, PGR, CHRM3, PPARG, CHRM2, BCL2, CASP3, and RELA, are related to more than five compounds, as shown in the compounds-targets network (Figure 4). These compounds and genes may be the key nodes in the network.Figure 3
The SM-RA mechanism network. The green octagon represents rheumatoid arthritis (RA), the purple hexagon represents the herbal medicine simiao pill (SM), while pink V’s represent compounds, and blue-purple ellipses represent genes.Table 2
Target genes interacting with compounds in the SM-RA network.
CompoundsTarget genesDelta 7-stigmastenolPGRPoriferast-5-en-3beta-olPGRCampesterolPGRNSC63551PGRBeta-daucosterol_qtPGRPoriferasta-7,22E-dien-3beta-olPGR28-Norolean-17-en-3-olPGRSpinasterolPGRSitosterol alpha1PGRSitosterolPGRCLRPGRChelerythrinePTGS1AstragalinPTGS1MandenolPTGS1Ferulic acidPTGS1, CHRM2phellamurin_qtESR1, NR3C1EpiberberineESR1, ARBerberinePTGS1, ESR1, ARCoptisinePTGS1, ESR1, ARWoreninePTGS1, ESR1, ARBerberrubinePTGS1, ESR1, ARThalifendinePTGS1, ESR1, ARFumarinePTGS1, CHRM3, HTR3ARutaecarpinePTGS1, AR, HTR3A, CYP3A4CavidinePTGS1, CHRM3, HTR3A, RXRBPalmatinePTGS1, ESR1, AR, ESR23β-acetoxyatractyloneCHRM3, AR, ACHE, CHRM2Dehydrotanshinone II ACHRM3, ESR1, AR, PPARG, ACHEInophyllum EPTGS1, ESR1, AR, ESR2, GSK3BStigmasterolPGR, PTGS1, AKR1B1, PLAU, CHRM3, CHRM2IsocorypalminePTGS1, CHRM3, HTR3A, CHRM2, DRD3, RXRB(S)-canadinePTGS1, CHRM3, HTR3A, CHRM2, DRD3, RXRBPaeonolPTGS1, CHRM2, RELA, BCL2, NFKBIA, ICAM1, TYRBeta-elemeneCHRM2, PTGS1, CHRM3, BCL2, RB1, TP63, CCNB1, RHOARutinRELA, IL6, CASP3, ALOX5, GSTP1, INS, FCER2, ITGB2Beta-sitosterolPGR, PTGS1, CHRM3, CHRM2, BCL2, CASP9, CASP3, CASP8, PRKCA, PON1EugenolPTGS1, CHRM3, CHRM2, PLAU, RELA, CYP1A1, ALOX5, AHR, ABCC2, MUC1WogoninPTGS1, ESR1, AR, PPARG, GSK3B, RELA, CCND1, BCL2, CASP9, IL6, CASP3, TP63, PTGER3, MCL1BaicaleinPTGS1, AR, RELA, VEGFA, BCL2, FOS, CASP3, TP63, HIF1A, FOSL1, CCNB1, AHR, IGF2, CYCS, NOX5KaempferolPTGS1, AR, PPARG, PGR, ACHE, CHRM2, RELA, IKBKB, BCL2, CASP3, MAPK8, PPARG, CYP3A4, CYP1A1, ICAM1, SELE, VCAM1, ALOX5, GSTP1, AHR, NR1I3, GSTM1QuercetinPTGS1, AR, PPARG, AKR1B1, ACHE, RELA, EGFR, VEGFA, CCND1, BCL2, FOS, CASP9, PLAU, RB1, IL6, CASP3, TP63, NFKBIA, CASP8,RAF1, PRKCA, HIF1A,ERBB2,PPARG,CYP3A4,CAV1,MYC,CYP1A1, ICAM1, SELE, VCAM1, PTGER3, BIRC5, NOS3, HSPB1, CCNB1, ALOX5, GSTP1, NFE2L2, NQO1, PARP1, AHR, COL3A1, DCAF5, NR1I3, HSF1, CRP, RUNX2, CTSD, IGFBP3, IGF2, IRF1, ERBB3, PON1, GSTM1Figure 4
The compounds-targets network. The pink V’s represent compounds and blue-purple ellipses represent target genes.


## 3.2.2. PPI Network

The PPI network is established by setting the confidence level of more than 0.40 and hiding the independent target protein nodes. The PPI network nodes represent proteins and edges represent protein-protein interactions. The network has 75 nodes and 1604 edges (Figure5). In addition, we analyzed the importance prioritization (adjacent nodes count of each protein) of proteins according to the network, and the leading 30 genes with higher connection were visualized by constructing a bar plot diagram (Figure 6), which indicates the 30 genes or proteins that may play a bridge role in connecting other nodes in the PPI network. These 30 genes or proteins include inflammation-associated genes (IL6 [36], NFKBIA [37]), cell proliferation-, differentiation-, and transformation-related genes (FOS [38], EGFR [39], MAPK8 [40], NR3C1 [41], RHOA [42], and PARP1 [43]), cell apoptosis-related genes (CASP3, CASP8 [44], CASP9 [45], MYC [46], CYCS [47], HIF1A [48], MCL1 [49], and GSK3B [50]), cell cycle-related gene (CCND1 [51]), hormone-related genes (INS [52], ESR1 [53], AR [54], and PGR [55]), angiogenesis-related gene (VEGFA [56]), and transcription factor (RELA [57]).Figure 5
The PPI network of SM-RA. Each node represents the relevant gene, and the edges represent protein-protein associations.Figure 6
Hub top 30 genes of the PPI network. They-axis displays significant top 30 genes, and the x-axis shows line counts of these genes.


## 3.3. Enrichment Analysis

### 3.3.1. GO Enrichment Analysis

GO analysis consisted of biological process (BP), cellular component (CC), and molecular function (MF). As showed in Figure7, the top 20 enrichment terms are visualized by the bar plot diagram. The results demonstrated that numerous targets are involved in various BPs associated with immune response and inflammation, such as the response to a steroid hormone, response to oxidative stress, and regulation of the apoptotic signaling pathway, which confirmed strongly the correlation with the pathogenesis in RA. The CC results showed that most of the targets are localized to the cellular membrane and nuclear chromatin part. The MF results indicated that many targets are associated with nuclear receptor activity and transcription factor activity.Figure 7
GO analysis of targets, the top 20 significant enrichment terms in BP (a), CC (b), and MF (c). They-axis shows significantly enriched biological process, cellular component, and molecular function categories of the target genes, respectively. The redder the color, the lower the P value. The x-axis displays the enrichment scores of these terms, and the length of the bar indicates the number of target genes in each pathway.
(a)(b)(c)


### 3.3.2. KEGG Enrichment Analysis

The KEGG pathways are applied to examine the function and signaling pathways of the identified target genes, with the top 20 of the potential pathways (P<0.05 and q<0.05) shown by a bar plot diagram (Figure 8) and visualized with the pathways-targets-compounds network (Figure 9). The results showed that numerous targets are associated with certain virus infections (such as Epstein-Barr virus infection) and cancer, which are associated with the onset and prognosis of RA.Figure 8
KEGG enrichment pathways (top 20). They-axis displays the top 20 significantly enriched KEGG pathways of the target genes. The redder the color, the smaller the P value. The x-axis represents the target genes counts, and the length of the bar indicates the number of target genes in each pathway.Figure 9
The pathways-targets-compounds network. The green diamonds represent pathways, the blue-purple ellipses represent genes, and the pink V’s represent compounds.


## 3.3.1. GO Enrichment Analysis

GO analysis consisted of biological process (BP), cellular component (CC), and molecular function (MF). As showed in Figure7, the top 20 enrichment terms are visualized by the bar plot diagram. The results demonstrated that numerous targets are involved in various BPs associated with immune response and inflammation, such as the response to a steroid hormone, response to oxidative stress, and regulation of the apoptotic signaling pathway, which confirmed strongly the correlation with the pathogenesis in RA. The CC results showed that most of the targets are localized to the cellular membrane and nuclear chromatin part. The MF results indicated that many targets are associated with nuclear receptor activity and transcription factor activity.Figure 7
GO analysis of targets, the top 20 significant enrichment terms in BP (a), CC (b), and MF (c). They-axis shows significantly enriched biological process, cellular component, and molecular function categories of the target genes, respectively. The redder the color, the lower the P value. The x-axis displays the enrichment scores of these terms, and the length of the bar indicates the number of target genes in each pathway.
(a)(b)(c)


## 3.3.2. KEGG Enrichment Analysis

The KEGG pathways are applied to examine the function and signaling pathways of the identified target genes, with the top 20 of the potential pathways (P<0.05 and q<0.05) shown by a bar plot diagram (Figure 8) and visualized with the pathways-targets-compounds network (Figure 9). The results showed that numerous targets are associated with certain virus infections (such as Epstein-Barr virus infection) and cancer, which are associated with the onset and prognosis of RA.Figure 8
KEGG enrichment pathways (top 20). They-axis displays the top 20 significantly enriched KEGG pathways of the target genes. The redder the color, the smaller the P value. The x-axis represents the target genes counts, and the length of the bar indicates the number of target genes in each pathway.Figure 9
The pathways-targets-compounds network. The green diamonds represent pathways, the blue-purple ellipses represent genes, and the pink V’s represent compounds.


## 3.4. Molecular Docking Analysis

The selected targets, including IL6, VEGFA, EGFR, and NFKBIA, play a significant role in the SM-RA network. The candidate compounds, including quercetin, kaempferol, baicalein, and wogonin, are the top 4 compounds (ranking by related target genes count) in the SM-RA network. These 4 target genes and 4 compounds are imported into MOE for molecular docking verification. The docking scores are shown in Table3. The action mode of NFKBIA and quercetin, kaempferol, baicalein, and wogonin and the action mode of wogonin and IL6, VEGFA, EGFR, and NFKBIA are shown in Figure 10.Table 3
Molecular docking scores.
IL6VEGFAEGFRNFKBIAQuercetin−4.7051−5.9131−6.0857−6.7291Kaempferol−4.9898−5.4844−5.7466−6.5524Baicalein−4.5010−5.5058−5.8372−6.2222Wogonin−4.6678−5.6466−6.1084−6.6169Figure 10
Molecular docking results. (a) The action mode of NFKBIA and quercetin, kaempferol, baicalein, and wogonin: (A) NFKBIA and quercetin; (B) NFKBIA and kaempferol; (C) NFKBIA and baicalein; (D) NFKBIA and wogonin. (b) The action mode of wogonin and IL6, VEGFA, EGFR, and NFKBIA: (A) wogonin and IL6; (B) wogonin and VEGFA; (C) wogonin and EGFR; (D) wogonin and NFKBIA.
(a)(b)


## 4. Discussion

In the present network pharmacological analysis, a total of 479 compounds were identified in the four herbal medicines of SM, and 72 compounds were yielded by ADME criteria screening. A total of 386 targets related to potential compounds and 3769 targets associated with RA were identified, and 77 target genes were obtained from the interaction of targets about SM identified compounds and RA. SM-RA network analysis visualized the interaction of multicomponents and multitargets about SM on RA. The compounds-targets network analysis indicated that the 6 compounds, including quercetin, kaempferol, baicalein, wogonin, beta-sitosterol, and eugenol, were linked to ≥10 target genes, and the 10 target genes (PTGS1, ESR1, AR, PGR, CHRM3, PPARG, CHRM2, BCL2, CASP3, and RELA) were core target genes in the network. GO enrichment analysis indicated that numerous targets are involved in response to a steroid hormone, oxidative stress, and regulation of the apoptotic signaling pathway in BP, are localized to the cellular membrane and nuclear chromatin part in CC, and are associated with nuclear receptor activity and transcription factor activity in MF. KEGG pathways analysis results indicated that numerous targets are associated with certain virus infections and cancer. Molecular docking showed that quercetin, kaempferol, baicalein, and wogonin have good binding activity with IL6, VEGFA, EGFR, and NFKBIA targets.About 72 identified compounds, particularly the 6 compounds, including quercetin, kaempferol, baicalein, wogonin, beta-sitosterol, and eugenol, were linked to more than 10 targets, indicating that these compounds might play a vital role in the process of RA treatment. Furthermore, certain compounds have exhibited the potential antirheumatic therapeutic activities except for wogonin (Table4). For instance, quercetin has been reported to decrease levels of tumor necrosis factor-α (TNF-α), interleukin-1β (IL-1β), interleukin-17 (IL-17), and monocyte chemotactic protein-1 (MCP-1) [58] and significantly reduced damage to interchondral joints, infiltration of inflammatory cells, and pannus formation [59]. Besides, kaempferol suppresses the proliferation and migration of RAFLS and the release of activated T-cell-mediated inflammatory cytokines and reduces osteoclast differentiation through targeting on the fibroblast growth factor receptor 3- (FGFR3-) ribosomal S6 kinase 2 (RSK2) signaling axis [60]. In addition, baicalein inhibits human rheumatoid arthritis fibroblast-like synoviocytes (RAFLS) proliferation involving suppression of nuclear factor kappa B (NF-κB) transcriptional activity and recombinant macrophage migration inhibitory factor- (MIF-) mediated signaling [61]. What’s more, β-Sitosterol could modulate the functions of macrophages and attenuates rheumatoid inflammation in CIA mice [62]. For eugenol, it is reported to be effective in ameliorating oxidative stress and inflammation in arthritic rats [25, 26]. Moreover, among the other 66 compounds, some articles previously reported the antirheumatic effect. For example, ferulic acid is reported to suppress osteoclast differentiation and bone erosion via the inhibition of receptor activator of nuclear factor кB ligand- (RANKL-) dependent NF-κB signaling pathway [63], and berberine could attenuate adjuvant-induced arthritic fibroblast-like synoviocytes (AA-FLS) proliferation and regulate the Th17/Treg imbalance [64]. Collectively, these active components exhibit antirheumatic effects from various aspects, including anti-inflammatory, immunoregulatory, reducing bone erosion and destruction, and attenuating oxidative stress. Therefore, these might indicate the collective effectiveness and diversity of constituents in SM for treating RA.Table 4
Potential anti-RA mechanisms of some compounds.
CompoundMechanismModelReferenceQuercetinDecreased TNF-α, IL-1β, IL-17, and MCP-1CIA miceHaleagrahara et al. [58]Decreased TNF-α in joints, reduced interchondral joints damage, inflammatory cells infiltration, and pannus formationCIA miceKawaguchi et al. [59]Promote RAFLS apoptosis by upregulating lncRNA metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) and inhibiting PI3K/AKT signal activation subsequentlyRAFLSPan et al. [72]Exerted anti-inflammatory, analgesic, and antioxidant effects by inhibiting NF-κB and regulating nuclear factor erythroid 2-related factor (Nrf2)/home oxygenase (HO-1) signalZymosan-induced arthritis miceGuazelli, et al. [73]Inhibited IL-17 and RANKL production, suppressed Th17 cellRAFLSKim HR, et al. 2019 [74]KaempferolInhibited RAFLS proliferation and migration, suppressed inflammatory cytokines (IL-17, IL-21, and TNF-α) by targeting FGFR3-RSK2 signalRAFLSLee, et al. [60]Inhibited RAFLS migration and invasion by blocking MAPK signalRAFLSPan et al. [75]Inhibited RAFLS proliferation, reduced MMPs, COX-2, and PGE2 production, inhibited NF-κB activationRAFLSYoon et al. [76]BaicaleinInhibited RAFLS proliferation by suppressing NF-κB activationRAFLSChen et al. [61]Beta-sitosterolInhibited inflammatory cytokines (iNOS, IL-1β), modulated macrophages functionsCIA miceLiu et al. [62]EugenolInhibited mononuclear infiltration, lowered TNF-α, TGF-β, and IFN-γCIA murineGrespan et al. [24]Reduced inflammatory cytokines (TNF-α, IL-6, and IL-10) and oxidative stressCIA ratMateen et al. [25]Reduced inflammatory cytokines (TNF-α, IL-6) and oxidative stressRA patientsMateen et al. [26]The anti-RA effect of identified compounds (quercetin, kaempferol, baicalein, beta-sitosterol, and eugenol) is partially associated with the potential target genes, including NFKBIA, IL6, and MAPK, and potential signals, including PI3K‐AKT, TNF, and IL‐17, indicating the interaction between multicomponents, multitargets, and multisignaling of SM treating RA.Among the main target genes (top 30) in the PPI network is INS, ranking first with the highest connection, which may affect the local inflammatory process of joint in RA [52], though research about the role of INS in RA is rare. IL6 is involved in the regulation of the immune response, inflammation, and hematopoiesis and confirmed the pathological roles in RA [65]. VEGFA contributes to promoting the angiogenic phenotype of RA [56]. EGFR is proved to be involved in the proliferation and cytokine production of synovial fibroblasts, the proliferation of endothelial cells, and the formation of osteoclasts [39]. CASP3, CASP8 [44], and CASP9 [45] are involved in the apoptosis of RA synoviocytes. NFKBIA is related to the inflammation of RA by regulating many genes for immune response, cell adhesion, differentiation, proliferation, angiogenesis, and apoptosis [37]. The antirheumatic effect of aforementioned baicalein, ferulic acid, etc. is partially associated with these target genes, indicating the interaction between multicomponents and multitargets of SM treating RA.KEGG pathway enrichment analysis indicated that certain types of virus infection and cancer might also be crucial in the network. The evidence that viral infection contributes to RA, such as Epstein-Barr virus infection [1], is strong, and RA is associated with an increased risk of cancer [66]. In addition, the KEGG pathway analysis also indicated that PI3K-Akt, TNF, and IL-17 signaling pathway may be a critical signaling pathway in the network pharmacology. The PI3K-Akt signaling pathway is involved in inflammatory cytokine production [67], proliferation and migration of RAFLS [68] and chondrocyte proliferation [69], and apoptosis and autophagy in RA [69]. Moreover, a pivotal role for the proinflammatory cytokines, including tumor necrosis factor (TNF) [70] and interleukin- (IL-) 17 [71], in RA joint pathology has been identified.Of the leading 30 target genes with a higher connection in the PPI network, IL6, VEGFA, EGFR, and NFKBIA play a critical role in the development of RA, which has been aforementioned. Besides, in the visualized pathways-targets network, IL6, VEGFA, EGFR, and NFKBIA are involved in numerous pathways, indicating that SM may exert anti-RA effects through multipathways and multitargets combined interaction. Furthermore, the molecular docking analysis was constructed to investigate the interaction of some candidate compounds and targets. For example, the absolute value of docking scores about NFKBIA and quercetin, kaempferol, baicalein, and wogonin is the highest in each group, indicating that NFKBIA has a higher binding affinity than other target genes. For wogonin, although there have been no relevant studies about the effect in RA, the docking results indicated that wogonin performed good binding activity with IL6, VEGFA, EGFR, and NFKBIA. In brief, the high binding affinities of these active components indicated that the therapeutic effects of SM treating RA were probably through the modulation of several related targets.As shown, the anti-RA effect of identified compounds (quercetin, kaempferol, baicalein, beta-sitosterol, and eugenol) is partially associated with the potential target genes, including NFKBIA, IL6, and MAPK, and potential signals, including PI3K-Akt, TNF, and IL-17, indicating the interaction between multicomponents, multitargets, and multisignaling of SM treating RA.


## 5. Limitation

This study has some limitations. It provides only a predictive overview of the pharmacological mechanisms of SM against RA based on the existing database, and further experiment verification in vivo and in vitro is necessary to ensure the reliability and reasonability of predicted results. First, posttranscriptional processing, translation regulation, and posttranslational processing and regulation play a critical role in gene expression regulation, and most of the research on mechanisms about SM treating RA is gene level in this study; therefore, an in-depth study needs to explore the related mechanism. Second, the key proteins and KEGG pathways need to be verified. Third, the anti-RA effect needs to be further verified in the animal model.Besides, for clinical applying of SM treating RA, owing to the ethnic, genetic, and possible etiological differences, potential mechanisms about the related therapeutic module coinciding with clinical applications are worthy of further experimental investigation. In addition, importantly, dosage exploration, oral bioavailability, water-solubility, pharmacokinetics, and potential side effects of SM will also need a thorough exploration.


## 6. Conclusion

In summary, a bioinformatics/topology-based strategy, including ADEM screening, bioinformatics, network topology, enrichment analysis, and molecular analysis, was applied for identification of the molecular mechanisms of SM against RA. The integrated strategy might make the decipherment of biological mechanisms more accurate and efficient. The SM-RA network, compounds-targets network, and pathways-targets network analysis visualized the interaction of multicomponents and multitargets about SM treating RA. In particular, quercetin, kaempferol, baicalein, wogonin, beta-sitosterol, and eugenol might be the candidate therapeutic agents, and PTGS1, ESR1, AR, PGR, CHRM3, PPARG, CHRM2, BCL2, CASP3, and RELA were identified as potential drug targets. The enrichment and PPI analysis revealed the biological functions of the grouping networks related to the pathogenesis of RA. The multicomponent cosynergism of the herbal combinations about SM was elaborated. The study also revealed the multifunctional synergetic mechanisms of SM, including certain virus infection and cancer, and PI3K-Akt, TNF, and IL-17 signaling pathway.



---

*Source: 5786053-2020-11-04.xml*