Next Article in Journal
Design and Synthesis of Thiazolo[5,4-f]quinazolines as DYRK1A Inhibitors, Part I
Next Article in Special Issue
Intramolecular Azide to Alkene Cycloadditions for the Construction of Pyrrolobenzodiazepines and Azetidino-Benzodiazepines
Previous Article in Journal
Five New Iridoids from Roots of Salvia digitaloides
Previous Article in Special Issue
Enantiomerically Pure Phosphonated Carbocyclic 2'-Oxa-3'-Azanucleosides: Synthesis and Biological Evaluation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unusual Structure-Energy Correlations in Intramolecular Diels–Alder Reaction Transition States

by
Justyna M. Żurek
,
Robert L. Rae
,
Martin J. Paterson
* and
Magnus W. P. Bebbington
*
Institute of Chemical Science, School of Engineering and Physical Sciences, William Perkin building, Heriot-Watt University, Edinburgh EH14 4AS, UK
*
Authors to whom correspondence should be addressed.
Molecules 2014, 19(10), 15535-15545; https://doi.org/10.3390/molecules191015535
Submission received: 23 July 2014 / Revised: 1 September 2014 / Accepted: 16 September 2014 / Published: 29 September 2014
(This article belongs to the Special Issue Cycloaddition Chemistry)

Abstract

:
Detailed analysis of calculated data from an experimental/computational study of intramolecular furan Diels–Alder reactions has led to the unusual discovery that the mean contraction of the newly forming C-C σ-bonds from the transition state to the product shows a linear correlation with both reaction Gibbs free energies and reverse energy barriers. There is evidence for a similar correlation in other intramolecular Diels–Alder reactions involving non-aromatic dienes. No such correlation is found for intermolecular Diels–Alder reactions.

Graphical Abstract

1. Introduction

The Diels–Alder reaction of furan has been widely used in the synthesis of complex targets and as a probe for the investigation of substituent effects [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]. We recently reported a comprehensive, computational and experimental study of halogenation effects in intramolecular Diels–Alder reactions [19]. The relative rates of many Diels–Alder reactions can be understood in terms of frontier orbital energy differences between the reaction partners [20,21,22,23,24,25]. Although furan has been widely used as a diene in both inter- and intra-molecular Diels–Alder (IMDAF) processes, and some reaction rates conform to the predictions of FMO theory, substituent effects are less clear cut than for conventional non-aromatic dienes. Other electronic effects, such as positive charge stabilisation, have been shown to be more important than FMO energies in some of these reactions [26,27,28]. Conformational preferences in the tether also play a major role and were the subject of considerable investigation in the 1980s and 1990s [2,29,30].
Our recent work [19] focused on halogenation effects in both diene and dienophile and resulted in the discovery of a C-X dipole-dipole interaction that contributed to the transition state energy. Calculated FMO energies did not correlate well with the experimentally observed reaction rates. Our tentative explanation for this was that the reactions under study had later transition states than are usual for Diels–Alder reactions, such that assuming a reactant-like transition state and using starting material FMO energies was a poor model. Evidence for the character of the transition state was obtained from calculated data, in which those reactions that were predicted to be more exergonic were found to have a greater contraction in the C-C distances in moving along the reaction coordinate from the transition state to the product (and therefore, an earlier transition state), consistent with the Hammond postulate (Scheme 1).
Scheme 1. Summary of computational results from the initial set of IMDAF substrates, 126.
Scheme 1. Summary of computational results from the initial set of IMDAF substrates, 126.
Molecules 19 15535 g005
We define the parameter, hereafter referred to as transition state σ-contraction, as follows:
%age TS σ-contraction = 100 * (Δa/a + Δb/b)/2;
where:
Δa = a − ap and Δb = b − bp
This set of reactions are themselves calculated to be much less exergonic than typical Diels–Alder reactions with non-aromatic dienes [28].

2. Results and Discussion

