Potential Therapeutic Mechanism of Radix Angelicae Biseratae and Dipsaci Radix Herb Pair against Osteoarthritis: Based on Network Pharmacology and Molecular Docking

Xi, Yujiang; Zhao, Ting; Shi, Mingqin; Zhang, Xiaoyu; Bao, Yanyuan; Gao, Jiamei; Shen, Jiayan; Wang, Hui; Xie, Zhaohu; Wang, Qi; Li, Zhaofu; Qin, Dongdong

doi:https://doi.org/10.1155/2023/2140327

Evidence-Based Complementary and Alternative Medicine

On this page

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

Special Issue

Real-World Big Data Processing and Analysis for Traditional Chinese Medicine 2022

View this Special Issue

Research Article | Open Access

Volume 2023 | Article ID 2140327 | https://doi.org/10.1155/2023/2140327

Potential Therapeutic Mechanism of Radix Angelicae Biseratae and Dipsaci Radix Herb Pair against Osteoarthritis: Based on Network Pharmacology and Molecular Docking

Yujiang Xi,¹Ting Zhao,¹Mingqin Shi,²Xiaoyu Zhang,¹Yanyuan Bao,¹Jiamei Gao,¹Jiayan Shen,¹Hui Wang,¹Zhaohu Xie,²Qi Wang,²Zhaofu Li,^1,2and Dongdong Qin²

Academic Editor: Xuezhong Zhou

Received21 Oct 2022

Revised13 Nov 2022

Accepted17 Mar 2023

Published14 Apr 2023

Abstract

Background. A major contributor to older disability is osteoarthritis. Radix Angelicae Biseratae (known as Duhuo in China, DH, the dried rhizome of Angelica pubescens) and Dipsaci Radix (known as Xuduan in China, XD, the dried rhizome of Dipsacus asper Wall) herb pair (DXHP) is widely used to treat osteoarthritis, but the underlying molecular mechanisms still have not been revealed. This research aimed to illustrate the therapeutic mechanism of DXHP against osteoarthritis through the techniques of network pharmacology and molecular docking. Methods. Gene targets for osteoarthritis and active ingredients for DXHP were screened based on the pharmacology public database and the gene-disease target database. The software program Cytoscape was used to visualize the active chemical target-disease gene network. The STRING biological information website was used to investigate protein interactions. On the Metascape bioinformatics website, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were carried out. The molecular docking of the important chemicals and primary targets identified by the aforementioned screening was performed using Autodock software. Results. Twenty-six active substances from the DXHP that had strong connections to 138 osteoarthritis-related targets were screened out. According to network analysis, TNF, GAPDH, IL-6, AKT-1, IL-1B, and VEGFA are prospective therapeutic targets, while osthole, cauloside A, ammidin, angelicone, beta-sitosterol, and asperosaponin VI may be significant active components. 1705 biological processes (BP), 155 molecular functions (MF), and 89 cellular components (CC) were identified by GO analysis. KEGG analysis indicated that IL-17, NF-kappa B, HIF-1, MAPK, and AGE-RAGE signaling pathways are potentially involved. Molecular docking showed that cauloside A, osthole, and β-sitosterol have excellent binding activity with main targets. Conclusions. This study comprehensively illuminated the active ingredients, potential targets, primary pharmacological effects, and relevant mechanisms of the DXHP in the treatment of OA. These findings provide fresh thoughts into the therapeutic mechanisms of the main active ingredients of DXHP and provide a reference for further exploration and clinical applications of DXHP.

1. Introduction

Destruction of cartilage, remodeling of the subchondral bone, the development of osteophytes, and synovial inflammation are all symptoms of osteoarthritis (OA), a degenerative joint disease [1–3]. OA is directly associated to heredity, age, obesity, injury, and chronic inflammation [4–7]. As populations of aging and obese individuals have grown in size, the prevalence of OA has also increased and is now considered a major public health problem worldwide [8]. According to recent studies, the global prevalence of OA grew by 48% between 1990 and 2019 [3, 9], and that OA affects more than 500 million persons, occurring most commonly in postmenopausal women over age 50 [10]. Some guidelines indicate that nonsteroidal anti-inflammatory drugs (NSAIDs) are the most crucial medications for OA [11, 12]. This class of drugs achieves its anti-inflammatory and analgesic effects by inhibiting cyclooxygenase (COX), a key enzyme in the metabolism of arachidonic acid (AA) [13]. However, the use of COX inhibitors alone is not ideal. Studies have shown that AA can also be metabolized by lipoxygenase (LOX) to produce inflammatory substances such as leukotrienes [14]. Therefore, single introduction of COX inhibitors cannot limit production of inflammatory mediators but instead stimulates and increases the release of inflammatory mediators in the LOX pathway, resulting in adverse reactions [15]. The current mainstream treatment options for OA have limited effects on delaying the development of joint inflammation. In the long term, such drugs cause serious damage to the functions of the digestive tract, liver, and kidneys [16, 17]. In addition, long-term use accelerates the progression of arthritis [18]. Therefore, research into alternative therapies that can contribute to the development of new drugs could improve treatment options for OA.

Traditional Chinese medicine (TCM), as a supplementary therapy, has accumulated rich theoretical knowledge that can potentially help advance clinical prevention and treatment strategies for OA [19]. The earliest research on OA in TCM can be traced back thousands of years to the ancient text “Huangdi Neijing [20].” In recent years, TCM has been suggested to have anti-inflammatory, anti-apoptotic, anti-oxidative, anti-metabolic, and proliferative effects in the treatment of OA diseases [21]. Herb pair refers to the combination of drugs that can be used simultaneously in TCM clinical to enhance efficacy or reduce adverse reactions. Compared with TCM formulas, herb pairs have a more defined synergistic mechanism [22]. Radix Angelicae Biseratae (known as Duhuo (DH), the dried rhizome of Angelica pubescens) and Dipsaci Radix (known as Xuduan (XD), the dried rhizome of Dipsacus asper Wall) are a classic herbal drug combination that has been widely used in TCM. Clinical statistical studies have shown that DXHP is one of the most frequently used Chinese medicine pairs in the treatment of OA [23]. In modern pharmacological research, the active ingredients in DH are linked to antioxidant, anti-inflammatory, and analgesic effects [24], while the active components of XD purportedly protect against bone loss, promote cartilage formation, and improve bone metabolism [25]. In addition, the pharmacological compatibility mechanism of DXHP is similarly reflected in the TCM theory. DXHP originated from Duhuo Xuduan decoction in the ancient Chinese medical book “Waitai Miyao.” Chinese Medicine Pharmacopoeia points out that DH helps relieve symptoms of rheumatism, particularly symptoms associated with dampness, swelling, and pain. XD supposedly tonifies the liver and kidneys, nourishes the muscles and bones, and repair fracture [26]. Hence, it is highly essential to research the pharmacological mechanism of DXHP therapy as a possible means of improving the treatment of OA. However, the molecular mechanism of DXHP in the treatment of OA remains unclear.

