Using structural bioinformatics to investigate the impact of non synonymous SNPs and disease mutations: scope and limitations
© Reumers et al; licensee BioMed Central Ltd. 2009
Published: 27 August 2009
Linking structural effects of mutations to functional outcomes is a major issue in structural bioinformatics, and many tools and studies have shown that specific structural properties such as stability and residue burial can be used to distinguish neutral variations and disease associated mutations.
We have investigated 39 structural properties on a set of SNPs and disease mutations from the Uniprot Knowledge Base that could be mapped on high quality crystal structures and show that none of these properties can be used as a sole classification criterion to separate the two data sets. Furthermore, we have reviewed the annotation process from mutation to result and identified the liabilities in each step.
Although excellent annotation results of various research groups underline the great potential of using structural bioinformatics to investigate the mechanisms underlying disease, the interpretation of such annotations cannot always be extrapolated to proteome wide variation studies. Difficulties for large-scale studies can be found both on the technical level, i.e. the scarcity of data and the incompleteness of the structural tool suites, and on the conceptual level, i.e. the correct interpretation of the results in a cellular context.
The molecular phenotype of a coding non synonymous SNP or disease associated mutation describes the functional and structural properties of a protein that are affected by a single amino acid substitution . In this study we want to address whether the concept of the in silico determined molecular phenotype can be employed for large-scale classification of SNPs and disease mutations. The attempt to classify a large set of mutations based on an incomplete molecular phenotype may seem naive at first glance, had it not been suggested that individual properties such as protein stability, the accessibility of the amino acid substitution site, and the location of variants in surface pockets are predictive determinants of the phenotypic effect of a variation [2–5]. A comparative study of protein stability predictors by Blundell and co-workers demonstrated that although protein stability changes caused by mutation can be relatively accurately estimated in silico, these predictions by themselves do not yield accuracy on large-scale classification between benign and disruptive mutations [6–8].
Furthermore, computational analyses rely heavily on the quality of the data under scrutiny and the computational methods used to evaluate these data. Before investigating 39 structural properties of proteins and amino acid substitutions for their predictive power regarding SNP classification, we have investigated what major liabilities are encountered when implementing an structural approach to SNP annotation and classification. The results are compared with those achieved by the best performers among the state-of-the-art tools.
Results and discussion
In this study we have identified the common issues that are encountered when performing large-scale analyses of structural properties of human coding variation. The first issue concerns the availability of structural data for nsSNPs and disease mutations, while the second involves the availability of computational tools to predict structural properties. The last issue concerns the quality of classification: are the training and evaluation data sets used in the analyses sufficient to extrapolate results for larger studies, and do the properties used have sufficient predictive power to separate the two data sets?
Structural coverage of human genetic variation
Despite structural genomics projects, the gap between sequence and structural information is still wide, and the coverage of variation data with structural data is estimated to be as low as 14% . We have investigated the boundaries of structural coverage by varying the quality requirements on the structural model (Supplementary Figure S1A in Additional file 1), the sequence identity between query sequence and modelled structure (Figure S1B), the percentage of the wild type sequence covered by the structural model (Figure S1C), and the length of the alignment between query and target (Figure S1D). Circa 12% of all nsSNPs present in the Ensembl Variation Database (release 44) can be mapped on a structural model, in accordance with the estimate cited previously. However, this percentage is valid only when no restrictions regarding sequence identity, sequence coverage or structure quality are applied.
In Figure S1A we see that the number of SNPs covered by structural data drops after 40% sequence identity. Requirements on sequence identity sufficient for prediction are different for various methods. Yue & Moult  found a sequence identity of 40% sufficient for accurate prediction, while Chasman & Adams  obtained the best results with identities higher than 60%. However, these methods do not use full atomic detail to assess the structural properties of an amino acid substitution, and thus to do not require high sequence identity to be able to model the substitution. We use the FoldX force field to model amino acid variation on structural models, which uses an all-atom representation of the structure. Although this introduces high accuracy of stability estimation [7, 9, 10], it also requires high quality structural models. Our standard restrictions on building high-confidence structural models using the FoldX force field are X-ray structures with a resolution lower than 2.5 Å and sequence identity higher than 80%. Applying these restrictions to the Ensembl data results in a data set of 5416 nsSNPs (circa 4% of the data, Figure S1B).
Summary of structural coverage of SNP data.
nsSNPs covered by high quality structural data
No additional criteria
Sequence coverage > 80 or alignment length > 100
Sequence identity > 80
Sequence coverage > 80 or alignment length > 100, and sequence identity > 80
Highly reliable nsSNPs covered by high quality structural data
Doublehit validation status, MAF > 0.01
Doublehit validation status, MAF > 0.01, sequence identity > 80
Doublehit validation status, MAF > 0.01, sequence coverage > 80 or alignment length > 100
Doublehit validation status, MAF > 0.01, sequence coverage > 80 or alignment length > 100, and sequence identity > 80
Predictability of structural properties
The second issue for a large-scale structural bioinformatics approach is the structural properties that are predictable with state of the art tools: how well can we describe the structural behaviour of a protein and its mutants? Previous structural studies have identified protein stability, aggregation and misfolding as determinants of correct functioning on the single protein level [8, 11, 12]. Mutations affecting the functional sites of a protein, such as DNA, ligand and protein interaction sites, are not considered within this scope, but the investigation of these sites will most certainly be of great importance to assess the impact of amino acid substitutions.
Tools have been developed that describe the structure and dynamics of a protein: stability, aggregation, amyloidosis, and folding. We have used computational methods that are capable of assessing the effects of a mutation on protein stability (FoldX), aggregation (Tango) and amyloidosis (Waltz). Although algorithms exist that can predict folding of small single domain proteins (e.g. Rosetta , FoldX , SimFold ), to date no computational method exists that can predict folding events on large multi-domain proteins, or that is applicable in genome wide studies.
Although we have limited ourselves to the analysis of structural features of single protein molecules, and have not investigated protein-protein interactions in this study, we have included an analysis of the binding of proteins to molecular chaperones, as it is directly related to correct folding of the protein. The high abundance of chaperones in the cell emphasises their crucial role in the protein quality control system , but this is not reflected in the availability of computational tools for chaperone binding. We have used the only available tool, the Hsp70 binding predictor Limbo , to assess chaperone binding variation caused by amino acid alteration.
The predictive power of structural properties
Following the recommendations of Care et al , we have used the SwissProt annotated disease and polymorphism data (SwissProt Variation Index release 52) as the evaluation data for our analyses. Mapping of these variants on high quality structural models (X-ray structures with resolution ≤ 2.5 Å, sequence identity with the model above 80%) yielded a data set of 240 positive (disease-associated) mutations and 400 negative variations (neutral nsSNPs) in 98 proteins. To ensure that the analyses are comparable, we applied the sequence based predictors to the same small data set as the predictors that use 3D structures or structural models.
Before we evaluated the discriminative power of the individual structural parameters, we wanted to assess whether our data showed distinguishable patterns for three important parameters. The first two criteria, stability difference and the degree of burial of the mutation site, have previously been identified as providing information about the severity of a mutation [5, 19]. The third criterion is difference in aggregation propensity, which has been cited as likely to be an important factor in disease susceptibility [12, 20] but thus far has not been applied in a proteome wide mutation analysis.
Predictive power of structural properties of the modeled variant proteins.
FoldX energy evaluation
Overall stability of residue
Backbone H bond
Sidechain H bond
Entropy side chain
Entropy main chain
Van der Waals contribution
Van der Waals clash
Side chain burial
Main chain burial
Entropy by sampling of possible side chain conformations
Entropy side chain
Predictive power of the differences between wild type and variant proteins for different structural properties.
FoldX energy evaluation
Overall stability difference
Overall stability diff. (surface)
Overall stability diff. (buried)
Backbone H bond
Sidechain H bond
Entropy main chain
Entropy side chain
Van der Waals contribution
Van der Waals clash
Entropy difference by sampling of possible side chain conformations
FoldX entropy difference
FoldX entropy diff. (buried)
FoldX entropy diff. (surface)
Tango (positive, more aggr.)
Tango (negative, less aggr.)
Waltz (positive, more aggr.)
Waltz (negative, less aggr.)
The results obtained from these detailed analyses are unanimous: none of the parameters evaluated can be used to separate the data. All MCC values are close to zero, and thus the predictions are no better than a random predictor would perform on the data. The high accuracy of FoldX for quantitative stability prediction has been proven in various studies [7, 9, 10], so we have high confidence in our stability estimations. In accordance with the analyses of Worth and co-workers , we find that high stability differences alone are no sufficient criterion to distinguish deleterious mutations and neutral variation. These results show that the dominant effect of for instance stability that was proposed in earlier large-scale studies [5, 21] can not always be generalised for other data.
The fact that none of the properties representing conformational differences between wild type and variant protein contain enough information to distinguish neutral and deleterious variation implies that large-scale classification based on singular structural properties is not feasible and requires a better understanding of how the complex interplay between biophysical and biochemical properties of a protein conspire to different tolerance for mutations in different proteins. Although we can predict the structural and functional impact of a mutation of a protein, it is not always feasible to translate this into a prediction of the overall phenotypic effect, i.e. will the mutation result in a disease phenotype or not.
Recent studies that combine structural and evolutionary information using machine learning techniques are able to classify relatively large data sets obtained for the SwissProt database successfully (summarised in Supplementary Table S2). Although the combination of these two types of information improves classification of disease mutations greatly, the incorporation of sequence conservation measures may obscure the mechanism underlying disease. Low frequency substitutions at conserved positions suggest that the mutation will not be tolerated, but will not teach us what the underlying reason of disease is. Although knowing that an amino acid is critical for correct function is of course useful, in a structural bioinformatics approach the focus is more on the molecular mechanism underlying disease.
A simple combination of the SNPeffect structural bioinformatics toolsuite on our evaluation data set showed that in our case, at least a linear combination of these methods is not sufficient to classify the data (TPR = 0.73, TNR = 0.27, MCC = 0). Although we have not fully explored the predictive power of the properties in a more sophisticated approach, such as machine learning techniques that use non-linear combinations (e.g. neural networks, support vector machines), the results obtained in the previous analyses have highlighted a major issue in disease mutation classification. The failure of the classification is mainly due to false positives, i.e. neutral mutations that are predicted to be deleterious. To assess the "predictiveness" of our data set, we applied the well-established evolutionary method SIFT  to our data and found that SIFT was also not able to classify effectively the data. In fact the results were even worse than our naive classifier (TPR = 0.69, TNR = 0.21, MCC = -0.12).
As an illustration of the influence of the data set used for evaluation on the performance of a predictor, we list the results for the variation in performance of SNP classification of SIFT, that uses evolutionary information to label SNPs (Supplementary Table S3). The Matthews correlation coefficient varies between -0.12 on our data set over 0.25 on human mutagenesis data, up to 0.59 on the HIV-1 protease mutagenesis set in the original SIFT paper . Although the methodology and underlying data (i.e. the BLOCKS database) is certainly sound, and there is no question that SIFT in most cases can be trusted to evaluate whether or not an amino acid change is tolerated in an evolutionary sense, this variability in classification success illustrates the importance of the choice of training and test data to build and evaluate predictors. SIFT was trained to classify mutations that disrupt the function of a protein, and may suffer from the same limitations as our structural approach. The ability to predict which mutations will affect function does not imply the ability to predict which mutations cause disease.
To date, we have not explored a machine learning technique that incorporates the functional effects predicted by our tools. Some of the problems above may be improved by using non-linear combinations of the structural properties. Since we pose great value in the interpretation of the classification, rule based techniques such as decision trees are our prime choice of machine learning technique. Studies using random forests to classify SNPs show results similar to other state-of-the-art classification techniques as support vector machines , and we plan to implement such techniques using the structural parameters described in this manuscript in the near future.
The concept of using the molecular phenotypic effect of a nsSNP to assess its effect on the structure and function of the protein it alters was first introduced by Bork and co-workers . The question has been raised to how much of this molecular phenotype is necessary to evaluate the contribution of a SNP to a disease phenotype: are there singular dominant properties that determine the impairment of structure and function, or do we need to consider the full ensemble of molecular properties to interpret the impact of the SNP? Other research groups have proposed that single properties such as stability  and solvent accessibility  can be used to classify SNPs. We have examined all the individual structural bioinformatics tools that were proposed in the SNPeffect toolsuite  for their ability to act as a binary classifier for deleterious and neutral SNPs. Neither of the individual properties that were examined could serve this purpose. Because several approaches were able to classify similar data sets as the one we have used, we applied the most used evolutionary method, SIFT , to our data set. As it was not able to classify our data set accurately, we argued that generalisation of the results presented by the state of the art classifiers might be an important issue. We illustrated this problem with the variability of performance of SIFT on 8 different data sets used in various analyses.
From these analyses we concluded that strict classification of SNPs is not feasible at the time, both because there are still many technical difficulties to overcome, and because the biological interpretation of the molecular phenotype in relation to a disease phenotype is a complex matter. Even at the single molecule level, we cannot assess how tolerant a specific protein is to structural variation. The inherent rigidity of a protein might influence the change in stability that is allowed before severe conformational changes are introduced. Furthermore, on the cellular level biological interpretation is even harder: we can not predict the role of the protein quality control system plays in this tolerance level, not all interactions are described at the molecular level, and much more. Even if we can predict the molecular effect accurately, this might not necessarily result in a disease phenotype because of functional redundancy of the protein.
However, not being able to classify human variation into disease mutations and neutral or beneficial variation does not mean that this approach or the methods developed are useless. By using high quality bioinformatics tools, we can select from a large pool of variations the candidates that are interesting for detailed investigation. This in itself is a valuable contribution, because the amount of variation data available is too massive to be investigated experimentally. In silico analyses can and will be used successfully as an addition to in vitro and in vivo studies.
Assembly of data sets
Statistics on the structural coverage and validation status of human non synonymous coding SNPs were performed on data from the Ensembl human variation database release 44, containing 12.2 million SNPs, of which 133698 cause an amino acid variation in a known transcript. The mapping of SNPs on protein structures was evaluated using the "ensppdbmapping" DAS service provided by the SPICE server . Positive and negative data sets (disease related mutations and polymorphisms) for the evaluation of SNP classification were designed with data from the SwissProt variation index  in the UniProt knowledge base (version 52.0, March 2007, ) that were mapped onto known PDB structures and high quality homologs thereof. The quality criteria described in the results section (models with resolution of 2.5 Å or higher, sequence identity of 80% or more) lead to structural models of 400 SNPs (negative) and 240 disease associated mutations (positive).
Structural bioinformatics tools
We have used the FoldX force field  for all mutant properties regarding structural location, protein stability and its various components, the Tango  and Waltz [, submitted] algorithms to assess the propensity for aggregation of wild type and variant proteins, and the Limbo algorithm [, submitted] to evaluate the chaperone-binding properties of amino acid sequences. A novel tool developed by Lenaerts et al (unpublished) was used to estimate the entropy of a specific amino acid site in a high-resolution structure.
The FoldX force field was developed for the fast and accurate estimation of the free change upon mutation on the stability of a protein or a protein complex [9, 14, 29, 32]. It uses an all-atom representation of these macromolecules, and has been validated on a test database of more than 1000 mutants from more than 20 different proteins. It currently yields a correlation of 0.78 with a standard deviation of 0.41 kcal/mol.
Modelling and evaluation of mutations in FoldX is performed with the BuildModel command. It is used first to model a homologous sequence on a structural model and to optimise the side chains to fit the new sequence, and then to evaluate the effect of a single amino acid variation. The Gibbs free energy of a protein is calculated with the Stability command. The various structural parameters used in the classification tests (backbone clash, backbone H bond formation, sidechain H bond formation, electrostatics, solvation of hydrophobic residues, solvation of polar residues, torsion clash, Van der Waals contribution, Van der Waals clash)
Entropy calculations based on side chain sampling
The method uses a rotamer database based on conditional statistics of dihedral angles derived from the WHAT IF data set . All amino acids from this data and their corresponding dihedral angles (10° bin) were used to derive the following probabilities: P(χ i ), P(χ i |χi-1) and P(χ i |χi-1, χi-2), except for χ1(P(χ1) and P(χ1|ϕ, ψ)). A set of n random rotamers can be derived from the probability distribution thus calculated. This will allow sampling of rotamers with greater resolution than classical rotamer libraries.
The free energy of each change is determined with FoldX.
The β-aggregation prediction algorithm Tango  uses a statistical mechanics approach to represent a competition between major conformational states: the random coil and the native conformations, as well as β-turn, α-helix and β-aggregate. Two windows of variable length slide over the sequence, and each such window can populate these conformational states according to a Boltzmann distribution. The frequency of population of each structural state for a given segment will be relative to its energy, which is derived from statistical and empirical parameters. To predict the β-aggregating segments of a peptide, Tango calculates the partition function of the phase space involving these conformational states. In our analysis we have used Tango to calculate the difference in aggregation tendency that results from an single amino acid variation.
Current methods for the prediction of the sequence determinants of amyloidosis suffer from two major problems: overpredicting amorphous cross β aggregates and missing amylogenic sequences that are enriched in the polar Q and N residues, such as the prion protein. The Waltz algorithm [, submitted] tackles these problems by taking into account amyloid hexapeptides from 48 new amyloid forming sequences, derived from 31 proteins. About half the proteins in this extended data set were not previously known to contain amyloidogenic sequences such as presenilin-2, titin and myosin. Waltz combines terms from amino acid sequence scoring in the learning set, physical property analysis and homology modelling. The method shows 84% sensitivity at 92% specificity on the AmylHex data set , and correctly identifies mutations in human proteins known to be associated with amyloid deposition.
Limbo is a Hsp70 binding site predictor that was built using a dual method combining sequence and structural information [, submitted]. Experimental DnaK binding data of 53 non-redundant peptide sequences was used to generate a sequence-based position-specific scoring matrix (PSSM) based on logarithm of the odds scores. Following an in silico alanine scan of the substrate peptide in the crystal structure of a DnaK-substrate complex (PDBID 1DKX) using FoldX, a structure-based PSSM that reflects the individual contribution of certain substrate residue types for DnaK binding was generated. The Limbo DnaK binding site predictor was obtained by combining the structure-based PSSM with a normalisation factor of 0.2 with the sequence-based PSSM. Limbo is able to correctly predict 89% of the true positives in a tested peptide set (high sensitivity), with a concurrent amount of only 5.9% false positives for a specific score threshold (high specificity). The robustness of the predictor was evaluated with a cross-validation test, resulting in a true positive rate of 72% true positives and a false positive rate of 5.9%. The predictor was able to identify an entire known DnaK binding site in the heat-shock promoter σ 32 . We have used Limbo to rank mutated proteins according to their DnaK binding affinity.
Evaluation of classification
Other measures, such as the proportion of correct predictions, are not useful when the two classes are of very different sizes. For example, assigning every object to the larger set achieves a high proportion of correct predictions, but is not generally a useful classification.
The MCC is generally regarded as a balanced measure which can be used even if the positive and negative classes are of very different sizes. It returns a value between -1 and +1. A coefficient of +1 represents a perfect prediction, 0 an average random prediction and -1 the worst possible prediction.
In a Receiver Operating Characteristic (ROC) curve the true positive rate (TPR) is plotted in function of the false positive rate (FPR) for different cut-off points. Each point on the ROC plot represents a sensitivity/specificity pair corresponding to a particular decision threshold. A test with perfect discrimination between positives and negatives has a ROC plot that passes through the upper left corner (100% sensitivity, 100% specificity). Therefore the closer the ROC plot is to the upper left corner, the higher the overall accuracy of the test .
Joke Reumers was supported by the Institute for the encouragement of Scientific Research and Innovation of Brussels (ISRIB), Belgium.
This article has been published as part of BMC Bioinformatics Volume 10 Supplement 8, 2009: Proceedings of the European Conference on Computational Biology (ECCB) 2008 Workshop: Annotation, interpretation and management of mutations. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/10?issue=S8.
- Sunyaev S, Lathe Wr, Bork P: Integration of genome data and protein structures: prediction of protein folds, protein interactions and "molecular phenotypes" of single nucleotide polymorphisms. Curr Opin Struct Biol 2001, 11: 125–130. 10.1016/S0959-440X(00)00175-5View ArticlePubMedGoogle Scholar
- Chasman D, Adams RM: Predicting the functional consequences of non-synonymous single nucleotide polymorphisms: Structure-based assessment of amino acid variation. J Mol Biol 2001, 307(2):683–706. 10.1006/jmbi.2001.4510View ArticlePubMedGoogle Scholar
- Ferrer-Costa C, Orozco M, de la Cruz X: Characterization of disease-associated single amino acid polymorphisms in terms of sequence and structure properties. J Mol Biol 2002, 315(4):771–786. 10.1006/jmbi.2001.5255View ArticlePubMedGoogle Scholar
- Stitziel NO, Tseng YY, Pervouchine D, Goddeau D, Kasif S, Liang J: Structural location of disease-associated single-nucleotide polymorphisms. J Mol Biol 2003, 327(5):1021–1030. 10.1016/S0022-2836(03)00240-7View ArticlePubMedGoogle Scholar
- Yue P, Li Z, Moult J: Loss of protein structure stability as a major causative factor in monogenic disease. J Mol Biol 2005, 353(2):459–473. 10.1016/j.jmb.2005.08.020View ArticlePubMedGoogle Scholar
- Worth CL, Burke DF, Blundell TL: Estimating the effects of single nucleotide polymorphisms on protein structure: how good are we at identifying likely disease associated mutations? Proceedings of Molecular Interactions – Bringing Chemistry to Life 2006.Google Scholar
- Burke DF, Worth CL, Priego EM, Cheng T, Smink LJ, Todd JA, Blundell TL: Genome bioinformatic analysis of nonsynonymous SNPs. BMC Bioinformatics 2007, 8: 301. 10.1186/1471-2105-8-301PubMed CentralView ArticlePubMedGoogle Scholar
- Worth CL, Bickerton GRJ, Schreyer A, Forman JR, Cheng TMK, Lee S, Gong S, Burke DF, Blundell TL: A structural bioinformatics approach to the analysis of nonsynonymous single nucleotide polymorphisms (nsSNPs) and their relation to disease. J Bioinform Comput Biol 2007, 5(6):1297–1318. 10.1142/S0219720007003120View ArticlePubMedGoogle Scholar
- Guerois R, Nielsen JE, Serrano L: Predicting changes in the stability of proteins and protein complexes: A study of more than 1000 mutations. J Mol Biol 2002, 320(2):369–387. 10.1016/S0022-2836(02)00442-4View ArticlePubMedGoogle Scholar
- Tokuriki N, Stricher F, Schymkowitz J, Serrano L, Tawfik DS: The stability effects of protein mutations appear to be universally distributed. J Mol Biol 2007, 369(5):1318–1332. 10.1016/j.jmb.2007.03.069View ArticlePubMedGoogle Scholar
- Steward RE, MacArthur MW, Laskowski RA, Thornton JM: Molecular basis of inherited diseases: a structural perspective. Trends Genet 2003, 19(9):505–513. 10.1016/S0168-9525(03)00195-1View ArticlePubMedGoogle Scholar
- DePristo M, Weinreich D, Hartl D: Missense meanderings in sequence space: A biophysical view of protein evolution. Nature Reviews Genetics 2005. AOP. AOP.Google Scholar
- Simons KT, Bonneau R, Ruczinski I, Baker D: Ab initio protein structure prediction of CASP III targets using ROSETTA. Proteins 1999., (Suppl 3):Google Scholar
- Serrano L, Guerois R: Fold-X: An algorithm to predict and engineer folding pathways. Abstr Pap Am Chem Soc 2001, 221: U395-U395.Google Scholar
- Fujitsuka Y, Chikenji G, Takada S: SimFold energy function for de novo protein structure prediction: consensus with Rosetta. Proteins 2006, 62(2):381–398. 10.1002/prot.20748View ArticlePubMedGoogle Scholar
- Soti C, Csermely P: Protein stress and stress proteins: implications in aging and disease. J Biosci 2007., 32:Google Scholar
- Van Durme J, Maurer-Stroh S, Wilkinson H, Rousseau F, Schymkowitz J: Accurate prediction of the sequence determinants of DnaK-peptide binding via a method that integrates homology modelling and experimental data. 2009, in press.Google Scholar
- Carvalho MA, Marsillac SM, Karchin R, Manoukian S, Grist S, Swaby R, Urmenyi T, Rondinelli E, Silva R, Gayol L, Baumbach L, Sutphen R, Pickard-Brzosowicz J, Nathanson K, Sali A, Goldgar D, Couch F, Radice P, Monteiro AN: Determination of cancer risk associated with germ line BRCA1 missense variants by functional analysis. Cancer Res 2007, 67(4):1494–1501. 10.1158/0008-5472.CAN-06-3297PubMed CentralView ArticlePubMedGoogle Scholar
- Ramensky V, Bork P, Sunyaev S: Human non-synonymous SNPs: server and survey. Nucleic Acid Res 2002, 30(17):3894–3900. 10.1093/nar/gkf493PubMed CentralView ArticlePubMedGoogle Scholar
- Worth CL, Blundell TL: Estimating the effects of SNPs on protein structure: loss of protein interactions and stability as indicators of mis-function and disease-association. Curr Top Biochem Res 2008, in press.Google Scholar
- Yue P, Melamud E, Moult J: SNPs3D: candidate gene and SNP selection for association studies. BMC Bioinformatics 2006, 7: 166. 10.1186/1471-2105-7-166PubMed CentralView ArticlePubMedGoogle Scholar
- Ng PC, Henikoff S: Predicting deleterious amino acid substitutions. Genome Res 2001, 11(5):863–874. 10.1101/gr.176601PubMed CentralView ArticlePubMedGoogle Scholar
- Bao L, Cui Y: Prediction of the phenotypic effects of non-synonymous single nucleotide polymorphisms using structural and evolutionary information. Bioinformatics 2005, 21(10):2185–2190. 10.1093/bioinformatics/bti365View ArticlePubMedGoogle Scholar
- Reumers J, Conde L, Medina I, Maurer-Stroh S, Van Durme J, Dopazo J, Rousseau F, Schymkowitz J: Joint annotation of coding and non-coding single nucleotide polymorphisms and mutations in the SNPeffect and PupaSuite databases. Nucleic Acids Res 2008, (36 Database):D825–9.Google Scholar
- Ng PC, Henikoff S: SIFT: predicting amino acid changes that affect protein function. Nucleic Acid Res 2003, 31(13):3812–3814. 10.1093/nar/gkg509PubMed CentralView ArticlePubMedGoogle Scholar
- Prlic A, Down TA, Hubbard TJ: Adding some SPICE to DAS. Bioinformatics 2005, 21(Suppl 2):ii40–1. 10.1093/bioinformatics/bti1106PubMed CentralView ArticlePubMedGoogle Scholar
- Yip YL, Famiglietti M, Gos A, Duek PD, David FPA, Gateau A, Bairoch A: Annotating single amino acid polymorphisms in the UniProt/Swiss-Prot knowledgebase. Hum Mutat 2008, 29(3):361–366. 10.1002/humu.20671View ArticlePubMedGoogle Scholar
- UniProt Consortium: The Universal Protein Resource (UniProt). Nucleic Acids Res 2007, (35 Database):D193–7. 10.1093/nar/gkl929Google Scholar
- Schymkowitz JWH, Rousseau F, Martins IC, Ferkinghoff-Borg J, Stricher F, Serrano L: Prediction of water and metal binding sites and their affinities by using the Fold-X force field. Proc Natl Acad Sci USA 2005, 102(29):10147–10152. 10.1073/pnas.0501980102PubMed CentralView ArticlePubMedGoogle Scholar
- Fernandez-Escamilla AM, Rousseau F, Schymkowitz J, Serrano L: Prediction of sequence-dependent and mutational effects on the aggregation of peptides and proteins. Nat Biotechnol 2004, 22(10):1302–1306. 10.1038/nbt1012View ArticlePubMedGoogle Scholar
- Maurer-Stroh S, Debulpaep M, Kuemmerer N, Lopez de la Paz M, Martins I, Reumers J, Copland A, Serpell L, Serrano L, Rousseau F, Schymkowitz J: An exploration of the sequence determinants of amyloid formation and the development of the WALTZ prediction algorithm. 2009, in press.Google Scholar
- Schymkowitz J, Borg J, Stricher F, Nys R, Rousseau F, Serrano L: The FoldX web server: an online force field. Nucleic Acid Res 2005, 33: W382-W388. 10.1093/nar/gki387PubMed CentralView ArticlePubMedGoogle Scholar
- Vriend G: What If – a molecular modeling and drug design program. J Mol Graph 1990, 8: 52. 10.1016/0263-7855(90)80070-VView ArticlePubMedGoogle Scholar
- Thompson MJ, Sievers SA, Karanicolas J, Ivanova MI, Baker D, Eisenberg D: The 3D profile method for identifying fibril-forming segments of proteins. Proc Natl Acad Sci USA 2006, 103: 4074–4078. 10.1073/pnas.0511295103PubMed CentralView ArticlePubMedGoogle Scholar
- McCarty JS, Rudiger S, Schonfeld HJ, Schneider-Mergener J, Nakahigashi K, Yura T, Bukau B: Regulatory region C of the E. coli heat shock transcription factor, sigma32, constitutes a DnaK binding site and is conserved among eubacteria. J Mol Biol 1996, 256(5):829–37. 10.1006/jmbi.1996.0129View ArticlePubMedGoogle Scholar
- Matthews BW: Comparison of the predicted and observed secondary structure of T4 phage lysozyme. Biochim Biophys Acta 1975, 405(2):442–451.View ArticlePubMedGoogle Scholar
- Zweig MH, Campbell G: Receiver-operating characteristic (ROC) plots: a fundamental evaluation tool in clinical medicine. Clin Chem 1993, 39(4):561–577.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.