Energy profiles for ketene cyclizations. interconversion of 1,3-oxazin-6-ones, mesoionic 1,3-oxazinium olates and acylketenes, imidoylketenes, oxoketenimines, and cyclization products

The energy surface connecting oxazinium olates 9, several possible conformers of ketenes 10 and 11, and the final cyclization products 12, 13 and 14, as well as the isomeric 1,3-oxazine-6-ones 15, ring opening of the latter to N-acylimidoylketenes 16, and subsequent rearrangement of 16 to oxoketenimines 17, azetinones 18, and the cyclization products 19 and 20 are evaluated computationally at the B3LYP/6-31G* and B3LYP/6-311+G*//B3LYP/6-31G* levels. The cyclizations of ketenes to oxazinium olates 9 and oxazines 15 have the characteristics of pseudopericyclic reactions. Plots of the energy vs internal reaction coordinate for the cyclization of transoid acylketenes such as 10 to 9 (via TS1) and 16 to 15 (via TS7) feature two inflection points and indicate that the part of the energy surface above the lower inflection points describe internal rotation of the acyl function in the ketene moiety, and the part below this point describes the cyclization of the cisoid ketene to the planar mesoionic oxazinium olate 9 or oxazinone 15. The 1,3-shifts of the OR group that interconvert ketenes 16 and ketenimines 17 via four-membered cyclic transition states TS8 behave similarly, the first portion (from the ketenimine side) of the activation barrier being due largely to internal rotation of substituents, and the top part being due to the 1,3-shift proper.

Leading coordinate analysis of reaction pathways in proton chain transfer: Application to a two-proton transfer model for the green fluorescent protein

The ‘leading coordinate’ approach to computing an approximate reaction pathway, with subsequent determination of the true minimum energy profile, is applied to a two-proton chain transfer model based on the chromophore and its surrounding moieties within the green fluorescent protein (GFP). Using an ab initio quantum chemical method, a number of different relaxed energy profiles are found for several plausible guesses at leading coordinates. The results obtained for different trial leading coordinates are rationalized through the calculation of a two-dimensional relaxed potential energy surface (PES) for the system. Analysis of the 2-D relaxed PES reveals that two of the trial pathways are entirely spurious, while two others contain useful information and can be used to furnish starting points for successful saddle-point searches. Implications for selection of trial leading coordinates in this class of proton chain transfer reactions are discussed, and a simple diagnostic function is proposed for revealing whether or not a relaxed pathway based on a trial leading coordinate is likely to furnish useful information.

New insights into NO-carbon and N2O-carbon reactions from quantum mechanical calculations

Density functional theory calculations were used to investigate the mechanisms of NO-carbon and N2O-carbon reactions. It was the first time that the importance of surface nitrogen groups was addressed in the kinetic behaviors of the NO-carbon reaction. It was found that the off-plane nitrogen groups that are adjacent to the zigzag edge sites and in-plane nitrogen groups that are located on the armchair sites make the bond energy of oxygen desorption even ca. 20% lower than that of the off-plane epoxy group adjacent to zigzag edge sites and in-plane o-quinone oxygen atoms on armchair sites; this may explain the reason why the experimentally obtained activation energy of the NO-carbon reaction is ca. 20% lower than that of the O-2-carbon reaction over 923 K. A higher ratio of oxygen atoms can be formed in the N2O-carbon reaction, because of the lower dissociation energy of N2O, which results in a higher ratio of off-plane epoxy oxygen atoms. The desorption energy of semiquinone with double adjacent off-plane oxygen groups is ca. 20% less than that of semiquinone with only one adjacent off-plane oxygen group. This may be the reason why the activation energy of N2O is also ca. 20% less than that of the O-2-carbon reaction. The new mechanism can also provide a good qualitative comparison for the relative reaction rates of NO-, N2O-, and O-2-carbon reactions. The anisotropic characters of these gas-carbon reactions can also be well explained.

Electron transfer in acetohydroxy acid synthase as a side reaction of catalysis. implications for the reactivity and partitioning of the carbanion/enamine form of (alpha-hydroxyethyl)thiamin diphosphate in a "nonredox" flavoenzyme

