Computational Exploration of Concerted and Zwitterionic Mechanisms of Diels\xe2\x80\x93Alder Reactions between 1,2,3-Triazines and Enamines and Acceleration by Hydrogen-Bonding Solvents
The mechanisms of Diels\xe2\x80\x93Alder reactions between 1,2,3-triazines and enamines have been explored with density functional theory computations. The focus of this work is on the origins of the different reactivities and mechanisms induced by substituents and by hexafluor-oisopropanol (HFIP) solvent. These inverse electron-demand Diels\xe2\x80\x93Alder reactions of triazines have wide applications in bioorthogonal chemistry and natural product synthesis. Both concerted and stepwise cycloadditions are predicted, depending on the nature of substituents and solvents. The nature of zwitterionic intermediates and the mechanism by which HFIP accelerates cycloadditions with enamines are characterized. Our results show the delicate nature of the concerted versus stepwise mechanism of inverse electron-demand Diels\xe2\x80\x93Alder reactions of 1,2,3-triazines, and that these mechanisms can be altered by electron-withdrawing substituents and hydrogen-bonding solvents.
INTRODUCTION<!>COMPUTATIONAL DETAILS<!>Uncatalyzed Diels\xe2\x80\x93Alder Reaction.<!>Solvent-Catalyzed Diels\xe2\x80\x93Alder Reaction.<!>FMO Analysis of the Triazine Diels\xe2\x80\x93Alder Reactions.<!>CONCLUSION
<p>The inverse electron-demand Diels–Alder reactions of heterocyclic azadienes are powerful approaches to construct six-membered heterocycles.1 These reactions have wide applications in bioorthogonal chemistry2 and natural product synthesis.3 Both Ess and Bickelhaupt and our group have systematically investigated the cycloadditions of benzene and ten different azabenzenes (pyridine, three diazines, three triazines, and three tetrazines) with the ethylene dienophile (Scheme 1a).4 Each nitrogen substitution causes higher Diels–Alder reactivity of azapolyenes due to more favorable orbital interactions and lower distortion energies.4b Through these computational studies, we now have a better under-standing of their intrinsic reactivities, especially on roles played by the number and position of the heteroatoms. While these simple reactions are all predicted to proceed via a concerted mechanism, substituents and solvents could change the mechanism and have various effects on their reactivity. For example, electron-withdrawing substituents on the C-5 position of 1,2,3-triazine have been found to increase its reactivity in inverse-electron demand Diels–Alder reactions.5b</p><p>Recently, the cycloadditions of a 1,2,3-triazine5 with an enamine was used as a key step in the concise total synthesis of the natural product methoxatin6 by the Boger group (Scheme 1b and 1c). It was discovered that non-nucleophilic fluorinated alcohol solvents, such as hexafluoroisopropanol (HFIP), accelerate these cycloadditions. The reaction of 1,2,3-triazine 1a with enamine 2 has a low yield in CHCl3 (Scheme 2), but has a much higher yield with HFIP as the solvent. Similarly, for the reaction of 1,2,3-triazine 1b with enamine 2, the use of HFIP as the solvent provides the cycloadduct 3b in good yield (86%). Typically, common Lewis acid catalysts are ineffective for inverse-electron demand Diels–Alder reaction because of the nonproductive consumption of the electron-rich dienophiles.7 Hydrogen-bond donation by HFIP appears to accelerate the reaction without destroying the enamine. Solvent-acceleration of inverse electron-demand Diels–Alder reactions with electron-deficient heterocyclic azadienes could be attributed to many factors, such as hydrogen-bonding and polarity of the solvent.8 We have used density functional theory (DFT) to explore the mechanisms of cycloadditions between several 1,2,3-triazines and enamines in different solvents.</p><!><p>All density functional theory (DFT) calculations were performed using Gaussian 09.9 Geometry optimizations and frequency calculations were performed at the M06–2X10/6–31G(d) 11 level of theory. Normal mode vibrational analysis confirmed that optimized structures are minima or transition structures. Single-point energies were computed on the optimized geometries with the conductor-like polarizable continuum model (CPCM)12 using chloroform (ε = 4.7) or HFIP (ε = 16.7) as solvent at the M06–2X/6–311+G(d,p) level of theory. An explicit HFIP was also included. Geometry optimizations and subsequent frequency calculations using M06–2X/6–311+G(d,p) with solvent produce qualitatively the same result, and were not carried out for all reactions.13 The frontier molecular orbitals (FMOs) and their energies were computed at the HF/6–311+G(d,p) level using the M06–2X/6–31G(d)-optimized geometries. All 3D renderings of stationary points were generated using CYLview.14</p><!><p>We first studied the uncatalyzed Diels–Alder reaction of parent 1,2,3-triazine 1a with enamine 2. Figure 1 shows the free energy profiles for the endo and exo pathways and the corresponding energetics of the reaction in CHCl3, while Figure 2 shows the transition structures. Both endo-TS1 and exo-TS1 are highly asynchronous concerted transition structures, each with a shorter forming C–C bond at 1.6–1.7 Å and a longer forming C–N bond at 2.3–2.4 Å. We showed earlier that C–C bond formation is always preferred over C–N bond formation in polyazine cycloadditions.4b The barriers for the Diels–Alder step through endo and exo pathways are 32.0 and 33.1 kcal/mol, respectively. The Diels–Alder adducts endo-1 and exo-1 are only about 7 kcal/mol more stable than endo-TS1 and exo-TS1, respectively, but the subsequent release of N2 is very exergonic, and has only low barriers (1.3 kcal/mol via exo-TS2 and 2.2 kcal/mol via endo-TS2). We earlier explored the rapid dynamics of N2 loss from similar intermediates.15 The low barrier implies a lifetime of a few 100 fs at most. The subsequent elimination of amine that leads to the aromatized compound 3a has a barrier of 38.0 kcal/mol (Figure S2). With the assistance of HFIP, the barrier is decreased to 24.4 kcal/mol (Figure S2).</p><p>The two electron-withdrawing groups (−CO2Et) at the 4- and 6-positions of 1,2,3-triazine 1b were modeled as CO2Me groups computationally. Figure 3 shows the free energy profiles for the endo and exo pathways of the Diels–Alder reaction of 1,2,3-triazine 1c with enamine 2, and Figure 4 shows transition structures and intermediates. Unlike the concerted Diels–Alder reaction of 1,2,3-triazine 1a with enamine 2, this reaction of 1,2,3-triazine 1c with enamine 2 undergoes a stepwise mechanism, forming the C–C bond first. The optimized transition structures for the first steps of endo and exo pathways have forming C–C bonds of ~1.8 Å, while the C–N distances are ~2.9 Å in TS3s. Shallow intermediates endo-3 and exo-3 are located on the potential energy surface for both the endo and the exo pathways, and they have comparable energies to the C–C bond forming transition structures endo-TS3 and exo-TS3. The forming C–C bonds are ~1.6 Å and the forming C–N bonds are ~2.3 Å in both second transition structures, endo-TS4 and exo-TS4. Interestingly, the C–N bond forming transition structures endo-TS4 and exo-TS4 are very similar to the concerted transition structures endo-TS1 and exo-TS1. The second barrier is only ~5 kcal/mol, and is the rate-determining step. Transition structure exo-TS4 is slightly favored over endo-TS4. The subsequent release of N2 is very exergonic, and the transition structures are lower in energy (22.5 kcal/mol for exo-TS5 and 26.0 kcal/mol for endo-TS5) than those of the second Diels–Alder step.</p><p>We have characterized the zwitterionic intermediates, endo-3 and exo-3. The charge separations on the triazine and enamine fragments, are 0.68 and 0.72 in the endo and exo species, respectively. The electrostatic potential (ESP) maps of endo-3 and exo-3, in Figure 4, show the charge separation graphically.</p><p>Compared with the reaction of concerted Diels–Alder reaction of 1,2,3-triazine 1a with enamine 2, the Diels–Alder reaction of 1,2,3-triazine 1c with enamine 2 undergoes a stepwise mechanism, with the C–N bond forming as the rate-limiting step. The overall barrier is lowered by 5.0 kcal/mol (27.0 kcal/mol vs 32.0 kcal/mol) by the two ester substituents. This is in agreement with the higher reactivity of 1,2,3-triazine 1b found in the experiment.6</p><!><p>We next explored the reaction of triazine and the diester-substituted triazine in HFIP. One explicit HFIP molecule was included. The energies and structures for the parent triazine are shown in Figures 5 and 6. The hydrogen-bonding of one molecule of HFIP solvent with the N-1 of 1,2,3-triazine 1a leads to complex 1a-HFIP. This association process is 5.1 kcal/mol endergonic but 8.1 kcal/mol exothermic (1a + HFIP → 1a-HFIP, ΔG = 5.1 kcal/mol, ΔH = −8.1 kcal/mol) when we use the normal standard state for gas-phase reactions (1 atm, 298.15K). Single-point energies with solvation were computed on the M06–2X/6–31G(d)-optimized structures in the gas phase. Correcting for 1 M standard state for solution-phase reactions, and the 9.5 M concentration of pure HFIP solvent,16 the ΔG becomes 1.9 kcal/mol. In addition, computed entropies of association in solution using harmonic assumption are frequently too negative.17 Considering all these factors, we believe the association process to form the complex 1a-HFIP is favorable in solution, and we have set the complex energy to zero.</p><p>Figure 5 shows the free energy profiles for the endo and exo pathways of the reaction. The HFIP-catalyzed Diels–Alder reaction of the parent triazine now proceeds via a stepwise mechanism very similar to the reaction of the substituent-activated triazine. As shown in Figure 6, the forming C–C bond is ~1.9 Å and the C–N distance is ~3.0 Å in both endo-TS6 and exo-TS6, and ~1.6 Å and ~2.1 Å in both endo-TS7 and exo-TS7, the second step of the reaction. The zwitterionic intermediates endo-6 and exo-6 are located on the potential energy surface for both the endo and the exo pathways, and they are ~3 kcal/mol more stable than endo-TS6 and exo-TS6. The barrier of reaction (the second step) is lowered by 4.5 kcal/mol versus the single TS for the uncatalyzed one. The hydrogen-bonding distance between H of HFIP and N1 of triazine 1a in zwitterionic intermediates endo-6 and exo-6 are shorter than in either complex 1a-HFIP (1.85 Å for H···N distance) or transition structures that forming C–C and C–N bonds. The barriers for forming the C–N bond are once again higher than that for forming the C–C bond. Transition structure exo-TS7 is very close to endo-TS7 in free energy. The subsequent release of N2 is rapid as always.</p><p>The hydrogen-bonding of one molecule of HFIP solvent with the N1 of the activated 4,6-dicarbomethoxy 1,2,3-triazine 1c leads to complex 1c-HFIP. Figure 7 shows the free energy profiles for the endo and exo pathways of the reaction, while Figure 8 shows transition structures and intermediates. The Diels–Alder reaction of 1c with enamine 2 undergoes stepwise mechanism in either CHCl3 or HFIP. The profile in Figure 7 is like the previous stepwise profile. The HFIP stabilizes the first transition state by 9 kcal/mol and the zwitterionic intermediate by 11 kcal/mol. The hydrogen-bonding distance between H of HFIP and N1 of triazine 1c in zwitterionic intermediates endo-9 and exo-9 (~1.6 Å for H···N distance) are shorter than in both complex 1c-HFIP (2.06 Å for H···N distance) and in transition structures that forming C–C and C–N bonds. The charge separation has increased to 0.82 in the favored intermediate exo-9, and the ESP map in Figure 8 reflects the close association of negative and positive regions of the intermediate. The subsequent step of N2 release now becomes the rate-limiting step, although it is still very facile. HFIP has lowered the rate-determining step by 6 kcal/mol, for an approximately 104 acceleration.</p><p>Because there are other hydrogen bond acceptors in 1c besides N1, we have also explored the possibility of HFIP binding to N2, N3, −CO2Me of C4, and −CO2Me of C6 of 1,2,3-triazine 1c (Figure 9). We also explored hydrogen-bonding to the nitrogen of the enamine, but that destabilize the transition state. The HFIP binding to N2, −CO2Me of C4, and −CO2Me of C6 of 1,2,3-triazine 1c is less favorable compared to binding to N1 of 1,2,3-triazine 1c. It is worth noting that transition structure for forming the C–C bond, endo-TS9C, with HFIP hydrogen bonding to N3 of 1,2,3-triazine 1c, is slightly favored over endo-TS9 by 0.7 kcal/mol. The hydrogen-bonding interaction with N3 of 1,2,3-triazine 1c is preferred. For the transition structure for forming the C–N bond, endo-TS10C, with HFIP hydrogen bonding to N3 of 1,2,3-triazine 1c, is disfavored by 3.1 kcal/mol compared to endo-TS10. There may be two or more HFIP molecules hydrogen bonding to the triazine, but it is difficult to evaluate the entropy change with association reaction involving multiple molecules. Because of the polarity of transition states and intermediates, we have also reoptimized them in HFIP solvent (see Figure S1 in Supporting Information). Only minor changes were observed.</p><!><p>Figure 10a shows the computed HOMO of enamine 2, the lowest vacant π orbitals (LUMO+1) of triazines 1a and 1c and HFIP-complexed triazines 1a-HFIP and 1c-HFIP. For all the four LUMO+1 vacant orbitals of triazines, there are two nodes in the triazine ring in each of them. Hydrogen bonding to an HFIP molecule lowers the energy of the LUMO+1 orbital by 0.33 eV (1a-HFIP). With two electron-withdrawing ester substituents at 4-and 6-positions, the energy of the LUMO+1 vacant orbital lowers by 0.60 eV (1c). The LUMO +1 vacant orbital of the HFIP-complexed complex 1c-HFIP is further decreased to an energy of 0.85 eV.</p><p>A plot of the highest Diels−Alder reaction barrier versus the HOMO–LUMO gap is shown in Figure 10b. There is a very good correlation between the barrier and the HOMO–LUMO gap, and the HFIP-complexed triazine 1c-HFIP is most reactive due to its lower energy of LUMO+1 compared to that of 1a, 1c and 1a-HFIP.</p><!><p>We have explored the mechanisms of Diels–Alder reactions between 1,2,3-triazines and enamine using DFT computations and investigated the origins of the different reactivities induced by substituents or the HFIP solvent. While the Diels–Alder reaction of unsubstituted 1,2,3-triazine and the enamine undergo a concerted mechanism, either electron-withdrawing groups or the strongly hydrogen-bonding solvent stabilize zwitterionic intermediates and lead to stepwise mechanisms. Our results imply the subtle nature of the concerted versus stepwise mechanism of inverse electron-demand Diels–Alder reactions of 1,2,3-triazines, which could be altered by electron-withdrawing substituents and hydrogen-bonding solvents. Systematic studies of the full suite of heterocyclic azadienes are underway in our research group and will be reported in due course. With a more complete picture of the mechanisms of these reactions, we envision that the power of tuning substituents and solvents could be the power of tuning substituents and solvents could be the power of tuning substituents and solvents could be harnessed more efficiently by experimentalists, to develop new variants and to discover new applications of inverse electron-demand Diels–Alder reactions.</p>