Association analysis of Vascular Endothelial Growth Factor-A (VEGF-A) polymorphism in rheumatoid arthritis using computational approaches.

Rheumatoid arthritis (RA), is marked by joint inflammation leading to pannus formation which results in cartilage destruction promoting bone erosion. The pathological hallmark of RA includes synovial hyperplasia and synovial angiogenesis. Active tissue neovascularization is observed in RA. Vascular endothelial Growth factor A (VEGFA), an endothelial cell-specific proangiogenic molecule is triggered by hypoxic cells and its levels are upregulated in RA. The aim of this study was to investigate functional and pathogenic VEGFA variants and to identify the impact of point mutation in VEGFA's interaction with VEGFR2 and how these polymorphisms affect the susceptibility and severity of RA. We investigated impact of these point mutations on the stability of VEGFA using various computational tools. These mutations were further identified by conservational profile as they are highly involved as structural and functional mutations. Furthermore, these selected variants were modelled and docked against targeted domain regions IGD2 and IGD3 of VEGFR2. Further molecular dynamic simulations were performed using Gromacs. Out of 168 nsSNPS, 19 were highlighted as highly pathogenic using insilico prediction tools. InterPro and ConSurf revealed domains and conserved variants respectively. After stability analysis, we concluded that almost all the mutations were responsible for decreasing the protein stability. HOPE predicted that all the selected damaging nsSNPs were present in the domain which is essential for the functioning of VEGFA protein. Constructed Ramachandran plot and ERRAT validated the quality of all the models. Based on the interactions predicted by STRING database, we performed Protein-Protein docking between VEGFA and VEGFR2. We found few conserved interactions and new polar contacts among wild-type and mutants with VEGFR2. From the simulations, we concluded that mutant R108Q was the most stabilizing mutant among all others whereas R82Q, C86Y, and R108W complexed with VEGFR2 were comparatively less stabilizing as compared to the wild type. This study provides insight into pathogenic nsSNPs that can affect VEGFA protein structure and function. These high-risk variants must be taken into consideration for genetic screening of patients suffering from RA.