Ab initio study of gas phase and water-assisted tautomerization of maleimide and formamide

Maleimide serves as an important starting material in the synthesis of drugs and enzyme inhibitors. In the present paper, knowing the importance of tautomerization in maleimide for its drug action, potential energy surface of maleimide is studied and its tautomerization has been discussed and compared with tautomerization of formamide. Gas phase tautomerization of maleimide requires large amount of energy (23·21 kcal/mol) in comparison to formamide (15·05 kcal/mol) at HF/6-31+G* level. Thus making the proton transfer reaction a difficult process in gas phase. Water molecule lowers the energy barrier of tautomerization thus facilitating the tautomerization of maleimide to 5-hydroxy-pyrrol-2-one. Water assisted tautomerization of maleimide requires 19·60 kcal/mol energy at HF/6-31+G* and 17·63 kcal/mol energy at B3LYP/6-31+G* level, a decrease of 3·61 and 5·96 kcal/mol over gas phase tautomerization. Whereas, tautomerization of formamide requires 14·16 and 12·84 kcal/mol energy, a decrease of 0·89 and 2·01 kcal/mol energy over gas phase tautomerization at HF/6-31+G* and B3LYP/6-31+G* level, respectively. Water-assisted tautomerization in maleimide and formamide showed that difference in energy barrier reduces to 2·83 kcal/mol from 10·41 kcal/mol (in gas phase) at B3LYP level, which resulted that maleimide readily undergoes tautomerization in water molecule.     

Theoretical study of 5-aminolevulinic acid tautomerization: a novel self-catalyzed mechanism

5-Aminolevulinic acid (5ALA) is the key synthetic building block in protoporphyrin IX (PpIX), the heme chromophore in mitochondria. In this study density functional theory calculations were performed on the tautomers of 5ALA and the tautomerization reaction mechanism from its enolic forms (5-amino-4-hydroxypent-3-enoic acid and 5-amino-4-hydroxypent-4-enoic acid) to the more stable 5ALA. The hydrated form 5-amino-4,4-dihydroxypentanoic acid was also studied. The lowest energy pathway of 5ALA tautomerization is by means of autocatalysis, in that an oxygen of the carboxylic group transfers the hydrogen atom as a "crane", with an activation energy of approximately 15 kcal/mol. This should be compared to the barriers of about 35 kcal/mol for water assisted tautomerization, and 60 kcal/mol for direct hydrogen transfer. For hydration of 5ALA, the water catalyzed activation barrier is found to be approximately 35 kcal/mol, approximately 5 kcal/mol lower than direct hydration.     

Dimers of and tautomerism between 2-pyrimidinethiol and 2(1H)-pyrimidinethione: a density functional theory (DFT) study

Density functional theory (BLYP, B3LYP, B3P86, B3PW91) with the 6-31+G(d,p), 6-311+G(d,p), and cc-pVTZ basis sets has been used to calculate structural parameters, relative energies, and vibrational spectra of 2-pyrimidinethiol (1) and 2(1H)-pyrimidinethione (2) and their hydrogen-bonded homodimers (C(2) 3, C(2h) [4](double dagger), C(2h) 5), monohydrates, and dihydrates and a heterodimer (6). Several transition state structures proposed for the tautomerization process have also been examined. At the B3PW91/6-311+G(d,p)//B3PW91/6-31+G(d,p) level of theory 2-pyrimidinethiol (1) is predicted to be 3.41 kcal/mol more stable (E(rel)) than 2(1H)-pyrimidinethione (2) in the gas phase and 2 is predicted to be 6.47 kcal/mol more stable than 1 in aqueous medium. An unfavorable planar intramolecular strained four center transition state (TS1) for the tautomerization of 1 and 2 in the gas-phase lies 29.07 kcal/mol higher in energy than 2-pyrimidinethiol (1). The C(2) 2-pyrimidinethiol dimer (3) is 6.84 kcal/mol lower in energy than the C(2) homodimer transition state structure ([11](double dagger)) that connects dimers 3 and 4. Transition state [11](double dagger) provides a facile pathway for tautomerization between 1 and 2 in the gas phase (monomer-dimer promoted tautomerization). The hydrogen bonded 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O structures are predicted to be 1.27 and 1.55 kcal/mol, respectively, higher in energy than 2(1H)-pyrimidinethione- - -H(2)O and 2(1H)-pyrimidinethione- - -2H(2)O. Water promoted tautomerization via cyclic transition states involving one water molecule (TS- - -H(2)O, [12](double dagger)) and two water molecules (TS- - -2H(2)O, [13](double dagger)) lie 11.42 and 11.44 kcal/mol, respectively, higher in energy than 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O. Thus, the hydrated transition states [12](double dagger) and [13](double dagger) are involved in the tautomerism between 1 and 2 in aqueous medium.     

