SABER: a computational method for identifying active sites for new reactions

A software suite, SABER (Selection of Active/Binding sites for Enzyme Redesign), has been developed for the analysis of atomic geometries in protein structures, using a geometric hashing algorithm (Barker and Thornton, Bioinformatics 2003;19:1644–1649). SABER is used to explore the Protein Data Bank (PDB) to locate proteins with a specific 3D arrangement of catalytic groups to identify active sites that might be redesigned to catalyze new reactions. As a proof‐of‐principle test, SABER was used to identify enzymes that have the same catalytic group arrangement present in o‐succinyl benzoate synthase (OSBS). Among the highest‐scoring scaffolds identified by the SABER search for enzymes with the same catalytic group arrangement as OSBS were L ‐Ala D/L ‐Glu epimerase (AEE) and muconate lactonizing enzyme II (MLE), both of which have been redesigned to become effective OSBS catalysts, demonstrated by experiments. Next, we used SABER to search for naturally existing active sites in the PDB with catalytic groups similar to those present in the designed Kemp elimination enzyme KE07. From over 2000 geometric matches to the KE07 active site, SABER identified 23 matches that corresponded to residues from known active sites. The best of these matches, with a 0.28 Å catalytic atom RMSD to KE07, was then redesigned to be compatible with the Kemp elimination using RosettaDesign. We also used SABER to search for potential Kemp eliminases using a theozyme predicted to provide a greater rate acceleration than the active site of KE07, and used Rosetta to create a design based on the proteins identified.     

Evolutionary Optimization of Computationally Designed Enzymes: Kemp Eliminases of the KE07 Series

Understanding enzyme catalysis through the analysis of natural enzymes is a daunting challenge—their active sites are complex and combine numerous interactions and catalytic forces that are finely coordinated. Study of more rudimentary (wo)man-made enzymes provides a unique opportunity for better understanding of enzymatic catalysis. KE07, a computationally designed Kemp eliminase that employs a glutamate side chain as the catalytic base for the critical proton abstraction step and an apolar binding site to guide substrate binding, was optimized by seven rounds of random mutagenesis and selection, resulting in a > 200-fold increase in catalytic efficiency. Here, we describe the directed evolution process in detail and the biophysical and crystallographic studies of the designed KE07 and its evolved variants. The optimization of KE07's activity to give a k_cat/K_M value of ∼ 2600 s^− 1 M^− 1 and an ∼ 10⁶-fold rate acceleration (k_cat/k_uncat) involved the incorporation of up to eight mutations. These mutations led to a marked decrease in the overall thermodynamic stability of the evolved KE07s and in the configurational stability of their active sites. We identified two primary contributions of the mutations to KE07's improved activity: (i) the introduction of new salt bridges to correct a mistake in the original design that placed a lysine for leaving-group protonation without consideration of its “quenching” interactions with the catalytic glutamate, and (ii) the tuning of the environment, the pK_a of the catalytic base, and its interactions with the substrate through the evolution of a network of hydrogen bonds consisting of several charged residues surrounding the active site.     

Optimization of the in-silico-designed kemp eliminase KE70 by computational design and directed evolution

Although de novo computational enzyme design has been shown to be feasible, the field is still in its infancy: the kinetic parameters of designed enzymes are still orders of magnitude lower than those of naturally occurring ones. Nonetheless, designed enzymes can be improved by directed evolution, as recently exemplified for the designed Kemp eliminase KE07. Random mutagenesis and screening resulted in variants with > 200-fold higher catalytic efficiency and provided insights about features missing in the designed enzyme. Here we describe the optimization of KE70, another designed Kemp eliminase. Amino acid substitutions predicted to improve catalysis in design calculations involving extensive backbone sampling were individually tested. Those proven beneficial were combinatorially incorporated into the originally designed KE70 along with random mutations, and the resulting libraries were screened for improved eliminase activity. Nine rounds of mutation and selection resulted in > 400-fold improvement in the catalytic efficiency of the original KE70 design, reflected in both higher k_cat values and lower K_m values, with the best variants exhibiting k_cat/K_m values of > 5 × 10⁴ s^− 1 M^− 1. The optimized KE70 variants were characterized structurally and biochemically, providing insights into the origins of the improvements in catalysis. Three primary contributions were identified: first, the reshaping of the active-site cavity to achieve tighter substrate binding; second, the fine-tuning of electrostatics around the catalytic His-Asp dyad; and, third, the stabilization of the active-site dyad in a conformation optimal for catalysis.     

