Skip to main content

Study on the potential active components and molecular mechanism of Xiao Huoluo Pills in the treatment of cartilage degeneration of knee osteoarthritis based on bioinformatics analysis and molecular docking technology

Abstract

Background

Knee osteoarthritis is a common joint degenerative disease. Xiao Huoluo Pills (XHLP) has been used to treat degenerative diseases such as osteoarthritis and hyperosteogeny. However, XHLP’s specific effective ingredients and mechanism of action against osteoarthritis have not been explored. Therefore, bioinformatics technology and molecular docking technology are employed in this study to explore the molecular basis and mechanism of XHLP in the treatment of knee osteoarthritis.

Methods

Public databases (TCMSP, Batman-TCM, HERB, DrugBank, and UniProt) are used to find the effective active components and corresponding target proteins of XHLP (screening conditions: OB > 30%, DL ≥ 0.18). Differentially expressed genes related to cartilage lesions of knee osteoarthritis are obtained based on the GEO database (screening conditions: adjust P value < 0.01, |log2 FC|≥1.0). The Venn package in R language and the BisoGenet plug-in in Cytoscape are adopted to predict the potential molecules of XHLP in the treatment of knee osteoarthritis. The XHLP-active component-target interaction network and the XHLP-knee osteoarthritis-target protein core network are constructed using Cytoscape software. Besides, GO/KEGG enrichment analysis on core genes is performed using the Bioconductor package and clusterProfiler package in the R language to explain the biological functions and signal pathways of the core proteins. Finally, molecular docking is performed through software such as Vina, LeDock, Discovery Studio 2016, PyMOL, AutoDockTools 1.5.6, so as to verify the binding ability between the active components of the drug and the core target protein.

Results

XHLP has been screened out of 71 potentially effective active compounds for the treatment of OA, mainly including quercetin, Stigmasterol, beta-sitosterol, Izoteolin, and ellagic acid. Knee osteoarthritis cartilage lesion sequencing data (GSE114007) was screened out of 1672 differentially expressed genes, including 913 upregulated genes and 759 downregulated genes, displayed as heat maps and volcano maps. Besides, 33 core target proteins are calculated by Venn data package in R and BisoGenet plug-in in Cytoscape. The enrichment analysis on these target genes revealed that the core target genes are mainly involved in biological processes such as response to oxygen levels, mechanical stimulus, vitamin, drug, and regulation of smooth muscle cell proliferation. These core target genes are involved in signaling pathways related to cartilage degeneration of knee osteoarthritis such as TNF signaling pathway and PI3K-Akt signaling pathway. Finally, the molecular docking verification demonstrates that some active components of the drug have good molecular docking and binding ability with the core target protein, further confirming that XHLP has the effect of inhibiting cartilage degeneration in knee osteoarthritis.

Conclusions

In this study, based on the research foundation of bioinformatics and molecular docking technology, the active components and core target molecules of XHLP for the treatment of cartilage degeneration of knee osteoarthritis are screened out, and the potential mechanism of XHLP inhibiting cartilage degeneration of knee osteoarthritis is deeply explored. The results provide theoretical basis and new treatment plan for XHLP in the treatment of knee osteoarthritis.

Introduction

The characteristic case changes of knee osteoarthritis (KOA) include articular cartilage degeneration, subchondral bone reactivity, osteophyte formation at the joint edge, synovial lesions, ligament laxity, and contracture [1, 2]. With the aging of the population, the incidence of knee osteoarthritis is increasing year by year. Globally, osteoarthritis affects approximately 250 million people [3]. Chronic pain and disability associated with osteoarthritis (OA) can lead to anxiety, depression, and suicidal emotions [4]. However, long-term use of non-steroidal anti-inflammatory drugs in OA patients is associated with gastrointestinal and cardiovascular side effects [5, 6]. Advanced KOA is mainly treated by joint replacement surgery, which is expensive and has many postoperative complications [7, 8]. Therefore, it is necessary to further explore the occurrence and development mechanism of KOA, and provide new theoretical basis and diagnosis and treatment guidelines for the precise treatment of KOA.

Traditional Chinese medicine has been used to treat KOA for many years. Starting from the local application of knee joint medicine and the whole body, it can effectively relieve pain symptoms, improve joint function, and protect joint structure [9, 10]. XHLP comes from the classic traditional Chinese medicine “Taiping Huimin He Jifang Prescription,” composed of six herbs: Chuanwu, Caowu, Dilong, Araceae, Frankincense, and Myrrh [11]. At present, it is mainly used to treat diseases such as rheumatoid arthritis, hyperosteosis, osteoarthritis, and stroke hemiplegia [12]. Modern pharmacological studies have confirmed that XHLP has the effects of anti-inflammatory, anti-oxidant damage, analgesia, blood vessel expansion, inhibiting platelet aggregation, and improving immunity [13]. Besides, the diester-type diterpene alkaloids of the active extracts of Radix Aconiti and Radix Aconitum play anti-inflammatory effects by inhibiting the process of prostaglandin metabolism [14] and chemotaxis of leukocytes mediated by chemokines [15]. Seo et al. [16] revealed that the antibacterial peptide lumbricusin of earthworm extract can significantly reduce the inflammatory reaction in vivo and in vitro caused by LPS such as COX-2, IL-6, and TNF-α [17]. By inhibiting the activation of neutrophil product NF-κB, frankincense extract can downregulate the expression of inflammatory factors TNF-α, IL-1, and INF-γ, alleviating inflammation [18, 19]. Meanwhile, it can inhibit the activity of collagenase MMP-1 and alleviate cartilage degeneration [20]. Therefore, XHLP is expected to have a strong anti-osteoarthritis effect. However, the potential pharmacological mechanism of XHLP and its interaction with osteoarthritis-related targets and pathways still need to be further investigated.

The traditional Chinese herbal medicine treatment KOA has passed the test of long-term medical practice. However, the “multi-component” and “multi-functional” of Chinese medicine has not been recognized by the world because its specific mechanism of action remains unclear [21]. Network pharmacology can reveal the complex network relationship between drugs, targets, and diseases through high-throughput screening, network visualization, and network topology analysis, and accurately predict and analyze the mechanism of action of traditional Chinese medicine compounds [22]. It has the characteristics of “multi-component, multi-target, multi-pathway,” turning the complex into simple. Additionally, potential differentially expressed genes are analyzed and discovered by using the genome sequencing data of KOA patients [23, 24], providing strong evidence for the biological targets and potential mechanisms of XHLP in the treatment of KOA. We speculate that the XHLP treatment of KOA through the characteristics of “multi-component, multi-target, multi-pathway.” This study aims to use bioinformatics analysis and molecular docking technology to clarify and verify the effective active components, targets, and potential mechanisms of XHLP in knee osteoarthritis (Fig. 1).

Fig. 1
figure 1

Schematic diagram for identifying the mechanism of XHLP anti-osteoarthritis by bioinformatics analysis

Methods

XHLP’s effective active components screening and identification

