A Dft Pursuit on Tandem Cycloaddition-Cycloreversion of 2- Pyrone and 1,4-Oxazinone with Strained Alkynes

Received Date: Jul 29, 2020 / Accepted Date: Aug 27, 2020 / Published Date: Aug 29, 2020 Abstract 2-Pyrone (2P) and 1,4-oxazinone (OXZ) react with strained alkynes viz., cyclohexyne (SA1), benzyne (SA2) and pyridyne (SA3) through concerted cycloaddition-cycloreversion reaction. All the reactions have been scanned and the activation parameters, FOE, deformation energies and variation in aromaticity computed at B3LYP/6-31g(d) level.The obtained results have been compared with that of acetylene. In SA1,faster cycloadditionis explained by its strain, in addition,aromatic gain is accountable to enhance the rate of cycloaddition in SA2 and SA3.


Introduction
Tandem cycloaddition-cycloreversion reactions are synthetically important to produce enormous natural products. One might reasonably expect that increased strain will lead to enhanced reactivity of a molecule and thus molecules become good candidates here to perform tandem reactions. Literature apparently shows the effect of ring size on the reactivity of cyclic alkynes [1][2][3][4][5][6][7]. Other strained alkynes (SAs) include arynes which are the reactive intermediates derived from substituted arenes that participate in some aromatic nucleophilic substitutions [8,9] and some cycloaddition reactions [10]. The most extensively studied aryne is benzyne which is known to undergo [4+2] cycloaddition with dienes like oxazole and thiophene derivatives [11][12][13]. Cyclohexyne has been reported to undergo tandem cycloaddition-cycloreversion with 2-pyrone to form tetrahydronaphthalene. This strategy when operated with some complex pyrones yielded analogues of planar frame work of dynemicin-A [14]. The effect of heteroatom on tandem cycloadditioncycloreversion of 2-pyrone and 1,4-oxazinone with hetero alkynes has also been explored theoretically [15,16]. On the other hand,addition of benzynes and pyridynes with 2-pyrone have so far not been reported Pyrone and 1,4-Oxazinone with Strained Alkynes DOI: https://doi.org/10.36811/jca.2020.110006 JCA: August-2020: Page No: [10][11][12][13][14][15][16][17][18][19][20][21] Page: 11 www.raftpubs.com experimentally. This is the stimulus to have the aim of this work i.e., (i) to monitor the effect of strained triple bond of cyclohexyne and arynes over the tandem cycloaddition and cycloreversion, (ii) the effect of strain and aromatic gain on the reactivity in benzyne and pyridyne and (iii) the effect of strain, aromatic gain and the effect of hetero atom in pyridyne reaction. This study discusses the computed results such as activation parameters, FO analyses, deformation energies, variation in aromaticity along the reaction coordinates of the transition states formed during the tandem cycloaddition-cycloreversion reaction of 2P and OXZ with strained alkynes viz., cyclohexyne (SA1), benzyne (SA2) and pyridyne (SA3) and compared with that of acetylene.

Computational Methods
All the reactions have been modelled and scanned using Gaussian 98 program [17] at B3LYP [18,19]/6-31G (d) level. The magnetic susceptibility isotropy has been calculated by computing the NMR shielding tensors at B3LYP/6-311+G (2d, p) using IGAIM method [20,21] which is a slight variation of CSGT [20,21] method. Natural bond orbital (NBO) calculations have been performed. For both the cycloaddition and cycloreversion processes, the variation in magnetic susceptibility isotropy (χiso) along the reaction coordinate has been monitored. Change in aromaticity during the reactions has been computed.

