Volume 13 Supplement 4
On the effect of protein conformation diversity in discriminating among neutral and disease related single amino acid substitutions
© Juritz et al; licensee BioMed Central Ltd. 2012
Published: 18 June 2012
Non-synonymous coding SNPs (nsSNPs) that are associated to disease can also be related with alterations in protein stability. Computational methods are available to predict the effect of single amino acid substitutions (SASs) on protein stability based on a single folded structure. However, the native state of a protein is not unique and it is better represented by the ensemble of its conformers in dynamic equilibrium. The maintenance of the ensemble is essential for protein function. In this work we investigated how protein conformational diversity can affect the discrimination of neutral and disease related SASs based on protein stability estimations. For this purpose, we used 119 proteins with 803 associated SASs, 60% of which are disease related. Each protein was associated with its corresponding set of available conformers as found in the Protein Conformational Database (PCDB). Our dataset contains proteins with different extensions of conformational diversity summing up a total number of 1023 conformers.
The existence of different conformers for a given protein introduces great variability in the estimation of the protein stability (ΔΔG) after a single amino acid substitution (SAS) as computed with FoldX. Indeed, in 35% of our protein set at least one SAS can be described as stabilizing, destabilizing or neutral when a cutoff value of ±2 kcal/mol is adopted for discriminating neutral from perturbing SASs. However, when the ΔΔG variability among conformers is taken into account, the correlation among the perturbation of protein stability and the corresponding disease or neutral phenotype increases as compared with the same analysis on single protein structures. At the conformer level, we also found that the different conformers correlate in a different way to the corresponding phenotype.
Our results suggest that the consideration of conformational diversity can improve the discrimination of neutral and disease related protein SASs based on the evaluation of the corresponding Gibbs free energy change.
Human single nucleotide polymorphisms (SNPs) are the most frequent type of genetic variation in humans. Less than 1% variations are associated with non-synonymous coding SNPs (nsSNPs). About 64,971 nsSNPs are presently listed as human polymorphisms and disease single amino acid substitutions, SASs, (http://www.uniprot.org/docs/humsavar) and approximately 40% of these SASs are disease related.
It has been documented that in proteins a single amino acid substitution (SAS), can produce the loss of function in different ways. Although the less frequently found , the most obvious mechanism at the disease origin is due to change of key residues participating directly in protein function. This is the case when residue substitution occurs at the active site or in binding-sites for substrate and/or allosteric regulators [2–4]. When the biological functional unit is a complex, SASs at the subunit interface may also hamper the activity [4, 5]. A second possible mechanism is related with the perturbation of protein stability. Residue substitution can indeed destabilize the native protein fold [1, 6]. Also stabilizing residue changes have been reported to be associated with diseases [7, 8]. Furthermore, related with protein stability alteration, the origin of pathogenesis was also related with anomalous post-translational modifications  and aggregation . The correlation among protein SASs and their involvement in human diseases has been proven to be moderate , suggesting that change in protein stability is not the only source of diseases.
Protein stability can be estimated measuring the variation of Gibbs free energy (ΔΔG) between the folded and unfolded state of the protein. Most of the experimental data reported in literature are contained in ProTherm , a thermodynamic database of proteins and their variation in different organisms. Alternatively, several computational methods have been developed to estimate stability changes caused by substitution of lateral side chains in proteins (ΔΔG=ΔGwild-ΔGmutated). Most of them rely on the analysis of the energetic and/or structural perturbation introduced by the variations in the protein native structure. Although computationally intensive, early methods used all atom models to estimate ΔΔG . Soon later, simplified potentials coupled with limited conformational searches [14, 15] and the use of different types of potentials, like those based on hydrophobic interactions , secondary structure , inter-residue contacts  and knowledge-based , allowed to study the effect of different mutations in proteins in a reasonable computational time. Recently machine learning based approaches have been implemented for the prediction of ΔΔG in proteins upon residue substitution taking as input either the protein structure or sequence (for a recent review see ). The discrimination among disease related and neutral SASs can be investigated by determining ΔΔG upon residue change in the protein. This analysis is based on the notion that most harmful SASs are related to protein stability perturbation above a certain threshold ΔΔG (±1 kcal/mol). Most methods suited to predict free energy changes have been recently benchmarked in relation to their ability to discriminate disease from neutral SASs based on the corresponding ΔΔG value and their performance has been proven to be rather poor .
In most methods, predictions of SASs effects are commonly estimated using a single structure of the corresponding protein [21–24]. This approach apparently underestimates the well established concept that the native state of a protein is better represented by an ensemble of conformers [25–27]. The conformational ensemble is a key concept to explain essential properties of proteins like function [28–30], enzyme and antibody promiscuity [31, 32], enzyme catalytic power , signal transduction , protein-protein recognition  and the origin of new functions . Conformers describing the native state of a protein exist in a dynamic equilibrium which changes in response to the presence of ligands such as substrate or allosteric modulators that shift the relative conformational population [37, 38]. From a practical perspective, conformational diversity could be described using experimentally available structures of the same protein obtained in alternative conditions. As these different structures for the same protein have been obtained under different conditions (for example presence of substrate, inhibitors or allosteric activators) they can be taken as snapshots of protein dynamics and then characterize putative conformers belonging to the native ensemble [39, 40]. In this way, the description of the native ensemble of the protein will be more or less complete depending on available experimental data. A way to describe the extension of conformational diversity could be the estimation of the maximum RMSD measured between the available conformers. The distribution of the structural diversity extension measured in this way in the protein space was recently studied . The analysis involves an all vs. all comparisons between structures of the same sequence deposited in PDB database . Conformer distribution exhibits a peak at 0.3 Å RMSD with a large skew that ends at about 24 Å RMSD.
In the present work we investigated how the presence and extension of conformational diversity affect the estimation of ΔΔG in proteins with neutral and disease related SASs. Like protein function relies on the existence and preservation of the dynamic ensemble of conformers, the study of the effect of a SAS in each conformer of a given protein could help to understand the loss of function. We used a set of 803 SASs (482 disease related and 323 neutral) in 119 proteins showing different extension of conformational diversity. These proteins were taken from the Protein Conformational Database (PCDB)  a redundant collection of protein structures linked with biological information. ΔΔG for each SAS in each conformer for a given protein was estimated using FoldX . We found that the ΔΔG estimated value for a SAS highly depends on the conformer used in the estimation. In 35% and 58% of the studied proteins we found that at least one SAS could be classified as neutral, stabilizing or destabilizing depending on a ΔΔG threshold value of ±2Kcal/mol and ±1 Kcal/mol respectively. We also found that the consideration of conformational diversity increases the performance of the prediction of disease related SASs based on ΔΔG analysis. Our results show that the different conformers correlate in different ways with the phenotype (disease or neutral) and that, in most cases, one conformer per protein correlates perfectly with the corresponding phenotype. Our results indicate that the use of conformational diversity may be important to understand the effects of neutral and disease related SASs on protein stability.
Extension of protein conformational diversity
Variation of ΔΔG estimation using conformational diversity
When each SAS is mapped into its corresponding protein, we found that in 35% of the cases there is at least one SAS that can be classified as neutral and stabilizing or neutral and destabilizing depending on the considered conformer. It is interesting to note that for this 35%, the average of the maximum RMSD between conformers is 3.78 Å compared to 2.38 Å for proteins without ambiguous predicted stability changes. The discriminative threshold of ΔΔG was set at ± 2 kcal/mol, as it was previously discussed by Worth and coworkers (Worth, Preissner, and Blundell 2011). However, other works have defined different thresholds to classify the stability changes of a SAS [47–49]. When we analyzed the data distribution for ΔΔG values using a threshold of ±1 kcal/mol, the number of proteins with at least one ambiguous prediction was even higher (58%) (these calculations can be done using the data and information included as additional file 1). Thus, the extension of conformational diversity measured by the maximum RMSD between conformers increases the variability of the stability prediction of SASs generating certain ambiguity in the prediction. It is noteworthy that this ambiguity is manifested in the same proportion respect to the total population of SASs in disease-associated SASs as in polymorphic SASs. That is, the proportion of neutral SASs with respect to the total number of SASs is similar to the proportion of neutral SASs with ambiguous prediction with respect to the total number of SASs with ambiguous predictions (almost 29%). From this we can conclude that the uncertainty in the evaluation of the thermodynamic effect of a SAS equally affects neutral and disease related SASs.
Summing up, all the observations indicate that there is a large spread of estimated ΔΔG values for the different conformers due to their structural differences. The spread and, eventually, the corresponding ambiguity in the prediction of ΔΔG, can blur the correlation among ΔΔG values and the SAS type (disease related or neutral) when only one protein structure is used, as it is routinely the case.
Prediction of disease related SASs using ΔΔG
Scoring the capability of discriminating among disease related and neutral SASs on different set of conformers.
The elucidation of the effect that a single amino acid substitution (SAS) has in a specific phenotype is a central problem in different areas of research. Particularly, in protein computational biology it is very challenging to understand the mechanisms that lead to disease . Estimation of protein stability and its perturbation after single SAS have played a key role to predict the effects of SASs [1, 6, 11, 51]. Notwithstanding important advances in the area, the estimation of stability in proteins is still problematic and its correlation with disease or neutral variants highly depends on the dataset considered . Here, we improved the correlation between the perturbation of protein stability and the classification of SASs as disease related or neutral. This improvement was obtained by considering maximum values of ΔΔG computed taking into account conformational diversity. As the native state is not unique, the conservation of the biological activity of a protein relies on the different conformers of the ensemble. It has been found that the different conformers for a given protein constraints in different way the substitution pattern of proteins . Then it is expected that the different conformers should have different robustness to SASs. Although the correlation is still far from being perfect, our results suggest that the effect of each SAS should be studied in all the protein conformers in order to obtain a better understanding of protein function perturbation and the disease origin. Our results also suggest that conformational diversity could add value to new computational tools for predicting SAS effects and play a key role to obtain a deeper understanding of the relationship among protein structure and function.
Data set collection
A list of proteins with disease and polymorphic SASs was extracted from http://www.uniprot.org/docs/humsavar. This list was linked with PCDB (http://pcdb.unq.edu.ar) database . PCDB contains a collection of redundant protein structures obtained in different conditions (for example presence of ligands, change in oligomerization state, etc.). These structures can be considered as snapshots of protein dynamism, assumption validated by previous works that have proved the correspondence between structural deformations detected under different crystallographic conditions and conformational changes related to the flexibility of the native state [39, 40]. The description of the conformational ensemble of the protein native state is limited to the information deposited in PCDB and PDB databases. The maximum RMSD (RMSDmax) between alpha carbon coordinates of the different conformers, calculated with MAMMOTH , is taken as a measure of the conformational diversity of the protein. From this cross linking, a dataset containing 119 proteins with different number of conformers was defined (8.6 conformers in average per protein, with a minimum of 2 and a maximum of 73). The length of all the structures in this dataset covers more than 70% of the length of the corresponding protein sequence. The dataset contains 803 SASs with 482 disease related and 323 neutral.
The presence of cognate ligands in each structure was determined using the Procognate database . The cognate ligands deposited in this database are those involved in the biological function of the proteins. After this filtering we found that 35 proteins have bound/unbound states in our dataset.
Stability and structural measurements
For each SAS and for each protein and conformation, we estimated ΔΔG values using the program FoldX . FoldX uses an empirical potential calibrated to fit in vitro ΔΔG values. Area exposed to solvent was calculated with Naccess program (http://www.bioinf.manchester.ac.uk/naccess/).
A value of MCC=1 defines the best possible prediction, while a value of MCC=-1 indicates the worst possible prediction. A value of MCC=0 corresponds to predictions made by chance.
RC thanks the following grants: PRIN 2009 project 009WXT45Y (Italian Ministry for University and Research: MIUR), COST BMBS Action TD1101 (European Union RTD Framework Programme), and PON project PON01_02249 (Italian Ministry for University and Research: MIUR).
MSF thanks Erasmus Mundus Action 2 Lot 13a – EADIC II proyect, funded by the European Commission and coordinated by the University of Bologna for the research fellowship. GP and MSF thank the following grants: PIP CONICET (112-200801-02849) and UNQ (1004/11). EJ has a Type II fellowship from CONICET.
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.
- Wang Z, Moult J: SNPs, protein structure, and disease. Hum Mutat. 2001, 17 (4): 263-270. 10.1002/humu.22.View ArticlePubMedGoogle Scholar
- Lofgren M, Banerjee R: Loss of allostery and coenzyme B12 delivery by a pathogenic mutation in adenosyltransferase. Biochemistry. 2011, 50 (25): 5790-5798. 10.1021/bi2006306.PubMed CentralView ArticlePubMedGoogle Scholar
- Shah GN, Bonapace G, Hu PY, Strisciuglio P, Sly WS: Carbonic anhydrase II deficiency syndrome (osteopetrosis with renal tubular acidosis and brain calcification): novel mutations in CA2 identified by direct sequencing expand the opportunity for genotype-phenotype correlation. Hum Mutat. 2004, 24 (3): 272-View ArticlePubMedGoogle Scholar
- Orosz F, Olah J, Ovadi J: Triosephosphate isomerase deficiency: new insights into an enigmatic disease. Biochim Biophys Acta. 2009, 1792 (12): 1168-1174. 10.1016/j.bbadis.2009.09.012.View ArticlePubMedGoogle Scholar
- Almeida-Souza L, Goethals S, de Winter V, Dierick I, Gallardo R, Van Durme J, Irobi J, Gettemans J, Rousseau F, Schymkowitz J, et al: Increased monomerization of mutant HSPB1 leads to protein hyperactivity in Charcot-Marie-Tooth neuropathy. J Biol Chem. 2010, 285 (17): 12778-12786. 10.1074/jbc.M109.082644.PubMed CentralView 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.020.View ArticlePubMedGoogle Scholar
- Ling SC, Albuquerque CP, Han JS, Lagier-Tourenne C, Tokunaga S, Zhou H, Cleveland DW: ALS-associated mutations in TDP-43 increase its stability and promote TDP-43 complexes with FUS/TLS. Proc Natl Acad Sci U S A. 2010, 107 (30): 13318-13323. 10.1073/pnas.1008227107.PubMed CentralView ArticlePubMedGoogle Scholar
- Seidle HF, Bieganowski P, Brenner C: Disease-associated mutations inactivate AMP-lysine hydrolase activity of Aprataxin. J Biol Chem. 2005, 280 (22): 20927-20931. 10.1074/jbc.M502889200.PubMed CentralView ArticlePubMedGoogle Scholar
- Alonso Adel C, Mederlyova A, Novak M, Grundke-Iqbal I, Iqbal K: Promotion of hyperphosphorylation by frontotemporal dementia tau mutations. J Biol Chem. 2004, 279 (33): 34873-34881. 10.1074/jbc.M405131200.View ArticleGoogle Scholar
- Dobson CM: The structural basis of protein folding and its links with human disease. Philos Trans R Soc Lond B Biol Sci. 2001, 356 (1406): 133-145. 10.1098/rstb.2000.0758.PubMed CentralView ArticlePubMedGoogle Scholar
- Casadio R, Vassura M, Tiwari S, Fariselli P, Luigi Martelli P: Correlating disease-related mutations to their effect on protein stability: a large-scale analysis of the human proteome. Hum Mutat. 2011, 32 (10): 1161-1170. 10.1002/humu.21555.View ArticlePubMedGoogle Scholar
- Gromiha MM, An J, Kono H, Oobatake M, Uedaira H, Sarai A: ProTherm: Thermodynamic Database for Proteins and Mutants. Nucleic Acids Res. 1999, 27 (1): 286-288. 10.1093/nar/27.1.286.PubMed CentralView ArticlePubMedGoogle Scholar
- Bash PA, Singh UC, Langridge R, Kollman PA: Free energy calculations by computer simulation. Science. 1987, 236 (4801): 564-568. 10.1126/science.3576184.View ArticlePubMedGoogle Scholar
- Lee C: Predicting protein mutant energetics by self-consistent ensemble optimization. J Mol Biol. 1994, 236 (3): 918-939. 10.1006/jmbi.1994.1198.View ArticlePubMedGoogle Scholar
- Lee C, Levitt M: Accurate prediction of the stability and activity effects of site-directed mutagenesis on a protein core. Nature. 1991, 352 (6334): 448-451. 10.1038/352448a0.View ArticlePubMedGoogle Scholar
- Koehl P, Delarue M: Polar and nonpolar atomic environments in the protein core: implications for folding and binding. Proteins. 1994, 20 (3): 264-278. 10.1002/prot.340200307.View ArticlePubMedGoogle Scholar
- Munoz V, Serrano L: Intrinsic secondary structure propensities of the amino acids, using statistical phi-psi matrices: comparison with experimental scales. Proteins. 1994, 20 (4): 301-311. 10.1002/prot.340200403.View ArticlePubMedGoogle Scholar
- Miyazawa S, Jernigan RL: Protein stability for single substitution mutants and the extent of local compactness in the denatured state. Protein Eng. 1994, 7 (10): 1209-1220. 10.1093/protein/7.10.1209.View ArticlePubMedGoogle Scholar
- Sippl MJ: Knowledge-based potentials for proteins. Curr Opin Struct Biol. 1995, 5 (2): 229-235. 10.1016/0959-440X(95)80081-6.View ArticlePubMedGoogle Scholar
- Khan S, Vihinen M: Performance of protein stability predictors. Hum Mutat. 2010, 31 (6): 675-684. 10.1002/humu.21242.View ArticlePubMedGoogle 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-166.PubMed CentralView ArticlePubMedGoogle Scholar
- Schymkowitz J, Borg J, Stricher F, Nys R, Rousseau F, Serrano L: The FoldX web server: an online force field. Nucleic Acids Res. 2005, 33 (Web Server issue): W382-388.PubMed CentralView ArticlePubMedGoogle Scholar
- Dehouck Y, Kwasigroch JM, Gilis D, Rooman M: PoPMuSiC 2.1 : a web server for the estimation of protein stability changes upon mutation and sequence optimality. BMC Bioinformatics. 12 (1): 151-
- Reumers J, Schymkowitz J, Rousseau F: Using structural bioinformatics to investigate the impact of non synonymous SNPs and disease mutations: scope and limitations. BMC Bioinformatics. 2009, 10 (Suppl 8): S9-10.1186/1471-2105-10-S8-S9.PubMed CentralView ArticlePubMedGoogle Scholar
- James LC, Tawfik DS: Conformational diversity and protein evolution--a 60-year-old hypothesis revisited. Trends Biochem Sci. 2003, 28 (7): 361-368. 10.1016/S0968-0004(03)00135-X.View ArticlePubMedGoogle Scholar
- Lange OF, Lakomek NA, Fares C, Schroder GF, Walter KF, Becker S, Meiler J, Grubmuller H, Griesinger C, de Groot BL: Recognition dynamics up to microseconds revealed from an RDC-derived ubiquitin ensemble in solution. Science. 2008, 320 (5882): 1471-1475. 10.1126/science.1157092.View ArticlePubMedGoogle Scholar
- Tsai CD, Ma B, Kumar S, Wolfson H, Nussinov R: Protein folding: binding of conformationally fluctuating building blocks via population selection. Crit Rev Biochem Mol Biol. 2001, 36 (5): 399-433. 10.1080/20014091074228.View ArticlePubMedGoogle Scholar
- del Sol A, Tsai CJ, Ma B, Nussinov R: The origin of allosteric functional modulation: multiple pre-existing pathways. Structure. 2009, 17 (8): 1042-1050. 10.1016/j.str.2009.06.008.PubMed CentralView ArticlePubMedGoogle Scholar
- Hilser VJ: Biochemistry. An ensemble view of allostery. Science. 2010, 327 (5966): 653-654. 10.1126/science.1186121.View ArticlePubMedGoogle Scholar
- Karplus M, Kuriyan J: Molecular dynamics and protein function. Proc Natl Acad Sci U S A. 2005, 102 (19): 6679-6685. 10.1073/pnas.0408930102.PubMed CentralView ArticlePubMedGoogle Scholar
- Khersonsky O, Roodveldt C, Tawfik DS: Enzyme promiscuity: evolutionary and mechanistic aspects. Curr Opin Chem Biol. 2006, 10 (5): 498-508. 10.1016/j.cbpa.2006.08.011.View ArticlePubMedGoogle Scholar
- James LC, Roversi P, Tawfik DS: Antibody multispecificity mediated by conformational diversity. Science. 2003, 299 (5611): 1362-1367. 10.1126/science.1079731.View ArticlePubMedGoogle Scholar
- Boehr DD, McElheny D, Dyson HJ, Wright PE: The dynamic energy landscape of dihydrofolate reductase catalysis. Science. 2006, 313 (5793): 1638-1642. 10.1126/science.1130258.View ArticlePubMedGoogle Scholar
- Smock RG, Gierasch LM: Sending signals dynamically. Science. 2009, 324 (5924): 198-203. 10.1126/science.1169377.PubMed CentralView ArticlePubMedGoogle Scholar
- Yogurtcu ON, Erdemli SB, Nussinov R, Turkay M, Keskin O: Restricted mobility of conserved residues in protein-protein interfaces in molecular simulations. Biophys J. 2008, 94 (9): 3475-3485. 10.1529/biophysj.107.114835.PubMed CentralView ArticlePubMedGoogle Scholar
- Tokuriki N, Tawfik DS: Protein dynamism and evolvability. Science. 2009, 324 (5924): 203-207. 10.1126/science.1169375.View ArticlePubMedGoogle Scholar
- Ma B, Shatsky M, Wolfson HJ, Nussinov R: Multiple diverse ligands binding at a single protein site: a matter of pre-existing populations. Protein Sci. 2002, 11 (2): 184-197.PubMed CentralView ArticlePubMedGoogle Scholar
- Gunasekaran K, Ma B, Nussinov R: Is allostery an intrinsic property of all dynamic proteins?. Proteins. 2004, 57 (3): 433-443. 10.1002/prot.20232.View ArticlePubMedGoogle Scholar
- Best RB, Lindorff-Larsen K, DePristo MA, Vendruscolo M: Relation between native ensembles and experimental structures of proteins. Proc Natl Acad Sci U S A. 2006, 103 (29): 10901-10906. 10.1073/pnas.0511156103.PubMed CentralView ArticlePubMedGoogle Scholar
- Zoete V, Michielin O, Karplus M: Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility. J Mol Biol. 2002, 315 (1): 21-52. 10.1006/jmbi.2001.5173.View ArticlePubMedGoogle Scholar
- Burra PV, Zhang Y, Godzik A, Stec B: Global distribution of conformational states derived from redundant models in the PDB points to non-uniqueness of the protein structure. Proc Natl Acad Sci U S A. 2009, 106 (26): 10505-10510. 10.1073/pnas.0812152106.PubMed CentralView ArticlePubMedGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28 (1): 235-242. 10.1093/nar/28.1.235.PubMed CentralView ArticlePubMedGoogle Scholar
- Juritz J, Fernandez-Alberti S, Parisi G: PCDB: A database of proteins with conformational diversity. Nucleic Acids Res. 2010, 39: D475-9.PubMed CentralView 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-4.View ArticlePubMedGoogle Scholar
- Kumar S, Ma B, Tsai CJ, Sinha N, Nussinov R: Folding and binding cascades: dynamic landscapes and population shifts. Protein Sci. 2000, 9 (1): 10-19.PubMed CentralView 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.069.View ArticlePubMedGoogle Scholar
- Worth CL, Preissner R, Blundell TL: SDM--a server for predicting effects of mutations on protein stability and malfunction. Nucleic Acids Res. 2011Google Scholar
- Calloni G, Zoffoli S, Stefani M, Dobson CM, Chiti F: Investigating the effects of mutations on protein aggregation in the cell. J Biol Chem. 2005, 280 (11): 10607-10613. 10.1074/jbc.M412951200.View ArticlePubMedGoogle Scholar
- Tokuriki N, Tawfik DS: Stability effects of mutations and protein evolvability. Curr Opin Struct Biol. 2009, 19 (5): 596-604. 10.1016/j.sbi.2009.08.003.View ArticlePubMedGoogle Scholar
- Thusberg J, Vihinen M: Pathogenic or not? And if so, then how? Studying the effects of missense mutations using bioinformatics methods. Hum Mutat. 2009, 30 (5): 703-714. 10.1002/humu.20938.View ArticlePubMedGoogle Scholar
- Worth CL, Bickerton GR, Schreyer A, Forman JR, Cheng TM, 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/S0219720007003120.View ArticlePubMedGoogle Scholar
- Juritz J, Palopoli N, Fornasari M, Fernandez-Alberti S, Parisi G: Protein conformational diversity modulates protein divergence. Accepted Mol Biol Evol 2012. 2012Google Scholar
- Ortiz AR, Strauss CE, Olmea O: MAMMOTH (matching molecular models obtained from theory): an automated method for model comparison. Protein Sci. 2002, 11 (11): 2606-2621.PubMed CentralView ArticlePubMedGoogle Scholar
- Bashton M, Nobeli I, Thornton JM: PROCOGNATE: a cognate ligand domain mapping for enzymes. Nucleic Acids Res. 2008, 36 (Database issue): D618-622.PubMed CentralPubMedGoogle Scholar