While these results were consistent with our original hypothesis that we were studying examples with unusually late transition states, we sought to clarify the matter further and now report the additional analysis of our original data and further high-level calculations that support our initial supposition, as well as uncovering an unusual and, so far as we are aware, unprecedented structure-energy correlation between reaction free energy and the contraction of the C-C distances in going from the transition state to the product.
A conventional method for determining the lateness of a transition state from thermochemical data is to construct a plot of activation Gibbs free energy vs. Gibbs free energy of the reaction based on the Leffler Equation (1) [31]:
G = α G°R
A plot of ΔG(activation) vs. ΔG(reaction) for these systems gives a straight line with gradient α = 0.55, indicating a transition state that is marginally product-like (Figure 1). This is in contrast to related calculations on reactions with conventional dienes, which suggest an enthalpic contribution of about 14% to the activation free energy [32], as well as data from KIE experiments on intermolecular reactions of butadiene and acrolein that also indicate early transition states [33].
Figure 1. Leffler plot of ∆G vs.G°R (at 383 K) for our original 26 IMDAF substrates, 126.
Figure 1. Leffler plot of ∆G vs.G°R (at 383 K) for our original 26 IMDAF substrates, 126.
Molecules 19 15535 g001
Examination of the original scatter plots [19] (reproduced in the Supplementary File) reveals a grouping of both activation energies and Gibbs free energies of the reaction according to the type of substrate present, with those containing halogenated dienophiles tending to have higher barriers and less exergonic free energies than those with halogenated furans. Substrates halogenated on both reaction partners were intermediate in both values, and similar trends were observed experimentally. The scatter plot of the average percentage C-C distance contraction from the transition state to the product was strikingly similar to the energy plots, prompting us to probe the nature of a possible correlation between the two. Somewhat to our surprise, a plot of reaction free energy versus average percentage C-C bond contraction showed a strong linear correlation (R2 = 0.93) (Figure 2a–c).
Perhaps less surprisingly, a linear correlation (R2 = 0.96) was also observed between the C-C distance contraction and the free energy change from the transition state to the product (i.e., the reverse Gibbs free energy barrier). Correlation with the activation energy was considerably weaker (R2 = 0.72). Data points include systems with diverse electronic properties and reaction asynchronicities, so the correlation occurs for a wide range of IMDAF substrates. The gradient of the best-fit line represents the sensitivity of the reaction to changes in substituents.
We next elected to probe the generality of the effect, with related calculations on representative sets of Diels–Alder reactions (see the Supplementary File for the details of the substrates). We chose to examine inter- and intra-molecular reactions of furan, cyclopentadiene and butadiene, including systems related to our original halogenated substrates and also those that exhibited widely varying dienophile electronics with both electron-withdrawing and electron-donating groups. This allowed us to examine whether or not the effects we were observing were confined to a specific set of substrates or to the use of aromatic or cyclic dienes. A combined plot of Gibbs’ free energy of the reaction against transition state contraction for all of the newly modelled intramolecular processes again showed a linear relationship, with a strong R2 value of 0.94. This correlation includes substrates with both three- and four-atom tethers, and the correlation appears relatively unaffected by tether length, although this is an issue that would merit further study. Two more plots corresponding to those constructed above showed a similar pattern, with a strong linear correlation between transition state contraction and the activation energy of the reverse reaction (R2 = 0.92), but a much weaker one (R2 = 0.16) between the contraction and the activation energy of the forward process (Figure 3a–c).
Figure 2. Plots of TS σ-contraction vs. key energetic parameters in our experimentally studied data. (a) TS σ-contraction vs. Gibbs free energies of the reaction for IMDAF reactions 126; (b) TS σ-contraction vs. reverse Gibbs free energy barrier for IMDAF reactions 126; (c) TS σ-contraction vs. Gibbs free energy of activation for IMDAF reactions 126.
Figure 2. Plots of TS σ-contraction vs. key energetic parameters in our experimentally studied data. (a) TS σ-contraction vs. Gibbs free energies of the reaction for IMDAF reactions 126; (b) TS σ-contraction vs. reverse Gibbs free energy barrier for IMDAF reactions 126; (c) TS σ-contraction vs. Gibbs free energy of activation for IMDAF reactions 126.
Molecules 19 15535 g002
Figure 3. Plots of TS σ-contraction vs. key energetic parameters in computed intramolecular DA reactions 3739, 5254 and 6568. (a) TS σ-contraction vs. Gibbs free energies of the reaction for new intramolecular DA reactions 3739, 5254 and 6568; (b) TS σ-contraction vs. reverse Gibbs free energy barrier for new intramolecular DA reactions 3739, 5254 and 6568; (c) TS σ-contraction vs. Gibbs free energy of activation for new intramolecular DA reactions 3739, 5254 and 6568.
Figure 3. Plots of TS σ-contraction vs. key energetic parameters in computed intramolecular DA reactions 3739, 5254 and 6568. (a) TS σ-contraction vs. Gibbs free energies of the reaction for new intramolecular DA reactions 3739, 5254 and 6568; (b) TS σ-contraction vs. reverse Gibbs free energy barrier for new intramolecular DA reactions 3739, 5254 and 6568; (c) TS σ-contraction vs. Gibbs free energy of activation for new intramolecular DA reactions 3739, 5254 and 6568.
Molecules 19 15535 g003
Notably, a Leffler plot using all of the new intramolecular Diels–Alder substrates does not show an obvious straight line. This presumably reflects the difference in the electronic properties of furan, cyclopentadiene and butadiene, which appears to limit our capacity to examine their transition state properties in this way. Nevertheless, a significant structure-energy correlation remains.
Examination of similar graphs for the intermolecular processes reveals that the patterns observed for the intramolecular reactions were not sustained. Correlations between the key energetic parameters and the transition state contraction were much weaker, for example, for a selection of intermolecular furan DA reactions (Figure 4a–c), although some evidence for the correlation remains for the reverse energy barrier.
The observed correlations were much weaker still for a sample set of intermolecular cyclopentadiene reactions (R2 = 0.04–0.6 for reactions 4251 and 5556) and butadiene reactions (R2 = 0.00–0.13 for reactions 5764), with a range of different dienophiles (see the Supplementary File). The difference in behaviour between the intermolecular and intramolecular reactions is perhaps related to the order imposed by the tethering of the reaction partners in the intramolecular processes. Indeed, the intermolecular processes we have examined appear to have more scope for asynchronicity, which could conceivably lead to a breakdown in the correlation between the energetics and the average value for the contraction of the two new σ-bond distances.
Figure 4. Plots of TS σ-contraction vs. key energetic parameters in intermolecular furan Diels–Alder reactions 2736, 4041. (a) TS σ-contraction vs. reverse Gibbs free energy of the reaction for intermolecular furan DA reactions 2736, 4041; (b) TS σ-contraction vs. reverse Gibbs free energy barrier for intermolecular furan DA reactions 2736, 4041; (c) TS σ-contraction vs. Gibbs free energy of activation for intermolecular furan Diels–Alder reactions 2736, 4041.
Figure 4. Plots of TS σ-contraction vs. key energetic parameters in intermolecular furan Diels–Alder reactions 2736, 4041. (a) TS σ-contraction vs. reverse Gibbs free energy of the reaction for intermolecular furan DA reactions 2736, 4041; (b) TS σ-contraction vs. reverse Gibbs free energy barrier for intermolecular furan DA reactions 2736, 4041; (c) TS σ-contraction vs. Gibbs free energy of activation for intermolecular furan Diels–Alder reactions 2736, 4041.
Molecules 19 15535 g004
It is apparent from the examination of the calculated transition state data that the correlations we have observed appear to hold for a wide range of different intramolecular Diels–Alder reactions, with very different reaction partners. It is unclear at present whether a simple mathematical expression exists that explains the observed correlations in terms of the underlying intermolecular forces, but in the opinion of the authors, it is likely that the correlation is represented by an approximation to linearity that covers a relatively small range of values for transition state bond contractions. However, the existence of the correlation for this series and, perhaps, as-yet undiscovered reactions may prove useful. For example, such predictable variation of the transition state structure would potentially be of use for those seeking to select reactions for the design of catalytic asymmetric processes, in which the perturbation of the transition state energy and, therefore, structure must be optimised to deliver high enantioselectivity.

