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.

Divergent Folding-Mediated Epistasis Among Unstable Membrane Protein Variants

Single Mutation Effects On Conformational Change And Membrane Deformation Of Influenza Hemagglutinin Fusion Peptides

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

Even With Nonnative Interactions, The Updated Folding Transition States Of The Homologs Proteins G & L Are Extensive And Similar

Solvation and Desolvation Effects in Protein Folding: Native Flexibility, Kinetic Cooperativity and Enthalpic Barriers under Isostability Conditions

The Folding Transition State of Protein L is Extensive with Nonnative Interactions (and Not Small and Polarized)

The Intrinsic Folding Mechanism Affects the Coupled Folding-Binding Process of Unfolded Proteins

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

The origin of mutational epistasis

Interplaying Roles of Native Topology and Chain Length in Marginally Cooperative and Noncooperative Folding of Small Protein Fragments

Viral Receptor-Binding Protein Evolves New Function through Mutations That Cause Trimer Instability and Functional Heterogeneity

Protein folding and binding can emerge as evolutionary spandrels through structural coupling

A Localized Specific Interaction Alters the Unfolding Pathways of Structural Homologues.

Synonymous mutations can alter protein dimerization through localized interface misfolding involving self-entanglements

The role of hydrophobic matching on transmembrane helix packing in cells

Entropic Folding Pathway of Human Epidermal Growth Factor Explored by Disulfide Scrambling and Amplified Collective Motion Simulations.

The curious case of A31P, a topology-switching mutant of the Repressor of Primer protein : A molecular dynamics study of its folding and misfolding

A biophysical protein folding model accounts for most mutational fitness effects in viruses

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

The crystal structure of a mutant protein with altered but improved hydrophobic core packing