A Comprehensive In Silico Analysis of the Functional and Structural Impact of Nonsynonymous SNPs in the ABCA1 Transporter Gene

Impact of highly deleterious non-synonymous polymorphisms on GRIN2A protein's structure and function

GRIN2A is a gene that encodes NMDA receptors found in the central nervous system and plays a pivotal role in excitatory synaptic transmission, plasticity and excitotoxicity in the mammalian central nervous system. Changes in this gene have been associated with a spectrum of neurodevelopmental disorders such as epilepsy. Previous studies on GRIN2A suggest that non-synonymous single nucleotide polymorphisms (nsSNPs) can alter the protein's structure and function. To gain a better understanding of the impact of potentially deleterious variants of GRIN2A, a range of bioinformatics tools were employed in this study. Out of 1320 nsSNPs retrieved from the NCBI database, initially 16 were predicted as deleterious by 9 tools. Further assessment of their domain association, conservation profile, homology models, interatomic interaction, and Molecular Dynamic Simulation revealed that the variant I463S is likely to be the most deleterious for the structure and function of the protein. Despite the limitations of computational algorithms, our analyses have provided insights that can be a valuable resource for further in vitro and in vivo research on GRIN2A-associated diseases.

Insilico prediction and functional analysis of nonsynonymous SNPs in human CTLA4 gene

The CTLA4 receptor is an immune checkpoint involved in the downregulation of T cells. Polymorphisms in this gene have been found to be associated with different diseases like rheumatoid arthritis, autosomal dominant immune dysregulation syndrome, juvenile idiopathic arthritis and autoimmune Addison's disease. Therefore, the identification of polymorphisms that have an effect on the structure and function of CTLA4 gene is important. Here we identified the most damaging missense or non-synonymous SNPs (nsSNPs) that might be crucial for the structure and function of CTLA4 using different bioinformatics tools. These in silico tools included SIFT, PROVEAN, PhD-SNP, PolyPhen-2 followed by MutPred2, I-Mutant 2.0 and ConSurf. The protein structures were predicted using Phyre2 and I-TASSER, while the gene-gene interactions were predicted by GeneMANIA and STRING. Our study identified three damaging missense SNPs rs1553657429, rs1559591863 and rs778534474 in coding region of CTLA4 gene. Among these SNPs the rs1553657429 showed a loss of potential phosphorylation site and was found to be highly conserved. The prediction of gene-gene interaction showed the interaction of CTlA4 with other genes and its importance in different pathways. This investigation of damaging nsSNPs can be considered in future while studying CTLA4 related diseases and can be of great importance in precision medicine.

A comprehensive in silico investigation into the nsSNPs of Drd2 gene predicts significant functional consequences in dopamine signaling and pharmacotherapy

DRD2 is a neuronal cell surface protein involved in brain development and function. Variations in the Drd2 gene have clinical significance since DRD2 is a pharmacotherapeutic target for treating psychiatric disorders like ADHD and schizophrenia. Despite numerous studies on the disease association of single nucleotide polymorphisms (SNPs) in the intronic regions, investigation into the coding regions is surprisingly limited. In this study, we aimed at identifying potential functionally and pharmaco-therapeutically deleterious non-synonymous SNPs of Drd2. A wide array of bioinformatics tools was used to evaluate the impact of nsSNPs on protein structure and functionality. Out of 260 nsSNPs retrieved from the dbSNP database, initially 9 were predicted as deleterious by 15 tools. Upon further assessment of their domain association, conservation profile, homology models and inter-atomic interaction, the mutant F389V was considered as the most impactful. In-depth analysis of F389V through Molecular Docking and Dynamics Simulation revealed a decline in affinity for its native agonist dopamine and an increase in affinity for the antipsychotic drug risperidone. Remarkable alterations in binding interactions and stability of the protein-ligand complex in simulated physiological conditions were also noted. These findings will improve our understanding of the consequence of nsSNPs in disease-susceptibility and therapeutic efficacy.

In Silico Tools and Approaches for the Prediction of Functional and Structural Effects of Single-Nucleotide Polymorphisms on Proteins: An Expert Review

Single-nucleotide polymorphisms (SNPs) are single-base variants that contribute to human biological variation and pathogenesis of many human diseases. Among all SNP types, nonsynonymous single-nucleotide polymorphisms (nsSNPs) can alter many structural, biochemical, and functional features of a protein such as folding characteristics, charge distribution, stability, dynamics, and interactions with other proteins/nucleotides. These modifications in the protein structure can lead nsSNPs to be closely associated with many multifactorial diseases such as cancer, diabetes, and neurodegenerative diseases. Predicting structural and functional effects of nsSNPs with experimental approaches can be time-consuming and costly; hence, computational prediction tools and algorithms are being widely and increasingly utilized in biology and medical research. This expert review examines the in silico tools and algorithms for the prediction of functional or structural effects of SNP variants, in addition to the description of the phenotypic effects of nsSNPs on protein structure, association between pathogenicity of variants, and functional or structural features of disease-associated variants. Finally, case studies investigating the functional and structural effects of nsSNPs on selected protein structures are highlighted. We conclude that creating a consistent workflow with a combination of in silico approaches or tools should be considered to increase the performance, accuracy, and precision of the biological and clinical predictions made in silico.

