How an enzyme might accelerate an intramolecular diels-alder reaction: Theozymes for the formation of salvileucalin B

Organic Letters

Volumn 12, Issue 6, 2010, Pages 1164-1167

  Tantillo, Dean J ^a  

^a UNIVERSITY OF CALIFORNIA   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

BIOLOGICAL PRODUCT; DITERPENE; ENZYME; SALVILEUCALIN B;

ARTICLE; BIOCATALYSIS; BIOSYNTHESIS; CHEMICAL STRUCTURE; CHEMISTRY; COMPUTER SIMULATION; ENZYME ACTIVE SITE; METABOLISM; QUANTUM THEORY; STEREOISOMERISM;

BIOCATALYSIS; BIOLOGICAL PRODUCTS; CATALYTIC DOMAIN; COMPUTER SIMULATION; DITERPENES; ENZYMES; MOLECULAR STRUCTURE; QUANTUM THEORY; STEREOISOMERISM;

EID: 77949867534     PISSN: 15237060     EISSN: 15237052     Source Type: Journal    
DOI: 10.1021/ol9028435     Document Type: Article

References (53)

1. (a) Aoyagi, Y.; Yamazaki, A.; Nakatsugawa, C.; Fukaya, H.; Takeya, K.; Kawaucji, S.; Izumi, H. Org. Lett. 2008, 10, 4429-4432.
2. (b) For a recent synthetic approach, to the salvileucalins, see: Levin, S.; Nani, R. R.; Reisman, S. E. Org. Lett. 2010, 12, 780-783.
3. Reviews: (a) Kelly, W. L. Org. Biomol. Chem. 2008, 6, 4483-4493.
4. (b) Stocking, E. M.; Williams, R. M. Aneew. Chem.. Int. Ed. 2003, 42, 3078-3115.
5. (C) Oikawa, H. Bull Chem. Soc. Jpn. 2005, 78, 537-554.
6. (d) Oikawa, H.; Tokiwano, T. Nat. Prod. Rep. 2004, 21, 321-352.
7. (e) Pohnert, G. ChemBioChem. 2001, 2, 873-875.
8. Macrophomate synthase was touted by some as the first wellcharacterized Diels-Alderase enzyme, but evidence has been accumulated that indicates that the cyclization catalyzed by this enzyme is stepwise, and therefore not truly pericyclic. See, for example: (a) Guimarães, C. R. W.; Udier-Blagovió, M.; Jorgensen, W. L. J. Am. Chem. Soc. 2005, 127, 3577-3588.
9. (b) Ose, T.; Watanabe, K.; Mie, T.; Honma, M.; Watanabe, H.; Yao, M.; Oikawa, H.; Tanaka, I. Nature 2003, 422, 185-189.
10. Keinan, W., Ed. Catalytic Antibodies; Wiley-VCH: Weinheim, Germany, 2005.
11. GAUSSIAN03 (Frisch, M. J. et al., GAUSSIAN03, Revision D.Ol; Gaussian, Inc., Pittsburgh, PA, 2003,
12. full reference in the Supporting Information) was employed for all calculations. All geometries were optimized at the B3LYP/6-31G(d) level of theory: Becke, A. D. J. Chem. Phys. 1993, 98, 5648-5652.
13. Becke, A. D. J. Chem. Phys. 1993, 98, 1372-1377.
14. Lee, C.; Yang, W.; Parr, R. G. Phys. Rev. B. Solid State 1988, 37, 785-789.
15. Stephens, P. J.; Devlin, F. J.; Chabalowski, C. F.; Frisch, M. J. J. Phys. Chem. 1994, 98, 11623-11627.
16. All stationary points were characterized as minima or transition state structures by analyzing their vibrational frequencies. Intrinsic reaction coordinate (IRC) calculations (Gonzalez, C.; Schlegel, H. B. J. Phys. Chem. 1990, 94, 5523-5527.
17. Fukui, K. Acc. Chem. Res. 1981, 14, 363-368)
18. 6 and MP2/6-31G(d) (Møller, C.; Plesset, M. S. Phys. Rev. 1934, 46, 618-622) were performed on some optimized B3LYP/6-31.G(d) geometries, as described in the text.
19. AU reported energies from B3LYP/6-31G(d) calculations include zero point energy corrections from frequency calculations, scaled by 0.9806: Scott, A. P.; Radom, L. J. Phys. Chem. 1996, 100, 16502-16513.
20. For some systems, the effects of solvent (CHCI3 [ε. 4.9]) were modeled using CPCM calculations (with UAKS radii), a self-consistent reaction field (SCRF) method: Barone, V.; Cossi, M. J. J. Phys. Chem. A 1998, 102, 1995-2001.
21. Barone, V.; Cossi, M.; Tomasi, J. J. Comput. Chem. 1998, 19, 404-417.
22. Takano, Y.; Houk, K. N. J. Chem. Theor. Comput. 2005, 1, 70-77.
23. estert⋯ H distance was 0.1.2 Å longer in solvent than in the gas phase; changes to the bonds that form and break during the Diels-Alder reaction, as well as changes to the O-H bonds in the water molecule, were all well less than 0.01 A. The barrier computed with structures optimized in solvent was 27.0 kcal/mol, ∼0.4 kcal/mol lower than that computed with single-point energies (see Figure 2). For system K, full optimization of the transition state structure in solvent failed, since, after an initial period of structural optimization, the calculation consistently oscillated between structures extremely close in geometry and differing in energy by less than 0.02 kcal/mol. Using the structure with lower energy of these two and the fully optimized reactant, an upper limit for the barrier of 25.8 kcal/mol (not zero-point energy corrected) was computed; this value is only ∼0.2 kcal/mol lower than that calculated with single-point energies (see Figure 2). Here, the largest geometric differences between the transition state structures optimized in the gas phase and in solvent were found in the O ⋯ H-N substructure: the O⋯ H distance was 0.16 Ä longer and the H-N distance was 0.06 Å shorter in solvent than in the gas phase; changes to the bonds that form and break during the Diels-Alder reaction were all less than 0.025 Å. Thus, using solvent single-point energies to compute barriers seems reasonable. The structural drawing in Figure 3 was produced using Ball & Stick (Müller, N.; Falk, A. Ball & Stick V.3.7.6, molecular graphics application for MacOS computers; Johannes Kepler University: Linz, Austria, 2000).
24. Matsuda, S. P. T.; Wilson, W. K.; Xiong, Q. Org. Biomol. Chem. 2006, 4, 530-543.
25. The furan ring was replaced by a hydrogen atom in order to reduce the time required for computations (this ring system would add five nonhydrogen atoms and some conformational flexibility). We feel that this approach is reasonable since the furan ring is not directly conjugated to any of the π-systems present in the salvileucalins.
26. (a) Ess, D. H.; Houk, K. N. J. Phys. Chem. A 2005, 109, 9542-9553.
27. (b) Guner, V. A.; Khuong, K. S.; Houk, K. N.; Chuma, A.; Pulay, P. J. Phys. Chem. A 2004, 108, 2959-2965.
28. (c) Guner, V.; Khuong, K. S.; Leach, A. G.; Lee, P. S.; Bartberger, M. D.; Houk, K. N. J. Phys. Chem. A 2003, 107, 11445-11459.
29. (d) Wiest, O.; Montiel, D. C.; Houk, K. N. J. Phys. Chem. A 1997, 101, 8378-8388.
30. For discussion and leading references, see: (a) Lee, J. K.; Tantillo, D. J. Adv. Phys. Org. Chem. 2003, 38, 183-218.
31. (b) Shutz, C. N.; Warshel, A. Proteins Struct. Funct. Genjet. 2001, 44, 400-417.
32. (c) Simonson, T.; Brooks, C. L., III. J. Am. Chem. Soc. 1996, 118, 8452-8458.
33. Note that all of the reactions in Figure 1 are predicted to be exothermic; trading π-bonds for σ-bonds counteracts the increase in strain upon cycloaddition.
34. For leading references on related examples, see: (a) Tantillo, D. J. Angew. Chem., Int. Ed. 2009, 48, 31-32. Selected studies involving [4 + 2] reactions:
35. (b) Linder, M.; Brinck, T. Org. Biomol. Chem. 2009, 7, 1304-1311.
36. (c) Anderson, C. D.; Dudding, T.; Gordillo, R.; Houk, K. N. Org. Lett. 2008, 10, 2749-2752.
37. (d) Zhang, X.; DeChancie, J.; Gunayadin, H.; Howdry, A. B.; Clemente, F. R.; Smith, A. J. T.; Handel, T. M.; Houk, K. N. J. Org. Chem. 2008, 73, 889-899.
38. (e) Cannizzaro, C. E.; Ashley, J. A.; Janda, K. D.; Houk, K. N. J. Am. Chem. Soc. 2003, 125, 2489-2506.
39. (f) Zhang, X.; Deng, Q.; Yoo, S. H.; Houk, K. N. J. Org. Chem. 2002, 67, 9043-9053.
40. For leading references, see: (a) Warshel, A.; Sharma, P. K.; Kato, M.; Xiang, Y.; Liu, H.; Olsson, M. H. M. Chem. ReV. 2006, 106, 3210-3235. Table 1 of this review lists enzymatic reactions for which free energies of activation (bound reactant(s) to bound rate-determining transition state structure) have been determined, and these range from 11 to 19 kcal/mol.
41. (b) Garcia-Viloca, M.; Gao, J.; Karplus, M.; Truhlar, D. G. Science 2004, 303, 186-195.
42. (c) Snider, M. G.; Temple, B. S.; Wolfenden, R. J. Phys. Org. Chem. 2004, 17, 586-591.
43. (d) Wolfenden, R.; Snider, M. J. Acc. Chem. Res. 2001, 34, 938-945.
44. Reviews: (a) Tantillo, D. J.; Chen, J.; Houk, K. N. Curr. Op. Chem. Biol. 1998, 2, 743-750.
45. (b) Tantillo, D. J.; Houk, K. N. Theozymes and Catalyst Design. In Stimulating Concepts in Chemistry; Wiley-VCH: Weinheim, Germany, 2000; pp 79-88.
46. (c) Sousa, F. P.; Fernandes, P. A.; Ramos, M. J. J. Phys. Chem. A 2009, 113, 14231-14236
47. For leading references, see: (a) Ordentlich, A.; Barak, D.; Kronman, C.; Ariel, N.; Segall, Y.; Velan, B.; Shafferman, A. J. Biol. Chem. 1998, 273, 19509-19517.
48. (b) Whiting, A. K.; Peticolas, W. L. Biochemistry 1994, 33, 552-561.
49. (a) Tantillo, D. J.; Houk, K. N. J. Comput. Chem. 2002, 23, 84-95.
50. (b) Tantillo, D. J.; Houk, K. N. Chem. Bial. 2001, 8, 535-545.
51. (c) Zhou, G. W.; Guo, J.; Huang, W.; Scanlan, T. S.; Fletterick, R. J. Science 1994, 265, 1059-1064.
52. Pauling, L. Nature 1948, 161, 707.
53. 11 all energies (B3LYP/ 6-31G(d)+0.9806ZPE) are in kcal/mol): "Figure Presented"