Results and Discussion
Cycloaddition and cycloreversion of 2P and OXZ with strained alkynes SA1, SA2 and SA3 form benzo fused aromatics which are illustrated in Scheme 1 and 2. The relative free energy profile for the sequence of 2P and OXZ with strained alkynes are represented in figure 1-6 along with the optimized geometries and notable bond distances. There is the possibility of two extruding groups from the cycloadduct viz., CO2 and HCN in the case of cycloreversion of strained alkyne with OXZ. When SA3 is the dienophile there are 2 regioisomeric adducts possible, denoted as Regio 1 and Regio 2 in Scheme 1(c)& 2(c).
The cycloaddition of the three SAs lead to highly stable and exothermic cycloadducts viz., A1, A2 and R1A3 & R2A3. The exothermicity of these cycloadducts are (ΔHr= -65.55, -80.03, -77.27 kcal mol -1 of A1, A2 and R2A3 respectively) two to three fold greater than that of simple acetylene (ΔHr= -29.38 kcal mol -1 ) due to two factors: (i) Release of strain in the triple bond (ii) Aromaticity gain in the adduct. It is revealed by the formation of strainless aromatic benzene and pyridine fused with bicyclic lactones A2 and A3 respectively, in 2P-SA2 and 2P-SA3 reactions. In contrast, there is only strain release, not the aromatic gain in 2P-SA1 cycloaddition. The barrier to cycloreversion shown in table 1. indicates that A1 extrudes CO2 faster and is comparable with that of 2P-Ace reaction. On the other hand, A2 and A3 (ΔG ‡ =15.95 kcal mol -1 and 16.91 kcal mol -1 ) have their cycloreversion barrier about 4-5 kcal mol -1 higher than 2P-SA1 reaction. This is ascertained to the excessive thermodynamic stability of the cycloadduct. Extra stabilization Page: 12 www.raftpubs.com in A2 and A3 due to aromatic gain also results in higher barrier to cycloreversion.

OXZ-SAs reactions
The barrier to cycloaddition is reduced for OXZ-SAs reactions also, as observed in OXZ-Ace. Reduction in barrier in OXZ-SAs reactions (ΔG ‡ =4.5 kcal mol -1 , 3.3 kcal mol -1 and 2.8 kcal mol -1 for SA1, SA2 and SA3 respectively) are less than that of simple alkyne case (ΔG ‡ =7 kcal mol -1 ). Nitrogen in the system facilitates the reaction due to fewer aromatic loss and its remote location,from the reacting triple bond in SA3, affects the barrier to the least extent. Now the reactivity order of the SAs with OXZ is SA2>SA3>SA1>Ace. The reduction in free energy of activation of the OXZ-SAs reaction presented in Table 2. It is ascribed to a much lesser activation enthalpy compared to 2P-Ace case shown in Table 1. It is notable that ΔH ‡ for the OXZ-SA2 and OXZ-SA3 reactions are slightly negative. ΔS ‡ for cycloaddition and cycloreversion resemble that of 2P-Ace and 2P-SAs reactions. In OXZ-SAs reactionsSA2 and SA3 form cycloadduct A4 and A5, whereas SA3 makes two regioisomeric adducts R1A6 and R2A6. Their stability order is A4<R1A6<A5<R2A6. Out of two pathways, similar to OXZ-Ace reaction, CO2 extrusion is facile than HCN elimination and the difference their barriers is nearly 15-20 kcal mol -1 . Cycloreversions of A4, A5, R1A6 and R2A6 (ΔG ‡ =14.33, 19.69, 19.37 and 21.34 kcal mol -1 respectively) are similar to 2P-Ace reaction.

Frontier Orbital Analyses
Strained alkynes have their HOMO raised and LUMO lowered than simple alkyne and hence show higher reactivity. FOE gaps in figure 7 show that the 2P-SA1 cycloaddition follows Inverse Electron Demand (IED), similar to 2P-Ace. Added unsaturations and hetero atom stabilize the FOs of SA2 and SA3 respectively. As their LUMOs are lowered, SA2 and SA3 react by normal electron demand (NED) type. These FOE results are supported by QCT at the dienophile in the 2P-SA2 and 2P-SA3 reactions whereas 2P-SA1 reaction contradicts the FOE prediction, as its QCT is negative. In OXZ-SAs reactions, cycloaddition of OXZ-SA1 is IED type, similar to 2P-SA1 and is supported by QCT. Reduction in FOE gap indicates fall in the activation barrier to OXZ-SA1 cycloaddition. OXZ-SA2 reaction is IED, in contrast QCT predicts it to be a NED type as -0.028 charge is on the dienophile. FOE gaps show OXZ-SA2 reaction is neutral electron demand type whereas its QCT value is -0.047, representing it to be NED. In OXZ-SA3 and 2P-SA1 reactions, FOE gaps predicts IED type while QCT values indicate as NED type reactions.

Effect of strain on reactivity
The barrier to cycloaddition for 2P-SAs is less than 2P-Ace reaction, recognizing the reactivity (deformation energy) promoted by the strained triple bond. The deformation energies of diene and dienophile (DE1and DE2 respectively), in the corresponding transition state are presented in table 3. In 2P-Ace reaction, both diene and dienophile are largely deformed whereas drastically reduced in 2P with strained alkynes which has the order as SA1>SA3>SA2. This aligns their barrier trend.Further, in OXZ reaction N in the embedded diene still reduces the deformation energy of both diene and dienophile.   Page: 13 www.raftpubs.com difference between the initial angles (viz., C8-C7-C12 and C7-C8-C9) in the alkyne and that of cycloadduct. Comparison of variation of C8-C7-C12 and C7-C8-C9 angles for the strained and simple alkyne with 2P along the reaction coordinate of cycloaddition TS is presented in figure 11. Greater the difference in angle, greater is the barrier and variation in angle is sharp in 2P-Ace whereas very slight in 2P-SAs reactions.

