Molecular Dynamic Simulations Unravel the Underlying Impact of Missense Mutation in Autoimmunity Gene PTPN22 on Predisposition to Rheumatoid Arthritis

Genetic mutations in various proteins have been implicated with increased risk or severity of rheumatoid arthritis (RA) in different population groups. In the present case-control study, we have investigated the risk association of single nucleotide mutations present in some of the highly reported anti-inflammatory proteins and/or cytokines, with RA susceptibility in the Pakistani subjects. The study involves 310 ethnically and demographically similar participants from whom blood samples were taken and processed for DNA extraction. Through extensive data mining, 5 hotspot mutations reported in 4 genes, that is, interleukin (IL)-4 (-590; rs2243250), IL-10 (-592; rs1800872), IL-10 (-1082; rs1800896), PTPN22 (C1858T; rs2476601), and TNFAIP3 (T380G; rs2230926), were selected for RA susceptibility analyses using genotyping assays. The results demonstrated the association of only 2 DNA variants [rs2243250 (odds ratio, OR = 2.025, 95% confidence interval, CI = 1.357-3.002, P = 0.0005 Allelic) and rs2476601 (OR = 4.25, 95% CI = 1.569-11.55, P = 0.004 Allelic)] with RA susceptibility in the local population. The former single nucleotide mutation was nonfunctional, whereas the latter, residing in the exonic region of a linkage-proven autoimmunity gene PTPN22, was involved in R620→W620 substitution. Comparative molecular dynamic simulations and free-energy calculations revealed a radical impact on the geometry/confirmation of key functional moieties in the mutant protein leading to a rather weak binding of W620 variant with the interacting receptor (SRC kinase). The interaction imbalance and binding instabilities provide convincing clues about the insufficient inhibition of T cell activation and/or ineffective clearance of autoimmune clones - a hallmark of several autoimmune disorders. In conclusion, the present research describes the association of 2 hotspot mutations in IL-4 promoter and PTPN22 gene with RA susceptibility in the Pakistani study cohort. It also details how a functional mutation in PTPN22 impacts the overall protein geometry, charge, and/or receptor interactions to contribute to RA susceptibility.