Computational study of the aminolysis of 2-benzoxazolinone

Three possible mechanisms (zwitterionic, neutral stepwise, and neutral concerted) of the ring-opening reaction of 2-benzoxazolinone (BO) upon aminolysis with methylamine were studied at the B3LYP/6-31G* level. In the gas phase, the neutral concerted mechanism is shown to be most favorable, which proceeds via a rate-determining barrier of 28-29 kcal/mol. The transition state, CTS, associated with this barrier is a four-centered one, where 1,2-addition of the N[bond]H of methylamine to the C[bond]O of BO ring occurs. The rate-determining barrier of the neutral stepwise pathway is found to be ca. 42 kcal/mol. The inclusion of solvent effects by a polarizable continuum model (PCM) does not change the conclusions based on the gas-phase study; the barrier at CTS is reduced to 20, 20, and 22 kcal/mol in water, ethanol, and acetonitrile, respectively.     

Aromatic dienophiles. 1. A theoretical study of an inverse-electron demand Diels-Alder reaction between 2-aminopyrrole and 1,3,5-triazine

This study is devoted to a detailed theoretical study of an inverse-electron demand Diels-Alder reaction (IDA) with 1,3,5-triazine as the diene and 2-aminopyrrole 1A(alpha) as the dienophile, which is a key step in a cascade reaction for the one-pot synthesis of purine analogues. Geometries were optimized with the B3LYP/6-31G* method and energies were evaluated with the MP2/6-311++G** method. This IDA reaction occurs through a stepwise mechanism, where the first step corresponds to the nucleophilic attack of 2-aminopyrrole to triazine to form a zwitterionic intermediate, which is in equilibrium with a neutral intermediate through a hydrogen transfer process, followed by a rate-determining ring-closure step. It is shown that the B3LYP method significantly overestimates the activation energy, whereas the MP2 method offers a reasonable activation barrier of 27.9 kcal/mol in the gas phase. The solvation effect has been studied by the PCM model. In DMSO, the calculated activation energy of the IDA reaction is decreased to 24.0 kcal/mol with a strong endothermicity of 17.4 kcal/mol due to the energy penalty of transforming two aromatic reactants into a nonaromatic IDA adduct. The possible stepwise [2+2] pathway is ruled out based on its higher activation and reaction energies than those of the [4+2] pathway. By comparing the IDA reactions of triazine to 2-aminopyrrole and pyrrole, we address two crucial roles of the alpha-amino substituent in lowering activation and reaction energies and controlling the reaction regiochemistry.     

Reactions of 1,3-cyclohexadiene with singlet oxygen. A theoretical study.

A thorough study of the reaction of singlet oxygen with 1,3-cyclohexadiene has been made at the B3LYP/6-31G(d) and CASPT2(12e,10o) levels. The initial addition reaction follows a stepwise diradical pathway to form cyclohexadiene endoperoxide with an activation barrier of 6.5 kcal/mol (standard level = CASPT2(12e,10o)/6-31G(d); geometries and zero-point corrections at B3LYP/6-31G(d)), which is consistent with an experimental value of 5.5 kcal/mol. However, as the enthalpy of the transition structure for the second step is lower than the diradical intermediate, the reaction might also be viewed as a nonsynchronous concerted reaction. In fact, the concertedness of the reaction is temperature dependent since entropy differences create a free energy barrier for the second step of 1.8 kcal/mol at 298 K. There are two ene reactions; one is a concerted mechanism (DeltaH(double dagger) = 8.8 kcal/mol) to 1-hydroperoxy-2,5-cyclohexadiene (5), while the other, which forms 1-hydroperoxy-2,4-cyclohexadiene (18), passes through the same diradical intermediate (9) as found on the pathway to endoperoxide. The major pathway from the endoperoxide is O-O bond cleavage (22.0 kcal/mol barrier) to form a 1,4-diradical (25), which is 13.9 kcal/mol less stable than the endoperoxide. From the diradical, two low-energy pathways exist, one to epoxyketone (29) and the other to the diepoxide (27), where both products are known to be formed experimentally with a product ratio sensitive to the nature of substitutents. A significantly higher activation barrier leads to C-C bond cleavage and direct formation of maleic aldehyde plus ethylene.     

