A novel computational and structural analysis of nsSNPs in CFTR gene

Computational biophysical characterization of the effect of gatekeeper mutations on the binding of ponatinib to the FGFR kinase

Fibroblast Growth Factor Receptor (FGFR) is connected to numerous downstream signalling cascades regulating cellular behavior. Any dysregulation leads to a plethora of illnesses, including cancer. Therapeutics are available, but drug resistance driven by gatekeeper mutation impedes the treatment. Ponatinib is an FDA-approved drug against BCR-ABL kinase and has shown effective results against FGFR-mediated carcinogenesis. Herein, we undertake molecular dynamics simulation-based analysis on ponatinib against all the FGFR isoforms having Val to Met gatekeeper mutations. The results suggest that ponatinib is a potent and selective inhibitor for FGFR1, FGFR2, and FGFR4 gatekeeper mutations. The extensive electrostatic and van der Waals interaction network accounts for its high potency. The FGFR3_VM mutation has shown resistance towards ponatinib, which is supported by their lesser binding affinity than wild-type complexes. The disengaged molecular brake and engaged hydrophobic spine were believed to be the driving factors for weak protein-ligand interaction. Taken together, the inhibitory and structural characteristics exhibited by ponatinib may aid in thwarting resistance based on Val-to-Met gatekeeper mutations at an earlier stage of treatment and advance the design and development of other inhibitors targeted at FGFRs harboring gatekeeper mutations.

Identification and In-Silico study of non-synonymous functional SNPs in the human SCN9A gene

Single nucleotide polymorphisms are the most common form of DNA alterations at the level of a single nucleotide in the genomic sequence. Genome-wide association studies (GWAS) were carried to identify potential risk genes or genomic regions by screening for SNPs associated with disease. Recent studies have shown that SCN9A comprises the NaV1.7 subunit, Na+ channels have a gene encoding of 1988 amino acids arranged into 4 domains, all with 6 transmembrane regions, and are mainly found in dorsal root ganglion (DRG) neurons and sympathetic ganglion neurons. Multiple forms of acute hypersensitivity conditions, such as primary erythermalgia, congenital analgesia, and paroxysmal pain syndrome have been linked to polymorphisms in the SCN9A gene. Under this study, we utilized a variety of computational tools to explore out nsSNPs that are potentially damaging to heath by modifying the structure or activity of the SCN9A protein. Over 14 potentially damaging and disease-causing nsSNPs (E1889D, L1802P, F1782V, D1778N, C1370Y, V1311M, Y1248H, F1237L, M936V, I929T, V877E, D743Y, C710W, D623H) were identified by a variety of algorithms, including SNPnexus, SNAP-2, PANTHER, PhD-SNP, SNP & GO, I-Mutant, and ConSurf. Homology modeling, structure validation, and protein-ligand interactions also were performed to confirm 5 notable substitutions (L1802P, F1782V, D1778N, V1311M, and M936V). Such nsSNPs may become the center of further studies into a variety of disorders brought by SCN9A dysfunction. Using in-silico strategies for assessing SCN9A genetic variations will aid in organizing large-scale investigations and developing targeted therapeutics for disorders linked to these variations.

Investigation of pathogenic germline variants in gastric cancer and development of "GasCanBase" database

BACKGROUND

Gastric cancer, which is also known as stomach cancer, can be influenced by both germline and somatic mutations. Non-synonymous Single Nucleotide Polymorphisms (nsSNPs) in germline have long been reported to play a pivotal role in cancer progression.

AIM

The aim of this study is to examine the nsSNP in GC-associated genes. The study also aims to develop a database with extensive information regarding the nsSNPs in the GC-associated genes and their impacts.

METHODS AND RESULTS

A total of 34,588 nsSNPs from 1,493,460 SNPs of the 40 genes were extracted from the available SNP database. Drug binding and energy minimization were examined by molecular docking and YASARA. To validate the existence of the germline CDH1 gene mutation (rs34466743) in the isolated blood DNA of gastric cancer (GC) patients, polymerase chain reaction (PCR) and DNA sequencing were performed. According to the results of the gene network analysis, 17 genes may interact with other types of cancer. A total of 11,363 nsSNPs were detected within the 40 GC genes. Among these, 474 nsSNPs were predicted to be damaging and 40 to be the most damaging. The SNPs in domain regions were thought to be strong candidates that alter protein functions. Our findings proposed that most of the selected nsSNPs were within the domains or motif regions. Free Energy Deviation calculation of protein structure pointed toward noteworthy changes in the structure of each protein that can demolish its natural function. Subsequently, drug binding confirmed the structural variation and the ineffectiveness of the drug against the mutant model in individuals with these germline variants. Furthermore, in vitro analysis of the rs34466743 germline variant from the CDH1 gene confirmed the strength and robustness of the pipeline that could expand the somatic alteration for causing cancer. In addition, a comprehensive gastric cancer polymorphism database named "GasCanBase" was developed to make data available to researchers.

CONCLUSION

The findings of this study and the "GasCanBase" database may greatly contribute to our understanding of molecular epidemiology and the development of precise therapeutics for gastric cancer. GasCanBase is available at: https://www.gascanbase.com/.

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.

In-silico assessment of high-risk non-synonymous SNPs in ADAMTS3 gene associated with Hennekam syndrome and their impact on protein stability and function