Computational SNP Analysis and Molecular Simulation Revealed the Most Deleterious Missense Variants in the NBD1 Domain of Human ABCA1 Transporter

The ATP-binding cassette transporter A1 (ABCA1) is a membrane-bound exporter protein involved in regulating serum HDL level by exporting cholesterol and phospholipids to load up in lipid-poor ApoA-I and ApoE, which allows the formation of nascent HDL. Mutations in the ABCA1 gene, when presents in both alleles, disrupt the canonical function of ABCA1, which associates with many disorders related to lipid transport. Although many studies have reported the phenotypic effects of a large number of ABCA1 variants, the pathological effect of non-synonymous polymorphisms (nsSNPs) in ABCA1 remains elusive. Therefore, aiming at exploring the structural and functional consequences of nsSNPs in ABCA1, in this study, we employed an integrated computational approach consisting of nine well-known in silico tools to identify damaging SNPs and molecular dynamics (MD) simulation to get insights into the magnitudes of the damaging effects. In silico tools revealed four nsSNPs as being most deleterious, where the two SNPs (G1050V and S1067C) are identified as the highly conserved and functional disrupting mutations located in the NBD1 domain. MD simulation suggested that both SNPs, G1050V and S1067C, changed the overall structural flexibility and dynamics of NBD1, and induced substantial alteration in the structural organization of ATP binding site. Taken together, these findings direct future studies to get more insights into the role of these variants in the loss of the ABCA1 function.

In silico analysis on the functional and structural impact of Rad50 mutations involved in DNA strand break repair

BACKGROUND

DNA double strand break repair is important to preserve the fidelity of our genetic makeup after DNA damage. Rad50 is one of the components in MRN complex important for DNA repair mechanism. Rad50 mutations can lead to microcephaly, mental retardation and growth retardation in human. However, Rad50 mutations in human and other organisms have never been gathered and heuristically compared for their deleterious effects. It is important to assess the conserved region in Rad50 and its homolog to identify vital mutations that can affect functions of the protein.

METHOD

In this study, Rad50 mutations were retrieved from SNPeffect 4.0 database and literature. Each of the mutations was analyzed using various bioinformatic analyses such as PredictSNP, MutPred, SNPeffect 4.0, I-Mutant and MuPro to identify its impact on molecular mechanism, biological function and protein stability, respectively.

RESULTS

We identified 103 mostly occurred mutations in the Rad50 protein domains and motifs, which only 42 mutations were classified as most deleterious. These mutations are mainly situated at the specific motifs such as Walker A, Q-loop, Walker B, D-loop and signature motif of the Rad50 protein. Some of these mutations were predicted to negatively affect several important functional sites that play important roles in DNA repair mechanism and cell cycle signaling pathway, highlighting Rad50 crucial role in this process. Interestingly, mutations located at non-conserved regions were predicted to have neutral/non-damaging effects, in contrast with previous experimental studies that showed deleterious effects. This suggests that software used in this study may have limitations in predicting mutations in non-conserved regions, implying further improvement in their algorithm is needed. In conclusion, this study reveals the priority of acid substitution associated with the genetic disorders. This finding highlights the vital roles of certain residues such as K42E, C681A/S, CC684R/S, S1202R, E1232Q and D1238N/A located in Rad50 conserved regions, which can be considered for a more targeted future studies.

In silico analysis of nonsynonymous single-nucleotide polymorphisms (nsSNPs) of the SMPX gene

Mutations in the SMPX gene can disrupt the regular activity of the SMPX protein, which is involved in the hearing process. Recent reports showing a link between nonsynonymous single-nucleotide polymorphisms (nsSNPs) in SMPX and hearing loss, thus classifying deleterious SNPs in SMPX will be an uphill task before designing a more extensive population study. In this study, damaging nsSNPs of SMPX from the dbSNP database were identified by using 13 bioinformatics tools. Initially, the impact of nsSNPs in the SMPX gene were evaluated through different in silico predictors; and the deleterious convergent changes were analyzed by energy-minimization-guided residual network analysis. In addition, the pathogenic effects of mutations in SMPX-mediated protein-protein interactions were also characterized by structural modeling and binding energy calculations. A total of four mutations (N19D, A29T, K54N, and S71L) were found to be highly deleterious by all the tools, which are located at highly conserved regions. Furthermore, all four mutants showed structural alterations, and the communities of amino acids for mutant proteins were readily changed, compared to the wild-type. Among them, A29T (rs772775896) was revealed as the most damaging nsSNP, which caused significant structural deviation of the SMPX protein, as a result reducing the binding affinity to other functional partners. These findings reflect the computational insights into the deleterious role of nsSNPs in SMPX, which might be helpful for subjecting wet-lab confirmatory analysis.

