Abstract
Keywords
Introduction
The incidence of osteoarthritis increases with advancing age, 1 and it exerts an economic burden on both patients and society. 2 Osteoarthritis is a degenerative joint disease that involves cartilage, subchondral bone, synovium, and systemic inflammatory response. 3 Furthermore, research suggests that osteoarthritis is a systemic disease that involves alterations in lipid metabolism. 4 Currently, therapeutic drugs are unable to stop the progression of osteoarthritis, so better tolerated and safer osteoarthritis remedy medications are essential. 5 Icariin is the primary active ingredient of epimedium (the related information for icariin is shown in Table 1). Numerous studies have shown the positive effects of icariin on osteoarthritis. Additionally, in vitro and in vivo experiments, the biological function of icariin in regulating lipid metabolism has been confirmed. 6 Therefore, icariin may play a crucial role in lipid metabolism in osteoarthritis.
Related Information About Icariin.
Network pharmacology, including systems biology, network analysis, connectivity, redundancy, and pleiotropy, provides a platform for drug discovery. 7 The TCM network pharmacology has been established by applying the network pharmacology methods in traditional Chinese medicine. It enables us to understand the molecular mechanisms of Chinese medicine better. Furthermore, it helps Chinese medicine transform from empirical therapy to evidence-based medicine. 8
Our study used network pharmacology to elaborate on the potential mechanism of icariin in regulating lipid metabolism in osteoarthritis and intuitively verified its effectiveness through molecular docking. We constructed the PPI network and screened the core targets. Based on the relevant literature available, we aimed to elaborate on the possible mechanisms by which icariin regulates lipid metabolism in osteoarthritis. These early results laid the foundation for subsequent studies.
Materials and Methods
Database and Software-Related Materials
The TCMSP database provides drug targets of active compounds and related diseases to facilitate the analysis of the mechanism of drug action. 9
Swiss Target Prediction server provides biologically active small molecules with ligand-based target prediction. 10
Pharm Mapper server performs potential drug target prediction for any small molecule and identifies potential candidates from the database within a few hours of run time. 11
UniProt database summarizes the information of experimentally validated or computationally predicted protein function. 12
OMIM database stores genotypic and phenotypic information. 13
The DrugBank database provides detailed drug data and comprehensive drug target information. 14
GeneCards database provides complete data on human genes that have been annotated and predicted. 15
STRING database allows rapid examination of the potential evidence supporting specific protein-protein association. 16
Cytoscape software combines biomolecular interaction networks and provides a unified framework for representing and incorporating network models. 17
AutoDock is a computational docking and virtual screening software for interactions between small molecules and receptors. 18
AutoDock Vina is a computational program for molecular docking and virtual screening. 19
The PDB database is an open-source database providing valuable biological resources. It holds archives of biomolecular structures. 20
PubChem is also an open-source database that provides necessary chemical information on cheminformatics, chemical biology, medicinal chemistry, and several other fields. 21
PyMoL is a biomolecular visualization program for molecular modeling. 22
Screen the Targets of Icariin, Osteoarthritis, and Lipid Metabolism
The potential targets of icariin were obtained from the TCMSP database, Swiss Target Prediction server, and Pharm Mapper server. The targets of osteoarthritis and lipid metabolism were obtained in OMIM, Drugbank, and GeneCards databases. Human genes were identified from the target information with the UniProt database. Finally, those targets were combined and deduplicated.
Construct the PPI Network and Screen Core Targets
To gain the common targets, we intersected the targets. We imported the common targets into the STRING database (Table 2). By using the STRING database, we constructed the PPI network, in which the targets separated from the networks were deleted. Cytoscape 3.8.2 software was used to construct the network diagram. The CytoNCA plug-in was used to analyze the topological property of the network, acquiring the important nodes in the PPI network. Nodes in the PPI network represented common targets; the larger the node, the higher the degree value.
Information About Core Targets.
GO and KEGG Enrichment Analyses
The clusterProfiler R package was used to conduct GO and KEGG enrichment analyses. The results of the GO and KEGG enrichment analyses were chosen based on
Construct the Drug-Target-Signaling Pathway-Disease Network
To visualize the drug-target-signaling pathway-disease network, drugs, common targets, signaling pathways, and disease-related data were imported into cytoscape3.8.2 software. By constructing the “drug-target-signaling pathway-disease” network, we can better analyze how the medicine affects lipid metabolism in osteoarthritis.
Molecular Docking
The targets with degree value greater than the median were chosen as the main targets. First, the 3D structure of the core target protein was obtained from the RCSB PDB database (information about the core target proteins is shown in Table 3), and PubChem was used to retrieve the two-dimensional molecular structures of icariin. Then the structure of icariin was energy minimized in chem3D. AutoDock 4.2 software was used to dehydrate and hydrogenate the protein structure of the core target. In this process, the small molecule ligand icariin was hydrogenated, and the ligand flexible bond was configured to be rotatable. Molecular docking simulations and binding energy calculations were performed using the AutoDock vina program. The number of docking poses was set to 10, and the other AutoDock Vina parameters were set to default. Finally, PyMoL software was used to visualize the results.
Information About the Core Target Proteins.
Literature Validation
To confirm the accuracy of our findings, we searched PubMed for relevant studies published over the previous 5 years using the keywords “icariin,” “osteoarthritis,” “lipid metabolism,” and “cytokines.” This approach was carried out independently by 2 researchers.
Results
Targets of Icariin, Osteoarthritis, and Lipid Metabolism
From the TCMSP database, Swiss Target Prediction server, and Pharm Mapper server, a total of 300 targets of icariin were identified. From the OMIM, Drugbank, and GeneCards databases, 1110 osteoarthritis targets were found and 1351 relevant targets for lipid metabolism were obtained.
The PPI Network Diagram and Core Targets
There were 12 common targets among icariin, osteoarthritis, and lipid metabolism, as depicted in Figure 1A. A PPI network with 11 nodes and 31 edges was constructed using the STRING database and the Cytoscape software (Figure 1B). The degree value is reflected in the nodes‘ color and size. The larger the node and the redder the color, the higher the degree value, which shows that it may be the core target. Information on core targets can be seen in Figure 1C and Table 2.

