Robust design and optimization of retroaldol enzymes.

Enzyme catalysts of a retroaldol reaction have been generated by computational design using a motif that combines a lysine in a nonpolar environment with water-mediated stabilization of the carbinolamine hydroxyl and ß-hydroxyl groups. Here, we show that the design process is robust and repeatable, with 33 new active designs constructed on 13 different protein scaffold backbones. The initial activities are not high but are increased through site-directed mutagenesis and laboratory evolution. Mutational data highlight areas for improvement in design. Different designed catalysts give different borohydride-reduced reaction intermediates, suggesting a distribution of properties of the designed enzymes that may be further explored and exploited.

Structural analyses of covalent enzyme-substrate analog complexes reveal strengths and limitations of de novo enzyme design.

We report the cocrystal structures of a computationally designed and experimentally optimized retro-aldol enzyme with covalently bound substrate analogs. The structure with a covalently bound mechanism-based inhibitor is similar to, but not identical with, the design model, with an RMSD of 1.4 Å over active-site residues and equivalent substrate atoms. As in the design model, the binding pocket orients the substrate through hydrophobic interactions with the naphthyl moiety such that the oxygen atoms analogous to the carbinolamine and ß-hydroxyl oxygens are positioned near a network of bound waters. However, there are differences between the design model and the structure: the orientation of the naphthyl group and the conformation of the catalytic lysine are slightly different; the bound water network appears to be more extensive; and the bound substrate analog exhibits more conformational heterogeneity than typical native enzyme-inhibitor complexes. Alanine scanning of the active-site residues shows that both the catalytic lysine and the residues around the binding pocket for the substrate naphthyl group make critical contributions to catalysis. Mutating the set of water-coordinating residues also significantly reduces catalytic activity. The crystal structure of the enzyme with a smaller substrate analog that lacks naphthyl ring shows the catalytic lysine to be more flexible than in the naphthyl-substrate complex; increased preorganization of the active site would likely improve catalysis. The covalently bound complex structures and mutagenesis data highlight the strengths and weaknesses of the de novo enzyme design strategy.

Kemp elimination catalysts by computational enzyme design.

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(-1)s(-1) 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.

De novo computational design of retro-aldol enzymes.

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.

New algorithms and an in silico benchmark for computational enzyme design.

The creation of novel enzymes capable of catalyzing any desired chemical reaction is a grand challenge for computational protein design. Here we describe two new algorithms for enzyme design that employ hashing techniques to allow searching through large numbers of protein scaffolds for optimal catalytic site placement. We also describe an in silico benchmark, based on the recapitulation of the active sites of native enzymes, that allows rapid evaluation and testing of enzyme design methodologies. In the benchmark test, which consists of designing sites for each of 10 different chemical reactions in backbone scaffolds derived from 10 enzymes catalyzing the reactions, the new methods succeed in identifying the native site in the native scaffold and ranking it within the top five designs for six of the 10 reactions. The new methods can be directly applied to the design of new enzymes, and the benchmark provides a powerful in silico test for guiding improvements in computational enzyme design.

Investigation of the mechanism of resistance to third-generation cephalosporins by class C beta-lactamases by using chemical complementation.

