Abstract:There is strong interest in accurate methods for predicting changes in protein stability resulting from amino acid mutations to the protein sequence. Recombinant proteins must often be stabilized to be used as therapeutics or reagents, and destabilizing mutations are implicated in a variety of diseases. Due to increased data availability and improved modeling techniques, recent studies have shown advancements in predicting changes in protein stability when a single point mutation is made. Less focus has been directed toward predicting changes in protein stability when there are two or more mutations, despite the significance of mutation clusters for disease pathways and protein design studies. Here, we analyze the largest available dataset of double point mutation stability and benchmark several widely used protein stability models on this and other datasets. We identify a blind spot in how predictors are typically evaluated on multiple mutations, finding that, contrary to assumptions in the field, current stability models are unable to consistently capture epistatic interactions between double mutations. We observe one notable deviation from this trend, which is that epistasis-aware models provide marginally better predictions on stabilizing double point mutations. We develop an extension of the ThermoMPNN framework for double mutant modeling as well as a novel data augmentation scheme which mitigates some of the limitations in available datasets. Collectively, our findings indicate that current protein stability models fail to capture the nuanced epistatic interactions between concurrent mutations due to several factors, including training dataset limitations and insufficient model sensitivity. Significance: Protein stability is governed in part by epistatic interactions between energetically coupled residues. Prediction of these couplings represents the next frontier in protein stability modeling. In this work, we benchmark protein stability models on a large dataset of double point mutations and identify previously overlooked limitations in model design and evaluation. We also introduce several new strategies to improve modeling of epistatic couplings between protein point mutations.

The Origin of Mutational Epistasis

The origin of mutational epistasis

Proteins Evolution Upon Point Mutations

Analysis of proteins in the light of mutations

How mutational epistasis impairs predictability in protein evolution and design

Physical principles underpinning molecular-level protein evolution

Addressing epistasis in the design of protein function

Emergent time scales of epistasis in protein evolution

Dynamic epistasis for different alleles of the same gene

The context-dependence of mutations: a linkage of formalisms

Biophysical inference of epistasis and the effects of mutations on protein stability and function

Coupled Molecular Dynamics Mediate Long- and Short-Range Epistasis Between Mutations That Affect Stability and Aggregation Kinetics

Protein stability models fail to capture epistatic interactions of double point mutations

Epistasis facilitates functional evolution in an ancient transcription factor

Divergent Folding-Mediated Epistasis Among Unstable Membrane Protein Variants

Importance of higher-order epistasis in large protein sequence-function relationships

Epistasis and Allele Specificity in the Emergence of a Stable Polymorphism in Escherichia coli

Exploring evolution to uncover insights into protein mutational stability

Global epistasis on fitness landscapes