Computational Molecular Phenotypic Analysis of PTPN22 (W620R), IL6R (D358A), and TYK2 (P1104A) Gene Mutations of Rheumatoid Arthritis

Abstract

Rheumatoid arthritis (RA) is a chronic autoimmune disorder of bone joints caused by the complex interplay between several factors like body physiology, the environment with genetic background. The recent meta-analysis of GWAS has expanded the total number of RA-associated loci to more than 100, out of which approximately ∼97% (98 variants) loci are located in non-coding regions, and the other ∼3% (3 variants) are in three different non-HLA genes, i.e., TYK2 (Prp1104Ala), IL6R (Asp358Ala), and PTPN22 (Trp620Arg). However, whether these variants prompt changes in the protein phenotype with regards to its stability, structure, and interaction with other molecules, remains unknown. Thus, we selected the three clinically pathogenic variants described above, as positive controls and applied diverse computational methods to scrutinize if those mutations cause changes in the protein phenotype. Both wild type and mutant protein structures of PTPN22 (W620R), IL6R (D358A), and TYK2 (P1104A) were modeled and studied for structural deviations. Furthermore, we have also studied the secondary structure characteristics, solvent accessibility and stability, and the molecular interaction deformities caused by the amino acid substitutions. We observed that simple nucleotide predictions of SIFT, PolyPhen, CADD and FATHMM yields mixed findings in screening the RA-missense variants which showed a ≥P-value threshold of 5 × 10−8 in genome wide association studies. However, structure-based analysis confirms that mutant structures shows subtle but significant changes at their core regions, but their functional domains seems to lose wild type like functional interaction. Our findings suggest that the multidirectional computational analysis of clinically potential RA-mutations could act as a primary screening step before undertaking functional biology assays.

Whole exome sequencing identifies rare biallelic ALMS1 missense and stop gain mutations in familial Alström syndrome patients

“…ANNOVAR software was used to annotate novel variants. SIFT, PolyPhen-2 and Mutation Tester algorithms were utilized to predict functional consequences of non-synonymous mutations in the exons (Shaik and Banaganapalli, 2019, Shaik et al, 2018).…”

Molecular insights into the coding region mutations of low‐density lipoprotein receptor adaptor protein 1 (LDLRAP1) linked to familial hypercholesterolemia

“…Up to 80% of single gene disease causing missense mutations are assumed to alter protein stability. In the present study, the DUET webserver predicted that W22S and R39W mutations located upstream of PTB domain, the V92M and R151W mutations located within the PTB domain, and the M281T mutation located in the carboxy terminal end of the LDLRAP1 protein are most deleterious as a result of their negative free energy values . The destabilization could be a result of the different physicochemical properties between the amino acids, such as its hydrophobic nature .…”

Unraveling the role of salt-sensitivity genes in obesity with integrated network biology and co-expression analysis

“…The PIN created from Bisogenet was optimized and imported to Cytoscape 3.2.1 in order to represent and measure the different parameters like DC and BC connected to network centrality of each individual protein in the biological network [20]. The Network Analyzer [21] Cytoscape plugin was deployed to monitor the network's local and global centrality parameters [14,[22][23][24].…”