Abstract
Osteoporosis (OP) is a degenerative skeletal disorder marked by a decrease in bone mineral density, which significantly raises the risk of fractures. The development of OP is closely associated with oxidative stress (OS), yet the specific molecular targets of OS in the context of OP, as well as the corresponding therapeutic strategies, are not well understood. This study aimed to systematically identify OS-related biomarkers for OP and explore potential therapeutic agents, providing new insights into the pathological mechanisms and precision therapy of OP. OP-related datasets (GSE56815 as the training set and GSE35958 as the validation set) were obtained from the gene expression omnibus database. Integrated with OS-related genes from GeneCards, 72 co-expressed differentially expressed genes associated with both OP and OS were identified. A protein-protein interaction network was constructed using search tool for the retrieval of interacting genes/proteins, and 9 hub genes were identified via Maximal Clique Centrality and MCODDE algorithms. Three core genes (SIRT2, neutrophil elastase [ELANE], and forkhead box O3 [FOXO3]) were further selected using random forest and support vector machine models. Diagnostic accuracy of key genes was evaluated by receiver operating characteristic curves. Immune infiltration analysis was performed with ImmuCellAI, while transcription factors were predicted through NetworkAnalyst to establish a miRNA-seed-transcription factor (TF) regulatory network. Potential drugs were screened from DGIdb, and molecular docking was conducted with CB-DOCK2 using ligand-protein structures from PubChem and UniProt. A FOXO3/sirtuin 1 (SIRT1) functional axis was subsequently constructed. The 72 differentially expressed genes showed enrichment in cell division, chromosome segregation, and cell cycle pathways. SIRT2, ELANE, and FOXO3 exhibited high diagnostic accuracy for OP (area under the curve > 0.7). Immune infiltration analysis revealed significant differences in key immune cells and subtypes between OP and controls. The miRNA-seed-TF network implicated miR-34a and miR-132. Molecular docking confirmed stable binding of 5 compounds to targets. Critically, the FOXO3-SIRT1 axis was identified as mediating oxidative stress in OP pathogenesis. FOXO3 and SIRT2 were upregulated; ELANE downregulated. This study identified SIRT2, ELANE, and FOXO3 as key OS-related biomarkers for OP and screened 5 potential therapeutic compounds. Integrating miRNA-TF network and drug sensitivity analysis, we propose the FOXO3/SIRT1 axis as central to OP oxidative stress, providing a foundation for mechanistic studies and targeted drug development.