The widespread use of antibiotics to treat bacterial infections has led to the continuing challenge of antibiotic resistance. For beta-lactam antibiotics, the most common form of resistance is the expression of beta-lactamase enzymes, which inactivate the antibiotics by cleavage of the beta-lactam core. In this study, chemical complementation, which is a general method to link the formation or cleavage of a chemical bond to the transcription of a reporter gene in vivo, was employed in combination with combinatorial mutagenesis to study the mechanism by which the class C beta-lactamase P99 might evolve resistance to the commonly administered third-generation cephalosporin cefotaxime. The chemical complementation system was first shown to be able to distinguish between the wild-type (wt) class C beta-lactamase P99 and the clinically isolated extended-spectrum class C beta-lactamase GC1 in the presence of cefotaxime. The system was then employed to evaluate the activity of mutants of wt P99 towards cefotaxime. A number of single-point mutations at position 221 (Tyr in wt P99) were identified that conferred resistance towards inhibition by cefotaxime, with as much as a 2000-fold increase in k(cat) and a 100-fold increase in k(cat)/K(M) (k(cat)=the rate of catalysis; K(M)=the Michaelis constant), as compared to those of the wt enzyme. Finally, the chemical complementation system was employed in a high-throughput screen to identify a number of mutants of P99 that have multiple mutations around the substrate-binding pocket that increase resistance towards cefotaxime inhibition. The catalytic turnover of cefotaxime by the most active mutant identified was 5500 times higher than that of the wt P99. The resistant mutants suggest a mechanism by which a number of mutations can confer resistance by increasing the flexibility of the Omega loop and altering the positioning of residue 221. Thus, as illustrated in this study, chemical complementation has the potential to be used as a high-throughput screen to study a wide range of enzyme-drug interactions.

Correlation between ligand-receptor affinity and the transcription readout in a yeast three-hybrid system.

The yeast two-hybrid assay has proven to be a powerful method to detect protein-protein interactions as well as to derive genome-wide protein interaction maps. More recently, three-hybrid assays have emerged as a means to detect both protein-RNA and protein-small molecule interactions. Despite the routine use of the two-hybrid assay and the potential of three-hybrid systems, there has been little quantitative characterization to understand how the strength of the protein interaction correlates with transcription activation. It is not known if the additional interaction in three-hybrid systems compromises the sensitivity of the system. Thus, here, we set out to determine the K(D) cutoff of a small molecule three-hybrid system and to determine if there is a correlation between the K(D) and the levels of transcription activation. A series of mutations to FK506-binding protein 12 (FKBP12) were designed to vary the affinity of this protein for the small molecule synthetic ligand for FK506-binding protein 12 (SLF). These FKBP12 variants were overexpressed and purified, and their K(D)'s for SLF were measured using a fluorescence polarization assay. Then the levels of transcription activation in a Mtx-DHFR yeast three-hybrid system were determined for these variants using a lacZ reporter gene. The K(D) cutoff of the Mtx yeast three-hybrid system is found to be ca. 50 nM. Further, the levels of transcription activation correlate with the strength of the binding interaction, though the dynamic range is only 1 order of magnitude. These results establish that the three-hybrid assay has the requisite sensitivity for drug discovery. However, the small dynamic range highlights a limitation to equilibrium-based assays for discriminating interactions based on affinity.

Receptor-dependence of the transcription read-out in a small-molecule three-hybrid system.

Small-molecule three-hybrid systems show promise as an in vivo alternative to affinity chromatography for detecting small-molecule-protein interactions. While several three-hybrid systems have been reported, little has been done to characterize these systems and, in particular, to test the assumption that the protein-small-molecule interaction can be varied without disrupting the transcription read-out. Recently we reported a dexamethasone-methotrexate chemical inducer of dimerization (CID) for use in the yeast three-hybrid system, based on the well-studied ligand-receptor pairs dexamethasone (Dex)-glucocorticoid receptor (GR) and methotrexate (Mtx)-dihydrofolate reductase (DHFR). Here we describe our first efforts to characterize this system, by focusing on a comparison of the activity of a bacterial and a mammalian DHFR as a test case of the influence of the ligand-receptor pair on the transcription read-out. By using a lacZ reporter gene, the activity of several GR and DHFR protein chimeras with different orientations and linker sequences and Dex-Mtx CIDs with different chemical linkers have been compared. In addition, Western analyses and in vivo biochemical assays have been carried out to confirm the integrity of the GR and DHFR protein chimeras. The transcription read-out is found to be much more sensitive to the structure of the protein chimeras than the CID. The most surprising result is that the levels of transcription activation are consistently higher with the bacterial than the mammalian DHFR, despite the fact that both proteins bind Mtx with an inhibition constant (K(I)) in the low pM range. These results set the stage for understanding three-hybrid systems at the biochemical level so that they can be used to detect ligand-receptor pairs with a range of structures and dissociation constants.