The design of new enzymes for reactions not catalysed by naturally occurring biocatalysts is a challenge for protein engineering and is a critical test of our understanding of enzyme catalysis. Here we describe the computational design of eight enzymes that use two different catalytic motifs to catalyse the Kemp elimination-a model reaction for proton transfer from carbon-with measured rate enhancements of up to 10 5 and multiple turnovers. Mutational analysis confirms that catalysis depends on the computationally designed active sites, and a high-resolution crystal structure suggests that the designs have close to atomic accuracy. Application of in vitro evolution to enhance the computational designs produced a .200-fold increase in k cat /K m (k cat /K m of 2,600 M 21 s 21 and k cat /k uncat of .10 6 ). These results demonstrate the power of combining computational protein design with directed evolution for creating new enzymes, and we anticipate the creation of a wide range of useful new catalysts in the future.Naturally occurring enzymes are extraordinarily efficient catalysts 1 . They bind their substrates in a well-defined active site with precisely aligned catalytic residues to form highly active and selective catalysts for a wide range of chemical reactions under mild conditions. Nevertheless, many important synthetic reactions lack a naturally occurring enzymatic counterpart. Hence, the design of stable enzymes with new catalytic activities is of great practical interest, with potential applications in biotechnology, biomedicine and industrial processes. Furthermore, the computational design of new enzymes provides a stringent test of our understanding of how naturally occurring enzymes work. In the past several years, there has been exciting progress in designing new biocatalysts 2,3 .Here we describe the use of our recently developed computational enzyme design methodology 4 to create new enzyme catalysts for a reaction for which no naturally occurring enzyme exists: the Kemp elimination 5,6 . The reaction, shown in Fig. 1a, has been extensively studied as an activated model system for understanding the catalysis of proton abstraction from carbon-a process that is normally restricted by high activation-energy barriers 7,8 . Computational design methodThe first step in our protocol for designing new enzymes is to choose a catalytic mechanism and then to use quantum mechanical transition state calculations to create an idealized active site with protein functional groups positioned so as to maximize transition state stabilization (Fig. 1b). The key step for the Kemp elimination is deprotonation of a carbon by a general base. We chose two different catalytic bases for this purpose: first, the carboxyl group of an aspartate or glutamate side chain, and, second, the imidazole of a histidine positioned and polarized by the carboxyl group of an aspartate or glutamate (we refer to this combination as a His-Asp dyad). The two choices have complementary strengths and weaknesses. The advantage of the carboxylate...

The creation of enzymes capable of catalyzing any desired chemical reaction is a grand challenge for computational protein design. Using new algorithms that rely on hashing techniques to construct active sites for multistep reactions, we designed retro-aldolases that use four different catalytic motifs to catalyze the breaking of a carbon-carbon bond in a nonnatural substrate. Of the 72 designs that were experimentally characterized, 32, spanning a range of protein folds, had detectable retro-aldolase activity. Designs that used an explicit water molecule to mediate proton shuffling were significantly more successful, with rate accelerations of up to four orders of magnitude and multiple turnovers, than those involving charged side-chain networks. The atomic accuracy of the design process was confirmed by the x-ray crystal structure of active designs embedded in two protein scaffolds, both of which were nearly superimposable on the design model.

The Diels-Alder reaction is a cornerstone in organic synthesis, forming two carbon-carbon bonds and up to four new stereogenic centers in one step. No naturally occurring enzymes have been shown to catalyze bimolecular Diels-Alder reactions. We describe the de novo computational design and experimental characterization of enzymes catalyzing a bimolecular Diels-Alder reaction with high stereoselectivity and substrate specificity. X-ray crystallography confirms that the structure matches the design for the most active of the enzymes, and binding site substitutions reprogram the substrate specificity. Designed stereoselective catalysts for carbon-carbon bond forming reactions should be broadly useful in synthetic chemistry.

We describe a computationally designed enzyme, formolase (FLS), which catalyzes the carboligation of three one-carbon formaldehyde molecules into one three-carbon dihydroxyacetone molecule. The existence of FLS enables the design of a new carbon fixation pathway, the formolase pathway, consisting of a small number of thermodynamically favorable chemical transformations that convert formate into a three-carbon sugar in central metabolism. The formolase pathway is predicted to use carbon more efficiently and with less backward flux than any naturally occurring one-carbon assimilation pathway. When supplemented with enzymes carrying out the other steps in the pathway, FLS converts formate into dihydroxyacetone phosphate and other central metabolites in vitro. These results demonstrate how modern protein engineering and design tools can facilitate the construction of a completely new biosynthetic pathway.computational protein design | pathway engineering | carbon fixation N ovel strategies are needed to address current challenges in energy storage and carbon sequestration. One approach is to engineer biological systems to convert one-carbon compounds into multicarbon molecules such as fuels and other high value chemicals. Many synthetic pathways to produce value-added chemicals from common feedstocks, such as glucose, have been constructed in organisms that lack one-carbon anabolic pathways, such as Escherichia coli or Saccharomyces cerevisiae (1-3); however, despite considerable effort, it has been difficult to introduce heterologous one-carbon fixing pathways into these organisms (4). Likely problems include the inherent complexity, environmental sensitivity, inefficiency, or unfavorable chemical driving force of naturally occurring one-carbon metabolic pathways (5).An optimal pathway for one-carbon utilization in common synthetic biology platforms would be (i) composed of a minimal number of enzymes, (ii) linear and disconnected from other metabolic pathways, (iii) thermodynamically favorable with a significant driving force at most or all steps, and (iv) capable of functioning in a robust manner under both aerobic and anaerobic conditions (5). A pathway with these properties could enable the assimilation of one-carbon molecules as the sole carbon source for the production of fuels and chemicals. Although no such pathway is known in nature, the established electrochemical reduction of carbon dioxide to formate under ambient temperatures and pressures in neutral aqueous solutions provides an attractive starting point for a onecarbon fixation pathway (5-8).We describe the computational design of an enzyme that catalyzes the carboligation of three one-carbon molecules into a single three-carbon molecule. This enzyme enables the construction of a new pathway, the formolase pathway, in which formate is converted into the central metabolite dihydroxyacetone phosphate (DHAP; Fig. 1). The use of computational protein design to reengineer catalytic activities opens up the pathway design space beyond that available based o...

The ability to redesign enzymes to catalyze non-cognate chemical transformations would have wide-ranging applications. We developed a computational method for repurposing the reactivity of active site functional groups of metalloenzymes to catalyze new reactions. Using this method, we engineered a zinc-containing murine adenosine deaminase to catalyze the hydrolysis of a model organophosphate with a catalytic efficiency kcat/Km ~104 M−1s−1 after directed evolution. In the high-resolution crystal structure of the enzyme, all but one of the designed residues adopt the designed conformation. The designed enzyme efficiently catalyzed the hydrolysis of the RP-isomer of a coumarinyl analog of the nerve agent cyclosarin, and showed striking substrate selectivity for coumarinyl leaving-groups. Computational redesign of native enzyme active sites complements directed evolution methods and offers a general approach for exploring their untapped catalytic potential for new reactivities.