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.

Computational analysis of sodium-dependent phosphate transporter SLC20A1/PiT1 gene identifies missense variations C573F, and T58A as high-risk deleterious SNPs.

SLC20A1/PiT1 is a sodium-dependent inorganic phosphate transporter, initially recognized as the retroviral receptor for Gibbon Ape Leukemia Virus in humans. SNPs in SLC20A1 is associated with Combined Pituitary Hormone Deficiency and Sodium Lithium Counter transport. Using in silico techniques, we have screened the nsSNPs for their deleterious effect on the structure and function of SLC20A1. Screening with sequence and structure-based tools on 430 nsSNPs, filtered 17 nsSNPs which are deleterious. To evaluate the role of these SNPs, protein modeling and MD simulations were performed. A comparative analysis of model generated with SWISS-MODEL and AlphaFold shows that many residues are in the disallowed region of Ramachandran plot. Since SWISS-MODEL structure has a 25-residue deletion, the AlphaFold structure was used to perform MD simulation for equilibration and structure refinement. Further, to understand perturbation of energetics, we performed in silico mutagenesis and ΔΔG calculation using FoldX on MD refined structures, which yielded SNPs that are neutral (3), destabilizing (12) and stabilizing (2) on protein structure. Furthermore, to elucidate the impact of SNPs on structure, we performed MD simulations to discern the changes in RMSD, Rg, RMSF and LigPlot of interacting residues. RMSF profiles of representative SNPs revealed that A114V (neutral) and T58A (positive) were more flexible & C573F (negative) was more rigid compared to wild type, which is also reflected in the changes in number of local interacting residues in LigPlot and ΔΔG. Taken together, our results show that SNPs can lead to structural perturbations and impact the function of SLC20A1 with potential implications for disease.Communicated by Ramaswamy H. Sarma.