Self-consistency test reveals systematic bias in programs for prediction change of stability upon mutation

Curating the gnomAD database: Report of novel variants in the globin-coding genes and bioinformatics analysis

Massive parallel sequencing technologies are facilitating the faster identification of sequence variants with the consequent capability of untangling the molecular bases of many human genetic syndromes. However, it is not always easy to understand the impact of novel variants, especially for missense changes, which can lead to a spectrum of phenotypes. This study presents a custom-designed multistep methodology to evaluate the impact of novel variants aggregated in the genome aggregation database for the HBB, HBA2, and HBA1 genes, by testing and improving its performance with a dataset of previously described alterations affecting those same genes. This approach scored high sensitivity and specificity values and showed an overall better performance than sequence-derived predictors, highlighting the importance of protein conformation and interaction specific analyses in curating variant databases. This study also describes the strengths and limitations of these structural studies and allows identifying residues in the globin chains more prone to tolerate substitutions.

DDGun: an untrained method for the prediction of protein stability changes upon single and multiple point variations

Background

Predicting the effect of single point variations on protein stability constitutes a crucial step toward understanding the relationship between protein structure and function. To this end, several methods have been developed to predict changes in the Gibbs free energy of unfolding (∆∆G) between wild type and variant proteins, using sequence and structure information. Most of the available methods however do not exhibit the anti-symmetric prediction property, which guarantees that the predicted ∆∆G value for a variation is the exact opposite of that predicted for the reverse variation, i.e., ∆∆G(A → B) = −∆∆G(B → A), where A and B are amino acids.

Results

Here we introduce simple anti-symmetric features, based on evolutionary information, which are combined to define an untrained method, DDGun (DDG untrained). DDGun is a simple approach based on evolutionary information that predicts the ∆∆G for single and multiple variations from sequence and structure information (DDGun3D). Our method achieves remarkable performance without any training on the experimental datasets, reaching Pearson correlation coefficients between predicted and measured ∆∆G values of ~ 0.5 and ~ 0.4 for single and multiple site variations, respectively. Surprisingly, DDGun performances are comparable with those of state of the art methods. DDGun also naturally predicts multiple site variations, thereby defining a benchmark method for both single site and multiple site predictors. DDGun is anti-symmetric by construction predicting the value of the ∆∆G of a reciprocal variation as almost equal (depending on the sequence profile) to -∆∆G of the direct variation. This is a valuable property that is missing in the majority of the methods.

Conclusions

Evolutionary information alone combined in an untrained method can achieve remarkably high performances in the prediction of ∆∆G upon protein mutation. Non-trained approaches like DDGun represent a valid benchmark both for scoring the predictive power of the individual features and for assessing the learning capability of supervised methods.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2923-1) contains supplementary material, which is available to authorized users.