Variation of aromaticity along the reaction coordinates
Aromatic gain during the reactions will definitely influence the rate. Among the three SAs involved reactions, SA2 and SA3 form benzo and pyrido-fused cycloadducts respectively, while SA1does not. The variation of aromaticity has been monitored by plotting the χiso along the reaction coordinate of both cycloaddition and cycloreversion of TSs. These plots of 2P-SA1 and 2P-SA2 reactions are collected in figures 14&15 respectively. Similar plots of 2P-SA3 and OXZ-SAs have not been done, as it resembles that of 2P-SA2.

2P-SA1 reaction
The change in χiso of TS1 (in the region -1 to +1 (amu) 1/2 .Bohr) along the reaction coordinate is shown in figure 12(a). The curve indicates a continuous aromatic enhancement during the reaction, at the TS there is an increase (more negative χiso value representing aromatic enhancement) while the rate of aromaticity gain is gradually decreased. Figure 6 (b) & (c) show the variation of χiso of pyrone and cyclohexyne ring, indicating a decrease in negative χiso till TS followed by a steeper increase. It explains both the rings lose aromaticity during cycloaddition and hence the continual aromaticity increase in plot 6(a) is due to the pericyclic closed curve formed. In this case, gain in aromaticity even after the TS can be explained as follows. The reacting LUMO of SA1 is perfectly oriented for good overlap with incoming diene. Also the change in angles are less along the reaction coordinate as the SA1is nearly close to its geometry in the adduct. This is reflected in deformation energy (DE2) TS being very low for the SAs in general. Therefore, even after the TS, deformation of SA1 will be necessarily low thereby increase in the aromaticity. During cycloreversion, except for the absolute value of χiso which is very high due to the fused benzene, the trend remains the same. Both plots 6 (d) & (e) show a linear increase in χiso owing to the formation of tetrahydronaphthalene (P1) after CO2 extrusion and hence the aromaticity of newly formed ring dominates the scenario.

2P-SA2 reaction
The variation of χiso during 2P-SA2 reaction shown in figure 7 indicates a profile similar to 2P-SA1 reaction. In this case, a continual linear increase in χiso even after the TS is observed. It is ascribed to the aromatic gain by the benzene ring (plot 7 (c)). The pyrone ring fragment (plot 7 (b)) loses aromaticity which is compared with a low DE1 in 2P-SA2 reaction. Hence the least deformation of the pyrone ring is reflected in aromaticity monitored as variation in χiso. During cycloreversion, negative value of χiso increases due to CO2 extrusion to form naphthelene. Therefore, unlike in a typical pericyclic reaction, there is a continual increase of aromaticity. Pyrone and 1,4-Oxazinone with Strained Alkynes DOI: https://doi.org/10.36811/jca.2020.110006 JCA: August-2020: Page No: 10-21 Page: 14 www.raftpubs.com   Page: 16 www.raftpubs.com

Conclusions
Tandem cycloaddition-cycloreversion of 2pyrone and 1,4-oxazinone with strained alkynes namely cyclohexyne(SA1), benzyne(SA2) and pyridine(SA3) occurs through concerted mechanism with slight asynchronicity due to the unsymmetric diene. Among the strained alkynes, benzyne reacts faster and pyridyne follows it, while cyclohexyne does slowly. All the three strained alkynes react extremely faster than Acetylene (Ace). Cycloreversion of 2P-SA1 adduct is comparable to 2P-Ace. In SA2 & SA3, cycloreversion is the rate determining step owing to the larger stability of adducts by strain release and aromatic gain by small ΔH ‡ and magnetic anisotropy variation. Aza diene cycloadds faster with strained alkynes whereas cycloreverts slower than the non aza diene. Deformation analysis explains extremely faster cycloaddition of SA1 with 2P & OXZ is due to its strain. In SA2 & SA3, strain and aromatic gain enhance the rate of cycloaddition and is demonstrated by variation of aromaticity. FO analysis shows the normal and inverse electron demands in these dienophiles especially in benzyne and pyridyne. Overall, the study proves that release of strain and aromaticity gain governs the tandem reaction of strained alkynes.