Theoretical studies of proton-transfer reactions of 2-hydroxypyridine--(H2O)n (n = 0-2) in the ground and excited states

The potential energy profiles for proton-transfer reactions of 2-hydroxypyridine and its complexes with water were determined by MP2, CASSCF and MR-CI calculations with the 6-31G** basis set. The tautomerization reaction between 2-hydroxypyridine (2HP) and 2-pyridone (2PY) does not take place at room temperature because of a barrier of approximately 35 kcal/mol for the ground-state pathway. The water-catalyzed enol-keto tautomerization reactions in the ground state proceed easily through the concerted proton transfer, especially for the two-water complex. The S1 tautomerization between the 2HP and 2PY monomers has a barrier of 18.4 kcal/mol, which is reduced to 5.6 kcal/mol for the one-water complex and 6.4 kcal/mol for the two-water complex. The results reported here predict that the photoinduced tautomerization reaction between the enol and keto forms involves a cyclic transition state having one or two water molecules as a bridge.     

2-Pyridinethiol/2-pyridinethione tautomeric equilibrium. A comparative experimental and computational study

The gas phase and solvent dependent preference of the tautomerization between 2-pyridinethiol (2SH) and 2-pyridinethione (2S) has been assessed using variable temperature Fourier transform infrared (FTIR) experiments, as well as ab initio and density functional theory computations. No spectroscopic evidence (nu(S)(-)(H) stretch) for 2SH was observed in toluene, C(6)D(6), heptane, or methylene chloride solutions. Although, C(s)() 2SH is 2.61 kcal/mol more stable than C(s)() 2S (CCSD(T)/cc-pVTZ//B3LYP/6-311+G(3df,2p)+ZPE), cyclohexane solvent-field relative energies (IPCM-MP2/6-311+G(3df,2p)) favor 2S by 1.96 kcal/mol. This is in accord with the FTIR observations and in quantitative agreement with the -2.6 kcal/mol solution (toluene or C(6)D(6)) calorimetric enthalpy for the 2S/2SH tautomerization favoring the thione. As the intramolecular transition state for the 2S, 2SH tautomerization (2TS) lies 25 (CBS-Q) to 30 kcal/mol (CCSD/cc-pVTZ) higher in energy than either tautomer, tautomerization probably occurs in the hydrogen bonded dimer. The B3LYP/6-311+G(3df,2p) optimized C(2) 2SH dimer is 10.23 kcal/mol + ZPE higher in energy than the C(2)(h)() 2S dimer and is only 2.95 kcal/mol + ZPE lower in energy than the C(2) 2TS dimer transition state. Dimerization equilibrium measurements (FTIR, C(6)D(6)) over the temperature range 22-63 degrees C agree: K(eq)(298) = 165 +/- 40 M(-)(1), DeltaH = -7.0 +/- 0.7 kcal/mol, and DeltaS = -13.4 +/- 3.0 cal/(mol deg). The difference between experimental and B3LYP/6-311+G(3df,2p) [-34.62 cal/(mol deg)] entropy changes is due to solvent effects. The B3LYP/6-311+G(3df,2p) nucleus independent chemical shifts (NICS) are -8.8 and -3.5 ppm 1 A above the 2SH and 2S ring centers, respectively, and the thiol is aromatic. Although the thione is not aromatic, it is stabilized by the thioamide resonance. In solvent, the large 2S dipole, 2-3 times greater than 2SH, favors the thione tautomer and, in conclusion, 2S is thermodynamically more stable than 2SH in solution.     

Solvent-mediated tautomerization of purine: single to quadruple proton transfer

We present calculations for the structures and the tautomerization reaction of purine and purine – (H₂O)_n (n=1–3) clusters. We find two pathways (via the carbene and the sp³-type intermediate) for the 9 ↔ 7 tautomerization of bare purine. The barrier heights for the 9 → 3 and 9 → 7 tautomerization of bare purine are calculated to be large (60–70 kcal/mol). Hydrogen bonding with the water molecule(s), however, dramatically lowers the 9 → 3 barrier by the concerted multiple proton transfer mechanism, favoring the formation of the conformer 3(H)- relative to the 7(H)-purine in the microsolvated environment, in contrast to the gas phase or the aqueous solution.     

Gas-phase thermochemical properties of pyrimidine nucleobases

The gas-phase acidity and proton affinity of thymine, cytosine, and 1-methyl cytosine have been examined using both theoretical (B3LYP/6-31+G*) and experimental (bracketing, Cooks kinetic) methods. This paper represents a comprehensive examination of multiple acidic sites of thymine and cytosine and of the acidity and proton affinity of thymine, cytosine, and 1-methyl cytosine. Thymine exists as the most stable "canonical" tautomer in the gas phase, with a DeltaH(acid) of 335 +/- 4 kcal mol(-1) (DeltaG(acid) = 328 +/- 4 kcal mol(-1)) for the more acidic N1-H. The acidity of the less acidic N3-H site has not, heretofore, been measured; we bracket a DeltaH(acid) value of 346 +/- 3 kcal mol(-1) (DeltaG(acid) = 339 +/- 3 kcal mol(-1)). The proton affinity (PA = DeltaH) of thymine is measured to be 211 +/- 3 kcal mol(-1) (GB = DeltaG = 203 +/- 3 kcal mol(-1)). Cytosine is known to have several stable tautomers in the gas phase in contrast to in solution, where the canonical tautomer predominates. Using bracketing methods in an FTMS, we measure a DeltaH(acid) for the more acidic site of 342 +/- 3 kcal mol(-1) (DeltaG(acid) = 335 +/- 3 kcal mol(-1)). The DeltaH(acid) of the less acidic site, previously unknown, is 352 +/- 4 kcal mol(-1) (345 +/- 4 kcal mol(-1)). The proton affinity is 228 +/- 3 kcal mol(-1) (GB = 220 +/- 3 kcal mol(-1)). Comparison of these values to calculations indicates that we most likely have a mixture of the canonical tautomer and two enol tautomers and possibly an imine tautomer under our conditions in the gas phase. We also measure the acidity and proton affinity of cytosine using the extended Cooks kinetic method. We form the proton-bound dimers via electrospray of an aqueous solution, which favors cytosine in the canonical form. The acidity of cytosine using this method is DeltaH(acid) = 343 +/- 3 kcal mol(-1), PA = 227 +/- 3 kcal mol(-1). We also examined 1-methyl cytosine, which has fewer accessible tautomers than cytosine. We measure a DeltaH(acid) of 349 +/- 3 kcal mol(-1) (DeltaG(acid) = 342 +/- 3 kcal mol(-1)) and a PA of 230 +/- 3 kcal mol(-1) (GB = 223 +/- 3 kcal mol(-1)). Our ultimate goal is to understand the intrinsic reactivity of nucleobases; gas-phase acidic and basic properties are of interest for chemical reasons and also possibly for biological purposes because biological media can be quite nonpolar.     

Spectroscopic evidence of keto-enol tautomerism in deliquesced malonic acid particles

Scanning transmission X-ray microscopy combined with near-edge X-ray absorption fine structure spectroscopy (STXM/NEXAFS) and optical microscopy coupled with Fourier transform infrared spectroscopy (micro-FTIR) have been applied to observe hygroscopic growth and chemical changes in malonic acid particles deposited on substrates. The extent of the hygroscopic growth of particles has been quantified in terms of the corresponding water-to-solute ratios (WSR) based on STXM/NEXAFS and micro-FTIR data sets. WSR values derived separately from two applied methods displayed a remarkable agreement with previous data reported in the literature. Comparison of NEXAFS and FTIR spectra acquired at different relative humidity (RH) shows efficient keto-enol tautomerization of malonic acid, with the enol form dominating at higher RH. The keto-enol equilibrium constants were calculated using relevant peak intensities in the carbon and oxygen K-edge NEXAFS spectra as a function of RH. We report strong dependence of the equilibrium constant on RH, with measured values of 0.18 ± 0.03, 1.11 ± 0.14, and 2.33 ± 0.37 corresponding to 2, 50, and 90% RH, respectively. Enols are important intermediates in aldol condensation reactions pertaining to formation and atmospheric aging of secondary organic aerosol (SOA). The present knowledge assumes that constituents of atmospheric deliquesced particles undergo aqueous chemistry with kinetic and equilibrium constants analogous to reactions in bulk solutions, which would estimate absolute dominance of the keto form of carboxylic acids. For instance, the keto-enol equilibrium constant of malonic acid in diluted aqueous solution is <10(-4). Our results suggest that in deliquesced micrometer-size particles, carboxylic acids may exist in predominantly enol forms that need to be explicitly considered in atmospheric aerosol chemistry.     

Overtone-induced decarboxylation: a potential sink for atmospheric diacids

Atmospheric photochemistry induced by solar excitation of vibrational overtone transitions has recently been demonstrated to be of importance in cleaving weak bonds (in HO(2)NO(2)) and inducing intramolecular rearrangement followed by reaction (in H(2)SO(4)). Here, we propose another potentially important process: the decarboxylation of organic acids. To demonstrate this possibility, we have calculated the decarboxylation pathways for malonic acid and its monohydrate. The barrier to the gas-phase decarboxylation was calculated to be in the range 26-28 kcal/mol at the B3LYP/6-311++G(3df,3pd) level of theory, in good agreement with previous results. The transition state is a six-membered ring structure which is accessed via concerted O-H and C-C stretches; excitation of v(OH) > or = 3 of either one of the OH stretching modes is sufficient to supply the energy needed for the decarboxylation. A low-energy isomer of the malonic acid-water complex forms an eight-membered, multiply hydrogen bonded structure, bound by 3-6 kcal/mol, somewhat less stable than the lowest energy, six-membered ring isomer. Decarboxylation of such complexes uses water as a catalyst; the water accepts an acidic proton from one malonic acid group and transfers a proton to the carbonyl of the other acid group. The barrier for this process is 20-22 kcal/mol, suggesting that complexes excited to v(OH) > or = 2 possess sufficient energy to react. Using estimated absorption cross sections for the OH overtone transitions, we suggest that the overtone-induced decarboxylation of malonic acid and its water complex is competitive with wet deposition of the acid and with gas-phase reaction with OH for removal of the acid.     

Theoretical studies of the tautomeric equilibria for five-member N-heterocycles in the gas phase and in solution

Tautomeric equilibria have been studied for five-member N-heterocycles and their methyl derivatives in the gas phase and in different solvents with dielectric constants of epsilon = 4.7-78.4. The free energy changes differently for tautomers upon solvation as compared to the gas phase, resulting in a shift of the equilibrium constant in solution. Solvents with increasing dielectric constant produce more negative solute-solvent interaction energies and increasing internal energies. The methyl-substituted imidazole and pyrrazole form delicate equilibria between two tautomeric forms. Depending on the solvent, the methyl-substituted triazoles and tetrazole have one or two major tautomers in solution. When estimating the relative solvation free energies by means of an explicit solvent model and using the FEP/MC method, one observes that the preferred tautomers differ in several cases from those predicted by the continuum solvent model. The 1,2-prototropic shift, as an intramolecular tautomerization path, requires about 50 kcal/mol activation energy for imidazole in the gas phase, and this route is also disfavored in a solution. The calculated activation free energy along the intramolecular path is 48-50 kcal/mol in chloroform and water as compared to a literature value of 13.6 kcal/mol for pyrrazole in DMSO. A molecular dynamics computer experiment favors the formation of an imidazole chain in chloroform, making the 1,3-tautomerization feasible along an intermolecular path in nonprotic solvents. In aqueous solution, one strong N-H...Ow hydrogen bond is formed for each species, whereas all other nitrogens in the ring form weaker, N...HwOw type hydrogen bonds. The tetrahydrofuran solvent acts as a hydrogen bond acceptor and forms N-H...Oether bonds. Molecules of the dichloromethane solvent are in favorable dipole-dipole interactions with the solute. The results obtained are useful in the design of N-heterocyclic ligands forming specified hydrogen bonds with protein side chains.     

The acidity and proton affinity of the damaged base 1,N6-ethenoadenine in the gas phase versus in solution: intrinsic reactivity and biological implications

1,N(6)-ethenoadenine (epsilonA) is a highly mutagenic lesion that is excised from human DNA by the enzyme alkyladenine DNA glycosylase (AAG). In an effort to understand the intrinsic properties of 1,N(6)-ethenoadenine, we examined its gas phase acidity and proton affinity using quantum mechanical calculations and mass spectrometric experimental methods. We measure two acidities for epsilonA: a more acidic site (DeltaH(acid) = 332 kcal mol(-1); DeltaG(acid) = 325 kcal mol(-1)) and a less acidic site (DeltaH(acid) = 367 kcal mol(-1); DeltaG(acid) = 360 kcal mol(-1)). We also find that the proton affinity of the most basic site of 1,N(6)-ethenoadenine is 232-233 kcal mol(-1) (GB = 224 kcal mol(-1)). These measurements, when compared to calculations, establish that, under our experimental conditions, we have only the canonical tautomer of 1,N(6)-ethenoadenine present. We also compare the gas phase acidic properties of epsilonA with that of the normal bases adenine and guanine and find that epsilonA is the most acidic. This supports the theory that AAG and other related enzymes may cleave damaged bases as anions. Furthermore, comparison of the gas phase and aqueous acidities indicates that the nonpolar environment of the enzyme enhances the acidity differences of epsilonA versus adenine and guanine.     

Ab initio study of the internal rotation barrier of formamide and the formamide–H2O complex

The barriers to internal rotation about the N? C bond of formamide and the formamide–H₂O complex have been examined by ab initio quantum chemistry methods. Both self-consistent field and correlated approaches have been utilized to determine the geometries, energies, and local harmonic vibrational frequencies of the minimum-energy and transition-state structures of these two systems. We find that formamide's rotation barrier of 14–15 kcal/mol is increased to 16–18 kcal/mol when a single H₂O molecule is attached. This result contrasts with the effect of a single H₂O solvent molecule on the barrier to tautomerization of formamide (to form formamidic acid HN? CHOH) for which a barrier lowering of ca. 20 kcal/mol has been observed. The rotation barriers obtained for formamide and for its H₂O complex are compared with barriers obtained experimentally in various solvents. © 1993 John Wiley & Sons, Inc.     

A direct ab initio dynamics study of the water-assisted tautomerization of formamide

Direct ab initio dynamics calculations based on a canonical variational transition-state theory with several multidimensional semiclassical tunneling approximations were carried out to obtain rate constants for the water-assisted tautomerization of formamide. The accuracy of the density functionals, namely, B-LYP, B3-LYP, and BH&H-LYP, were examined. We found that the BH&H-LYP method yields the most accurate transition-state properties when comparing it to ab initio MP2 and QCISD results, whereas B-LYP and B3-LYP methods predict barrier heights too low. Reaction path information was calculated at both the MP2 and nonlocal hybrid BH&H-LYP levels using the 6–31G(d,p) basis set. At the BH&H-LYP level, we found that the zero-point energy motion lowers the barrier to tautomerization in the formamide-water complex by 3.6 kcal/mol. When tunneling is considered, the activation energy at the BH&H-LYP level at 300 K is 17.1 kcal/mol. This is 3.4 kcal/mol below the zero-point-corrected barrier and 7.0 kcal/mol below the classical barrier. Excellent agreement between BH&H-LYP and MP2 rate constants further supports the use of BH&H-LYP for rate calculations of large systems. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 861–874, 1997     

水溶液中亚硝基甲烷与甲醛肟之间互变异构反应的理论研究

采用超分子-连续介质(PCM)模型,在密度泛函B3LYP/6-311++G^**水平上对水溶液中亚硝基甲烷异构化反应的机理进行了理论研究.结果表明,在水溶液中亚硝基甲烷可以通过两条反应途径(Ⅰ和Ⅱ)经质子迁移得到更稳定的重排产物--反式甲醛肟,但优势反应途径与在气相反应不同.在水溶液中亚硝基甲烷异构化反应最有可能的途径Ⅰ是通过氢迁移先生成顺式甲醛肟,然后绕N-O键旋转生成更稳定的反式甲醛肟.并且由于水分子的催化作用使得反应活化能从气相中240.6和196.2 kJ/mol分别降低至水溶液中的61.7和92.1 kJ/mol.     

22-hydroxybenziporphyrin: switching of antiaromaticity by phenol-keto tautomerization

22-hydroxybenziporphyrin, a porphyrin analogue containing a phenol moiety, has been shown to exist as an equilibrium mixture of two distinctly different tautomers. One of them actually contains the hydroxy group and shows the local [6]annulene aromaticity in the phenol fragment. The other tautomer contains a keto group and exhibits a [20]annulenoid structure characterized by macrocyclic antiaromaticity. The tautomerization process has been investigated in detail using variable-temperature 1H NMR spectroscopy. The process is very fast, with an estimated activation energy of ca. 30 kJ/mol. Further insight into the energetics of the tautomerization is obtained from density functional (DFT) calculations. Surprisingly, the estimated energy of the antiaromatic keto species is 3-5 kcal/mol lower than the energy of the phenolic tautomer. The geometric and magnetic manifestations of aromaticity and antiaromaticity in the two tautomers are probed using a number of computational devices, including Wiberg bond indices, resonace weights derived from the harmonic oscillator model, and nucleus-independent chemical shifts. It is shown that mixing of phenolic and keto contributions in both tautomers is stronger than that in related tautomers of phenol. This effect is caused by extensive conjugation with the tripyrrolic unit of 22-hydroxybenziporphyrin and, to a lesser extent, by intramolecular hydrogen bonding.     

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.     

Standard enthalpies of the formation of malonic acid and products of its dissociation in an aqueous solution

Enthalpies of the dissolution of malonic acid in aqueous solutions of perchloric acid and sodium perchlorate were measured at ionic strength I = 1.0; 1.5; 2.0; 2.5 mol/l and T = 298.15 K by calorimetry. The standard enthalpies of the formation of malonic acid and the products of its dissociation in a aqueous solution were calculated.     

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.     

Reaction mechanism and tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution: a combined Monte Carlo and quantum mechanics study

A combined Monte Carlo and quantum mechanical study was carried out to analyze the tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution. Second- and fourth-order M?ller-Plesset perturbation theory calculations indicate that in the gas phase thiol (Pym-SH) is more stable than the thione (Pym-NH) by ca. 8 kcal/mol. In aqueous solution, thermodynamic perturbation theory implemented on a Monte Carlo NpT simulation indicates that both the differential enthalpy and Gibbs free energy favor the thione form. The calculated differential enthalpy is DeltaH(SH)(-->)(NH)(solv) = -1.7 kcal/mol and the differential Gibbs free energy is DeltaG(SH)(-->)(NH)(solv) = -1.9 kcal/mol. Analysis is made of the contribution of the solute-solvent hydrogen bonds and it is noted that the SH group in the thiol and NH group in the thione tautomers act exclusively as a hydrogen bond donor in aqueous solution. The proton transfer reaction between the tautomeric forms was also investigated in the gas phase and in aqueous solution. Two distinct mechanisms were considered: a direct intramolecular transfer and a water-assisted mechanism. In the gas phase, the intramolecular transfer leads to a large energy barrier of 34.4 kcal/mol, passing through a three-center transition state. The proton transfer with the assistance of one water molecule decreases the energy barrier to 17.2 kcal/mol. In solution, these calculated activation barriers are, respectively, 32.0 and 14.8 kcal/mol. The solvent effect is found to be sizable but it is considerably more important as a participant in the water-assisted mechanism than the solvent field of the solute-solvent interaction. Finally, the calculated total Gibbs free energy is used to estimate the equilibrium constant.     

Dimerization and double proton transfer-induced tautomerism of 4(3H)-pyrimidinone in solution studied by IR spectroscopy and quantum chemical calculations

The tautomerism and dimerization of 4(3H)-pyrimidinone (4(3H)Pyr) in carbon tetrachloride (CCl4) and chloroform (CHCl3) solutions were investigated using IR spectroscopy and quantum chemical calculations. The observed IR spectra in the NH and OH stretching regions clearly revealed the predominance of the keto tautomer in both solvent systems. The enol form only exists in a very small proportion in the CCl4 solution. The tautomeric constant for the two monomers KT[OH/NH] = 0.012 and DeltaE = 2.62 kcal/mol were estimated at 25 degrees C. This result was supported by the self-consistent reaction field/polarizable continuum (SCRF/PCM) calculation at the MP4(full, SDQ)/aug-cc-pVDZ level, which predicted DeltaE = 3.06 kcal/mol in CCl4. In the C=O stretching region, two bands were observed, suggesting the coexistence of two keto structures at equilibrium. The calculated IR spectra indicated that the bands at 1711 and 1675 cm(-1) arise from the keto monomer and keto-keto (KK) ring dimer, respectively. At elevated temperature, the populations of both the keto and enol monomers increased for the CCl4 solution. The present study revealed that the keto <--> enol tautomerization does not occur in the isolated monomer molecule. The double proton transfer (DPT) reaction in the KK ring dimer presumably plays a substantial role in the population increase of the enol monomer. To our knowledge, this may be the first observation of the tautomerization in a model base pair via the temperature-induced ground-state DPT reaction under a nonpolar liquid environment reported so far. This tautomerism can serve as a mimic circumstance for the spontaneous mutations induced by proton transfer in the DNA base pairs.