Network Pharmacology and Molecular Docking-Based Prediction of the Mechanism of Qianghuo Shengshi Decoction against Rheumatoid Arthritis

Zeng, Zhihao; Hu, Jiaoting; Jiang, Jieyi; Xiao, Guanlin; Yang, Ruipei; Li, Sumei; Li, Yangxue; Huang, Huajing; Zhong, Huixian; Bi, Xiaoli

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

BioMed Research International

On this page

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

Research Article | Open Access

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

Network Pharmacology and Molecular Docking-Based Prediction of the Mechanism of Qianghuo Shengshi Decoction against Rheumatoid Arthritis

Zhihao Zeng,¹Jiaoting Hu,²Jieyi Jiang,^3,4Guanlin Xiao,^3,4Ruipei Yang,¹Sumei Li,^3,4Yangxue Li,^3,4Huajing Huang,¹Huixian Zhong,¹and Xiaoli Bi^1,3,4

Academic Editor: Mohammad Hassan Baig

Received17 Nov 2020

Revised18 Jul 2021

Accepted16 Aug 2021

Published06 Sept 2021

Abstract

Qianghuo Shengshi decoction (QHSSD) is a classical Chinese medicine formula, which is used in clinical practice for the treatment of rheumatoid arthritis (RA) in China. However, the pharmacological mechanism of QHSSD on RA has remained unclear by now. We collected and screened active compounds and its potential targets by the pharmacology platform of Chinese herbal medicines. In addition, the therapeutic targets of RA were obtained and selected from databases. Network construction analyzed that 128 active compounds may act on 87 candidate targets and identified a total of 18 hub targets. GO annotation and KEGG enrichment investigated that the action mechanism underlying the treatment of RA by QHSSD might be involved in cell proliferation, angiogenesis, anti-inflammation, and antioxidation. Finally, molecular docking verification showed that TP53, VEGFA, TNF, EGFR, and NOS3 may be related to the RA treatment and molecular dynamics simulation showed the stability of protein-ligand interactions. In this work, QHSSD might exert therapeutic effect through a multicomponent, multitarget, and multipathway in RA from a holistic aspect, which provides basis for its mechanism of action and subsequent experiments.

1. Introduction

Rheumatoid arthritis (RA) is a chronic autoimmune disease with the following symptoms: tender, warm, swollen, and stiff joints in the morning and as RA progresses, resulting in bone erosion, joint deformity, and physical disabilities [1–4]. RA affects approximately 1% of the worldwide population and mostly occurs in middle-aged women [5]. The drug treatment of RA mainly focuses on nonsteroidal anti-inflammatory drugs (NSAIDs), steroids, disease-modifying antirheumatic drugs (DMARDs), and glucocorticoids. However, these drugs can only delay the progression of RA, which may also increase the risk of infection. Long-term use of these drugs will cause drug dependence and a cascade of side effects, such as digestive tract irritation reaction, rashes, liver toxicity, and hair loss [2, 6, 7].

There is a rising trend of the application of traditional Chinese medicine in the treatment of RA, and the RA patients are more inclined to choose complementary medicine as a medical treatment [8–10]. Traditional Chinese medicine (TCM) has a long history of clinical use in China, used as complementary and alternative medicines in treatment [11]. Qianghuo Shengshi decoction (QHSSD) is a classic clinical prescription from internal and external causes (Nei Wai Shang Bian Huo Lun). Accordingly, QHSSD can dispel wind and dampness and sweat and alleviate pain, consisting of seven herbs: Notopterygii rhizoma et radix (Qianghuo, QH), Angelicae pubescentis radix (Duhuo, DH), Saposhnikoviae radix (Fangfeng FF), Ligustici rhizoma et radix (Gaoben, GB), Chuanxiong rhizoma (Chuangxiong, CX), Viticis fructus (Manjingzi, MJZ), and Glycyrrhizae radix et rhizoma (Gancao, GC). QHSSD has been used to treat influenza, allergic purpura, and inflammatory diseases, such as RA [12–15]. Unlike the single-target agent with a clear mechanism, however, the mechanism of QHSSD underlying the treatment of RA still remains a blur.

Network pharmacology is a novel approach to understand drug’s pharmacological mechanism in the network perspective and to find out the interaction between the drug and the body on the basis of biological networks [16]. There is a growing number of researches that applied the network pharmacology to elucidate the mechanism between TCM/herbs and disease treatment, involving multicompound and multitarget ones [17–24]. We therefore tried to explore the mechanism of action of QHSSD on RA by network pharmacology and to further predict the interaction modes between QHSSD and the predicted targets by molecular docking and molecular dynamics. The flowchart of our study is presented in Figure 1.

2. Materials and Methods

2.1. Collecting and Screening of Active Components of QHSSD

