Identification and prediction of the consequences of nonsynonymous SNPs in glyceraldehyde 3-phosphate dehydrogenase (GAPDH) gene of zebrafish Danio rerio

Identification and prediction of the consequences of nonsynonymous SNPs in glyceraldehyde 3-phosphate dehydrogenase (GAPDH) gene of zebrafish Danio rerio

The recent progress in DNA sequencing and computational algorithms dramatically heightened the value of single nucleotide polymorphism (SNP) databases. The experimental studies of mutation for the particular gene of interest are laborious and time-consuming. In this study, we used several computational algorithms to investigate the structural and functional consequences of SNPs on the glyceraldehyde 3-phosphate dehydrogenase (GAPDH) gene of zebrafish (Danio rerio). The GAPDH gene is involved in multiple cellular activities, in addition to its association with classical glycolytic pathway. It has most commonly been used as the housekeeping gene for differential gene expression profiling, including western blots. The computational analysis, which used sequence-based tools such as SIFT, PANTHER, PROVEAN, and I-Mutant2.0, predicted deleterious nonsynonymous SNPs (nsSNPs) at p.W194R in the GAPDH gene. The SWISS Model tool generated a homology-modeled 3D structure of GAPDH by applying homologous protein structures retrieved from the RCSB PDB databank (PDB ID: 1U8F_O, 1ZNQ_O, 1J0X_O, 3GPD_R, and 1SZJ_G). The sequence identities of the above selected models were 86%, 86%, 86%, 83%, and 75%, respectively. The Procheck quality assessment in the Ramachandran plot showed that the incorporated mutation (p.W194R) led to a shifting of amino acid residues from the most favored regions towards the disallowed region, as compared to the native counterpart. The lowered ERRAT and PROSA Z scores in the mutant confirmed the overall deteriorating quality of mutant GAPDH. We also predicted an effect of mutation on increased H-bonding patterns due to mutant residue. The molecular docking analysis, while binding GAPDH to its ligand (NAD+), also revealed elevated global energy in the mutant. Altogether, this study demonstrated the impact of nsSNPs on the protein structure and function. The results of work provide an insight/roadmap for future investigations with regard to the pathological consequences of nsSNPs using in vivo methods