(A) Common targets of icariin, osteoarthritis, and lipid metabolism, (B) PPI network diagram, and (C) the network map of core targets.
GO and KEGG Enrichment Analyses
In total, 810 biological processes, 17 cellular components, and 63 molecular functions were obtained from the GO enrichment analysis. The results showed that the common targets were primarily related to biological processes, including “tissue remodeling” and “steroid hormone response.” Figure 2 shows the top 10 biological processes, cellular components, and molecular functions. Through KEGG enrichment analysis, we discovered 62 signaling pathways. The top 20 KEGG pathways were intercepted from smallest to greatest based on

GO analyses: top 10 biological processes, cellular components, and molecular functions.

KEGG enrichment analyses: top 20 signaling pathways.
The Drug-Target-Signaling Pathway-Disease Network Diagram
The network diagram of “drug-target-signaling pathway-disease” (Figure 4) was constructed by Cytoscape 3.8.2 software, in which the nodes are represented by different colors and shapes; orange represents common targets, green represents the signaling pathways, red represents the drug, and the purple represents the disease.

Drug-target-signaling pathway-disease network diagram. The red node represents icariin, the yellow nodes represent gene targets, the green nodes represent signaling pathways, and the purple node represents disease.
Molecular Docking
Icariin was docked ten times with TNF, PTGS2, CCND1, MMP2, and ESR1. All docking binding free energies were smaller than or equal to −7 kcal/mol. It demonstrated that icariin could dock with TNF, PTGS2, CCND1, MMP2, and ESR1 without any external force and develop the optimal chimeric structure. Table 4 shows the minimal free energy of docking binding of icariin to TNF, PTGS2, CCND1, MMP2, and ESR1. Figure 5 shows the molecular docking results.

