Integrated Molecular Docking with Network Pharmacology to Reveal the Molecular Mechanism of Simiao Powder in the Treatment of Acute Gouty Arthritis

Fan, Yihua; Liu, Wei; Jin, Yue; Hou, Xu; Zhang, Xuewu; Pan, Hudan; Lu, Hang; Guo, Xiaojing

doi:https://doi.org/10.1155/2021/5570968

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Methods Results Discussion Conclusion Abbreviations Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Treatment of Musculoskeletal Diseases with Chinese Medicine

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 5570968 | https://doi.org/10.1155/2021/5570968

Integrated Molecular Docking with Network Pharmacology to Reveal the Molecular Mechanism of Simiao Powder in the Treatment of Acute Gouty Arthritis

Yihua Fan,^1,2Wei Liu ,^1,2Yue Jin,^1,2Xu Hou,³Xuewu Zhang,⁴Hudan Pan,⁵Hang Lu,^1,2and Xiaojing Guo^1,2

Academic Editor: Jun Jiang

Received14 Jan 2021

Revised12 Mar 2021

Accepted16 Apr 2021

Published28 Apr 2021

Abstract

Background. The incidence of gout has been rapidly increasing in recent years with the changing of diet. At present, modern medications used in the clinical treatment of gout showed several side effects, such as gastrointestinal damage and the increased risk of cardiovascular disease. The traditional Chinese prescription Simiao Powder (SMP) has a long history in the treatment of acute gouty arthritis (AGA) and has a good curative effect. However, the mechanism and target of its therapeutic effects are still not completely understood. Methods. Potential active compounds (PACs) and targets of SMP were found in the TCMSP database, and the disease target genes related to AGA were obtained by searching CTD, DisGeNET, DrugBank, GeneCards, TTD, OMIM, and PharmGKB disease databases with “acute gouty arthritis” and “Arthritis, Gouty” as keywords, respectively. The network of “Traditional Chinese medicine (TCM)-PACs-potential targets of acute gouty arthritis” was constructed with the Cytoscape 3.7.2 software, and the target genes of acute gouty arthritis were intersected with genes regulated by active compounds of SMP. The resultant common gene targets were input into Cytoscape 3.7.2 software, and the BisoGenet plug-in was used to construct a PPI network. The GO functional enrichment analysis and KEGG pathway enrichment analysis of the intersecting target proteins were performed using R software and corresponding program packages. The molecular docking verification was carried out between the potentially active compounds of SMP and the core target at the same time. Results. 40 active components and 203 targets were identified, of which 95 targets were common targets for the drugs and diseases. GO function enrichment analysis revealed that SMP regulated several biological processes, such as response to lipopolysaccharide and oxidative stress, RNA polymerase II transcription regulator complex, protein kinase complex, and other cellular and molecular processes, including DNA-binding transcription factor binding. Results of KEGG pathway analysis showed that SMP was associated with AGA-related pathways such as interleukin-17 (IL-17), tumor necrosis factor (TNF), p53, and hypoxia-inducible factor 1 (HIF-1) signaling pathways. The results of molecular docking showed that active compounds in SMP exhibited strong binding to five core protein receptors (TP53, FN1, ESR1, CDK2, and HSPA5). Conclusions. Active components of SMP, such as quercetin, kaempferol, wogonin, baicalein, beta-sitosterol, and rutaecarpine, showed therapeutic effects on AGA. These compounds were strongly associated with core target proteins (such as TP53, FN1, ESR1, CDK2, and HSPA5). This study reveals that IL-17, TNF, p53, and HIF-1 signaling pathways mediate the therapeutic effects of SMP on AGA. These findings expand our understanding of the mechanism of SMP in the treatment of AGA.

1. Introduction

Gouty arthritis (GA) is characterized by acute onset of single joint redness, swelling, heat, and pain at night or in the morning. The pain progressively worsens, and in severe cases, patients often suffer from cutting pain, which often reaches a peak within 24 hours [1]. The number of gout patients has increased drastically due to dietary changes and high intake of high-purine foods and beer [2]. According to the gout epidemic survey in Hong Kong, the prevalence of gout increased from 1.56% in 2006 to −2.92% in 2016 [3]. In the United States, the prevalence of gout was about 3.9% in 2016, with an estimated 9.2 million gout patients [4]. Acute gouty arthritis (AGA) is the first symptom of gout. Severe pain can reduce the quality of life of patients. Recurrent acute arthritis can destroy joints and damage internal organs [5]. Gout treatment guidelines issued in 2012 by the American College of Rheumatology (ACR) recommend nonsteroidal anti-inflammatory drugs (NSAIDs) or oral colchicine for acute attacks of gout [6]. It has been reported that NSAIDs are harmful to the gastrointestinal tract, liver, and kidney, central nervous system, etc. [7]. Colchicine is an anti-inflammatory drug commonly used in the treatment of acute gout. However, when taken in large doses, it easily causes gastrointestinal tract discomfort, liver and kidney damage, or heart disease. Meanwhile, it may lead to poisoning in elderly patients [8]. Therefore, safer and more effective drugs for the treatment of AGA are urgently needed.

Simiao Powder (SMP) is a well-known prescription used in ancient China for the treatment of joint pain [9]. This prescription originated from Chengfang Biandu in the Qing Dynasty (1904 A.D.). SMP is composed of Rhizoma Atractylodis (Cang Zhu), Cortex Phellodendri (Huang Bo), Radix Vladimiriae (Niu Xi), and Semen Coicis (Yi Yi Ren), and it could clear heat, dry dampness, and promote blood circulation to invigorate the tendons. In the gouty arthritis model, Lin et al. revealed that SMP treatment could effectively relieve joint symptoms and improve laboratory indexes (e.g., serum uric acid content or inflammatory cytokines IL-1b, IL-9, IFN-g, and so on). Moreover, the anti-inflammatory effect of SMP was stronger than febuxostat [10]. A systematic review and meta-analysis demonstrated that SMP possesses anti-inflammatory, analgesic, and uric acid-lowering effects and hence can treat AGA. Jiawei reported that SMP is superior to other anti-inflammatory drugs or combined uric acid-lowering drugs. Compared to modern medications, SMP is associated with fewer adverse events [11]. However, the mechanism and the therapeutic target of SMP in the management of AGA are vague and are the main limiting factors for its widespread clinical application [11].