Acetohydroxy acid synthases (AHAS) are thiamin diphosphate- (ThDP-) and FAD-dependent enzymes that catalyze the first common step of branched-chain amino acid biosynthesis in plants, bacteria, and fungi. Although the flavin cofactor is not chemically involved in the physiological reaction of AHAS, it has been shown to be essential for the structural integrity and activity of the enzyme. Here, we report that the enzyme-bound FAD in AHAS is reduced in the course of catalysis in a side reaction. The reduction of the enzyme-bound flavin during turnover of different substrates under aerobic and anaerobic conditions was characterized by stopped-flow kinetics using the intrinsic FAD absorbance. Reduction of enzyme-bound FAD proceeds with a net rate constant of k' = 0.2 s(-1) in the presence of oxygen and approximately 1 s(-1) under anaerobic conditions. No transient flavin radicals are detectable during the reduction process while time-resolved absorbance spectra are recorded. Reconstitution of the binary enzyme-FAD complex with the chemically synthesized intermediate 2-(hydroxyethyl)-ThDP also results in a reduction of the flavin. These data provide evidence for the first time that the key catalytic intermediate 2-(hydroxyethyl)ThDP in the carbanionic/enamine form is not only subject to covalent addition of 2-keto acids and an oxygenase side reaction but also transfers electrons to the adjacent FAD in an intramolecular redox reaction yielding 2-acetyl-ThDP and reduced FAD. The detection of the electron transfer supports the idea of a common ancestor of acetohydroxy acid synthase and pyruvate oxidase, a homologous ThDP- and FAD-dependent enzyme that, in contrast to AHASs, catalyzes a reaction that relies on intercofactor electron transfer.

Converged quantum calculations of HO2 bound states and resonances for J=6 and 10

Bound and resonance states of HO2 are calculated quantum mechanically using both the Lanczos homogeneous filter diagonalization method and the real Chebyshev filter diagonalization method for nonzero total angular momentum J=6 and 10, using a parallel computing strategy. For bound states, agreement between the two methods is quite satisfactory; for resonances, while the energies are in good agreement, the widths are in general agreement. The quantum nonzero-J specific unimolecular dissociation rates for HO2 are also calculated. (C) 2004 American Institute of Physics.

Facile ketene-ketene and ketene-ketenimine rearrangements: A study of the 1,3-migration of alpha-substituents interconverting alpha-imidoylketenes and alpha-oxoketenimines, a pseudopericyclic reaction

Theoretical calculations (B3LYP/6-311+G(3df,2p)//B3LYP/6-31G*) of the 1,3 migration of NR2 transforming alpha-oxoketenimines 1 to alpha-imidoylketenes 3 and vice versa indicate that this process is a pseudo-pericyclic reaction with a low activation energy (NH2 97 kJ mol(-1), N(CH3)(2) 62 kJ mol(-1)). The oxoketenimines were found to be more stable (by 18-35 kJ mol(-1)) which is in line with experimental observations. The hindered amine rotation in the amide and amidine moieties adjacent to the cumulenes are important in the migration of the NR2 group, as one of the rotation transition states is close to the 1,3 migration pathway. This gives an interesting potential energy surface with a valley-ridge inflection (VRI) between the orthogonal hindered amine rotation and 1,3 migration transition states. The imidoylketene may also undergo ring closure to an azetinone 5; however, this is metastable, and under the conditions that allow the 1,3-migration, the oxoketenimine 1 will be favored. The imine NH E/Z-interconversion of the ketenimine group takes place by inversion and has a low activation barrier (similar to40 kJ mol(-1)). In all the amidines examined the E/Z-interconversion of the imine function was predicted to be by rotation with a high barrier (>80 kJ mol(-1)), in contrast to all other reported imine E/Z-interconversions which are by inversion.

Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Unimolecular Rovibrational Bound and Resonance States for Large Angular Momentum: J=20 Calculations for HO2

We explore the calculation of unimolecular bound states and resonances for deep-well species at large angular momentum using a Chebychev filter diagonalization scheme incorporating doubling of the autocorrelation function as presented recently by Neumaier and Mandelshtam [Phys. Rev. Lett. 86, 5031 (2001)]. The method has been employed to compute the challenging J=20 bound and resonance states for the HO2 system. The methodology has firstly been tested for J=2 in comparison with previous calculations, and then extended to J=20 using a parallel computing strategy. The quantum J-specific unimolecular dissociation rates for HO2-> H+O-2 in the energy range from 2.114 to 2.596 eV have been reported for the first time, and comparisons with the results of Troe and co-workers [J. Chem. Phys. 113, 11019 (2000) Phys. Chem. Chem. Phys. 2, 631 (2000)] from statistical adiabatic channel method/classical trajectory calculations have been made. For most of the energies, the reported statistical adiabatic channel method/classical trajectory rate constants agree well with the average of the fluctuating quantum-mechanical rates. Near the dissociation threshold, quantum rates fluctuate more severely, but their average is still in agreement with the statistical adiabatic channel method/classical trajectory results.