We collected the chemical compounds of QHSSD from the traditional Chinese medicine systems pharmacology database and analysis platform (TCMSP) (https://tcmspw.com/tcmsp.php), a common tool for network pharmacology research of TCM [25]. Oral bioavailability (OB) and druglikeness (DL) were set as the main parameters to screen the active components according to the absorption, distribution, metabolism, and excretion (ADME) process. We set and as the screening threshold values to select the candidate components through literature search [23, 26, 27].

2.2. Potential Targets of QHSSD

Considering that the target screening criteria differ between different databases (including TCMSP, PubChem, and Swiss Target Prediction databases), therefore, only the TCMSP database was applied in this work. We predicted relevant targets (DrugBank, https://go.drugbank.com/) of active components in QHSSD on the TCMSP and then transformed the target name to a standard gene name on Uniport (https://www.uniprot.org/) database and removed the duplications.

2.3. Identification of Associated Targets of RA

RA-related genes were collected from the CTD (http://ctdbase.org/) database [28] and GeneCards (https://www.genecards.org/) database [29], with the keyword “rheumatoid arthritis.” The genes from the above databases merged and removed the duplications.

2.4. Compound-Target Network Construction and Analysis

Venny 2.1.0 (https://bioinfogp.cnb.csic.es/tools/venny/) [30] was used to find out the overlapping targets between compound targets and disease targets. To explore the relationship between compounds and disease more reasonably, Cytoscape 3.7.2 (https://cytoscape.org/) [31] with a visualized tool and all node degrees of networks calculated was used to construct a compound-potential target network.

2.5. Protein-Protein Interaction (PPI) Network Construction and Analysis

The Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) (https://string-db.org/) was used to predict the interaction between proteins [32]. According to the previous steps, drug-associated genes and disease targets were analyzed on STRING to create a protein-protein interaction network (PPI network) by importing overlapping targets. Cytoscape was then used to construct a core PPI network.

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

We imported the hub targets on the Bioconductor clusterProfiler org.Hs.eg.db and DOSE packages [33, 34] of R 4.0.2 (https://cran.r-project.org/src/base/R-4/) software to create the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and construct Gene Ontology (GO) biological function enrichment analysis with and as the thresholds.

2.7. Compound-Target Molecular Docking

According to previous studies [19, 20], compounds with a higher degree value and potential targets were evaluated by molecular docking. The 3D structure format file of the active compounds screened from the previous step and the hub targets were downloaded from the PubChem (https://pubchem.ncbi.nlm.nih.gov/) database and the RCSB PDB (https://www.rcsb.org/), respectively. The active compounds are ligands and the targets are receptors. AutoDockTools 1.5.6 software was used to remove the water molecules, isolate proteins, add nonpolar hydrogen, and calculate Gasteiger charges for the structure. Both ligands and receptors were kept as PDBQT format docked through AutoDock Vina [35]. The Grid Box size and position of macromolecule were predicted by POCASA 1.1 (http://g6altair.sci.hokudai.ac.jp/g6/service/pocasa/). Then, the protein-ligand interaction was visualized with PyMOL 2.3.0 software.

2.8. Molecular Dynamics Simulations

Molecular dynamics (MD) simulations were used to evaluate the structural characteristics of the protein-ligand system [36]. In this work, the GROMACS Program package [37] and CHARMM36 force field [38] were used to simulate the molecular dynamics of the complexes in the TIP3P water model. Ions (Na⁺ or Cl⁻) were added to neutralize charges. The energy of systems was minimized through a 4000-step steepest descent method and a 2000-step conjugate gradient method to eliminate high-energy collisions between molecules. Then, the systems were heated from 0 K to 300 K within 100 ps and simulated as an NPT ensemble with normal temperature (300 K) and normal pressure (101 kPa) for 100 ps; unconstrained MD simulation was carried out for 10 ns after the system reached dynamic equilibrium. Periodic boundary conditions and particle mesh Ewald (PME) methods were used to deal with long-range electrostatic interactions. The short-range Coulomb interaction was 0.9 nm, and the SHAKE algorithm was used to restrict the bond to the hydrogen atom. The temperature and pressure of the system were controlled by V-rescale and Parrinello-Rahman. The stability of systems was evaluated using the root-mean-square deviation (RMSD) of the aligned protein-ligand coordinate set calculated against the initial frame.

3. Results

3.1. Compound-Target Network and Analysis

We collected 153 compounds and 152 putative targets of QHSSD (Table S1 and Figure S1) on the screening criteria from the TCMSP. Among these compounds, we excluded 25 compounds without RA-related targets and 128 active compounds were obtained in this work. 1846 targets were selected from 29368 genes related to RA by the screening threshold of inference on the CTD database. Meanwhile, 306 genes were obtained from 4465 targets with a relevance on GeneCards. Targets from the above database merged and removed the overlapping parts, and at last, we received a total of 1992 RA-related targets. Based on the Venn diagram, we obtained 87 overlapping targets from 152 putative targets of QHSSD and 1992 RA-related targets. Consequently, 87 targets (Figure 2 and Table 1) were considered as potential targets for the treatment. To further investigate the interaction between the compounds and targets, we constructed a compound-target network (Figure 3) with 215 nodes and 1239 edges. The network results manifested that the average degree value of the compounds was 9.68; we received 66 active compounds. Five compounds were considered to have high connection with putative targets of RA: quercetin (), kaempferol (), luteolin (), wogonin (), and isorhamnetin () (Table 2).

3.2. PPI Network and Analysis

Based on the previous results, 87 potential targets were imported in STRING with the organism set as “homo sapiens.” Those data with a confidence level higher than 0.4 constructed a complete PPI network including 87 nodes and 1062 edges. We chose three important parameters, including “degree (DC),” “betweenness centrality (BC),” and “closeness centrality (CC),” as the thresholds to screen the key genes to establish major hub nodes by Cytoscape. First, we set the medium of DC, BC, and CC as thresholds, , , and , respectively, and obtained 37 nodes and 472 edges. Then, the 37 key nodes were further filtered by the second threshold of , , and and we got 18 nodes and 153 edges at last (Figure 4). The hub targets of QHSSD in RA therapy are as follows: IL-6, TP53, VEGFA, JUN, TNF, MAPK1, MAPK8, EGF, PTGS2, MAPK3, EGFR, ESR1, IL-1B, CAT, CCL2, MAPK14, MMP2, and NOS3.

3.3. GO and KEGG Pathway Enrichment Analysis

We performed GO and KEGG pathways to further elucidate the biological functions of 18 hub targets systematically. There were 1534 GO items (1474 terms of biological process, 40 terms of molecular function, and 20 terms of cellular component) and 142 pathways of KEGG enrichment. The top 20 GO items and KEGG pathways are displayed in Figure 5. The biological processes of targets were mainly enriched in oxidative response/stress (GO: 0034599, GO: 0006979, GO: 0000302, and GO: 0034614), inflammatory response (GO: 0032496, GO: 0031663, and GO: 0071222), and biochemical stimulation (GO: 0062197, GO: 0002237, GO: 0071216, and GO: 0071219). Additionally, the KEGG pathways were most enriched in immune- and inflammatory-associated pathways, including the MAPK signaling pathway (hsa04010), IL-17 signaling pathway (hsa04657), TNF signaling pathway (hsa04668), Toll-like receptor signaling pathway (hsa04620), and C-type lectin receptor signaling pathway (hsa04625), and involved in angiogenesis and atherosclerosis pathways: fluid shear stress and atherosclerosis (hsa05418) and relaxin signaling pathway (hsa04926). The compounds of QHSSD may act on these signaling pathways and yield a therapeutic effect.

(a)

(b)

(c)

(d)

3.4. Molecular Docking Analysis

To further elucidate the mechanism of action, 18 hub targets with their related compounds were docked and screened by the binding affinity on AutoDock Vina based on the compound-target network. We investigated a total of 39 pairs of ligand-protein interaction (Table 3), with binding as the screening threshold. The lower the value of the docking affinity, the stronger the binding ability between the compound and the active site of the target. Accordingly, JUN and TNF get a low binding ability with their related compounds, as well as JUN-quercetin, JUN-kaempferol, and TNF-quercetin. Therefore, 3 pairs of combinations were excluded, and at last, we retrieved 47 compound-target pairs. In addition, cognate docking was used to verify the rationality of molecular docking. In order to verify the rationality of the docking for the TP53 and quercetin system, we investigated the docking of IRAK4 (PDB ID: 2a9i) and quercetin; the binding affinity of IRAK4-quercetin was −6.6 kcal/mol. Small-molecule ligand potentially could fit into the interface port formed by the interaction of amino acid residues in protein. Overall, these complexes provide further evidence that these proteins could act as therapeutic target for the compound of QHSSD.

3.5. Molecular Dynamics Simulations

To further research the stability of the molecular dynamics trajectory, we took a complex (TP53-quercetin) as an example to analyze the root mean square deviation (RMSD) of the TP53 backbone atoms in this complex, as shown in Figure 6. The RMSD of the TP53 backbone atom began to show an upward trend, whose trajectory was stable and reached an equilibrium state. The RMSD remained stable at about 0.2 nm after 4 ns, but in 7 ns, the trajectory fluctuated slightly. In addition, after 2 ns of the molecular dynamics simulation, the potential energy of the system also tends to be the minimum and stabilizes. The total interaction energy value was . The 4 ns simulated trajectory was extracted by GROMACS as the representative frame of the interaction between quercetin and TP53 (Figure 7); we could find that the amino acid ARG15 of TP53 and quercetin have hydrophobic interactions, resulting in quercetin and TP53 more closely. Moreover, quercetin and TP53 could form three hydrogen bonds (the hydrogen bond between the oxygen atom of quercetin and the hydrogen atom of amino acid residue GLY204 and the hydrogen bond between the hydrogen atom of quercetin and the nitrogen atom of amino acid residues GLU209 and GLN28) and caused the complex to bind more tightly.

4. Discussion

There has been increased researches on how traditional Chinese herbal medicine or formulae work on various diseases based on network pharmacology [21, 22, 24]. In our study, we first collected the putative targets of compounds in QHSSD from the database. And then, the compound-target network was constructed, which revealed that one herb could correspond with multiple components and multiple targets. The traditional Chinese medicine integrative database (TCMID) and a Bioinformatics Analysis Tool for Molecular Mechanism of Traditional Chinese Medicine (BATMAN-TCM) database were used to predict the targets of Chinese herbal medicine [27, 28], whereas we found that the results of these databases were different because of the screening criteria. To ensure the consistency of data, we only used the TCMSP database to collect active compounds and related targets. In compound-target network analysis, the top five compounds (quercetin, kaempferol, luteolin, wogonin, and isorhamnetin) were considered to be crucial compounds in QHSSD. Surprisingly, these compounds are both flavonoids which have antioxidant and anti-inflammatory effects [39–42]. For instance, researches indicated that quercetin mediate the anti-inflammatory mechanism of reducing acute laryngitis by restoring the Th17/Treg balance and activating heme oxygenase 1 and isorhamnetin increased the levels of TNF-α, IL-1β, IL-6, IL-17A, and IL-17F as well as decreased levels of IL-35 and IL-10 in the CIA mice, resulting in an anti-inflammatory effect [43, 44]. Luteolin prevented cartilage destruction in OA rats in vivo through inhibiting IL-1β-induced phosphorylation of NF-κB [45]. In conclusion, these active compounds are crucial material basis of the potential therapeutic effect of QHSSD on RA.

Further, we selected 18 hub targets from the PPT network, as follows: IL-6, TP53, VEGFA, JUN, TNF, MAPK1, MAPK8, EGF, PTGS2, MAPK3, EGFR, ESR1, IL-1B, CAT, CCL2, MAPK14, MMP2, and NOS3. Combined with GO annotation and pathway enrichment analysis, we preliminarily considered that the treatment of RA by QHSSD may participate in biological processes of proliferation, inflammatory response, and angiogenesis. Previous studies indicated that IL-6 and TNF antagonists decrease the inflammatory response of RA patients and IL-1B stimulates inflammation and degradation of the bone and cartilage [46–49]. EGF is an important growth factor, regulating the proliferation of cells. MAPK enables the growth and proliferation of cells by binding with other protein kinases or factors [50]. In addition, the formation of new blood vessels supplied nutrients and oxygen to the augmented inflammatory cell mass, contributing to perpetuation of joint disease [51]. VEGFA, a proangiogenic molecule, promotes the angiogenic phenotype of RA and is upregulated in RA [52, 53]. These literatures are consistent with our prediction. Therefore, these hub targets are considered as crucial proteins in RA therapy.

To further verify the hypothesis of action mechanism, we performed molecular docking on the hub target binding with crucial compounds. Surprisingly, we found that quercetin, luteolin, and wogonin could both fitly bind with TP53, a transcription factor, regulating cell cycle. Moreover, quercetin and luteolin also act on VEGFA. Previous studies revealed that TP53 mutations are linked to the VEGF pathway and showed that the transcriptional target of TP53 inhibited proliferation and angiogenesis via VEGFA [54–56]. We therefore hypothesized that quercetin and luteolin positively worked on TP53 to regulate cell cycle, and meanwhile act on VEGFA to suppress vascular formation. Additionally, luteolin could simultaneously act on EGFR and TNF, decreasing cell proliferation and the level of proinflammatory cytokines. This is in line with the researches before; the result of which manifested that luteolin inhibited IL-1β-induced expression of NO, PGE2, TNF-α, and MMP-2 [45, 50, 57, 58]. On the one hand, the synergistic action of quercetin and luteolin produces the effect of inhibiting cell proliferation and angiogenesis by multiple proteins and also has an anti-inflammatory effect. On the other hand, kaempferol acting on NOS3 may increase NO levels and prevent oxidative damage [59]. From a holistic aspect, QHSSD could inhibit the cell proliferation and angiogenesis and yield anti-inflammatory and antioxidant effects through regulating TP53, VEGFA, EGFR, TNF, and NOS3.

MD is crucial to the discovery of new drug [60] and is a powerful method for studying protein-ligand interactions [61]. In contrast to experimental methods alone, MD simulations can address diseases caused by protein misfolding and virtual screening; it also identifies the stability of protein-ligand complexes and ligand binding kinetics [62]. Previous studies found that quercetin can play an important role in inflammation by acting on the TP53 target [54–56]. In our study, MD simulation was used to investigate the interactions of TP53-quercetin. This complex was simulated within 10 ns and reached an equilibrium state. We found that ligand has three hydrogen bonds with protein; a study [63] had shown that hydrogen bonding has an advantage over shapeless forces if it can form hydrogen bonds. In this study, MD simulation revealed that the structure of TP53-quercetin was very stable when the simulation temperature was 300 K. It is important to explore the application of new drug by molecular docking and MD simulations [61]; however, it is still necessary to perform in vitro experiment to verify the prediction.

5. Conclusion

To conclude, we could learn that multicompounds act on multitargets in TCM therapy. We speculated that the therapeutic effect on the treatment of RA by QHSSD may be focused on antiangiogenesis, meanwhile, anti-inflammation, especially TP53 and VEGFA targets, plays a vital role during the treatment. This study was based on the combination of network pharmacology and molecular docking to predict and identify the mechanism of QHSSD in the treatment of RA, which provides the foundation for subsequent experimental investigation and offers ideas for the multidimensional and multilevel research of TCM formulae.

Data Availability

The data used to support the findings of this study are included within the supplementary information files.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this article.

Authors’ Contributions

Xiaoli Bi designed this study and provided technical support toward the study. Zhihao Zeng, Guanlin Xiao, Ruipei Yang, Huajing Huang, and Huixian Zhong collected the data. Zhihao Zeng and Jiaoting Hu analyzed the data and wrote the manuscript. Jieyi Jiang, Sumei Li, and Yangxue Li revised the manuscript.

Acknowledgments

Thanks are due for Warren L. DeLano, the original PyMOL author of the open-source software. Our research was supported by the “Significant News Drugs Creation” Science and Technology Major Project of China (grant number 2018ZX09721005).

Supplementary Materials

Table S1: information of active compounds of QHSSD. Figure S1: herb-compound-target network of QHSSD. (Supplementary Materials). (Supplementary Materials)

References

T.-p. Li, A.-h. Zhang, J.-h. Miao et al., “Applications and potential mechanisms of herbal medicines for rheumatoid arthritis treatment: a systematic review,” RSC Advances, vol. 9, no. 45, pp. 26381–26392, 2019.
View at: Publisher Site | Google Scholar
T. D. Wilsdon and C. L. Hill, “Managing the drug treatment of rheumatoid arthritis,” Australian Prescriber, vol. 40, no. 2, pp. 51–58, 2017.
View at: Publisher Site | Google Scholar
F. Angelotti, A. Parma, G. Cafaro, R. Capecchi, A. Alunno, and I. Puxeddu, “One year in review 2017: pathogenesis of rheumatoid arthritis,” Clinical and Experimental Rheumatology, vol. 35, 2017.
View at: Google Scholar
A. Gibofsky, “Overview of epidemiology, pathophysiology, and diagnosis of rheumatoid arthritis,” American Journal of Managed Care, vol. 18, pp. 295–302, 2012.
View at: Google Scholar
M. V. Nemtsova, D. V. Zaletaev, I. V. Bure et al., “Epigenetic Changes in the Pathogenesis of Rheumatoid Arthritis,” Frontiers in Genetics, vol. 10, 2019.
View at: Publisher Site | Google Scholar
J. M. Kahlenberg and D. A. Fox, “Advances in the medical treatment of rheumatoid arthritis,” Hand Clinics, vol. 27, no. 1, pp. 11–20, 2011.
View at: Publisher Site | Google Scholar
P. Conigliaro, P. Triggianese, E. D. Martino, M. S. Chimenti, and R. Perricone, “Challenges in the treatment of rheumatoid arthritis,” Autoimmunity Reviews, vol. 18, no. 7, pp. 706–713, 2019.
View at: Publisher Site | Google Scholar
H.-D. Pan, Y. Xiao, W.-Y. Wang, R.-T. Ren, E. L.-H. Leung, and L. Liu, “Traditional Chinese Medicine as a Treatment for Rheumatoid Arthritis: From Empirical Practice to Evidence-Based Therapy,” Engineering, vol. 5, no. 5, pp. 895–906, 2019.
View at: Publisher Site | Google Scholar
Q. Xing, L. Fu, Z. Yu, X. Zhou, and A. Vannacci, “Efficacy and Safety of Integrated Traditional Chinese Medicine and Western Medicine on the Treatment of Rheumatoid Arthritis: A Meta-Analysis,” Evidence-Based Complementary and Alternative Medicine, vol. 2020, Article ID 4348709, 15 pages, 2020.
View at: Publisher Site | Google Scholar
H. Chen, T. Pan, P. Liu, P. Wang, S. Xu, and H. Yamada, “Baihu Jia Guizhi Decoction Improves Rheumatoid Arthritis Inflammation by Regulating Succinate/SUCNR1 Metabolic Signaling Pathway,” Evidence-Based Complementary and Alternative Medicine, vol. 2019, Article ID 3258572, 11 pages, 2019.
View at: Publisher Site | Google Scholar
T.-T. Sham, C.-O. Chan, Y.-H. Wang, J.-M. Yang, D. K.-W. Mok, and S.-W. Chan, “A review on the traditional Chinese medicinal herbs and formulae with hypolipidemic effect,” BioMed Research International, vol. 2014, Article ID 925302, 21 pages, 2014.
View at: Publisher Site | Google Scholar
J.-J. Wang, “Treatment of arthralgia by Qianghuo Shengshi decoction and Huangqi Guizhi Wuwu decoction (report of 300 cases),” Journal of Hubei University of Science and Technology(Medical Sciences), vol. 28, p. 150, 2014.
View at: Google Scholar
M.-Z. Li, “Clinical Observation on Treatment of Rheumatoid Arthritis with Modification of Qianghuo Shengshi Decoction Combined with Traditional Chinese Medicine,” Mod. Health Care, pp. 187-188, 2017.
View at: Google Scholar
Y. Li, H. Liu, X. Deng, G. Ding, and S. Lin, “The clinical application status of Qianghuo Shengshi decoction,” Journal of Jiangxi College of Traditional Chinese Medicine, vol. 29, pp. 113-114, 2017.
View at: Google Scholar
F. Hongxia, W. Chunhao, P. Lijuan et al., “Study on the Anti-Inflammatory and Analgesic Effects of Qianghuo Shengshi Decoction,” in The 9th Annual Meeting of Pharmaceutical Toxicology-New Era·New Technology·New Strategy·New HealthWuhan, p. 1, Hubei, China, 2019.
View at: Google Scholar
W. Zhang, Network pharmacology: a further description, Network Pharmacology, 2016.
J. Zhao, J. Yang, S. Tian, and W. Zhang, “A survey of web resources and tools for the study of TCM network pharmacology,” Frontiers of Electrical and Electronic Engineering in China, vol. 7, pp. 17–29, 2019.
View at: Google Scholar
R. Zhang, X. Zhu, H. Bai, and K. Ning, “Network pharmacology databases for traditional Chinese medicine: review and assessment,” Frontiers in Pharmacology, vol. 10, 2019.
View at: Google Scholar
Y. Zhang, X. Li, C. Guo, J. Dong, and L. Liao, “Mechanisms of Spica Prunellae against thyroid-associated ophthalmopathy based on network pharmacology and molecular docking,” BMC Complementary Medicine and Therapies, vol. 20, no. 1, p. 229, 2020.
View at: Publisher Site | Google Scholar
J. Zhang, X. Liu, W. Zhou et al., “A bioinformatics investigation into molecular mechanism of Yinzhihuang granules for treating hepatitis B by network pharmacology and molecular docking verification,” Scientific Reports, vol. 10, 2020.
View at: Google Scholar
H.-M. An, D.-R. Huang, H. Yang et al., “Comprehensive chemical profiling of Jia-Wei-Qi-Fu-Yin and its network pharmacology-based analysis on Alzheimer's disease,” Journal of Pharmaceutical and Biomedical Analysis, vol. 189, article 113467, 2020.
View at: Publisher Site | Google Scholar
D. Gao, S.-N. Wu, C.-E. Zhang et al., “Exploration in the mechanism of rhubarb for the treatment of hyperviscosity syndrome based on network pharmacology,” Journal of Ethnopharmacology, vol. 261, article 113078, 2020.
View at: Publisher Site | Google Scholar
W. Tao, X. Xu, X. Wang et al., “Network pharmacology-based prediction of the active ingredients and potential targets of Chinese herbal _Radix Curcumae_ formula for application to cardiovascular disease,” Journal of Ethnopharmacology, vol. 145, no. 1, pp. 1–10, 2013.
View at: Publisher Site | Google Scholar
Y. Zhou, R. Wu, F.-F. Cai et al., “Xiaoyaosan decoction alleviated rat liver fibrosis via the TGFβ/Smad and Akt/FoxO3 signaling pathways based on network pharmacology analysis,” Journal of Ethnopharmacology, vol. 264, article 113021, 2021.
View at: Publisher Site | Google Scholar
J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, no. 1, p. 13, 2014.
View at: Publisher Site | Google Scholar
X. Xu, W. Zhang, C. Huang et al., “A novel chemometric method for the prediction of human oral bioavailability,” International Journal of Molecular Ences, vol. 13, no. 6, pp. 6964–6982, 2012.
View at: Publisher Site | Google Scholar
S. Tian, J. Wang, Y. Li, D. Li, L. Xu, and T. Hou, “The application of _in silico_ drug-likeness predictions in pharmaceutical research,” Advanced Drug Delivery Reviews, vol. 86, pp. 2–10, 2015.
View at: Publisher Site | Google Scholar
C. J. Grondin, A. P. Davis, T. C. Wiegers, J. A. Wiegers, and C. J. Mattingly, “Accessing an expanded exposure science module at the comparative toxicogenomics database,” Environmental Health Perspectives, vol. 126, no. 1, p. 014501, 2018.
View at: Publisher Site | Google Scholar
M. Rebhan, V. Chalifa-Caspi, J. Prilusky, and D. Lancet, “GeneCards: a novel functional genomics compendium with automated data mining and query reformulation support,” Bioinformatics, vol. 14, no. 8, pp. 656–664, 1998.
View at: Publisher Site | Google Scholar
J. C. O. Venny, “An interactive tool for comparing lists with Venn’s diagrams,” 2007, https://bioinfogp.cnb.csic.es/tools/venny/index.html.
View at: Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., Cytoscape: A software environment for integrated models of biomolecular interaction networks, Routledge, 2003.
S. Damian, A. L. Gable, L. David et al., “STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, pp. D607–D613, 2019.
View at: Publisher Site | Google Scholar
J. B. Zou, H. B. Chai, X. F. Zhang, D. Y. Guo, and Y. J. Shi, “Reconstruction of the lncRNA-miRNA-mRNA network based on competitive endogenous RNA reveal functional lncRNAs in cerebral infarction,” Scientific Reports, vol. 9, no. 1, p. 12176, 2019.
View at: Publisher Site | Google Scholar
Y. Guangchuang, W. Li-Gen, Y. Guang-Rong, and H. Qing-Yu, “DOSE: an R/Bioconductor package for disease ontology semantic and enrichment analysis,” Bioinformatics, vol. 31, pp. 608-609, 2015.
View at: Publisher Site | Google Scholar
O. Trott and A. J. Olson, “AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading,” Journal of Computational Chemistry, vol. 31, 2009.
View at: Publisher Site | Google Scholar
T. Marco, V. Alessia, and A. A. Maria, “Reverse screening on Indicaxanthin from Opuntia ficus-indica as natural chemoactive and chemopreventive agent,” Journal of Theoretical Biology, vol. 455, pp. 147–160, 2018.
View at: Google Scholar
C. Kutzner, S. Páll, M. Fechner, A. Esztermann, B. L. De Groot, and H. Grubmüller, “More bang for your buck: improved use of GPU nodes for GROMACS 2018,” Journal of Computational Chemistry, vol. 40, no. 27, pp. 2418–2431, 2019.
View at: Publisher Site | Google Scholar
J. Huang, S. Rauscher, G. Nawrocki et al., “Mackerell, CHARMM36m: an improved force field for folded and intrinsically disordered proteins,” Biophysical Journal, vol. 112, pp. 175a–176a, 2017.
View at: Google Scholar
S. Martinez-Florez, J. González-Gallego, J. M. Culebras, and M. J. Tuñón, “Flavonoids: properties and anti-oxidizing action,” Nutricion Hospitalaria, vol. 17, 2002.
View at: Google Scholar
L. Zhang, A. S. Ravipati, S. R. Koyyalamudi et al., “Antioxidant and anti-inflammatory activities of selected medicinal plants containing phenolic and flavonoid compounds,” Journal of Agricultural and Food Chemistry, vol. 59, no. 23, pp. 12361–12367, 2011.
View at: Publisher Site | Google Scholar
M. H. Johnson, E. G. de Mejia, J. Fan, M. A. Lila, and G. G. Yousef, “Anthocyanins and proanthocyanidins from blueberry–blackberry fermented beverages inhibit markers of inflammation in macrophages and carbohydrate-utilizing enzymes in vitro,” Molecular Nutrition & Food Research, vol. 57, no. 7, pp. 1182–1197, 2013.
View at: Publisher Site | Google Scholar
Y. S. Yi, “Regulatory roles of flavonoids on inflammasome activation during inflammatory responses,” Molecular Nutrition & Food Research, vol. 62, article e1800147, 2018.
View at: Google Scholar
L. Li, S. Lin, D. Zhang, H. Xie, and Z. Wei, “Quercetin restores Th17/Treg balance and activates heme oxygenase 1 mediated anti-inflammatory effect on the mechanism of reducing acute laryngitis,” Hebei Medical Journal, vol. 42, pp. 2730–2734, 2020.
View at: Google Scholar
Y. Huang, L. Guo, R. Chitti et al., “Wogonin ameliorate complete Freund's adjuvant induced rheumatoid arthritis via targeting NF-κB/MAPK signaling pathway,” BioFactors (Oxford, England), vol. 46, no. 2, pp. 283–291, 2020.
View at: Publisher Site | Google Scholar
J. Fei, B. Liang, C. Jiang, H. Ni, and L. Wang, “Luteolin inhibits IL-1β-induced inflammation in rat chondrocytes and attenuates osteoarthritis progression in a rat model,” Biomedicine & Pharmacotherapy = Biomedecine & Pharmacotherapie, vol. 109, pp. 1586–1592, 2019.
View at: Publisher Site | Google Scholar
P. P. Tak and J. R. Kalden, “Advances in rheumatology: new targeted therapeutics,” Arthritis Research & Therapy, vol. 13, Supplement 1, p. S5, 2011.
View at: Publisher Site | Google Scholar
J. S. Smolen, R. B. M. Landewé, J. W. J. Bijlsma et al., “EULAR recommendations for the management of rheumatoid arthritis with synthetic and biological disease-modifying antirheumatic drugs: 2019 update,” Annals of the Rheumatic Diseases, vol. 79, no. 6, pp. 685–699, 2020.
View at: Publisher Site | Google Scholar
I. B. McInnes and G. Schett, “The pathogenesis of rheumatoid arthritis,” The New England Journal of Medicine, vol. 365, no. 23, pp. 2205–2219, 2011.
View at: Publisher Site | Google Scholar
S. Mateen, A. Zafar, S. Moin, A. Q. Khan, and S. Zubair, “Understanding the role of cytokines in the pathogenesis of rheumatoid arthritis,” Clinica Chimica Acta; International Journal of Clinical Chemistry, vol. 455, pp. 161–171, 2016.
View at: Publisher Site | Google Scholar
M. F. Rimawi, P. B. Shetty, H. L. Weiss et al., “Epidermal growth factor receptor expression in breast cancer association with biologic phenotype and clinical outcomes,” Cancer, vol. 116, no. 5, pp. 1234–1242, 2010.
View at: Publisher Site | Google Scholar
P. C. Taylor, “VEGF and imaging of vessels in rheumatoid arthritis,” Arthritis Research, vol. 4, Supplement 3, 2002.
View at: Google Scholar
A. Paradowska-Gorycka, A. Pawlik, K. Romanowska-Prochnicka et al., “Relationship between VEGF gene polymorphisms and serum VEGF protein levels in patients with rheumatoid arthritis,” PLoS One, vol. 11, no. 8, article e0160769, 2016.
View at: Publisher Site | Google Scholar
J. P. Yi, Y. Z. Wu, N. Yu, Z. W. Yu, F. Y. Xie, and Q. Yuan, “VEGF gene polymorphisms affect serum protein levels and alter disease activity and synovial lesions in rheumatoid arthritis,” Medical Science Monitor : International Medical Journal of Experimental and Clinical Research, vol. 22, pp. 316–324, 2016.
View at: Publisher Site | Google Scholar
A. M. Li, A. Boichard, and R. Kurzrock, “MutatedTP53is a marker of increasedVEGFexpression: analysis of 7,525 pan-cancer tissues,” Cancer Biology & Therapy, vol. 21, no. 1, pp. 95–100, 2020.
View at: Publisher Site | Google Scholar
X. Chen, K. Zeng, M. Xu et al., “P53-induced miR-1249 inhibits tumor growth, metastasis, and angiogenesis by targeting VEGFA and HMGA2,” Cell Death & Disease, vol. 10, no. 2, p. 131, 2019.
View at: Publisher Site | Google Scholar
J. J. Wheler, F. Janku, A. Naing et al., “TP53 alterations correlate with response to VEGF/VEGFR inhibitors: implications for targeted therapeutics,” Molecular Cancer Therapeutics, vol. 15, no. 10, pp. 2475–2485, 2016.
View at: Publisher Site | Google Scholar
B. C. Zhang, Z. Li, W. Xu, C. H. Xiang, and Y. F. Ma, “Luteolin alleviates NLRP3 inflammasome activation and directs macrophage polarization in lipopolysaccharide-stimulated RAW264.7 cells,” American Journal of Translational Research, vol. 10, no. 1, pp. 265–273, 2018.
View at: Google Scholar
D. M. Anson, R. M. Wilcox, E. D. Huseman et al., “Luteolin decreases epidermal growth factor receptor-mediated cell proliferation and induces apoptosis in glioblastoma cell lines,” Basic & Clinical Pharmacology & Toxicology, vol. 123, no. 6, pp. 678–686, 2018.
View at: Publisher Site | Google Scholar
A. A. Olabiyi, F. B. Carvalho, N. B. Bottari et al., “Tiger nut and walnut extracts modulate extracellular metabolism of ATP and adenosine through the NOS/cGMP/PKG signalling pathway in kidney slices,” Phytomedicine : International Journal of Phytotherapy and Phytopharmacology, vol. 43, pp. 140–149, 2018.
View at: Publisher Site | Google Scholar
S. A. Hollingsworth and R. O. Dror, “Molecular dynamics simulation for all,” Neuron, vol. 99, no. 6, pp. 1129–1143, 2018.
View at: Publisher Site | Google Scholar
C. J. Morris and D. D. Corte, “Using molecular docking and molecular dynamics to investigate protein-ligand interactions,” Modern Physics Letters B, vol. 35, no. 8, p. 2130002, 2021.
View at: Publisher Site | Google Scholar
X. Liu, D. Shi, S. Zhou, H. Liu, H. Liu, and X. Yao, “Molecular dynamics simulations and novel drug discovery,” Expert Opin Drug Discov, vol. 13, no. 1, pp. 23–37, 2018.
View at: Publisher Site | Google Scholar
Z. Cao, Z. Song, F. Liang, X. An, and Y. Yang, “Hydrogen bonding sewing interface,” RSC Advances, vol. 10, no. 30, pp. 17438–17443, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Zhihao Zeng 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

887

Downloads

1256

Citations