Disproportionation of Gold(II) complexes. A density functional study of ligand and solvent effects

A computational study of gold(II) disproportionation is presented for the atomic ion as well as complexes with chloride and neutral ligands. The Au2+ atomic ion is stable to disproportionation, but the barrier is more than halved to 119 kcal/mol in an aqueous environment vs 283 kcal/mol in the gas phase. For dissociative disproportionation of chloride complexes, the loss of chlorine, either as an atom (Delta G(aq) = +20 kcal/mol) or as an anion (Delta G(aq) = +15 kcal/mol) represents the largest calculated barrier. The calculated transition state for associative disproportionation is only 9 kcal/mol above separated Au(II)Cl3(-) anions. For the disproportionation of Au(II)L3 complexes with neutral ligands, disproportionation is highly endergonic in the gas phase. Calculations imply that for synthesis of a monometallic Au(II) complex, a nonpolar solvent is preferred. With the exception of [Au(CO)3]2+, disproportionation of Au(II)L3 complexes to Au(I)L and Au(III)L3 is exergonic in solution phase for the ligands investigated. The driving force is provided by the very favorable solvation free energy of the trivalent gold complex. The solvation free energy contribution to the reaction (Delta G(solv)) is very large for small and polar ligands such as ammonia and water. Furthermore, calculations imply that choosing ligands that would yield neutral species upon disproportionation may provide an effective route to thwart this decomposition pathway for Au(II) complexes. Likewise, bulkier ligands that yield larger, more weakly solvated complex ions would appear to be desirable.     

Degenerate electron exchange reaction of n-alkane radical cations in solution

The degenerate electron exchange (DEE) reaction involving radical cations (RCs) of n-nonane, n-dodecane, and n-hexadecane in n-hexane solution was studied over the temperature range 253-313 K using the method of time-resolved magnetic field effect in recombination fluorescence of spin-correlated radical ion pairs. In the dilute solutions the rate constant of DEE was found to be 200 times slower than the diffusion limit. Using n-nonane as an example, we showed that two reasons are responsible for the low value of the RC self-exchange rate: (1) conformational variability of molecules and RCs and (2) the activation barrier of DEE reaction. The calculations of the reaction enthalpy performed by the B3LYP/6-31G(d) method indicated that electron transfer can be effective only upon collision of RC with a neutral molecule either in the all-trans conformation or in the conformation differing from the latter by rotation of the end ethyl fragment. The activation barrier of the DEE reaction was estimated using the reorganization energy of the internal degrees of freedom calculated at the B3LYP level and was found to be about 6 kcal/mol. A possible influence of the interaction between RC and a neutral molecule in an encounter complex on DEE rate constant is also discussed.     

Fixation of heterocumulenes. 1. A theoretical study on the irreversible reaction of CO(2) with a 2-lithio-1,3-dithiane

The mechanism of electrophilic addition of CO(2) to 2-lithio-2-phenyl-1, 3-dithiane-tetrahydrofuran-tetramethylethylendiamine, a Corey-Seebach "umpolungs" reagent for nucleophilic acylation, was investigated at the B3YLP/6-311+G//HF/6-31+G level of theory with specific solvation effects being included in the study. The overall reaction is exothermic by 18.8 kcal/mol and proceeds via precomplexation of the CO(2) (necessary for activation). The reaction barrier is calculated to be 11.9 kcal/mol (relative to the precomplex 2a) and represents a lower boundary for the activation energy. The reaction barrier is shown to originate from specific solvent effects. A detailed survey of the electronic effects governing the reaction is given. The ab initio results were then compared with semiempirical PM3 calculations, which were extended to include the heterocumulenes COS and CS(2). As expected, the reduced electrophilicity of these molecules leads to a higher activation barrier by retention of the same reaction mechanism.     

