Mechanism and structural aspects of thermal Curtius rearrangement. Quantum chemical study

The electronic and geometric structures of formyl, acetyl, and benzoyl azides were studied and fragments of the potential surfaces for the thermal Curtius rearrangement of these azides into the corresponding isocyanates were calculated by density functional theory at the PBE/TZ2P level. Acyl azides adopt two stable, conformations syn and anti, with respect to the C-N bond. The syn conformers are more stable than their anti analogs. The activation energies of the syn-anti isomerization in the series under study are 9.4, 7.0, and 9.2 kcal mol⁻¹, respectively, and the activation energies of the reverse reaction are 8.5, 6.1, and 2.5 kcal mol⁻¹. The rearrangement of syn-acyl azides is a one-step process, in which elimination of N₂ occurs synchronously with the rearrangement of atoms and bonds to form isocyanates. The activation energies of the rearrangements of syn-HC(O)N₃, syn-MeC(O)N₃, and syn-PhC(O)N₃ are 28.0, 32.9, and 34.5 kcal mol⁻¹, respectively. The rearrangement of the anti conformers of the above-mentioned azides involves the formation of singlet acylnitrene. The activation energies of the latter process are 34.6, 32.9, and 32.3 kcal mol⁻¹, respectively. The activation energies of the rearrangement of acylnitrenes into isocyanates are 20.9, 18.9, and 13.6 kcal mol⁻¹, respectively. The energy characteristics of the process and the structural data for the starting compounds, final products, and transition states provide evidence that the thermal Curtius rearrangement occurs predominantly by a concerted mechanism. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 10, pp. 2200–2209, October, 2005.     

A comparison of acetyl- and methoxycarbonylnitrenes by computational methods and a laser flash photolysis study of benzoylnitrene

Density functional theory (DFT), CCSD(T), and CBS-QB3 calculations were performed to understand the chemical and reactivity differences between acetylnitrene (CH(3)C(=O)N) and methoxycarbonylnitrene (CH(3)OC(=O)N) and related compounds. CBS-QB3 theory alone correctly predicts that acetylnitrene has a singlet ground state. We agree with previous studies that there is a substantial N-O interaction in singlet acetylnitrene and find a corresponding but weaker interaction in methoxycarbonylnitrene. Methoxycarbonylnitrene has a triplet ground state because the oxygen atom stabilizes the triplet state of the carbonyl nitrene more than the corresponding singlet state. The oxygen atom also stabilizes the transition state of the Curtius rearrangement and accelerates the isomerization of methoxycarbonylnitrene relative to acetylnitrene. Acetyl azide is calculated to decompose by concerted migration of the methyl group along with nitrogen extrusion; the free energy of activation for this concerted process is only 27 kcal/mol, and a free nitrene is not produced upon pyrolysis of acetyl azide. Methoxycarbonyl azide, on the other hand, does have a preference for stepwise Curtius rearrangement via the free nitrene. The bimolecular reactions of acetylnitrene and methoxycarbonylnitrene with propane, ethylene, and methanol were calculated and found to have enthalpic barriers that are near zero and free energy barriers that are controlled by entropy. These predictions were tested by laser flash photolysis studies of benzoyl azide. The absolute bimolecular reaction rate constants of benzoylnitrene were measured with the following substrates: acetonitrile (k = 3.4 x 10(5) M(-1) (s-1)), methanol (6.5 x 10(6) M(-1) s(-1)), water (4.0 x 10(6) M(-1) s(-1)), cyclohexane (1.8 x 10(5) M(-1) s(-1)), and several representative alkenes. The activation energy for the reaction of benzoylnitrene with 1-hexene is -0.06 +/- 0.001 kcal/mol. The activation energy for the decay of benzoylnitrene in pentane is -3.20 +/- 0.02 kcal/mol. The latter results indicate that the rates of reactions of benzoylnitrene are controlled by entropic factors in a manner reminiscent of singlet carbene processes.     

