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 unusually strong reversible binding of biotin by avidin and streptavidin has been investigated by density functional and MP2 ab initio quantum mechanical methods. The solvation of biotin by water has also been studied through QM/MM/MC calculations. The ureido moiety of biotin in the bound state hydrogen bonds to five residues, three to the carbonyl oxygen and one for each -NH group. These five hydrogen bonds act cooperatively, leading to stabilization that is larger than the sum of individual hydrogen-bonding energies. The charged aspartate is the key residue that provides the driving force for cooperativity in the hydrogen-bonding network for both avidin and streptavidin by greatly polarizing the urea of biotin. If the residue is removed, the network is disrupted, and the attenuation of the energetic contributions from the neighboring residues results in significant reduction of cooperative interactions. Aspartate is directly hydrogen-bonded with biotin in streptavidin and is one residue removed in avidin. The hydrogen-bonding groups in streptavidin are computed to give larger cooperative hydrogen-bonding effects than avidin. However, the net gain in electrostatic binding energy is predicted to favor the avidin-bicyclic urea complex due to the relatively large penalty for desolvation of the streptavidin binding site (specifically expulsion of bound water molecules). QM/MM/MC calculations involving biotin and the ureido moiety in aqueous solution, featuring PDDG/PM3, show that water interactions with the bicyclic urea are much weaker than (strept)avidin interactions due to relatively low polarization of the urea group in water.

Glutamate transporters regulate excitatory amino acid neurotransmission across neuronal and glial cell membranes by coupling the translocation of their substrate (aspartate or glutamate) into the intracellular (IC) medium to the energetically favorable transport of sodium ions or other cations. The first crystallographically resolved structure of this family, the archaeal aspartate transporter, GltPh, has served as a structural paradigm for elucidating the mechanism of substrate translocation by these transporters. Two helical hairpins, HP2 and HP1, at the core domains of the three subunits that form this membrane protein have been proposed to act as the respective extracellular and IC gates for substrate intake and release. Molecular dynamics simulations using the outward-facing structure have confirmed that the HP2 loop acts as an EC gate. The mechanism of substrate release at atomic scale, however, remained unknown due to the lack of structural data until the recent determination of the inward-facing structure of GltPh. In the present study, we use this recently resolved structure to simulate the release of substrate to the cytoplasm and the roles of HP1 and HP2 in this process. The highly flexible HP2 loop is observed to serve as an activator (or initiator) prompting the release of a gatekeeper Na+ to the cytoplasm and promoting the influx of water molecules from the cytoplasm, which effectively disrupt substrate–protein interactions and drive the dislodging of the substrate from its binding site. The completion of substrate release and exit, however, entails the opening of the highly stable HP1 loop as well. Overall, the unique conformational flexibility of the HP2 loop, the dissociation of a Na+, the hydration of binding pocket, and final yielding of the HP1 loop 3-Ser motif emerge as the successive events controlling the release of the bound substrate to the cell interior by glutamate transporters.

The design of active sites has been carried out using quantum mechanical calculations to predict the rate-determining transition state of a desired reaction in presence of the optimal arrangement of catalytic functional groups (theozyme). Eleven versatile reaction targets were chosen, including hydrolysis, dehydration, isomerization, aldol, and Diels-Alder reactions. For each of the targets, the predicted mechanism and the rate-determining transition state (TS) of the uncatalyzed reaction in water is presented. For the rate-determining TS, a catalytic site was designed using naturalistic catalytic units followed by an estimation of the rate acceleration provided by a reoptimization of the catalytic site. Finally, the geometries of the sites were compared to the X-ray structures of related natural enzymes. Recent advances in computational algorithms and power, coupled with successes in computational protein design, have provided a powerful context for undertaking such an endeavor. We propose that theozymes are excellent candidates to serve as the active site models for design processes.