Identification of osteoporosis ferroptosis-related markers and potential therapeutic compounds based on bioinformatics methods and molecular docking technology
Authors: Long Shi-Wei , Li Shi-Hong , Li Jian , He Yang , Tan Bo , Jing Hao-Han , Zheng Wei , Wu Juan 
Published: April 2024

Abstract
Research background and purpose
Osteoporosis (OP) is one of the most common bone diseases worldwide, characterized by low bone mineral density and susceptibility to pathological fractures, especially in postmenopausal women and elderly men. Ferroptosis is one of the newly discovered forms of cell death regulated by genes in recent years. Many studies have shown that ferroptosis is closely related to many diseases. However, there are few studies on ferroptosis in osteoporosis, and the mechanism of ferroptosis in osteoporosis is still unclear. This study aims to identify biomarkers related to osteoporosis ferroptosis from the GEO (Gene Expression Omnibus) database through bioinformatics technology, and to mine potential therapeutic small molecule compounds through molecular docking technology, trying to provide a basis for the diagnosis and treatment of osteoporosis in the future.

Materials and methods
We downloaded the ferroptosis-related gene set from the FerrDb database (http://www.zhounan.org/ferrdb/index.html), downloaded the data sets GSE56815 and GSE7429 from the GEO database, and used the R software “limma” package to screen differentially expressed genes (DEGs) from GSE56815, and intersected with the ferroptosis gene set to obtain ferroptosis-related DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed by the R software “clusterProfiler” package. The random forest model was further screened to obtain essential ferroptosis genes. R software “corrplot” package was used for correlation analysis of essential ferroptosis genes, and the Wilcox test was used for significance analysis. The lncRNA-miRNA-mRNA-TF regulatory network was constructed using Cytoscape software. The least absolute shrinkage and selection operator (LASSO) was used to construct a disease diagnosis model, and a Receiver operating characteristic (ROC) curve was drawn to evaluate the diagnostic performance, and then GSE7429 was used to verify the reliability of the diagnosis model. Molecular docking technology was used to screen potential small molecule compounds from the Drugbank database. Finally, a rat osteoporosis model was constructed, and peripheral blood mononuclear cells were extracted for qRT-PCR detection to verify the mRNA expression levels of crucial ferroptosis genes.

Result
Six DEGs related to ferroptosis were initially screened out. GO function and KEGG pathway enrichment analysis showed that ferroptosis-related DEGs were mainly enriched in signaling pathways such as maintenance of iron ion homeostasis, copper ion binding function, and ferroptosis. The random forest model identified five key ferroptosis genes, including CP, FLT3, HAMP, HMOX1, and SLC2A3. Gene correlation analysis found a relatively low correlation between these five key ferroptosis genes. The lncRNA-miRNA-mRNA-TF regulatory network shows that BAZ1B and STAT3 may also be potential molecules. The ROC curve of the disease diagnosis model shows that the model has a good diagnostic performance. Molecular docking technology screened out three small molecule compounds, including NADH, Midostaurin, and Nintedanib small molecule compounds. qRT-PCR detection confirmed the differential expression of CP, FLT3, HAMP, HMOX1 and SLC2A3 between OP and normal control group.

Conclusion
This study identified five key ferroptosis genes (CP, FLT3, HAMP, HMOX1, and SLC2A3), they were most likely related to OP ferroptosis. In addition, we found that the small molecule compounds of NADH, Midostaurin, and Nintedanib had good docking scores with these five key ferroptosis genes. These findings may provide new clues for the early diagnosis and treatment of osteoporosis in the future.