Multi-omics Data Integration Analysis Identified Therapeutic Targets and Potential Reuse Drugs for Osteoporosis

Aims: To facilitate drug discovery and development for the treatment of osteoporosis Background: With global aging, osteoporosis has become a common problem threatening the health of the elderly. It is of important clinical value to explore new targets for drug intervention and develop promising drugs for the treatment of osteoporosis. Objective: To understand the major molecules that mediate the communication between the cell populations of bone marrow-derived mesenchymal stem cells (BM-MSCs) in osteoporosis and osteoarthritis patients and identify potential reusable drugs for the treatment of osteoporosis. Methods: Single-cell RNA sequencing (scRNA-seq) data of BM-MSCs in GSE147287 dataset were classified using the Seurat package. CellChat was devoted to analyzing the ligand-receptor pairs (LR pairs) contributing to the communication between BM-MSCs subsets. The LR pairs that were differentially expressed between osteoporosis samples and control samples and significantly correlated with immune score were screened in the GSE35959 dataset, and the differentially expressed gene in both GSE35959 and GSE13850 data sets were identified as targets from a single ligand or receptor. The therapeutic drugs for osteoporosis were screened by network proximity method, and the top-ranked drugs were selected for molecular docking and molecular dynamics simulation with the target targets. Results: Twelve subsets of BM-MSCs were identified, of which CD45-BM-MSCS_4, CD45-BM- MSCS_5, and CD45+ BM-MSCs_5 subsets showed significantly different distributions between osteoporosis samples and osteoarthritis samples. Six LR pairs were identified in the bidirectional communication between these three BM-MSCs subsets and other BM-MSCs subsets. Among them, MIF-CD74 and ITGB2-ICAM2 were significantly correlated with the immune score. CD74 was identified as the target, and a total of 48 drugs targeting CD47 protein were identified. Among them, DB01940 had the lowest free energy binding score with CD74 protein and the binding state was very stable. Conclusion: This study provided a new network-based framework for drug reuse and identified initial insights into therapeutic agents targeting CD74 in osteoporosis, which may be meaningful for promoting the development of osteoporosis treatment.