___

  • ReferencesAlanazi M, Abduljaleel Z, Khan W, Warsy AS, Elrobh M, Khan Z, Al Amri A, Bazzi MD (2011). In silico analysis of single nucleotide polymorphism (SNPs) in human β-globin gene. PLoS ONE 6: e25876.Barber RD, Harmer DW, Coleman RA, Clark BJ (2005). GAPDH as a housekeeping gene: analysis of GAPDH mRNA expression in a panel of 72 human tissues. Physiol Genomics 21: 389–395.Butterfield DA, Hardas SS, Lange ML (2010). Oxidatively modified glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and Alzheimer’s disease: many pathways to neurodegeneration. J Alzheimers Dis 20: 369–393.Capriotti E, Fariselli P, Casadio R (2005). I-Mutant2.0: predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res 33: W306–310.Chakravarti R, Aulak KS, Fox PL, Stuehr DJ (2010). GAPDH regulates cellular heme insertion into inducible nitric oxide synthase. P Natl Acad Sci USA 107: 18004–18009.Chen L, Li Q, Li L (2014). The modeled structures of Deg5 and Deg8 proteases in Arabidopsis thaliana. Turk J Biol 38: 168–176.Choi Y, Sims GE, Murphy S, Miller JR, Chan AP (2012). Predicting the functional effect of amino acid substitutions and indels. PLoS ONE 7: e46688.De Alencar SA, Lopes JC (2010). A comprehensive in silico analysis of the functional and structural impact of SNPs in the IGF1R gene. J Biomed Biotechnol 2010: 715139. De Cristofaro R, Carotti A, Akhavan S, Palla R, Peyvandi F, Altomare C, Mannucci PM (2006). The natural mutation by deletion of Lys9 in the thrombin A-chain affects the pKa value of catalytic residues, the overall enzyme’s stability and conformational transitions linked to Na+ binding. FEBS J 273: 159–169.Eshkoor SA, Ismail P, Rahman SA (2014). Does CYP1A1 gene polymorphism affect cell damage biomarkers and ageing? Turk J Biol 38: 219–225.Filiz E, Tombuloğlu H (2015). Genome-wide distribution of superoxide dismutase (SOD) gene families in Sorghum bicolor. Turk J Biol 39: 49–59.George Priya Doss C, Rajasekaran R, Sudandiradoss C, Ramanathan K, Purohit R, Sethumadhavan R (2008). A novel computational and structural analysis of nsSNPs in CFTR gene. Genomic Med 2: 23–32.George Priya Doss C, Rajith B (2012). Computational refinement of functional single nucleotide polymorphisms associated with AT M gene. PLoS ONE 7: e34573.Gökgöz NB, Yalaz S, Avcı NG, Buldum G, Ölmez EO, Akbulut BS (2015). Investigation of the in vivo interaction between β-lactamase and its inhibitor protein. Turk J Biol 39: 485–492.Han JH, Kerrison N, Chothia C, Teichmann SA (2006). Divergence of interdomain geometry in two-domain proteins. Structure 14: 935–945.Hardt M, Laine RA (2004). Mutation of active site residues in the chitin-binding domain ChBDChiA1 from chitinase A1 of Bacillus circulans alters substrate specificity: use of a green fluorescent protein binding assay. Arch Biochem Biophys 426: 286–297.Honda JY, Kobayashi I, Kiyoshima T, Yamaza H, Xie M, Takahashi K, Enoki N, Nagata K, Nakashima A, Sakai H (2008). Glycolytic enzyme Pgk1 is strongly expressed in the developing tooth germ of the mouse lower first molar. Histol Histopathol 23: 423–432.Hussain MRM, Shaik NA, Al-Aama JY, Asfour HZ, Khan FS, Masoodi TA, Khan MA, Shaik NS (2012). In silico analysis of Single Nucleotide Polymorphisms (SNPs) in human BRAFgene. Gene 508: 188–196.Johnson MM, Houck J, Chen C (2005). Screening for deleterious nonsynonymous single-nucleotide polymorphisms in genes involved in steroid hormone metabolism and response. Cancer Epidem Biomar 14: 1326–1329.Jones R, Ruas M, Gregory F, Moulin S, Delia D, Manoukian S, Rowe J, Brookes S, Peters G (2007). A CDKN2A mutation in familial melanoma that abrogates binding of p16INK4a to CDK4 but not CDK6. Cancer Res 67: 9134–9141.Kamaraj B, Purohit R (2013). In silico screening and molecular dynamics simulation of disease-associated nsSNP in TYRP1 gene and its structural consequences in OCA3. Biomed Res Int 2013: 697051.Khan S, Vihinen M (2010). Performance of protein stability predictors. Hum Mutat 31: 675–684.Kumar P, Henikoff S, Ng PC (2009). Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc 4: 1073–1081.Luduena RF (1993). Are tubulin isotypes functionally significant? Mol Biol Cell 4: 445–457.Mazzola JL, Sirover MA (2001). Reduction of glyceraldehyde-3-phosphate dehydrogenase activity in Alzheimer’s disease and in Huntington’s disease fibroblasts. J Neurochem 76: 442–449.Mohamoud HS, Hussain MR, El-Harouni AA, Shaik NA, Qasmi ZU, Merican AF, Baig M, Anwar Y, Asfour H, Bondagji N et al. (2014). First comprehensive in silico analysis of the functional and structural consequences of SNPs in human GalNAc-T1gene. Comput Math Method M 2014: 904052.Mohapatra C, Barman HK, Panda RP, Kumar S, Das V, Mohanta R, Mohapatra SD, Jayasankar P (2010). Cloning of cDNA and prediction of peptide structure of Plzf expressed in the spermatogonial cells of Labeo rohita. Mar Genomics 3: 157–163.Ng PC, Henikoff S (2003). SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res 31: 3812–3814.Ng PC, Henikoff S (2006). Predicting the effects of amino acid substitutions on protein function. Annu Rev Genom Hum G 7: 61–80.Ng SY, Gunning P, Eddy R, Ponte P, Leavitt J, Shows T, Kedes L (1985). Evolution of the functional human beta-actin gene and its multi-pseudogene family: conservation of noncoding regions and chromosomal dispersion of pseudogenes. Mol Cell Biol 5: 2720–2732.Ode H, Matsuyama S, Hata M, Neya S, Kakizawa J, Sugiura W, Hoshino T (2007). Computational characterization of structural role of the non-active site mutation M36I of human immunodeficiency virus type 1 protease. J Mol Biol 370: 598–607.Robbins AR, Ward RD, Oliver C (1995). A mutation in glyceraldehyde 3-phosphate dehydrogenase alters endocytosis in CHO cells. J Cell Biol 130: 1093–1104.Sarropoulou E, Nousdili D, Kotoulas G, Magoulas A (2011). Functional divergences of GAPDH isoforms during early development in two perciform fish species. Mar Biotechnol 13: 1115–1124.Schaefer C, Rost B (2012). Predict impact of single amino acid change upon protein structure. BMC Genomics 13: S4.Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, Diemer K, Muruganujan A, Narechania A (2003). PANTHER: a library of protein families and subfamilies indexed by function. Genome Res 13: 2129–2141.Thomas PD, Kejariwal A (2004). Coding single-nucleotide polymorphisms associated with complex vs. Mendelian disease: evolutionary evidence for differences in molecular effects. P Natl Acad Sci USA 101: 15398–15403.Tisdale EJ, Artalejo CR (2007). A GAPDH mutant defective in Src-dependent tyrosine phosphorylation impedes Rab2-mediated events. Traffic 8: 733–741.Wang Z, Gerstein M, Snyder M (2009). RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 10: 57–63.Wang Z, Moult J (2001). SNPs, protein structure, and disease. Hum Mutat 17: 263–270.Wiederstein M, Sippl MJ (2007). ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res 35: W407–410.