Cited by 202Machine Learning in BioinformaticsProtein Structure and DynamicsGenomics and Rare Diseases

MutaBind2: predicting the impacts of single and multiple mutations on protein-protein interactions

Open access
Zhang, Ning, Chen, Yuting, Lu, Haoyu, Zhao, Feiyang, Alvarez, Roberto Vera, Goncearenco, Alexander, Panchenko, Anna R., Li, Minghui·2020·iScience

Abstract

Missense mutations may affect proteostasis by destabilizing or over-stabilizing protein complexes and changing the pathway flux. Predicting the effects of stabilizing mutations on protein-protein interactions is notoriously difficult because existing experimental sets are skewed toward mutations reducing protein-protein binding affinity and many computational methods fail to correctly evaluate their effects. To address this issue, we developed a method MutaBind2, which estimates the impacts of single as well as multiple mutations on protein-protein interactions. MutaBind2 employs only seven features, and the most important of them describe interactions of proteins with the solvent, evolutionary conservation of the site, and thermodynamic stability of the complex and each monomer. This approach shows a distinct improvement especially in evaluating the effects of mutations increasing binding affinity. MutaBind2 can be used for finding disease driver mutations, designing stable protein complexes, and discovering new protein-protein interaction inhibitors.

AI snapshot

Written from the paper’s full text.

TL;DR

MutaBind2 is a structure-based method that predicts the effects of single and multiple mutations on protein-protein binding affinity using a seven-term scoring function. It improves performance especially for mutations that increase binding affinity, while requiring a 3D structure of the protein-protein complex.

Key findings

  • MutaBind2 uses only seven features, with the most important describing protein-solvent interactions, evolutionary conservation of the site, and thermodynamic stability of the complex and each monomer.
  • In cross-validation, MutaBind2 achieved Pearson R = 0.76 for single mutations and R = 0.74 for multiple mutations under leave-one-complex-out (CV4); under leave-one-binding-site-out (CV5), R = 0.69 for single and R = 0.71 for multiple mutations.
  • On the independent S1748 set, retraining on the SKEMPI set enriched with reverse mutations improved R from 0.38 to 0.63 and reduced RMSE from 1.51 to 1.25 kcal mol−1; for highly increasing mutations, MutaBind2 reached AUC = 0.86 versus MutaBind 0.65, BeAtMuSiC 0.56, and FoldX 0.74.
  • For multiple mutations, the worst performance was for three or more mutations on the same chain (R = 0.61 in CV4), and the best was for double mutations on multiple chains (R = 0.85 in CV4).

Methods

  • MutaBind2 employs a new minimization protocol and a scoring function with seven terms, trained on a larger SKEMPI2 dataset containing 1.7 times more mutations and 3.3 times more complexes than SKEMPI, and enriched with structural models of reverse mutations to balance mutations increasing binding affinity.
  • Performance was assessed with five cross-validation procedures (CV1–CV5), independent test sets including S1748, S4191, M1707, S487, and S8338, and comparisons with BeAtMuSiC, FoldX, iSEE, and mCSM-PPI2.
  • The paper also evaluated predictions against yeast two-hybrid data covering 4,797 mutation-interaction pairs, of which 451 mutations could be mapped to known 3D structures and 147 were interaction-disruptive; main metrics included Pearson R, RMSE, slope, ROC/AUC, MCC, sensitivity, and specificity.

Limitations

  • Stated: MutaBind2 requires the 3D structure of a protein-protein complex; six of seven features are calculated from that structure, limiting application to mutations that cannot be mapped to a structural complex.
  • Stated: multiple-mutation instances with more than 10 mutations are not supported because such cases are rare in the training data and prediction accuracy is low.
  • Evident from the text: outliers are more likely in complexes with higher protein-protein binding affinity and at mutated sites with more hydrogen bonds, and performance improves significantly after removing 5% outliers.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1016/j.isci.2020.100939

Cites

40 of 40 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Claims the paper makes

