Impact of in vivo protein folding probability on local fitness landscapes

Abstract It is incompletely understood how biophysical properties like protein stability impact molecular evolution and epistasis. Epistasis is defined as specific when a mutation exclusively influences the phenotypic effect of another mutation, often at physically interacting residues. By contrast, nonspecific epistasis results when a mutation is influenced by a large number of non-local mutations. As most mutations are pleiotropic, the in vivo folding probability – governed by basal protein stability – is thought to determine activity-enhancing mutational tolerance, which implies that nonspecific epistasis is dominant. However, evidence exists for both specific and nonspecific epistasis as the prevalent factor, with limited comprehensive datasets to validate either claim. Using deep mutational scanning we probe how in vivo enzyme folding probability impacts local fitness landscapes. We computationally designed two different variants of the amidase AmiE in which catalytic efficiencies are statistically indistinguishable but the enzyme variants have lower probabilities of folding in vivo. Local fitness landscapes show only slight alterations among variants, with essentially the same global distribution of fitness effects. However, specific epistasis was predominant for the subset of mutations exhibiting positive sign epistasis. These mutations mapped to spatially distinct locations on AmiE near the initial mutation or proximal to the active site. Intriguingly, the majority of specific epistatic mutations were codon-dependent, with different synonymous codons resulting in fitness sign reversals. Together, these results offer a nuanced view of how protein folding probability impacts local fitness landscapes, and suggest that transcriptional-translational effects are an equally important determinant as stability in determining evolutionary outcomes.