Isomerization pathways from the norbornadiene to the cycloheptatriene radical cation by opening a bridgehead-methylene bond: a theoretical investigation

Three skeletal rearrangement channels for the norbornadiene (N*+) to the 1,3,5-cycloheptatriene (CHT*+) radical cation conversion, initialized by opening a bridgehead-methylene bond in N*+, are investigated using the quantum chemical B3LYP, MP2 and CCSD(T) methods in conjunction with the 6-311 +G(d,p) basis set. Two of the isomerizations proceed through the norcaradiene radical cation (NCD*+), either through a concerted path (N*+ - NCD*+), or by a stepwise mechanism via a stable intermediate (N*+ - I1 - NCD*+). At the CCSD(T)/6-311 +G(d,p)//B3LYP/6-311 +G(d,p) level, the lowest activation energy, 28.9 kcal mol(-1), is found for the concerted path whereas the stepwise path is found to be 2.3 kcal mol(-1) higher. On both pathways, NCD*+ rearranges further to CHT*+ with significantly less activation energy. The third channel proceeds from N*+ through I1 and then directly to CHT*+, with an activation energy of 37.1 kcal mol(-1). The multi-step channel reported earlier by our group, which proceeds from N*+ to CHT*+ via the quadricyclane and the bicyclo[2.2.1]hepta-2-ene-5-yl-7-ylium radical cations, is 4.6 kcal mol(-1) lower than the most favorable path of the present study. If the methylene group is substituted with C(CH3)2, however, the concerted path is estimated to be 5.6 kcal mol(-1) lower than the corresponding substituted multi-step path at the B3LYP/6-311+(d,p) level. This shows that substitution of particular positions can have dramatic effects on altering reaction barriers in the studied rearrangements. We also note that identical energies are computed for CHT*+ and NCD*+ whereas, in earlier theoretical investigations, the former was reported to be 6-17 kcal mol(-1) more stable than the latter. Finally, a bent geometry is obtained for CHT*+ with MP2/6-311 +G(d,p) in contradiction with the planar conformation reported for this cation in earlier computational studies.     

Theoretical study of the electrocyclic ring closure of hydroxypentadienyl cations

At the 6-311G* level of theory, DFT methods predict that the rearrangement of 1,4-dihydroxy-5-methylpentadienyl cation 1 (R = Me) to protonated trans-3-hydroxy-2-methylcyclopent-4-en-1-one 2, an intermediate step in the Piancatelli reaction or rearrangement of furfuryl carbinols to trans-2-alkyl(aryl)-3-hydroxycyclopent-4-en-1-one, is a concerted electrocyclic process. Energetic, magnetic, and stereochemical criteria are consistent with a conrotatory electrocyclic ring closure of the most stable out,out-1 isomer to afford trans-2. Although the out,in-1 isomer is thermodynamically destabilized by 6.84 kcal mol(-1), the activation energy for its cyclization is slightly lower (5.29 kcal mol(-1) versus 5.95 kcal mol(-1)). The cyclization of the isomers of 1 with the C1-hydroxy group inwards showed considerably higher activation energies than their outwards counterparts. in,out-1, although close in energy to out,out-1 (difference of 1.57 kcal mol(-1)) required about 10 kcal mol(-1) more to reach the corresponding transition structure. The value measured for the activation energy of in,in-1 (17.32 kcal mol(-1)) eliminates the alternative conrotatory electrocyclization of this isomer to provide trans-2. Geometric scrambling by isomerization of the terminal C1--C2 bond of 1 is also unlikely to compete with electrocyclization. The possibility to interpret the 1-->2 reaction as a nonpericyclic cationic cyclization was also examined through NBO analysis, and the study of bond lengths and atomic charges. It was found that the 1-->2 concerted rearrangement benefits from charge separation at the cyclization termini, an effect not observed in related concerted electrocyclic processes, such as the classical Nazarov reaction 3-->4 or the cyclization of the isomeric 2-hydroxypentadienyl cation 5.     