0 of 16 cited sources checked · 0 supported · 0 flagged

  • Findingown data

    MutaBind2 shows a distinct improvement especially in evaluating the effects of mutations increasing binding affinity.

    “This approach shows a distinct improvement especially in evaluating the effects of mutations increasing binding affinity.”
    • Also stated in: § Validation of MutaBind2 on Independent Test Sets · boosted
      “predictions for mutations increasing binding affinity were significantly improved without compromising the accuracy for predicting mutations decreasing binding affinity (Table 1).”
    Quote in the paperSubstantial claim
  • Findingcitation

    The outliers usually correspond to mutations with extreme experimental values.

    “Consistent with the observation of Pires et al. (Pires et al., 2016), the outliers usually correspond to those mutations with extreme experimental values (highly decreasing/increasing binding affinity)”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites Pires et al., 2016 source

      Not checked in the free run.

  • Findingreasoning

    MutaBind2 employs only seven features.

    “MutaBind2 employs only seven features”
    • Also stated in: § Results and Discussion · plain
      “we designed a scoring function with seven terms instead of 10 used in the previous MutaBind version”
    Quote in the paperSubstantial claim
  • Findingown data

    MutaBind and MutaBind2 show excellent performance in distinguishing interaction-disruptive from other mutations.

    “Figures 3 and S5 show excellent performance of MutaBind and MutaBind2 in distinguishing interaction-disruptive from other mutations.”
    Quote in the paperSubstantial claim
  • Findingown data

    MutaBind2's performance is improved significantly after removing 5% outliers.

    “Figure S11A shows that the performance is improved significantly after removing 5% outliers.”
    Quote in the paperNot a substantial claim
  • May combine two claims. 14
  • May combine two claims.

    Findingown data

    The correlation coefficients of each cross-validation round exceed 0.80 for CV1 and CV2 and about 0.70 for CV3 for both single and multiple mutations.

    “The correlation coefficients of each cross-validation round exceed 0.80 for “CV1” and “CV2” and about 0.70 for “CV3” cross-validation for both single and multiple mutations.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingreasoning

    MutaBind and MutaBind2, unlike other methods, estimate interactions of a protein with the solvent, which is one of the most important terms together with the site evolutionary conservation and thermodynamic stability of a protein complex and each binding partner.

    “Importantly, MutaBind and MutaBind2, unlike other methods mentioned earlier, estimate interactions of a protein with the solvent, which is one of the most important terms together with the site evolutionary conservation and thermodynamic stability of a protein complex and each binding partner (Table S2).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    Outliers are more likely to appear in complexes with higher protein-protein binding affinity and at mutated sites with higher number of hydrogen bonds.

    “we found that the outliers are more likely to appear in complexes with higher protein-protein binding affinity and at the mutated sites with higher number of hydrogen bonds (Figure S11B).”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    The Pearson correlation coefficient between experimental and computed ΔΔG values using CV4 reaches 0.76 for single mutations and 0.74 for multiple mutations.

    “The Pearson correlation coefficient between experimental and computed ΔΔG values using the “leave-one-complex-out” (“CV4”) procedure reaches 0.76 for single mutations and 0.74 for multiple mutations (Table 1 and Figure 2).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    MutaBind2 performed well for almost all categories of multiple mutations; the worst performance was for three or more mutations on the same chain (R = 0.61 in CV4), and the best was for double mutations on multiple chains (R = 0.85 in CV4).

    “We found that MutaBind2 performed well for almost all categories; the worst performance was observed for three or more mutations introduced on the same chain (R = 0.61 in CV4 validation), and the best performance was achieved for double mutations on multiple chains (R = 0.85 in CV4 validation).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The CV5 correlation coefficient is 0.69 for single mutations and 0.71 for multiple mutations.

    “the correlation coefficient remained statistically significant, being equal to 0.69 for single mutations and 0.71 for multiple mutations (Table 1 and Figure S4D).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    MutaBind2 improves the performance for both interfacial and non-interfacial mutations compared with the previous version and outperforms other methods on the S1748 set.

    “MutaBind2 further improves the performance for both interfacial and non-interfacial mutations compared with the previous version and outperforms other methods on the S1748 set (Figure S7).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    MutaBind2 was trained on a larger dataset from SKEMPI2, encompassing 1.7 times more mutations and 3.3 times more complexes compared with SKEMPI.

    “we trained our models on a much larger dataset from SKEMPI2, which encompassed 1.7 times more mutations and 3.3 times more complexes compared with SKEMPI used for training of MutaBind.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The most important features of MutaBind2 describe interactions of proteins with the solvent, evolutionary conservation of the site, and thermodynamic stability of the complex and each monomer.

    “the most important of them describe interactions of proteins with the solvent, evolutionary conservation of the site, and thermodynamic stability of the complex and each monomer.”
    • Also stated in: § Prediction of Mutations Highly Decreasing and Increasing Binding Affinity · plain
      “MutaBind and MutaBind2, unlike other methods mentioned earlier, estimate interactions of a protein with the solvent, which is one of the most important terms together with the site evolutionary conservation and thermodynamic stability of a protein complex and each binding partner (Table S2).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    Enriching the training set with mutations increasing binding affinity improved performance on an independent S1748 set, with correlation coefficient increasing from 0.38 to 0.63 and RMSE decreasing from 1.51 to 1.25 kcal mol−1.

    “the performance on an independent S1748 set improved considerably (the correlation coefficient increased from 0.38 to 0.63, and root-mean-square error decreased from 1.51 to 1.25 kcal mol−1).”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The methodologies of training procedures employed by many computational methods rely on experimental sets enriched with mutations that decrease binding.
  • May combine two claims.

    Findingown data

    Predictions for mutations increasing binding affinity were significantly improved without compromising the accuracy for predicting mutations decreasing binding affinity.

    “Moreover, predictions for mutations increasing binding affinity were significantly improved without compromising the accuracy for predicting mutations decreasing binding affinity (Table 1).”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The methodologies of training procedures employed by many computational methods rely on experimental sets enriched with mutations that decrease binding.
  • May combine two claims.

    Premisecitation

    Missense mutations may affect proteostasis by destabilizing or over-stabilizing protein complexes and changing the pathway flux.

    “Missense mutations may affect proteostasis by destabilizing or over-stabilizing protein complexes and changing the pathway flux.”
    • Also stated in: § Introduction · hedged
      “However, the latter may also have important consequences leading to disruption of proteostasis by over-stabilizing transient protein complexes (Nishi et al., 2013, Stefl et al., 2013, Rutten et al., 2018, Shoichet et al., 1995, Nagatani et al., 2007, Witham et al., 2011, Jubb et al., 2016) and changing the pathway flux.”
    Quote in the paperSubstantial claim0 of 3 sources checked
    • cites Nishi et al., 2013

      Not checked in the free run.

    • cites Stefl et al., 2013

      Not checked in the free run.

    • cites Rutten et al., 2018

      Not checked in the free run.

  • May combine two claims.

    Interpretationreasoning

    MutaBind2 can be used for finding disease driver mutations, designing stable protein complexes, and discovering new protein-protein interaction inhibitors.

    “MutaBind2 can be used for finding disease driver mutations, designing stable protein complexes, and discovering new protein-protein interaction inhibitors.”
    • Also stated in: § Introduction · plain
      “MutaBind2 can be applied to a large number of tasks, including, but not limited to, finding disease driver mutations and understanding their molecular mechanisms, assessing the effects of sequence variants on protein fitness, structural modeling of mutant complexes, and designing protein interaction inhibitors (Goncearenco et al., 2017).”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites Goncearenco et al., 2017

      Not checked in the free run.

    • Rests on an assumption: Large-scale rapid experimental assays that would allow the assessment of thousands of variants are still limited.
  • May combine two claims.

    Premisecitation

    The previous version MutaBind could predict single mutations highly decreasing binding affinity relatively well but failed to annotate mutations highly increasing affinity.

    “The previous version, MutaBind, could predict single mutations highly decreasing binding affinity relatively well but failed to annotate mutations highly increasing affinity.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites Li et al., 2016b

      Not checked in the free run.

    • Rests on an assumption: The methodologies of training procedures employed by many computational methods rely on experimental sets enriched with mutations that decrease binding.