Identification and in silico analysis of functional SNPs of human TAGAP protein: A comprehensive study

Genetic polymorphisms in TAGAP gene have been associated with many diseases including rheumatoid arthritis, multiple sclerosis and other autoimmune disorders. Identifying functional SNPs in such disease associated genes is an uphill task hence before planning larger population study, it is better to scrutinize putative functional SNPs. In this study we used various computational approaches to identify nsSNPs which are deleterious to the structure and/or function of TAGAP protein that might be causing these diseases. Computational analysis was performed by five different in silico tools including SIFT, PROVEAN, PolyPhen-2, PhD-SNP and SNPs&GO. The study concludes that mutations of Glycine → Glutamic Acid at position 120, Glycine → Tryptophan at position 141 and Valine → Methionine at position 151 are major mutations in native TAGAP protein which might contribute to its malfunction and ultimately causing disease. The study also proposed 3D structures of native TAGAP protein and its three mutants. Future studies should consider these nsSNPs as main target mutations in various diseases involving TAGAP malfunction. This is the first comprehensive study, where TAGAP gene variants were analyzed using in silico tools hence will be of great help while considering large scale studies and also in developing precision medicines for cure of diseases related to these polymorphisms. Furthermore, animal models of various autoimmune diseases and having these mutations might be of help in exploring their precise roles.

Functional and Structural Impact of ATP-Binding Cassette Transporter A1 R219K and I883M Gene Polymorphisms in Obese Children and Adolescents

INTRODUCTION

Obesity is a serious medical condition that affects children and adolescents. ATP-binding cassette transporter A1 (ABCA1) protein is known to mediate the transport of intracellular cholesterol and phospholipids across the cell membranes. Thus, we aimed to investigate the association between ABCA1 gene polymorphisms and overweight/obesity risk, and to evaluate their relation to the lipid profile.

MATERIALS AND METHODS

The study included in silico analysis of ABCA1 gene and protein. Two genetic variants in ABCA1 gene-R219K (rs2230806; G/A) and I883M (rs2066714; A/G)-were genotyped in 128 normal weight and 128 overweight/obese subjects using polymerase chain reaction-restriction fragment length polymorphism technology. Anthropometric and biochemical assessments were performed.

RESULTS

Our findings suggest that the heterozygote GA genotype of R219K polymorphism increased susceptibility to obesity under the heterozygous model (odds ratio 2.75, 95 % CI 1.01-6.12; p = 0.014) compared with the control group. This susceptibility could be gender-specific, with higher risk among females. In addition, the A variant was associated with a higher degree of obesity (p < 0.001). On the other hand, individuals with the G variant of I883M polymorphism showed lower susceptibility to obesity under all genetic models (allelic, homozygote, heterozygote, dominant, and recessive models; p < 0.05), with no observed association with body mass index or degree of obesity. However, both single nucleotide polymorphisms (SNPs) showed significant differences in lipid levels among patients with different genotypes.

CONCLUSIONS

The study results suggest that R219K and I883M SNPs of the ABCA1 gene may play a role in susceptibility to obesity in our Egyptian population; the former increases susceptibility and phenotype severity, and the latter is protective. Larger epidemiological studies are needed for validation of the results.

The Bioinformatics Report of Mutation Outcome on NADPH Flavin Oxidoreductase Protein Sequence in Clinical Isolates of H. pylori

frxA gene has been implicated in the metronidazole nitro reduction by H. pylori. Alternatively, frxA is expected to contribute to the protection of urease and to the in vivo survival of H. pylori. The aim of present study is to report the mutation effects on the frxA protein sequence in clinical isolates of H. pylori in our community. Metronidazole resistance was proven in 27 of 48 isolates. glmM and frxA genes were used for molecular confirmation of H. pylori isolates. The primer set for detection of whole sequence of frxA gene for the effect of mutation on protein sequence was used. DNA and protein sequence evaluation and analysis were done by blast, Clustal Omega, and T COFFEE programs. Then, FrxA protein sequences from six metronidazole-resistant clinical isolates were analyzed by web-based bioinformatics tools. The result of six metronidazole-resistant clinical isolates in comparison with strain 26695 showed ten missense mutations. The result with the STRING program revealed that no change was seen after alterations in these sequences. According to consensus data involving four methods, residue substitutions at 40, 13, and 141 increase the stability of protein sequence after mutation, while other alterations decrease. Residue substitutions at 40, 43, 141, 138, 169, and 179 are deleterious, while, V7I, Q10R, V34I, and V96I alterations are neutral. As FrxA contribute to survival of bacterium and in regard to the effect of mutations on protein function, it might affect the survival and bacterium phenotype and it need to be studied more. Also, none of the stability prediction tool is perfect; iStable is the best predictor method among all methods.