Quadricyclane radical cation rearrangements: a computational study of the transformations to 1,3,5-cycloheptatriene and norbornadiene

An alternative skeletal rearrangement of the quadricyclane radical cation (Q*+) explains the side products formed in the one-electron oxidation to norbornadiene. First, the bicyclo[2.2.1]hepta-2-ene-5-yl-7-ylium radical cation, with an activation energy of 14.9 kcal mol(-1), is formed. Second, this species can further rearrange to 1,3,5-cycloheptatriene through two plausible paths, that is, a multistep mechanism with two shallow intermediates and a stepwise path in which the bicyclo[3.2.0]hepta-2,6-diene radical cation is an intermediate. The multistep rearrangement has a rate-limiting step with an estimated activation energy of 16.5 kcal mol(-1), which is 2.8 kcal mol(-1) lower in energy than the stepwise mechanism. However, the lowest activation energy is found for the Q*+ cycloreversion to norbornadiene that has a transition structure, in close correspondence with earlier studies, and an activation energy of 10.1 kcal mol(-1), which agrees well with the experimental estimate of 9.3 kcal mol(-1). The computational estimates of activation energies were done using the CCSD(T)/6-311+G(d,p) method with geometries optimized on the B3LYP/6-311+G(d,p) level, combined with B3LYP/6-311+G(d,p) frequencies.     

5-羟基色胺受体激动剂(1R,2S)-(-)-2-(2-羟基苯基)-N,N-二丙基环丙胺的不对称合成

利用手性双唑啉与三氟甲磺酸亚铜催化2-甲氧基苯乙烯与重氮乙酸二环己基甲酯的不对称环丙烷化反应合成了手性环丙烷羧酸酯,用氢氧化钠对其进行选择性水解得到全反式环丙烷羧酸,其ee值经GC测定为88%.进一步经过Curtius重排、烷基化等反应及重结晶等步骤合成了光学纯的具有生理活性的环丙胺化合物1a和1b.     

Computational study of the Curtius-like rearrangements of phosphoryl, phosphinyl, and phosphinoyl azides and their corresponding nitrenes

The free energies of reaction (DeltaG) and activation (DeltaG) were determined for the Curtius-like rearrangement of dimethylphosphinoyl, dimethylphosphinyl, and dimethylphosphoryl azides as well as the corresponding singlet and triplet nitrenes by CBS-QB3 and B3LYP computational methods. From CASSCF calculations, it was established that the closed-shell configuration was the lower energy singlet state for each of these nitrenes. The triplet states of dimethylphosphinyl- and dimethylphosphorylnitrene are the preferred ground states. However, the closed-shell singlet state is the ground state for dimethylphosphinoylnitrene. The CBS-QB3 DeltaG values for the Curtius-like rearrangements of dimethylphosphinyl and dimethylphosphoryl azides were 45.4 and 47.0 kcal mol-1, respectively. For the closed-shell singlet dimethylphosphinyl- and dimethylphosphorylnitrene, the CBS-QB3 DeltaG values for the rate-limiting step of the Curtius-like rearrangement were 22.9 and 18.0 kcal mol-1, respectively. It is unlikely that the nitrenes will undergo a Curtius-like rearrangement because of competing bimolecular reactions that have lower activation barriers. The pharmacology of weaponized organophosphorus compounds can be investigated using phosphorylnitrenes as photoaffinity labels. Dominant bimolecular reactivity is a desirable quality for a photoaffinity label to possess, and thus, the resistance of phosphorylnitrenes to intramolecular Curtius-like rearrangements increases their usefulness as photoaffinity labels.     

A novel class of compounds: synthesis of 5,5′-carbonyl-bis(5,6-dihydro-4H-furo- and thieno-[2,3-c]pyrrol-4-ones)