Reactivity of alkyl versus silyl peroxides. The consequences of 1, 2-silicon bridging on the epoxidation of alkenes with silyl hydroperoxides and bis(trialkylsilyl)peroxides

The bond dissociation energies for a series of silyl peroxides have been calculated at the G2 and CBS-Q levels of theory. A comparison is made with the O-O BDE of the corresponding dialkyl peroxides, and the effect of the O-O bond strength on the activation barrier for oxygen atom transfer is discussed. The O-O bond dissociation enthalpies (DeltaH(298)) for bis (trimethylsilyl) peroxide (1) and trimethylsilyl hydroperoxide (2) are 54.8 and 53.1 kcal/mol, respectively at the G2 (MP2) and CBS-Q levels of theory. The O-O bond dissociation energies computed at G2 and G2(MP2) levels for bis(tert-butyl) peroxide and tert-butyl hydroperoxide are 45.2 and 48.3 kcal/mol, respectively. The barrier height for 1,2-methyl migration from silicon to oxygen in trimethylsilyl hydroperoxide is 47.9 kcal/mol (MP4//MP2/6-31G). The activation energy for the oxidation of trimethylamine to its N-oxide by bis(trimethylsilyl) peroxide is 28.2 kcal/mol (B3LYP/6-311+G(3df,2p)// B3LYP/6-31G(d)). 1,2-Silicon bridging in the transition state for oxygen atom transfer to a nucleophilic amine results in a significant reduction in the barrier height. The barrier for the epoxidation of E-2-butene with bis(dimethyl(trifluoromethyl))silyl peroxide is 25.8 kcal/mol; a reduction of 7.5 kcal/mol relative to epoxidation with 1. The activation energy calculated for the epoxidation of E-2-butene with F(3)SiOOSiF(3) is reduced to only 2.2 kcal/mol reflecting the inductive effect of the electronegative fluorine atoms.     

Investigation of ethylene dimerization over faujasite zeolite by the ONIOM method.

Ethylene dimerization was investigated by using an 84T cluster of faujasite zeolite modeled by the ONIOM3(MP2/6-311++G(d,p):HF/6-31G(d):UFF) method. Concerted and stepwise mechanisms were evaluated. In the stepwise mechanism, the reaction proceeds by protonation of ethylene to form the surface ethoxide and then C--C bond formation between the ethoxide and the second ethylene molecule to give the butoxide product. The first step is rate-determining and has an activation barrier of 30.06 kcal mol(-1). The ethoxide intermediate is rather reactive and readily reacts with another ethylene molecule with a smaller activation energy of 28.87 kcal mol(-1). In the concerted mechanism, the reaction occurs in one step of simultaneous protonation and C--C bond formation. The activation barrier is calculated to be 38.08 kcal mol(-1). Therefore, the stepwise mechanism should dominate in ethylene dimerization.     

Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction

To investigate fundamental features of enzyme catalysis, there is a need for high-level calculations capable of modelling crucial, unstable species such as transition states as they are formed within enzymes. We have modelled an important model enzyme reaction, the Claisen rearrangement of chorismate to prephenate in chorismate mutase, by combined ab initio quantum mechanics/molecular mechanics (QM/MM) methods. The best estimates of the potential energy barrier in the enzyme are 7.4-11.0 kcal mol(-1)(MP2/6-31+G(d)//6-31G(d)/CHARMM22) and 12.7-16.1 kcal mol(-1)(B3LYP/6-311+G(2d,p)//6-31G(d)/CHARMM22), comparable to the experimental estimate of Delta H(++)= 12.7 +/- 0.4 kcal mol(-1). The results provide unequivocal evidence of transition state (TS) stabilization by the enzyme, with contributions from residues Arg90, Arg7, and Arg63. Glu78 stabilizes the prephenate product (relative to substrate), and can also stabilize the TS. Examination of the same pathway in solution (with a variety of continuum models), at the same ab initio levels, allows comparison of the catalyzed and uncatalyzed reactions. Calculated barriers in solution are 28.0 kcal mol(-1)(MP2/6-31+G(d)/PCM) and 24.6 kcal mol(-1)(B3LYP/6-311+G(2d,p)/PCM), comparable to the experimental finding of Delta G(++)= 25.4 kcal mol(-1) and consistent with the experimentally-deduced 10(6)-fold rate acceleration by the enzyme. The substrate is found to be significantly distorted in the enzyme, adopting a structure closer to the transition state, although the degree of compression is less than predicted by lower-level calculations. This apparent substrate strain, or compression, is potentially also catalytically relevant. Solution calculations, however, suggest that the catalytic contribution of this compression may be relatively small. Consideration of the same reaction pathway in solution and in the enzyme, involving reaction from a 'near-attack conformer' of the substrate, indicates that adoption of this conformation is not in itself a major contribution to catalysis. Transition state stabilization (by electrostatic interactions, including hydrogen bonds) is found to be central to catalysis by the enzyme. Several hydrogen bonds are observed to shorten at the TS. The active site is clearly complementary to the transition state for the reaction, stabilizing it more than the substrate, so reducing the barrier to reaction.     

Ab initio and density function theory computational studies of the CH4 + H → CH3 + H2 reaction

The activation barrier for the CH₄ + H → CH₃ + H₂ reaction was evaluated with traditional ab initio and Density Functional Theory (DFT) methods. None of the applied ab initio and DFT methods was able to reproduce the experimental activation barrier of 11.0-12.0 kcal/mol. All ab initio methods (HF, MP2, MP3, MP4, QCISD, QCISD(T), G1, G2, and G2MP2) overestimated the activation energy. The best results were obtained with the G2 and G2MP2 ab initio computational approaches. The zero-point corrected energy was 14.4 kcal mol⁻¹. Some of the exchange DFT methods (HFB) computed energies which were similar to the highly accurate ab initio methods, while the B3LYP hybrid DFT methods underestimated the activation barrier by 3 kcal mol⁻¹. Gradient-corrected DFT methods underestimated the barrier even more. The gradient-corrected DFT method that incorporated the PW91 correlational functional even generated a negative reaction barrier. The suitability of some computational methods for accurately predicting the potential energy surface for this hydrogen radical abstraction reaction was discussed.     

Theoretical study of the mechanism and rate constant of the B + CO2 reaction

The different stationary points on the potential energy surface relative to the title reaction have been reinvestigated at the B3LYP/aug-cc-pVDZ level with relative energies computed at the CCSD(T)/aug-cc-pVTZ level with B3LYP/aug-cc-pVDZ optimized geometries and by using the G3B3 composite method. Two entrance channels have been identified. The first one corresponds to boron addition at one of the oxygen atoms of the CO 2 molecule leading to trans-BOCO, which is found to be about 27 kcal/mol exothermic with a potential energy barrier of 16.4 kcal/mol (G3B3). The second channel, which has not been identified in previous theoretical works, corresponds to a direct insertion of the boron atom into a CO bond and leads to OBCO. The B + CO 2 --> OBCO step is found to be about 84 kcal/mol exothermic and needs to overcome a potential energy barrier of only 3.6 kcal/mol (G3B3). The rate constant at 300 K of the insertion step, calculated by using TST theory with G3B3 calculated activation energy value, is 5.4 10 (-14) cm (3) molecule (-1) s (-1), in very good agreement with the experimental data ((7.0 +/- 2.8) 10 (-14) cm (3) molecule (-1) s (-1), DiGiuseppe, T. G.; Davidovits, P. J. Chem. Phys. 1981, 74, 3287). The one corresponding to the addition process is found to be several orders of magnitude smaller because of a much higher potential energy barrier. The addition channel would not contribute to the title reaction even at high temperature. A modified Arrhenius equation has been fitted in the 300-1000 K temperature range, which might be useful for chemical models.     

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.     

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.     

Multilayer formulation of the fragment molecular orbital method (FMO)

The fragment molecular orbital method (FMO) has been generalized to allow for multilayer structure. Fragments are assigned to layers, and each layer can be described with a different basis set and/or level of electron correlation. Interlayer boundaries are treated in the general spirit of the FMO method since they also coincide with some interfragment boundaries. The question of the one- and two-layer FMO accuracy dependence upon the fragmentation scheme is also addressed. The new method has been applied to predict the reaction barrier and the reaction heat for the Diels-Alder reaction with a representative set of reactants based on dividing fragments in two layers. The 6-31G* basis set has been used for the active site and the 6-31G*, 6-31G, 3-21G, and STO-3G basis sets have been used for the substituents. Different levels of electron correlation (RHF, B3LYP, and MP2) have been applied to layers in systematic fashion. The one-layer FMO errors in the reaction barrier and the reaction heat were 2.0 kcal/mol or less for all levels applied (RHF, B3LYP, and MP2), relative to full ab initio methods. For the two-layer method the error was found to be several kcal/mol. Benchmark calculations of the activation barrier for the decarboxylation of phenylcyanoacetate by beta-cyclodextrin demonstrated that the two-layer calculations are efficient, being 36 times faster than the regular DFT, as well as accurate, with the error being 1.0 kcal/mol.     

Docking, triggering, and biological activity of dynemicin A in DNA: a computational study

The triggering and biological activity of the naturally occurring enediyne dynemicin A (1) was investigated, both inside and outside the minor groove of the duplex 10-mer B-DNA sequence d(CTACTACTGG).d(CCAGTAGTAG), using density functional theory (B3LYP with the 3-21G and 6-31G(d) basis set), BD(T)/cc-pVDZ (Brueckner doubles with a perturbative treatment of triple excitations), and the ONIOM approach. Enediyne 1 is triggered by NADPH in a strongly exothermic reaction (-88 kcal/mol), which involves a number of intermediate steps. Untriggered 1 has a high barrier for the Bergman cyclization (52 kcal/mol) that is lowered after triggering to 16.7 kcal/mol due to an epoxide opening and the accompanying strain relief. The Bergman reaction of triggered 1 is slightly exothermic by 2.8 kcal/mol. The singlet biradical formed in this reaction is kinetically stable (activation enthalpies of 19.5 and 21.8 kcal/mol for retro-Bergman reactions) and is as reactive as para-benzyne. The activity-relevant docking mode is an edge-on insertion into the minor groove, whereas the intercalation between base pairs, although leading to larger binding energies, excludes a triggering of 1 and the development of its biological activity. Therefore, an insertion-intercalation model is developed, which can explain all known experimental observations made for 1. On the basis of the insertion-intercalation model it is explained why large intercalation energies suppress the biological activity of dynemicin and why double-strand scission can be achieved only in a two-step mechanism that involves two enediyne molecules, explaining thus the high ratio of single-strand to double-strand scission observed for 1.     

Relative reactivity of peracids versus dioxiranes (DMDO and TFDO) in the epoxidation of alkenes. A combined experimental and theoretical analysis

Comparative analysis of the calculated gas-phase activation barriers (DeltaE++) for the epoxidation of ethylene with dimethyldioxirane (DMDO) and peroxyformic acid (PFA) [15.2 and 16.4 kcal/mol at QCISD(T)// QCISD/6-31+G(d,p)] and E-2-butene [14.3 and 13.2 kcal/mol at QCISD(T)/6-31G(d)//B3LYP/6-311+G(3df,2p)] suggests similar oxygen atom donor capacities for both oxidants. Competition experiments in CH(2)Cl(2) solvent reveal that DMDO reacts with cyclohexene much faster than peracetic acid/acetic acid under scrupulously dried conditions. The rate of DMDO epoxidation is catalyzed by acetic acid with a reduction in the classical activation barrier of 8 kcal/mol. In many cases, the observed increase in the rate for DMDO epoxidation in solution may be attributed to well-established solvent and hydrogen-bonding effects. This predicted epoxidative reactivity for DMDO is not consistent with what has generally been presumed for a highly strained cyclic peroxide. The strain energy (SE) of DMDO has been reassessed and its moderated value (about 11 kcal/mol) is now more consistent with its inherent gas-phase reactivity toward alkenes in the epoxidation reaction. The unusual thermodynamic stability of DMDO is largely a consequence of the combined geminal dimethyl- and dioxa-substitution effects and unusually strong C-H and C-CH(3) bonds. Methyl(trifluoromethyl)dioxirane (TFDO) exhibits much lower calculated activation barriers than DMDO in the epoxidation reaction (the average DeltaDeltaE++ values are about 7.5 kcal/mol). The rate increase relative to DMDO of approximately 10(5), while consistent with the higher strain energy for TFDO (SE approximately 19 kcal/mol) is attributed largely to the inductive effect of the CF(3) group. We have also examined the effect of alkene strain on the rate of epoxidation with PFA. The epoxidation barriers are only slightly higher for the strained alkenes cyclopropene (DeltaE++ = 14.5 kcal/mol) and cyclobutene (DeltaE++ = 13.7 kcal/mol) than for cyclopentene (DeltaE++ = 12.1 kcal/mol), reflecting the fact there is little relief of strain in the transition state. Alkenes strained by twist or pi-bond torsion do exhibit much lower activation barriers.     

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.     

Understanding the participation of quadricyclane as nucleophile in polar [2sigma + 2sigma + 2pi] cycloadditions toward electrophilic pi molecules

The formal [2sigma + 2sigma + 2pi] cycloaddition of quadricyclane, 1, with dimethyl azodicarboxylate, 2, in water has been studied using DFT methods at the B3LYP/6-31G** and MPWB1K/6-31G** levels. In the gas phase, the reaction of 1 with 2 has a two-stage mechanism with a large polar character and an activation barrier of 23.2 kcal/mol. Inclusion of water through a combined discrete-continuum model changes the mechanism to a two-step model where the first nucleophilic attack of 1 to 2 is the rate-limiting step with an activation barrier of 14.7 kcal/mol. Analysis of the electronic structure of the transition state structures points out the large zwitterionic character of these species. A DFT analysis of the global electrophilicity and nucleophilicity of the reagents provides a sound explanation about the participation of 1 as a nucleophile in these cycloadditions. This behavior is reinforced by a further study of the reaction of 1 with 1,1-dicyanoethylene.     

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.     

A theoretical examination of the acid-catalyzed and noncatalyzed ring-opening reaction of an oxirane by nucleophilic addition of acetate. Implications to epoxide hydrolases

Ab initio and density functional calculations have been performed to gain a better understanding of the epoxide ring-opening reaction catalyzed by epoxide hydrolase. The S(N)2 reaction of acetate with 1S,2S-trans-2-methylstyrene oxide to provide the corresponding diol acetate ester was studied with and without general-acid catalysis. MP2 and DFT (B3LYP) calculations predict, for the noncatalyzed reaction, a central barrier of approximately 20-21 kcal/mol separating the reactants from products depending on which carbon center in the epoxide is undergoing attack. From these gas-phase reactions the immediate alkoxide products are not energetically far below their associated transition states such that the reaction is predicted to be endothermic. Inclusion of aqueous solvation effects via a polarizable continuum model predicts the activation barrier to increase by almost 10 kcal/mol due to the solvation of the acetate ion nucleophile. The activation barrier for the epoxide ring-opening reaction is reduced to approximately 10 kcal/mol when phenol, as the general-acid catalyst, is included in the gas-phase calculations. This is due to the immediate product being the neutral ester rather than the corresponding alkoxide. The transition state in the general-acid-catalyzed reaction is earlier than that for the noncatalyzed reaction and the reaction is highly exothermic. Molecular mechanics calculations of 1S,2S-trans-2-methylstyrene oxide in the active site of murine epoxide hydrolase show two possible binding conformations. Both conformers have the epoxide oxygen forming hydrogen bonds with the acidic hydrogens of the catalytic tyrosines (Tyr381 and Tyr465). These two conformations likely lead to different products since the nucleophile (Asp333-CO(2)(-)) is positioned to react with either carbon center in the epoxide.