XHLP is composed of six traditional Chinese medicines: Chuanwu, Caowu, Dilong, Araceae, Frankincense, and Myrrh, with the composition ratio of 1.2:1.2:1.2:1.2:1.0:1.0. XHLP’s constituent medicines “Chuanwu,” “Caowu,” “Dilong,” “Araceae,” “Frankincense,” and “Myrrh” were entered into the TCM pharmacology system database [25] (TCMSP: traditional Chinese medicine system pharmacology, https://tcmspw.com/tcmsp.php) in turn for search. Oral bioavailability (OB) refers to the speed and degree of absorption of oral drugs into the circulation of the body (including the bioavailability and utilization rate of the drug). Drug likeness (DL) is a concept based on the physical and chemical properties and molecular structure of existing drugs, used to evaluate the similarity between a compound and a known drug. According to the common screening criteria of the TCMSP database and previous literature results [26, 27], the active components with OB > 30% and DL ≥ 0.18 were selected for subsequent analysis and research.

Target protein screening and drug-active component-target network construction

The TCMSP database is used to search for the target of the active components of the drugs obtained by screening, and all the target protein names obtained from the screening of the TCMSP database are imported into the DrugBank database (https://www.drugbank.ca/) and the UniProt database (https://www.uniprot.org). The species is set as Homo sapiens. Comparison and gene name correction are performed on the obtained target protein, converting into the target gene name and deleting the target protein without corresponding gene name. The obtained active component of the medicine and the active component establish a corresponding relationship with the target. The drug-active component-target interaction network model is obtained by Cytoscape 3.7.2 software [28]. The “Networkanalyze” function in Cytoscape 3.7.2 software [29] is employed to perform network topology analysis, screening out core active components and potential target genes.

Obtaining gene expression profile of human osteoarthritis cartilage tissue

The high-throughput sequencing data of cartilage tissue genes related to osteoarthritis are downloaded in the GEO (gene expression omnibus) database (https://www.ncbi.nlm.nih.gov/geo/). The screening conditions are (1) osteoarthritis, (2) human, and (3) cartilage tissue. The whole-genome sequencing data GSE114007 that meets the screening criteria are obtained, including 18 normal knee cartilage tissue samples from K-L I patients and 20 knee cartilage tissue samples from K-L IV KOA patients. All the above sequencing data sets are downloaded based on the GEO database of GPL11154 and GPL18573 platforms. The general clinical information of patients is presented in Table 1.

Table 1 Characteristics of study subjects

Acquisition and analysis of differential genes in cartilage tissue of osteoarthritis

The original sequencing data of GSE114007 are downloaded from the GEO database. The FastQC software is adopted to perform quality assessment and quality control processing on the original data. Besides, the Hisat2 software is used to select a reference genome to compare and analyze the clean data obtained after quality control processing. Then, the number of gene reads is quantitatively analyzed using the HTSeq Python package to obtain the count file. Next, the counts data obtained by sorting are filtered. Filter conditions are (1) variance filtering and (2) low abundance filtering. Besides, the data are standardized using the Log2-counts per million method. The DESeq2 package in R language is employed to perform difference analysis and calculate the difference genes between the two groups. Set P < 0.01. adjust P value < 0.01 indicates that the change range is greater than or equal to twice |log2 FC|≥1.0 as the criterion for screening differential genes. Among them, log2FC≥1.0 and log2FC≤-1.0 represent upregulation and downregulation of gene expression, respectively. Finally, the differentially expressed genes in cartilage tissue of the osteoarthritis group and the healthy control group were derived, namely KOA differentially expressed genes (DEGs). The online analysis website ClustVi s[30] (https://biit.cs.ut.ee/clustvis/) is adopted to draw the heat map and cluster analysis of the differential genes obtained by screening. Then, the adjust P value in the processed data is converted to -log10. According to log2 FC, -log10 (adjust P value) is divided into upregulated genome, downregulated genome, no statistical difference, and non-expressed difference genome. The processed data are imported into GraphPad Prism 7 to draw a volcano map.

Intersection and interactive network construction of active component of medicine and osteoarthritis target

The Venn diagram of the common target genes in the target of the active component of the drug and the target of the disease is displayed through the online Venn diagram production website InteractiVenn (http://www.interactivenn.net/) [31]. The obtained common target genes are imported into the online analysis website STRING [32] to construct protein-protein interaction (PPI). Besides, the BisoGenet plug-in in the Cytoscape 3.7.2 software [33] is used to construct the target PPI of the drug active component and the PPI network of the osteoarthritis differential gene. The “Networkanalyze” function in Cytoscape 3.7.2 software is used for network topology analysis. The PPI and KOA differential gene PPI networks of the active pharmaceutical ingredient are combined to obtain the PPI network shared by the active pharmaceutical ingredient PPI and the KOA differential gene PPI network. The CytoNCA plug-in in the Cytoscape 3.7.2 software is adopted to calculate the multi-center network topology for the PPI shared by the target of the active component of the drug and the KOA differential gene. Then, screening is conducted according to the median value of Degree, BetweennessCentrality (BC), ClosenessCentrality (CC), LAC, and NeighborhoodConnectivity (NC) values.

Enrichment analysis of target gene by GO/KEGG

To further explore the potential mechanism of XHLP active components on KOA and the function of common target genes, we use the Bioconductor package and clusterProfiler package [34] in the R language to perform gene ontology (GO) and KEGG (Kyoto Encyclopedia of Genes and Genomes) signal pathway analysis on the target genes obtained in the previous step. The GO analysis mainly includes cellular component (CC), molecular function (MF), and biological process (BP) of potential target genes. Meanwhile, the Pathview package is used to draw the corresponding signal path diagram.

Molecular docking verification of the binding of the active component of the drug to the target protein

The molecular docking verification of XHLP and potential targets is conducted to explore the interaction between them. The structural formula of the active components of XHLP is downloaded from the TCMSP database (https://tcmspw.com/tcmsp.php). The Chem3D software is used to make 3D structures of effective active components. Then, the 3D structure of the core protein gene is downloaded from the PDB database (http://www.rcsb.org/). The PyMOL software is employed to perform operations such as dewatering and phosphate removal of proteins. Besides, the AutoDockTools 1.5.6 software [35] is adopted to convert the pdb format of active components and core protein gene files into pdbqt format and search for active pockets. Finally, the Vina script [36] and LeDock software [37] are applied to calculate the molecular binding energy and display the results of molecular docking. Simultaneously, Discovery Studio 2016 is used to find the docking site and calculate the LibDockScore [38, 39]. Next, the output molecular docking results are imported into PyMOL software for molecular docking verification display. Binding energy < 0 indicates that the ligand and the receptor can bind spontaneously. The binding energy ≤− 5.0 kcal mol−1 is selected as the screening basis to evaluate the reliability of biological information analysis and prediction [40].

Results

Screening of effective active components of XHLP

A total of 568 active components of XHLP were retrieved through TCMSP database and Batman-TCM database, including 12 species of Radix Aconiti, 27 species of Radix Aquilinum, 123 species of Araceae, 127 species of Frankincense, and 276 species of Myrrh. However, Earth Dragon failed to retrieve the active components in the TCMSP database and Batman-TCM database. Therefore, three active components of Earthworm were obtained through literature search and HERB (A high-throughput experiment- and reference-guided database of traditional Chinese medicine, http://herb.ac.cn/) identification by the herbal group. Under the screening conditions of OB ≥ 30% and DL ≥ 0.18, 71 effective active compounds were selected from the active components of XHLP, including 3 species of Radix Aconiti and 8 species of Radix Aconiti, 3 species of Earth Dragon, 7 species of Araceae, 8 species of Frankincense, and 45 species of Myrrh; besides, there is 1 common active component of Radix Aconiti and Radix Aquilinum, and 2 common active components of Araceae and Myrrh (Table 2).

Table 2 XHLP effective active component information

Construction and analysis of XHLP-active component-target network

The 71 active components in XHLP were entered into the TCMSP database and Batman-TCM database to find their corresponding target proteins. The 206 target proteins obtained were imported into the DrugBank database and the UniProt database for comparison and correction, and 170 standard gene names were output. Next, the XHLP composition and effective active components and the corresponding relationship between the active components and the target was obtained. The XHLP-active component-target network was constructed using Cytoscape 3.7.2 software. The network has 248 nodes, including 6 drug nodes, 71 active compound nodes, and 170 target nodes, with a total of 506 edges (Fig. 2). Each edge in the network represents the active component contained in the drug and the interaction between the active component and the target gene. The “Networkanalyze” function in the Cytoscape 3.7.2 software is employed to perform topology analysis and calculation on the network. Among them, Degree value and BetweennessCentrality value (BC) are the key parameters to measure the nodes in the network. The active components in this network have an average relationship with 6 targets, and each target has an interconnection relationship with an average of 2.5 active components. Therefore, the active components in XHLP may be related to multiple targets while one target may be related to multiple active components. It can be revealed by analyzing the Degree value and BC value of the active component node and target gene node in the network that the top 5 active molecules are MOY27-quercetin, B2-Stigmasterol, B1-beta-sitosterol, CAW1-Izoteolin, and MOY2-ellagic acid, which can be connected to 136, 53, 53, 18, and 17 target proteins, respectively; the top six targets are PGR, NCOA2, NR3C2, PTGS2, PTGS1, and RXRA, which can interact with 25, 24, 18, 18, 12, and 12 active components, respectively (Table 3).

Fig. 2
figure 2

“XHLP-Active Component-Target” Network: the purple circles in the network represent the drugs in XHLP, the orange hexagons represent the active components, and the green diamonds represent the target genes

Table 3 Degree and BC values of the top 8 effective active components and targets of XHLP

Acquisition and analysis of differentially expressed genes in KOA cartilage tissue

Our research is based on 20 cases of cartilage tissue samples from KOA patients (8 male cases, 12 female cases, aged 66.20 ± 7.16, and K-L IV patients accounted for 20 cases) and 18 cases of normal cartilage tissue samples (aged 36.61 ± 13.08, 13 male cases and 5 female cases) from the sequencing data set GSE114007. The general clinical characteristics of the patients are presented in Table 1. The filter conditions are variance filter = 15.0 and low abundance = 4.0; the standardization method is Log2-counts per millio. Box diagram and PCA diagram of the data filtering and before and after standardization are illustrated in Fig. 3A, B. In the GSE114007 data set, 1672 differentially expressed genes were screened, including 913 upregulated genes and 759 downregulated genes. A heat map was drawn for the DEGs selected in GSE114007 and the top 50 DEGs with the most significant differences in the adjust P Value selection (Fig. 3C, D). Among them, red and green represent upregulation and downregulation of gene expression. Then, -log10 conversion was performed on the adjust P value of the gene after the differential analysis process. According to log2 FC, -log10 (adjust P value) is divided into upregulated genome, downregulated genome, and no statistically different genome). The results were imported into GraphPad Prism 7.0 to draw a volcano map (Fig. 3E).

Fig. 3
figure 3

Bioinformatics analysis of GSE114007. A Box plots of samples before and after normalization, representing the counts distribution of each sample. B PCA plot before and after normalization; each point represents a sample; the farther the distance between the two samples, the greater the difference between the two samples. C Heat map of total 1672 DEGs in GSE114007. D Heat map of top 50 DEGs in GSE114007. Tissue samples are presented as columns; individual genes are presented as rows. Red indicates upregulated genes and green indicates downregulated genes in patients with KOA. The patients in the top rows in pink indicates columns were normals; the patients in the top rows in blue indicates columns were patients with KOA. E The volcano plot of total 1672 DEGs in GSE114007. The vertical blue lines correspond to Log2FC up and down and the horizontal orange line represents an adjusted P value of 0.01. The red point represents the differentially upregulated genes with statistical significance. The green point represents the differentially downregulated genes with statistical significance.

XHLP effective active component-KOA intersection target analysis and mutual PPI network construction

The 170 drug active component targets obtained above and the 1672 differentially expressed genes obtained from the differential analysis of the sequencing data set GSE114007 were introduced into the online Venn diagram production website InteractiVenn. Then, a total of 33 potential targets for XHLP treatment of KOA cartilage changes were obtained by matching mapping. The Venn diagram of the effective active component-KOA intersection target is exhibited in Fig. 4A. Besides, 33 potential target genes were imported into STRING online analysis website (https://string-db.org/); hide unconnected targets were set, and “medium confidence > 0.400” in the lowest interaction score was determined; the result data of protein-protein interaction were exported (Fig. 4B, C). Next, the PPI network of potential target genes was obtained by Cytoscape 3.7.2 software, and the cytoHubba plug-in was adopted to screen Hub genes [41]. Two different algorithms, Degree and Betweenness, were employed to the top 6 potential core target genes: VEGFA, CCND1, MYC, JUN, MMP9, and MMP2 (Fig. 4D–F). Meanwhile, the BisoGenet plug-in in the Cytoscape 3.7.2 software was applied to construct the PPI network of the drug active component target and the KOA differential gene PPI network (Fig. 5A, B). The two PPI networks were merged and mapped to an intersection network (Fig. 5C). The network topology of the intersection network was analyzed and calculated using CytoNCA to obtain parameter values such as Degree, Betweenness, BetweennessCentrality (BC), ClosenessCentrality (CC), NeighborhoodConnectivity (NC), and LAC. The first screening threshold was Degree > 66 (that is, twice the median value of Degree). The result revealed that there are 1170 nodes and 52,739 edges in the secondary network graph (Fig. 5D). The second screening thresholds were Degree > 109, Betweenness > 460.887, BC > 0.00043, CC > 0.4421, NC > 146.685, and LAC > 16.999. The results indicated that there are 144 nodes and 1723 edges in the three-level network graph, which also included genes such as VEGFA, CCND1, MYC, JUN, MMP9, and MMP2 (Fig. 5E).

Fig. 4
figure 4

“XHLP-KOA-Common Targets” network analysis. A Venn diagram of KOA-related targets of XHLP into the active components. B PPI network generated on XHLP-KOA mapped targets. C Top 30 candidate targets extracted from PPI according to the node degree. D PPI network diagram of KOA joint targets treated with XHLP. E PPI network of the top 10 hub proteins of co-targets based on degree centrality. F PPI network of the top 10 hub proteins of co-targets based on Betweenness centrality

Fig. 5
figure 5

Process of topological screening for the PPI network. A PPI network of XHLP targets. B PPI network of KOA differentially expressed genes. C PPI network of key targets of merging mapping. D PPI network of key genes with a screening threshold of Degree>66 (2 times the median value of Degree). E PPI network of key genes with a second screening threshold of Degree>109, Betweenness>460.887, BC>0.00043, CC>0.4421, NC>146.685和LAC>16.999 (the median value of all the above)

Enrichment analysis of intersection target by GO/KEGG

The Bioconductor package and clusterProfiler package in R language were used to perform GO and KEGG pathway enrichment analysis on 33 XHLP active components-KOA intersection targets. The GO analysis of 33 potential target genes demonstrated that the biological process (BP) mainly focuses on response to oxygen levels, mechanical stimulus, vitamin, drug, and regulation of smooth muscle cell proliferation; cellular component (CC) is mainly concentrated in collagen-containing extracellular matrix, endoplasmic reticulum lumen, and fibrillar collagen trimer; molecular function (MF) is mainly manifested in activating transcription factor binding, growth factor binding, and core promoter sequence-specific DNA binding (Table 4, Fig. 6A–C). The enrichment analysis indicated that the KEGG pathway mainly focuses on AGE-RAGE signaling pathway in diabetic complications, Relaxin signaling pathway, TNF signaling pathway, PI3K-Akt signaling pathway, fluid shear stress, and atherosclerosis signaling pathways. The pathview package was adopted to display the signal pathway diagram related to cartilage (Fig 6D–G).

Table 4 Annotation of intersection target gene GO/KEGG enrichment analysis
Fig. 6
figure 6

GO and KEGG enrichment analyses of 33 candidate targets. A Three types of GO terms for target genes and gene products. B The Bubble chart of the top 20 terms extracted according to the p.adjust value based on GO enrichment analysis. C The Bar chart of the top 20 terms extracted according to the p.adjust value based on GO enrichment analysis. D The Bubble chart of the top 20 pathways extracted according to the p.adjust value based on KEGG enrichment analysis. E The Bar graph of the top 20 pathways extracted according to the p.adjust value based on KEGG enrichment analysis. F Top 5 significantly enriched BP, CC, and MF categories based on gene ontology, and the Top 8 significantly enriched BP, CC, and MF categories. G TNF signaling pathway map of KOA targets treated with XHLP. H PI3K-Akt signaling pathway map of KOA targets treated with XHLP. (the red nodes represent the most potential targets, arrows represent the activation effect, T arrows represent the inhibition effect, and segments show the activation effect or inhibition effect)

Molecular docking verification

The top 5 active molecules in the 2.2 results (quercetin, Stigmasterol, beta-sitosterol, Izoteolin, and ellagic acid) and the top 6 potential core target genes in the 2.4 results (VEGFA, CCND1, MYC, JUN, MMP9, and MMP2) were selected for molecular docking verification. The Chem3D software was employed to draw the corresponding 3D structure according to the structural formulas of the 5 effective active components, and output in mol*2 format. The 3D structures of 6 core proteins were downloaded from the PDB database and output in pdb format. The AutoDockTools 1.5.6 software was used to convert active components and core proteins into pdbqt format to find active pockets, that is, the ligand is combined with one or more amino acid residues to form H bonds, H-π bonds, H-π bonds, or π-π bonds and other active sites. Vina script, LeDock, and Discovery Studio 2016 software were operated to calculate the binding energy of ligand and receptor. The results suggested that quercetin can form a stable docking model with CCND1, JUN, MMP9, and MMP2, respectively; beta-sitosterol can dock with JUN and MMP2 protein ligands; izoteolin and ellagic acid can form stable docking with JUN, MMP9, and MMP2 protein ligands (Table 5). When the receptor-ligand binding energy calculated by Vina and LeDock was ≤− 5.0 kcal mol−1 and DS can find the docking site, the results of the active component receptor output by the LeDock software were imported into Pymol, and the 3D molecular docking display was performed with the protein ligand (Fig. 7). Simultaneously, the Discovery Studio 2016 software was adopted to display the results of molecular docking that meet the requirements in two dimensions (Fig. 8).

Table 5 Comparison of binding energy of molecular docking Vina, LeDock, and DS
Fig. 7
figure 7

3D Molecular docking model of the 5 key active ingredients with the 6 hub-targets. A The molecular docking of quercetin with CCND1(2W96). B The molecular docking of quercetin with JUN(3V3V). C The molecular docking of quercetin with MMP9 (6ESM). D The molecular docking of quercetin with MMP2(4WK7). E The molecular docking of beta-sitosterol with JUN(3V3V). F The molecular docking of beta-sitosterol with MMP2(4WK7). G The molecular docking of Izoteolin with JUN(3V3V). H The molecular docking of Izoteolin with MMP9(6ESM). I The molecular docking of Izoteolin with MMP2(4WK7). J The molecular docking of ellagic acid with JUN(3V3V). K The molecular docking of ellagic acid with MMP9(6ESM). L The molecular docking of ellagic acid with MMP2(4WK7)

Fig. 8
figure 8

2D virtual molecular docking of bioactive ingredients from XHLP for knee osteoarthritis targets. A The molecular docking of quercetin with CCND1(2W96). B The molecular docking of quercetin with JUN(3V3V). C The molecular docking of quercetin with MMP9(6ESM). D The molecular docking of quercetin with MMP2(4WK7). E The molecular docking of beta-sitosterol with JUN(3V3V). F The molecular docking of beta-sitosterol with MMP2(4WK7). G The molecular docking of Izoteolin with JUN(3V3V). H The molecular docking of Izoteolin with MMP9(6ESM). I The molecular docking of Izoteolin with MMP2(4WK7). J The molecular docking of ellagic acid with JUN(3V3V). K The molecular docking of ellagic acid with MMP9(6ESM). L The molecular docking of ellagic acid with MMP2(4WK7)

Discussion

Bioinformatics was first proposed by Dr. Li Huaan in 1978. People gradually realized the importance of combining computer science with biology. With the continuous surge of experimental data and biological information (about doubling every 15 months), information management, and analysis has become an urgent problem to be solved. Therefore, data mining can help us find valuable information from massive amounts of data and convert it into meaningful knowledge. At present, life science research is transforming from experimental analysis combined with data accumulation to data analysis to guide experimental verification [42]. Network pharmacology and bioinformatics analysis can be based on massive experimental data and clinical trial results, contributing to exploring the potential mechanism of TCM treatment and prevention of diseases, especially for the research on the complex chemical components of Chinese herbal medicine and the multi-target and multi-mechanism of diseases [43]. Studies have discovered that some of the extracts in the XHLP composition of drugs play an essential role in the treatment of KOA. However, the specific treatment mechanism has not been thoroughly studied. Therefore, based on bioinformatics analysis and mining, this study explored the potential effective active components, molecular targets, and mechanism of action of XHLP in the treatment of KOA cartilage degeneration, and further used molecular docking technology to verify the binding ability of the effective active components in XHLP with molecular targets. It fully verified the molecular basis of XHLP treatment of KOA and provided a theoretical basis for clinical treatment of KOA. Besides, bioinformatics analysis demonstrated that the potential mechanism of XHLP in the treatment of KOA cartilage degeneration has the features of “multi-component, multi-target, multi-pathway.”

In this study, the XHLP-active component-target network was constructed. The network topology analysis revealed that the potential active components of XHLP treatment KOA mainly include quercetin, Stigmasterol, beta-sitosterol, Izoteolin, and ellagic acid. Hu et al. [44] discovered that quercetin can inhibit the inflammation and apoptosis of chondrocytes and provide cartilage cells with an environment promoting cartilage formation to enhance cartilage repair in the OA environment. Gabay [45] observed that Stigmasterol is a plant sterol that binds to chondrocyte membranes, can reduce the expression of MMP related to cartilage degradation and the pro-inflammatory mediator PGE2 by inhibiting the NF-κB pathway, and has the properties of anti-inflammatory and inhibiting cartilage matrix catabolism. KIM et al. [46] proposed that in high-fat diet-induced intestinal inflammation, beta-sitosterol can inactivate NF-κB by interfering with the binding between LPS and TLR4 to reduce the production of pro-inflammatory cytokines. Besides, Lin et al. [47] demonstrated that ellagic acid can inhibit inflammation caused by IL-1β stimulation; reduce the expression of MMP-13, ADAMTS5, COX-2, and INOS in human OA chondrocytes; and weaken IL-1β-induced inflammation by inhibiting the NF-κB pathway. There is no related research on Izoteolin in KOA. Therefore, the effective active components in XHLP may inhibit cartilage cell apoptosis, reduce cartilage matrix catabolism, and impede the NF-κB pathway to lower the production of pro-inflammatory cytokines, playing a synergistic effect in slowing down the progression of KOA cartilage degeneration.

Through the mutual mapping of the XHLP target and the differentially expressed genes of KOA, 33 potential molecular targets were obtained. Furthermore, Cytohubba, BisoGenet, and CytoNCA plug-in in Cytoscape software were used to obtain 6 potential therapeutic targets for XHLP treatment of KOA, namely VEGFA, CCND1, MYC, JUN, MMP9, and MMP2. A large number of studies have revealed that the expression of VEGFA is elevated in patients with advanced OA [48,49,50]. VEGFA can upregulate the expression of MMP-1, MMP-3, and MMP-13 in chondrocytes and downregulate the expression of key extracellular matrix components (including aggrecan and type II collagen) to accelerate the process of KOA cartilage degeneration. Zan et al. [51] discovered that CCND1 gene silencing can enhance the inhibitory effect on IL-1β-induced proliferation of OA chondrocytes. Zou et al. [52] put forward that silencing the C-myc gene can reduce the effect of IL-1β on the cycle progression and apoptosis of chondrocytes in OA rats. Rhee et al. [53] found that IL-1β activates the expression of JUN in articular chondrocytes, inhibiting the expression of SOX9 and type II collagen in chondrocytes and the anabolic reaction of the extracellular matrix, accelerating the damage of arthritic chondrocytes. Studies have suggested that MMP-9 exacerbates local inflammation by increasing the expression of pro-inflammatory mediators such as TNF-a and IL-1, leading to accelerated degeneration of cartilage and overgrowth of subchondral bone [54, 55]. Besides, MMP-2 can effectively lyse and inactivate monocyte chemotactic protein 3 (MCP-3). This not only prevents the inflammatory response from starting but also removes inflammatory factors such as TNF-a and IL-17 in the body, preventing the further development of inflammation [56]. Therefore, the above targets may all play a role in the development of KOA cartilage degeneration and are potential targets for XHLP to treat KOA.

The enrichment analysis of target genes GO/KEGG indicated that “oxygen levels,” “mechanical stimulus,” “collagen-containing extracellular matrix,” “TNF signaling pathway,” “PI3K-Akt signaling pathway,” and “Fluid shear stress” may be a potential way of XHLP treatment of KOA. Many studies [57,58,59] have revealed that the mechanical stress environment of joints is a significant factor affecting the activity of cartilage cells in the body, and accelerates the catabolism of osteoarthritis cartilage by interacting with pro-inflammatory factors. TNF-α can inhibit the synthesis of matrix macromolecules; stimulate the synthesis and secretion of cartilage degrading enzymes, prostaglandin (PGE2); and promote the synthesis of inflammatory factors, leading to cartilage destruction, increased inflammation, and finally degeneration of KOA cartilage [60, 61]. Moreover, the activation of PI3K/AKT signaling pathway can inhibit chondrocyte apoptosis and promote the proliferation of chondrocytes to slow down the process of KOA cartilage degeneration [62, 63]. Meanwhile, the effective active components in XHLP have the characteristics of inhibiting chondrocyte apoptosis, reducing cartilage matrix catabolism, and eliminating the production of pro-inflammatory cytokines. Therefore, XHLP can play a role in inhibiting chondrocyte apoptosis, anti-inflammatory, and other mechanisms to treat or delay the degeneration of KOA cartilage.

Finally, the binding ability of the active components of XHLP with molecular target domains was verified in this study through molecular docking technology. The results demonstrated that quercetin can form a stable docking model with CCND1, JUN, MMP9, and MMP2, respectively; beta-sitosterol can dock with JUN and MMP2 protein ligands; Izoteolin and ellagic acid can form stable docking with JUN, MMP9, and MMP2 protein ligands. According to the calculations with Vina, LeDock, and DS software, quercetin can form the most stable molecular docking with the 6ESM ligand of MMP9 protein; MMP9 is a pro-inflammatory protein, and can exert a direct anti-inflammatory effect to treat KOA cartilage degeneration when combined with quercetin, which is the effective active component in XHLP. However, this study still has some limitations. First, the results of this study lack cell and animal experimental verification. Another limitation is the false negative due to the targets rooted from a different database, which may have biased impact because of different experimental conditions. We are also aware of the public databases that have limited information. Moreover, Earth Dragon failed to retrieve the active components in the TCMSP database and Batman-TCM database.

Conclusion

In this study, the potentially effective active components, molecular targets, and action pathways of XHLP in the treatment of KOA cartilage degeneration were analyzed using bioinformatics. Besides, 71 effective active components were discovered for potential treatment of KOA; 1672 differentially expressed genes in cartilage degeneration with KOA were discovered through bio-information analysis. Mutual mapping was used to obtain 33 potential targets for the XHLP treatment of KOA. This study discovered 33 molecular targets such as VEGFA, CCND1, MYC, JUN, MMP9, and MMP2. These molecular targets are associated with the functions of “inflammatory response,” “mechanical stimulus,” “TNF signaling pathway,” and “PI3K-Akt signaling pathway,” involved in signaling pathways related to inflammation and chondrocyte apoptosis. Moreover, molecular docking verification was conducted to verify that the active components of XHLP can form a stable docking model with potential molecular targets. To sum up, the XHLP treatment of KOA plays a synergistic role in anti-inflammatory and inhibiting chondrocyte apoptosis through the characteristics of “multi-component, multi-target, multi-pathway.”

Availability of data and materials

All the data will be available upon motivated request to the corresponding author of the present paper.

Abbreviations

OA:

Osteoarthritis

KOA:

Knee osteoarthritis

XHLP:

Xiao huoluo pills

TCMSP:

Traditional Chinese medicine system pharmacology

OB:

Oral bioavailability

DL:

Drug-likeness

GEO:

Gene expression omnibus

DEGs:

Differentially expressed genes

PPI:

Protein-protein interaction

BC:

Betweenness centrality

NC:

Neighborhood connectivity

GO:

Gene ontology

KEGG:

Kyoto encyclopedia of genes and genomes

CC:

Cellular component

MF:

Molecular function

BP:

Biological process

References

  1. Cross M, Smith E, Hoy D, Nolte S, Ackerman I, Fransen M, et al. The global burden of hip and knee osteoarthritis: estimates from the global burden of disease 2010 study. Ann Rheum Dis. 2014;73(7):1323–30. https://doi.org/10.1136/annrheumdis-2013-204763.

    Article  PubMed  Google Scholar 

  2. Kim JR, Yoo JJ, Kim HA. Therapeutics in Osteoarthritis Based on an Understanding of Its Molecular  Pathogenesis. Int J Mol Sci. 2018;19(3).

  3. Tang X, Wang S, Zhan S, Niu J, Tao K, Zhang Y, et al. The prevalence of symptomatic knee osteoarthritis in China: results from the China Health and Retirement Longitudinal Study. Arthritis Rheumatol. 2016;68(3):648–53. https://doi.org/10.1002/art.39465.

    Article  PubMed  Google Scholar 

  4. Sayre EC, Esdaile JM, Kopec JA, Singer J, Wong H, Thorne A, et al. Specific manifestations of knee osteoarthritis predict depression and anxiety years in the future: Vancouver Longitudinal Study of Early Knee Osteoarthritis. BMC Musculoskelet Disord. 2020;21(1):467. https://doi.org/10.1186/s12891-020-03496-8.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Latourte A, Kloppenburg M, Richette P. Emerging pharmaceutical therapies for osteoarthritis. Nat Rev Rheumatol. 2020;16(12):673–88. https://doi.org/10.1038/s41584-020-00518-6.

    Article  PubMed  Google Scholar 

  6. UK NCGC. Osteoarthritis: Care and Management in Adults. London: National Institute for Health and Care Excellence (UK); 2014.

  7. Glyn-Jones S, Palmer AJ, Agricola R, Price AJ, Vincent TL, Weinans H, et al. Osteoarthritis. Lancet. 2015;386(9991):376–87. https://doi.org/10.1016/S0140-6736(14)60802-3.

    Article  CAS  PubMed  Google Scholar 

  8. Hochberg MC, Altman RD, April KT, Benkhalti M, Guyatt G, McGowan J, et al. American College of Rheumatology 2012 recommendations for the use of nonpharmacologic and pharmacologic therapies in osteoarthritis of the hand, hip, and knee. Arthritis Care Res (Hoboken). 2012;64(4):465–74. https://doi.org/10.1002/acr.21596.

    Article  CAS  Google Scholar 

  9. Yang M, Jiang L, Wang Q, Chen H, Xu G. Traditional Chinese medicine for knee osteoarthritis: an overview of systematic review. Plos One. 2017;12(12):e189884.

    Google Scholar 

  10. Wang L, Zhang XF, Zhang X, Guo DY, Duan YW, Wang ZC, et al. Evaluation of the therapeutic effect of traditional Chinese medicine on osteoarthritis: a systematic review and meta-analysis. Pain Res Manag. 2020;2020:5712187.

    PubMed  PubMed Central  Google Scholar 

  11. The State Pharmacopoeia Committee of People’s Republic of China. Pharmacopoeia of People’s Republic of China, vol. 1. Beijing, China: Chemical Industry Press; 2020. p. 40–1, 58,127,193,233,247–248.

    Google Scholar 

  12. Zhou Y, Cui Y, Liu Y. Shen G: [Bioavailability study on xiaohuoluo pills]. Zhongguo Zhong Yao Za Zhi. 1995;20(3):159–61 191.

    PubMed  Google Scholar 

  13. Pan J. XLZD: Immunosuppresive, antioxidative, anti-inflammatory and analgesic effect of Xiahuoluo pills. Chin J Clin Rehabil. 2006;2006(47):183–8.

    Google Scholar 

  14. Cui P, Han H, Wang R, Yang L. Identification and determination of Aconitum alkaloids in Aconitum herbs and Xiaohuoluo pill using UPLC-ESI-MS. Molecules. 2012;17(9):10242–57. https://doi.org/10.3390/molecules170910242.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kiss T, Orvos P, Bánsághi S, Forgo P, Jedlinszki N, Tálosi L, et al. Identification of diterpene alkaloids from Aconitum napellus subsp. firmum and GIRK channel activities of some Aconitum alkaloids. Fitoterapia. 2013;90:85–93. https://doi.org/10.1016/j.fitote.2013.07.010.

    Article  CAS  PubMed  Google Scholar 

  16. Seo M, Lee JH, Baek M, Kim MA, Ahn MY, Kim SH, et al. A novel role for earthworm peptide Lumbricusin as a regulator of neuroinflammation. Biochem Biophys Res Commun. 2017;490(3):1004–10. https://doi.org/10.1016/j.bbrc.2017.06.154.

    Article  CAS  PubMed  Google Scholar 

  17. Li C, Chen M, Li X, Yang M, Wang Y, Yang X. Purification and function of two analgesic and anti-inflammatory peptides from coelomic fluid of the earthworm, Eisenia foetida. Peptides. 2017;89:71–81. https://doi.org/10.1016/j.peptides.2017.01.016.

    Article  CAS  PubMed  Google Scholar 

  18. Ge CY, Zhang JL. Bioactive sesquiterpenoids and steroids from the resinous exudates of Commiphora myrrha. Nat Prod Res. 2019;33(3):309–15. https://doi.org/10.1080/14786419.2018.1448811.

    Article  CAS  PubMed  Google Scholar 

  19. Liang Y. LPHQ: Acetyl-11-keto-beta-boswellic acid regulates the activities of matrix metalloproteinases-1,-2,-9. Chin J Pathophysiol. 2009;25(10):2004–11.

    CAS  Google Scholar 

  20. Ammon HP. Modulation of the immune system by Boswellia serrata extracts and boswellic acids. Phytomedicine. 2010;17(11):862–7. https://doi.org/10.1016/j.phymed.2010.03.003.

    Article  PubMed  Google Scholar 

  21. Li S, Zhang B. Traditional Chinese medicine network pharmacology: theory, methodology and application. Chin J Nat Med. 2013;11(2):110–20. https://doi.org/10.1016/S1875-5364(13)60037-0.

    Article  PubMed  Google Scholar 

  22. Lee WY, Lee CY, Kim YS, Kim CE. The methodological trends of traditional herbal medicine employing network Pharmacology. Biomolecules. 2019;9(8).

  23. Santolini M, Barabási AL. Predicting perturbation patterns from the topology of biological networks. Proc Natl Acad Sci U S A. 2018;115(27):E6375–83. https://doi.org/10.1073/pnas.1720589115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Chen S, Kang J, Xing Y, Zhao Y, Milton DK. Estimating large covariance matrix with network topology for high-dimensional biomedical data. Comput Stat Data Anal. 2018;127:82–95. https://doi.org/10.1016/j.csda.2018.05.008.

    Article  Google Scholar 

  25. Ru J, Li P, Wang J, Zhou W, Li B, Huang C, et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. 2014;6(1):13. https://doi.org/10.1186/1758-2946-6-13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Wan Y, Xu L, Liu Z, Yang M, Jiang X, Zhang Q, et al. Utilising network pharmacology to explore the underlying mechanism of Wumei Pill in treating pancreatic neoplasms. BMC Complement Altern Med. 2019;19(1):158. https://doi.org/10.1186/s12906-019-2580-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Xu X, Zhang W, Huang C, Li Y, Yu H, Wang Y, et al. A novel chemometric method for the prediction of human oral bioavailability. Int J Mol Sci. 2012;13(6):6964–82. https://doi.org/10.3390/ijms13066964.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Assenov Y, Ramírez F, Schelhorn SE, Lengauer T, Albrecht M. Computing topological parameters of biological networks. Bioinformatics. 2008;24(2):282–4. https://doi.org/10.1093/bioinformatics/btm554.

    Article  CAS  PubMed  Google Scholar 

  30. Metsalu T, Vilo J. ClustVis: a web tool for visualizing clustering of multivariate data using principal component analysis and heatmap. Nucleic Acids Res. 2015;43(W1):W566–70. https://doi.org/10.1093/nar/gkv468.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Heberle H, Meirelles GV, Da SF, Telles GP, Minghim R. InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams. BMC Bioinformatics. 2015;16(1):169. https://doi.org/10.1186/s12859-015-0611-3.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13. https://doi.org/10.1093/nar/gky1131.

    Article  CAS  PubMed  Google Scholar 

  33. Martin A, Ochagavia ME, Rabasa LC, Miranda J, Fernandez-de-Cossio J, Bringas R. BisoGenet: a new tool for gene network building, visualization and analysis. BMC Bioinformatics. 2010;11(1):91. https://doi.org/10.1186/1471-2105-11-91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Yu G, Wang LG, Han Y. He QY: clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010;31(2):455–61. https://doi.org/10.1002/jcc.21334.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Nguyen NT, Nguyen TH, Pham T, Huy NT, Bay MV, Pham MQ, et al. Autodock Vina adopts more accurate binding poses but Autodock4 forms better binding affinity. J Chem Inf Model. 2020;60(1):204–11. https://doi.org/10.1021/acs.jcim.9b00778.

    Article  CAS  PubMed  Google Scholar 

  37. Cao H, Sun Y, Wang L, Zhao C, Fu J, Zhang A. Understanding the microscopic binding mechanism of hydroxylated and sulfated polybrominated diphenyl ethers with transthyretin by molecular docking, molecular dynamics simulations and binding free energy calculations. Mol Biosyst. 2017;13(4):736–49. https://doi.org/10.1039/C6MB00638H.

    Article  CAS  PubMed  Google Scholar 

  38. Zhang L, Shi X, Huang Z, Mao J, Mei W, Ding L, et al. Network pharmacology approach to uncover the mechanism governing the effect of Radix achyranthis Bidentatae on osteoarthritis. BMC Complement Med Ther. 2020;20(1):121. https://doi.org/10.1186/s12906-020-02909-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Zhang J, Zhang Q, Chen X, Liu Y, Xue J, Dahan A, et al. Revealing synergistic mechanism of multiple components in Gandi capsule for diabetic nephropathy therapeutics by network pharmacology. Evid Based Complement Alternat Med. 2018;2018:6503126.

    PubMed  PubMed Central  Google Scholar 

  40. Veber DF, Johnson SR, Cheng HY, Smith BR, Ward KW, Kopple KD. Molecular properties that influence the oral bioavailability of drug candidates. J MED CHEM. 2002;45(12):2615–23. https://doi.org/10.1021/jm020017n.

    Article  CAS  PubMed  Google Scholar 

  41. Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):S11.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Samara B, Deliorman M, Sukumar P, Qasaimeh MA. Cryopreservable arrays of paper-based 3D tumor models for high throughput drug screening. Lab Chip. 2021;21(5):844–54. https://doi.org/10.1039/D0LC01300E.

    Article  CAS  PubMed  Google Scholar 

  43. Huang C, Zheng C, Li Y, Wang Y, Lu A, Yang L. Systems pharmacology in drug discovery and therapeutic insight for herbal medicines. Brief Bioinform. 2014;15(5):710–33. https://doi.org/10.1093/bib/bbt035.

    Article  PubMed  Google Scholar 

  44. Hu Y, Gui Z, Zhou Y, Xia L, Lin K, Xu Y. Quercetin alleviates rat osteoarthritis by inhibiting inflammation and apoptosis of chondrocytes, modulating synovial macrophages polarization to M2 macrophages. Free Radic Biol Med. 2019;145:146–60. https://doi.org/10.1016/j.freeradbiomed.2019.09.024.

    Article  CAS  PubMed  Google Scholar 

  45. Gabay O, Sanchez C, Salvat C, Chevy F, Breton M, Nourissat G, et al. Stigmasterol: a phytosterol with potential anti-osteoarthritic properties. Osteoarthritis Cartilage. 2010;18(1):106–16. https://doi.org/10.1016/j.joca.2009.08.019.

    Article  CAS  PubMed  Google Scholar 

  46. Kim KA, Lee IA, Gu W, Hyam SR, Kim DH. β-Sitosterol attenuates high-fat diet-induced intestinal inflammation in mice by inhibiting the binding of lipopolysaccharide to toll-like receptor 4 in the NF-κB pathway. Mol Nutr Food Res. 2014;58(5):963–72. https://doi.org/10.1002/mnfr.201300433.

    Article  CAS  PubMed  Google Scholar 

  47. Lin Z, Lin C, Fu C, Lu H, Jin H, Chen Q, et al. The protective effect of Ellagic acid (EA) in osteoarthritis: an in vitro and in vivo study. Biomed Pharmacother. 2020;125:109845. https://doi.org/10.1016/j.biopha.2020.109845.

    Article  CAS  PubMed  Google Scholar 

  48. Zupan J, Vrtačnik P, Cör A, Haring G, Weryha G, Visvikis-Siest S, et al. VEGF-A is associated with early degenerative changes in cartilage and subchondral bone. Growth Factors. 2018;36(5-6):263–73. https://doi.org/10.1080/08977194.2019.1570926.

    Article  CAS  PubMed  Google Scholar 

  49. Hamilton JL, Nagao M, Levine BR, Chen D, Olsen BR, Im HJ. Targeting VEGF and its receptors for the treatment of osteoarthritis and associated Pain. J Bone Miner Res. 2016;31(5):911–24. https://doi.org/10.1002/jbmr.2828.

    Article  PubMed  Google Scholar 

  50. Vadalà G, Russo F, Musumeci M, Giacalone A, Papalia R, Denaro V. Targeting VEGF-A in cartilage repair and regeneration: state of the art and perspectives. J Biol Regul Homeost Agents. 2018;32(6 Suppl. 1):217–24.

    PubMed  Google Scholar 

  51. Zan PF, Yao J, Wu Z, Yang Y, Hu S, Li GD. Cyclin D1 gene silencing promotes IL-1β-induced apoptosis in rat chondrocytes. J Cell Biochem. 2018;119(1):290–9. https://doi.org/10.1002/jcb.26172.

    Article  CAS  PubMed  Google Scholar 

  52. Zou J, Li XL, Shi ZM, Xue JF. Effects of C-myc gene silencing on interleukin-1β-induced rat chondrocyte cell proliferation, apoptosis and cytokine expression. J Bone Miner Metab. 2018;36(3):286–96. https://doi.org/10.1007/s00774-017-0845-4.

    Article  CAS  PubMed  Google Scholar 

  53. Rhee J, Park SH, Kim SK, Kim JH, Ha CW, Chun CH, et al. Inhibition of BATF/JUN transcriptional activity protects against osteoarthritic cartilage destruction. Ann Rheum Dis. 2017;76(2):427–34. https://doi.org/10.1136/annrheumdis-2015-208953.

    Article  CAS  PubMed  Google Scholar 

  54. Teng P, Liu Y, Dai Y, Zhang H, Liu WT, Hu J. Nicotine attenuates osteoarthritis pain and matrix metalloproteinase-9 expression via the α7 Nicotinic Acetylcholine Receptor. J Immunol. 2019;203(2):485–92. https://doi.org/10.4049/jimmunol.1801513.

    Article  CAS  PubMed  Google Scholar 

  55. Jackson MT, Moradi B, Smith MM, Jackson CJ, Little CB. Activation of matrix metalloproteinases 2, 9, and 13 by activated protein C in human osteoarthritic cartilage chondrocytes. Arthritis Rheumatol. 2014;66(6):1525–36. https://doi.org/10.1002/art.38401.

    Article  CAS  PubMed  Google Scholar 

  56. Xue M, McKelvey K, Shen K, Minhas N, March L, Park SY, et al. Endogenous MMP-9 and not MMP-2 promotes rheumatoid synovial fibroblast survival, inflammation and cartilage degradation. Rheumatology (Oxford). 2014;53(12):2270–9. https://doi.org/10.1093/rheumatology/keu254.

    Article  CAS  Google Scholar 

  57. Eskelinen A, Tanska P, Florea C, Orozco GA, Julkunen P, Grodzinsky AJ, et al. Mechanobiological model for simulation of injured cartilage degradation via pro-inflammatory cytokines and mechanical stimulus. Plos Comput Biol. 2020;16(6):e1007998. https://doi.org/10.1371/journal.pcbi.1007998.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Zevenbergen L, Gsell W, Chan DD, Vander SJ, Himmelreich U, Neu CP, et al. Functional assessment of strains around a full-thickness and critical sized articular cartilage defect under compressive loading using MRI. Osteoarthritis Cartilage. 2018;26(12):1710–21. https://doi.org/10.1016/j.joca.2018.08.013.

    Article  CAS  PubMed  Google Scholar 

  59. Myller K, Korhonen RK, Töyräs J, Salo J, Jurvelin JS, Venäläinen MS. Computational evaluation of altered biomechanics related to articular cartilage lesions observed in vivo. J Orthop Res. 2019;37(5):1042–51. https://doi.org/10.1002/jor.24273.

    Article  CAS  PubMed  Google Scholar 

  60. Wang T, He C. Pro-inflammatory cytokines: the link between obesity and osteoarthritis. Cytokine Growth Factor Rev. 2018;44:38–50. https://doi.org/10.1016/j.cytogfr.2018.10.002.

    Article  CAS  PubMed  Google Scholar 

  61. Liao CR, Wang SN, Zhu SY, Wang YQ, Li ZZ, Liu ZY, et al. Advanced oxidation protein products increase TNF-α and IL-1β expression in chondrocytes via NADPH oxidase 4 and accelerate cartilage degeneration in osteoarthritis progression. Redox Biol. 2020;28:101306. https://doi.org/10.1016/j.redox.2019.101306.

    Article  CAS  PubMed  Google Scholar 

  62. Sun K, Luo J, Guo J, Yao X, Jing X, Guo F. The PI3K/AKT/mTOR signaling pathway in osteoarthritis: a narrative review. Osteoarthritis Cartilage. 2020;28(4):400–9. https://doi.org/10.1016/j.joca.2020.02.027.

    Article  CAS  PubMed  Google Scholar 

  63. Ke H, Mou X, Xia Q. Remifentanil repairs cartilage damage and reduces the degradation of cartilage matrix in post-traumatic osteoarthritis, and inhibits IL-1β-induced apoptosis of articular chondrocytes via inhibition of PI3K/AKT/NF-κB phosphorylation. Ann Transl Med. 2020;8(22):1487. https://doi.org/10.21037/atm-20-6000.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Grateful acknowledgement is made to my supervisor Professor Liang Du and Wang haibin who gave me considerable help by means of suggestion, comments, and criticism. In addition, I deeply appreciate the contribution to this thesis made in various ways by my friends and classmates.

Funding

This research was financially supported by the National Natural Science Foundation of China (No.81774339;No.82074462) to Wang Haibin.

Author information

Authors and Affiliations

Authors

Contributions

C.W.J and L.T.Y were responsible for the design and implementation of the study presented. C.W.J and H.Q were responsible for the inclusions and literature review. Y.P, Z.G.Y, H.F.Y,W.Z.H,P.H, and L.B.L were responsible for the acquisition of the data. C.W.J and L.T.Y prepared the initial draft of the manuscript. L.D and W.H.B gave critical feedback during the study and critically revised the submitted manuscript for important intellectual content. All authors have read and approved the final manuscript to be submitted.

Corresponding authors

Correspondence to Du Liang or Haibin Wang.

Ethics declarations

Ethical approval and consent to participate

This study was conducted in agreement with the Declaration of Helsinki and its later amendments or comparable ethical standards and had been approved by the ethics board of The First Affiliated Hospital of Guangzhou University of Chinese Medicine (No: Y《2020》119).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, W., Lin, T., He, Q. et al. Study on the potential active components and molecular mechanism of Xiao Huoluo Pills in the treatment of cartilage degeneration of knee osteoarthritis based on bioinformatics analysis and molecular docking technology. J Orthop Surg Res 16, 460 (2021). https://doi.org/10.1186/s13018-021-02552-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13018-021-02552-w

Keywords