We hereby report the first synthesis of novel class of compounds, 5,5′-carbonyl-bis(5,6-dihydro-4H-thieno- and furo-[2,3-c]pyrrol-4-one starting from methyl 2-(2-methoxy-2-oxoethyl) thiophene- and furan-3-carboxylate, respectively. The ester functionalities connected to methylene group were regiospecifically converted to the desired monoacyl azides. Curtius rearrangement of acyl azides followed by hydrolysis of the formed isocyanates gave the symmetrical urea derivatives. Cyclization of the ester groups provided the target compounds.     

Tuning the Bergman cyclization by introduction of metal fragments at various positions of the enediyne. Metalla-Bergman cyclizations

We expand the scope of the Bergman cyclization by exploring computationally the rearrangement of two osmaenediynes and one rhodaenediyne. The three hypothetical metallaenediynes are constructed by substituting the 14-electron Os(PH3)3 fragment for the C fragment, or the 15-electron Os(PH3)3H or Rh(PH3)3 fragments for the sp2 CH fragment, of 3-ene-1,5-diyne. This replacement is guided by the isolobal analogy and previous metallabenzene chemistry. The rearrangement of osmaenediyne with an Os(PH3)3 fragment in place of C is exothermic by 3 kcal/mol (the parent Bergman reaction is computed to be endothermic by 5 kcal/mol) and associated with a significant decrease in the barrier to rearrangement to 13 kcal/mol (the Ea of the parent reaction computed at the same level of theory is 33 kcal/mol). The replacement of a CH by the isolobal analogue Os(PH3)3H reduces the energy of activation for the rearrangement to 23 kcal/mol and produces a corresponding metalladiradical that is 8 kcal/mol less stable that the corresponding osmaenediyne. The activation energy corresponding to the rearrangement of the rhodaenediyne is the same as that of the organic parent enediyne. Interesting polytopal rearrangements of metallaenediynes and the diradical nature of the resulting intermediates are also explored.     

Regioselective synthesis of the 3,4-dihydrofuro[3,2-d]pyrimidin-2(1H)-one skeleton: a new class of compound

We hereby report the first preparation of the 3,4-dihydrofuro[3,2-d]pyrimidin-2(1H)-one skeleton formed by two controlled Curtius rearrangements of the corresponding acyl azides, prepared from 2-(2-methoxy-2-oxoethyl)furan-3-carboxylate via the hydrazide. Rearrangement of the acyl azides followed by trapping by nucleophiles and intramolecular trapping provided the target compounds.     

Thermolysis of alkyl sulfoxides and derivatives: a comparison of experiment and theory.

Gas-phase activation data were obtained for model sulfoxide elimination reactions. The activation enthalpy for methyl 3-phenylpropyl sulfoxide is 32.9 +/- 0.9 kcal/mol. Elimination by methyl vinyl sulfoxide to form acetylene has an enthalpic barrier of 41.6 +/- 0.8 kcal/mol and that of 3-phenylpropyl methanesulfinate to form hydrocinnamaldehyde is 34.6 +/- 0.6 kcal/mol. Calculations at the MP2/6-311+G(3df,2p)//MP2/6-31G(d,p) level for simplified models of these reactions provide barriers of 32.3, 40.3, and 32.7 kcal/mol, respectively. A series of other compounds are examined computationally, and it is shown that the substituent effects on the sulfoxide elimination reaction are much more straightforward to interpret if DeltaH data are available in addition to the usually determined DeltaH++. The activation enthalpy of the reverse addition reaction is also subject to structural variation and can usually be rationalized on the basis of nucleophilicity of the sulfur or polarity matching between the sulfenic acid and olefin derivative.     

McLafferty rearrangement of the radical cations of butanal and 3-fluorobutanal: a theoretical investigation of the concerted and stepwise mechanisms

The stepwise and concerted pathways for the McLafferty rearrangement of the radical cations of butanal (Bu(+)) and 3-fluorobutanal (3F-Bu(+)) are investigated with density functional theory (DFT) and ab initio methods in conjunction with the 6-311+G(d,p) basis set. A concerted transition structure (TS) for Bu(+), (H), is located with a Gibbs barrier height of 37.7 kcal/mol as computed with CCSD(T)//BHandHLYP. Three pathways for the stepwise rearrangement of Bu(+) have been located, which are all found to involve different complexes. The barrier height for the H(gamma) transfer is found to be 2.2 kcal/mol, while the two most favorable TSs for the C(alpha)-C(beta) cleavage are located 8.9 and 9.2 kcal/mol higher. The energies of the 3F-Bu(+) system have been calculated with the promising hybrid meta-GGA MPWKCIS1K functional of DFT. Interestingly, the fluorine substitution yields a barrier height of only 20.5 kcal/mol for the concerted TS, (3F-H). A smaller computed dipole moment, 12.1 D, for (3F-H) compared with 103.2 D for (H) might explain the stabilization of the substituted TS. The H(gamma) transfer, with a barrier height of 4.9 kcal/mol, is found to be rate-determining for the stepwise McLafferty rearrangement of 3F-Bu(+), in contrast to the unsubstituted case. By inspection of the spin and charge distributions of the stationary points, it is noted that the bond cleavages in the concerted rearrangements are mainly of heterolytic nature, while those in the stepwise channels are found to be homolytic.     

DFT studies on the acid-catalyzed Curtius reaction: the Schmidt reaction

The Schmidt reaction is the acid-catalyzed analogue of the Curtius reaction and is extensively used in organic synthesis. In this work, the mechanism of this reaction has been explored using DFT calculations at the B3LYP/6-311+G(d,p) level. Protonated formyl azide may undergo rearrangement to the product, protonated isocyanic acid, with simultaneous extrusion of molecular nitrogen (concerted mechanism), or undergo rearrangement to the anti conformer, followed by removal of nitrogen to form the nitrenium ion, which then rearranges to the final product, protonated isocyanic acid (step-wise mechanism). Like the Curtius reaction, it is found that the concerted pathway is definitely preferred. The key role of acidification in decreasing the overall energy barrier is more highlighted in case of phenyl substitution, with negligible effect on the lower homologues. For methoxy and amine substituents, there is very little difference in the activation energies of the concerted and step-wise reactions, with the former being still slightly preferred. Unlike the parent compound, the rearrangement of substituted nitrenium ion in some cases involves side reactions like C-H insertion and cyclization.     

Structure and photochemical behaviour of 3-azido-acrylophenones: a matrix isolation infrared spectroscopy study

(Z)-3-Azido-3-methoxycarbonyl-2-chloro-acrylophenone (MACBP) has been synthesized, isolated in low temperature argon and xenon matrices and studied by FTIR spectroscopy, complemented by DFT(B3LYP)/6-311++G(d,p) calculations. The molecule was characterized both structurally and spectroscopically, and its photochemistry used to probe the mechanism of photo-induced conversion of 3-azido-acrylophenones into oxazoles. In situ UV irradiation (λ = 235 nm) of matrix-isolated MACBP yielded as primary photoproduct a 2H-azirine, which undergoes subsequent photoisomerization to methyl 4-chloro-5-phenyl-1,3-oxazole-2-carboxylate. In a competitive process, a ketenimine is also formed upon photolysis of MACBP. The reported results indicate that this ketenimine must be formed from the starting 3-azido-acrylophenone via a Curtius type concerted rearrangement.     

Thermal isomerization of acetylnitrene: a quantum-chemical study

The electronic structure and pathways of thermal isomerization of formylnitrene and acetylnitrene were studied by the B3LYP/6-311G(d,p) density functional method and ab initio G2(MP2,SVP) computational procedure using the geometries obtained from B3LYP calculations. According to G2 calculations, both nitrenes have singlet ground states while the energies of the corresponding triplet states are 2.8 and 5.7 kcal mol^–1 higher. For acetylnitrene, the activation barrier to the nitrene isocyanate isomerization was estimated at 28.9 kcal mol^–1 (G2). Calculations revealed no pathway for single-step isomerization of nitrene into cyanate in both systems. The formation of methyl cyanate from isocyanate is thermodynamically unfavorable (E = 26.5 kcal mol^–1) and requires a high activation barrier (89.4 kcal mol^–1) should be overcome. Based on the results obtained, the pathways of transformation of nitrene formed in thermal decomposition of acetyl azide (Curtius rearrangement) were analyzed.     

Mechanisms of the Knoevenagel hetero Diels-Alder sequence in multicomponent reactions to dihydropyrans: experimental and theoretical investigations into the role of water

The role of water in a multicomponent domino reaction (MCR) involving styrene, 2,4-pentanedione, and formaldehyde was studied. Whereas anhydrous conditions produced no reaction, the MCR successfully proceeded in the presence of water, affording the targeted dihydropyran derivatives with good yield. The mechanism of this MCR (Knoevenagel hetero Diels-Alder sequence) was studied with and without explicit water molecules using the SMD continuum solvation model in combination with the B3LYP density functional and the 6-311++G** basis set to compute the water and acetone (aprotic organic solvent) solution Gibbs free energies. In the Knoevenagel step, we found that water acted as a proton relay to favor the formation of more flexible six-membered ring transition state structures both in concerted (direct H(2)O elimination) and stepwise (keto-enol tautomerization and dehydration) pathways. The inclusion of a water molecule in our model resulted in a significant decrease (-8.5 kcal mol(-1)ΔG(water)(?)) of the direct water elimination activation barrier. Owing to the presence of water, all chemical steps involved in the MCR mechanism had activation free energies barriers lower than 39 kcal mol(-1) at 25 °C in aqueous solvent (<21 kcal mol(-1) ZPE corrected electronic energies barriers). Consequently, the MCR proceeded without the assistance of any catalyst.     

Theoretical study on the reaction pathways of HFCO + H2O

For the reaction of methanoyl fluoride with water, both optimized structures and vibrational wavenumbers of reaction intermediates, transition structures and product complexes were calculated and characterized with theory at the MP2/6-311++G(d,p) level. Including a catalytic path and concerted and stepwise hydrolysis paths, possible reaction mechanisms were also investigated. The catalytic reaction of HFCO yielding HF and CO has the smallest activation barrier, 29.6 kcal/mol, whereas for the concerted hydrolysis 33.0 kcal/mol is required to overcome the barrier to form transoid HCOOH + HF, which is less than for the stepwise counterpart, 42.0 kcal/mol.     

Quantum mechanical study of the potential energy surface of the ClO + NO2 reaction

In the study of the reaction pathways of the ClO + NO2 reaction including reliable structures of the reactants, products, intermediates, and transition states as well as energies the MP2/6-311G(d), B3LYP/6-311G(d), and G2(MP2) methods have been employed. Chlorine nitrate, ClONO2, is formed by N-O association without an entrance barrier and is stabilized by 29.8 kcal mol(-1). It can undergo either a direct 1,3 migration of Cl or OCl rotation to yield an indistinguishable isomer. The corresponding barriers are 45.8 and 7.1 kcal mol(-1), respectively. ClONO2 can further decompose into NO3 + Cl with an endothermicity of 46.4 kcal mol(-1). The overall endothermicity of the NO2 + ClO --> NO3 + Cl reaction is calculated to be 16.6 kcal mol(-1). The formation of cis-perp and trans-perp conformer of chlorine preoxynitrite, ClOONO(cp) and ClOONO(tp), are exothermic by 5.4 and 3.8 kcal mol(-1), respectively. Calculations on the possible reaction pathways for the isomerization of ClOONO to ClONO2 showed that the activation barriers are too high to account for appreciable nitrate formation from peroxynitrite isomerization. All quoted relative energies are related to G2(MP2) calculations.     

Role of the somersault rearrangement in the oxidation step for flavin monooxygenases (FMO). A comparison between FMO and conventional xenobiotic oxidation with hydroperoxides

Model quantum mechanical calculations presented for C-4a-flavin hydroperoxide (FlHOOH) at the B3LYP/6-311+G(d,p) level suggest a new mechanism for flavoprotein monooxygenase (FMO) oxidation involving a concerted homolytic O-O bond cleavage in concert with hydroxyl radical transfer from the flavin hydroperoxide rather than an S(N)2-like displacement by the substrate on the C-4a-hydroperoxide OOH group. Homolytic O-O bond cleavage in a somersault-like rearrangement of hydroperoxide C-4a-flavinhydroperoxide (1) (FLHO-OH → FLHO···HO) produces an internally hydrogen-bonded HO(?) radical intermediate with a classical activation barrier of 27.0 kcal/mol. Model hydroperoxide 1 is used to describe the transition state for the key oxidation step in the paradigm aromatic hydroxylase, p-hydroxybenzoate hydroxylase (PHBH). A comparison of the electron distribution in the transition structures for the PHBH hydroxylation of p-hydroxybenzoic acid (ΔE(?) = 23.0 kcal/mol) with that of oxidation of trimethylamine (ΔE(?) = 22.3 kcal/mol) and dimethyl sulfide (ΔE? = 14.1 kcal/mol) also suggests a mechanism involving a somersault mechanism in concert with transfer of an HO(?) radical to the nucleophilic heteroatom center with a hydrogen transfer back to the FLH-O residue after the barrier is crossed to produce the final product, FLH-OH. In each case the hydroxylation barrier was less than that of the O-O rearrangement barrier in the absence of a substrate supporting an overall concerted process. All three transition structures bear a resemblance to the TS for the comparable hydroxylation of isobutane (ΔE(?) = 29.2 kcal/mol) and for simple Fenton oxidation by aqueous iron(III) hydroperoxides. To our surprise the oxidation of N- and S-nucleophiles with conventional oxidants such as alkyl hydroperoxides and peracids also proceeds by HO(?) radical transfer in a manner quite similar to that for tricyclic hydroperoxide 1. Stabilization of the developing oxyradical produced by somersault rearrangement for concerted enzymatic oxidation with tricyclic hydroperoxide 1 results in a reduced overall activation barrier.     

Quantum mechanical study of the ring-closing reaction of the hexatriene radical cation

The ring-closing reaction of hexatriene radical cation 1(*)(+) to 1,3-cyclohexadiene radical cation 2(*)(+) was studied computationally at the B3LYP/6-31G* and QCISD(T)/6-311G*//QCISD/6-31G* levels of theory. Both, concerted and stepwise mechanisms were initially considered for this reaction. Upon evaluation at the B3LYP level of theory, three of the possible pathways-a concerted C(2)-symmetric via transition structure 3(*)(+) and stepwise C(1)-symmetric pathways involving three-membered ring intermediate 5(*)(+) and four-membered ring intermediate 6(*)(+)-were rejected due to high-energy stationary points along the reaction pathway. The two remaining pathways were found to be of competing energy. The first proceeds through the asymmetric, concerted transition structure 4(*)(+) with an activation barrier E(a) = 16.2 kcal/mol and an overall exothermicity of -23.8 kcal/mol. The second pathway, beginning from the cis,cis,trans rotamer of 1(*)(+), proceeds by a stepwise pathway to the cyclohexadiene product with an overall exothermicity of -18.6 kcal/mol. The activation energy for the rate-determining step in this process, the formation of the intermediate bicyclo[3.1.0]hex-2-ene via transition structure 9(*)(+), was found to be 20.4 kcal/mol. More rigorous calculations of a smaller subsection of the potential energy hypersurface at the QCISD(T)//QCISD level confirmed these findings and emphasized the importance of conformational control of the reactant.