The entire human body is a complex network, and the occurrence and progression of a disease are the results of the comprehensive action of multiple factors. Single-targeted therapy cannot meet the needs of multitarget and multipathway treatment of diseases. The network pharmacology concept put forward by Hopkins in 2007 [12] provides a “multiway” and “multitarget” approach to drug analysis. It is similar to the thoughts of TCM compound prescription, as well as the characteristics of the overall TCM compatibility and comprehensive and multipathway treatment of diseases.

The present study employed network pharmacology and molecular docking techniques to establish the active components of SMP in AGA treatment and predicted the potential target and drug action pathways of SMP. Collectively, we aimed to explore the potential mechanism of action of SMP and provide a reference for pharmacological research and clinical application of this prescription. The specific flow chart is shown in Figure 1.

2. Methods

2.1. Screening of Drug Components and SMP Targets

Rhizoma Atractylodis (Cang Zhu), Cortex Phellodendri (Huang Bo), Radix Vladimiriae (Niu Xi), and Semen Coicis (Yi Yi Ren), SMP components, were acquired from the traditional Chinese medicine systems pharmacology database and analysis platform (TCMSP, http://tcmspw.com/tcmsp.php), capturing all the potential active components of drugs in SMP. We screened the active compounds in SMP following the standards of oral bioavailability (OB) ≥30% and drug-likeness (DL) ≥0.18. Then, a search was undertaken for corresponding targets of active compounds on the TCMSP platform. The identified drug targets were input into the UniProt database (https://www.uniprot.org/). To obtain the abbreviation of the human gene for each target, we set the protein type as “Homo sapiens.”

2.2. Search and Screening of AGA-Related Genes

The search for target genes related to AGA was achieved using the terms “acute gouty arthritis” and “Arthritis, Gouty” in the Human Gene Database (GeneCards, https://www.genecards.org/), Online Mendelian Inheritance in Man (OMIM, https://www.omim.org/), the Pharmacogenomics Knowledge Base database (PharmGKB, https://www.pharmgkb.org/), DrugBank database (https://www.drugbank.ca/), DisGeNET database (https://www.disgenet.org/), the Therapeutic Target Database (TTD, http://db.idrblab.net/ttd/), and Comparative Toxicogenomics Database (CTD, http://ctdbase.org/). The results of the above databases were integrated to obtain the target genes of AGA after deleting the repeated genes.

2.3. Potential SMP Targets in the Treatment of AGA

We intersected the targets regulated by active compounds in SMP with the disease target genes of AGA. After that, common targets were obtained. The R 4.0.2 software and its Venn diagram package were employed to generate a Venn diagram of the SMP-AGA-related target.

2.4. Network Construction of Potential Active Compound-Potential Target

The potential targets and active compounds corresponding to the potential targets of SMP in AGA treatment were introduced into Cytoscape 3.7.2 software. Then, we constructed the network of “TCM-active compound-potential targets” and conducted the visualization display.

2.5. Construction of the Protein-Protein Interaction (PPI) Network and Screening of Core Targets

The intersection target was imported into Cytoscape 3.7.2 software. The PPI network was constructed via its plug-in BisoGenet with the following data sources: DIP, BIOGRID, HPRD, BIND, MINT, and INTACT. The plug-in CytoNCA was employed for topology analysis. Then, the degree centrality (DC) betweenness centrality (BC) of each node was calculated. A larger DC and BC of the node at the protein location implied that the protein was more important in the constructed PPI network [13]. Because we screened out various important target proteins, other protein interaction parameters, such as closeness centrality (CC) and eigenvector centrality (EC), were further applied to screen core target proteins.

2.6. Gene Ontology (GO) Functional Enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Enrichment Analyses

The ID of the intersection target was obtained using R4.2.0 software and its “org.Hs.eg.db” package. Then, the “clusterProfiler” and other program packages were adopted for GO functional enrichment analysis of the target. A histogram was constructed after screening out the first 10 functional categories of biological process (BP), cellular component (CC), and molecular function (MF). The target proteins were subjected to the KEGG enrichment analysis, whereby the first 30 KEGG pathways were selected to construct the histogram. Information of the first 20 KEGG pathways was introduced into Cytoscape 3.7.2 to draw the KEGG pathway-gene network map. After that, we selected a signal pathway closely related to GA for visual display.

2.7. SMP Potential Active Compounds and Protein Molecular Docking

The 2D structure of the potential compound was retrieved from PubChem (https://pubchem.ncbi.nlm.nih.gov/). For small molecules, their energy was minimized using ChemBio3D 2014 software. The 2D structure was converted into a 3D structure. We downloaded the 3D structure of the core target protein from the Protein Data Bank (PDB) database (http://www.rcsb.org). Using the Pymol software, the target protein receptor molecules were processed, including dehydration and removal of ligand small molecules. The target protein receptor molecule was hydrogenated using Autodock Tools 4.2.6 software. The center coordinate and size of the box were set based on the position of the active site of the protein molecule and the area where it potentially acted on the ligand small molecule. Molecular docking was achieved using AutoDock Vina, whereby lower binding energy depicted a better affinity between the receptor and the ligand. The binding energy ≤0 kcal/mol indicated that the compound could bind and interact with the target, whereas the binding energy <−5 kcal/mol demonstrated a very strong binding force. The docking results with the best binding force between each core target protein and active compound were simultaneously presented in 2D images and 3D structures.

3. Results

3.1. Target Screening of SMP and AGA

In total, 403 compounds of SMP were obtained from the TCMSP database, including 49 from Rhizoma Atractylodis (Cang Zhu), 140 from Cortex Phellodendri (Huang Bo), 176 from Radix Vladimiriae (Niu Xi), and 38 from Semen Coicis (Yi Yi Ren). Besides, 52 active compounds were screened according to OB ≥ 30% and DL ≥ 0.18, including 4 from Rhizoma Atractylodis (Cang Zhu), 25 from Cortex Phellodendri (Huang Bo), 17 from Radix Vladimiriae (Niu Xi), and 6 from Semen Coicis (Yi Yi Ren). After eliminating repeated active compounds, 40 active compounds were obtained. From the TCMSP platform, 916 targets of the 40 active compounds were obtained, including 60 from Rhizoma Atractylodis (Cang Zhu), 430 from Cortex Phellodendri (Huang Bo), 384 from Radix Vladimiriae (Niu Xi), and 42 from Semen Coicis (Yi Yi Ren). After eliminating repeated targets, 203 targets and their abbreviations were obtained. Information of active compounds and the number of action targets of TCM in SMP are outlined in Table 1. Through a search in GeneCards, OMIM, PharmGKB, DrugBank, DisGeNET, TTD, and CTD databases, 1,204 target genes associated with AGA were obtained (Figure 2).

3.2. Potential Targets and Corresponding Active Compounds of SMP in the Treatment of AGA

The drug-disease intersection target (the potential target of SMP in AGA treatment) was obtained by intersecting the 203 drug targets with 1204 disease targets using R software, with 95 potential targets (Figure 3). Notably, 33 active compounds corresponded to the 95 potential targets. They were the active compounds of SMP in AGA treatment.

3.3. Construction of TCM-PACs-Potential Target of AGA Network

We drew and analyzed the potential targets of SMP in the treatment of AGA using Cytoscape 3.7.2 network drawing software. Then, a PAC-potential target network diagram of SMP for AGA treatment was constructed (Figure 4). The network diagram comprised 133 nodes (including 95 target genes, 33 drug active components, 4 TCM, and 1 disease name) and 250 edges. We further calculated the degree of the compound in the figure. Notably, a higher degree implied that the compound played a more critical role in the network. The top six active compounds included MOL000098-quercetin, MOL000422-kaempferol, MOL000173-wogonin, MOL002714-baicalein, MOL000358beta-sitosterol, and MOL002662-rutaecarpine, with 79, 30, 21, 18, 12, and 7 as the corresponding target numbers, respectively.

Figure 4

A network of the drug-PAC-potential target of AGA. The left circle represents the active ingredients. Different colors distinguish drug sources of each ingredient; blue represents the active ingredient from Cortex Phellodendri (Huang Bo), yellow represents the active ingredient from Rhizoma Atracylodis (Cang Zhu), orange represents the active ingredient from Semen Coicis (Yi Yi Ren), green represents the active ingredient from Radix Vladimiriae (Niu Xi), and purple represents the common ingredient of various drugs. The square node on the right denotes the target.

3.4. Construction of PPI Network, Topological Analysis, and Determination of Core Targets

After obtaining the potential target genes, they were imported into Cytoscape 3.7.2 software. We then drew the PPI network diagram using the Bisogent plug-in and obtained 5,264 nodes and 133,011 connections (Figure 5(a)). For the first topology analysis, DC ≥ 74 and 982 nodes and 43,726 connections were obtained (Figure 5(b)). The second topology analysis was conducted with BC ≥ 925.611 (average), CC ≥ 0.581 (average), and EC ≥ 0.096 (average) to identify the crucial target genes; eventually, 54 nodes and 829 lines were obtained. Target proteins with the top five degree values were selected as the core target proteins in the protein network (Figure 5(c)), where TP53 protein, FN1 protein, ESR1 protein, CDK2 protein, and HSPA5 protein corresponded to 946 edges, 781 edges, 1,762 edges, 687 edges, and 445 edges, respectively.

3.5. GO Functional Enrichment and KEGG Pathway Enrichment Analyses

For GO analysis, BP, CC, and MF were included, and 1,938 items were obtained via BP enrichment analysis (including response to lipopolysaccharide, oxidative stress, reactive oxygen species (ROS), etc.); 36 items via CC enrichment analysis (including RNA polymerase II transcription regulator complex, protein kinase complex, and membrane raft); 145 items via MF enrichment analysis (including DNA-binding transcription factor binding, heme binding, and RNA polymerase II-specific DNA-binding transcription factor binding). To draw the GO function histogram of SMP in AGA treatment, we selected the top 10 BP, CC, and MF results at (Figure 6).

To further elucidate the pathways regulated by the therapeutic target genes, we did a KEGG pathway analysis. Results revealed that these target genes were distributed in 155 pathways. To construct a histogram, the pathways with the top 30 enriched genes were selected (Figure 7). GA-related pathways were mainly associated with IL-17, TNF, p53, and HIF-1 signaling pathways. We also analyzed the association of the top 20 pathways with their corresponding target genes by constructing a KEGG pathway-gene network diagram (Figure 8). There were 83 nodes (including 20 pathways and 63 target genes) and 403 lines. Based on the number of regulatory pathways, the top five genes were RELA, TP53, MAPK1, IKBKB, and CHUK with 16, 15, 15, 15, and 15 regulatory pathways, respectively. A strong relationship existed between GA and inflammation (the IL-17 pathway is closely related to inflammation). Therefore, we constructed a diagram to demonstrate the association of the regulatory target of SMP with the IL-17 pathway (Figure 9).

3.6. Molecular Docking

Using the AutoDock Vina software, PACs in the network map of “TCM-PAC-potential targets of AGA” were docked with five core target proteins (TP53, FN1, ESR1, CDK2, and HSPA5) in PPI. A higher binding activity between the compound and the target protein receptor implied lower binding energy. Results demonstrated that the docking of 33 active compounds with the 5 core target protein receptors was less than 0 kcal/mol (a majority was less than −5 kcal/mol) (Table 2). The molecular docking results of the compounds with the highest binding energy corresponding to the five targets are displayed in 2D images and 3D structures (Figure 10).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

4. Discussion

The repeated attacks and protracted course of acute gouty arthritis can lead to joint swelling, deformity, and destruction of bone and joint cartilage. Tophi can be formed in many parts of the body, and it could seriously invade joints, kidneys, and even aorta and heart valves in severe cases [14–16]. For the treatment of gout, anti-inflammatory analgesics are mainly used in the acute stage. After the pain of patients is completely relieved, the use of uric acid-lowering drugs to correct the serum uric acid level is adopted [7, 17, 18]. Studies have found that although febuxostat and allopurinol can effectively reduce blood uric acid levels, they can also increase the risk of cardiovascular disease [19]. Traditional Chinese medicine compounds could function in anti-inflammatory analgesics and lowering uric acid at the same time, which has certain advantages compared with anti-inflammatory analgesics alone [20–22]. Besides, its relatively fewer adverse reactions and lower price have attracted more and more scholars’ attention.

According to Figure 4, the key active components of SMP in the treatment of GA were quercetin, kaempferol, wogonin, baicalein, beta-sitosterol, and rutaecarpine. Quercetin is the active compound of Cortex Phellodendri (Huang Bo) and Radix Vladimiriae (Niu Xi), kaempferol is the active compound of Radix Vladimiriae (Niu Xi), wogonin is the active compound of Rhizoma Atractylodis (Cang Zhu) and Radix Vladimiriae (Niu Xi), baicalein is the active compound of Radix Vladimiriae (Niu Xi), beta-sitosterol is the active compound of Cortex Phellodendri (Huang Bo) and Radix Vladimiriae (Niu Xi), and rutaecarpine is the active compound of Cortex Phellodendri (Huang Bo). Among them, quercetin and wogonin are the common active components of two traditional Chinese medicines, which is similar to the compatibility and synergism of TCM compounds. Considering Semen Coicis (Yi Yi Ren) as an indispensable TCM in this prescription, the active ingredient stigmasterol, characterized by a great degree value in Figure 4, was established as one of the crucial active compounds.

These PACs in SMP have a certain curative effect on inhibiting inflammation and lowering blood uric acid level [23]. Among them, quercetin, baicalein, wogonin, and rutaecarpine showed strong anti-inflammatory effects [19–21]. Baicalein could inhibit the activity of NF-κB and some inflammatory factors, such as monocyte chemotactic protein-1 (MCP-1), interleukins, tumor necrosis factor, and cellular adhesion molecules. In addition, it also clears the reactive oxygen species (ROS) to show anti-inflammatory and antioxidant effects [24]. Moreover, wogonin, a flavonoid from medicinal plants, has shown various biological activities. It potentially inhibits the production of inflammatory mediators by macrophages and lymphocytes [25]. Inhibition of NF-κB and mitogen-activated protein kinase (MAPK) pathways in rutaecarpine-derivative R3 inhibits the NLRP3 inflammasome activation, which consequently lowers the expression of proinflammatory cytokine IL-1β [26, 27]. The pathogenesis of GA is highly associated with the inflammatory cascade induced by the activation of the inflammatory bodies of NLRP3 [2]. Thus, the inhibitory effect of rutaecarpine on the inflammatory body activation of NLRP3 may influence GA treatment. Besides, quercetin and kaempferol may significantly elevate the total antioxidant capacity and inhibit the xanthine oxidase (XOD) effect in hyperuricemia rats [28]. In another study, inhibition of XOD activity could reduce the serum uric acid level and the formation of ROS, consequently reducing the deposition of gout urate crystals [29]. Additional studies found that quercetin could reduce serum uric acid levels in hyperuricemia mice and impede renal insufficiency. This was achieved via the regulation of renal organic ion transport protein and uromodulin, and demonstrated fewer side effects than allopurinol [30, 31]. Beta-sitosterol (SIT), the main compound among various phytosterols, can reduce the expression of chemokines and inflammatory cytokines. It then regulates various biological functions, including anti-inflammatory, analgesic, immunomodulatory, and antimicrobial functions [32, 33]. Also, stigmasterol can inhibit the inflammation induced by sodium urate crystals and XOD activity in mice, a phenomenon that validates its antigout effect [34].

Herein, to obtain the core target proteins of SMP in the treatment of GA disease, the PPI network of intersecting target proteins was constructed. Also, we performed two topological screens on the SMP target for GA treatment based on DC, BC, CC, and EC. TP53, FN1, ESR1, CDK2, and HSPA5 were the five core target proteins. Based on the existing reports, TP53, also known as p53, contributes to transcriptional regulation of the cell cycle, DNA repair, cell survival, and cell metabolism [35]. Inhibition of p53 exerts a substantial effect on promoting the senescence of IL-1β-induced chondrocytes [36]. The potential role of TP53 in the regulation of inflammation has been described, whereby it significantly inhibits the production of proinflammatory factor IL-6 [37]. The overexpression of FN1 can activate TGF-Akt/PI3K/Akt signal pathway to promote cell viability and differentiation capacity [38]. Cyclin-dependent kinase (CDK) is the main regulator of cell division. It can be potentially triggered by different cyclins at different cell cycle stages [39], and it can regulate the secretion of inflammatory cytokines in macrophages [40] and inhibit CDK2. In consequence, a reduction in chemotaxis of primary neutrophils outlines the role of CDK2 in regulating inflammation [39]. Estrogen receptor 1 (ESR1) is a subtype of human ESR, expressed in chondrocytes, stromal cells, and osteoblasts [41]. In particular, ESR can induce the transcription of related target genes to promote the proliferation and differentiation of tissue cells as they bind to estrogen [42]. Heat shock protein family A member 5 (HSPA5) gene can also participate in the inflammatory process of various types of osteoarthritis as it encodes the binding immunoglobulin protein (BiP) of the endoplasmic reticulum Hsp70 family [43].

Following the KEGG signal pathway analysis, the mechanisms of SMP in the treatment of AGA were mainly anti-inflammation and protection of the cartilage. The critical inflammatory signaling pathways included TNF and IL-17, whereas the apoptosis signaling pathways included p53 and HIF.

Current pieces of research on the mechanism of gout have reported that the phagocytosis of MSU crystal triggers the change of inflammatory cell state of synovial cells in GA pathogenesis [14]. For instance, the deposition of MSU crystals promotes the local release of TNF-α and IL-1β, which causes persistent gout episodes [44]. IL-17 also stimulates macrophages and monocytes to produce proinflammatory factors, including IL-1β, TNF-α, and IFN-γ (Figure 9). As a powerful proinflammatory cytokine, IL-17 can increase serum IL-17 levels 8 hours following the acute attack of gout [45]. Moreover, in several arthritic animal models, the inhibition of IL-17 expression was found to limit inflammation and reduce joint erosion [46]. The TNF signaling pathway has also proved to be a crucial inflammatory signal pathway. Signal transduction is mediated by TNF receptor 1 (TNFR1) and receptor-2 (TNFR2). TNFR2 is, in most cases, expressed in immune cells. It binds to TNF-α and TNF-β to regulate immune response [47]. TNF-α is a proinflammatory cytokine produced by different cell types, such as macrophages, lymphocytes, and fibroblasts, and reflects inflammation, infection, and other environmental stresses [48]. Previous reports showed that inhibition of TNF-α expression in synovial cells during AGA could effectively inhibit local inflammation [49, 50].

The physiological dose of soluble uric acid exerts a specific protective effect on cartilage [51]. However, following repeated attacks by GA, the articular cartilage edge is damaged, the articular surface becomes irregular, the joint space becomes narrow, and MSU deposition occurs [52]. A study found that reduced articular cartilage matrix secretion in patients with GA was associated with chondrocyte apoptosis [53]. p53 can regulate cell differentiation and senescence. The inhibition of p53 expression potentially improves senescence and affects the destruction of IL-1β-induced chondrocytes [36]. Compared to the inflammatory signaling pathway, the HIF-1 signaling pathway is slightly less in gout studies. As a transcription factor, activation of HIF-1α expression in synovial cells may be induced during hypoxia. Synovial hyperplasia follows the repeated infiltration of inflammatory cells in the joint. On the other hand, hypoxia is associated with immune cell infiltration and synovial hyperplasia. Therefore, the upregulation of synovium-derived HIF-1α occurs in the pathogenesis of rheumatoid arthritis (RA) and osteoarthritis (OA), which are hypoxia-associated diseases [54, 55]. HIF-1α exerts a regulatory role in the function of immune cells. Activated LPS-induced macrophages can also express HIF-1α, which is crucial in glycolysis and the induction of proinflammation [56]. The therapeutic effect of SMP on AGA may be derived through the regulation of the inflammatory signal pathway, p53 apoptosis signal, and HIF-1signaling pathway. And other top-ranked signaling pathways, such as the AGE-RAGE signaling pathway in diabetic complications, fluid shear stress, and atherosclerosis, etc., still lack their gout-related studies, which are likely to be the future direction of gout mechanism research.

According to the space matching and energy matching relationship between molecules, the core target protein of SMP in AGA treatment was docked with active compounds. This provides a reference for the subsequent research and development of targeted drugs. The 2D image of molecular docking depicts the binding mode between the core target protein and the compound as well as the interaction with the surrounding amino acid residues. Hydrophobic interaction best describes the main association of target protein TP53, FN1, and ESR1 with Inophyllum E. The target protein CDK2 residue Leu83 (A) forms a hydrogen bond with Inophyllum E, whereas the residues Asp34 (A), Lys96 (A), Thr37 (A), and Gly227 (A) of the target protein HSPA5 form five hydrogen bonds with baicalein. The main forces between ligand and protein are hydrophobic force and hydrogen bond, which both are chemical bonds with strong binding force. According to Table 2, key components of SMP have strong binding activity with core target proteins, such as TP53, FN1, ESR1, CDK2, and HSPA5, suggesting that SMP may play its pharmacological role by regulating these key targets. However, there are still some shortcomings in this study, and the conclusions need to be further verified by in vivo, in vitro, and clinical trials.

5. Conclusion

The present study systematically expounds the core compounds and molecular action mechanism of SMP in the treatment of AGA. The core compounds of SMP, including quercetin, kaempferol, wogonin, baicalein, and beta-sitosterol, were screened via network pharmacology, whereas the core targets of AGA including TP53, FN1, ESR1, CDK2, and HSPA5 were screened out via PPI. These core compounds are characterized by stable binding activity to the core targets. Also, SMP could regulate signal pathways related to AGA disease, for example, TNF, IL-17, p53, and HIF signaling pathways. These findings suggest the potential synergistic role of SMP in treating AGA via multiple components, multitargets, multiple biological functions, and multiple signal pathways.

Abbreviations

SMP:	Simiao Powder
AGA:	Acute gouty arthritis
TCM:	Traditional Chinese medicine
IL-17:	Interleukin-17
TNF:	Tumor necrosis factor
HIF-1:	Hypoxia-inducible factor 1
GA:	Gouty arthritis
ACR:	American College of Rheumatology
NSAIDs:	Nonsteroidal anti-inflammatory drugs
TCMSP:	Traditional Chinese medicine systems pharmacology database and analysis platform
OB:	Oral bioavailability
DL:	Drug-likeness
OMIM:	Online Mendelian Inheritance in Man
PharmGKB:	Pharmacogenomics Knowledge Base database
TTD:	Therapeutic Target Database
CTD:	Comparative Toxicogenomics Database
PACs:	Potential active compounds
PPI:	Protein-protein interaction
DC:	Degree centrality
BC:	Betweenness centrality
CC:	Closeness centrality
EC:	Eigenvector centrality
GO:	Gene ontology
KEGG:	Kyoto Encyclopedia of Genes and Genomes
BP:	Biological process
CC:	Cellular component
MF:	Molecular function
PDB:	Protein Data Bank
ROS:	Reactive oxygen species
MSU:	Monosodium urate
NF-κB:	Nuclear factor-kappa B
MAPK:	Mitogen-activated protein kinase
XOD:	Xanthine oxidase
SIT:	Beta-sitosterol
XO:	Xanthine oxidase
CDK:	Cyclin-dependent kinase
ESR1:	Estrogen receptor 1
HSPA5:	Heat shock protein family A (Hsp70) member 5
BiP:	Binding immunoglobulin protein
TNFR1:	TNF receptor 1
TNFR2:	TNF receptor-2
RA:	Rheumatoid arthritis
OA:	Osteoarthritis.

Data Availability

The data for this study can be provided by the corresponding author (Wei Liu: [email protected]).

Conflicts of Interest

There are no conflicts of interest between these authors.

Authors’ Contributions

Yihua Fan, Wei Liu, and Yue Jin have contributed equally to this work and are co-first authors.

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 (Nos. 81673927 and 82074377), the National Key Research and Development Project (No. 2018YFC1705201), the Tianjin Science and Technology Commission Program (No. 2019ZD12), Tianjin Key Specialty Program (20200602-1), the Special Scientific Research for Traditional Chinese Medicine (No. 201507001-07), and Qihuang Scholar of the “Hundred and Thousand” Talent Project of the Inheritance and Innovation of Traditional Chinese Medicine (Qihuang Project) (No. 201901011w).

References

R. Terkeltaub, “What makes gouty inflammation so variable?” BMC Medicine, vol. 15, no. 1, p. 158, 2017.
View at: Publisher Site | Google Scholar
T. Pascart and F. Lioté, “Gout: state of the art after a decade of developments,” Rheumatology (Oxford, England), vol. 58, no. 1, pp. 27–44, 2019.
View at: Publisher Site | Google Scholar
M. Tsoi, “Epidemiology of gout in Hong Kong: a population-based study from 2006 to 2016,” Arthritis Research & Therapy, vol. 22, no. 1, p. 204, 2020.
View at: Publisher Site | Google Scholar
M. Chen-Xu, “Contemporary prevalence of gout and hyperuricemia in the United States and decadal trends: the national health and nutrition examination survey, 2007–2016,” Arthritis & Rheumatology (Hoboken, N.J.), vol. 71, no. 6, pp. 991–999, 2019.
View at: Google Scholar
L. Staurengo-Ferrari, “Trans-chalcone attenuates pain and inflammation in experimental acute gout arthritis in mice,” Frontiers in Pharmacology, vol. 9, p. 1123, 2018.
View at: Publisher Site | Google Scholar
D. Khanna, P. P. Khanna, J. D. Fitzgerald et al., “2012 American College of Rheumatology guidelines for management of gout. Part 2: therapy and antiinflammatory prophylaxis of acute gouty arthritis,” Arthritis Care & Research, vol. 64, no. 10, pp. 1447–1461, 2012.
View at: Publisher Site | Google Scholar
S. Bindu, S. Mazumder, and U. Bandyopadhyay, “Non-steroidal anti-inflammatory drugs (NSAIDs) and organ damage: a current perspective,” Biochemical Pharmacology, vol. 180, Article ID 114147, 2020.
View at: Publisher Site | Google Scholar
Y. Jiang, “Colchicine: risk of severe toxicity and death from drug overdose,” Chinese Journal of Pharmacoepidemiology, vol. 20, no. 2, p. 102, 2011.
View at: Google Scholar
M. A.-O. Tang, “Integrating network pharmacology with molecular docking to unravel the active compounds and potential mechanism of simiao pill treating rheumatoid arthritis,” Evidence-Based Complementary and Alternative Medicine, vol. 2020, Article ID 5786053, 16 pages, 2020.
View at: Publisher Site | Google Scholar
X. Lin, “Simiao decoction alleviates gouty arthritis by modulating proinflammatory cytokines and the gut ecosystem,” Frontiers in Pharmacology, vol. 11, 2020.
View at: Publisher Site | Google Scholar
Y.-F. Liu, Y. Huang, C.-Y.-Z. Wen et al., “The effects of modified simiao decoction in the treatment of gouty arthritis: a systematic review and meta-analysis,” Evidence-Based Complementary and Alternative Medicine, vol. 2017, p. 12, 2017.
View at: Publisher Site | Google Scholar
A. L. Hopkins, “Network pharmacology,” Nature Biotechnology, vol. 25, no. 10, pp. 1110-1111, 2007.
View at: Publisher Site | Google Scholar
Y. Wang, “Drug target prediction based on the herbs components: the study on the multitargets pharmacological mechanism of qishenkeli acting on the coronary heart disease,” Evidence-Based Complementary and Alternative Medicine, vol. 2012, Article ID 698531, 10 pages, 2012.
View at: Publisher Site | Google Scholar
B. N. Cronstein and P. Sunkureddi, “Mechanistic aspects of inflammation and clinical management of inflammation in acute gouty arthritis,” JCR: Journal of Clinical Rheumatology, vol. 19, no. 1, pp. 19–29, 2013.
View at: Publisher Site | Google Scholar
M. H. Pillinger and B. F. Mandell, “Therapeutic approaches in the treatment of gout,” Seminars in Arthritis and Rheumatism, vol. 50, no. 3, pp. S24–S30, 2020.
View at: Publisher Site | Google Scholar
C. Yokose, N. Lu, H. Xie et al., “Heart disease and the risk of allopurinol-associated severe cutaneous adverse reactions: a general population-based cohort study,” Canadian Medical Association Journal, vol. 191, no. 39, pp. E1070–E1077, 2019.
View at: Publisher Site | Google Scholar
Y.-F. Huang, C. Bai, F. He, Y. Xie, and H. Zhou, “Review on the potential action mechanisms of Chinese medicines in treating coronavirus disease 2019 (COVID-19),” Pharmacological Research, vol. 158, Article ID 104939, 2020.
View at: Publisher Site | Google Scholar
X. Zheng, F. Wu, X. Lin, L. Shen, and Y. Feng, “Developments in drug delivery of bioactive alkaloids derived from traditional Chinese medicine,” Drug Delivery, vol. 25, no. 1, pp. 398–416, 2018.
View at: Publisher Site | Google Scholar
L. Li, “Insights into the action mechanisms of traditional Chinese medicine in osteoarthritis,” Evidence-based Complementary and Alternative Medicine, vol. 2017, Article ID 5190986, 13 pages, 2017.
View at: Publisher Site | Google Scholar
X. He, X. Wang, J. Fang et al., “The genus Achyranthes: a review on traditional uses, phytochemistry, and pharmacological activities,” Journal of Ethnopharmacology, vol. 203, pp. 260–278, 2017.
View at: Publisher Site | Google Scholar
H.-C. Ko, Y.-H. Wang, K.-T. Liou et al., “Anti-inflammatory effects and mechanisms of the ethanol extract of Evodia rutaecarpa and its bioactive components on neutrophils and microglial cells,” European Journal of Pharmacology, vol. 555, no. 2-3, pp. 211–217, 2007.
View at: Publisher Site | Google Scholar
W. Yan-Hang and K.-W. Zeng, “Natural products as a crucial source of anti-inflammatory drugs: recent trends and advancements,” Traditional Medicine Research, vol. 4, no. 5, p. 257, 2019.
View at: Google Scholar
H. Wang, H. Duan, S. Chen et al., “Chinese herbal medicine si-miao-san decoction for acute gouty arthritis: a protocol for systematic review and meta-analysis of randomized controlled trials,” Medicine, vol. 99, no. 32, Article ID e21510, 2020.
View at: Publisher Site | Google Scholar
L. Baandrup, S. Dinda, S. DasSharma, R. Banik, A. Chakraborty, and M. Dinda, “Therapeutic potentials of baicalin and its aglycone, baicalein against inflammatory disorders,” European Journal of Medicinal Chemistry, vol. 131, pp. 68–80, 2017.
View at: Google Scholar
D. L. Huynh, T. H. Ngau, N. H. Nguyen, G.-B. Tran, and C. T. Nguyen, “Potential therapeutic and pharmacological effects of Wogonin: an updated review,” Molecular Biology Reports, vol. 47, no. 12, pp. 9779–9789, 2020.
View at: Google Scholar
Q. Qin, H. Liu, J. Shou, Y. Jiang, H. Yu, and X. Wang, “The inhibitor effect of RKIP on inflammasome activation and inflammasome-dependent diseases,” Cellular & Molecular Immunology, vol. 18, p. 992, 2021.
View at: Google Scholar
J. Luo, X. Wang, X. Jiang et al., “Rutaecarpine derivative R3 attenuates atherosclerosis via inhibiting NLRP3 inflammasome‐related inflammation and modulating cholesterol transport,” The FASEB Journal, vol. 34, no. 1, pp. 1398–1411, 2020.
View at: Publisher Site | Google Scholar
F. Haidari, “Effects of parsley (Petroselinum crispum) and its flavonol constituents, kaempferol and quercetin, on serum uric acid levels, biomarkers of oxidative stress and liver xanthine oxidoreductase aactivity in Oxonate-induced hyperuricemic rats,” Iranian Journal of Pharmaceutical Research: IJPR, vol. 10, no. 4, pp. 811–819, 2011.
View at: Google Scholar
S. Ahmad, “Therapeutic efficacy of urinile against gouty arthritis,” Dose-Response, vol. 18, no. 4, 2020.
View at: Publisher Site | Google Scholar
Q.-H. Hu, X. Zhang, X. Wang, R.-Q. Jiao, and L.-D. Kong, “Quercetin regulates organic ion transporter and uromodulin expression and improves renal function in hyperuricemic mice,” European Journal of Nutrition, vol. 51, no. 5, pp. 593–606, 2012.
View at: Publisher Site | Google Scholar
J. X. Zhu, Y. Wang, L. D. Kong, C. Yang, and X. Zhang, “Effects of Biota orientalis extract and its flavonoid constituents, quercetin and rutin on serum uric acid levels in oxonate-induced mice and xanthine dehydrogenase and xanthine oxidase activities in mouse liver,” Journal of Ethnopharmacology, vol. 93, no. 1, pp. 133–140, 2004.
View at: Publisher Site | Google Scholar
Q. Yang, D. Yu, and Y. Zhang, “β-sitosterol attenuates the intracranial aneurysm growth by suppressing TNF-α-mediated mechanism,” Pharmacology, vol. 104, pp. 303–311, 2019.
View at: Google Scholar
S. Babu and S. Jayaraman, “An update on β-sitosterol: a potential herbal nutraceutical for diabetic management,” Biomedicine & Pharmacotherapy, vol. 131, Article ID 110702, 2020.
View at: Publisher Site | Google Scholar
C. Martins de Sá Müller, G. B. Coelho, M. Carolina de Paula Michel Araújo, and D. A. Saúde-Guimarães, “Lychnophora pinaster ethanolic extract and its chemical constituents ameliorate hyperuricemia and related inflammation,” Journal of Ethnopharmacology, vol. 242, Article ID 112040, 2019.
View at: Publisher Site | Google Scholar
E. M. Pinto and G. P. Zambetti, “What 20 years of research has taught us about the TP53 p.R337H mutation,” Cancer, vol. 126, no. 21, pp. 4678–4686, 2020.
View at: Google Scholar
H. Lu, D. Wang, H. Li et al., “GPR39 agonist TC-G 1008 ameliorates IL-1β-induced chondrocyte senescence,” Artificial Cells, Nanomedicine, and Biotechnology, vol. 47, no. 1, pp. 2612–2617, 2019.
View at: Publisher Site | Google Scholar
T. Zhang, “p53 predominantly regulates IL-6 production and suppresses synovial inflammation in fibroblast-like synoviocytes and adjuvant-induced arthritis,” Arthritis Research & Therapy, vol. 18, no. 1, p. 271, 2016.
View at: Publisher Site | Google Scholar
H. Zhang, X. Chen, P. Xue, X. Ma, and J. Li, “FN1 promotes chondrocyte differentiation and collagen production via TGF-β/PI3K/Akt pathway in mice with femoral fracture,” Gene, Article ID 145253, 2020.
View at: Google Scholar
A. Y. Hsu, D. Wang, S. Liu et al., “Phenotypical microRNA screen reveals a noncanonical role of CDK2 in regulating neutrophil migration,” Proceedings of the National Academy of Sciences, vol. 116, no. 37, pp. 18561–18570, 2019.
View at: Publisher Site | Google Scholar
P. Laphanuwat and S. Jirawatnotai, “Immunomodulatory roles of cell cycle regulators,” Frontiers in Cell and Developmental Biology, vol. 7, p. 23, 2019.
View at: Publisher Site | Google Scholar
Y. Oshima, K.-i. Matsuda, A. Yoshida, N. Watanabe, M. Kawata, and T. Kubo, “Localization of estrogen receptors alpha and beta. In the articular surface of the rat femur,” Acta Histochemica et Cytochemica, vol. 40, no. 1, pp. 27–34, 2007.
View at: Publisher Site | Google Scholar
W. Liu, F. M. Shao, L. Yan, H. X. Cao, and D. Qiu, “Polymorphisms in the gene encoding estrogen receptor alpha are associated with osteoarthritis in Han Chinese women,” International Journal of Clinical and Experimental Medicine, vol. 7, no. 12, pp. 5772–5777, 2014.
View at: Google Scholar
J. Wang, J. Lee, D. Liem, and P. Ping, “HSPA5 Gene encoding Hsp70 chaperone BiP in the endoplasmic reticulum,” Gene, vol. 618, pp. 14–23, 2017.
View at: Publisher Site | Google Scholar
Y. Zamudio-Cuevas, J. Fernández-Torres, G. A. Martínez-Nava et al., “Phagocytosis of monosodium urate crystals by human synoviocytes induces inflammation,” Experimental Biology and Medicine, vol. 244, no. 5, pp. 344–351, 2019.
View at: Publisher Site | Google Scholar
Y. Liu, Q. Zhao, Y. Yin, M. A. McNutt, T. Zhang, and Y. Cao, “Serum levels of IL-17 are elevated in patients with acute gouty arthritis,” Biochemical and Biophysical Research Communications, vol. 497, no. 3, pp. 897–902, 2018.
View at: Publisher Site | Google Scholar
M. I. Koenders, E. Lubberts, B. Oppers-Walgreen et al., “Blocking of interleukin-17 during reactivation of experimental arthritis prevents joint inflammation and bone erosion by decreasing RANKL and interleukin-1,” The American Journal of Pathology, vol. 167, no. 1, pp. 141–149, 2005.
View at: Publisher Site | Google Scholar
M. K. McCoy and M. G. Tansey, “TNF signaling inhibition in the CNS: implications for normal brain function and neurodegenerative disease,” Journal of Neuroinflammation, vol. 5, no. 1, p. 45, 2008.
View at: Publisher Site | Google Scholar
S. Abramson and Y. Yazici, “Biologics in development for rheumatoid arthritis: relevance to osteoarthritis,” Advanced Drug Delivery Reviews, vol. 58, no. 2, pp. 212–225, 2006.
View at: Publisher Site | Google Scholar
Z. Wang, X. Fan, Y. Xu, K. Chen, X. Yu, and G. Sun, “Efficacy of Xixiancao (Herba Siegesbeckiae Orientalis) on interactions between nuclear factor kappa-B and inflammatory cytokines in inflammatory reactions of rat synovial cells induced by sodium urate,” Journal of Traditional Chinese Medicine = Chung I Tsa Chih Ying Wen pan, vol. 40, no. 5, pp. 774–781, 2020.
View at: Publisher Site | Google Scholar
A. Liu, “The effects of carvacrol on inflammatory cytokines and the apoptosis of cartilage in rats with gouty arthritis,” Journal of Clinical and Experimental Medicine, vol. 16, no. 6, pp. 531–534, 2017.
View at: Google Scholar
J. Lai, “Physiological concentrations of soluble uric acid are chondroprotective and anti-inflammatory,” Scientific Reports, vol. 7, no. 1, p. 2359, 2017.
View at: Publisher Site | Google Scholar
F. Oliviero, A. Scanu, and L. Punzi, “Metabolism of crystals within the joint,” Reumatismo, vol. 63, no. 4, pp. 221–229, 2012.
View at: Google Scholar
F. Zhu, L. Yin, and L. Ji, “Suppressive effect of sanmiao formula on experimental gouty arthritis by inhibiting cartilage matrix degradation: an in vivo and in vitro study,” pp. 36–42, 2016.
View at: Google Scholar
G. George, G. L. Shyni, and K. G. Raghu, “Current and novel therapeutic targets in the treatment of rheumatoid arthritis,” Inflammopharmacology, vol. 28, no. 6, pp. 1457–1476, 2020.
View at: Publisher Site | Google Scholar
Y. Yu, W. Cai, J. Zhou et al., “Anti-arthritis effect of berberine associated with regulating energy metabolism of macrophages through AMPK/HIF-1α pathway,” International Immunopharmacology, vol. 87, Article ID 106830, 2020.
View at: Publisher Site | Google Scholar
M. Wu, “Chaetocin attenuates gout in mice through inhibiting HIF-1α and NLRP3 inflammasome-dependent IL-1β secretion in macrophages,” Archives of Biochemistry and Biophysics, vol. 670, pp. 94–103, 2019.
View at: Google Scholar

Copyright

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

PDF Download Citation

Download other formats

Order printed copies

Views

1129

Downloads

1237

Citations