Experimental and theoretical study of the OH vibrational spectra and overtone chemistry of gas-phase vinylacetic acid

In this study we present the gas-phase vibrational spectrum of vinylacetic acid with a focus on the nu = 1-5 vibrational states of the OH stretching transitions. Cross sections for nu = 1, 2, 4 and 5 of the OH stretching vibrational transitions are derived on the basis of the vapor pressure data obtained for vinylacetic acid. Ab initio calculations are used to assist in the band assignments of the experimental spectra, and to determine the threshold for the decarboxylation of vinylacetic acid. When compared to the theoretical energy barrier to decarboxylation, it is found that the nu OH = 4 transition with thermal excitation of low frequency modes or rotational motion and nu OH = 5 transitions have sufficient energy for the reaction to proceed following overtone excitation.     

Electronic Factors Influencing the Decarboxylation of beta-Keto Acids. A Model Enzyme Study

A theoretical study of the mechanism of decarboxylation of beta-keto acids is described. A cyclic transition structure was found with essentially complete proton transfer from the carboxylic acid to the beta-carbonyl group. The activation barrier for decarboxylation of formylacetic acid is predicted to be 28.6 kcal/mol (MP4SDTQ/6-31G//MP2/6-31G) while loss of CO(2) from its anion exhibits a barrier of only 20.6 kcal/mol (MP4SDTQ/6-31+G//MP2/6-31+G). Barrier heights of decarboxylation of malonic acid and alpha,alpha-dimethylacetoacetic acid are predicted to be 33.2 and 26.7 kcal/mol, respectively. Model enzyme studies using a thio methyl ester of malonate anion suggests that the role of malonyl-CoA is to afford a polarizable sulfur atom to stabilize the developing enolate anion in the transition structure for decarboxylation. Adjacent positively charged ammonium ions are also observed to stabilize the loss of CO(2) from a carboxylate anion by through-bond Coulombic stabilization of the transition structure.     

Energetics and mechanism of the decomposition of trifluoromethanol

The thermal instability of alpha-fluoroalcohols is generally attributed to a unimolecular 1,2-elimination of HF, but the barrier to intramolecular HF elimination from CF3OH is predicted to be 45.1 +/- 2 kcal/mol. The thermochemical parameters of trifluoromethanol were calculated using coupled-cluster theory (CCSD(T)) extrapolated to the complete basis set limit. High barriers of 42.9, 43.1, and 45.0 kcal/mol were predicted for the unimolecular decompositions of CH2FOH, CHF2OH, and CF3OH, respectively. These barriers are lowered substantially if cyclic H-bonded dimers of CF3OH with complexation energies of approximately 5 kcal/mol are involved. A six-membered ring dimer has an energy barrier of 28.7 kcal/mol and an eight-membered dimer has an energy barrier of 32.9 kcal/mol. Complexes of CF3OH with HF lead to strong H-bonded dimers, trimers and tetramers with complexation energies of approximately 6, 11, and 16 kcal/mol, respectively. The dimer, CH3OH:HF, and the trimers, CF3OH:2HF and (CH3OH)2:HF, have decomposition energy barriers of 26.7, 20.3, and 22.8 kcal/mol, respectively. The tetramer (CH3OH:HF)2 gives rise to elimination of two HF molecules with a barrier of 32.5 kcal/mol. Either CF3OH or HF can act as catalysts for HF-elimination via an H-transfer relay. Because HF is one of the decomposition products, the decomposition reactions become autocatalytic. If the energies due to complexation for the CF3OH-HF adducts are not dissipated, the effective barriers to HF elimination are lowered from approximately 20 to approximately 9 kcal/mol, which reconciles the computational results with the experimentally observed stabilities.     

Self-catalyzed keto-enol tautomerization of malonic acid

We demonstrate through quantum chemical calculations that the keto-enol tautomerization of malonic acid can be catalyzed by the two tautomers of malonic acid itself. This self-catalyzed process proceeds with a relatively low barrier (Gibbs energy ca. 13 kcal/mol in gas phase, 20 kcal/mol in aqueous phase), and involves the concerted transfer of two protons between the substrate and the carboxylic acid functionality of the malonic acid catalyst. This mechanism is expected to compete with the proton relay mechanism currently favored to explain the tautomerization of malonic acid in aqueous media. Malonic acid is an important constituent of secondary organic aerosol where the present chemistry may play a role in determining chemical composition.     

Proton migration and tautomerism in protonated triglycine.

Proton migration in protonated glycylglycylglycine (GGG) has been investigated by using density functional theory at the B3LYP/6-31++G(d,p) level of theory. On the protonated GGG energy hypersurface 19 critical points have been characterized, 11 as minima and 8 as first-order saddle points. Transition state structures for interconversion between eight of these minima are reported, starting from a structure in which there is protonation at the amino nitrogen of the N-terminal glycyl residue following the migration of the proton until there is fragmentation into protonated 2-aminomethyl-5-oxazolone (the b(2) ion) and glycine. Individual free energy barriers are small, ranging from 4.3 to 18.1 kcal mol(-)(1). The most favorable site of protonation on GGG is the carbonyl oxygen of the N-terminal residue. This isomer is stabilized by a hydrogen bond of the type O-H.N with the N-terminal nitrogen atom, resulting in a compact five-membered ring. Another oxygen-protonated isomer with hydrogen bonding of the type O-H.O, resulting in a seven-membered ring, is only 0.1 kcal mol(-)(1) higher in free energy. Protonation on the N-terminal nitrogen atom produces an isomer that is about 1 kcal mol(-)(1) higher in free energy than isomers resulting from protonation on the carbonyl oxygen of the N-terminal residue. The calculated energy barrier to generate the b(2) ion from protonated GGG is 32.5 kcal mol(-)(1) via TS(6-->7). The calculated basicity and proton affinity of GGG from our results are 216.3 and 223.8 kcal mol(-)(1), respectively. These values are 3-4 kcal mol(-)(1) lower than those from previous calculations and are in excellent agreement with recently revised experimental values.     

Density functional calculations of ATP systems. 2. ATP hydrolysis at the active site of actin

The hydrolysis of adenosine 5'-triphosphate (ATP) at the active site of actin has been studied using density functional calculations. The active site is modeled by the triphosphate tail of ATP, an Mg cation, surrounding water molecules, and the nearby protein residues. Four reaction paths have been followed by constraining coordinates that allow phosphate stretching, nucleophilic attack of the catalytic water, and OH(-) formation via water deprotonation. The lowest-energy barrier (21.0 kcal/mol) is obtained for a dissociative reaction where the terminal phosphate breaks on approaching the catalytic water, followed by proton release via a proton wire mechanism. A higher barrier (39.6 kcal/mol) results for an associative reaction path where OH(-) is formed first, with a pentacoordinated phosphorus atom (P-O distances 2.1 A). Stretching the terminal bridging P-O bond results in bond rupture at 2.8 A with an energy barrier of 28.8 kcal/mol. The residues Gln137 and His161 are not important in the reactions, but insight into their roles in vivo has been obtained. The favored coordination of the end products H(2)PO(4)(-) and ADP(3-) includes a hydrogen bond and an O-Mg-O bridge between the phosphates as well as a hydrogen bond between H(2)PO(4)(-) and the Ser14 side chain. The total energy is 2.1 kcal/mol lower than in the initial reactants. Classical simulations of ATP- and ADP.P(i)-actin show few hydrolysis-induced differences in the protein structure, indicating that phosphate migration is necessary for a change in conformation.     

Stability of the hydrogen trioxy radical via infrared action spectroscopy

The hydrogen trioxy radical (HO3) has been proposed as an intermediate in several important chemical reactions and relaxation processes involving OH in the atmosphere. In this work, the gas-phase infrared action spectrum of HO3 is obtained in the OH overtone region, along with the product state distribution of the OH fragment following dissociation. The highest observed OH product channel sets an upper limit for the HO-O2 binding energy of 6.12 kcal mol(-1). The experimental stability of HO3 and derived equilibrium constant imply that up to 66% of atmospheric OH may be converted into HO3 in the tropopause region.     

Second OH overtone excitation and statistical dissociation dynamics of peroxynitrous acid

The second OH overtone transition of the trans-perp conformer of peroxynitrous acid (tp-HOONO) is identified using infrared action spectroscopy. HOONO is produced by the recombination of photolytically generated OH and NO(2) radicals, and then cooled in a pulsed supersonic expansion. The second overtone transition is assigned to tp-HOONO based on its vibrational frequency (10 195.3 cm(-1)) and rotational band contour, which are in accord with theoretical predictions and previous observations of the first overtone transition. The transition dipole moment associated with the overtone transition is rotated considerably from the OH bond axis, as evident from its hybrid band composition, indicating substantial charge redistribution upon OH stretch excitation. The overtone band exhibits homogeneous line broadening that is attributed to intramolecular vibrational redistribution, arising from the coupling of the initially excited OH stretch to other modes that ultimately lead to dissociation. The quantum state distributions of the OH X (2)Pi (nu=0) products following first and second OH overtone excitation of tp-HOONO are found to be statistical by comparison with three commonly used statistical models. The product state distributions are principally determined by the tp-HOONO binding energy of 16.2(1) kcal mol(-1). Only a small fraction of the OH products are produced in nu=1 following the second overtone excitation, consistent with statistical predictions.     

Is the 1,3,5-tridehydrobenzene triradical a cyclopropenyl radical analogue?

The gas-phase heat of formation (DeltaH(f,298)) of the 1,3,5-tridehydrobenzene triradical has been determined by using a negative ion thermochemical cycle. The first three measurements carried out were of the gas-phase acidity of 3,5-dichlorobenzoic acid, the enthalpy for decarboxylation of 3,5-dichlorobenzoate, and the enthalpy for chloride loss from 3,5,-dichlorophenide and constitute the measurement of the heat of formation for 5-chloro-m-benzyne. The last two measurements, the electron affinity of 5-chloro-m-benzyne, and the threshold for chloride loss from 5-chloro-m-benzyne, when combined with DeltaH(f,298) of 5-chloro-m-benzyne, give the heat of formation of the triradical. The 5-chloro-m-benzyne heat of formation is 116.2 +/- 3.7 kcal/mol. The heat of formation of the 1,3,5-tridehydrobenzene triradical measured in this work is 179.1 +/- 4.6 kcal/mol. This heat of formation was used to derive the bond dissociation energy (BDE) at the 5-position of m-benzyne, a third BDE in benzene. The BDE, at 109.2 +/- 5.6 kcal/mol, is ca. 4 kcal/mol lower than the first BDE in benzene (112.9 kcal/mol) and significantly higher than the BDE of phenyl radical at the meta position. The agreement between the first and third BDEs implies that the triradical is best described as a phenyl radical that interacts little with a m-benzyne moiety. The experimentally measured BDE is in good agreement with multireference configuration interaction calculations, which predict a (2)A(1) ground state for the Jahn-Teller distorted triradical. The trends in the first, second, and third BDEs of benzene are similar to those found for cyclopropane, suggesting a cyclopropenyl-like electronic structure within the six-membered ring of the 1,3,5-benzene triradical.     

Gas-phase proton-transfer pathways in protonated histidylglycine

Pathways for proton transfer in the histidylglycine cation are examined in the gas-phase environment with the goal of understanding the mechanism by which protons may become mobile in proteins with basic amino acid residues. An extensive search of the potential energy surface is performed using density functional theory (DFT) methods. After corrections for zero-point energy are included, it is found that all the lowest energy barriers for proton transfer between the N-terminus and the imidazole ring have heights of only a few kcal/mol, while those between the imidazole ring and the backbone amide oxygen have heights of approximately 15 kcal/mol when the proton is moving from the ring to the backbone and only a few kcal/mol when moving from the backbone to the imidazole ring. In mass spectrometric techniques employing collision-induced dissociation to dissociate protein complex ions or to fragment peptides, these barriers can be overcome, and the protons mobilized. Copyright (c) 2008 John Wiley & Sons, Ltd.     

Dissociation dynamics of thiolactic acid at 193 nm: detection of the nascent OH product by laser-induced fluorescence

Electronically excited thiolactic acid (2-mercaptopropionic acid), H(3)C-CH(SH)-COOH, undergoes the C-OH bond cleavage on excitation to the S(2) state at 193 nm, generating the primary product OH (v,J), which is detected by laser-induced fluorescence technique in a collisionless condition of flow system. The partitioning of the available energy between vibrational, rotational, and translational degrees of freedom of nascent photofragments is obtained from relative intensities of ro-vibronic lines in laser-induced fluorescence spectrum of OH, and their Doppler profiles. The rotational population of OH (v(")=0) is characterized by rotational temperature of 408+/-25 K. OH is produced in a vibrationally cold state, i.e., mostly in v(")=0. The average translational energy of OH (v(")=0,J(")) is found to be 21.5+/-2.0 kcal/mol, which implies 25.6 kcal/mol of energy in relative translation of photoproducts corresponding to the f(t) value of approximately 0.6. The observed high translational energy is due to the presence of a barrier in the exit channel, implying that the C-OH bond scission takes place on an electronically excited potential energy surface. The observed partitioning of the available energy between various degrees of the photofragments is theoretically modeled, and the hybrid model, with 26.0 kcal/mol of barrier in the exit channel, is found to explain the measured data quite well. The experimental results are also supported with ab initio molecular orbital calculations for both the ground and the excited electronic states. Time-dependent density functional theory is used to understand the nature of various electronic transitions connecting the lower excited states. Potential energy curves as a function of the C-OH bond length of thiolactic acid suggest distinct exit barriers in the S(1), T(1), and T(2) states. But, we could locate the transition state structure for OH formation in the S(1) state alone. Thus, although thiolactic acid is excited to the S(2) state at 193 nm, it undergoes internal conversion to S(1) where it dissociates to yield OH. In addition to the OH channel from excited electronic states, we studied theoretically all probable dissociation channels occurring on the ground electronic state of thiolactic acid.     

Effects of torsion on O-H stretch overtone spectra and direct overtone photolysis of methyl hydroperoxide

We use laser photoacoustic spectroscopy to obtain overtone spectra at three through six quanta of O-H stretch excitation (3nu(OH)-6nu(OH)) for methyl hydroperoxide (MeOOH). Extending the spectral regions beyond our previous work reveals new features that can be attributed to transitions involving torsion about the O-O bond. Experimental spectral profiles (3nu(OH)-6nu(OH)) and cross sections (3nu(OH)-5nu(OH)) at room temperature show a good agreement with the simulated spectra that we obtain from ab initio calculations employing a vibration-torsion model at 298 K. A Birge-Sponer analysis yields experimental values for the O-H stretch frequency (omega=3773+/-15 cm(-1)) and anharmonicity (omegax=94+/-3 cm(-1)). We also detect OH radicals by laser-induced fluorescence and present photodissociation action spectra of MeOOH in the regions of 4nu(OH) and 5nu(OH). While the spectral profile at 5nu(OH) mimics the photoacoustic spectrum, the peak intensity for transitions to torsionally excited states is relatively more intense in the action spectrum at 4nu(OH), reflecting the fact that the 4nu(OH) excitation energy is below the literature dissociation energy (D0=42.6+/-1 kcal mol(-1)) so that features in the action spectrum come from thermally populated excited states. Finally, we use our calculations to assign contributions to individual peaks in the room-temperature spectra and relate our findings to a recent dynamics study in the literature.     

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

Overtone-induced degradation of perfluorinated alcohols in the atmosphere

Perfluorinated alcohols (PFOHs) are thermally unstable and degrade via loss of HF, ultimately forming perfluorocarboxylic acids. Experiments and calculations of the high activation barrier for the decomposition of CF3OH suggest that the reaction occurs exclusively heterogeneously, perhaps with the involvement of water. Here, we use density functional theory calculations to show that PFOHs may readily complex with water and are expected to be present as a few percent of the total PFOH concentration under ambient atmospheric conditions. The presence of water lowers the HF elimination barrier sufficiently that excitation to low-lying O-H vibrational overtone levels (vOH >or= 3) may cause reaction. Photolysis rate constants for CF3OH x H2O and CF3CF2OH x H2O were estimated to be 6.1 x 10(-8) and 5.6 x 10(-8) s(-1), respectively. PFOH-water complexes should undergo degradation much faster than the corresponding gas-phase unimolecular decomposition of PFOHs, which requires excitation into the vOH = 5 or 6 vibrational level. Overtone-driven gas-phase reactions of PFOH-water complexes could be moderately competitive with heterogeneous reactions with water in dry locations. Overtone-induced degradation of PFOHs is likely a modest atmospheric source of PFCAs to the environment.     

Near infrared photochemistry of pyruvic acid in aqueous solution

Recent experimental and theoretical results have suggested that organic acids such as pyruvic acid, can be photolyzed in the ground electronic state by the excitation of the OH stretch vibrational overtone. These overtones absorb in the near-infrared and visible regions of the spectrum where the solar photons are plentiful and could provide a reaction pathway for the organic acids and alcohols that are abundant in the earth's atmosphere. In this paper the overtone initiated photochemistry of aqueous pyruvic acid is investigated by monitoring the evolution of carbon dioxide. In these experiments CO(2) is being produced by excitation in the near-infrared, between 850 nm and ～1150 nm (11,765-8696 cm(-1)), where the second OH vibrational overtone (Δν = 3) of pyruvic acid is expected to absorb. These findings show not only that the overtone initiated photochemical decarboxylation reaction occurs but also that in the aqueous phase it occurs at a lower energy than was predicted for the overtone initiated reaction of pyruvic acid in the gas phase (13,380 cm(-1)). A quantum yield of (3.5 ± 1.0) × 10(-4) is estimated, suggesting that although this process does occur, it does so with a very low efficiency.     

DFT study of the glutathione peroxidase-like activity of phenylselenol incorporating solvent-assisted proton exchange

Modeling of the glutathione peroxidase-like activity of phenylselenol has been accomplished using density-functional theory and solvent-assisted proton exchange (SAPE). SAPE is a modeling technique intended to mimic solvent participation in proton transfer associated with chemical reaction. Within this method, explicit water molecules incorporated into the gas-phase model allow relay of a proton through the water molecules from the site of protonation in the reactant to that in the product. The activation barriers obtained by SAPE for the three steps of the GPx-like mechanism of PhSeH fall within the limits expected for a catalytic system at physiological temperatures (DeltaG(1)++ = 19.1 kcal/mol; DeltaG(2)++= 6.6 kcal/mol; G(3)++ = 21.7 kcal/mol) and are significantly lower than studies which require direct proton transfer. The size of the SAPE network is also considered for the model of the reduction of the selenenic acid, step 2 of the GPx-like cycle. Use of a four-water network better accommodates the reaction pathway and reduces the activation barrier by 5 kcal/mol over the two-water model.     

Infrared overtone spectroscopy and unimolecular decay dynamics of peroxynitrous acid

Peroxynitrous acid (HOONO) is generated in a pulsed supersonic expansion through recombination of photolytically generated OH and NO(2) radicals. A rotationally resolved infrared action spectrum of HOONO is obtained in the OH overtone region at 6971.351(4) cm(-1) (origin), providing definitive spectroscopic identification of the trans-perp (tp) conformer of HOONO. Analysis of the rotational band structure yields rotational constants for the near prolate asymmetric top, the ratio of the a-type to c-type components of the transition dipole moment for the hybrid band, and a homogeneous linewidth arising from intramolecular vibrational energy redistribution and/or dissociation. The quantum state distribution of the OH (nu=0,J(OH)) products from dissociation is well characterized by a microcanonical statistical distribution constrained only by the energy available to products, 1304+/-38 cm(-1). This yields a 5667+/-38 cm(-1) [16.2(1) kcal mol(-1)] binding energy for tp-HOONO. An equivalent available energy and corresponding binding energy are obtained from the highest observed OH product state. Complementary high level ab initio calculations are carried out in conjunction with second-order vibrational perturbation theory to predict the spectroscopic observables associated with the OH overtone transition of tp-HOONO including its vibrational frequency, rotational constants, and transition dipole moment. The same approach is used to compute frequencies and intensities of multiple quantum transitions that aid in the assignment of weaker features observed in the OH overtone region, in particular, a combination band of tp-HOONO involving the HOON torsional mode.     

Measurement of the fourth O-H overtone absorption cross section in acetic acid using cavity ring-down spectroscopy

We report the absolute absorption cross sections of the fourth vibrational O-H (5ν(OH)) overtone in acetic acid using cavity ring-down spectroscopy. For compounds that undergo photodissociation via overtone excitation, such intensity information is required to calculate atmospheric photolysis rates. The fourth vibrational overtone of acetic acid is insufficiently energetic to effect dissociation, but measurement of its cross section provides a model for other overtone transitions that can affect atmospheric photochemistry. Though gas-phase acetic acid exists in equilibrium with its dimer, this work shows that only the monomeric species contributes to the acetic acid overtone spectrum. The absorption of acetic acid monomer peaks at ～615 nm and has a peak cross section of 1.84 × 10(-24) cm(2)·molecule(-1). Between 612 and 620 nm, the integrated cross section for the acetic acid monomer is (5.23 ± 0.73) × 10(-24) cm(2)·nm·molecule(-1) or (1.38 ± 0.19) × 10(-22) cm(2)·molecule(-1)·cm(-1). This is commensurate with the integrated cross section values for the fourth O-H overtone of other species. Theoretical calculations show that there is sufficient energy for hydrogen to transition between the two oxygen atoms, which results in an overtone-induced conformational change.     

Heats of formation and singlet-triplet separations of hydroxymethylene and 1-hydroxyethylidene

Thermochemical parameters of hydroxymethylene (HC:OH) and 1-hydroxyethylidene (CH3C:OH) were evaluated by using coupled-cluster, CCSD(T), theory, in conjunction with the augmented correlation consistent, aug-cc-pVnZ, basis sets, with n = D, T, Q, and 5, extrapolated to the complete basis set limit. The predicted value at 298 K for Delta Hf(CH2O) is -26.0 +/- 1 kcal/mol, as compared to an experimental value of -25.98 +/- 0.01 kcal/mol, and for Delta Hf(CH:OH) it is 26.1 +/- 1 kcal/mol. The hydroxymethylene-formaldehyde energy gap is 52.1 +/- 0.5 kcal/mol, the singlet-triplet separation of hydroxymethylene is Delta E(ST)(HC:OH) = 25.3 +/- 0.5 kcal/mol, the proton affinity is PA(HC:OH) = 222.5 +/- 0.5 kcal/mol, and the ionization energy is IEa(HC:OH) = 8.91 +/- 0.04 eV. The predicted value at 298 K for Delta Hf(CH3CHO) is -39.1 +/- 1 kcal/mol as compared to an experimental value of -40.80 +/- 0.35 kcal/mol, and for Delta Hf(CH3C:OH) it is 11.2 +/- 1 kcal/mol. The hydroxyethylidene-acetaldehyde energy gap is 50.6 +/- 0.5 kcal/mol, the singlet-triplet separation of 1-hydroxyethylidene is Delta E(ST)(CH3C:OH) = 30.5 +/- 0.5 kcal/mol, the proton affinity is PA(CH3C:OH) = 234.7 +/- 0.5 kcal/mol, and the ionization energy is IEa(CH3C:OH) = 8.18 +/- 0.04 eV. The calculated energy differences between the carbene and aldehyde isomers, and, thus, the heats of formation of the carbenes, differ from the experimental values by 2.5 kcal/mol.     

Computational mutation analysis of hydrogen abstraction and radical rearrangement steps in the catalysis of coenzyme B12-dependent diol dehydratase

A mutation analysis of the catalytic functions of active-site residues of coenzyme B(12)-dependent diol dehydratase in the conversion of 1,2-propanediol to 1,1-propanediol has been carried out by using QM/MM computations. Mutants His143Ala, Glu170Gln, Glu170Ala, and Glu170Ala/Glu221Ala were considered to estimate the impact of the mutations of His143 and Glu170. In the His143Ala mutant the activation energy for OH migration increased to 16.4 from 11.5 kcal mol(-1) in the wild-type enzyme. The highest activation energy, 19.6 kcal mol(-1), was measured for hydrogen back-abstraction in this reaction. The transition state for OH migration is not sufficiently stabilized by the hydrogen-bonding interaction formed between the spectator OH group and Gln170 in the Glu170Gln mutant, which demonstrates that a strong proton acceptor is required to promote OH migration. In the Glu170Ala mutant, a new strong hydrogen bond is formed between the spectator OH group and Glu221. A computed activation energy of 13.6 kcal mol(-1) for OH migration in the Glu170Ala mutant is only 2.1 kcal mol(-1) higher than the corresponding barrier in the wild-type enzyme. Despite the low activation barrier, the Glu170Ala mutant is inactive because the subsequent hydrogen back-abstraction is energetically demanding in this mutant. OH migration is not feasible in the Glu170Ala/Glu221Ala mutant because the activation barrier for OH migration is greatly increased by the loss of COO(-) groups near the spectator OH group. This result indicates that the effect of partial deprotonation of the spectator OH group is the most important factor in reducing the activation barrier for OH migration in the conversion of 1,2-propanediol to 1,1-propanediol catalyzed by diol dehydratase.