Sandbox Reserved 1852

Introduction

The Diels-Alderase aims to catalyze the Diels-Alder reaction [Fig. 1A] between 4-carboxybenzyl-trans-1,3-butadiene-1-carbamate and N,N-dimethylacrylamide [Fig. 1B] for use in synthetic organic chemistry. Specifically, the enzyme surpasses uncatalyzed reactions by generating a product that is entirely stereoselective for the 3R,4S endo form [Fig.1C]. The Diels-Alderase was built using de novo enzyme design, using computational modeling and refinement through collaborative problem-solving from online users. The first generation Diels-Alderase was made using the Rosetta computational design program, where a potential active site was built and tested against a library of scaffold proteins. Later, as the active site was perfected, future generations of the Diels-Alderase were made using an online protein folding game called Foldit, where players competed to improve binding efficiency by completing various challenges.^[1]

Figure 1. A) Example mechanism of a simple Diels-Alder reaction. B) Diels-Alderase substrates. Diene is 4-carboxybenzyl trans-1,3-butadiene-1-carbamate; dienophile is N,N- dimethylacrylamide. C) Illustration of 3R, 4S endo stereoisomerism, which the Diels-Alderase is selective for.

The Diels-Alderase was designed to connect a diene and dienophile in a Diels-Alder reaction. It accomplishes this by decreasing the energy gap between the dienophile’s lowest unoccupied molecular orbital (LUMO) and the diene’s highest occupied molecular orbital (HOMO) in the transition state.^[2] The current most active form of the Diels-Alderase is modelled under the PDB code 4o5t.

The binding pocket of 4o5t is selective for two substrates, 4-carboxybenzyl trans-1,3-butadiene-1-carbamate (diene) and N,N- dimethylacrylamide (dienophile). These substrates are shown as a single, combined ligand, 4-{[2-(phosphonooxy)ethyl]carbamoyl}benzyl [(1R,6S)-6-(dimethylcarbamoyl)cyclohex-2-en-1-yl]carbamate, in the protein model. The binding site contains a hydrogen bond donor (Tyr134) which lowers the LUMO energy and stabilizes the negative charge on the dienophile and a hydrogen bond acceptor (Glu208) that increases the HOMO energy and stabilizes the positive charge on the diene.^[2] Both of these H-bonding interactions work to stabilize the transition state, while also orienting the substrates in optimal conformations for reacting.

Overall, the Diels-Alderase stimulates improvement in synthetic laboratories and demonstrates early success in the now-prominent world of computational enzyme design.

General Structure

Figure 2. Binding pocket and substrate. Shown is the binding pocket of the enzyme shown as surface, highlighting the electrostatics of the two catalytic residues, Tyr134 and Glu208. The ligand is color coded based on original structure: the dieophile is in yellow and the diene is in green. The reaction proceeds via attack of the C6 on the C5, shifting electron density to C2, which attacks C1.

Scaffold

After early Rosetta computational modelling, an ideal protein was found in the 6-bladed beta-propeller of Loligo vulgalis, or the European Squid. ^[2][3] The protein is relatively simple, with only one chain, one unit, 324 residues, and no extra ligands, metal ions, or small molecules bound.

Active Site

In the designed active site, stabilize the transition state of the Diels-Alder reaction. The Tyr134 acts as a to the oxygen on the dienophile (see Fig. 2). Q208 acts as a to the nitrogen on the diene. These interactions help reduce the energetic gap between orbitals allowing the reaction to proceed.

Helix Cap

In the evolution process, a 16-residue alpha-helix to the top of the binding site. The hydrophobic helix “functions as a lid to constrain the substrates in a productive orientation for reaction,” decreasing the K_m of the enzyme and increasing the catalytic efficiency, as seen in the measured kinetics of the enzyme.^[1]

Mechanism

Figure 3. Active site mechanism

The key to this enzyme’s success as a catalyst lies in its ability to lower the energy gap between reactants. To accomplish this, the two active site residues, Tyr134 and Glu208, use hydrogen bonding to assist the reaction in a variety of ways. First, it allows specific binding of the ligand in the active site, selecting for molecules with certain stereochemistry at and around the catalytic residues, specifically the carbamate and carbonyl of the diene and dienophile, respectively. This promotes the reaction by stabilizing the molecules in close proximity to one another, which also promotes the characteristic stereoselectivity.^[2] Second, the bonds affect the energetics of the molecules. By donating a hydrogen to the carbonyl of dienophile, Tyr134 helps to decrease the electron density around the molecule, lowering the energy of the lowest unoccupied molecular orbital (LUMO). Conversely, by abstracting the hydrogen from the carbamate of the dienophile, Glu208 increases the electron density and thus the energy of the highest occupied molecular orbital (HOMO). By closing the gap between these orbitals, the enzyme lowers the activation energy required to allow the orbitals to react. Finally, these interactions help to stabilize the accumulated charges in the transition state. By decreasing electron density in the dienophile, Tyr134 helps to stabilize the accumulated negative charge in the transition state. The Glu208, then, helps stabilize the accumulated positive charge by increasing the electron density of the diene. Calculations predict that this helps to stabilize the transition state by nearly 5 kcal/mol.^[2] All together, these interactions make it much easier for the reaction to proceed in a very stereoselective and favorable manner.

Development and Evolution

DA_20_00

