Abstract
Rheumatoid arthritis (RA) is a chronic autoimmune disorder characterized by inflammation and joint destruction. Various genetic studies have identified 100 susceptible genes linked with RA. Among them, four SNPs (rs6671847, rs1801274, rs17400517, and rs6668534) in the FCGR2A gene have been shown to be associated with RA susceptibility. The functional consequences of these pathogenic variants on the protein structure, stability, and interaction with other molecules remain unknown. In this study, we used multiple in silico methods to examine the effects of these mutations on the FCGR2A gene. We annotated the gene structurally, functionally, and Predicted physiochemical characterization, gene expression profiling, and Post-translational modification (PTM) analysis. We found that the (rs1801274 H167R) mutation in FCGR2A had a significant effect on the protein structure and interactions, whereas the other three SNPs (rs6671847, rs17400517, and rs6668534) had no significant impact on protein function. We also performed protein-protein interaction and molecular docking studies to identify potential therapeutic targets for RA. Our structure-based analysis showed that these mutations affect the protein core regions and functional domains, leading to altered protein-protein interactions both directly and indirectly. We identified three novel potential target sites for the development of drugs to treat RA by developing a drugability profile. These findings may have significant implications for the pathogenesis of RA and can be validated further through in vitro studies. Our study demonstrates the potential utility of multidirectional computational analysis for screening clinically relevant mutations in RA before undertaking biological assays. The findings highlight the importance of understanding the structural and functional consequences of pathogenic variants for developing effective therapeutic strategies. Overall, the insights gained from our in-silico approaches could aid in the development of novel therapies for RA and could pave the way for modified medicine in the treatment of this complex disease.
References
Halushka MK, Fan J-B, Bentley K, Hsie L, Shen N, Weder A, et al. Patterns of single-nucleotide polymorphisms in candidate genes for blood-pressure homeostasis. Nature genetics. 1999;22(3):239-47.
Angelotti F, Parma A, Cafaro G, Capecchi R, Alunno A, Puxeddu I. One year in review 2017: pathogenesis of rheumatoid arthritis. Clin Exp Rheumatol. 2017;35(3):368-78.
Yau AC, Holmdahl R. Rheumatoid arthritis: identifying and characterising polymorphisms using rat models. Disease models & mechanisms. 2016;9(10):1111-23.
Van der Woude D, Houwing‐Duistermaat JJ, Toes RE, Huizinga TW, Thomson W, Worthington J, et al. Quantitative heritability of anti–citrullinated protein antibody–positive and anti–citrullinated protein antibody–negative rheumatoid arthritis. Arthritis & Rheumatism: Official Journal of the American College of Rheumatology. 2009;60(4):916-23.
Sode J, Vogel U, Bank S, Andersen PS, Hetland ML, Locht H, et al. Genetic variations in pattern recognition receptor loci are associated with anti-TNF response in patients with rheumatoid arthritis. PloS one. 2015;10(10):e0139781.
Lenz TL, Deutsch AJ, Han B, Hu X, Okada Y, Eyre S, et al. Widespread non-additive and interaction effects within HLA loci modulate the risk of autoimmune diseases. Nature genetics. 2015;47(9):1085-90.
Stahl EA, Wegmann D, Trynka G, Gutierrez-Achury J, Do R, Voight BF, et al. Bayesian inference analyses of the polygenic architecture of rheumatoid arthritis. Nature genetics. 2012;44(5):483-9.
Kim K, Bang S-Y, Lee H-S, Bae S-C. Update on the genetic architecture of rheumatoid arthritis. Nature Reviews Rheumatology. 2017;13(1):13-24.
Pandey JP. Genomewide association studies and assessment of risk of disease. The New England journal of medicine. 2010;363(21):2076-7.
Ravetch JV, Bolland S. Igg fc receptors. Annual review of immunology. 2001;19(1):275-90.
Yang Y, Peng L, He C, Xing S, Li D, Jin T, et al. The role of genetic variants in FCGR2A on the risk of rheumatoid arthritis in the Han Chinese population. 2020.
McWilliam H, Li W, Uludag M, Squizzato S, Park YM, Buso N, et al. Analysis tool web services from the EMBL-EBI. Nucleic acids research. 2013;41(W1):W597-W600.
Garg VK, Avashthi H, Tiwari A, Jain PA, Ramkete PW, Kayastha AM, et al. MFPPI–multi FASTA ProtParam interface. Bioinformation. 2016;12(2):74.
Julenius K. NetCGlyc 1.0: prediction of mammalian C-mannosylation sites. Glycobiology. 2007;17(8):868-76.
Kiemer L, Lund O, Brunak S, Blom N. Coronavirus 3CL pro proteinase cleavage sites: Possible relevance to SARS virus pathology. BMC bioinformatics. 2004;5(1):1-9.
Johansen MB, Kiemer L, Brunak S. Analysis and prediction of mammalian protein glycation. Glycobiology. 2006;16(9):844-53.
Gupta R, Jung E, Brunak S. Prediction of N-glycosylation sites in human proteins. 2004.
Steentoft C, Vakhrushev SY, Joshi HJ, Kong Y, Vester‐Christensen MB, Schjoldager KTB, et al. Precision mapping of the human O‐GalNAc glycoproteome through SimpleCell technology. The EMBO journal. 2013;32(10):1478-88.
Blom N, Sicheritz‐Pontén T, Gupta R, Gammeltoft S, Brunak S. Prediction of post‐translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics. 2004;4(6):1633-49.
Duckert P, Brunak S, Blom N. Prediction of proprotein convertase cleavage sites. Protein Engineering Design and Selection. 2004;17(1):107-12.
Guex N, Peitsch MC. SWISS‐MODEL and the Swiss‐Pdb Viewer: an environment for comparative protein modeling. electrophoresis. 1997;18(15):2714-23.
Lovell SC, Davis IW, Arendall III WB, De Bakker PI, Word JM, Prisant MG, et al. Structure validation by Cα geometry: ϕ, ψ and Cβ deviation. Proteins: Structure, Function, and Bioinformatics. 2003;50(3):437-50.
Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, et al. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic acids research. 2010;39(suppl_1):D561-D8.
Bairoch A, Apweiler R. The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic acids research. 2000;28(1):45-8.
Perez-Vilar J, Randell SH, Boucher RC. C-Mannosylation of MUC5AC and MUC5B Cys subdomains. Glycobiology. 2004;14(4):325-37.
Hou L, Zhao H. A review of post-GWAS prioritization approaches. Frontiers in genetics. 2013;4:280.
Cooper GM, Shendure J. Needles in stacks of needles: finding disease-causal variants in a wealth of genomic data. Nature Reviews Genetics. 2011;12(9):628-40.
Al-Abbasi FA, Mohammed K, Sadath S, Banaganapalli B, Nasser K, Shaik NA. Computational protein phenotype characterization of IL10RA mutations causative to early onset inflammatory bowel disease (IBD). Frontiers in genetics. 2018;9:146.
Khan S, Vihinen M. Spectrum of disease-causing mutations in protein secondary structures. BMC Structural Biology. 2007;7(1):1-18.
Miao Z, Cao Y. Quantifying side-chain conformational variations in protein structure. Scientific reports. 2016;6(1):1-10.
Hijikata A, Tsuji T, Shionyu M, Shirai T. Decoding disease-causing mechanisms of missense mutations from supramolecular structures. Scientific reports. 2017;7(1):1-8.
Rao VS, Srinivas K, Sujini G, Kumar G. Protein-protein interaction detection: methods and analysis. International journal of proteomics. 2014;2014.
Van Zundert G, Rodrigues J, Trellet M, Schmitz C, Kastritis P, Karaca E, et al. The HADDOCK2. 2 web server: user-friendly integrative modeling of biomolecular complexes. Journal of molecular biology. 2016;428(4):720-5.

This work is licensed under a Creative Commons Attribution 4.0 International License.
Copyright (c) 2023 Pak-Euro Journal of Medical and Life Sciences