3. Experimental Section

All calculations were performed with the Gaussian09 program (versions a, c and d) [34]. Geometries of all reactants, products and transition states were initially optimised with DFT theory, the B3LYP functional and the 6–31G(d) basis set. Transition states were obtained via optimisation using the Berny algorithm. The complete basis set (CBS) method of Petersson and co-workers was then used on these geometries to obtain accurate thermochemical data [35,36,37,38]. The specific extrapolation used here is the CBS-QB3 one, the details of which are given in Paragraph 1 of the Supplementary Information. All calculations were performed in the gas phase. Free energies with the CBS-QB3 model were calculated for two temperatures: 298 K and 383 K (the temperature of the experimental studies) for reactions 126. There was only a minimal influence of the temperature on the optimised bond lengths for substrates and transition states, and so, the calculations for reactions 2768 were run only at a temperature of 298 K.

4. Conclusions

An in-depth analysis of calculated transition state structural data and reaction energetics has revealed an unusual linear correlation between reaction free energy and transition state contraction in intramolecular Diels–Alder reactions. This correlation appears not to be restricted to intermolecular reactions of furan, but does not hold for intermolecular processes. Further investigation is underway to probe the generality of the effect with regard to other types of reaction and to examine its relevance to methodological development.

Supplementary Materials

