A novel computational and structural analysis of nsSNPs in CFTR gene.
TL;DR: The genetic variations that can alter the expression and function of the CFTR gene responsible for causing cystic fibrosis are analyzed using computational methods to identify potential candidates for future studies on CFTR mutations.
Abstract: Single Nucleotide Polymorphisms (SNPs) are being intensively studied to understand the biological basis of complex traits and diseases. The Genetics of human phenotype variation could be understood by knowing the functions of SNPs. In this study using computational methods, we analyzed the genetic variations that can alter the expression and function of the CFTR gene responsible candidate for causing cystic fibrosis. We applied an evolutionary perspective to screen the SNPs using a sequence homology-based SIFT tool, which suggested that 17 nsSNPs (44%) were found to be deleterious. The structure-based approach PolyPhen server suggested that 26 nsSNPS (66%) may disrupt protein function and structure. The PupaSuite tool predicted the phenotypic effect of SNPs on the structure and function of the affected protein. Structure analysis was carried out with the major mutation that occurred in the native protein coded by CFTR gene, and which is at amino acid position F508C for nsSNP with id (rs1800093). The amino acid residues in the native and mutant modeled protein were further analyzed for solvent accessibility, secondary structure and stabilizing residues to check the stability of the proteins. The SNPs were further subjected to iHAP analysis to identify htSNPs, and we report potential candidates for future studies on CFTR mutations.
Citations
More filters
TL;DR: WNK1-SPAK pathway-activated increases in CFTR bicarbonate permeability are altered by CFTRBD variants through multiple mechanisms, and molecular dynamics simulations suggest physical restriction of the CFTR channel and altered dynamic channel regulation.
Abstract: 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.
152 citations
TL;DR: 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.
Abstract: 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.
111 citations
TL;DR: A comprehensive analysis of the functional and structural impact of all known SNPs in this gene using publicly available computational prediction tools shows that a mutation from arginine to cysteine at position 1216 on the surface of the protein caused the greatest impact on stability.
Abstract: 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.
50 citations
TL;DR: This study provides a significant insight into the underlying molecular mechanism involved in albinism associated with OCA1A and determines that certain mutations can affect the dynamic properties of protein and can lead to disease conditions.
49 citations
Cites methods from "A novel computational and structura..."
...Our group also tried to evaluate the accuracy of SIFT and PolyPhen based predictions on different oncogens (Doss et al., 2008a; 2008b; Kumar and Purohit, 2012a; 2012b; Rajasekaran et al., 2008a; 2008b)....
[...]
Journal Article•
TL;DR: A total of 44 novel mtDNA variations were identified and 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.
Abstract: Purpose
To screen mitochondrial DNA (mtDNA) for nucleotide variations in primary congenital glaucoma (PCG).
48 citations
Cites background from "A novel computational and structura..."
...SIFT is a sequence homology-based tool that sorts intolerant from tolerant amino acid substitutions and predicts whether an amino acid substitution in a protein will have a phenotypic effect [19-22]....
[...]
References
More filters
TL;DR: The results of an international collaboration to produce and make freely available a draft sequence of the human genome are reported and an initial analysis is presented, describing some of the insights that can be gleaned from the sequence.
Abstract: The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.
22,269 citations
"A novel computational and structura..." refers background or methods in this paper
...The recent sequencing of the human genome (Venter et al. 2001; Lander et al. 2001) C. George Priya Doss R. Rajasekaran C. Sudandiradoss K. Ramanathan R. Purohit R. Sethumadhavan (&) Bioinformatics Division, School of Biotechnology, Chemical and Biomedical Engineering, Vellore Institute of…...
[...]
...The recent sequencing of the human genome (Venter et al. 2001; Lander et al. 2001) C....
[...]
TL;DR: A set of simple and physically motivated criteria for secondary structure, programmed as a pattern‐recognition process of hydrogen‐bonded and geometrical features extracted from x‐ray coordinates is developed.
Abstract: For a successful analysis of the relation between amino acid sequence and protein structure, an unambiguous and physically meaningful definition of secondary structure is essential. We have developed a set of simple and physically motivated criteria for secondary structure, programmed as a pattern-recognition process of hydrogen-bonded and geometrical features extracted from x-ray coordinates. Cooperative secondary structure is recognized as repeats of the elementary hydrogen-bonding patterns “turn” and “bridge.” Repeating turns are “helices,” repeating bridges are “ladders,” connected ladders are “sheets.” Geometric structure is defined in terms of the concepts torsion and curvature of differential geometry. Local chain “chirality” is the torsional handedness of four consecutive Cα positions and is positive for right-handed helices and negative for ideal twisted β-sheets. Curved pieces are defined as “bends.” Solvent “exposure” is given as the number of water molecules in possible contact with a residue. The end result is a compilation of the primary structure, including SS bonds, secondary structure, and solvent exposure of 62 different globular proteins. The presentation is in linear form: strip graphs for an overall view and strip tables for the details of each of 10.925 residues. The dictionary is also available in computer-readable form for protein structure prediction work.
14,077 citations
"A novel computational and structura..." refers methods in this paper
...We obtained the information about secondary structures of the proteins using the program DSSP (Kabsch and Sander 1983)....
[...]
...Therefore, in this work, the computational algorithms namely SIFT, PolyPhen, PupaSuite, FASTSNP, ASA View, DSSP and SRide were used to identify the deleterious nsSNPs that are likely to affect the function and structure of the protein and showed the htSNPs which are in the haplotype blocks using iHAP analysis....
[...]
TL;DR: Comparative genomic analysis indicates vertebrate expansions of genes associated with neuronal function, with tissue-specific developmental regulation, and with the hemostasis and immune systems are indicated.
Abstract: A 2.91-billion base pair (bp) consensus sequence of the euchromatic portion of the human genome was generated by the whole-genome shotgun sequencing method. The 14.8-billion bp DNA sequence was generated over 9 months from 27,271,853 high-quality sequence reads (5.11-fold coverage of the genome) from both ends of plasmid clones made from the DNA of five individuals. Two assembly strategies-a whole-genome assembly and a regional chromosome assembly-were used, each combining sequence data from Celera and the publicly funded genome effort. The public data were shredded into 550-bp segments to create a 2.9-fold coverage of those genome regions that had been sequenced, without including biases inherent in the cloning and assembly procedure used by the publicly funded group. This brought the effective coverage in the assemblies to eightfold, reducing the number and size of gaps in the final assembly over what would be obtained with 5.11-fold coverage. The two assembly strategies yielded very similar results that largely agree with independent mapping data. The assemblies effectively cover the euchromatic regions of the human chromosomes. More than 90% of the genome is in scaffold assemblies of 100,000 bp or more, and 25% of the genome is in scaffolds of 10 million bp or larger. Analysis of the genome sequence revealed 26,588 protein-encoding transcripts for which there was strong corroborating evidence and an additional approximately 12,000 computationally derived genes with mouse matches or other weak supporting evidence. Although gene-dense clusters are obvious, almost half the genes are dispersed in low G+C sequence separated by large tracts of apparently noncoding sequence. Only 1.1% of the genome is spanned by exons, whereas 24% is in introns, with 75% of the genome being intergenic DNA. Duplications of segmental blocks, ranging in size up to chromosomal lengths, are abundant throughout the genome and reveal a complex evolutionary history. Comparative genomic analysis indicates vertebrate expansions of genes associated with neuronal function, with tissue-specific developmental regulation, and with the hemostasis and immune systems. DNA sequence comparisons between the consensus sequence and publicly funded genome data provided locations of 2.1 million single-nucleotide polymorphisms (SNPs). A random pair of human haploid genomes differed at a rate of 1 bp per 1250 on average, but there was marked heterogeneity in the level of polymorphism across the genome. Less than 1% of all SNPs resulted in variation in proteins, but the task of determining which SNPs have functional consequences remains an open challenge.
12,098 citations
"A novel computational and structura..." refers background or methods in this paper
...The recent sequencing of the human genome (Venter et al. 2001; Lander et al. 2001) C. George Priya Doss R. Rajasekaran C. Sudandiradoss K. Ramanathan R. Purohit R. Sethumadhavan (&) Bioinformatics Division, School of Biotechnology, Chemical and Biomedical Engineering, Vellore Institute of…...
[...]
...The recent sequencing of the human genome (Venter et al. 2001; Lander et al. 2001) C....
[...]
TL;DR: The dbSNP database is a general catalog of genome variation to address the large-scale sampling designs required by association studies, gene mapping and evolutionary biology, and is integrated with other sources of information at NCBI such as GenBank, PubMed, LocusLink and the Human Genome Project data.
Abstract: In response to a need for a general catalog of genome variation to address the large-scale sampling designs required by association studies, gene mapping and evolutionary biology, the National Center for Biotechnology Information (NCBI) has established the dbSNP database [S.T.Sherry, M.Ward and K.Sirotkin (1999) Genome Res., 9, 677–679]. Submissions to dbSNP will be integrated with other sources of information at NCBI such as GenBank, PubMed, LocusLink and the Human Genome Project data. The complete contents of dbSNP are available to the public at website: http://www.ncbi.nlm.nih.gov/SNP. The complete contents of dbSNP can also be downloaded in multiple formats via anonymous FTP at ftp:// ncbi.nlm.nih.gov/snp/.
6,449 citations
"A novel computational and structura..." refers background in this paper
...human population (Sherry et al. 2001; Hinds et al. 2005; The International Hapmap Consortium 2003) opens the way for the development of a detailed understanding of the mechanisms by which genetic variation results in phenotype variation....
[...]
...…Vellore, Tamil Nadu 632014, India e-mail: rsethumadhavan@vit.ac.in together with the large number of SNPs present in the human population (Sherry et al. 2001; Hinds et al. 2005; The International Hapmap Consortium 2003) opens the way for the development of a detailed understanding of the…...
[...]
Baylor College of Medicine1, Chinese Academy of Sciences2, Chinese National Human Genome Center3, University of Hong Kong4, The Chinese University of Hong Kong5, Hong Kong University of Science and Technology6, Illumina7, McGill University8, Washington University in St. Louis9, University of California, San Francisco10, Wellcome Trust Sanger Institute11, Beijing Normal University12, Health Sciences University of Hokkaido13, Shinshu University14, University of Tsukuba15, Howard University16, University of Ibadan17, Case Western Reserve University18, University of Utah19, Cold Spring Harbor Laboratory20, Johns Hopkins University21, University of Oxford22, North Carolina State University23, National Institutes of Health24, Massachusetts Institute of Technology25, Chinese Academy of Social Sciences26, Kyoto University27, Nagasaki University28, Wellcome Trust29, Genome Canada30, Foundation for the National Institutes of Health31, University of Maryland, Baltimore32, Vanderbilt University33, Stanford University34, University of California, Berkeley35, New York University36, University of Oklahoma37, University of New Mexico38, Université de Montréal39, University of California, Los Angeles40, University of Michigan41, University of Wisconsin-Madison42, London School of Economics and Political Science43, Genetic Alliance44, GlaxoSmithKline45, University of Washington46, Harvard University47, University of Chicago48, Fred Hutchinson Cancer Research Center49, University of Tokyo50
TL;DR: The HapMap will allow the discovery of sequence variants that affect common disease, will facilitate development of diagnostic tools, and will enhance the ability to choose targets for therapeutic intervention.
Abstract: The goal of the International HapMap Project is to determine the common patterns of DNA sequence variation in the human genome and to make this information freely available in the public domain. An international consortium is developing a map of these patterns across the genome by determining the genotypes of one million or more sequence variants, their frequencies and the degree of association between them, in DNA samples from populations with ancestry from parts of Africa, Asia and Europe. The HapMap will allow the discovery of sequence variants that affect common disease, will facilitate development of diagnostic tools, and will enhance our ability to choose targets for therapeutic intervention.
5,926 citations