Hennekam Lymphangiectasia-Lymphedema Syndrome 3 (HKLLS3) is a rare genetical disorder caused by mutations in a few genes including ADAMTS3. It is characterized by lymphatic dysplasia, intestinal lymphangiectasia, severe lymphedema and distinctive facial appearance. Up till now, no extensive studies have been conducted to elucidate the mechanism of the disease caused by various mutations. As a preliminary investigation of HKLLS3, we sorted out the most deleterious nonsynonymous single nucleotide polymorphisms (nsSNPs) that might affect the structure and function of ADAMTS3 protein by using a variety of in silico tools. A total of 919 nsSNPs in the ADAMTS3 gene were identified. 50 nsSNPs were predicted to be deleterious by multiple computational tools. 5 nsSNPs (G298R, C567Y, A370T, C567R and G374S) were found to be the most dangerous and can be associated with the disease as predicted by different bioinformatics tools. Modelling of the protein shows it can be divided into segments 1, 2 and 3, which are connected by short loops. Segment 3 mainly consists of loops without substantial secondary structures. With prediction tools and molecular dynamics simulation, some SNPs were found to significantly destabilize the protein structure and disrupt the secondary structures, especially in segment 2. The deleterious effects of mutations in segment 1 are possibly not from destabilization but from other factors such as the change in phosphorylation as suggested by post-translational modification (PTM) studies. This is the first-ever study of ADAMTS3 gene polymorphism, and the predicted nsSNPs in ADAMST3, some of which have not been reported yet in patients, will serve for diagnostic purposes and further therapeutic implications in Hennekam syndrome, contributing to better diagnosis and treatment.

In Silico Analysis of nsSNPs of Human KRAS Gene and Protein Modeling Using Bioinformatic Tools

The KRAS gene belongs to the RAS family and codes for 188 amino acid residues of KRAS protein, with a molecular mass of 21.6 kD. Non-synonymous single-nucleotide polymorphisms (nsSNPs) have been identified within the coding region in which some are associated with different diseases. However, structural changes are not well defined yet. In this study, we first categorized SNPs in the KRAS coding area and then used computational methods to determine their impact on the protein structure and stability. In addition, the three-dimensional model of KRAS was taken from the Protein Data Bank for structural modeling. Furthermore, genomic data were extracted from a variety of sources, including the 1000 Genome Project, dbSNPs, and ENSEMBLE, and assessed through in silico methods. Based on various tools used in this study, 10 out of 48 missense SNPs with rsIDs were found deleterious. The substitution of alanine for proline at position 146 pushed several residues toward the center of the protein. Arginine instead of leucine has a minor effect on protein structure and stability. In addition, the substitution of proline for leucine at the 34th position disrupted the structure and led to a bigger size than the wild-type protein, hence interrupting the protein interaction. Using the well-intended computational approach and applying several bioinformatic tools, we characterized and identified most damaging nsSNPs and further explored the structural dynamics and stability of KRAS protein.

Role of Mitochondrial Mutations in Ocular Aggregopathy

Background Mitochondria are essential cellular organelles that are responsible for oxidative stress-induced damage in age-dependent neurodegenerations such as glaucoma. Previous studies have linked mitochondrial DNA (mtDNA) mutations to cellular energy shortages that result in eye degeneration. Methodology To look for nucleotide variations in mtDNA in exfoliation syndrome/glaucoma (XFS/XFG), we performed a polymerase chain reaction (PCR) to amplify the entire coding region of the mitochondrial genome from peripheral blood of XFS/XFG (n = 25) patients and controls (n = 25). Results This study identified a total of 65 variations in XFS/XFG patients, of which 25 (38%) variations were non-synonymous single-nucleotide polymorphism (nsSNPs). Out of 25 nsSNPs, seven (five nsSNP in MT-ND4 and two in MT-ATP6 gene) were predicted as pathogenic using four different software, namely, SIFT, Polyphene2, mutation taster, and MutPred2. The pathogenic nsSNPs were then subjected to structural change analysis using online tools. Conclusions The pathogenic nsSNPs were found in both proteins' transmembrane domains and were expected to be conserved, but with lower protein stability (ΔΔG <- 0.5), indicating a possibly harmful effect in exfoliation. However, three-dimensional protein analysis indicated that the predicted mutations in MT-ND4 and MT-ATP6 were unlikely to alter the protein function.

Novel Disease-Associated Missense Single-Nucleotide Polymorphisms Variants Predication by Algorithms Tools and Molecular Dynamics Simulation of Human TCIRG1 Gene Causing Congenital Neutropenia and Osteopetrosis

