Bioinformatics analysis combined with molecular dynamics simulation validation to elucidate the potential molecular mechanisms of Jianshen Decoction for treatment of osteoporotic fracture

Abstract

Osteoporotic fracture (OPF) is a prevalent skeletal disease in the middle-aged and elderly. In clinical practice, Jianshen Decoction (JSD) has been used to treat OPFs. However, the specific effective components and mechanisms of JSD on OPF have not been explored. Therefore, this study used bioinformatics analysis combined with molecular dynamics simulation validation to explore the molecular mechanism of JSD treatment of OPF. Public databases (TCMSP, Batman TCM) were used to find the effective active components and corresponding target proteins of JSD (screening conditions: OB ≥ 30%, drug-likeness ≥ 0.18, half-life ≥ 4). Differentially expressed genes (DEGs) related to OPF lesions were obtained based on the gene expression omnibus database (screening conditions: adjustPvalue < .01, | log2FC | ≥ 1.0). The BisoGenet plug-in and the CytoNCA plug-in of Cytoscape were used to derive the potential core target proteins of JSD in the treatment of OPF. The JSD active ingredient target interaction network and the JSD-OPF target protein core network were constructed using the Cytoscape software. In addition, the R language Bioconductor package and clusterProfiler package were used to perform gene ontology (GO)/Kyoto Encylopedia Of Genes And Genome (KEGG) enrichment analysis on core genes to explain the biological functions and signal pathways of core proteins. Finally, molecular docking and molecular dynamics simulations were carried out through PyMOL, AutoDockTools 1.5.6, Vina, LeDock, Discovery Studio (DS) 2019, and other software to verify the binding ability of drug active ingredients and core target proteins. A total of 245 targets and 70 active components were identified. Through protein-protein interaction (PPI) network construction, 39 core targets were selected for further research. GO/KEGG enrichment analysis showed that the DNA-binding transcription factor binding, RNA polymerase II-specific DNA-binding transcription factor binding, MAPK signaling pathway, and ErbB signaling pathway were mainly involved. The results of molecular docking and molecular dynamics simulations supported the good interaction between MYC protein and Quercetin/Stigmasterol. In this study, bioinformatics, molecular docking, and molecular dynamics simulations were used for the first time to clarify the active components, molecular targets, and key biological pathways of JSD in the treatment of OPF, providing a theoretical basis for further research.