In recent years, network pharmacology, which combines multidisciplinary knowledge and methods, has the characteristics of being systematic and holistic. Molecular docking can calculate the binding energies between ligands and receptors and predict reasonable binding mode and has been used to explore the molecular process of drug active ingredients acting on the human body. Therefore, the aim of this research is to clarify the molecular underpinnings behind the effects of DXHP on OA using network pharmacology and molecular docking technologies. The specific technical route adopted in this study is shown in Figure 1.

2. Materials and Methods

2.1. Acquisition of Bioactive Compounds of DXHP

To obtain more comprehensive compositional data, the active ingredients of DH and XD were retrieved from three herbal pharmacological databases, including the TCMSP data platform (https://tcmsp-e.com/) [27], the ETCM data platform (https://www.tcmip.cn/ETCM/index.html) [28], and the SymMap data platform (https://www.symmap.org/) [29]. According to previous studies and recommended standards [30], oral bioavailability (OB) ≥30% is regarded as having an excellent absorption. As the selection criteria for “drug-like” compounds in traditional Chinese medicinal materials, DL ≥0.18 is thought to be appropriate for drug design. Therefore, OB and DL criteria were chosen to investigate compounds that meet basic pharmacokinetic criteria for effective drug development [31]. Finally, we searched PubMed and CNKI databases for compounds that did not meet the criteria but had been proved to have significant therapeutic effects on OA.

2.2. Target Prediction for Compounds

The active ingredients obtained after the aforementioned screening were queried by PubChem (database https://pubchem.ncbi.nlm.nih.gov/) for SMILES codes of their potential ingredients, which were entered into Swiss Target Prediction (https://www.swisstargetprediction.ch/) for further prediction of active ingredient targets, setting the screening parameter criteria to the probability value >0.1. To normalize gene nomenclature and species and avoid over-annotation of homologous proteins, these target genes were transformed into matching gene symbols that corresponded to the “Homo sapiens” species using the UniProt information system (https://www.uniprot.org/) [32].

2.3. Collecting Herb-Disease Genes

We retrieved OA-related genes from the GEO data platform (https://www.ncbi.nlm.nih.gov/geo/) [33] database using “osteoarthritis” as a keyword search. The datasets were processed using the robust multiarray average algorithm for background correction and matrix data normalization. The screening conditions for significantly differential genes were and |log FC| >1.5. The differentially expressed genes for OA were obtained. Then, we search the OA-related genes in the following databases: GeneCards (https://www.genecards.org) [34], the TTD information system (https://db.idrblab.net/ttd) [35], and the OMIM information system (https://www.omim.org) [36]. For these searches, the species source was limited to “Homo sapiens.” The above-obtained genes were merged, and the duplicated genes were deleted to obtain the defined disease targets of osteoarthritis. Subsequently, the OA-related targets were matched with the active compound targets of DXHP to obtain the target of DXHP in the treatment of OA. The Venn diagram was drawn using the Venny2.1 plotter (https://bioinfogp.cnb.csic.es/tools/venny/index.html).

2.4. Protein-Protein Interaction (PPI) Network

Based on the STRING Bioinformatics system (https://www.string-db.org/), a protein-protein interaction (PPI) network was crafted for the common targets of compounds and diseases [37]. The choice of species was human, and the confidence level was set to >0.4. Others parameters keep the default settings. To identify the primary target genes for OA therapy that have a strong correlation with DXHP, the topological properties of common targets were analyzed using the CytoNCA topology analyzer (A plugin for Cytoscape). The network analyzer tool was applied for topology analysis, referring to three parameters, namely the degree centrality value (DC), the betweenness centrality value (BC), and the closeness centrality value (CC). The greater the DC of the node, the higher the importance of the node in a PPI network. Targets with all three parameters above average were selected and ordered by the degree centrality value.

2.5. Protein-Protein Interaction Enrichment MCODE Analysis

MCODE, as a clustering algorithm, helps to capture modules with high-quality biological processes in a large PPI spherical network and helps to discover a subset of targets that are closely related to that functional module. MCODE subcluster enrichment analysis was performed on the PPI network. PPI enrichment analysis used the following databases: STRING, InWebIM, OmniPath, and BioGRID. Only the STRING Bioinformatics system (Physical Score >0.132) and BioGRID data are used in the Physical Interaction module. A subset of the proteins in the resultant network physically interacts with at least one other item on the list. The MCODE methodology was used to find components of highly linked networks in subnetworks with 3 to 500 proteins.

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

The Metascape database (https://metascape.org/) is a powerful broad-coverage and fast-updating gene function annotation analysis tool that can analyze a large number of gene or protein functions. In order to effectively research the biological ontology of DXHP in the regulation of OA and to clarify the biological process of each core target protein and its function in signaling pathway transduction, GO functional annotation and KEGG pathway enrichment analysis were performed based on the Metascape biological system. The common targets of DXHP and OA obtained in the above screening were imported into the target gene list, with confined to just human species and correction of all target genes to their recognized gene symbols. Multiple comparisons were performed using Benjamini–Hochberg’s FDR correction to avoid false positives. Bonferroni-corrected values <0.01 for GO and KEGG terms were considered significant. The minimal count for the KEGG analysis was 3, and the enrichment factor was more than 1.5. Finally, the top 20 items were picked, and the annotated chart was formed on the Bioinformatics (https://www.bioinformatics.com/) platform for visualization.

2.7. Molecular Docking Verification of Compound Target

In network pharmacology, molecular docking is a critical tool for verifying compound-target interactions. It works by combining proteins of known targets obtained from network pharmacology with small compounds of active ingredients in natural drugs and then evaluating the strength and activity of the binding [38]. The mol2 format files of the main active components of DXHP were searched and obtained from the PubChem database in advance. Next, the selected protein’s best resolution 3D structure was acquired in the PDB database (https://www.rcsb.org/) and their PDBID files were downloaded. For these 3D structures, the species source was limited to “Homo sapiens” and contains the crystal structure of the complete pocket. MGLtool 1.5.7 was used to process the protein by adding hydrogen, calculating the charge, merging the nonpolar hydrogen, and saved as a PDBQT format file as a docking ligand. Gridbox coordinates and docking box sizes were set, and molecular docking was performed using Autodock Vina1.1.2. The lowest binding energy score conformation was selected, with lower binding energy scores indicating better docking activity and strength. The molecular docking structures were demonstrated using Pymol 2.3.

3. Results

3.1. Screening of DXHP Active Compounds

We obtained 106 active components of DH and 63 active components of XD. After screening by OB and DL standards, 20 active ingredients with good ADME properties were obtained for the subsequent study. In addition to this, DH and XD were found to contain some pharmacologically active ingredients, which were excluded because their OB and DL values were less than the screening criteria. Thus, according to previous literature reports, we included a total of six active ingredients such as osthole, columbianadin, umbelliferone, asperosaponin VI, ursolic acid, and loganin [39, 40]. A total of 26 DXHP active compounds and their corresponding 463 targets were screened (Table 1). The results show that a single compound can regulate multiple targets, indicating that DXHP has multicomponent and multitarget components. The active compounds of DXHP and corresponding genes are shown in Figure 2.

3.2. Drug-Disease Core Target Acquisition

In the GEO database, we identified the dataset that met the criteria as GSE169077. The chip contained six osteoarthritis cartilage tissue samples as the experimental group and five normal knee cartilage tissue samples as the control group. The results of differential expression gene analysis of osteoarthritis showed that 792 genes were found to be differentially expressed. The heatmap and Limma package in the R language were used to draw the dataset volcano map of differential genes (Figure 3). The GeneCards database yielded 3647 OA-related genes in total, and 885 OA-related genes were obtained using a relevance score >1 as the screening criterion. The TTD database attained 33 OA-related genes in total, whereas the OMIM database attained 30 OA-related genes. The above gene sets were combined, and duplicate values were removed to obtain a total of 1778 OA-related genes. The screened drug targets and OA-related genes were input into Venny 2.1.0 plotter, a Venn figure was constructed, and a total of 138 common targets were obtained (Figure 4).

3.3. PPI Network Construction and Centiscape Analysis

The herb-disease common targets were imported into the STRING Bioinformatics system. Species option was selected as “Homo sapiens,” the PPI network was constructed (Figure 5), and the CSV file was entered into Cytoscape software 3.9.1 for network topology analysis using the network Stats tool and the centiscape plugin for the construction of the three parameter correlations such as DC, BC, and CC. 138 nodes and 1599 edges made up the network, which had an average connectedness of 23.70. A total of 25 core targets were selected through the centiscape plugin by selecting genes with greater than average DC, BC, and CC parameters as key targets and ranking them by the degree centrality value. The 10 highest ranked core targets are listed in Table 2. DH and XD have 155 overlapping targets, 69 of which are associated with OA (Figure 6).

3.4. GO Functional Annotation

GO functional annotation analysis showed that the enrichment results included 1705 biological processes (BPs), 155 molecular functions (MFs), and 89 cellular components (CCs). Most GO annotations of BP were contained to regulation of inflammatory response, response to hormone, inflammatory response, positive regulation of locomotion, cellular response to lipid, cellular response to organic cyclic compound, positive regulation of cell death, protein phosphorylation, positive regulation of MAPK cascade, and cell activation. MF annotation was mainly involved in protein kinase activity, endopeptidase activity, transcription factor binding, nuclear receptor activity, cytokine receptor binding, protease binding, protein homodimerization activity, transmembrane receptor protein kinase activity, oxidoreductase activity, and phosphatase binding. The majority of CC annotations were included to the extracellular matrix, membrane raft, receptor complex, cell body, organelle outer membrane, ficolin-1-rich granule lumen, endolysosome, protein kinase complex, side of the membrane, and endoplasmic reticulum lumen. Based on their Q-value (the Q-value was used for multiple testing, it was calculated using the Benjamini–Hochberg procedure, and higher Q-values showed the greater GO term enrichment), the top 10 terms of BP, CC, and MF were rated. The terms are presented in Figure 7.

3.5. KEGG Pathway Analysis

KEGG pathway analysis enriched 179 signaling pathways after FDR correction. The top 20 pathways included pathways in cancer, the AGE-RAGE signaling pathway in diabetic complications, the IL-17 signaling pathway, the NF-kappa B signaling pathway, lipid and atherosclerosis, the TNF signaling pathway, the MAPK signaling pathway, hepatitis B, the toll-like receptor signaling pathway, prostate cancer, the HIF-1 signaling pathway, the PI3K-Akt signaling pathway, microRNAs in cancer, apoptosis, fluid shear stress and atherosclerosis, chagas disease, the VEGF signaling pathway, human cytomegalovirus infection, the FoxO signaling pathway, and osteoclast differentiation. Based on their Q-value, the top 20 critical signaling pathways are presented in Figure 8. The figures of the two inflammation-related signaling pathways with the highest number of core genes enriched are shown in Figures 9 and 10.

3.6. PPI Enrichment MCODE Analysis

Core genes were identified through network construction and MCODE analysis using the complete datasets for the independent enrichment analysis of gene clusters. The enrichment analysis of biological processes was used for each MCODE component (Figure 11 and Table 3). MCODE analysis showed that the core genes of AKT1, IL1B, TLR9, VEGFA, MMP2, PTGS1, NOS2, IL6, etc., are mainly involved in the biological process of regulating the inflammatory response in OA. The core genes of ESR1, MTOR, MAPK1, CASP3, CCL2, PGR, etc., are mainly involved in the biological process of regulating defense response in OA. The three terms with the highest scores based on the values were used as the functional description of the corresponding components (Table 4).

3.7. Molecular Docking Verification

To further validate the reliability of the binding of key targets and components screened by the network analysis, molecular docking verification of some key targets and important active components was carried out using Autodock Vina1.2.0 software. It is largely accepted that the steadier the binding structure, the smaller the binding energy of the receptor-ligand docking and therefore the higher the likelihood of interaction between the two, with a binding energy < −5.0 kcal/mol as the screening criterion. Binding energies < −5.0 kcal/mol indicate potential activity, and docking with binding energies < −7.0 kcal/mol is extremely stable [41]. The top six active ingredients in the component-target network were selected for molecular docking with the top six targets in the PPI network. The results showed that about 83% of targets and active components exhibited binding ability and 16% exhibited extremely strong binding ability. These molecular docking findings align with earlier network screening conclusions, which indirectly validate the treatment ability of DXHP on OA and demonstrate the reliability of network pharmacology applied to this study. The docking results of binding affinity and detailed compound-target interactions are presented in Table 5 and Figure 12.

(a)

(b)

(c)

(d)

4. Discussion

OA has a high prevalence in the elderly population, and modern medical treatment options are primarily symptomatic, with no options available to curb disease progression [42]. The advantages of DXHP in relieving OA symptoms as well as delaying OA disease progression have been demonstrated in previous studies [43], but its mechanism of action remains incompletely elucidated. Therefore, we investigated the mechanism of action of DSHP on OA in a more systematic way using a network pharmacology and the molecular docking approach.

According to this investigation, DXHP’s primary therapeutic ingredients for OA include substances such as osthole, asperosaponin VI, angelicone, beta-sitosterol, ammidin, and cauloside A, among which β-sitosterol is the overlapping component of DH and XD. In recent years, the pharmacological effects of asperosaponin VI in anti-inflammatory analgesia, prevention of osteoporosis, neuroprotection, and anti-apoptosis have attracted the attention of many scholars [44]. Many in vivo and cell experiments have shown that asperosaponin VI can phosphorylate ERK/2 protein to promote the expression of osteogenic genes such as ALP, OCN, COL1 and RUNX2, which is closely related to the PI3K/Akt signaling pathway and the MAPK signaling pathway [45]. In addition, asperosaponin VI could significantly reduce MDA, TNF-α, IL6, and IL10 by significantly inhibiting oxidative stress and inflammatory response in tissues [46]. Osthole, one of the main active components of DH, is thought to improve bone metabolism and promote osteoblast activation [47]. Recent studies on the anti-inflammatory mechanism of osthole have found that osthole can inhibit the formation and resorption activity of osteoclasts by inhibiting the activation of NF-κB and NFATc1 and reducing the expression of osteoclast-specific genes such as CTSK, MMP-9, TRAP, integrin β3, C-SRC, and NFATc1 [48]. β-sitosterol has a strong down-regulate effect on proinflammatory factors, such as IL-1β, IL-6, and TNF-α and quenches ROS produced by the human body through antioxidant activity. It can also alleviate inflammation through eosinophil percolation [49]. Chen et al. assessed the knee joints of OA rabbits by morphological and histological methods and found that β-sitosterol could significantly inhibit the secretion of matrix metalloproteinases (MMPs) and inhibit the degradation of cartilage [50].

The observations of the coincide target PPI analysis demonstrated that the core targets of DXHP in the healing of OA encompassed multiple targets such as TNF, GAPDH, IL6, AKT1, IL1B, VEGFA, CASP3, STAT3, MMP9, and PTGS2. TNF-α has been demonstrated in many diseases, including osteoarthritis, autoimmune diseases, ankylosing spondylitis, insulin resistance, psoriasis, nephropathy, and cancer [51]. TNF-α is associated with many cytokines, and it has been found that TNF-α implicated in angiogenesis by synergistic induction of VEGF production with IL-1β and IL-6 [52]. Nie et al. found that TNF-α could regulate the harmony between Treg cells and TH17 and TH1 in the joint synovium through FOXP3 dephosphorylation [51]. When cells are stimulated with NO, GAPDH will be nitrosylated, bound to E3 ubiquitin ligase Siah1, and undergo nuclear translocation and apoptosis [53]. The results of molecular docking suggest that DXHP may achieve the purpose of treating OA by binding to GAPDH. Studies have shown that [54] knockdown of the IL-6 gene in an OA rat model can lead to inhibition of MMP13 expression and secretion, which may be related to the inhibition of c-fos/ap-1-mediated inflammatory stimulation in OA chondrocytes. IL-6 can also cause an increase in the expression of MMP9, prevents the formation of type II collagen and proteoglycans, accelerates the degradation of extracellular matrix, and influences bone resorption by activating osteoclasts [55]. Shahine and Elhadidi found that the expression amount of the IL-1β gene in OA organisms was positively correlated with the pain index [56]. Correspondingly, high expression of IL-1β was detected in the synovial membrane and fluid of OA organisms and was positively correlated with OA disease [57]. Apoptotic chondrocytes are essential for the progression of OA. AKT1 is a key down-regulate gene kinase involved in the PI3K pathway. Researchers have demonstrated that the PI3K/AKT pathway that involves AKT1 prevents chondrocyte apoptosis [58]. TP53 inhibits DNA replication, induces apoptosis, and accelerates cartilage degradation [59]. In addition, researchers have demonstrated that the VEGFA is closely associated with many pathological responses, such as osteophyte formation and cartilage degeneration in OA [60]. The expression of VEGFA stimulates the division of vascular endothelial cells to promote angiogenesis, accelerate the exchange of nutrients in the knee joint and the metabolism of inflammatory products, and promote the growth of cartilage synovial cells.

The KEGG enrichment analysis showed that the main enriched inflammation-related signaling pathways included AGE-RAGE, IL-17, NF-kappa B, MAPK, TLR, HIF-1, PI3K-Akt, and the VEGF signaling pathway. The other signaling pathways enriched included pathways in cancer, lipid and atherosclerosis, hepatitis B, prostate cancer, microRNAs in cancer, apoptosis, fluid shear stress and atherosclerosis, chagas disease, human cytomegalovirus infection, FoxO signaling pathway, and osteoclast differentiation. The MCODE analysis revealed that the regulatory response to inflammation was the most significant biological process in treating OA by DXHP. The NF-κb pathway is closely related to cartilage destruction in OA and targeted therapy of OA [61]. Previous research has showed that the secretion level of NF-κb in joint synovial fluid and peripheral blood of OA patients is increased, and the degree of up-regulation of the NF-κb pathway is positively correlated with the degree of cartilage erosion and destruction [62]. After phosphorylation, IκBα dissociates from NF-κB, and NF-κB is activated and enters the nucleus, promoting the synthesis and secretion of TNF-α and IL-6 [63, 64] and eventually leading to the degeneration of articular cartilage [65, 66]. In addition, IL-17 can trigger the release of chemokines, cytokines, antimicrobial peptides, and matrix metalloproteinases from mesenchymal and bone marrow cells [67]. Bai et al. [68] found that the level of IL-17 in the serum of OA patients increased and was positively correlated with the severity of OA. Due to the lack of capillaries in articular cartilage, the cartilage microenvironment is essentially a hypoxic environment. The main hypoxia-inducible factor (HIF) of articular chondrocytes is a transcription factor, which is also the main mediator of homeostatic response that enables cells to survive under hypoxic conditions [69]. In OA, the expression of HIF-1α decreases [70], and the loss of HIF-1α can up-regulate the expression of MMP13, degrade col2A1 and ACAN, promote chondrocyte degradation, and promote the development of OA [71]. The advanced glycation end product receptor (RAGE) is secreted in a variety of cells, including macrophages and mast cells [72]. AGE is a ligand of RAGE and is the end product of glycosylation of proteins and sugar [73]. RAGE is essential for the induction of several inflammatory genes as well as important signaling pathways linked to proinflammatory responses. Although AGE induces inflammation by exciting NF-κB and MAPK in a variety of cells including osteocytes [74], it can also up-regulate the expression of PGE2 and NO via the MAPK pathway and induces the chondrocytes inflammatory [75].

As an herb pair, DH and XD have many overlapping targets and active ingredients, and the compatibility mechanism of DXHP in treating OA may be related to this. β-sitosterol, as its overlapping active ingredient, is associated with anti-inflammatory and antioxidant effects, immune regulation, and bone metabolic balance [76]. Furthermore, the aforementioned KEGG and MCODE analyses revealed that the overlapping targets such as IL6, IL-1β, AKT1, VEGFA, MMP13, and STAT3 are vital links in inflammatory signalling pathways such as AGE-RAGE, IL-17, MAPK, and NF-kappa B and are extensively involved in the biological processes of regulating the inflammatory response in OA. This finding may help us to better understand the compatibility mechanism of DXHP from a molecular perspective.

Molecular docking tests were carried to further confirm the molecular mechanism of DXHP in the treatment of OA. The results revealed that osthole, angelicone, cauloside A, and β-sitosterol have excellent binding activity with multiple key targets, such as TNF, GAPDH, IL-6, AKT1, IL-1β, and VEGFA. This also explains the therapeutic mechanism of DXHP from another perspective. These results offer insight for the subsequent application of network pharmacology methods to improve the efficiency of natural medicine ingredient development, which could facilitate the development of new high-efficiency, low-toxicity, multi-target OA drugs capable of improving symptoms and delaying disease progression.

Compared to recently published network pharmacology studies of similar diseases [77], we have expanded the source of the drug and disease database and introduced differential genes from the GEO microarray, making the source more comprehensive. In order to capture modules with high-quality biological processes in large PPI networks, we apply the MCODE clustering algorithm, which allows for a further in-depth interpretation of the target network. In addition, we have investigated the overlapping active ingredients and targets of this herb pair to elucidate the compatibility mechanisms of DXHP from a new perspective.

Despite the study’s advantages, several limitations should be mentioned. Firstly, the sources of component and target data are scattered and vary widely between databases, resulting in insufficiently accurate data sources [78]. Secondly, the validation of the findings is only based on computer simulations of molecular docking, which has limitations in terms of convincingness [79]. Therefore, the key pathways and targets identified in this study need to be validated in more in-depth in vivo and in vitro experiments.

5. Conclusion

In this study, we found that DXHP is a valuable TCM herb pair for treating patients with OA through multiple components, targets, and pathways. Its pharmacological mechanism could be through the AGE-RAGE, IL-17, NF-kappa B, MAPK, TLR, HIF-1, PI3K-Akt, and VEGF signaling pathways to alleviate OA. Hopefully, our research may provide a scientific basis for the prodrug discovery of its natural ingredients and the identification of therapeutic targets in the future.

Data Availability

The original contributions presented in the study are included within the article/supplementary material. The data supporting the findings of the current study are available from the corresponding authors upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Yujiang Xi, Ting Zhao, and Mingqin Shi contributed equally to this work.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (31960178, 82160923, 81960870, and 82260928); Applied Basic Research Programs of Science and Technology Commission Foundation of Yunnan Province (2019FA007); Key Laboratory of Traditional Chinese Medicine for Prevention and Treatment of Neuropsychiatric Diseases, Yunnan Provincial Department of Education; Scientific Research Projects for High-level Talents of Yunnan University of Chinese Medicine (2019YZG01); Young Top-Notch Talent in 10,000 Talent Program of Yunnan Province (YNWR-QNBJ-2019-235); National Science and Technology Innovation 2030 Major Program (2021ZD0200900); Construction Project of National Traditional Chinese Medicine Clinical Research Base (2018 No. 131); Clinical Cooperative Project of Chinese and Western Medicine for Major and Knotty Diseases; the Yunnan Provincial Key Laboratory Construction Project Funding; Yunnan Provincial Key Laboratory of Chinese Medicine Rheumatology and Immunology; the Yunnan Provincial Ten Thousands Program Famous Doctor Special; Yunnan Province Qingguo Wang Expert Workstation Construction Project (202005AF150017); Yunnan Applied Basic Research Projects-Union Foundation (2019FF002(-031)); and Scientific Research Fund Project of the Yunnan Provincial Department of Education (2021Y461).

References

L. J. Wang, N. Zeng, Z. P. Yan, J. T. Li, and G. X. Ni, “Post-traumatic osteoarthritis following ACL injury,” Arthritis Research and Therapy, vol. 22, no. 1, p. 57, 2020.
View at: Publisher Site | Google Scholar
D. Bryazka, M. B. Reitsma, M. G. Griswold et al., “Population-level risks of alcohol consumption by amount, geography, age, sex, and year: a systematic analysis for the Global Burden of Disease Study 2020,” The Lancet, vol. 400, no. 10347, pp. 185–235, 2022.
View at: Publisher Site | Google Scholar
D. J. Hunter and S. Bierma-Zeinstra, The Lancet, vol. 393, no. 10182, pp. 1745–1759, 2019.
View at: Publisher Site
E. Macfarlane, M. J. Seibel, and H. Zhou, “Arthritis and the role of endogenous glucocorticoids,” Bone Res, vol. 8, no. 1, p. 33, 2020.
View at: Publisher Site | Google Scholar
D. A. Young, M. J. Barter, and J. Soul, “Osteoarthritis year in review: genetics, genomics, epigenetics,” Osteoarthritis and Cartilage, vol. 30, no. 2, pp. 216–225, 2022.
View at: Publisher Site | Google Scholar
P. Kulkarni, A. Martson, R. Vidya, S. Chitnavis, and A. Harsulkar, “Pathophysiological landscape of osteoarthritis,” Advances in Clinical Chemistry, vol. 100, pp. 37–90, 2021.
View at: Publisher Site | Google Scholar
T. A. Perry, M. J. Parkes, R. J. Hodgson, D. T. Felson, N. K. Arden, and T. W. O'Neill, “Association between Bone marrow lesions & synovitis and symptoms in symptomatic knee osteoarthritis,” Osteoarthritis and Cartilage, vol. 28, no. 3, pp. 316–323, 2020.
View at: Publisher Site | Google Scholar
G. Aubourg, S. J. Rice, P. Bruce-Wootton, and J. Loughlin, “Genetics of osteoarthritis,” Osteoarthritis and Cartilage, vol. 30, no. 5, pp. 636–649, 2022.
View at: Publisher Site | Google Scholar
D. Veret, C. Jorgensen, and J. M. Brondello, “Osteoarthritis in time for senotherapeutics,” Joint Bone Spine, vol. 88, no. 2, Article ID 105084, 2021.
View at: Publisher Site | Google Scholar
I. A. Szilagyi, J. H. Waarsing, D. Schiphof, J. B. J. van Meurs, and S. M. A. Bierma-Zeinstra, “Towards sex-specific osteoarthritis risk models: evaluation of risk factors for knee osteoarthritis in males and females,” Rheumatology, vol. 61, no. 2, pp. 648–657, 2022.
View at: Publisher Site | Google Scholar
Q. Y. Shi, X. Y. Tan, and W. Y. Li, “Joint surgery group of orthopedic branch of Chinese medical as- sociation, orthopedic Expert committee of Wu jieping medical foundation. Consensus of four-stepladder program of knee osteoarthritis (2018),” Chin J Joint Surg (Electronic Edition), vol. 13, no. 1, pp. 124–130, 2019.
View at: Google Scholar
D. Huang, “Professional Committee of Orthopedics and Traumatology of Chi- nese Society of Integrative Medicine. Guidelines for diagnosis and treatment of knee osteoarthritis with integrated traditional Chinese and western medicine,” Nature Medicine J, vol. 98, no. 45, pp. 3653–3658, 2018.
View at: Google Scholar
A. A. El-Malah, M. M. Gineinah, P. K. Deb et al., “Selective COX-2 inhibitors: road from success to controversy and the quest for repurposing,” Pharmaceuticals, vol. 15, no. 7, p. 827, 2022.
View at: Publisher Site | Google Scholar
R. Dey, S. Dey, A. Samadder, A. Saxena, and S. Nandi, “Natural inhibitors against potential targets of cyclooxygenase, lipoxygenase and leukotrienes,” Comb Chem High Throughput Screen, vol. 25, no. 14, pp. 2341–2357, 2021.
View at: Publisher Site | Google Scholar
M. Sisa, M. Dvorakova, V. Temml, V. Jarosova, T. Vanek, and P. Landa, “Synthesis, inhibitory activity and in silico docking of dual COX/5-LOX inhibitors with quinone and resorcinol core,” European Journal of Medicinal Chemistry, vol. 204, Article ID 112620, 2020.
View at: Publisher Site | Google Scholar
L. Brosseau, G. A. Wells, G. P. Kenny et al., “The implementation of a community-based aerobic walking program for mild to moderate knee osteoarthritis (OA): a knowledge translation (KT) randomized controlled trial (RCT): Part I: the Uptake of the Ottawa Panel clinical practice guidelines (CPGs),” BMC Public Health, vol. 12, no. 1, p. 871, 2012.
View at: Publisher Site | Google Scholar
J. Y. Reginster, A. Neuprez, M. P. Lecart, N. Sarlet, and O. Bruyere, “Role of glucosamine in the treatment for osteoarthritis,” Rheumatology International, vol. 32, no. 10, pp. 2959–2967, 2012.
View at: Publisher Site | Google Scholar
D. Dragos, M. Gilca, L. Gaman et al., “Phytomedicine in joint disorders,” Nutrients, vol. 9, no. 1, p. 70, 2017.
View at: Publisher Site | Google Scholar
C. Wang, Y. Gao, Z. Zhang et al., “Safflower yellow alleviates osteoarthritis and prevents inflammation by inhibiting PGE2 release and regulating NF-κB/SIRT1/AMPK signaling pathways,” Phytomedicine, vol. 78, Article ID 153305, 2020.
View at: Publisher Site | Google Scholar
Q. L. Yuan, P. Wang, L. Liu et al., “Acupuncture for musculoskeletal pain: a meta-analysis and meta-regression of sham-controlled randomized clinical trials,” Scientific Reports, vol. 6, no. 1, Article ID 30675, 2016.
View at: Publisher Site | Google Scholar
L. Li, H. Liu, W. Shi et al., “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
S. Feng, T. Wang, L. Fan et al., “Exploring the potential therapeutic effect of Eucommia ulmoides-Dipsaci Radix herbal pair on osteoporosis based on network pharmacology and molecular docking technology,” RSC Advances, vol. 12, no. 4, pp. 2181–2195, 2022.
View at: Publisher Site | Google Scholar
C. Li, Research on Medication Rules and Systematic Evaluation of Osteoarthritis Based on Data Mining, Shandong University of Chinese Medicine, Shandong, China, 2020.
L. Yang, A. Hou, S. Wang et al., “A review of the botany, traditional use, phytochemistry, analytical methods, pharmacological effects, and toxicity of angelicae pubescentis Radix,” Evidence-based Complementary and Alternative Medicine, vol. 2020, Article ID 7460781, 28 pages, 2020.
View at: Publisher Site | Google Scholar
J. Zhao, G. Liang, and Y. Han, “Research progress of dipsoponin VI from Dipsophidium dipsophidium in the prevention and treatment of osteoporosis,” Chinese Journal of Osteoporosis, vol. 26, no. 05, pp. 755–759, 2020.
View at: Google Scholar
Pharmacopoeia Commission, Pharmacopoeia of the People’s Republic of China, China Medical Science Press, Beijing, China, 2015.
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
H. Y. Xu, Y. Q. Zhang, Z. M. Liu et al., “ETCM: an encyclopaedia of traditional Chinese medicine,” Nucleic Acids Research, vol. 47, no. D1, pp. D976–D982, 2019.
View at: Publisher Site | Google Scholar
Y. Wu, F. Zhang, K. Yang et al., “SymMap: an integrative database of traditional Chinese medicine enhanced by symptom mapping,” Nucleic Acids Research, vol. 47, no. 1, pp. D1110–D1117, 2019.
View at: Publisher Site | Google Scholar
Y. Wen, L. He, R. Peng et al., “A novel strategy to evaluate the quality of herbal products based on the chemical profiling, efficacy evaluation and pharmacokinetics,” Journal of Pharmaceutical and Biomedical Analysis, vol. 161, pp. 326–335, 2018.
View at: Publisher Site | Google Scholar
S. K. Yan, J. Zhao, S. Dou, P. Jiang, R. H. Liu, and W. D. Zhang, “Methodology of modernization research in traditional Chinese medicine based on systems biology and network biology,” Chinese Journal of Natural Medicines, vol. 7, no. 4, pp. 249–259, 2009.
View at: Publisher Site | Google Scholar
H. Wang, H. Wang, J. Zhang et al., “Molecular mechanism of crataegi folium and alisma rhizoma in the treatment of dyslipidemia based on network pharmacology and molecular docking,” Evidence-based Complementary and Alternative Medicine, vol. 2022, Article ID 4891370, 13 pages, 2022.
View at: Publisher Site | Google Scholar
T. Barrett, S. E. Wilhite, P. Ledoux et al., “NCBI GEO: archive for functional genomics data sets--update,” Nucleic Acids Research, vol. 41, no. D1, pp. D991–D995, 2012.
View at: Publisher Site | Google Scholar
G. Stelzer, N. Rosen, I. Plaschkes et al., “The GeneCards suite: from gene data mining to disease genome sequence analyses,” Current Protocols in Bioinformatics, vol. 54, pp. 1.30.1–1.30.33, 2016.
View at: Publisher Site | Google Scholar
Y. Zhou, Y. Zhang, X. Lian et al., “Therapeutic target database update 2022: facilitating drug discovery with enriched comparative data of targeted agents,” Nucleic Acids Research, vol. 50, no. D1, pp. D1398–D1407, 2022.
View at: Publisher Site | Google Scholar
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, pp. 1.2.1–1.2.12, 2017.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, D. Lyon et al., “STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–D613, 2019.
View at: Publisher Site | Google Scholar
J. K. Maier and P. Labute, “Assessment of fully automated antibody homology modeling protocols in molecular operating environment,” Proteins, vol. 82, no. 8, pp. 1599–1610, 2014.
View at: Publisher Site | Google Scholar
C. Zhang, C. Yang, and C. Xiao, “Studies on chemical constituents of Dipsaci Radix,” Shandong Science, vol. 35, no. 05, pp. 1–9, 2022.
View at: Google Scholar
L. Zhou and J. Zen, “Research progress on chemical composition and pharmacological activity of Radix Angelicae Biseratae,” Chinese Modern Traditional Medicine, vol. 21, no. 12, pp. 1739–1748, 2019.
View at: Google Scholar
K. Y. Hsin, S. Ghosh, and H. Kitano, “Combining machine learning systems and multiple docking simulation packages to improve docking prediction reliability for network pharmacology,” PLoS One, vol. 8, no. 12, Article ID e83922, 2013.
View at: Publisher Site | Google Scholar
S. T. Skou and E. M. Roos, “Physical therapy for patients with knee and hip osteoarthritis: supervised, active treatment is current best practice,” Clinical & Experimental Rheumatology, vol. 120, no. 5, pp. 112–117, 2019.
View at: Google Scholar
H. Zhu, A Retrospective Study of Osteoarthritis Treated by Master of Traditional Chinese Medicine Zhou Zhongying and Discussion of Clinical Experience, Nanjing University of Traditional Chinese Medicine, Nanjing, China, 2022.
H. Tian, F. Zhao, and Y. Li, “Research progress of asperosaponin VI,” Chinese Journal of Experimental Medicine, vol. 24, no. 05, pp. 226–234, 2018.
View at: Google Scholar
K. Ke, Q. Li, X. Yang et al., “Asperosaponin VI promotes bone marrow stromal cell osteogenic differentiation through the PI3K/AKT signaling pathway in an osteoporosis model,” Scientific Reports, vol. 6, no. 1, Article ID 35233, 2016.
View at: Publisher Site | Google Scholar
C. Li, Y. Gao, J. Tian, Y. Xing, H. Zhu, and J. Shen, “Long-term oral Asperosaponin VI attenuates cardiac dysfunction, myocardial fibrosis in a rat model of chronic myocardial infarction,” Food and Chemical Toxicology, vol. 50, no. 5, pp. 1432–1438, 2012.
View at: Publisher Site | Google Scholar
J. Zhang, X. Xie, and J. Li, “Research overview of osthole extract in prevention and treatment of primary osteoporosis,” Chinese Journal of Osteoporosis, vol. 26, no. 10, pp. 1546–1549, 2020.
View at: Google Scholar
Y. Ma, L. Wang, S. Zheng et al., “Osthole inhibits osteoclasts formation and bone resorption by regulating NF-κB signaling and NFATc1 activations stimulated by RANKL,” Journal of Cellular Biochemistry, vol. 120, no. 9, pp. 16052–16061, 2019.
View at: Publisher Site | Google Scholar
J. N. Choi, Y. H. Choi, J. M. Lee et al., “Anti-inflammatory effects of β-sitosterol-β-D-glucoside from Trachelospermum jasminoides (Apocynaceae) in lipopolysaccharide-stimulated RAW 264.7 murine macrophages,” Natural Product Research, vol. 26, no. 24, pp. 2340–2343, 2012.
View at: Publisher Site | Google Scholar
W. P. Chen, C. Yu, P. F. Hu, J. P. Bao, J. L. Tang, and L. D. Wu, “Stigmasterol blocks cartilage degradation in rabbit model of osteoarthritis,” Acta Biochimica Polonica, vol. 59, no. 4, pp. 537–541, 2012.
View at: Publisher Site | Google Scholar
H. Nie, Y. Zheng, R. Li et al., “Phosphorylation of FOXP3 controls regulatory T cell function and is inhibited by TNF-α in rheumatoid arthritis,” Nature Medicine, vol. 19, no. 3, pp. 322–328, 2013.
View at: Publisher Site | Google Scholar
H. Nakahara, J. Song, M. Sugimoto et al., “Anti-interleukin-6 receptor antibody therapy reduces vascular endothelial growth factor production in rheumatoid arthritis,” Arthritis & Rheumatism, vol. 48, no. 6, pp. 1521–1529, 2003.
View at: Publisher Site | Google Scholar
K. Li, M. Huang, P. Xu et al., “Microcystins-LR induced apoptosis via S-nitrosylation of GAPDH in colorectal cancer cells,” Ecotoxicology and Environmental Safety, vol. 190, Article ID 110096, 2020.
View at: Publisher Site | Google Scholar
A. Haseeb, M. Y. Ansari, and T. M. Haqqi, “Harpagoside suppresses IL-6 expression in primary human osteoarthritis chondrocytes,” Journal of Orthopaedic Research, vol. 35, no. 2, pp. 311–320, 2017.
View at: Publisher Site | Google Scholar
T. Wang and C. He, “Pro-inflammatory cytokines: the link between obesity and osteoarthritis,” Cytokine & Growth Factor Reviews, vol. 44, pp. 38–50, 2018.
View at: Publisher Site | Google Scholar
E. Shahine and A. Elhadidi, “AB0776 efficacy of glucosamine sulphate in lowering serum level of interleukin-1β in symptomatic primary knee osteoarthritis: clinical and laboratory study,” Annals of the Rheumatic Diseases, vol. 73, no. 2, pp. 1061.1–1061, 2014.
View at: Publisher Site | Google Scholar
T. Mabey and S. Honsawek, “Cytokines as biochemical markers for knee osteoarthritis,” World Journal of Orthopedics, vol. 6, no. 1, pp. 95–105, 2015.
View at: Publisher Site | Google Scholar
K. M. Tong, D. C. Shieh, C. P. Chen et al., “Leptin induces IL-8 expression via leptin receptor, IRS-1, PI3K, Akt cascade and promotion of NF-κB/p300 binding in human synovial fibroblasts,” Cellular Signalling, vol. 20, no. 8, pp. 1478–1488, 2008.
View at: Publisher Site | Google Scholar
M. Lin, Y. Lin, X. Li et al., “Warm sparse-dense wave inhibits cartilage degradation in papain-induced osteoarthritis through the mitogen-activated protein kinase signaling pathway,” Experimental and Therapeutic Medicine, vol. 14, no. 4, pp. 3674–3680, 2017.
View at: Publisher Site | Google Scholar
N. Saetan, S. Honsawek, A. Tanavalee et al., “Relationship of plasma and synovial fluid vascular endothelial growth factor with radiographic severity in primary knee osteoarthritis,” International Orthopaedics, vol. 38, no. 5, pp. 1099–1104, 2014.
View at: Publisher Site | Google Scholar
E. Jimi, F. Huang, and C. Nakatomi, “NF-κB signaling regulates physiological and pathological chondrogenesis,” International Journal of Molecular Sciences, vol. 20, no. 24, p. 6275, 2019.
View at: Publisher Site | Google Scholar
M. Ostojic, A. Zevrnja, K. Vukojevic, and V. Soljic, “Immunofluorescence analysis of NF-kB and iNOS expression in different cell populations during early and advanced knee osteoarthritis,” International Journal of Molecular Sciences, vol. 22, no. 12, p. 6461, 2021.
View at: Publisher Site | Google Scholar
M. C. Choi, J. Jo, J. Park, H. K. Kang, and Y. Park, “NF-κB signaling pathways in osteoarthritic cartilage destruction,” Cells, vol. 8, no. 7, p. 734, 2019.
View at: Publisher Site | Google Scholar
T. Tanaka, M. Narazaki, and T. Kishimoto, “IL-6 in inflammation, immunity, and disease,” Cold Spring Harbor Perspectives in Biology, vol. 6, no. 10, Article ID a016295, 2014.
View at: Publisher Site | Google Scholar
B. B. Aggarwal, S. C. Gupta, and J. H. Kim, “Historical perspectives on tumor necrosis factor and its superfamily: 25 years later, a golden journey,” Blood, vol. 119, no. 3, pp. 651–665, 2012.
View at: Publisher Site | Google Scholar
P. Wang, P. P. Guan, C. Guo, F. Zhu, K. Konstantopoulos, and Z. Y. Wang, “Fluid shear stress-induced osteoarthritis: roles of cyclooxygenase-2 and its metabolic products in inducing the expression of proinflammatory cytokines and matrix metalloproteinases,” The FASEB Journal, vol. 27, no. 12, pp. 4664–4677, 2013.
View at: Publisher Site | Google Scholar
Y. Quan, B. Zhou, Y. Wang et al., “Association between IL17 polymorphisms and risk of cervical cancer in Chinese women,” Clinical and Developmental Immunology, vol. 2012, Article ID 258293, 6 pages, 2012.
View at: Publisher Site | Google Scholar
Y. Bai, S. Gao, Y. Liu, S. Jin, H. Zhang, and K. Su, “Correlation between Interleukin-17 gene polymorphism and osteoarthritis susceptibility in Han Chinese population,” BMC Medical Genetics, vol. 20, no. 1, p. 20, 2019.
View at: Publisher Site | Google Scholar
C. Maes, G. Carmeliet, and E. Schipani, “Hypoxia-driven pathways in bone development, regeneration and disease,” Nature Reviews Rheumatology, vol. 8, no. 6, pp. 358–366, 2012.
View at: Publisher Site | Google Scholar
J. Fernández-Torres, Y. Zamudio-Cuevas, G. A. Martínez-Nava, and A. G. López-Reyes, “Hypoxia-Inducible Factors (HIFs) in the articular cartilage: a systematic review,” European Review for Medical and Pharmacological Sciences, vol. 21, no. 12, pp. 2800–2810, 2017.
View at: Google Scholar
W. Bouaziz, J. Sigaux, D. Modrowski et al., “Interaction of HIF1α and β-catenin inhibits matrix metalloproteinase 13 expression and prevents cartilage damage in mice,” Proceedings of the National Academy of Sciences of the U S A, vol. 113, no. 19, pp. 5453–5458, 2016.
View at: Publisher Site | Google Scholar
J. H. Rosenberg, V. Rai, M. F. Dilisio, and D. K. Agrawal, “Damage-associated molecular patterns in the pathogenesis of osteoarthritis: potentially novel therapeutic targets,” Molecular and Cellular Biochemistry, vol. 434, no. 1-2, pp. 171–179, 2017.
View at: Publisher Site | Google Scholar
A. Suzuki, A. Yabu, and H. Nakamura, “Advanced glycation end products in musculoskeletal system and disorders,” Methods, vol. 203, pp. 179–186, 2022.
View at: Publisher Site | Google Scholar
J. Xie, J. D. Méndez, V. Méndez-Valenzuela, and M. M. Aguilar-Hernández, “Cellular signalling of the receptor for advanced glycation end products (RAGE),” Cellular Signalling, vol. 25, no. 11, pp. 2185–2197, 2013.
View at: Publisher Site | Google Scholar
Z. Rasheed and T. M. Haqqi, “Endoplasmic reticulum stress induces the expression of COX-2 through activation of eIF2α, p38-MAPK and NF-κB in advanced glycation end products stimulated human chondrocytes,” Biochimica et Biophysica Acta (BBA) - Molecular Cell Research, vol. 1823, no. 12, pp. 2179–2189, 2012.
View at: Publisher Site | 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
X. Q. Hao, Y. Q. Kou, J. W. Xie Xj, J. B. Lv, and J. Su, “Network pharmacology-based study of the anti-oxidative mechanism of san miao wan in treatment of arthritis,” World J Tradit Chin Med, vol. 8, pp. 100–109, 2022.
View at: Google Scholar
H. Wang, J. Liu, J. He et al., “Potential mechanisms underlying the therapeutic roles of sinisan formula in depression: based on network pharmacology and molecular docking study,” Frontiers in Psychiatry, vol. 13, Article ID 1063489, 2022.
View at: Publisher Site | Google Scholar
S. Li, M. Niu, K. Yang, B. Zhang, H. Xu, and M. Yang, “Network pharmacology evaluation method guidance-draft,” World J Tradit Chin Med, vol. 7, pp. 146–154, 2021.
View at: Google Scholar

Copyright

Copyright © 2023 Yujiang Xi 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

505

Downloads

445

Citations