Supplementary Figure can be accessed at https://www.mdpi.com/1420-3049/19/10/15535/s1.

Supplementary Files

Supplementary File 1

Acknowledgments

We thank the EPSRC (for funding to Justyna M. Żurek and Robert L. Rae) and the EU for an ERC starter grant (to Martin J. Paterson).

Author Contributions

Justyna M. Żurek and Martin J. Paterson carried out the high level calculations. Robert L. Rae, Magnus W. P. Bebbington and Martin J. Paterson analysed the data. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ciganek, E. The Intramolecular Diels–Alder Reaction. In Organic Reactions; John Wiley & Sons: New York, NY, USA, 2004. [Google Scholar]
  2. Kappe, C.O.; Murphree, S.S.; Padwa, A. Synthetic Applications of Furan Diels–Alder chemistry. Tetrahedron 1997, 53, 14179–14233. [Google Scholar]
  3. Vogel, P.; Cossy, J.; Plumet, J.; Arjona, O. Derivatives of 7-oxabicyclo[2.2.1]heptanes in nature and as useful synthetic intermediates. Tetrahedron 1999, 55, 13521–13642. [Google Scholar]
  4. Padwa, A.; Flick, A.C. Intramolecular Diels–Alder reactions for natural product synthesis. Adv. Heterocycl. Chem. 2013, 110, 1–41. [Google Scholar]
  5. Cooper, H.D.; Wright, D.L. The tandem ring opening/ring closing metathesis route to oxaspirocycles. Molecules 2013, 18, 2438–2448. [Google Scholar]
  6. Pitsinos, E.N.; Athinaios, N.; Veroniki, V.P. Enantioselective synthesis of (−)-Laurenditerpenol. Org. Lett. 2012, 14, 4666–4669. [Google Scholar]
  7. Wang, C.-C.; Li, W.-D.Z. A Convergent and Stereocontrolled Cycloaddition Strategy toward Eudesmane Sesquiterpenoid: Total Synthesis of (±)-6β,14-Epoxyeudesm-4(15)-en-1β-ol. J. Org. Chem. 2012, 77, 4217–4225. [Google Scholar]
  8. Fischer, P.; Gruner, M.; Jaeger, A.; Kataeva, O.; Metz, P. Total Synthesis of Pamamycin-6498B. Chem. Eur. J. 2011, 17, 13334–13340. [Google Scholar]
  9. Tanino, K.; Takahashi, M.; Tomata, Y.; Tokura, H.; Uehara, T.; Takashi, N.; Miyashita, M. Total synthesis of solanoeclipon A. Nat. Chem. 2011, 3, 484–488. [Google Scholar]
  10. Petronijevic, F.R.; Wipf, P. Total Synthesis of (+/−)-Cycloclavine and (+/−)-5-epi-Cycloclavine. J. Am. Chem. Soc. 2011, 133, 7704–7707. [Google Scholar]
  11. Xue, Y.-P.; Li, W.-D.Z. Formal Total Synthesis of (+/−)-Estrone via the Furano Diene Approach. J. Org. Chem. 2011, 76, 57–64. [Google Scholar]
  12. O’Keefe, M.B.; Mans, D.M.; Kaelin, D.E.; Martin, S.F. Total Synthesis of Isokidamycin. J. Am. Chem. Soc. 2010, 132, 15528–15530. [Google Scholar]
  13. Mitsuru, S.; Hayashi, Y. Chemistry of Epoxyquinols A, B, and C and Epoxytwinol A. Eur. J. Org. Chem. 2007, 23, 3783–3800. [Google Scholar]
  14. Sparks, S.M.; Chen, C.-L.; Martin, S.F. Tandem Intramolecular Benzyne-furan cycloadditions. Total synthesis of vineomycinone B2 methyl ester. Tetrahedron 2007, 63, 8619–8635. [Google Scholar]
  15. Morton, G.E.; Barrett, A.G.M. Interative Benzyne-Furan Cycloaddition Reactions: Studies toward the Total Synthesis of ent-Sch 47554 and ent-Sch 47555. Org. Lett. 2006, 8, 2859–2861. [Google Scholar]
  16. Padwa, A.; Ginn, J.D. Studies on the Synthesis of (+/−)-Stenine: A Combined Intramolecular [4 + 2]- Cycloaddition/Rearrangement Cascade. J. Org. Chem. 2005, 70, 5197–5206. [Google Scholar]
  17. Takadoi, M.; Katoh, T.; Tadashi, I.; Terashima, S. Synthetic studies of himbacine, a potent agonist of the muscarinic M2 subtype receptor 1. Stereoselective total synthesis and antagonistic activity of enantiomeric pairs of himbacine and (2'S,6'R)-diepihimbacine, 4-epihimbacine, and novel himbacine congeners. Tetrahedron 2002, 58, 9903–9923. [Google Scholar]
  18. Mitsuru, S.; Yamaguchi, J.; Kakeya, H.; Osada, H.; Hayashi, Y. Total Synthesis of (+)-Epoxyquinols A and B. Angew. Chem. Int. Ed. 2002, 41, 3192–3194. [Google Scholar]
  19. Rae, R.L.; Żurek, J.M.; Paterson, M.J.; Bebbington, M.W.P. Halogenation effects in intramolecular furan Diels–Alder reactions: Broad scope synthetic and computational studies. Org. Biomol. Chem. 2013, 11, 7946–7952. [Google Scholar]
  20. Fleming, I. Molecular Orbitals and Organic Chemical Reactions; Wiley: Chichester, UK, 2010. [Google Scholar]
  21. Houk, K.N. The frontier molecular orbital theory of cycloaddition reaction. Acc. Chem. Res. 1975, 8, 361–369. [Google Scholar]
  22. Sauer, J.; Sustmann, R. Mechanistic aspects of Diels–Alder reactions: A critical survey. Angew. Chem. Int. Ed. Engl. 1980, 19, 779–807. [Google Scholar]
  23. Brown, F.K.; Houk, K.N.; Burnell, D.J.; Valenta, Z. Modeling Control of Facial Stereoselectivity. Diels–Alder Cycloadditions of Unsymmetrically Substituted Cyclopentadienes. J. Org. Chem. 1987, 52, 3050–3059. [Google Scholar]
  24. Robiette, R.; Marchand-Brynaert, J.; Peeters, D. Influence of Diene Substitution on Transition State Stabilization in Diels–Alder Reaction. J. Org. Chem. 2002, 67, 6823–6826. [Google Scholar]
  25. James, N.C.; Um, J.M.; Padias, A.B.; Hall, H.K., Jr.; Houk, K.N. Computational investigation of the competition between the concerted Diels–Alder reaction and formation of diradicals in reactions of acrylonitrile with non-polar dienes. J. Org. Chem. 2013, 78, 6582–6592. [Google Scholar]
  26. Crawford, K.R.; Bur, S.K.; Straub, C.S.; Padwa, A. Intramolecular Cyclization Reactions of 5-Halo and 5-Nitro-Substituted Furans. Org. Lett. 2003, 5, 3337–3340. [Google Scholar]
  27. Padwa, A.; Crawford, K.R.; Straub, C.S.; Pieniazek, S.N.; Houk, K.N. Halo Substituent Effects on Intramolecular Cycloadditions involving Furanyl Amides. J. Org. Chem. 2006, 71, 5432–5439. [Google Scholar]
  28. Pienazek, S.N.; Houk, K.N. The Origin of the Halogen Effect on Reactivity and Reversibility of Diels–Alder Cycloadditions involving Furan. Angew. Chem. Int. Ed. 2006, 45, 1442–1445. [Google Scholar]
  29. Nakamura, N.; Takahasi, I.; Yamada, S.; Dobashi, Y.; Kitagawa, O. Intramolecular Diels–Alder reaction of N-allyl 2-furoyl amides: Effect of strain and amide rotational isomerisation. Tetrahedron Lett. 2011, 52, 53–55. [Google Scholar]
  30. Tantillo, D.J.; Houk, K.N.; Jung, M.E. Origins of Stereoselectivity in Intramolecular cycloadditions of dienes and dienophiles linked by ester and amide tethers. J. Org. Chem. 2001, 66, 1938–1940. [Google Scholar]
  31. Williams, A. Free Energy Relationships in Organic and Bioorganic Chemistry; RSC Publishing: Cambridge, UK, 2003. [Google Scholar]
  32. Huisgen, R.; Schug, R. New Mechanistic Criterion for Early and Late Transition States. J. Am. Chem. Soc. 1976, 98, 7819–7821. [Google Scholar]
  33. Gajewski, J.J.; Peterson, K.B.; Kagel, J.R. Transition State Structure Variation in the Diels–Alder reaction from Secondary Deuterium Kinetic Isotope effects: The reaction of a nearly symmetrical diene and dienophile is nearly synchronous. J. Am. Chem. Soc. 1987, 109, 5545–5546. [Google Scholar]
  34. Frisch, M.J.T.G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09; Revision D.01; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  35. Nyden, M.R.; Petersson, G.A. Complete Basis Set Correlation Energies 1. The Asymptotic Convergence of Pair Natural Orbital Expansions. J. Chem. Phys. 1981, 75, 1843–1862. [Google Scholar]
  36. Ochterski, J.W.; Petersson, G.A.; Montgomery, J.A. A complete basis set model chemistry 5. Extensions to six or more heavy atoms. J. Chem. Phys. 1996, 104, 2598–2619. [Google Scholar]
  37. Montgomery, J.A.; Frisch, M.J.; Ochterski, J.W.; Petersson, G.A. A complete basis set model chemistry. VI. Use of density functional geometries and frequencies. J. Chem. Phys. 1999, 110, 2822–2827. [Google Scholar]
  38. Montgomery, J.A.; Frisch, M.J.; Ochterski, J.W.; Petersson, G.A. A complete basis set model chemistry. VII. Use of the minimum population localization method. J. Chem. Phys. 2000, 112, 6532–6542. [Google Scholar]
  • Sample Availability: Samples of the compounds are not available from the authors.