T Cell Immune Regulator 1, ATPase H + Transporting V0 Subunit A3 (TCIRG1 gene provides instructions for making one part, the a3 subunit, of a large protein complex known as a vacuolar H + -ATPase (V-ATPase). V-ATPases are a group of similar complexes that act as pumps to move positively charged hydrogen atoms (protons) across membranes. Single amino acid changes in highly conserved areas of the TCIRG1 protein have been linked to autosomal recessive osteopetrosis and severe congenital neutropenia. We used multiple computational approaches to classify disease-prone single nucleotide polymorphisms (SNPs) in TCIRG1. We used molecular dynamics analysis to identify the deleterious nsSNPs, build mutant protein structures, and assess the impact of mutation. Our results show that fifteen nsSNPs (rs199902030, rs200149541, rs372499913, rs267605221, rs374941368, rs375717418, rs80008675, rs149792489, rs116675104, rs121908250, rs121908251, rs121908251, rs149792489 and rs116675104) variants are likely to be highly deleterious mutations as by incorporating them into wild protein they destabilize the wild protein structure and function. They are also located in the V-ATPase I domain, which may destabilize the structure and impair TCIRG1 protein activation, as well as reduce its ATPase effectiveness. These mutants have not yet been identified in patients suffering from CN and osteopetrosis while (G405R, R444L, and D517N) reported in our study are already associated with osteopetrosis. Mutation V52L reported in our study was identified in a patient suspected for CN. Finally, these mutants can help to further understand the broad pool of illness susceptibilities associated with TCIRG1 catalytic kinase domain activation and aid in the development of an effective treatment for associated diseases.

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.

Prediction of Deleterious Non-synonymous SNPs of Human STK11 Gene by Combining Algorithms, Molecular Docking, and Molecular Dynamics Simulation

Serine-threonine kinase11 (STK11) is a tumor suppressor gene which plays a key role in regulating cell growth and apoptosis. It is widely known as a multitasking kinase and engaged in cell polarity, cell cycle arrest, chromatin remodeling, energy metabolism, and Wnt signaling. The substitutions of single amino acids in highly conserved regions of the STK11 protein are associated with Peutz-Jeghers syndrome (PJS), which is an autosomal dominant inherited disorder. The abnormal function of the STK11 protein is still not well understood. In this study, we classified disease susceptible single nucleotide polymorphisms (SNPs) in STK11 by using different computational algorithms. We identified the deleterious nsSNPs, constructed mutant protein structures, and evaluated the impact of mutation by employing molecular docking and molecular dynamics analysis. Our results show that W239R and W308C variants are likely to be highly deleterious mutations found in the catalytic kinase domain, which may destabilize structure and disrupt the activation of the STK11 protein as well as reduce its catalytic efficiency. The W239R mutant is likely to have a greater impact on destabilizing the protein structure compared to the W308C mutant. In conclusion, these mutants can help to further realize the large pool of disease susceptibilities linked with catalytic kinase domain activation of STK11 and assist to develop an effective drug for associated diseases.

Computational analysis for the determination of deleterious nsSNPs in human MTHFD1 gene

Single nucleotide polymorphisms (SNPs) are the most common genetic polymorphisms and play a major role in many inherited diseases. Methylenetetrahydrofolate dehydrogenase 1 (MTHFD1) is one of the enzymes involved in folate metabolism. In the present study, the functional and structural consequences of nsSNPs of human MTHFD1 gene was analyzed using various computational tools like SIFT, PolyPhen2, PANTHER, PROVEAN, SNAP2, nsSNPAnalyzer, PhD-SNP, SNPs&GO, I-Mutant, MuPro, ConSurf, InterPro, NCBI Conserved Domain Search tool, ModPred, SPARKS-X, RAMPAGE, FT Site and PyMol. Out of 327 nsSNPs form human MTHFD1 gene, total 45 SNPs were predicted as functionally most significant SNPs, among which 17 were highly conserved and functional, 17 were highly conserved and structural residues. Among 45 most significant SNPs, 15 were predicted to be involved in post translational modifications. The p.Gly165Arg may interfere in homodimer interface formation. The p.Asn439Lys and p.Asp445Asn may interfere in binding interactions of MTHFD1 protein with cesium cation and potassium. The two SNPs (p.Asp562Gly and p.Gly637Cys) might interfere in interactions of MTHFD1 with ligand.

In Silico Analysis of nsSNPs of Carp TLR22 Gene Affecting its Binding Ability with Poly I:C

Immune response mediated by toll-like receptor 22 (TLR22), only found in teleost/amphibians, is triggered by double-stranded RNA binding to its LRR (leucine-rich repeats) ecto-domain. Accumulated evidences suggested that missense mutations in TLR genes affect its function. However, information on mutation linked pathogen recognition for TLR22 was lacking. The present study was commenced for predicting the effect of non-synonymous single-nucleotide polymorphisms (nsSNPs) on the pathogen recognizable LRR domain of TLR22 of farmed carp, Labeo rohita. The sequence-based algorithms (SIFT, PROVEAN and I-Mutant2.0) indicated that three SNPs (out of 27) such as p.L159F (rs76759876) and p.L529P (rs749355507) of LRR, and p.I836M (rs750758397) of intracellular motifs could potentially disrupt protein function. The 3D structure was generated using MODELLER 9.13 and further validated by SAVEs server. The simulated molecular docking of native TLR22 and mutants with poly I:C ligand indicated that mutations positioned at p.L159F and p.L529P of the LRR region affects the binding affinity significantly. This is the first kind of study of predicting nsSNPs of teleost TLR22 with disturbed ligand binding affinity with its extra-cellular LRR domain and thereby likely hindrance in subsequent signal transduction. This study serves as a guide for in vivo evaluation of impact of mutation on immune response mediated by teleost TLR22 gene.

Identification of Deleterious Mutations in Myostatin Gene of Rohu Carp (Labeo rohita) Using Modeling and Molecular Dynamic Simulation Approaches

The myostatin (MSTN) is a known negative growth regulator of skeletal muscle. The mutated myostatin showed a double-muscular phenotype having a positive significance for the farmed animals. Consequently, adequate information is not available in the teleosts, including farmed rohu carp, Labeo rohita. In the absence of experimental evidence, computational algorithms were utilized in predicting the impact of point mutation of rohu myostatin, especially its structural and functional relationships. The four mutations were generated at different positions (p.D76A, p.Q204P, p.C312Y, and p.D313A) of MSTN protein of rohu. The impacts of each mutant were analyzed using SIFT, I-Mutant 2.0, PANTHER, and PROVEAN, wherein two substitutions (p.D76A and p.Q204P) were predicted as deleterious. The comparative structural analysis of each mutant protein with the native was explored using 3D modeling as well as molecular-dynamic simulation techniques. The simulation showed altered dynamic behaviors concerning RMSD and RMSF, for either p.D76A or p.Q204P substitution, when compared with the native counterpart. Interestingly, incorporated two mutations imposed a significant negative impact on protein structure and stability. The present study provided the first-hand information in identifying possible amino acids, where mutations could be incorporated into MSTN gene of rohu carp including other carps for undertaking further in vivo studies.

In silico Analysis Revealed High-risk Single Nucleotide Polymorphisms in Human Pentraxin-3 Gene and their Impact on Innate Immune Response against Microbial Pathogens

Pentraxin-3 (PTX-3) protein is an evolutionary conserved protein that acts as a soluble pattern-recognition receptor for pathogens and plays important role in innate immune response. It recognizes various pathogens by interacting with extracellular moieties such as glactomannan of conidia (Aspergillus fumigatus), lipopolysaccharide of Pseudomonas aeruginosa, Streptococcus pneumonia and Salmonella typhimurium. Thus, PTX-3 protein helps to clear these pathogens by activating downstream innate immune process. In this study, computational methods were used to analyze various non-synonymous single nucleotide polymorphisms (nsSNPs) in PTX-3 gene. Three different databases were used to retrieve SNP data sets followed by seven different in silico algorithms to screen nsSNPs in PTX-3 gene. Sequence homology based approach was used to identify nsSNPs. Conservation profile of PTX-3 protein amino acid residues were predicted by ConSurf web server. In total, 10 high-risk nsSNPs were identified in pentraxin-domain of PTX-3 gene. Out of these 10 high-risk nsSNPs, 4 were present in the conserved structural and functional residues of the pentraxin-domain, hence, selected for structural analyses. The results showed alteration in the putative structure of pentraxin-domain. Prediction of protein–protein interactions analysis showed association of PTX-3 protein with C1q component of complement pathway. Different functional and structural residues along with various putative phosphorylation sites and evolutionary relationship were also predicted for PTX-3 protein. This is the first extensive computational analyses of pentraxin protein family with nsSNPs and will serve as a valuable resource for future population based studies.

Identification of hub glycogenes and their nsSNP analysis from mouse RNA-Seq data

Glycogenes regulate a large number of biological processes such as cancer and development. In this work, we created an interaction network of 923 glycogenes to detect potential hubs from different mouse tissues using RNA-Seq data. DAVID functional cluster analysis revealed enrichment of immune response, glycoprotein and cholesterol metabolic processes. We also explored nsSNPs that may modify the expression and function of identified hubs using computational methods. We observe that the number of nsSNPs predicted by any two methods to affect protein function is 4, 7 and 2 for FLT1, NID2 and TNFRSF1B. Residues in the native and mutant proteins were analyzed for solvent accessibility and secondary structure change. Analysis of hubs can help in determining their degree of conservation and understanding their functions in biological processes. The nsSNPs proposed in this work may be further targeted through experimental methods for understanding structural and functional relationships of hub mutants.

Molecular Dynamics: New Frontier in Personalized Medicine

The field of drug discovery has witnessed infinite development over the last decade with the demand for discovery of novel efficient lead compounds. Although the development of novel compounds in this field has seen large failure, a breakthrough in this area might be the establishment of personalized medicine. The trend of personalized medicine has shown stupendous growth being a hot topic after the successful completion of Human Genome Project and 1000 genomes pilot project. Genomic variant such as SNPs play a vital role with respect to inter individual's disease susceptibility and drug response. Hence, identification of such genetic variants has to be performed before administration of a drug. This process requires high-end techniques to understand the complexity of the molecules which might bring an insight to understand the compounds at their molecular level. To sustenance this, field of bioinformatics plays a crucial role in revealing the molecular mechanism of the mutation and thereby designing a drug for an individual in fast and affordable manner. High-end computational methods, such as molecular dynamics (MD) simulation has proved to be a constitutive approach to detecting the minor changes associated with an SNP for better understanding of the structural and functional relationship. The parameters used in molecular dynamic simulation elucidate different properties of a macromolecule, such as protein stability and flexibility. MD along with docking analysis can reveal the synergetic effect of an SNP in protein-ligand interaction and provides a foundation for designing a particular drug molecule for an individual. This compelling application of computational power and the advent of other technologies have paved a promising way toward personalized medicine. In this in-depth review, we tried to highlight the different wings of MD toward personalized medicine.

In silico analysis of consequences of non-synonymous SNPs of Slc11a2 gene in Indian bovines

The aim of our study was to analyze the consequences of non-synonymous SNPs in Slc11a2 gene using bioinformatic tools. There is a current need of efficient bioinformatic tools for in-depth analysis of data generated by the next generation sequencing technologies. SNPs are known to play an imperative role in understanding the genetic basis of many genetic diseases. Slc11a2 is one of the major metal transporter families in mammals and plays a critical role in host defenses. In this study, we performed a comprehensive analysis of the impact of all non-synonymous SNPs in this gene using multiple tools like SIFT, PROVEAN, I-Mutant and PANTHER. Among the total 124 SNPs obtained from amplicon sequencing of Slc11a2 gene by Ion Torrent PGM involving 10 individuals of Gir cattle and Murrah buffalo each, we found 22 non-synonymous. Comparing the prediction of these 4 methods, 5 nsSNPs (G369R, Y374C, A377V, Q385H and N492S) were identified as deleterious. In addition, while tested out for polar interactions with other amino acids in the protein, from above 5, Y374C, Q385H and N492S showed a change in interaction pattern and further confirmed by an increase in total energy after energy minimizations in case of mutant protein compared to the native.

•

22 nsSNPs were predicted to decrease the stability of protein based on I-Mutant.

•

From these SNPs, 5 was identified as deleterious by SIFT, PROVEAN, and PANTHER.

•

Y374C, Q385H and N492S were found to be damaging.

Mechanisms of CFTR functional variants that impair regulated bicarbonate permeation and increase risk for pancreatitis but not for cystic fibrosis

CFTR is a dynamically regulated anion channel. Intracellular WNK1-SPAK activation causes CFTR to change permeability and conductance characteristics from a chloride-preferring to bicarbonate-preferring channel through unknown mechanisms. Two severe CFTR mutations (CFTRsev) cause complete loss of CFTR function and result in cystic fibrosis (CF), a severe genetic disorder affecting sweat glands, nasal sinuses, lungs, pancreas, liver, intestines, and male reproductive system. We hypothesize that those CFTR mutations that disrupt the WNK1-SPAK activation mechanisms cause a selective, bicarbonate defect in channel function (CFTRBD) affecting organs that utilize CFTR for bicarbonate secretion (e.g. the pancreas, nasal sinus, vas deferens) but do not cause typical CF. To understand the structural and functional requirements of the CFTR bicarbonate-preferring channel, we (a) screened 984 well-phenotyped pancreatitis cases for candidate CFTRBD mutations from among 81 previously described CFTR variants; (b) conducted electrophysiology studies on clones of variants found in pancreatitis but not CF; (c) computationally constructed a new, complete structural model of CFTR for molecular dynamics simulation of wild-type and mutant variants; and (d) tested the newly defined CFTRBD variants for disease in non-pancreas organs utilizing CFTR for bicarbonate secretion. Nine variants (CFTR R74Q, R75Q, R117H, R170H, L967S, L997F, D1152H, S1235R, and D1270N) not associated with typical CF were associated with pancreatitis (OR 1.5, p = 0.002). Clones expressed in HEK 293T cells had normal chloride but not bicarbonate permeability and conductance with WNK1-SPAK activation. Molecular dynamics simulations suggest physical restriction of the CFTR channel and altered dynamic channel regulation. Comparing pancreatitis patients and controls, CFTRBD increased risk for rhinosinusitis (OR 2.3, p<0.005) and male infertility (OR 395, p<<0.0001). WNK1-SPAK pathway-activated increases in CFTR bicarbonate permeability are altered by CFTRBD variants through multiple mechanisms. CFTRBD variants are associated with clinically significant disorders of the pancreas, sinuses, and male reproductive system.

Identification of functional SNPs in BARD1 gene and in silico analysis of damaging SNPs: based on data procured from dbSNP database

BACKGROUND

The BARD1 gene encodes for the BRCA1-associated RING domain (BARD1) protein. Germ line and somatic mutations in BARD1 are found in sporadic breast, ovarian and uterine cancers. There is a plethora of single nucleotide polymorphisms (SNPs) which may or may not be involved in the onset of female cancers. Hence, before planning a larger population study, it is advisable to sort out the possible functional SNPs. To accomplish this goal, data available in the dbSNP database and different computer programs can be used. To the best of our knowledge, until now there has been no such study on record for the BARD1 gene. Therefore, this study was undertaken to find the functional nsSNPs in BARD1.

RESULT

2.85% of all SNPs in the dbSNP database were present in the coding regions. SIFT predicted 11 out of 50 nsSNPs as not tolerable and PolyPhen assessed 27 out of 50 nsSNPs as damaging. FastSNP revealed that the rs58253676 SNP in the 3' UTR may have splicing regulator and enhancer functions. In the 5' UTR, rs17489363 and rs17426219 may alter the transcriptional binding site. The intronic region SNP rs67822872 may have a medium-high risk level. The protein structures 1JM7, 3C5R and 2NTE were predicted by PDBSum and shared 100% similarity with the BARD1 amino acid sequence. Among the predicted nsSNPs, rs4986841, rs111367604, rs13389423 and rs139785364 were identified as deleterious and damaging by the SIFT and PolyPhen programs. Additionally, I-Mutant showed a decrease in stability for these nsSNPs upon mutation. Finally, the ExPASy-PROSIT program revealed that the predicted deleterious mutations are contained in the ankyrin ring and BRCT domains.

CONCLUSION

Using the available bioinformatics tools and the data present in the dbSNP database, the four nsSNPs, rs4986841, rs111367604, rs13389423 and rs139785364, were identified as deleterious, reducing the protein stability of BARD1. Hence, these SNPs can be used for the larger population-based studies of female cancers.

Screening and structural evaluation of deleterious Non-Synonymous SNPs of ePHA2 gene involved in susceptibility to cataract formation

Age-related cataract is clinically and genetically heterogeneous disorder affecting the ocular lens, and the leading cause of vision loss and blindness worldwide. Here we screened nonsynonymous single nucleotide polymorphisms (nsSNPs) of a novel gene, EPHA2 responsible for age related cataracts. The SNPs were retrieved from dbSNP. Using I-Mutant, protein stability change was calculated. The potentially functional nsSNPs and their effect on protein was predicted by PolyPhen and SIFT respectively. FASTSNP was used for functional analysis and estimation of risk score. The functional impact on the EPHA2 protein was evaluated by using SWISSPDB viewer and NOMAD-Ref server. Our analysis revealed 16 SNPs as nonsynonymous out of which 6 nsSNPs, namely rs11543934, rs2291806, rs1058371, rs1058370, rs79100278 and rs113882203 were found to be least stable by I-Mutant 2.0 with DDG value of > -1.0. nsSNPs, namely rs35903225, rs2291806, rs1058372, rs1058370, rs79100278 and rs113882203 showed a highly deleterious tolerance index score of 0.00 by SIFT server. Four nsSNPs namely rs11543934, rs2291806, rs1058370 and rs113882203 were found to be probably damaging with PSIC score of ≥ 2. 0 by Polyp hen server. Three nsSNPs namely, rs11543934, rs2291806 and rs1058370 were found to be highly polymorphic with a risk score of 3-4 with a possible effect of Non-conservative change and splicing regulation by FASTSNP. The total energy and RMSD value was higher for the mutant-type structure compared to the native type structure. We concluded that the nsSNP namely rs2291806 as the potential functional polymorphic that is likely to have functional impact on the EPHA2 gene.

Phenotype-optimized sequence ensembles substantially improve prediction of disease-causing mutation in cystic fibrosis

Cystic fibrosis transmembrane conductance regulator (CFTR) mutation is associated with a phenotypic spectrum that includes cystic fibrosis (CF). The disease liability of some common CFTR mutations is known, but rare mutations are seen in too few patients to categorize unequivocally, making genetic diagnosis difficult. Computational methods can predict the impact of mutation, but prediction specificity is often below that required for clinical utility. Here, we present a novel supervised learning approach for predicting CF from CFTR missense mutation. The algorithm begins by constructing custom multiple sequence alignments called phenotype-optimized sequence ensembles (POSEs). POSEs are constructed iteratively, by selecting sequences that optimize predictive performance on a training set of CFTR mutations of known clinical significance. Next, we predict CF disease liability from a different set of CFTR mutations (test-set mutations). This approach achieves improved prediction performance relative to popular methods recently assessed using the same test-set mutations. Of clinical significance, our method achieves 94% prediction specificity. Because databases such as HGMD and locus-specific mutation databases are growing rapidly, methods that automatically tailor their predictions for a specific phenotype may be of immediate utility. If the performance achieved here generalizes to other systems, the approach could be an excellent tool to help establish genetic diagnoses.

Screening and Evaluation of Deleterious SNPs in APOE Gene of Alzheimer's Disease

Introduction. Apolipoprotein E (APOE) is an important risk factor for Alzheimer's disease (AD) and is present in 30–50% of patients who develop late-onset AD. Several single-nucleotide polymorphisms (SNPs) are present in APOE gene which act as the biomarkers for exploring the genetic basis of this disease. The objective of this study is to identify deleterious nsSNPs associated with APOE gene. Methods. The SNPs were retrieved from dbSNP. Using I-Mutant, protein stability change was calculated. The potentially functional nonsynonymous (ns) SNPs and their effect on protein was predicted by PolyPhen and SIFT, respectively. FASTSNP was used for functional analysis and estimation of risk score. The functional impact on the APOE protein was evaluated by using Swiss PDB viewer and NOMAD-Ref server. Results. Six nsSNPs were found to be least stable by I-Mutant 2.0 with DDG value of >−1.0. Four nsSNPs showed a highly deleterious tolerance index score of 0.00. Nine nsSNPs were found to be probably damaging with position-specific independent counts (PSICs) score of ≥2.0. Seven nsSNPs were found to be highly polymorphic with a risk score of 3-4. The total energies and root-mean-square deviation (RMSD) values were higher for three mutant-type structures compared to the native modeled structure. Conclusion. We concluded that three nsSNPs, namely, rs11542041, rs11542040, and rs11542034, to be potentially functional polymorphic.

Path to facilitate the prediction of functional amino acid substitutions in red blood cell disorders--a computational approach

Background

A major area of effort in current genomics is to distinguish mutations that are functionally neutral from those that contribute to disease. Single Nucleotide Polymorphisms (SNPs) are amino acid substitutions that currently account for approximately half of the known gene lesions responsible for human inherited diseases. As a result, the prediction of non-synonymous SNPs (nsSNPs) that affect protein functions and relate to disease is an important task.

Principal Findings

In this study, we performed a comprehensive analysis of deleterious SNPs at both functional and structural level in the respective genes associated with red blood cell metabolism disorders using bioinformatics tools. We analyzed the variants in Glucose-6-phosphate dehydrogenase (G6PD) and isoforms of Pyruvate Kinase (PKLR & PKM2) genes responsible for major red blood cell disorders. Deleterious nsSNPs were categorized based on empirical rule and support vector machine based methods to predict the impact on protein functions. Furthermore, we modeled mutant proteins and compared them with the native protein for evaluation of protein structure stability.

Significance

We argue here that bioinformatics tools can play an important role in addressing the complexity of the underlying genetic basis of Red Blood Cell disorders. Based on our investigation, we report here the potential candidate SNPs, for future studies in human Red Blood Cell disorders. Current study also demonstrates the presence of other deleterious mutations and also endorses with in vivo experimental studies. Our approach will present the application of computational tools in understanding functional variation from the perspective of structure, expression, evolution and phenotype.

Cataloguing functionally relevant polymorphisms in gene DNA ligase I: a computational approach

A computational approach for identifying functionally relevant SNPs in gene LIG1 has been proposed. LIG1 is a crucial gene which is involved in excision repair pathways and mutations in this gene may lead to increase sensitivity towards DNA damaging agents. A total of 792 SNPs were reported to be associated with gene LIG1 in dbSNP. Different web server namely SIFT, PolyPhen, CUPSAT, FASTSNP, MAPPER and dbSMR were used to identify potentially functional SNPs in gene LIG1. SIFT, PolyPhen and CUPSAT servers predicted eleven nsSNPs to be intolerant, thirteen nsSNP to be damaging and two nsSNPs have the potential to destabilize protein structure. The nsSNP rs11666150 was predicted to be damaging by all three servers and its mutant structure showed significant increase in overall energy. FASTSNP predicted twenty SNPs to be present in splicing modifier binding sites while rSNP module from MAPPER server predicted nine SNPs to influence the binding of transcription factors. The results from the study may provide vital clues in establishing affect of polymorphism on phenotype and in elucidating drug response.

Mitochondrial complex 1 gene analysis in keratoconus

PURPOSE

Keratoconus is characterized by the thinning of corneal stroma, resulting in reduced vision. The exact etiology of keratoconus (KC) is still unknown. The involvement of oxidative stress (OS) in this disease has been reported. However, the exact mechanism of OS in keratoconus is still unknown. Thus we planned this study to screen mitochondrial complex I genes for sequence changes in keratoconus patients and controls, as mitochondrial complex I is the chief source of reactive oxygen species (ROS) production.

METHODS

A total of 20 keratoconus cases and 20 healthy controls without any ocular disorder were enrolled in this study. Mitochondrial complex I genes (ND1, 2, 3, 4, 4L, 5, and 6) were amplified in all patients and controls using 12 pairs of primers by PCR. After sequencing, DNA sequences were analyzed against the mitochondrial reference sequence NC_012920. Haplogroup frequency based Principle Component Analysis (PCA) was constructed to determine whether the gene pool of keratoconus patients is closer to major populations in India.

RESULTS

DNA sequencing revealed a total 84 nucleotide variations in patients and 29 in controls. Of 84 nucleotide changes, 18 variations were non-synonymous and two novel frame-shift mutations were detected in cases. Non-synonymous mtDNA sequence variations may account for increased ROS and decreased ATP production. This ultimately leads to OS; which is a known cause for variety of corneal abnormalities. Haplotype analysis showed that most of the patients were clustered under the haplogroups: T, C4a2a, R2'TJ, M21'Q1a, M12'G2a2a, M8'CZ and M7a2a, which are present as negligible frequency in normal Indian population, whereas only few patients were found to be a part of the other haplogroups like U7 (Indo-European), R2 and R31, whose origin is contentious.

CONCLUSIONS

Mt complex I sequence variations are the main cause of elevated ROS production which leads oxidative stress. This oxidative stress then starts a cascade of events which ultimately can lead to keratoconus. Prompt antioxidant therapy should be initiated in keratoconus patients to minimize ROS related damage.

In silico searching for disease-associated functional DNA variants

Experimental analyses of disease-associated DNA variants have provided significant insights into the functional implications of sequence variation. However, such experiment-based approaches for identifying functional DNA variants from a pool with a large number of neutral variants are challenging. Computational biology has the opportunity to play an important role in the identification of functional DNA variants in large-scale genotyping studies, ultimately yielding new drug targets and biomarkers. This chapter outlines in silico methods to predict disease-associated functional DNA variants so that the number of DNA variants screened for association with disease can be reduced to those that are most likely to alter gene function. To explore possible relationships between genetic mutations and phenotypic variation, different computational methods like Sorting Intolerant from Tolerant (SIFT, an evolutionary-based approach), Polymorphism Phenotyping (PolyPhen, a structure-based approach) and PupaSuite are discussed for prioritization of high-risk DNA variants. The PupaSuite tool aims to predict the phenotypic effect of DNA variants on the structure and function of the affected protein as well as the effect of variants in the non-coding regions of the same genes. To further investigate the possible causes of disease at the molecular level, deleterious nonsynonymous variants can be mapped to 3D protein structures. An analysis of solvent accessibility and secondary structure can also be performed to understand the impact of a mutation on protein function and stability. This chapter demonstrates a 'real-world' application of some existing bioinformatics tools for DNA variant analysis.

Evaluation of the disease liability of CFTR variants

Over 1600 novel sequence variants in the CFTR gene have been reported to the CF Mutation Database (http://www.genet.sickkids.on.ca/cftr/Home.html). While about 25 mutations are well characterized by clinical studies and functional assays, the disease liability of most of the remaining mutations is either unclear or unknown. This gap in knowledge has implications for diagnosis, therapy selection, and counseling for patients and families carrying an uncharacterized CFTR mutation. This chapter will describe a critical approach to assessing the disease implications of CFTR mutations utilizing clinical data, literature review, functional testing, and bioinformatic in silico methods.

Axenfeld-Rieger Syndrome Associated with Congenital Glaucoma and Cytochrome P4501B1 Gene Mutations

Developmental anomalies of the ocular anterior chamber angle may lead to an incomplete development of the structures that form the conventional aqueous outflow pathway. Thus, disorders that present with such dysfunction tend to be associated with glaucoma. Among them, Axenfeld-Rieger (ARS) malformation is a rare clinical entity with an estimated prevalence of one in every 200,000 individuals. The changes in eye morphogenesis in ARS are highly penetrant and are associated with 50% risk of development of glaucoma. Mutations in the cytochrome P4501B1 (CYP1B1) gene have been reported to be associated with primary congenital glaucoma and other forms of glaucoma and mutations in pituitary homeobox 2 (PITX2) gene have been identified in ARS in various studies. This case was negative for PITX2 mutations and compound heterozygote for CYP1B1 mutations. Clinical manifestations of this patient include bilateral elevated intraocular pressure (>40 mmHg) with increased corneal diameter (>14 mm) and corneal opacity. Patient also had iridocorneal adhesions, anteriorly displaced Schwalbe line, anterior insertion of iris, broad nasal bridge and protruding umbilicus. This is the first study from north India reporting CYP1B1 mutations in Axenfeld-Rieger syndrome with bilateral buphthalmos and early onset glaucoma. Result of this study supports the role of CYP1B1 as a causative gene in ASD disorders and its role in oculogenesis.

A comprehensive in silico analysis of the functional and structural impact of SNPs in the IGF1R gene

Insulin-like growth factor 1 receptor (IGF1R) acts as a critical mediator of cell proliferation and survival. Many single nucleotide polymorphisms (SNPs) found in the IGF1R gene have been associated with various diseases, including both breast and prostate cancer. The genetics of these diseases could be better understood by knowing the functions of these SNPs. In this study, we performed a comprehensive analysis of the functional and structural impact of all known SNPs in this gene using publicly available computational prediction tools. Out of a total of 2412 SNPs in IGF1R retrieved from dbSNP, we found 32 nsSNPs, 58 sSNPs, 83 mRNA 3' UTR SNPs, and 2225 intronic SNPs. Among the nsSNPs, a total of six missense nsSNPs were found to be damaging by both a sequence homology-based tool (SIFT) and a structural homology-based method (PolyPhen), and one nonsense nsSNP was found. Further, we modeled mutant proteins and compared the total energy values with the native IGF1R protein, and showed that a mutation from arginine to cysteine at position 1216 (rs61740868) on the surface of the protein caused the greatest impact on stability. Also, the FASTSNP tool suggested that 31 sSNPs and 3 intronic SNPs might affect splicing regulation. Based on our investigation, we report potential candidate SNPs for future studies on IGF1R mutations.

Mitochondrial DNA analysis in primary congenital glaucoma

PURPOSE

To screen mitochondrial DNA (mtDNA) for nucleotide variations in primary congenital glaucoma (PCG).

METHODS

The entire coding region of the mitochondrial genome was amplified by polymerase chain reaction from 35 PCG patients and 40 controls. The full mtDNA genome except the D-loop was sequenced. All sequences were analyzed against mitochondrial reference sequence NC_012920.

RESULTS

MtDNA sequencing revealed a total of 132 and 58 nucleotide variations in PCG and controls, respectively. Of 132 nucleotide variations, 42 (31.81%) were non-synonymous and 82 (62.12%) were synonymous changes, and 8 were in RNA genes. The highest number of nucleotide variations were recorded in complex I followed by complex IV, then complex V. Eight patients (22.85%) had potentially pathogenic mtDNA nucleotide changes and twenty (57.14%) had mtDNA sequence changes associated with elevated reactive oxygen species (ROS) production. Mitochondria not only constitute the energy-generating system in the cell, but are also critically involved in calcium signaling and apoptosis. Mitochondrial function can be affected by mutations in mitochondrial and nuclear DNA, chemical insults to components of the electron transport chain, and a lack of substrates such as oxygen. Mitochondrial dysfunction results in an excessive generation of free radicals and reduced mitochondrial respiration. Developing trabecular meshwork (TM) is deficient in antioxidant enzymes, and thus is more susceptible to oxidative stress (OS) induced damage. Previous studies have documented certain mtDNA sequence variations associated with elevated ROS levels and OS. Three such changes (G10398A, A12308G, and G13708A) were present in our patients. Elevated ROS may cause OS. This OS may further damage mtDNA and may cause decreased mitochondrial respiration. This may lead to impaired growth, development and differentiation of TM and consequently trabecular-dysgenesis, which is a characteristic feature of PCG. OS affects both TM and retinal ganglion cells (RGCs) and may be involved in the neuronal death affecting the optic nerve in glaucoma. There are several studies which point to mitochondrial dysfunction in different types of glaucoma and critically participate in RGC death. Recent studies also implicate mitochondrial dysfunction-associated OS as a risk factor for glaucoma patients. It has been reported that elevated hydrostatic pressure causes breakdown of the mitochondrial network by mitochondrial fission and induce cristae depletion and cellular ATP reduction in differentiated RGC-5 cells in vitro as well as in vivo.

CONCLUSIONS

A total of 44 novel mtDNA variations were identified in this study. Non-synonymous mtDNA variations may adversely affect respiratory chain, impair OXPHOS pathway result in low ATP production, high ROS production and impair growth, development and differentiation of TM lead to trabecular-dysgenesis and consequently RGC's death. Such cases with mtDNA variations and consequent OS may benefit by early diagnosis and prompt management by antioxidant therapy. This may delay OS induced injury to TM and RGCs and hence improve visual prognosis.

Identification of four novel cytochrome P4501B1 mutations (p.I94X, p.H279D, p.Q340H, and p.K433K) in primary congenital glaucoma patients

PURPOSE

Primary congenital glaucoma (PCG) is an autosomal recessive eye disorder that is postulated to result from developmental defects in the anterior eye segment. Mutations in the cytochrome P4501B1 (CYP1B1) gene are a predominant cause of congenital glaucoma. In this study we identify CYP1B1 mutations in PCG patients.

METHODS

Twenty-three unrelated PCG patients and 50 healthy controls were enrolled in the study. CYP1B1 was screened for mutations by PCR and DNA sequencing.

RESULTS

DNA sequencing revealed a total of 15 mutations. Out of these, four (p.I94X, p.H279D, p.Q340H, and p.K433K) were novel mutations and five were known pathogenic mutations. Five coding single nucleotide polymorphisms and one intronic single nucleotide polymorphism (rs2617266) were also found. Truncating mutations (p.I94X and p.R355X) were associated with the most severe disease phenotype. It is possible that patients with two null alleles with no catalytic activity may present with a more severe phenotype of the disease compared to patients with one null allele (heterozygous). The disease phenotype of patients with CYP1B1 mutations was more severe compared with the clinical phenotype of patients negative for CYP1B1 mutations.

CONCLUSION

Mutations in CYP1B1 are a major cause for PCG in our patients. Identifying mutations in subjects at risk of developing glaucoma, particularly among relatives of PCG patients, is of clinical significance. These developments may help in reducing the disease frequency in familial cases. Such studies will be of benefit in the identification of pathogenic mutations in different populations and will enable us to develop simple and rapid diagnostic tests for analyzing such cases.