14 of 14
  • The paper relies on this without support.

    Proteins and their binding interfaces tend to be optimized in evolution with regard to their stability.

    “proteins and their binding interfaces tend to be optimized in evolution with regard to their stability”
  • The paper relies on this without support.

    The methodologies of training procedures employed by many computational methods rely on experimental sets enriched with mutations that decrease binding.

    “the methodologies of training procedures employed by many computational methods rely on these experimental sets”
  • The paper relies on this without support.

    Large-scale rapid experimental assays that would allow the assessment of thousands of variants are still limited.

    “large-scale rapid experimental assays that would allow the assessment of thousands of variants are still limited”
  • The paper relies on this without support.

    The majority of mutations produce rather subtle effects with unclear clinical significance.

    “the majority of mutations produce rather subtle effects with unclear clinical significance”
4 of 4

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Cites

Flagged first, then by citations

Cited by

Sorted by citations

Authors
Zhang, N., Chen, Y., Lu, H., Zhao, F., Alvarez, R. V., Goncearenco, A., Panchenko, A. R., & Li, M.
Affiliations
  • Soochow University (CN)
  • National Institutes of Health (US)
Corresponding author
  • Panchenko, Anna R.
  • Li, Minghui
Venue
iScience
Publisher
Elsevier BV
Published
2020-03
Volume
23
Issue
3
Pages
100939
ISSN
2589-0042 (ISSN-L)
OpenAlex
W3007500173
Cited by
199
References
40
Citations by year
2026: 15 · 2025: 35 · 2024: 41 · 2023: 45 · 2022: 29 · 2021: 25 · 2020: 8 · 2012: 1
FWCI
8.68
Type
journal article
Language
en
Open access
Locations
  • iScience (publishedVersion, cc-by-nc-nd)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by-nc-nd)
  • PubMed Central (publishedVersion, cc-by-nc-nd)
  • PubMed (publishedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
Funders
  • Natural Science Foundation of Jiangsu Province
  • National Natural Science Foundation of China
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Machine Learning in Bioinformatics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Protein Structure and Dynamics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genomics and Rare Diseases
Keywords
protein-protein interaction, missense mutations, binding affinity prediction, stabilizing mutations, multiple mutations, protein complex stability, evolutionary conservation, solvent accessibility, protein-protein interaction inhibitors, proteostasis

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Claims read by deepseek-v4.1-flash under definition version 3 (extract-v3) on 2026-10-10.

Classified by jev-1.13-20260917 (gold-guideline-v0.2) on 2026-10-10.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review