(A) Docking of ICA and TNF. (B) Docking of ICA and PTGS2. (C) Docking of ICA and CCND1. (D) Docking of ICA and MMP2. (E) Docking of ICA and ESR1.
The Minimum Free Energy of Docking Binding.
Results of Literature Validation
We searched the literature and found 38 papers that were relevant to our research. There were 24 studies that involved icariin in the treatment of osteoarthritis; 12 studies that involved the pathogenesis of osteoarthritis and lipid metabolism; and 2 studies that involved icariin in regulating lipid metabolism. TNFα and IL-1β were the most researched, receiving 6 and 7 mentions, respectively. Figure 6 depicts the findings of the literature validation. Among these targets, TNFα, MMP2, ESR1, and PTGS2 were consistent with our findings, implying that our findings were correct. Icariin may be beneficial in treating lipid metabolism in osteoarthritis by targeting TNFα, MMP2, ESR1, and PTGS2.

(A) Studies in the literature connected to icariin and osteoarthritis are represented in blue, research related to icariin and lipid metabolism is represented in orange, and studies related to osteoarthritis and lipid metabolism are represented in green. (B) The proportion of various targets specified in the relevant studies.
Discussion
Icariin is a Potential Modulator of Lipid Metabolism in Osteoarthritis
Icariin has anti-inflammatory properties. 23 Icariin was discovered to suppress chondrocyte autophagy and apoptosis. 24 Furthermore, a study on high cholesterol diet-induced atherosclerotic rats revealed that icariin could reduce TC, TG, and LDL-C while increasing HDL-C, significantly improving hyperlipidemia. 25 It is reported that icariin can lower blood TC and LDL-C levels among rabbits on a high-cholesterol diet. 26 One of the risk factors for the development of osteoarthritis is high serum cholesterol levels. 27 Statin lipid-lowering medicines can help slow down the course of osteoarthritis. 28 Based on these findings, it is not difficult to draw the following conclusions. Firstly, icariin has excellent anti-inflammatory and lipid-regulating abilities, which may positively affect osteoarthritis and lipid metabolism. Secondly, the development of osteoarthritis is related to lipid metabolism, and the regulation of lipid metabolism is beneficial to the treatment of osteoarthritis. Thirdly, icariin may be a potential modulator of lipid metabolism in osteoarthritis. After further validation of our study with the literature, we suggest that TNFα, PTGS2, MMP2, and ESR1 are the 4 possible drug targets of icariin in regulating lipid metabolism in osteoarthritis.
Icariin may Regulate Lipid Metabolism in Osteoarthritis by Inhibiting the Inflammatory Response
TNFα is a serum endotoxin-inducing factor, cachexia, and differentiation-inducing factor, mediating inflammatory responses and innate immunity. 29 The significance of TNFα in the pathophysiology of osteoarthritis has been demonstrated through its detection in the synovial fluid of osteoarthritis patients. 30 The level of TNFα also increases as the osteoarthritis score increases. 31 In addition, TNFα levels that are abnormally high in people with osteoarthritis are believed to be the leading cause of cartilage degradation. 32 COX-2 is another name for PTGS2 and is the crucial enzyme that converts arachidonic acid into prostaglandins. During inflammation, COX-2 levels increase, leading to increased prostaglandin production, and subsequent inflammation was significantly reduced with the administration of COX-2-specific inhibitors. 33 COX-2 is highly expressed in synovial vessels, lining cells, and fibroblast-like cells in synovial tissues of patients with various kinds of arthritis. 34 Inflammatory cytokines TNFα and COX-2 have been discovered to be inhibited by icariin. 35 These results suggested that icariin might exert an anti-inflammatory effect by inhibiting TNFα and COX-2.
Obesity is caused primarily by a disorder in lipid metabolism. Moreover, obesity is associated with an increased rate of osteoarthritis. This suggests that disorders of lipid metabolism further increase the risk of osteoarthritis. Higher levels of adipokines, particularly TNFα, will be found in Obesity-related osteoarthritis people. 36 TNFα also stimulates the synthesis of COX-2 and PGE-2, 37 aggravating the inflammatory response and speeding up the osteoarthritis process. Thus, TNF and COX-2 are important inflammatory factors in the lipid metabolism disorder of osteoarthritis. Their high expression in lipid metabolism disorders further aggravates the inflammatory response and accelerates the progression of osteoarthritis. Icariin can regulate osteoarthritis lipid metabolism by inhibiting the inflammatory response induced by TNFα and COX-2. It's interesting to note that TNFα is a key player in this chain of events, as it activates COX-2 and other enzymes to speed up the inflammatory response. TNFα is probably the primary pharmacological target.
Icariin may Regulate Lipid Metabolism in Osteoarthritis by Inhibiting Extracellular Matrix Breakdown
Adipose tissue secretes leptin which is an adipokine. Leptin is involved in the pathology of osteoarthritis, resulting in chondrocyte enlargement, cartilage degradation, and ossificariintion. 36 Leptin levels are greater than usual in people with osteoarthritis of the knee. 38 Leptin increases MMP-2, MMP-9, and collagen II levels by promoting matrix metalloproteinase expression. 39 Furthermore, higher mechanical stresses produce increased loading on cartilage and subchondral bone in patients with obesity-related osteoarthritis, 40 and the increased joint loading limits cartilage matrix synthesis. 41 MMP-2, a matrix metalloproteinase, is involved in extracellular matrix remodeling and protein breakdown. 42 Thus, MMP-2 has an irreplaceable role in lipid metabolism in osteoarthritis, as it aggravates cartilage wear by promoting extracellular matrix’s degradation, while MMP-2 levels are further increased in lipid metabolism disorders, which can further accelerate the osteoarthritis’ process. It is reported that icariin can suppress the expression of MMP-2. 43 As a result, icariin may modulate lipid metabolism in osteoarthritis by inhibiting extracellular matrix disintegration.
Estrogen Receptors: New Research Directions in Future
Estrogen receptors were discovered as one of the potential therapeutic targets. Estrogen works by binding to estrogen receptors in the body. Estrogen is one of the first hormones discovered and has an essential role in female reproduction. In addition, estrogen is involved in male reproduction and other systems, including neuroendocrine, vascular, skeletal, and immune systems in both males and females. 44 The incidence of osteoarthritis increases with age and is higher in women than in men, particularly after 50 years, indicating that estrogen plays an important role in the development of osteoarthritis. 45 Estrogen can also regulate lipogenesis and lipolysis in adipose tissue. Furthermore, it inhibits lipid synthesis and promotes lipolysis. 46 The estrogen receptor may be a new direction of research in the regulation of lipid metabolism in osteoarthritis by icariin due to the role of estrogen in osteoarthritis and lipid metabolism.
Conclusion
We investigated the potential molecular mechanism of icariin while regulating lipid metabolism in osteoarthritis. Five core targets, including TNF, PTGS2, CCND1, MMP2, and ESR1 were identified. Further literature research suggests that TNFα, MMP2, PTGS2, and ESR1 may be significant. We suggest that icariin may be a potential therapeutic agent for regulating lipid metabolism in osteoarthritis. Icariin could regulate lipid metabolism in osteoarthritis by inhibiting the inflammatory response. And icariin may regulate lipid metabolism in osteoarthritis via inhibiting extracellular matrix breakdown. Meanwhile, we suggested that ESR1 may be a significantly pharmacological target, and our findings serve as a reference for understanding the therapeutic mechanism.
