Research Article | Open Access
Pan Kang, Zhiming Wu, Yue Zhong, Zihao Wang, Chi Zhou, Shaochuan Huo, Hai Guo, Songtao Li, Kun Xu, Lingyun Liu, Shuai Chen, Hongyu Tang, Haibin Wang, "A Network Pharmacology and Molecular Docking Strategy to Explore Potential Targets and Mechanisms Underlying the Effect of Curcumin on Osteonecrosis of the Femoral Head in Systemic Lupus Erythematosus", BioMed Research International, vol. 2021, Article ID 5538643, 14 pages, 2021. https://doi.org/10.1155/2021/5538643
A Network Pharmacology and Molecular Docking Strategy to Explore Potential Targets and Mechanisms Underlying the Effect of Curcumin on Osteonecrosis of the Femoral Head in Systemic Lupus Erythematosus
Background. Systemic lupus erythematosus (SLE) is a refractory immune disease, which is often complicated with osteonecrosis of the femoral head (ONFH). Curcumin, the most active ingredient of Curcuma longa with a variety of biological activities, has wide effects on the body system. The study is aimed at exploring the potential therapeutic targets underlying the effect of curcumin on SLE-ONFH by utilizing a network pharmacology approach and molecular docking strategy. Methods. Curcumin and its drug targets were identified using network analysis. First, the Swiss target prediction, GeneCards, and OMIM databases were mined for information relevant to the prediction of curcumin targets and SLE-ONFH-related targets. Second, the curcumin target gene, SLE-ONFH shared gene, and curcumin-SLE-ONFH target gene networks were created in Cytoscape software followed by collecting the candidate targets of each component by R software. Third, the targets and enriched pathways were examined by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Eventually, a gene-pathway network was constructed and visualized by Cytoscape software; key potential central targets were verified and checked by molecular docking and literature review. Results. 201 potential targets of curcumin and 170 related targets involved in SLE-ONFH were subjected to network analysis, and the 36 intersection targets indicated the potential targets of curcumin for the treatment of SLE-ONFH. Additionally, for getting more comprehensive and accurate candidate genes, the 36 potential targets were determined to be analyzed by network topology and 285 candidate genes were obtained finally. The top 20 biological processes, cellular components, and molecular functions were identified, when corrected by a value ≤ 0.05. 20 related signaling pathways were identified by KEGG analysis, when corrected according to a Bonferroni value ≤ 0.05. Molecular docking showed that the top three genes (TP53, IL6, VEGFA) have good binding force with curcumin; combined with literature review, some other genes such as TNF, CCND1, CASP3, and MMP9 were also identified. Conclusion. The present study explored the potential targets and signaling pathways of curcumin against SLE-ONFH, which could provide a better understanding of its effects in terms of regulating cell cycle, angiogenesis, immunosuppression, inflammation, and bone destruction.
Systemic lupus erythematosus (SLE) is a chronic autoimmune disease, which prevalence ranges from 20 to 150 cases per 100,000 population and appears to be increasing as the disease is recognized more readily and survival increases . SLE can cause multiple organ dysfunction such as the lung, kidney, blood vessel, and joint; the diverse clinical manifestations of SLE present a challenge to the clinician . Osteonecrosis (ON) remains a serious complication in SLE with prevalence ranging from 10% to 50%, and the femoral head is the most common site [3–9]. The pathogenesis of SLE-ONFH has not yet been fully elucidated. Except for joint arthroplasty in the late inactive phase, the active phase usually requires effective drugs or other measures for SLE itself, which is more complicated. To date, the drug for SLE is mainly glucocorticoids with immunosuppressive effect, which have achieved good therapeutic effects. However, some side effects have resulted because of their extensive use. A large number of studies have reported that the use of glucocorticoids can cause femoral head necrosis (ONFH) [10–12] and long-term high-dose cyclophosphamide can induce infection, reproductive toxicity, and other adverse reactions [13, 14]. Obviously, none of these seems to be an ideal treatment for SLE-ONFH in the future. Traditional Chinese medicine (TCM) has the characteristics of multitarget, multicomponent, multifunction, and generally low toxicity or side effects [15–18]. It has been widely used in the treatment of various diseases including immune diseases [19, 20] and has attracted more and more attention in recent years. The Chinese medicine arsenic being found to be effective in treating leukemia is one of the best examples . Curcumin, an effective monomeric component mainly extracted from the traditional Chinese medicine turmeric, has been reported to have protective effects on SLE and ONFH [22–25]. However, its potential targets and mechanisms for the treatment of SLE-ONFH are not still very clear, and there is no systematic and comprehensive understanding of the relationship between therapeutic targets and pathways. With the improvement of bioinformatics and related databases, network pharmacology and molecular docking strategies are widely used in drug research, which revealed the interrelationship between drug-target-pathway-disease. In this study, we constructed an interaction network between SLE and ONFH, as well as their interaction network with curcumin’s targets. And molecular docking technology was used to verify the potential targets of curcumin in the treatment of SLE-ONFH. In addition, relevant signal pathways are obtained by bioinformatics analysis to reveal the underlying mechanisms. The entire research workflow is shown in Figure 1.
2. Materials and Methods
2.1. Data Preparation
2.1.1. Chemical Structures
The 2D chemical structure of curcumin (Figure 2) was obtained from PubChem (http://pubchem.ncbi.nlm.nih.gov), which is an open chemical database that provides information on compound structures and descriptive data.
2.1.2. Prediction of Curcumin Targets
The Swiss target prediction database (http://www.swisstargetprediction.ch/) is a web-based tool, on-line since 2014, to perform ligand-based target prediction for any bioactive small molecule . To get curcumin relevant target genes, the “SMILES” format of curcumin was obtained and inputted into the Swiss target prediction database by searching curcumin in the PubChem database. In addition, the HERB database (http://herb.ac.cn/) was used to get curcumin targets for a comprehensive collection of drug targets, which is a high-throughput experiment- and reference-guided database of traditional Chinese medicine . The targets are mainly from TCMID (http://www.megabionet.org/tcmid/) and TCMSP (https://tcmspw.com/tcmsp.php) databases [28, 29]. Finally, each of the predicted targets was inputted into the UniProt database (https://www.uniprot.org/) for screening. In order to increase credibility, we only need human targets which have been annotated and reviewed, and the potential targets of curcumin were identified by eliminating duplicate and nonstandard targets.
2.1.3. Predict Targets of Curcumin against SLE-ONFH
The SLE and ONFH-related targets were identified in the GeneCards database (https://www.genecards.org/) using the phrase “Systemic lupus erythematosus”, “SLE”, “Osteonecrosis of Femoral Head”, and “ONFH” as the keyword. The targets were also supplemented with the OMIM database (https://omim.org/) . A Venn script was installed in the R software 4.0.2 software to get targets of SLE-ONFH, which was the intersection genes of two diseases. To determine the predicted targets of curcumin relevant to SLE-ONFH, the potential targets of curcumin were taken intersection with the relevant targets involved in SLE-ONFH by R software.
2.2. Network Construction
2.2.1. Curcumin Target Network and SLE-ONFH Target Network
The curcumin target network and SLE-ONFH network were constructed and visualized using Cytoscape 3.8.0 software [31–33]; they were uploaded to the String database (https://string-db.org/), which is including direct and indirect interactions between proteins and calculates a confidence score for all protein interactions [34, 35]. The score is proportional to the confidence in the protein interaction, ranged from 0 to 1, and we set the score greater than 0.9 to be high confidence [36, 37].
2.2.2. Curcumin-SLE-ONFH Target PPI Network
Shared targets of curcumin, SLE, and ONFH were obtained by running the Venn script as described above, which was installed in the R 4.0.2 software. They were uploaded to Cytoscape software. The PPI data were obtained from the Database of Interacting Proteins (DIP™), Biological General Repository for Interaction Datasets (Bio GRID), Human Protein Reference Database (HPRD), IntAct Molecular Interaction Database (IntAct), Molecular INTeraction database (MINT), and biomolecular interaction network database (BIND) using the plugin BisoGenet of Cytoscape software .
2.3. Topology Analysis
After uploading the common gene above, the nodes with topological importance in the interaction network were screened by calculating Degree Centrality (DC), betweenness centrality (BC), closeness centrality (CC), Eigenvector Centrality (EC), local average connectivity-based method (LAC), and Network Centrality (NC) with the Cytoscape plugin CytoNCA. At the heart of Cytoscape is a network. A simple network diagram consists of nodes and edges. Each node represents a gene. The node-node connection represents the interaction between these nodes. DC believes that the greater the number of neighbors of a node, the greater its impact. CC was used to calculate the importance of nodes. BC is a kind of centrality measure. EC is an indicator in ranking the importance of nodes in the network. The closer the EC of a node is to the network radius, the closer the node is to the network center. LAC is often used to determine the importance of proteins by assessing the relationship between proteins and their neighbors. NC is a method based on the edge clustering coefficient, which considers not only the centrality of nodes but also the relationship between nodes and adjacent nodes [39–41]. These parameters represent the topological importance, and they have been reported about their definitions and computational formulas. They were all used in network pharmacology and systems pharmacology .
2.4. Enrichment Analysis
GO analysis with the biological process, cellular component, and molecular function was carried out by R software. The result is obtained by running the script installed on R software including “colorspace,” “stringi,” “ggplot2,” “BiocManager,” “clusterProfiler,” and “enrichplot.” The top 20 items were selected and visualized; a Bonferroni-corrected value ≤ 0.05 was set as the default option. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis method is similar to the GO analysis above.
2.5. Gene-Pathway Network Analysis
The gene-pathway network was constructed based on the significantly enriched pathways with genes that regulated these pathways. The topological analysis of 20 pathways and 285 genes was carried out with BC. The squares represented target genes, and the V-shapes represented pathways in the network.
2.6. Molecular Docking
Molecular docking simulations were used to explore the potential interaction between the top 3 genes of the network above. The SDF format of curcumin was downloaded from the PubChem database; Chem 3D software was used to convert the SDF format into mol2 format file; then, the RESB database was used to get the PDB format structure of the top 3 target proteins. After that, solvent molecules and ligands were removed by Pymol software and saved as pdbqt format. Docking simulations were performed in Autodock1.1.2 software, and Discovery Studio 2020 was used to visually analyze the docking conformation at last [43, 44].
3.1. Target Network Analysis
3.1.1. Curcumin Target Network
100, 152, and 11 target genes of curcumin were identified by searching the Swiss, HERB, and STITCH databases, respectively, and 240 genes were selected as curcumin target candidate genes through running R software to remove duplicate genes. In order to ensure credibility, the 240 target genes were inputted into the STRING database; by setting the highest confidence (0.9000) and selecting the condition “Homo sapiens,” 201 targets were finally identified after hiding disconnected nodes in the network. The PPI network of them was saved as “tsv” file and visualized by Cytoscape software according to the degree (Figure 3). This network represented the regulatory relationship of the curcumin targets. In the figure, the nodes represented the target points, and the edges indicated the relationships between the targets. Four important values (Table 1) average shortest path length (ASPL), betweenness centrality (BC), closeness centrality (CC), and clustering coefficient were analyzed; ASPL represented the smallest number of links that a node needed to connect the whole network. In a network, the smaller the ASPL, the more efficient the signal conduction through a node. BC represented the probability that the signal passed through the nodes, and the higher the value of BC, the more important the node. The closeness centrality (CC) and the clustering coefficient could reflect the network tightness, and a tighter network was of higher efficiency [45, 46]. The degree of a node referred to the number of edges between nodes in a network. The greater a node’s degree is, the more nodes are connected to it and the more important it is in the network [47, 48]. We used the value of the degree to analyze the node’s importance in the network and screening the top 30 genes, which is more than twice the median of the target genes at least . It indicated that these targets might be more effective in the treatment of SLE-ONFH.
3.1.2. Targets and Network of SLE-ONFH
After removing the duplication, we obtained 4338 SLE-related targets and 300 ONFH-related targets from the GeneCards and OMIM databases, respectively. In order to improve the comprehensiveness and credibility of the data, these targets were all included in the current study. Finally, 170 intersection genes related to SLE-ONFH were obtained by running the Venn map script installed in R software, and 36 genes related to SLE, ONFH, and curcumin were also received, as shown in Figure 4. The STRING database and Cytoscape software were used to construct and visualize the PPI network of 170 genes (Figure 5); the method is similar to the description in 3.1.1.
3.1.3. Curcumin-SLE-ONFH Target Network
According to the network pharmacology analysis method, 36 overlapping genes for curcumin, SLE, and ONFH were defined as potential genes (Table 2) for curcumin against SLE-ONFH. They were imported into the STRING database, and the results were saved as a “tsv” file. Then, the file was imported into Cytoscape software, and the curcumin-SLE-ONFH target network diagram was generated (Figure 6). We found that the number of nodes was 36, the number of edges was 406, the average node degree was 22.55, and the genes with a degree greater than 22.55 were represented by an ellipse in the network; they may be important, among which TP53, VEGFA, IL6, TNF, EGFR, CASP3, MMP9, etc.
3.2. Identification of Candidate Targets for Curcumin against SLE-ONFH Based on Topology Analysis
To reveal the mechanisms underlying curcumin’s effects on SLE-ONFH and obtain candidate genes more comprehensively and accurately, topological analysis of 36 potential genes was performed through the plugin BisoGenet and CytoNCA of Cytoscape; a network consisting of 3653 nodes and 87728 edges is presented in Figure 7(a). The median degree of all nodes was 28, and the nodes with more than 56 degrees were identified as significant targets according to the previous research . A network of significant targets for curcumin against SLE-ONFH was constructed, and it contained 945 nodes and 39,757 edges (Figure 7(b)). The median values of DC, BC, CC, EC, LAC, and NC were 117, 6,403.125, 0.318309, 0.024411, 19.61376, and 27.984, respectively. The candidate targets were further screened, and 285 targets with were identified (Figure 7(c)).
3.3. GO and KEGG Enrichment Analysis
3.3.1. GO Analysis
GO analysis of 285 candidate targets was analyzed based on biological process (BP), cellular component (CC), and molecular function (MF). 2,364 GO terms were finally enriched, 1,973 in biological process, 200 in cellular component, and 191 in molecular function. The top 20 terms are shown in Figure 8 ( value < 0.05). The results demonstrated many targets involved in multiple BPs associated with mRNA catabolic process, regulation of DNA-binding transcription factor activity, and regulation of cell cycle phase transition, which confirmed the correlation with the gene transcriptional activity in SLE-ONFH. The CC results suggested that most targets are localized to the focal adhesion, cell-substrate junction, and nuclear chromatin. The MF results showed that many targets are associated with ubiquitin protein ligase binding and DNA-binding transcription factor binding.
3.3.2. KEGG Pathway Analysis
The pathways that are significantly influenced by curcumin in the process of treating SLE-ONFH were identified by KEGG pathway analysis; the analysis method is similar to the GO analysis by running the corresponding script in R software. The top 20 significant pathways ( value < 0.05) are identified and shown in Figure 9(a), including cell cycle, viral carcinogenesis, hepatitis B, hepatitis C, chronic myeloid leukemia, and human T-cell leukemia virus 1 infection pathways; like the GO analysis above, they were closely related to cell proliferation, differentiation, senescence, and apoptosis. The most significant signaling pathway is the cell cycle pathway, in which the central genes such as TP53 play the pivotal role as shown in Figure 9(b). The data of analysis are shown in Table 3.
3.4. Gene-Pathway Network Analysis
The gene-pathway network was constructed based on the significantly enriched pathways and genes that regulated these pathways, which are presented in Figure 10. Topological analysis of 20 pathways and 285 candidate genes was carried out with BC. The squares represented target genes, and the V shapes represented pathways in the network. The network diagram suggested that TP53 has the most maximum BC and was the core target gene. Other several genes also had larger BC, such as HSP90AA1, GAPDH, MYC, and RPS27A. They might be the key target genes for the regulatory pathways of curcumin against SLE-ONFH.
3.5. Molecular Docking
Molecular docking was applied to validate the binding action mode of the top three genes (TP53, VEGFA, and IL6), according to the network analysis. The results revealed that they can interact with curcumin; the combination of curcumin with them is shown in Figure 11. For curcumin and TP53, a ring forms a hydrophobic interaction with ILe232, a hydrogen bond and a hydrophobic interaction with HIS233, and a hydrogen bond with GLY199; the carbonyl group between the a and b rings forms a hydrogen bond with GLU224 and PRO223 (Figure 11(a)). For curcumin and VEGFA, a ring forms a hydrogen bond with GLU93, and the carbonyl group between the a and b rings forms a hydrogen bond with SER95; the b ring forms a hydrophobic interaction with MET94 and a hydrogen bond with ARG82 (Figure 11(b)). For curcumin and IL6, a ring forms a hydrophobic interaction with LEU178 and ARG182 and forms a hydrogen bond with ASP26, the carbonyl group between the a and b rings forms a hydrogen bond with ARG30, and the b ring forms a hydrogen bond with ARG179 (Figure 11(c)). The Affinity (kcal/mol) of curcumin and TP53, VEGFA, and IL6 was -6.5, -5, and -5.9, respectively, which means they can interact well indeed.
4.1. Role of Network Pharmacology and Molecular Docking Strategy in the Research of Autoimmune Diseases
SLE is an autoimmune disease with multiple organ involvement and multiple antibodies. It is an urgent problem in the field of medicine. High disease activity of SLE is closely related to osteonecrosis especially ONFH, which has become the main cause of disability in SLE patients. In the whole process of disease development, immunity, inflammation, apoptosis, and angiogenesis are important factors [50–54]. Because of the complexity of its treatment, it is more worthwhile to find the appropriate drugs and explore the mechanism(s) than the final operation. Traditional Chinese medicine (TCM) has the characteristics of multitargets and multipathway effect, which may be more effective for complex diseases. As a drug research strategy, network pharmacology has been widely applied in drug research, which gives us a better understanding of the role of drug and their compounds . In addition, molecular docking can help reveal the relationship between the drug and the target. Hence, in the present study, combined with topology analysis and molecular docking strategy, we used network pharmacology to investigate the potential targets involved in curcumin against SLE-ONFH by constructing and analyzing target networks, performing enrichment analysis with targets, and revealing the potential mechanism of curcumin against SLE-ONFH.
4.2. Analysis of Central Genes and Significant Pathways of Curcumin against SLE-ONFH
Curcumin has been found to be useful in the treatment of SLE and ONFH, but the central genes and significant pathways of curcumin against SLE-ONFH are not clear. Through the target interaction network, we found 36 potential targets as described in Results; there were 16 targets (TP53, VEGFA, IL6, TNF, EGFR, CASP3, ESR1, MMP9, CCND1, TGFB1, IL1B, CXCL8, FGF2, SRC, MMP2, and SERPINE1) with connectivity greater than the median degree. The TP53 gene is at the core of the whole network, which means that it is a central target. p53 protein is a tumor suppressor that inhibits the growth of aberrant cells; functional p53 is believed to sense DNA damage and, subsequently, to induce DNA repair, growth arrest, or apoptosis of the aberrant cell [56, 57]. On the one hand, as a cell cycle protein, p53 inhibits the proliferation of bone marrow-derived mesenchymal stem cells (BM-MSCs) of SLE patients. Previous studies have shown that p53 protein is highly expressed in SLE patients compared with the control group, and the activation of the p53 mediated signaling pathway was closely associated with the senescence of BM-MSCs from SLE patients . On the other hand, it is wildly known that the function of BM-MSCs is closely related to the pathology of ONFH . One study shows that p53 and Parkin that coregulate mitophagy in BM-MSCs could promote the repair of early steroid-induced ONFH . Hence, it is at least indicated that p53 and its associated cell cycle signaling pathways are associated with BM-MSCs in SLE-ONFH patients. What is more, as an effective drug for the treatment of SLE and ONFH, curcumin could improve the function of BM-MSCs at the same time and is more likely to play a role by targeting p53 and related cell cycle signal pathway proteins. Of course, it also needs to be further determined by in vivo or in vitro experiments. In the network above, the second central gene is VEGFA. VEGFA is a highly specific vascular endothelial growth factor, which is involved in the occurrence and development of many angiogenesis-dependent diseases. A study found that the serum VEGF level in the active SLE group was significantly higher than that in the control and inactive SLE groups . Moreover, it is reported that SNP rs2010963 of VEGFA is significantly associated with the risk of ONFH . It is worth mentioning that VEGFA is also significantly associated with other autoimmune diseases such as rheumatoid arthritis . IL6 and TNF are important cytokines involved in the immunopathological mechanism of SLE [63, 64]; they mediate the cellular inflammatory process. EGFR, epidermal growth factor receptor, which is Bsr I polymorphism is associated with SLE , and it also influences bone formation . The process of ONFH is manifested by the apoptosis of bone cells, CASP3, as an extremely important molecule for apoptosis, which participated in the apoptosis process of bone-related cells . MMP2 and MMP9 are members of the matrix metalloproteinase family. Studies have found that they are susceptibility genes for ONFH. MMP2 is known to be increased in osteoarthritic cartilage while MMP-9 regulates apoptosis of hypertrophic chondrocytes and is a key regulator of growth plate angiogenesis, thus rendering MMP-9 essential for normal bone development and remodeling . ESR-1 and IL1B are crucial metabolism and inflammatory molecules; they are involved in the process of osteoblast-osteoclast metabolism and inflammation [69, 70]. Previous studies have also shown that they are associated with susceptibility to SLE [71–73]. CCND1, a cyclin indispensable for cell proliferation, participates in cell cycle regulation . TGFB1, transforming growth factor-β, is associated with angiogenesis and vascularization and facilitates inflammation and joint damage . CXCL8 is a member of the CXC family of chemokines, which is expressed in the kidneys of SLE patients. A meta-analysis showed that SLE patients have higher circulating IL-8 levels than normal controls and IL-8 levels in SLE patients are influenced by age, region, and disease duration . FGF2 is a possible marker of kidney damage in SLE. SRC, the first identified viral oncogene, encoding tyrosine-protein kinase, has been reported to cause autoinflammatory bone disease . SERPINE1 is also known as the plasminogen activator inhibitor-1; previous studies have demonstrated that decreased fibrinolytic activity due to elevated plasminogen activator inhibitor-1 (PAI-1) levels correlates with ONFH pathogenesis .
4.3. Research Value and Bright Prospects for Revealing the Targets and Mechanism of Curcumin against SLE-ONFH
At present, the treatment of SLE depends on the organ involved in the disease. Patients with mild symptoms can use low-dose corticosteroids; however, moderate and severe SLE may require higher doses of corticosteroids or other immunosuppressive agents. The adverse effects and limitations of hormones and immunosuppressive agents make TCM become a candidate therapeutic drug because of their unique advantages. The characteristic of SLE is that the complement system is activated; curcumin could inhibit the complement cascade. Studies have found that curcumin inhibits the increase of matrix protein, glial fibrillary acidic protein, and vimentin in the hippocampus of lupus mice . As mentioned above, we found that curcumin may be able to treat SLE-ONFH by regulating TP53, involving the cell cycle pathway. As we all know, p53 protein plays a very important role in controlling the cell cycle, apoptosis, and DNA repair. There are abnormalities in BM-MSCs in SLE patients. Curcumin could have a positive effect on the treatment of SLE-ONFH by influencing the fate of BM-MSCs in SLE-ONFH patients. This provides an idea for the study of small molecular compounds specifically targeting p53 protein in SLE-ONFH disease. We initially confirmed this through molecular docking strategy. Combined with literature review, we found that curcumin could stabilize p53 by interaction with NAD(P)H: quinone oxidoreductase 1 in tumor-derived cell lines. Above all, according to our analysis, curcumin regulates p53 and influences the cell cycle. It may also regulate VEGFA, IL6, ERK1/2 (MAPK1/MAPK3), AKT1, MYC, CDKN1A, NF-KB, and other genes to effectively against SLE-ONFH because they have an important position in the gene-pathway network; this requires further research in the future, and we will continue to explore their roles in the process in SLE-ONFH.
The data that support the findings of this study are available from the corresponding authors upon request.
Conflicts of Interest
The authors report no conflicts of interest regarding this work.
Pan Kang worked on the investigation, methodology, writing (original draft), resources, and software. Zhiming Wu performed the investigation and methodology. Haibin Wang and Hongyu Tang performed the funding acquisition, project administration, and writing (review and editing). Yue Zhong, Zihao Wang, and Kun Xu worked on picture and software. Chi Zhou, Shaochuan Huo, Hai Guo, Songtao Li, Lingyun Liu, and Shuai Chen worked on data curation and advice. Pan Kang and Zhiming Wu contributed equally to this work.
This study was supported by the National Natural Science Foundation of China (Nos.: 81774339 and 82074462) and project funded by Guangzhou Science and Technology Plan Project (No.: 201707010319). Thanks are due to Mr. Wang Haibin, for the care and help he gave to the first author in the study and work.
- G. C. Tsokos, “Systemic lupus erythematosus,” The New England Journal of Medicine, vol. 365, no. 22, pp. 2110–2121, 2011.
- G. Yaniv, G. Twig, D. B.-A. Shor et al., “A volcanic explosion of autoantibodies in systemic lupus erythematosus: a diversity of 180 different antibodies found in SLE patients,” Autoimmunity Reviews, vol. 14, no. 1, pp. 75–79, 2015.
- R. A. Asherson, “Avascular necrosis of bone and antiphospholipid antibodies in systemic lupus erythematosus,” The Journal of Rheumatology, vol. 20, no. 2, pp. 284–288, 1993.
- R. Hisada, M. Kato, N. Ohnishi et al., “Antiphospholipid score is a novel risk factor for idiopathic osteonecrosis of the femoral head in patients with systemic lupus erythematosus,” Rheumatology (Oxford), vol. 58, no. 4, pp. 645–649, 2019.
- C. C. Mok, C. S. Lau, and R. W. Wong, “Risk factors for avascular bone necrosis in systemic lupus erythematosus,” British Journal of Rheumatology, vol. 37, no. 8, pp. 895–900, 1998.
- M. Y. Mok, V. T. Farewell, and D. A. Isenberg, “Risk factors for avascular necrosis of bone in patients with systemic lupus erythematosus: is there a role for antiphospholipid antibodies?” Annals of the Rheumatic Diseases, vol. 59, no. 6, pp. 462–467, 2000.
- M. Petri, “Musculoskeletal complications of systemic lupus erythematosus in the Hopkins Lupus Cohort: an update,” Arthritis Care and Research, vol. 8, no. 3, pp. 137–145, 1995.
- A. Rascu, K. Manger, H.-G. Kraetsch, J. R. Kalden, and B. Manger, “Osteonecrosis in systemic lupus erythematosus, steroid-induced or a lupus-dependent manifestation?” Lupus, vol. 5, no. 4, pp. 323–327, 1996.
- M. Sayarlioglu, N. Yuzbasioglu, M. Inanc et al., “Risk factors for avascular bone necrosis in patients with systemic lupus erythematosus,” Rheumatology International, vol. 32, no. 1, pp. 177–182, 2012.
- L. Cui, Q. Zhuang, J. Lin et al., “Multicentric epidemiologic study on six thousand three hundred and ninety five cases of femoral head osteonecrosis in China,” International Orthopaedics, vol. 40, no. 2, pp. 267–276, 2016.
- T. Kubo, K. Ueshima, M. Saito, M. Ishida, Y. Arai, and H. Fujiwara, “Clinical and basic research on steroid-induced osteonecrosis of the femoral head in Japan,” Journal of Orthopaedic Science, vol. 21, no. 4, pp. 407–413, 2016.
- Y. Sun, C. Q. Zhang, S. B. Chen, J. G. Sheng, D. X. Jin, and B. F. Zeng, “Treatment of femoral head osteonecrosis in patients with systemic lupus erythematosus by free vascularised fibular grafting,” Lupus, vol. 18, no. 12, pp. 1061–1065, 2009.
- L. Durcan, T. O'Dwyer, and M. Petri, “Management strategies and future directions for systemic lupus erythematosus in adults,” Lancet, vol. 393, no. 10188, pp. 2332–2343, 2019.
- L. Lyons, M. ELBeltagy, G. Bennett, and P. Wigmore, “The effects of cyclophosphamide on hippocampal cell proliferation and spatial working memory in rat,” PLoS One, vol. 6, no. 6, p. e21445, 2011.
- S.-H. Liu, W.-C. Chuang, W. Lam, Z. Jiang, and Y.-C. Cheng, “Safety surveillance of traditional Chinese medicine: current and future,” Drug Safety, vol. 38, no. 2, pp. 117–128, 2015.
- H. Tang, S. He, X. Zhang et al., “A network pharmacology approach to uncover the pharmacological mechanism of XuanHuSuo powder on osteoarthritis,” Evidence-based Complementary and Alternative Medicine, vol. 2016, Article ID 3246946, 10 pages, 2016.
- S. H. Venkatesha, B. M. Berman, and K. D. Moudgil, “Herbal medicinal products target defined biochemical and molecular mediators of inflammatory autoimmune arthritis,” Bioorganic & Medicinal Chemistry, vol. 19, no. 1, pp. 21–29, 2011.
- X.-X. Zhao, C. Peng, H. Zhang, and L.-P. Qin, “Sinomenium acutum: a review of chemistry, pharmacology, pharmacokinetics, and clinical use,” Pharmaceutical Biology, vol. 50, no. 8, pp. 1053–1061, 2012.
- Y. Tang, Y. Zhang, L. Li, Z. Xie, C. Wen, and L. Huang, “Kunxian capsule for rheumatoid arthritis: inhibition of inflammatory network and reducing adverse reactions through drug matching,” Frontiers in Pharmacology, vol. 11, 2020.
- Q. Xing, L. Fu, Z. Yu, and X. Zhou, “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.
- L. Chen, H.-M. Zhu, Y. Li et al., “Arsenic trioxide replacing or reducing chemotherapy in consolidation therapy for acute promyelocytic leukemia (APL2012 trial),” Proceedings of the National Academy of Sciences of the United States of America, vol. 118, no. 6, 2021.
- E. L. Dent, E. B. Taylor, H. R. Turbeville, and M. J. Ryan, “Curcumin attenuates autoimmunity and renal injury in an experimental model of systemic lupus erythematosus,” Physiological Reports, vol. 8, no. 13, p. e14501, 2020.
- S. Jin, C. Meng, Y. He et al., “Curcumin prevents osteocyte apoptosis by inhibiting M1-type macrophage polarization in mice model of glucocorticoid-associated osteonecrosis of the femoral head,” Journal of Orthopaedic Research, vol. 38, no. 9, pp. 2020–2030, 2020.
- A. Türkoğlu, “Effect of intraperitoneal curcumin instillation on postoperative peritoneal adhesions,” Medical Principles and Practice, vol. 24, no. 2, pp. 153–158, 2015.
- J. Zhao, J. Wang, M. Zhou, M. Li, M. Li, and H. Tan, “Curcumin attenuates murine lupus via inhibiting NLRP3 inflammasome,” International Immunopharmacology, vol. 69, pp. 213–216, 2019.
- A. Daina, O. Michielin, and V. Zoete, “SwissTargetPrediction: updated data and new features for efficient prediction of protein targets of small molecules,” Nucleic Acids Research, vol. 47, no. W1, pp. W357–w364, 2019.
- S. S. Fang, L. Dong, L. Liu et al., “HERB: a high-throughput experiment- and reference-guided database of traditional Chinese medicine,” Nucleic Acids Research, vol. 49, no. D1, pp. D1197–d1206, 2021.
- 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, 2014.
- R. Xue, Z. Fang, M. Zhang, Z. Yi, C. Wen, and T. Shi, “TCMID: traditional Chinese medicine integrative database for herb molecular mechanism analysis,” Nucleic Acids Research, vol. 41, no. Database issue, pp. D1089–D1095, 2013.
- J. S. Amberger and A. Hamosh, “Searching Online Mendelian Inheritance in Man (OMIM): a knowledgebase of human genes and genetic phenotypes,” Current Protocols in Bioinformatics, vol. 58, no. 1, 2017.
- S. Wei, X. Zhou, M. Niu et al., “Network pharmacology exploration reveals the bioactive compounds and molecular mechanisms of Li-Ru-Kang against hyperplasia of mammary gland,” Molecular Genetics and Genomics, vol. 294, no. 5, pp. 1159–1171, 2019.
- J. Zhao, C. Lv, Q. Wu et al., “Computational systems pharmacology reveals an antiplatelet and neuroprotective mechanism of Deng-Zhan-Xi-Xin injection in the treatment of ischemic stroke,” Pharmacological Research, vol. 147, p. 104365, 2019.
- J. Zhou, Q. Wang, Z. Xiang et al., “Network pharmacology analysis of traditional Chinese medicine formula Xiao Ke Yin Shui treating type 2 diabetes mellitus,” Evidence-based Complementary and Alternative Medicine, vol. 2019, Article ID 4202563, 15 pages, 2019.
- J. M. Ortiz-Rodriguez, C. Ortega-Ferrusola, M. C. Gil et al., “Transcriptome analysis reveals that fertilization with cryopreserved sperm downregulates genes relevant for early embryo development in the horse,” PLoS One, vol. 14, no. 6, p. e0213420, 2019.
- D. Szklarczyk, A. Franceschini, S. Wyder et al., “STRING v10: protein-protein interaction networks, integrated over the tree of life,” Nucleic Acids Research, vol. 43, no. D1, pp. D447–D452, 2015.
- J. L. Deng, Y. H. Xu, and G. Wang, “Identification of potential crucial genes and key pathways in breast cancer using bioinformatic analysis,” Frontiers in Genetics, vol. 10, 2019.
- D. Szklarczyk, J. H. Morris, H. Cook et al., “The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible,” Nucleic Acids Research, vol. 45, no. D1, pp. D362–d368, 2017.
- A. Martin, M. E. Ochagavia, L. C. Rabasa, J. Miranda, J. Fernandez-de-Cossio, and R. Bringas, “BisoGenet: a new tool for gene network building, visualization and analysis,” BMC Bioinformatics, vol. 11, no. 1, p. 91, 2010.
- X. Liu, Z. Hong, J. Liu et al., “Computational methods for identifying the critical nodes in biological networks,” Briefings in Bioinformatics, vol. 21, no. 2, pp. 486–497, 2020.
- M. Li, J. Wang, X. Chen, H. Wang, and Y. Pan, “A local average connectivity-based method for identifying essential proteins from the network level,” Computational Biology and Chemistry, vol. 35, no. 3, pp. 143–150, 2011.
- J. Wang, M. Li, H. Wang, and Y. Pan, “Identification of essential proteins based on edge clustering coefficient,” IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol. 9, no. 4, pp. 1070–1080, 2012.
- Y. Tang, M. Li, J. Wang, Y. Pan, and F.-X. Wu, “CytoNCA: a cytoscape plugin for centrality analysis and evaluation of protein interaction networks,” Biosystems, vol. 127, pp. 67–72, 2015.
- M. A. Alamri and M. A. Alamri, “Pharmacophore and docking-based sequential virtual screening for the identification of novel Sigma 1 receptor ligands,” Bioinformation, vol. 15, no. 8, pp. 586–595, 2019.
- D. Seeliger and B. L. de Groot, “Ligand docking and binding site analysis with PyMOL and Autodock/Vina,” Journal of Computer-Aided Molecular Design, vol. 24, no. 5, pp. 417–422, 2010.
- J. Lv, Y.-X. Zhu, Y.-Q. Liu, and X. Xue, “Distinctive pathways characterize A. actinomycetemcomitans and P. gingivalis,” Molecular Biology Reports, vol. 42, no. 2, pp. 441–449, 2015.
- G. A. Pavlopoulos, P. I. Kontou, A. Pavlopoulou, C. Bouyioukos, E. Markou, and P. G. Bagos, “Bipartite graphs in systems biology and medicine: a survey of methods and applications,” Gigascience, vol. 7, no. 4, pp. 1–31, 2018.
- D. Bu, Z. Su, J. Zou, M. Meng, and C. Wang, “Study of the mechanism underlying therapeutic effect of Compound Longmaining on myocardial infarction using a network pharmacology-based approach,” Biomedicine & Pharmacotherapy, vol. 118, p. 109234, 2019.
- J. Y. Zhang, C. L. Hong, H. S. Chen et al., “Target identification of active constituents of Shen Qi Wan to treat kidney yang deficiency using computational target fishing and network pharmacology,” Frontiers in Pharmacology, vol. 10, 2019.
- Y. Zhang, Z. Li, M. Yang et al., “Identification of GRB2 and GAB1 coexpression as an unfavorable prognostic factor for hepatocellular carcinoma by a combination of expression profile and network analysis,” PLoS One, vol. 8, no. 12, article e85170, 2013.
- P. M. Chen and G. C. Tsokos, “T cell abnormalities in the pathogenesis of systemic lupus erythematosus: an update,” Current Rheumatology Reports, vol. 23, no. 2, 2021.
- J. Liu, X. Wang, X. Yang, Q. Yan, S. Wang, and W. Han, “Investigating the role of angiogenesis in systemic lupus erythematosus,” Lupus, vol. 24, no. 6, pp. 621–627, 2015.
- T. Liu, M. Son, and B. Diamond, “HMGB1 in systemic lupus erythematosus,” Frontiers in Immunology, vol. 11, 2020.
- J.-M. Wang, W.-D. Xu, Z.-C. Yuan, Q. Wu, J. Zhou, and A.-F. Huang, “Serum levels and gene polymorphisms of angiopoietin 2 in systemic lupus erythematosus patients,” Scientific Reports, vol. 11, no. 1, 2021.
- B. Yousefi, M. Rastin, M. R. Hatef, J. Shariati, R. Alimohammadi, and M. Mahmoudi, “In vitro modulatory effect of dehydroepiandrosterone sulfate on apoptosis and expression of apoptosis-related genes in patients with systemic lupus erythematosus,” Journal of Cellular Physiology, vol. 234, no. 8, pp. 12676–12684, 2019.
- B. Boezio, K. Audouze, P. Ducrot, and O. Taboureau, “Network-based approaches in pharmacology,” Molecular Informatics, vol. 36, no. 10, 2017.
- S. Lee, B. Elenbaas, A. Levine, and J. Griffith, “p53 and its 14 kDa C-terminal domain recognize primary DNA damage in the form of insertion/deletion mismatches,” Cell, vol. 81, no. 7, pp. 1013–1020, 1995.
- L. J. Ko and C. Prives, “p53: puzzle and paradigm,” Genes & Development, vol. 10, no. 9, pp. 1054–1072, 1996.
- Z. Gu, J. Jiang, W. Tan et al., “p53/p21 Pathway Involved in Mediating Cellular Senescence of Bone Marrow- Derived Mesenchymal Stem Cells from Systemic Lupus Erythematosus Patients,” Clinical & Developmental Immunology, vol. 2013, article 134243, pp. 1–13, 2013.
- S. Xiang, Z. Li, and X. Weng, “The role of lncRNA RP11-154D6 in steroid-induced osteonecrosis of the femoral head through BMSC regulation,” Journal of Cellular Biochemistry, vol. 120, no. 10, pp. 18435–18445, 2019.
- F. Zhang, W. Peng, J. Zhang et al., “P53 and Parkin co-regulate mitophagy in bone marrow mesenchymal stem cells to promote the repair of early steroid-induced osteonecrosis of the femoral head,” Cell Death & Disease, vol. 11, no. 1, 2020.
- W. Ma, K. Xin, K. Chen et al., “Relationship of common variants in _VEGFA_ gene with osteonecrosis of the femoral head: A Han Chinese population based association study,” Scientific Reports, vol. 8, no. 1, p. 16221, 2018.
- J. Ramírez-Bello, D. Cadena-Sandoval, J. M. Fragoso et al., “The VEGFA -1154G/A polymorphism is associated with reduced risk of rheumatoid arthritis but not with systemic lupus erythematosus in Mexican women,” The Journal of Gene Medicine, vol. 20, no. 6, p. e3024, 2018.
- J. N. Stannard and J. M. Kahlenberg, “Cutaneous lupus erythematosus: updates on pathogenesis and associations with systemic lupus,” Current Opinion in Rheumatology, vol. 28, no. 5, pp. 453–459, 2016.
- J. Wenzel, “Cutaneous lupus erythematosus: new insights into pathogenesis and therapeutic strategies,” Nature Reviews Rheumatology, vol. 15, no. 9, pp. 519–532, 2019.
- C.-M. Huang, C.-H. Tsai, C.-L. Chen, C.-P. Chang, C.-C. Lai, and F.-J. Tsai, “Epidermal growth factor receptor (EGFR) gene Bsr I polymorphism is associated with systemic lupus erythematosus,” Lupus, vol. 13, no. 10, pp. 773–776, 2004.
- F. Yang, Z. W. Lin, T. Y. Huang et al., “Ligustilide, a major bioactive component of Angelica sinensis, promotes bone formation via the GPR30/EGFR pathway,” Scientific Reports, vol. 9, no. 1, 2019.
- H. Yu, P. Liu, W. Zuo et al., “Decreased angiogenic and increased apoptotic activities of bone microvascular endothelial cells in patients with glucocorticoid-induced osteonecrosis of the femoral head,” BMC Musculoskeletal Disorders, vol. 21, no. 1, p. 277, 2020.
- S. Grässel, “Expression profile of matrix metalloproteinase-2 and -9 and their endogenous tissue inhibitors in osteonecrotic femoral heads,” International Journal of Molecular Medicine, vol. 26, no. 1, pp. 127–133, 2010.
- A. B. Khalid and S. A. Krum, “Estrogen receptors alpha and beta in bone,” Bone, vol. 87, pp. 130–135, 2016.
- Y. K. Jung, Y. M. Kang, and S. Han, “Osteoclasts in the inflammatory arthritis: implications for pathologic osteolysis,” Immune Netw, vol. 19, no. 1, 2019.
- M. N. Drehmer, D. Andrade, I. A. Pereira et al., “Estrogen receptor alpha gene (ESR1) polymorphism can contribute to clinical findings in systemic lupus erythematosus patients,” Lupus, vol. 26, no. 3, pp. 294–298, 2017.
- G. G. Song, J.-H. Kim, Y. H. Seo, S. J. Choi, J. D. Ji, and Y. H. Lee, “Associations between interleukin 1 polymorphisms and susceptibility to systemic lupus erythematosus: a meta-analysis,” Human Immunology, vol. 75, no. 1, pp. 105–112, 2014.
- Y. Rupasree, “Association of estrogen receptor 1 (ESR1) haplotypes with risk for systemic lupus erythematosus among south Indians,” Indian Journal of Experimental Biology, vol. 53, no. 11, pp. 714–718, 2015.
- H. W. Yang, M. Chung, T. Kudo, and T. Meyer, “Competing memories of mitogen and p53 signalling control cell-cycle entry,” Nature, vol. 549, no. 7672, pp. 404–408, 2017.
- J. Belasco and N. Wei, “Psoriatic arthritis: what is happening at the joint?” Rheumatology and Therapy, vol. 6, no. 3, pp. 305–315, 2019.
- Y.-M. Mao, C.-N. Zhao, L.-N. Liu et al., “Increased circulating interleukin-8 levels in systemic lupus erythematosus patients: a meta-analysis,” Biomarkers in Medicine, vol. 12, no. 11, pp. 1291–1302, 2018.
- K. Abe, A. Cox, N. Takamatsu et al., “Gain-of-function mutations in a member of the Src family kinases cause autoinflammatory bone disease in mice and humans,” Proceedings of the National Academy of Sciences of the United States of America, vol. 116, no. 24, pp. 11872–11877, 2019.
- H. Kim, “Significant associations of PAI-1 genetic polymorphisms with osteonecrosis of the femoral head,” BMC Musculoskeletal Disorders, vol. 12, no. 1, p. 160, 2011.
- S. Foxley, M. Zamora, B. Hack et al., “Curcumin aggravates CNS pathology in experimental systemic lupus erythematosus,” Brain Research, vol. 1504, pp. 85–96, 2013.
Copyright © 2021 Pan Kang 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.