- Open Access
Disease-related mutations predicted to impact protein function
BMC Genomics volume 13, Article number: S11 (2012)
Non-synonymous single nucleotide polymorphisms (nsSNPs) alter the protein sequence and can cause disease. The impact has been described by reliable experiments for relatively few mutations. Here, we study predictions for functional impact of disease-annotated mutations from OMIM, PMD and Swiss-Prot and of variants not linked to disease.
Most disease-causing mutations were predicted to impact protein function. More surprisingly, the raw predictions scores for disease-causing mutations were higher than the scores for the function-altering data set originally used for developing the prediction method (here SNAP). We might expect that diseases are caused by change-of-function mutations. However, it is surprising how well prediction methods developed for different purposes identify this link. Conversely, our predictions suggest that the set of nsSNPs not currently linked to diseases contains very few strong disease associations to be discovered.
Firstly, annotations of disease-causing nsSNPs are on average so reliable that they can be used as proxies for functional impact. Secondly, disease-causing nsSNPs can be identified very well by methods that predict the impact of mutations on protein function. This implies that the existing prediction methods provide a very good means of choosing a set of suspect SNPs relevant for disease.
Evolution leads to genetic diversity
The selection of survival under changing conditions guides the cell’s genetic makeup (“genotype”) that is dynamically fit for retaining important cellular functions (“phenotype”). Today’s genetic landscape represents the current state of a sampling process that continuously creates new phenotypes. This process yields genetic variation across and within species. In human, single nucleotide polymorphisms (SNPs) are essential for genetic diversity [1, 2]. Non-synonymous SNPs (nsSNPs) alter the amino acid sequence. Some of these mutations affect protein structure and/or function and could increase susceptibility to disease.
Do disease-causing mutations impact protein function?
Disease-causing mutations occur often inside the protein (buried) and at hydrogen-bonding residues [3–5]. Protein function is often associated with evolutionarily conserved residues [4, 6–9]. Most known disease-related nsSNPs in proteins of known 3D (three-dimensional) structure appear to affect structurally important residues and sites relevant for function . For instance, disease-associated mutations can affect protein interactions . In protein kinases, they have been shown to cluster into the functionally important catalytic core [11, 12]. The above trends confirm the expectation that mutations cause disease because they damage important proteins.
Experts have established the above trends by laboriously inspecting small sets of well-curated proteins. Could less well-versed experts with better algorithms have established valid trends about disease-causing mutations for large data set by automatically extracting data set of disease-related mutations and their predicted functional effects? At OMIM’s infancy, a few years ago, we failed to accomplish this; i.e. observed trends did not differ much from random. This has changed. Here, we provide data that strongly suggest an affirmative answer to the question and demonstrate that we have a large repository of disease-causing mutations. To pick the most important practical result of our work: today’s disease-causing mutations can serve as an excellent proxy for “change of function”.
We used SNPdbe  as the underlying source for amino acid substitutions, functional effect annotations and disease relations. This comprehensive new resource integrates variants from dbSNP , Swiss-Prot , PMD , and OMIM  and annotations of functional effects (from Swiss-Prot and PMD) and disease (from SwissVar , PMD and OMIM). The term ‘genetic disease’ is rather heterogeneous, covering Mendelian, monogenic disorders and polygenic diseases, exhibiting more complex genotypic patterns. Here, we do not differentiate between the different disease-types. Instead we aim at analyzing all disease-causing mutations.
We created the following five subsets from SNPdbe (Additional file 2). (1) Set of disease-related + observed effect mutations: We collected 1,105 human nsSNPs (from 217 proteins) that were annotated to be both disease-causing and functionally non-neutral. (2) Set of disease-related mutations: We obtained a set of amino acid substitutions in human proteins with disease-association. We extracted 26,404 mutations (3,419 proteins) with disease annotations but no annotated functional effect. (3) Set of observed effect mutations: We collected 36,317 mutants in 3,790 proteins with experimentally observed effect. We excluded mutations with disease associations. This set constitutes a part of the “functional effects” sets annotated in PMD; it served as the positive training set for SNAP . Note that after our filtering the resulting set of mutations with observed effect and the set of disease-related mutants did NOT overlap. (4) Set of mutations with unknown disease relation : We extracted 251,414 variants (28,913 proteins) without known disease associations. (5) Set of random mutations: We randomly selected one mutation in each of the 28,913 proteins from the set of mutants of unknown disease relation such that the mutated position was maximally distant from any other mutation observed in the given protein.
Prediction of effect
For the vast majority of point mutants (single amino acid changes or nsSNPs) in human, the impact on protein function remains unknown. For all mutations in the above four data sets (disease-causing, disease-relation unknown, observed function-changing, and random), we predicted their effects on function with SNAP  and SIFT . Both methods provide binary classifications (effect/neutral) along with a more detailed score. SNAP scores range from -100 (strongly predicted as neutral) to 100 (strongly predicted to change function); the distance from the binary decision boundary (0) measures the reliability of the effect. Essentially, stronger predictions are also more reliable, i.e. the higher the score, the more likely the mutation impacts function [19, 21, 22]. For a small data set, we previously established that SNAP scores correlate with the severity of change; i.e. high (positive) SNAP scores relate to more severe functional effects [19, 21, 22].
SIFT  scores range from 0 to 1 and aim at characterizing the normalized probability of tolerable amino acid substitution. Values < 0.05 imply prediction of functional change; all other values are considered neutral. As with many other prediction methods, the distance to the decision boundary (0.05) reflects the reliability of a particular prediction . For many prediction methods developed in our group (protein-protein binding [24–26], protein-DNA binding , backbone flexibility ), the strength of an effect correlated with prediction strength, e.g. ISIS predicted binding hot spots stronger than other residues involved in the interaction . Although we never used the strength of an effect to train our methods, this correlation is intuitive: stronger effects are more consistent and therefore become stronger carved into the machine-learning model. Similarly, SIFT scores could be used to prioritize amino acid substitutions . In this perspective, we consider the distance from the default decision boundary (0.05) as the magnitude of the effect.
SNAP and SIFT aspire to solve the same problem with different means. SNAP was trained on literature-derived  mutants that are either functionally similar to the wild-type (neutral) or alter function (effect) in either direction (decrease and increase of function). SIFT on the other hand infers probabilities of functional change from residue conservation in alignments of evolutionarily related proteins. While SNAP operates on an experimentally substantiated definition of change, SIFT uses conservation scores of amino acids as a proxy for functional change. Although both methods largely capture the underlying biological meaning of functional change, their predictions disagree often. Thus, the methods are likely orthogonal, picking up different aspects of protein function.
In addition, we applied PhD-SNP  to predict whether mutations in all five sets are disease-causing or neutral. PhD-SNP offers several modes striking different balances between runtime and performance. We used the most accurate mode that uses both sequence and evolutionary profiles.
We represented our resulting distributions using box plots [30, 31]. The lower and upper box edges depict the first and third quartiles of the distributions, respectively. The length of the box is the interquartile range of the distribution. The bold bar inside the box represents the median, while dashed lines reach to the most extreme data points, that are no more than 1.5 times the interquartile range away from the upper or lower box edge. Note that each box covers half the distribution.
Results and discussion
Disease-causing mutations strongly predicted to change protein function
We applied SNAP and SIFT to the 26,404 annotated disease related mutants (Methods). At the default threshold, SNAP predicted over 86% of the disease related mutations to impact function (Fig. 1A, B, 2) and SIFT ~59% (Fig. 2, Additional file 1). SNAP predictions were very strong: about half of the effect predictions had levels of severity of >40 (Fig. 1B, dashed black curve).
In our experience, SNAP scores >40 are exceptional when applying the method to new data. To clarify this point, the observed effect mutations were the very same data set that trained SNAP. We ascertained that this set had no overlap with the disease related mutations (Methods). Usually, machine-learning methods perform much better on the training than on the testing set. This also holds for SNAP; hence, the distribution of SNAP scores for the training set of observed effect mutants is expected to be closer to ‘more effect’ than for any other data set. We observed the opposite (Fig. 1B: red vs. dashed black lines): effect predictions were stronger for the disease related mutations than for our observed effect training set, e.g. while just over 40% of the training set reached a score >40, 47% of the disease related mutations did. A difference of seven percentage points might not be perceived as high, but the effect is significantly higher for comparison to testing on the training set. SIFT overall also predicted the disease related mutations stronger than the observed effect data, but the difference was not significant (Additional file 1).
Do disease-related mutations with an observed effect alter function even more? We analyzed the predicted functional effect of disease-associated mutations with observed effect (disease-related+observed effect). About 90% were predicted to impact function (4% more than for disease related), while over 53% had SNAP scores higher than 40 (6% more than for disease related; Fig. 1A, B solid black line, Fig. 2). SIFT showed a similar trend: 66% in the set of disease related+observed effect compared to 59% in disease related mutations (Fig. 2, Additional file 1). This suggests that the most reliable source of impact mutations is by connecting disease relations and independent experimental observations.
As negative control, the predictions differed greatly for the 251,414 mutants with unknown disease relation. First, only about 51% of those were predicted to have an effect by SNAP (Fig. 1A, B, 2), and only 39% by SIFT (Fig. 2, Additional file 1). Second, only 12% of those had a SNAP score larger than 40 (Fig. 1B, dashed green curve).
Many mutations with unknown effect predicted to alter function
SNAP and SIFT predicted much more effect for disease related mutations than in mutants with unknown disease relation. Still, many of those mutations were predicted to change protein function. However, much fewer mutants with unknown disease relation were predicted to significantly change function than the disease related mutations (Fig. 1B: strong effect for 14% of mutants unknown disease relation - dashed green line - vs. 48% of disease related mutations - dashed black line). Comparing the prediction trends between the two data sets suggests that the mutations of unknown disease relation will never become a ‘disease-rich’ set (i.e. through newly discovered disease associations). Random mutations were even less often predicted to have strong effect (~7%, Fig. 1B, dashed blue line). This result suggests that many experimental annotations of ‘functional impact’ remain to be determined/observed for the set of mutations with unknown disease relation (roughly > 7%-14%).
Same trend found in predicted disease mutations
If disease related can serve as a good proxy for (strong) functional impact, then a method trained to predict disease-causing mutations should reveal the reverse and thus confirm the same: predicted disease is expected to be enriched in observed effect compared to mutations of unknown disease relation. We analyzed the fraction of predicted disease by applying PhD-SNP (Methods) to our five data sets. PhD-SNP predicted >64% of the observed effect mutations as disease related (Fig. 2), while only 26% of mutations with unknown disease relation were predicted to be disease associated. Furthermore, we confirmed the other observations already found in functional impact predictions: Random mutations appear to have the lowest impact on disease (only 22%, Fig. 2).
PhD-SNP predicted both disease-related sets to contain most disease mutants (86% in disease related+observed effect and 74% in disease related, Fig. 2). This was expected due to the important overlap between our data and the training set of PhD-SNP . Nonetheless, the increase in predicted-disease mutations of 12% once again suggested that observed effect mutants play a major role in disease.
Our findings show that if a mutation leads to disease then a change in function plays a major role in explaining the cause (59%-86%). This finding cannot be inverted due to the overlap of score distributions of disease related mutants and mutants with unknown disease annotation (Fig. 1A, Additional file 1); i.e. strong effect on function does not imply disease.
Our comparison between mutations annotated as disease related and those experimentally annotated function changing (observed effect) does not imply that there is anything special about disease-causing mutations. Instead, our findings highlight differences in the severity of functional effect. That is, on average, assuming that a disease causing mutation has a functional effect is more reliable than experimentally evaluating functional change.
We compared disease-associated single point mutations (nsSNPs) predicted to change protein function with those of unknown disease-association. Implicitly, we tested the reliability of annotations that link mutations to disease and the extent to which predictions of functional effect overlap with disease causation.
As opposed to other studies addressing this question [3–6, 10–12], we used predictions of functional effect to determine the fraction of deleterious point mutations in two different populations of human variants: disease related (or disease-causing) mutations and mutations without any knowledge of phenotypic effect. The major findings were: (1) annotations of disease-causation provide a good approximation of functional effect. (2) Methods developed to predict the impact of mutations onto protein function clearly identify disease-causing mutations as those that change function. In other words, their predictions provide a valuable first step towards the study of the molecular impact of disease.
This work was supported by a grant from the Alexander von Humboldt foundation through the German Ministry for Research and Education (BMBF: Bundesministerium fuer Bildung und Forschung); YB was supported by the SEBS, Rutgers, New Brunswick startup funds.
Consortium GP: A map of human genome variation from population-scale sequencing. Nature. 2010, 467 (7319): 1061-1073. 10.1038/nature09534.
Collins FS, Brooks LD, Chakravarti A: A DNA polymorphism discovery resource for research on human genetic variation. Genome Res. 1998, 8 (12): 1229-1231.
Gong S, Blundell TL: Structural and functional restraints on the occurrence of single amino acid variations in human proteins. PLoS ONE. 2010, 5 (2): e9186-10.1371/journal.pone.0009186.
Sunyaev S, Ramensky V, Bork P: Towards a structural basis of human non-synonymous single nucleotide polymorphisms. Trends Genet. 2000, 16 (5): 198-200. 10.1016/S0168-9525(00)01988-0.
Wang Z, Moult J: SNPs, protein structure, and disease. Hum Mutat. 2001, 17 (4): 263-270. 10.1002/humu.22.
Talavera D, Taylor MS, Thornton JM: The (non)malignancy of cancerous amino acidic substitutions. Proteins. 2010, 78 (3): 518-529.
Lichtarge O, Bourne HR, Cohen FE: An evolutionary trace method defines binding surfaces common to protein families. J Mol Biol. 1996, 257 (2): 342-358. 10.1006/jmbi.1996.0167.
Rausell A, Juan D, Pazos F, Valencia A: Protein interactions and ligand binding: from protein subfamilies to functional specificity. Proc Natl Acad Sci U S A. 2010, 107 (5): 1995-2000. 10.1073/pnas.0908044107.
Landgraf R, Xenarios I, Eisenberg D: Three-dimensional cluster analysis identifies interfaces and functional residue clusters in proteins. J Mol Biol. 2001, 307 (5): 1487-1502. 10.1006/jmbi.2001.4540.
Schuster-Bockler B, Bateman A: Protein interactions in human genetic diseases. Genome Biol. 2008, 9 (1): R9-10.1186/gb-2008-9-1-r9.
Torkamani A, Schork NJ: Distribution analysis of nonsynonymous polymorphisms within the human kinase gene family. Genomics. 2007, 90 (1): 49-58. 10.1016/j.ygeno.2007.03.006.
Torkamani A, Verkhivker G, Schork NJ: Cancer driver mutations in protein kinase genes. Cancer Lett. 2009, 281 (2): 117-127. 10.1016/j.canlet.2008.11.008.
Schaefer C, Meier A, Rost B, Bromberg Y: SNPdbe: constructing an nsSNP functional impacts database. Bioinformatics. 2012, 28 (4): 601-602. 10.1093/bioinformatics/btr705.
Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K: dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001, 29 (1): 308-311. 10.1093/nar/29.1.308.
Boeckmann B, Bairoch A, Apweiler R, Blatter MC, Estreicher A, Gasteiger E, Martin MJ, Michoud K, O'Donovan C, Phan I, et al: The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003, 31 (1): 365-370. 10.1093/nar/gkg095.
Kawabata T, Ota M, Nishikawa K: The Protein Mutant Database. Nucleic Acids Res. 1999, 27 (1): 355-357. 10.1093/nar/27.1.355.
Amberger J, Bocchini CA, Scott AF, Hamosh A: McKusick's Online Mendelian Inheritance in Man (OMIM). Nucleic Acids Res. 2009, 37 (Database issue): D793-796.
Mottaz A, David FP, Veuthey AL, Yip YL: Easy retrieval of single amino-acid polymorphisms and phenotype information using SwissVar. Bioinformatics. 2010, 26 (6): 851-852. 10.1093/bioinformatics/btq028.
Bromberg Y, Rost B: SNAP: predict effect of non-synonymous polymorphisms on function. Nucleic Acids Res. 2007, 35 (11): 3823-3835. 10.1093/nar/gkm238.
Ng PC, Henikoff S: Predicting deleterious amino acid substitutions. Genome Res. 2001, 11 (5): 863-874. 10.1101/gr.176601.
Bromberg Y, Overton J, Vaisse C, Leibel RL, Rost B: In silico mutagenesis: a case study of the melanocortin 4 receptor. FASEB J. 2009, 23 (9): 3059-3069. 10.1096/fj.08-127530.
Bromberg Y, Rost B: Comprehensive in silico mutagenesis highlights functionally important residues in proteins. Bioinformatics. 2008, 24 (16): i207-212. 10.1093/bioinformatics/btn268.
Ng PC, Henikoff S: SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003, 31 (13): 3812-3814. 10.1093/nar/gkg509.
Ofran Y, Rost B: Analysing six types of protein-protein interfaces. J Mol Biol. 2003, 325: 377-387. 10.1016/S0022-2836(02)01223-8.
Ofran Y, Rost B: ISIS: interaction sites identified from sequence. Bioinformatics. 2007, 23 (2): e13-16. 10.1093/bioinformatics/btl303.
Ofran Y, Rost B: Protein-protein interaction hot spots carved into sequences. PLoS Computational Biology. 2007, 3 (7): e119-10.1371/journal.pcbi.0030119.
Ofran Y, Mysore V, Rost B: Prediction of DNA-binding residues from sequence. Bioinformatics. 2007, 23 (13): i347-353. 10.1093/bioinformatics/btm174.
Schlessinger A, Yachdav G, Rost B: PROFbval: predict flexible and rigid residues in proteins. Bioinformatics. 2006, 22 (7): 891-893. 10.1093/bioinformatics/btl032.
Capriotti E, Calabrese R, Casadio R: Predicting the insurgence of human genetic diseases associated to single point protein mutations with support vector machines and evolutionary information. Bioinformatics. 2006, 22 (22): 2729-2734. 10.1093/bioinformatics/btl423.
McGill R, Tukey JW, Larsen WA: Variations of Box Plots. The American Statistician. 1978, 32 (1): 12-16.
Tukey JW: Exploratory data analysis. 1977, Reading, Mass.: Addison-Wesley Pub. Co.
Special thanks to Laszlo Kajan (TUM), Guy Yachdav (TUM/Columbia University), and Tim Karl (TUM) for help with software and hardware; to Marlena Drabik (TUM) for administrative support. Thanks to Rolf Apweiler (UniProt, EBI, Hinxton), Amos Bairoch (CALIPHO, SIB, Geneva), Ioannis Xenarios (Swiss-Prot, SIB, Geneva), their crews, and those from OMIM, PMD, SwissVar and dbSNP for maintaining excellent databases. Last, but not least, thanks to all experimentalists who enabled this analysis by making their data publicly available.
This article has been published as part of BMC Genomics Volume 13 Supplement 4, 2012: SNP-SIG 2011: Identification and annotation of SNPs in the context of structure, function and disease. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/13/S4.
CS participated in the design of the study, performed the data analysis and helped to draft the manuscript. YB participated in the design of the study and helped draft the manuscript. DA participated in the design of the study. BR participated in the coordination and design of the study and helped to draft the manuscript.
The authors declare they have no competing interests.
Electronic supplementary material
About this article
Cite this article
Schaefer, C., Bromberg, Y., Achten, D. et al. Disease-related mutations predicted to impact protein function. BMC Genomics 13, S11 (2012). https://doi.org/10.1186/1471-2164-13-S4-S11
- Functional Effect
- Decision Boundary
- Disease Association
- Functional Impact
- Impact Function