Abstract
To investigate the potential of ginsenosides in treating osteoporosis, ginsenoside compound K (GCK) was selected to explore the potential targets and mechanism based on network pharmacology (NP). Based on text mining from public databases, 206 and 6590 targets were obtained for GCK and osteoporosis, respectively, in which 138 targets were identified as co-targets of GCK and osteoporosis using intersection analysis. Five central gene clusters and key genes (STAT3, PIK3R1, VEGFA, JAK2 and MAP2K1) were identified based on Molecular Complex Detection (MCODE) analysis through constructing a protein-protein interaction network using the STRING database. Gene Ontology (GO) analysis implied that phosphatidylinositol-related biological process, molecular modification and function may play an important role for GCK in the treatment of osteoporosis. Function and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis suggested that the c-Fms-mediated osteoclast differentiation pathway was one of the most important mechanisms for GCK in treating osteoporosis. Meanwhile, except for being identified as key targets based on cytoHubba analysis using Cytoscape software, MAPK and PI3K-related proteins were enriched in the downstream of the c-Fms-mediated osteoclast differentiation pathway. Molecular docking further confirmed that GCK could interact with the cavity on the surface of a c-Fms protein with the lowest binding energy (-8.27 Kcal/moL), and their complex was stabilized by hydrogen bonds (Thr578 (1.97 Å), Leu588 (2.02 Å, 2.18 Å), Ala590 (2.16 Å, 2.84 Å) and Cys 666 (1.93 Å)), van der Waals and alkyl hydrophobic interactions. Summarily, GCK could interfere with the occurrence and progress of osteoporosis through the c-Fms-mediated MAPK and PI3K signaling axis regulating osteoclast differentiation.