Challenges and advances in validating enzyme design proposals: the case of kemp eliminase catalysis

One of the fundamental challenges in biotechnology and biochemistry is the ability to design effective enzymes. Despite recent progress, most of the advances on this front have been made by placing the reacting fragments in the proper places, rather than by optimizing the preorganization of the environment, which is the key factor in enzyme catalysis. Thus, rational improvement of the preorganization would require approaches capable of evaluating reliably the actual catalytic effect. This work considers the catalytic effects in different Kemp eliminases as a benchmark for a computer-aided enzyme design. It is shown that the empirical valence bond provides a powerful screening tool, with significant advantages over current alternative strategies. The insights provided by the empirical valence bond calculations are discussed with an emphasis on the ability to analyze the difference between the linear free energy relationships obtained in solution and those found in the enzymes. We also point out the trade-off between the reliability and speed of the calculations and try to determine what it takes to realize reliable computer-aided screening.     

How protein stability and new functions trade off

Numerous studies have noted that the evolution of new enzymatic specificities is accompanied by loss of the protein's thermodynamic stability (DeltaDeltaG), thus suggesting a tradeoff between the acquisition of new enzymatic functions and stability. However, since most mutations are destabilizing (DeltaDeltaG>0), one should ask how destabilizing mutations that confer new or altered enzymatic functions relative to all other mutations are. We applied DeltaDeltaG computations by FoldX to analyze the effects of 548 mutations that arose from the directed evolution of 22 different enzymes. The stability effects, location, and type of function-altering mutations were compared to DeltaDeltaG changes arising from all possible point mutations in the same enzymes. We found that mutations that modulate enzymatic functions are mostly destabilizing (average DeltaDeltaG = +0.9 kcal/mol), and are almost as destabilizing as the "average" mutation in these enzymes (+1.3 kcal/mol). Although their stability effects are not as dramatic as in key catalytic residues, mutations that modify the substrate binding pockets, and thus mediate new enzymatic specificities, place a larger stability burden than surface mutations that underline neutral, non-adaptive evolutionary changes. How are the destabilizing effects of functional mutations balanced to enable adaptation? Our analysis also indicated that many mutations that appear in directed evolution variants with no obvious role in the new function exert stabilizing effects that may compensate for the destabilizing effects of the crucial function-altering mutations. Thus, the evolution of new enzymatic activities, both in nature and in the laboratory, is dependent on the compensatory, stabilizing effect of apparently "silent" mutations in regions of the protein that are irrelevant to its function.     

Directed evolution: an approach to engineer enzymes

Directed evolution is being used increasingly in industrial and academic laboratories to modify and improve commercially important enzymes. Laboratory evolution is thought to make its biggest contribution in explorations of non-natural functions, by allowing us to distinguish the properties nurtured by evolution. In this review we report the significant advances achieved with respect to the methods of biocatalyst improvement and some critical properties and applications of the modified enzymes. The application of directed evolution has been elaborately demonstrated for protein solubility, stability and catalytic efficiency. Modification of certain enzymes for their application in enantioselective catalysis has also been elucidated. By providing a simple and reliable route to enzyme improvement, directed evolution has emerged as a key technology for enzyme engineering and biocatalysis.     

Directed evolution of a new catalytic site in 2-keto-3-deoxy-6-phosphogluconate aldolase from Escherichia coli

BACKGROUND: Aldolases are carbon bond-forming enzymes that have long been identified as useful tools for the organic chemist. However, their utility is limited in part by their narrow substrate utilization. Site-directed mutagenesis of various enzymes to alter their specificity has been performed for many years, typically without the desired effect. More recently directed evolution has been employed to engineer new activities onto existing scaffoldings. This approach allows random mutation of the gene and then selects for fitness to purpose those proteins with the desired activity. To date such approaches have furnished novel activities through multiple mutations of residues involved in recognition; in no instance has a key catalytic residue been altered while activity is retained. RESULTS: We report a double mutant of E. coli 2-keto-3-deoxy-6-phosphogluconate aldolase with reduced but measurable enzyme activity and a synthetically useful substrate profile. The mutant was identified from directed-evolution experiments. Modification of substrate specificity is achieved by altering the position of the active site lysine from one beta strand to a neighboring strand rather than by modification of the substrate recognition site. The new enzyme is different to all other existing aldolases with respect to the location of its active site to secondary structure. The new enzyme still displays enantiofacial discrimination during aldol addition. We have determined the crystal structure of the wild-type enzyme (by multiple wavelength methods) to 2.17 A and the double mutant enzyme to 2.7 A resolution. CONCLUSIONS: These results suggest that the scope of directed evolution is substantially larger than previously envisioned in that it is possible to perturb the active site residues themselves as well as surrounding loops to alter specificity. The structure of the double mutant shows how catalytic competency is maintained despite spatial reorganization of the active site with respect to substrate.     

