LIANG Pengchen1,ZHOU Ziyan2,3, CHANG Qing2, YI Qingqing2, SUN Miaomiao2,3, TANG Yeling2,3, CAO Liou4,5, YANG Jie6
Journal of Jiangsu University(Medicine Edition).
2022, 32(1):
1-7,12.
Objective:The L1000 micro-interference data set was used to screen compounds for the treatment of osteoporosis, and the evaluation and cytological activity verification by multiple biological information methods were conducted, aiming to discover new osteoporosis treatment drugs.Methods:The data related to osteoporosis(OP)were retrieved in the GEO(Gene Expression Omnibus)database, and the difference in gene expression levels on the datasets was analyzed using limma package in R language. Pathway enrichment analysis of differential genes was performed using MetaScape data platform. The resulting differential OP gene expression profiles were analyzed to match potential OP therapeutic compounds in the L1000 dataset of the Connectivity Map database. The molecular docking study of the matching compounds and bone metabolism targets was conducted by Autodock_Vina software, and the ADMET(absorption, distribution, metabolism, excretion, and toxicity)characteristics of the matching compounds were calculated and analyzed by using pkCSM database. The activity of the candidate compound was studied by a cell proliferative activity experiment.Results:After analyzing the dataset of osteoporosis expression profile, 195 differential genes related to osteoporosis were obtained, including 127 genes with significantly up-regulated expression and 68 genes with significantly downregulated expression. There were 146 GO Biological Processes with differential gene enrichment;There were nine enriched KEGG signaling pathways and nine enriched Reactome Gene Sets. A total of 10 potential therapeutic compounds were matched in the L1000 dataset: diethylstilbestrol, PLX-4720, mezlocillin, KU-C103428N, rhapontin, maravir, quinine, MST-312, risperidone, and CO-102862. Through molecular docking studies, we found that the matching compounds have multi-target characteristics, and the main docking activities with the compounds are CAⅡ, PGR, ERβ, BMP2 and other OP targets. The ADMET properties of potential therapeutic compounds were calculated using the ADMET calculation module of the pkCSM database, which showed that mezlocillin had good ADME performance and low toxicity, with good comprehensive evaluation. Cell proliferation experiment showed that mezlocillin at 20 μg/mL had the most significant effect on the proliferation of osteoblasts.Conclusion:The candidate OP therapeutic compound, mezlocillin, was screened from the L1000 data set with good ADMET characteristics and bone cell proliferation promoting activity.