Share and Cite

MDPI and ACS Style

Żurek, J.M.; Rae, R.L.; Paterson, M.J.; Bebbington, M.W.P. Unusual Structure-Energy Correlations in Intramolecular Diels–Alder Reaction Transition States. Molecules 2014, 19, 15535-15545. https://doi.org/10.3390/molecules191015535

AMA Style

Żurek JM, Rae RL, Paterson MJ, Bebbington MWP. Unusual Structure-Energy Correlations in Intramolecular Diels–Alder Reaction Transition States. Molecules. 2014; 19(10):15535-15545. https://doi.org/10.3390/molecules191015535

Chicago/Turabian Style

Żurek, Justyna M., Robert L. Rae, Martin J. Paterson, and Magnus W. P. Bebbington. 2014. "Unusual Structure-Energy Correlations in Intramolecular Diels–Alder Reaction Transition States" Molecules 19, no. 10: 15535-15545. https://doi.org/10.3390/molecules191015535

APA Style

Żurek, J. M., Rae, R. L., Paterson, M. J., & Bebbington, M. W. P. (2014). Unusual Structure-Energy Correlations in Intramolecular Diels–Alder Reaction Transition States. Molecules, 19(10), 15535-15545. https://doi.org/10.3390/molecules191015535

Article Metrics

Back to TopTop