During initial computer modelling, over one million potential Diels-Alderase active sites were matched to potential protein scaffolds.^[2] Computer optimization narrowed this down to 84 potential models on various scaffolds, and researchers attempted to grow and purify those proteins within an E. coli host. Of the 50 proteins that were successfully purified, only 2 proteins proved to be sufficiently active after LC-MS screening. DA_20_00, which used a beta-propeller scaffold, had the most success in further mutations and therefore became the Diels-Alderase of choice.^[2] However, this initial enzyme's active site had very little catalytic activity, seen in its low catalytic efficiency after kinetic screening.^[2][4]

DA_20_10

DA_20_10 provided key mutations in and around the active site that increased the hydrophobicity, provided structural stability, and increased interactions between the ligand and surrounding residues.

Q162R

Glu 162, a , resides near the top of the binding site, and is about than 3Å from the ligand in most models on the enzyme. It can act as a hydrogen bond donor to the terminal phosphate on the ligand when in proximity. To increase this interaction, Glu162 was mutated to an , which decreased the length of the potential hydrogen bond to within 2.5Å in most models, increasing the strength of the interaction.^[2]

S284A

Ser 284 resides deep within the binding pocket of the enzyme. Choosing a to mutation increases the hydrophobicity of the binding pocket and reduce reactivity, without also changing any steric characteristics in the region unintentionally near the catalytic residues.^[2]

A285N

, as follows, is also buried within the binding pocket. Introducing this mutation increases steric hindrance with the catalytic tyrosine, reducing the number of rotamers the residue has to increase the reactivity of the enzyme by lowering the distance between Y134 and the ligand.^[2]

CE6

The DA_20_10 model of the Diels Alderase was further enhanced by players of the online game "Foldit." Building on preliminary early data, players were asked to optimize various helical structures that would surround and support the ligand. After over 100,000 designs were tested, the top-scoring CE6 model was finalized, containing as that favorably constrains ligand orientation. This "cap" consists of two helices--helix one spans from residues 36-44, and helix two spans from residues 48-56.^[1]

CE20

In this generation, it was found that the most catalytically efficient models had mutated T34, P48, and R56 to . These mutations further tightened the binding pocket and create a more hydrophobic environment.^[4]

Kinetics

Figure 4. A) Catalytic efficiencies of key Diels-Alderase generations. Kinetic data was measured at 25°C, in PBS, at pH 7.4. B) Improvement of catalytic efficiency across generations.^[4]

Classic Michaelis-Menten kineticswere determined for each generation of the enzyme. As the Diels-Alderase relies on a catalyzed interaction between both the diene and dienophile, a Michaelis binding constant (K_m value) was determined for each substrate separately before catalytic efficiency was calculated. The CE20 model of the enzyme is over 300-fold more efficient than the first enzyme model due to increasing active site specificity.^[4]

Applications

The CE20 model is the most efficient Diels-Alderase yet, surpassing many other biological (antibody) and artificial (ribozyme, metalloenzyme) attempts at catalyzing the Diels-Alder reaction. Even then, the CE20 model has a catalytic efficiency value at least 4 orders of magnitude lower than those seen in other moderately-efficient natural enzymes, demonstrating the innate slowness of the Diels-Alder reaction.^[4]

Though the rate of product formation using this enzyme is not significantly different from that found when reactants reflux free in solution, the Diels-Alderase shows a vast improvement in product stereoselectivity. When refluxed in a room temperature aqueous solution containing the necessary substrates, the enzyme catalyzed an over 90% conversion rate, producing only the 3R,4S endo cyclohexane product isomer. By comparison, refluxing the substrates free in solution for a similar duration of time yields a racemic (66:34) mixture of endo and exo products. It is primarily for these stereoselective benefits that this enzyme is valuable for synthetic purposes. ^[4]

Future improvement of the Diels-Alderase will likely revolve around the improvement of catalytic efficiency, further constriction of the active site, and selective production of varying stereoisomers.

References

1. ↑ ^{1.0 1.1 1.2} Eiben CB, Siegel JB, Bale JB, Cooper S, Khatib F, Shen BW, Players F, Stoddard BL, Popovic Z, Baker D. Increased Diels-Alderase activity through backbone remodeling guided by Foldit players. Nat Biotechnol. 2012 Jan 22;30(2):190-2. doi: 10.1038/nbt.2109. PMID:22267011 doi:http://dx.doi.org/10.1038/nbt.2109
2. ↑ ^{2.00 2.01 2.02 2.03 2.04 2.05 2.06 2.07 2.08 2.09 2.10} Siegel JB, Zanghellini A, Lovick HM, Kiss G, Lambert AR, St Clair JL, Gallaher JL, Hilvert D, Gelb MH, Stoddard BL, Houk KN, Michael FE, Baker D. Computational design of an enzyme catalyst for a stereoselective bimolecular Diels-Alder reaction. Science. 2010 Jul 16;329(5989):309-13. PMID:20647463 doi:329/5989/309
3. ↑ Scharff EI, Koepke J, Fritzsch G, Lucke C, Ruterjans H. Crystal structure of diisopropylfluorophosphatase from Loligo vulgaris. Structure. 2001 Jun;9(6):493-502. PMID:11435114
4. ↑ ^{4.0 4.1 4.2 4.3 4.4 4.5} Preiswerk N, Beck T, Schulz JD, Milovnik P, Mayer C, Siegel JB, Baker D, Hilvert D. Impact of scaffold rigidity on the design and evolution of an artificial Diels-Alderase. Proc Natl Acad Sci U S A. 2014 May 20. pii: 201401073. PMID:24847076 doi:http://dx.doi.org/10.1073/pnas.1401073111

Student Contributors

Taylor Donahue, Kate Thuma, Micah Zile