In-silico analysis predicts disruption of normal angiogenesis as a causative factor in osteoporosis pathogenesis

Angiogenesis-osteogenesis coupling is critical for proper functioning and maintaining the health of bones. Any disruption in this coupling, associated with aging and disease, might lead to loss of bone mass. Osteoporosis (OP) is a debilitating bone metabolic disorder that affects the microarchitecture of bones, gradually leading to fracture. Computational analysis revealed that normal angiogenesis is disrupted during the progression of OP, especially postmenopausal osteoporosis (PMOP). The genes associated with OP and PMOP were retrieved from the DisGeNET database. Hub gene analysis and molecular pathway enrichment were performed via the Cytoscape plugins STRING, MCODE, CytoHubba, ClueGO and the web-based tool Enrichr. Twenty-eight (28) hub genes were identified, eight of which were transcription factors (HIF1A, JUN, TP53, ESR1, MYC, PPARG, RUNX2 and SOX9). Analysis of SNPs associated with hub genes via the gnomAD, I-Mutant2.0, MUpro, ConSurf and COACH servers revealed the substitution F201L in IL6 as the most deleterious. The IL6 protein was modeled in the SWISS-MODEL server and the substitution was analyzed via the YASARA FoldX plugin. A positive ΔΔG (1.936) of the F201L mutant indicates that the mutated structure is less stable than the wild-type structure is. Thirteen hub genes, including IL6 and the enriched molecular pathways were found to be profoundly involved in angiogenesis/endothelial function and immune signaling. Mechanical loading of bones through weight-bearing exercises can activate osteoblasts via mechanotransduction leading to increased bone formation. The present study suggests proper mechanical loading of bone as a preventive strategy for PMOP, by which angiogenesis and the immune status of the bone can be maintained. This in silico analysis could be used to understand the molecular etiology of OP and to develop novel therapeutic approaches.