Icariin: A Potential Molecule for Treatment of Knee Osteoarthritis

Abstract

Background: Knee osteoarthritis (KOA) is a degenerative disease that develops over time. Icariin (ICA) has a positive effect on KOA, although the mechanism is unknown. To investigate drug-disease connections and processes, network pharmacology is commonly used. The molecular mechanisms of ICA for the treatment of KOA were investigated using network pharmacology, molecular docking and literature research approaches in this study.

Methods: We gathered KOA-related genes using the DisGeNET database, the OMIM database, and GEO microarray data. TCMSP database, Pubchem database, TTD database, SwissTargetPrediction database, and Pharmmapper database were used to gather ICA-related data. Following that, a protein-protein interaction (PPI) network was created. Using the Metascape database, we performed GO and KEGG enrichment analyses. After that, we built a targets-pathways network. Furthermore, molecular docking confirms the prediction. Finally, we looked back over the last 5 years of literature on icariin for knee osteoarthritis to see if the findings of this study were accurate.

Results: core targets relevant to KOA treatment include TNF, IGF1, MMP9, PTGS2, ESR1, MMP2 and so on. The main biological process involved regulation of inflammatory response, collagen catabolic process, extracellular matrix disassembly and so on. The most likely pathways involved were the IL-17 signaling pathway, TNF signaling pathway, Estrogen signaling pathway.

Conclusion: ICA may alleviate KOA by inhibiting inflammation, cartilage breakdown and extracellular matrix degradation. Our study reveals the molecular mechanism of ICA for the treatment of KOA, demonstrating its potential value for further research and as a new drug.

Introduction

Knee osteoarthritis (KOA) is a chronic progressive disease with pain, swelling, and even deformity of the knee joint as the main clinical manifestation (Martel-Pelletier et al., 2016; Primorac et al., 2020). OA is the 11th leading cause of disability worldwide, which not only affects patients physically and mentally, but also represents a serious social burden (Palazzo et al., 2016). An epidemiological survey showed that the number of KOA patients in the United States doubled from the 1900’s to the 2000’s (Wallace et al., 2017). According to research in 2010 (Cross et al., 2014), the global standard age prevalence of KOA was 3.8%, and the prevalence increased sharply with age. KOA pathology mainly involves articular cartilage and synovium (Michael et al., 2010). In the physiological state, the cartilage matrix is in a dynamic balance of cytokine-controlled production and breakdown. When mechanically altered or stimulated by inflammation, cartilage degrades and progresses to KOA if the compensatory capacity of the joint is exceeded (Michael et al., 2010). The pathogenesis of KOA, however, is not yet clear. To better guide the clinical treatment of KOA, its pathogenesis needs to be studied urgently.

Icariin (ICA) is one of the main components of the Chinese medicine Epimedium. Epimedium was first recorded in the Yellow Emperor’s Inner Canon in 100 B.C. Traditional Chinese medicine believed that Epimedium could strengthen muscles and bones, and it was widely used in clinical practice. Due to the development of modern medicine, its active ingredients have been further studied. Modern research has discovered that ICA, one of its active ingredients, has strong biological activity on various systems, such as the nervous system (Jin et al., 2019), reproductive system (Zhang et al., 2021), skeletal system (Pan et al., 2017; Wang et al., 2018) and so on. ICA was found to improve bone loss caused by estrogen deficiency through the IGF-1 pathway, the ERK pathway and the JNK pathway (Song et al., 2013; Zhou et al., 2021). It also has the effect of preventing cartilage degradation and promoting cartilage repair, as proved by some in vivo and in vitro experiments (Wei et al., 2016; Zhang et al., 2019). The most widely used drugs in non-surgical treatment are NSAIDs such as cyclooxygenase-2 (COX-2) inhibitors (Bannuru et al., 2019). However, the long-term use of NSAIDs leads to gastrointestinal complications such as indigestion, gastritis and heartburn (Curtis et al., 2019). Since there is no specific treatment for KOA, there is an urgent need to find a safe and reliable treatment method.

Network pharmacology analysis is one of the most commonly used methods to study the pharmacology of traditional Chinese medicine (TCM). By searching the databases, constructing the networks, enrichment analysis and other steps, we can get the components, potential targets, correlation with diseases and other information of TCM (Yuan et al., 2017). This leads to a more in-depth study of the possible mechanisms of TCM in the treatment of diseases. Molecular docking technology can verify the binding between TCM molecules and predicted targets to validate the accuracy of these targets (Pinzi and Rastelli 2019). Therefore, this study explores the molecular mechanism of ICA for KOA based on network pharmacology, molecular docking, and literature research and the flow chart is shown in Figure 1.

FIGURE 1

Materials and Methods

Collection of Osteoarthritis Related Genes