Directed Evolution: An Approach to Engineer Enzymes

ABSTRACT

Directed evolution is being used increasingly in industrial and academic laboratories to modify and improve commercially important enzymes. Laboratory evolution is thought to make its biggest contribution in explorations of non-natural functions, by allowing us to distinguish the properties nurtured by evolution. In this review we report the significant advances achieved with respect to the methods of biocatalyst improvement and some critical properties and applications of the modified enzymes. The application of directed evolution has been elaborately demonstrated for protein solubility, stability and catalytic efficiency. Modification of certain enzymes for their application in enantioselective catalysis has also been elucidated. By providing a simple and reliable route to enzyme improvement, directed evolution has emerged as a key technology for enzyme engineering and biocatalysis.     

Improving upon nature: active site remodeling produces highly efficient aldolase activity toward hydrophobic electrophilic substrates

The substrate specificity of enzymes is frequently narrow and constrained by multiple interactions, limiting the use of natural enzymes in biocatalytic applications. Aldolases have important synthetic applications, but the usefulness of these enzymes is hampered by their narrow reactivity profile with unnatural substrates. To explore the determinants of substrate selectivity and alter the specificity of Escherichia coli 2-keto-3-deoxy-6-phosphogluconate (KDPG) aldolase, we employed structure-based mutagenesis coupled with library screening of mutant enzymes localized to the bacterial periplasm. We identified two active site mutations (T161S and S184L) that work additively to enhance the substrate specificity of this aldolase to include catalysis of retro-aldol cleavage of (4S)-2-keto-4-hydroxy-4-(2'-pyridyl)butyrate (S-KHPB). These mutations improve the value of k(cat)/K(M)(S-KHPB) by >450-fold, resulting in a catalytic efficiency that is comparable to that of the wild-type enzyme with the natural substrate while retaining high stereoselectivity. Moreover, the value of k(cat)(S-KHPB) for this mutant enzyme, a parameter critical for biocatalytic applications, is 3-fold higher than the maximal value achieved by the natural aldolase with any substrate. This mutant also possesses high catalytic efficiency for the retro-aldol cleavage of the natural substrate, KDPG, and a >50-fold improved activity for cleavage of 2-keto-4-hydroxy-octonoate, a nonfunctionalized hydrophobic analogue. These data suggest a substrate binding mode that illuminates the origin of facial selectivity in aldol addition reactions catalyzed by KDPG and 2-keto-3-deoxy-6-phosphogalactonate aldolases. Furthermore, targeting mutations to the active site provides a marked improvement in substrate selectivity, demonstrating that structure-guided active site mutagenesis combined with selection techniques can efficiently identify proteins with characteristics that compare favorably to those of naturally occurring enzymes.     

Catalytic and binding poly-reactivities shared by two unrelated proteins: The potential role of promiscuity in enzyme evolution

It is generally accepted that enzymes evolved via gene duplication of existing proteins. But duplicated genes can serve as a starting point for the evolution of a new function only if the protein they encode happens to exhibit some activity towards this new function. Although the importance of such catalytic promiscuity in enzyme evolution has been proposed, little is actually known regarding how common promiscuous catalytic activities are in proteins or their origins, magnitudes, and potential contribution to the survival of an organism. Here we describe a pattern of promiscuous activities in two completely unrelated proteins-serum albumins and a catalytic antibody (aldolase antibody 38C2). Despite considerable structural dissimilarities-in the shape of the cavities and the position of catalytic lysine residues-both active sites are able to catalyze the Kemp elimination, a model reaction for proton transfer from carbon. We also show that these different active sites can bind promiscuously an array of hydrophobic negatively charged ligands. We suggest that the basic active-site features of an apolar pocket and a lysine residue can act as a primitive active site allowing these promiscuous activities to take place. We also describe, by modelling product formation at different substrate concentrations, how promiscuous activities of this kind- inefficient and rudimentary as they are-can provide a considerable selective advantage and a starting point for the evolution of new functions.     

Rational engineering of enzyme allosteric regulation through sequence evolution analysis

Control of enzyme allosteric regulation is required to drive metabolic flux toward desired levels. Although the three-dimensional (3D) structures of many enzyme-ligand complexes are available, it is still difficult to rationally engineer an allosterically regulatable enzyme without decreasing its catalytic activity. Here, we describe an effective strategy to deregulate the allosteric inhibition of enzymes based on the molecular evolution and physicochemical characteristics of allosteric ligand-binding sites. We found that allosteric sites are evolutionarily variable and comprised of more hydrophobic residues than catalytic sites. We applied our findings to design mutations in selected target residues that deregulate the allosteric activity of fructose-1,6-bisphosphatase (FBPase). Specifically, charged amino acids at less conserved positions were substituted with hydrophobic or neutral amino acids with similar sizes. The engineered proteins successfully diminished the allosteric inhibition of E. coli FBPase without affecting its catalytic efficiency. We expect that our method will aid the rational design of enzyme allosteric regulation strategies and facilitate the control of metabolic flux.     

Enhancing the Promiscuous Phosphotriesterase Activity of a Thermostable Lactonase (GkaP) for the Efficient Degradation of Organophosphate Pesticides

The phosphotriesterase-like lactonase (PLL) enzymes in the amidohydrolase superfamily hydrolyze various lactones and exhibit latent phosphotriesterase activities. These enzymes serve as attractive templates for in vitro evolution of neurotoxic organophosphates (OPs) with hydrolytic capabilities that can be used as bioremediation tools. Here, a thermostable PLL from Geobacillus kaustophilus HTA426 (GkaP) was targeted for joint laboratory evolution with the aim of enhancing its catalytic efficiency against OP pesticides. By a combination of site saturation mutagenesis and whole-gene error-prone PCR approaches, several improved variants were isolated. The most active variant, 26A8C, accumulated eight amino acid substitutions and demonstrated a 232-fold improvement over the wild-type enzyme in reactivity (k(cat)/K(m)) for the OP pesticide ethyl-paraoxon. Concomitantly, this variant showed a 767-fold decrease in lactonase activity with δ-decanolactone, imparting a specificity switch of 1.8 × 10(5)-fold. 26A8C also exhibited high hydrolytic activities (19- to 497-fold) for several OP pesticides, including parathion, diazinon, and chlorpyrifos. Analysis of the mutagenesis sites on the GkaP structure revealed that most mutations are located in loop 8, which determines substrate specificity in the amidohydrolase superfamily. Molecular dynamics simulation shed light on why 26A8C lost its native lactonase activity and improved the promiscuous phosphotriesterase activity. These results permit us to obtain further insights into the divergent evolution of promiscuous enzymes and suggest that laboratory evolution of GkaP may lead to potential biological solutions for the efficient decontamination of neurotoxic OP compounds.     

The catalytic effect of dihydrofolate reductase and its mutants is determined by reorganization energies

The effect of distant mutations on the catalytic reaction of dihydrofolate reductase (DHFR) is reexamined by empirical valence bond simulations. The simulations reproduce for the first time the changes in the observed rate constants (without the use of adjustable parameters for this purpose) and show that the changes in activation barriers are strongly correlated with the corresponding changes in the reorganization energy. The preorganization of the polar groups of enzymes is the key catalytic factor, and anticatalytic mutations destroy this preorganization. Some anticatalytic mutations in DHFR also increase the distance between the donor and acceptor, but this effect is not directly related to catalysis since the native enzyme and the uncatalyzed reaction in water have similar average donor-acceptor distances. Insight into the effect of a mutation is provided by constructing the relevant free energy surfaces in terms of the generalized solute-solvent coordinates. It is shown how the mutations change the reaction coordinate and the activation barrier, and it is clarified that the corresponding changes do not reflect dynamical effects. It is also pointed out that all reactions in a condensed phase involve correlated motions (both in enzymes and in solution) and that the change of such motions upon mutations is a result of the change in the shape of the multidimensional reaction path on the solute-solvent surface, rather than the reason for the change in rate constant. Thus, as far as catalysis is concerned, the change in the activation barrier is due to the change in the electrostatic preorganization energy.     

Rationally selected single-site mutants of the Thermoascus aurantiacus endoglucanase increase hydrolytic activity on cellulosic substrates

Variants of the Thermoascus aurantiacus Eg1 enzyme with higher catalytic efficiency than wild-type were obtained via site-directed mutagenesis. Using a rational mutagenesis approach based on structural bioinformatics and evolutionary analysis, two positions (F16S and Y95F) were identified as priority sites for mutagenesis. The mutant and parent enzymes were expressed and secreted from Pichia pastoris and the single site mutants F16S and Y95F showed 1.7- and 4.0-fold increases in k(cat) and 1.5- and 2.5-fold improvements in hydrolytic activity on cellulosic substrates, respectively, while maintaining thermostability. Similar to the parent enzyme, the two variants were active between pH 4.0 and 8.0 and showed optimal activity at temperature 70°C at pH 5.0. The purified enzymes were active at 50°C for over 12 h and retained at least 80% of initial activity for 2 h at 70°C. In contrast to the improved hydrolysis seen with the single mutation enzymes, no improvement was observed with a third variant carrying a combination of both mutations, which instead showed a 60% reduction in catalytic efficiency. This work further demonstrates that non-catalytic amino acid residues can be engineered to enhance catalytic efficiency in pretreatment enzymes of interest.     

Force-clamp spectroscopy detects residue co-evolution in enzyme catalysis

Understanding how the catalytic mechanisms of enzymes are optimized through evolution remains a major challenge in molecular biology. The concept of co-evolution implicates that compensatory mutations occur to preserve the structure and function of proteins. We have combined statistical analysis of protein sequences with the sensitivity of single molecule force-clamp spectroscopy to probe how catalysis is affected by structurally distant correlated mutations in Escherichia coli thioredoxin. Our findings show that evolutionary anti-correlated mutations have an inhibitory effect on enzyme catalysis, whereas positively correlated mutations rescue the catalytic activity. We interpret these results in terms of an evolutionary tuning of both the enzyme-substrate binding process and the chemistry of the active site. Our results constitute a direct observation of distant residue co-evolution in enzyme catalysis.     

The evolution of catalytic efficiency and substrate promiscuity in human theta class 1-1 glutathione transferase

Theta class glutathione transferases (GST) from various species exhibit markedly different catalytic activities in conjugating the tripeptide glutathione (GSH) to a variety of electrophilic substrates. For example, the human theta 1-1 enzyme (hGSTT1-1) is 440-fold less efficient than the rat theta 2-2 enzyme (rGSTT2-2) with the fluorogenic substrate 7-amino-4-chloromethyl coumarin (CMAC). Large libraries of hGSTT1-1 constructed by error-prone PCR, DNA shuffling, or saturation mutagenesis were screened for improved catalytic activity towards CMAC in a quantitative fashion using flow cytometry. An iterative directed evolution approach employing random mutagenesis in conjunction with homologous recombination gave rise to enzymes exhibiting up to a 20,000-fold increase in k(cat)/K(M) compared to hGSTT1-1. All highly active clones encoded one or more mutations at residues 32, 176, or 234. Combinatorial saturation mutagenesis was used to evaluate the full complement of natural amino acids at these positions, and resulted in the isolation of enzymes with catalytic rates comparable to those exhibited by the fastest mutants obtained via directed evolution. The substrate selectivities of enzymes resulting from random mutagenesis, DNA shuffling, and combinatorial saturation mutagenesis were evaluated using a series of distinct electrophiles. The results revealed that promiscuous substrate activities arose in a stochastic manner, as they did not correlate with catalytic efficiency towards the CMAC selection substrate. In contrast, chimeric enzymes previously constructed by homology-independent recombination of hGSTT-1 and rGSTT2-2 exhibited very different substrate promiscuity profiles, and showed a more defined relationship between evolved and promiscuous activities.     

Strategy and success for the directed evolution of enzymes

Directed evolution is widely used to improve enzymes, particularly for industrial biocatalytic processes. Molecular biology advances present many new strategies for directed evolution. Commonly used techniques have led to many successful examples of enzyme improvement, yet there is still a need to improve both the efficiency and capability of directed evolution. Recent strategies aimed at making directed evolution faster and more efficient take better advantage of available structural and sequence information. The underlying principles that lead to early dead-ends for directed evolution experiments are also discussed along with recent strategies designed to by-pass them. Several emerging methods for creating novel enzymes are also discussed including examples of catalytic activity for which there is no precedent in nature. Finally, the combined use of several strategies is likely to be required in practice to improve multiple target properties of an enzyme, as successfully shown by a recent industrial example.