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.

How mutational epistasis impairs predictability in protein evolution and design

Addressing epistasis in the design of protein function

The origin of mutational epistasis

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

Dynamic epistasis for different alleles of the same gene

Epistasis arises from shifting the rate-limiting step during enzyme evolution of a β-lactamase

Network of epistatic interactions in an enzyme active site revealed by large-scale deep mutational scanning

Emergent time scales of epistasis in protein evolution

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

Epistasis facilitates functional evolution in an ancient transcription factor

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

Strong Selection Significantly Increases Epistatic Interactions in the Long-Term Evolution of a Protein

Higher-order epistasis and phenotypic prediction

Modular epistasis and the compensatory evolution of gene deletion mutants

The simplicity of protein sequence-function relationships

A mechanistic view of enzyme evolution

Reducing Epistasis and Pleiotropy Can Avoid the Survival of the Flattest Tragedy