Firstly, microarray information on differential DNA expression between normal subjects and KOA patients was collected from the GEO database (https://www.ncbi.nlm.nih.gov/geo/) (Edgar et al., 2002), Series GSE169077. The microarray data was normalized with log₂ (fold change) > 1 or log₂ (fold change) < -1 as the standard, which was considered to have significant differences in the expression. The KOA related genes were collected by integrating databases through searches with the key words “knee osteoarthritis” and “Homo sapiens,” followed by taking the intersection. Two databases include: DisGeNET database (https://www.disgenet.org/web/DisGeNET/menu/home) (Piñero et al., 2017) and OMIM database (https://omim.org/) (Amberger et al., 2015). Next, a visual volcano map and a gene expression heat map were created using Graphpad Prism 8.0.2 (GraphPad Software, San Diego, California United States, www.graphpad.com) and TBtools (https://github.com/CJ-Chen/TBtools) (Chen et al., 2020).

Potential Target Genes of Icariin in the Treatment of Knee Osteoarthritis

To evaluate the biological information of ICA, it is chemical name, molecular weight, oral bioavailability (OB) and drug likeness (DL) were obtained from the TCMSP database (https://tcmspw.com/tcmsp.php) (Ru et al., 2014). Studies have shown that OB≥40% is considered to have good bioavailability; DL ≥0.18 is considered to be a drug-like compound (Liu et al., 2013). Then the 2D and 3D molecular structures of ICA were obtained through the Pubchem database (https://pubchem.ncbi.nlm.nih.gov/) (Kim et al., 2021).

Next, the potential targets of ICA were predicted using the SwissTargetPrediction database (http://swisstargetprediction.ch/) (Gfeller et al., 2013), Therapeutic Target Database (http://db.idrblab.net/ttd/) (Wang Y. et al., 2020) and Pharmmapper database (http://www.lilab-ecust.cn/pharmmapper/) (Wang et al., 2017), followed by taking intersection set. After that, all targets were converted into gene symbols standardized through the Uniprot database (https://www.uniprot.org/) (Morgat et al., 2020).

Construction of Protein Interaction Network

The potential targets of ICA were matched with KOA-related genes to obtain the core targets. The venn diagram was drawn through venny 2.1.0 (https://bioinfogp.cnb.csic.es/tools/venny/) (Oliveros 2007-2015). Next, the targets were analyzed by the String database (https://string-db.org/) (Szklarczyk et al., 2019). After that, the PPI network was edited with Cytoscape 3.7.2 software (http://www.cytoscape.org) and the degree and edge betweenness of the nodes were calculated by the NetworkAnalyzer tool.

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Enrichment Analysis

Based on the above data, GO and KEGG enrichment analysis was performed in the Metascape database (http://metascape.org/gp/index.html) (Zhou et al., 2019). The Gene Ontology (GO) enrichment analysis includes cellular composition, biological process, and molecular function, whereas the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis can profile metabolic pathways of genes in the cell as well as their systemic function. Unlike the last update of the David database in 2016, the Metascape database was updated on 1 February 2021. After this, a visual bubble diagram was drawn using the Bioinformatic online tool (http://bioinformatics.com.cn/).

Construction of Targets-Pathways Network

To investigate the interactions between key genes and pathways of ICA and KOA, the ICA-KOA-targets-pathways network was constructed using the results of metascape database analysis. The degree and edge betweenness were analyzed by topological analysis tools in the cytoscape 3.7.2 software.

Molecular Docking

Proteins corresponding to the top ten key genes were selected to dock with the ICA molecule to verify their affinity. Crystal structures of proteins were downloaded from the PDB database (https://www.rcsb.org/) (Burley et al., 2021) with the search criteria: Homo sapiens, refine resolution<3.0 Å and release date 2015–2021. then we imported the protein structures into PyMOL 2.2.0 software (https://pymol.org) (Schrodinger 2015) for modification including removal of water molecules, separation of ligands and addition of hydrogen. Using AutoDockTools 1.5.6 software (http://autodock.scripps.edu/) (Goodsell and Olson 1990) to add charge to protein molecule and set up a docking grid box at the center of the molecule (Seeliger and de Groot 2010). Molecular docking was performed using AutoDock Vina 1.1.2 software (http://vina.scripps.edu/) (Trott and Olson 2010). We chose it as the software for molecular docking because it is efficient, accurate, and has a new way of evaluation. Finally, after analyzing the binding energy of the molecule, choosing the conformation with the lowest binding energy, and observing the formation of hydrogen bonds, we drew a binding diagram with PyMOL 2.2.0.

Search and Analysis of Literature

To validate the predictions, we also searched for relevant studies over the last 5 years, including in vivo studies, in vitro studies and reviews. The search terms we used were “icariin”, “osteoarthritis”, “knee osteoarthritis”, “cartilage”, “chondrocytes”, and “extracellular matrix”. To ensure the relevance of the literature to this study, the results were collated by two independent researchers.

Result

Acquisition of Knee Osteoarthritis Related Genes

A total of 1105 differentially expressed genes were identified by analysis of gene microarray Series GSE169077, which was used to draw a volcano may and heat map (Figure 2). The results showed a large difference in gene expression between OA patients and normal subjects. Then 720 out of 9329 genes from the Genecards database (relevance score >4.26), 368 genes from the Disgenet database, 56 genes from the OMIM database, and finally 398 genes were obtained.

FIGURE 2

Analysis of Icariin and Collection of its Potential Targets

Molecular formula, molecular weight, 2D structure, 3D structures, and ADME information of ICA were collected from the TCMSP database and the Pubchem database (Table 1). The OB (41.58%) and DL (0.61) of ICA proved its high bioavailability and drug-like property, which can be used as a drug molecule for the next step of research. A total of 334 potential targets were obtained by combining and removing duplicate ones from the SwissTargetPrediction database, the Pharmmapper database, and the Therapeutic Target Database. Next, the ICA targets were intersected with the KOA genes to obtain the key genes and plotted a Venn diagram (Figure 3).

TABLE 1

Molecular Name	Molecular Formula	Molecular Weight	2D Structure	3D Structure
Icariin	C33H40O15	676.66 g/mol

Basic information of ICA.

FIGURE 3

Protein-Protein Interaction Network Analysis

The obtained 21 common genes were uploaded to the String database to produce the PPI network. After that, the results were imported into cytoscape 3.7.2 software and the topological parameters of the network was analyzed. The result showed that the network had a total of 18 nodes and 78 edges (Figure 4A). According to the analysis results, TNF was the most important target in the network (Degree 17, BC 0.24, CC 0.95). To find the key genes of ICA in treatment of KOA, a core network with 5 nodes and 10 edges was filtered by Degree>11 (Figure 4B). The other 4 genes were MMP9, IGF1, PTGS2 and ESR1. Finally, the 5 core genes screened by the NetworkAnalyzer tool and their Degree values and Betweenness Centrality values were shown in Table 2.

FIGURE 4

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Enrichment Analysis

351 GO items were collected through the analysis of the Metascape database (p 0.01, count 3), while the KEGG entries were 112 (p 0.01, count 3). The top 15 most significant entries were filtered in ascending p-value order and displayed in a visual bubble diagram, and all results were shown if less than 15 (Figure 5). The results show that the most significant molecular function included serine-type endopeptidase activity, serine-type peptidase activity, metallopeptidase activity, nuclear receptor activity and so on. While the biological process was primarily involved in processes such as collagen catabolism, extracellular matrix disassembly, collagen metabolic process, and other extracellular matrix-related processes. Another aspect is the processes associated with inflammation, such as regulation of inflammatory responses, regulation of neuroinflammatory responses and so on. KEGG results showed that genes shared by ICA and KOA were mainly enriched in the IL-17 signaling pathway, TNF signaling pathway and Estrogen signaling pathway, excluding other disease pathways (Figure 5D; Table 3). Then a targets-pathways network was constructed based on the results (Figure 6).

FIGURE 5

FIGURE 6

Molecular Docking Analysis

Based on the data above, the 21 core target genes were individually imported into the Autodock Vina 1.1.2 software to dock with the ICA. Affinity was the score used by the software to calculate the binding ability. It is generally believed that affinity < −7 kcal/mol indicated a stronger binding activity (Ferreira et al., 2015; Gromiha et al., 2017; Cui et al., 2020). After calculation, it was found that 17 out of 21 genes met the criteria of affinity < −7 kcal/mol (Table 4). The docking diagrams of ICA and genes were drawn using Pymol software (Figure 7). In addition, we found the hydrogen bonds formed between molecules and represented by yellow dotted lines in the diagrams.

FIGURE 7

Results of the Literature Collection

Through our study of the literature, we screened twenty-eight articles from the last 5 years relevant to this study. Of these, 8 papers mentioned that icariin inhibited inflammation, 6 papers indicated that icariin reduced chondrocyte apoptosis and 4 papers indicated that icariin inhibited extracellular matrix degradation. In addition, 9 papers mentioned that icariin promotes chondrocyte proliferation, differentiation, cellular autophagy or cartilage repair (Figure 8A). In the literature on the underlying mechanisms, 3 papers suggested that icariin can treat OA via the NF-κB signaling pathway, 3 papers mentioned that icariin inhibited the expression of TNF-α, 3 papers mentioned that icariin downregulated the expression of MMPs, while 2 papers mentioned that icariin reduced the inflammatory response caused by IL-β (Figure 8B).

FIGURE 8

Discussion

Despite the significant developments in the diagnosis and treatment of KOA, single drug development for a single target is still slow (Wang et al., 2015). This is because KOA is a complex pathological process involving multiple substances and regulated by multiple pathways (Abramoff and Caldera 2020). Coincidentally, traditional Chinese medicine has a multi-target and multi-pathway biological effect. As a traditional Chinese medicine, epimedium is widely used in clinical practice. It has also been widely used in the treatment of knee osteoarthritis from ancient times to the present day with good healing properties. As research progressed, ICA, one of the main components of Epimedium, was found to have strong biological activity. Due to its low price and easy availability, it has a high potential value in the treatment of KOA. In the present study, based on network pharmacology, we found that icariin has good bioavailability. Meanwhile, we assembled a number of key targets and pathways, as well as biological functions, and reviewed the literature over the last 5 years to validate them.

Icariin May be Useful for the Treatment of Knee Osteoarthritis by Inhibiting Inflammation

Inflammation is a major factor in the development and progression of KOA. On the one hand, inflammatory mediators promote cartilage apoptosis and extracellular matrix degradation. On the other hand, the products of chondrocyte apoptosis and matrix breakdown induce the production of inflammatory substances that further aggravate the KOA condition (Martel-Pelletier et al., 2016). These inflammatory mediators include tumour necrosis factor (TNF), interleukin (IL), nitric oxide (NO), prostaglandin E2 (PGE₂) and so on (Fei et al., 2019; Li et al., 2020).

Our study found that the key targets of icariin for the treatment of KOA include TNF, PTGS2 as well as RELA (NF-κB). TNF is a powerful pro-inflammatory factor that binds to two types of receptors, including TNFR1 and TNFR2, which are expressed on the surface of both chondrocytes and synovial cells (Ansari et al., 2020) (Figure 9). TNF binding to the receptor activates the downstream NF-κB signaling pathway, thereby up-regulating the synthesis of inflammatory mediators (PTGS2, NO, PGE₂, IL-6, etc.) and matrix degrading enzymes (MMPs, ADAMTS). In addition, it can inhibit the synthesis of type II collagen (Wojdasiewicz et al., 2014). A study reported that icariin can inhibit NF-κB activation and nuclear translocation triggered by TNF-α, as well as reduce chondrocyte apoptosis and inflammation (Mi et al., 2018). In addition, according to the study of Pengzhen Wang et al. (2020), icariin inhibited the activation of NF-κB and HIF-2α in the chondrocytes of OA mice and attenuated the inflammatory response driven by TNF-α. In addition, the expression of MMPs was inhibited. These findings confirm that icariin reduces KOA inflammation by down-regulating the TNF signalling pathway through inhibition of TNF-α, NF-κB. In addition to the above proven targets, we also identified PTGS2 as a potential target for icariin. This could be an area for further research.

FIGURE 9

Icariin may contribute to the treatment of KOA by promoting cartilage repair and attenuating cartilage extracellular matrix degradation.

In KOA, chondrocyte apoptosis and extracellular matrix degradation are severe pathological changes (Choi et al., 2019). The extracellular matrix is mostly made up of type 2 collagenase and aggrecan, which can be degraded by matrix-degrading enzymes like MMPs and ADAMTS (Guilak et al., 2018). And this is associated with the accumulation of inflammatory substances in KOA joints including TNF-α, IL-6, IL-1β (Krishnan and Grodzinsky 2018). Degradation of the extracellular matrix promotes chondrocyte apoptosis at the same time (Komori 2016).

Our research has found substantial evidence that icariin can promote cartilage repair in KOA over the last 5 years (Wang et al., 2016; Liu et al., 2019; Zhang et al., 2019). Other studies have shown that icariin regulates the TDP-43 signaling pathway and NF-κB signaling pathway to reduce chondrocyte apoptosis (Mi et al., 2018; Huang et al., 2019). In addition, Tang Y, et al. showed that icariin mediated the PI3K/AKT/mTOR pathway to activate chondrocyte autophagy (Tang et al., 2021).

In addition to the above findings, we found that insulin-like growth factor 1 (IGF-1) is also one of the key targets of icariin. It is closely related to the proliferation, differentiation and maintenance of the phenotype of cartilage (Oh and Chun 2003). One study found that icariin enhances IGF-1 signalling, promotes osteogenic differentiation and has estrogen-like effects (Zhou et al., 2021). This suggests that icariin promotes chondrocyte proliferation and differentiation via IGF-1, thereby repairing KOA-damaged cartilage.

We also found that icariin also binds to the estrogen receptors and activates estrogen signaling pathway (Figure 9). The estrogen signaling pathway also plays a significant role in the skeletal system. In particular, post-menopausal women have a high incidence of arthritis and cardiovascular disease. Estrogen was found to have good anti-inflammatory effects and reduce cartilage damage (Martín-Millán and Castañeda 2013). This may be one of the reasons why there are more KOA patients among females than males over the age of 50.

Network pharmacology as a new technology has a wide range of applications in exploring the molecular mechanisms of Chinese medicine. In this study, we predicted the molecular mechanism of ICA for the treatment of KOA and obtained 21 relevant targets. Further analysis obtained the core targets as TNF, MMP9, IGF1, PTGS2, ESR1 and so on. After GO and KEGG analysis, we found that the key pathway for ICA treatment of KOA is the IL-17 signaling pathway, TNF signaling pathway and Estrogen signaling pathway with a regulatory effect on the extracellular matrix degradation and inflammation. Molecular docking verified the affinity of key targets. Finally, we reviewed the literature related to icariin for treatment of KOA in the last 5 years, and the results are generally consistent with the prediction.

Conclusion

Taken together, ICA is able to inhibit KOA inflammation, chondrocyte apoptosis, matrix degradation, and promote cartilage repair through multiple pathways at multiple targets. Through network pharmacology and molecular docking technology, we obtained the interaction relationship between genes, protein molecules and revealed the molecular mechanism of ICA for the treatment of KOA, making it possible as a new drug.

References

• 1

  AbramoffB.CalderaF. E. (2020). Osteoarthritis: Pathology, Diagnosis, and Treatment Options. Med. Clin. North. Am.104 (2), 293–311. 10.1016/j.mcna.2019.10.007

  • CrossRef
  • Google Scholar

• 2

  AmbergerJ. S.BocchiniC. A.SchiettecatteF.ScottA. F.HamoshA. (2015). OMIM.org: Online Mendelian Inheritance in Man (OMIM®), an Online Catalog of Human Genes and Genetic Disorders. Nucleic Acids Res.43 (Database issue), D789–D798. 10.1093/nar/gku1205

  • CrossRef
  • Google Scholar

• 3

  AnsariM. Y.AhmadN.HaqqiT. M. (2020). Oxidative Stress and Inflammation in Osteoarthritis Pathogenesis: Role of Polyphenols. Biomed. Pharmacother.129, 110452. 10.1016/j.biopha.2020.110452

  • CrossRef
  • Google Scholar

• 4

  BannuruR. R.OsaniM. C.VaysbrotE. E.ArdenN. K.BennellK.Bierma-ZeinstraS. M. A.et al (2019). OARSI Guidelines for the Non-surgical Management of Knee, Hip, and Polyarticular Osteoarthritis. Osteoarthritis Cartilage27 (11), 1578–1589. 10.1016/j.joca.2019.06.011

  • CrossRef
  • Google Scholar

• 5

  BurleyS. K.BhikadiyaC.BiC.BittrichS.ChenL.CrichlowG. V.et al (2021). RCSB Protein Data Bank: Powerful New Tools for Exploring 3D Structures of Biological Macromolecules for Basic and Applied Research and Education in Fundamental Biology, Biomedicine, Biotechnology, Bioengineering and Energy Sciences. Nucleic Acids Res.49 (D1), D437–D51. 10.1093/nar/gkaa1038

  • CrossRef
  • Google Scholar

• 6

  ChenC.ChenH.ZhangY.ThomasH. R.FrankM. H.HeY.et al (2020). TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant13 (8), 1194–1202. 10.1016/j.molp.2020.06.009

  • CrossRef
  • Google Scholar

• 7

  ChoiM. C.JoJ.ParkJ.KangH. K.ParkY. (2019). NF-κB Signaling Pathways in Osteoarthritic Cartilage Destruction. Cells8 (7). 10.3390/cells8070734

  • CrossRef
  • Google Scholar

• 8

  CrossM.SmithE.HoyD.NolteS.AckermanI.FransenM.et al (2014). The Global burden of Hip and Knee Osteoarthritis: Estimates from the Global burden of Disease 2010 Study. Ann. Rheum. Dis.73 (7), 1323–1330. 10.1136/annrheumdis-2013-204763

  • CrossRef
  • Google Scholar

• 9

  CuiQ.ZhangY. L.MaY. H.YuH. Y.ZhaoX. Z.ZhangL. H.et al (2020). A Network Pharmacology Approach to Investigate the Mechanism of Shuxuening Injection in the Treatment of Ischemic Stroke. J. Ethnopharmacol257, 112891. 10.1016/j.jep.2020.112891

  • CrossRef
  • Google Scholar

• 10

  CurtisE.FuggleN.ShawS.SpoonerL.NtaniG.ParsonsC.et al (2019). Safety of Cyclooxygenase-2 Inhibitors in Osteoarthritis: Outcomes of a Systematic Review and Meta-Analysis. Drugs Aging36 (Suppl. 1), 25–44. 10.1007/s40266-019-00664-x

  • CrossRef
  • Google Scholar

• 11

  EdgarR.DomrachevM.LashA. E. (2002). Gene Expression Omnibus: NCBI Gene Expression and Hybridization Array Data Repository. Nucleic Acids Res.30 (1), 207–210. 10.1093/nar/30.1.207

  • CrossRef
  • Google Scholar

• 12

  FeiJ.LiangB.JiangC.NiH.WangL. (2019). Luteolin Inhibits IL-1β-induced Inflammation in Rat Chondrocytes and Attenuates Osteoarthritis Progression in a Rat Model. Biomed. Pharmacother.109, 1586–1592. 10.1016/j.biopha.2018.09.161

  • CrossRef
  • Google Scholar

• 13

  FerreiraL. G.Dos SantosR. N.OlivaG.AndricopuloA. D. (2015). Molecular Docking and Structure-Based Drug Design Strategies. Molecules20 (7), 13384–13421. 10.3390/molecules200713384

  • CrossRef
  • Google Scholar

• 14

  GfellerD.MichielinO.ZoeteV. (2013). Shaping the Interaction Landscape of Bioactive Molecules. Bioinformatics29 (23), 3073–3079. 10.1093/bioinformatics/btt540

  • CrossRef
  • Google Scholar

• 15

  GoodsellD. S.OlsonA. J. (1990). Automated Docking of Substrates to Proteins by Simulated Annealing. Proteins8 (3), 195–202. 10.1002/prot.340080302

  • CrossRef
  • Google Scholar

• 16

  GromihaM. M.YugandharK.JemimahS. (2017). Protein-protein Interactions: Scoring Schemes and Binding Affinity. Curr. Opin. Struct. Biol.44, 31–38. 10.1016/j.sbi.2016.10.016

  • CrossRef
  • Google Scholar

• 17

  GuilakF.NimsR. J.DicksA.WuC. L.MeulenbeltI. (2018). Osteoarthritis as a Disease of the Cartilage Pericellular Matrix. Matrix Biol.71-72, 40–50. 10.1016/j.matbio.2018.05.008

  • CrossRef
  • Google Scholar

• 18

  HuangH.ZhangZ. F.QinF. W.TangW.LiuD. H.WuP. Y.et al (2019). Icariin Inhibits Chondrocyte Apoptosis and Angiogenesis by Regulating the TDP-43 Signaling Pathway. Mol. Genet. Genomic Med.7 (4), e00586. 10.1002/mgg3.586

  • CrossRef
  • Google Scholar

• 19

  JinJ.WangH.HuaX.ChenD.HuangC.ChenZ. (2019). An Outline for the Pharmacological Effect of Icariin in the Nervous System. Eur. J. Pharmacol.842, 20–32. 10.1016/j.ejphar.2018.10.006

  • CrossRef
  • Google Scholar

• 20

  KimS.ChenJ.ChengT.GindulyteA.HeJ.HeS.et al (2021). PubChem in 2021: New Data Content and Improved Web Interfaces. Nucleic Acids Res.49 (D1), D1388–D95. 10.1093/nar/gkaa971

  • CrossRef
  • Google Scholar

• 21

  KomoriT. (2016). Cell Death in Chondrocytes, Osteoblasts, and Osteocytes. Int. J. Mol. Sci.17 (12), 2045. 10.3390/ijms17122045

  • CrossRef
  • Google Scholar

• 22

  KrishnanY.GrodzinskyA. J. (2018). Cartilage Diseases. Matrix Biol.71-72, 51–69. 10.1016/j.matbio.2018.05.005

  • CrossRef
  • Google Scholar

• 23

  LiL.LiZ.LiY.HuX.ZhangY.FanP. (2020). Profiling of Inflammatory Mediators in the Synovial Fluid Related to Pain in Knee Osteoarthritis. BMC Musculoskelet. Disord.21 (1), 99. 10.1186/s12891-020-3120-0

  • CrossRef
  • Google Scholar

• 24

  LiuH.WangJ.ZhouW.WangY.YangL. (2013). Systems Approaches and Polypharmacology for Drug Discovery from Herbal Medicines: an Example Using Licorice. J. Ethnopharmacol146 (3), 773–793. 10.1016/j.jep.2013.02.004

  • CrossRef
  • Google Scholar

• 25

  LiuY.YangJ.LuoZ.LiD.LuJ.WangQ.et al (2019). Development of an Injectable Thiolated Icariin Functionalized Collagen/hyaluronic Hydrogel to Promote Cartilage Formation In Vitro and In Vivo. J. Mater. Chem. B7 (17), 2845–2854. 10.1039/c9tb00211a

  • CrossRef
  • Google Scholar

• 26

  Martel-PelletierJ.BarrA. J.CicuttiniF. M.ConaghanP. G.CooperC.GoldringM. B.et al (2016). Osteoarthritis. Nat. Rev. Dis. Primers2, 16072. 10.1038/nrdp.2016.72

  • CrossRef
  • Google Scholar

• 27

  Martín-MillánM.CastañedaS. (2013). Estrogens, Osteoarthritis and Inflammation. Jt. bone Spine80 (4), 368–373. 10.1016/j.jbspin.2012.11.008

  • CrossRef
  • Google Scholar

• 28

  MiB.WangJ.LiuY.LiuJ.HuL.PanayiA. C.et al (2018). Icariin Activates Autophagy via Down-Regulation of the NF-Κb Signaling-Mediated Apoptosis in Chondrocytes. Front. Pharmacol.9, 605. 10.3389/fphar.2018.00605

  • CrossRef
  • Google Scholar

• 29

  MichaelJ. W.Schlüter-BrustK. U.EyselP. (2010). The Epidemiology, Etiology, Diagnosis, and Treatment of Osteoarthritis of the Knee. Dtsch Arztebl Int.107 (9), 152–162. 10.3238/arztebl.2010.0152

  • CrossRef
  • Google Scholar

• 30

  MorgatA.LombardotT.CoudertE.AxelsenK.NetoT. B.GehantS.et al (2020). Enzyme Annotation in UniProtKB Using Rhea. Bioinformatics36 (6), 1896–1901. 10.1093/bioinformatics/btz817

  • CrossRef
  • Google Scholar

• 31

  OhC. D.ChunJ. S. (2003). Signaling Mechanisms Leading to the Regulation of Differentiation and Apoptosis of Articular Chondrocytes by Insulin-like Growth Factor-1. J. Biol. Chem.278 (38), 36563–36571. 10.1074/jbc.M304857200

  • CrossRef
  • Google Scholar

• 32

  OliverosJ. C. (2007-2015). Venny. An Interactive Tool for Comparing Lists with Venn's Diagrams. Available at: https://bioinfogp.cnb.csic.es/tools/venny/index.html.

  • Google Scholar

• 33

  PalazzoC.NguyenC.Lefevre-ColauM. M.RannouF.PoiraudeauS. (2016). Risk Factors and burden of Osteoarthritis. Ann. Phys. Rehabil. Med.59 (3), 134–138. 10.1016/j.rehab.2016.01.006

  • CrossRef
  • Google Scholar

• 34

  PanL.ZhangY.ChenN.YangL. (2017). Icariin Regulates Cellular Functions and Gene Expression of Osteoarthritis Patient-Derived Human Fibroblast-like Synoviocytes. Int. J. Mol. Sci.18 (12). 10.3390/ijms18122656

  • CrossRef
  • Google Scholar

• 35

  PiñeroJ.BravoÀ.Queralt-RosinachN.Gutiérrez-SacristánA.Deu-PonsJ.CentenoE.et al (2017). DisGeNET: a Comprehensive Platform Integrating Information on Human Disease-Associated Genes and Variants. Nucleic Acids Res.45 (D1), D833–D9. 10.1093/nar/gkw943

  • CrossRef
  • Google Scholar

• 36

  PinziL.RastelliG. (2019). Molecular Docking: Shifting Paradigms in Drug Discovery. Int. J. Mol. Sci.20 (18). 10.3390/ijms20184331

  • CrossRef
  • Google Scholar

• 37

  PrimoracD.MolnarV.RodE.JelečŽ.ČukeljF.MatišićV.et al (2020). Knee Osteoarthritis: A Review of Pathogenesis and State-Of-The-Art Non-operative Therapeutic Considerations. Genes (Basel)11 (8), 854. 10.3390/genes11080854

  • CrossRef
  • Google Scholar

• 38

  RuJ.LiP.WangJ.ZhouW.LiB.HuangC.et al (2014). TCMSP: a Database of Systems Pharmacology for Drug Discovery from Herbal Medicines. J. Cheminform6, 13. 10.1186/1758-2946-6-13

  • CrossRef
  • Google Scholar

• 39

  SchrodingerL. L. C. (2015). The PyMOL Molecular Graphics System. Version 1.8.

  • Google Scholar

• 40

  SeeligerD.de GrootB. L. (2010). Ligand Docking and Binding Site Analysis with PyMOL and Autodock/Vina. J. Comput. Aided Mol. Des.24 (5), 417–422. 10.1007/s10822-010-9352-6

  • CrossRef
  • Google Scholar

• 41

  SongL.ZhaoJ.ZhangX.LiH.ZhouY. (2013). Icariin Induces Osteoblast Proliferation, Differentiation and Mineralization through Estrogen Receptor-Mediated ERK and JNK Signal Activation. Eur. J. Pharmacol.714 (1-3), 15–22. 10.1016/j.ejphar.2013.05.039

  • CrossRef
  • Google Scholar

• 42

  SzklarczykD.GableA. L.LyonD.JungeA.WyderS.Huerta-CepasJ.et al (2019). STRING V11: Protein-Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-wide Experimental Datasets. Nucleic Acids Res.47 (D1), D607–D13. 10.1093/nar/gky1131

  • CrossRef
  • Google Scholar

• 43

  TangY.LiY.XinD.ChenL.XiongZ.YuX. (2021). Icariin Alleviates Osteoarthritis by Regulating Autophagy of Chondrocytes by Mediating PI3K/AKT/mTOR Signaling. Bioengineered12 (1), 2984–2999. 10.1080/21655979.2021.1943602

  • CrossRef
  • Google Scholar

• 44

  TrottO.OlsonA. J. (2010). AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading. J. Comput. Chem.31 (2), 455–461. 10.1002/jcc.21334

  • CrossRef
  • Google Scholar

• 45

  WallaceI. J.WorthingtonS.FelsonD. T.JurmainR. D.WrenK. T.MaijanenH.et al (2017). Knee Osteoarthritis Has Doubled in Prevalence since the Mid-20th century. Proc. Natl. Acad. Sci. U S A.114 (35), 9332–9336. 10.1073/pnas.1703856114

  • CrossRef
  • Google Scholar

• 46

  WangK.XuJ.HunterD. J.DingC. (2015). Investigational Drugs for the Treatment of Osteoarthritis. Expert Opin. Investig. Drugs24 (12), 1539–1556. 10.1517/13543784.2015.1091880

  • CrossRef
  • Google Scholar

• 47

  WangP.MengQ.WangW.ZhangS.XiongX.QinS.et al (2020). Icariin Inhibits the Inflammation through Down-Regulating NF-Κb/hif-2α Signal Pathways in Chondrocytes. Biosci. Rep.40 (11). 10.1042/BSR20203107

  • CrossRef
  • Google Scholar

• 48

  WangP.ZhangF.HeQ.WangJ.ShiuH. T.ShuY.et al (2016). Flavonoid Compound Icariin Activates Hypoxia Inducible Factor-1α in Chondrocytes and Promotes Articular Cartilage Repair. PloS one11 (2), e0148372. 10.1371/journal.pone.0148372

  • CrossRef
  • Google Scholar

• 49

  WangX.ShenY.WangS.LiS.ZhangW.LiuX.et al (2017). PharmMapper 2017 Update: a Web Server for Potential Drug Target Identification with a Comprehensive Target Pharmacophore Database. Nucleic Acids Res.45 (W1), W356–W60. 10.1093/nar/gkx374

  • CrossRef
  • Google Scholar

• 50

  WangY.ZhangS.LiF.ZhouY.ZhangY.WangZ.et al (2020). Therapeutic Target Database 2020: Enriched Resource for Facilitating Research and Early Development of Targeted Therapeutics. Nucleic Acids Res.48 (D1), D1031–D41. 10.1093/nar/gkz981

  • CrossRef
  • Google Scholar

• 51

  WangZ.WangD.YangD.ZhenW.ZhangJ.PengS. (2018). The Effect of Icariin on Bone Metabolism and its Potential Clinical Application. Osteoporos. Int.29 (3), 535–544. 10.1007/s00198-017-4255-1

  • CrossRef
  • Google Scholar

• 52

  WeiC. C.PingD. Q.YouF. T.QiangC. Y.TaoC. (2016). Icariin Prevents Cartilage and Bone Degradation in Experimental Models of Arthritis. Mediators Inflamm.2016, 9529630. 10.1155/2016/9529630

  • CrossRef
  • Google Scholar

• 53

  WojdasiewiczP.PoniatowskiŁ. A.SzukiewiczD. (2014). The Role of Inflammatory and Anti-inflammatory Cytokines in the Pathogenesis of Osteoarthritis. Mediators Inflamm.2014, 561459. 10.1155/2014/561459

  • CrossRef
  • Google Scholar

• 54

  YuanH.MaQ.CuiH.LiuG.ZhaoX.LiW.et al (2017). How Can Synergism of Traditional Medicines Benefit from Network Pharmacology?Molecules22 (7). 10.3390/molecules22071135

  • CrossRef
  • Google Scholar

• 55

  ZhangJ.ZhangC.LiuA.JiQ.RenL.MaC.et al (2021). Synthesis of Icariin-Zinc and its Protective Effect on Exercise Fatigue and Reproductive System Related Glands in Male Rats. Front. Pharmacol.12, 611722. 10.3389/fphar.2021.611722

  • CrossRef
  • Google Scholar

• 56

  ZhangJ.ZhangD.WuC.LiuA.ZhangC.JiaoJ.et al (2019). Icariin-conditioned Serum Engineered with Hyaluronic Acid Promote Repair of Articular Cartilage Defects in Rabbit Knees. BMC Complement. Altern. Med.19 (1), 155. 10.1186/s12906-019-2570-0

  • CrossRef
  • Google Scholar

• 57

  ZhouL.PoonC. C.-W.WongK.-Y.CaoS.DongX.ZhangY.et al (2021). Icariin Ameliorates Estrogen-Deficiency Induced Bone Loss by Enhancing IGF-I Signaling via its Crosstalk with Non-genomic ERα Signaling. Phytomedicine82, 153413. 10.1016/j.phymed.2020.153413

  • CrossRef
  • Google Scholar

• 58

  ZhouY.ZhouB.PacheL.ChangM.KhodabakhshiA. H.TanaseichukO.et al (2019). Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat. Commun.10 (1), 1523. 10.1038/s41467-019-09234-6

  • CrossRef
  • Google Scholar