A combinatorially complete epistatic fitness landscape in an enzyme active site

Protein engineering often targets binding pockets or active sites which are enriched in epistasis-- non-additive interactions between amino acid substitutions--and where the combined effects of multiple single substitutions are difficult to predict. Few existing sequence-fitness datasets capture epistasis at large scale, especially for enzyme catalysis, limiting the development and assessment of model-guided enzyme engineering approaches. We present here a combinatorially complete, 160,000-variant fitness landscape across four residues in the active site of an enzyme. Assaying the native reaction of a thermostable {beta}-subunit of tryptophan synthase (TrpB) in a non-native environment yielded a landscape characterized by significant epistasis and many local optima. These effects prevent simulated directed evolution approaches from efficiently reaching the global optimum. There is nonetheless wide variability in the effectiveness of different directed evolution approaches, which together provide experimental benchmarks for computational and machine learning workflows. The most-fit TrpB variants contain a substitution that is nearly absent in natural TrpB sequences--a result that conservation-based predictions would not capture. Thus, although fitness prediction using evolutionary data can enrich in more-active variants, these approaches struggle to identify and differentiate among the most-active variants, even for this near-native function. Overall, this work presents a new, large-scale testing ground for model-guided enzyme engineering and suggests that efficient navigation of epistatic fitness landscapes can be improved by advances in both machine learning and physical modeling. Significance statementPredictive models for protein engineering seek to capture the relationship between protein sequence and function. While many methods and datasets exist for predicting the effects of single substitutions across a range of protein functions, fewer capture interactions among substitutions, which are much more difficult to predict. Even fewer do this comprehensively for a catalytic function. To provide a testbed for evaluating predictive models for enzyme engineering, we constructed and analyzed a 160,000-member enzyme sequence-fitness dataset at four interacting residues near the active site of tryptophan synthase, capturing significant non-additive effects of substitutions on catalytic function. It is necessary to predict and understand such interactions in order to efficiently traverse evolutionary landscapes and build machine learning models that accelerate protein engineering.