Double proton transfer behavior and one-electron oxidation effect in double H-bonded glycinamide-formic acid complex

The behavior of double proton transfer occurring in a representative glycinamide-formic acid complex has been investigated at the B3LYP/6-311 + + G( * *) level of theory. Thermodynamic and, especially, kinetic parameters, such as tautomeric energy, equilibrium constant, and barrier heights have been discussed, respectively. The relevant quantities involved in the double proton transfer process, such as geometrical changes, interaction energies, and intrinsic reaction coordinate calculations have also been studied. Computational results show that the participation of a formic acid molecule favors the proceeding of the proton transfer for glycinamide compared with that without mediate-assisted case. The double proton transfer process proceeds with a concerted mechanism rather than a stepwise one since no ion-pair complexes have been located during the proton transfer process. The calculated barrier heights are 11.48 and 0.85 kcal/mol for the forward and reverse directions, respectively. However, both of them have been reduced by 2.95 and 2.61 kcal/mol to 8.53 and -1.76 kcal/mol if further inclusion of zero-point vibrational energy corrections, where the negative barrier height implies that the reverse reaction should proceed with barrierless spontaneously, analogous to that occurring between glycinamide and formamide. Furthermore, solvent effects on the thermodynamic and kinetic processes have also been predicted qualitatively employing the isodensity surface polarized continuum model within the framework of the self-consistent reaction field theory. Additionally, the oxidation process for the double H-bonded glycinamide-formic acid complex has also been investigated. Contrary to that neutral form possessing a pair of two parallel intermolecular H bonds, only a single H bond with a comparable strength has been found in its ionized form. The vertical and adiabatic ionization potentials for the neutral complex have been determined to be about 9.40 and 8.69 eV, respectively, where ionization is mainly localized on the glycinamide fragment. Like that ionized glycinamide-formamide complex, the proton transfer in the ionized complex is characterized by a single-well potential, implying that the proton initially attached to amide N4 in the glycinamide fragment cannot be transferred to carbonyl O13 in the formic acid fragment at the geometry of the optimized complex.     

2-hydroxypyridine <--> 2-pyridone tautomerization: catalytic influence of formic acid

A 1:1 hydrogen-bonded complex between 2-pyridone and formic acid has been characterized using laser-induced-fluorescence excitation and dispersed fluorescence spectroscopy in a supersonic jet expansion. Under the same expansion condition, the fluorescence signal of the tautomeric form of the complex (2-hydroxypyridine...formic acid) is absent, although both the bare tautomeric molecules exhibit well-resolved laser-induced-fluorescence spectra. Quantum chemistry calculation at the DFT/B3LYP/6-311++G** level predicts that in the ground electronic state the activation barrier for tautomerization from hydroxy to keto form in bare molecules is very large (approximately 34 kcal/mol). However, the process turns out to be nearly barrierless when assisted by formic acid, and double proton transfer occurs via a concerted mechanism.     

Fragmentation and deformation mechanism of glycine isomers in gas phase: Investigations of charge effect

The structural parameters, relative stability, proton transfer energy barriers of four typical and life related isomers and conformers of different charged (n=0,+/-1,+/-2) glycine species have been investigated using B3LYP, BHLYP, and CCSD(T) methods. Results indicate that those neutral and (+/-1)-charged species are stable. For the (+2)-charged cases, all four triplet-state glycine species and only the singlet-state zwitterionic one are stable. On the other hand, only the singlet-state zwtterionic glycine ((1)GlyZW(-2)) and the corresponding neutral form counterpart ((1)Gly(-2)) are stable for the (-2)-charged cases. Either of the two stable structures holds a proton lying in the position (2-3 A) of being separated from its corresponding parental species. Those unstable divalent glycine species are dissociated into different smaller species spontaneously according to the characters of their different structures and electron spins. The presented fragmentation and deformation mechanisms can effectively predict and satisfactorily explain some experimental phenomena, which had been puzzling the mass spectrometry chemists. Also, the mechanisms should be suitable for any other similar molecule systems. Comparisons of the relative energies of the four (+1)-charged glycine species show that doublet-state glycine III ((2)GlyIII1) is more stable in energy by 12.1 kcal/mol than the (+1)-charged glycine Gly ((2)Gly1). This is consistent with the energy ordering of their corresponding mono-valence metal ion-bound derivatives. In addition, calculations show that an intramolecular proton transfer of (2)Gly(-1) to become its zwitterionic counterpart is preferred due to its least activation energy barrier (5.8 kcal/mol) among four discussed processes.     

Double proton transfer in the isolated and DNA-embedded guanine-cytosine base pair

The energetics and dynamics of double proton transfer (DPT) is investigated theoretically for the Watson-Crick conformation of the guanine-cytosine (GC) base pair. Using semiempirical density functional theory the isolated and DNA-embedded GC pair is considered. Differences in the energetics and dynamics of DPT thus addresses the question of how relevant studies of isolated base pairs are for the understanding of processes occurring in DNA. Two-dimensional potential energy surfaces involving the transferring hydrogen atoms and the proton donors and acceptors are presented for both systems. The DPT reaction is accompanied by a contraction of the distance between the two bases with virtually identical energetic barriers being 18.8 and 18.7 kcal/mol for the isolated and DNA-embedded system, respectively. However, the transition state for DPT in the DNA-embedded GC pair is offset by 0.1 A to larger N-H separation compared to the isolated GC pair. Using activated ab initio molecular dynamics, DPT is readily observed for the isolated base pair with a minimal amount of 21.4 kcal/mol of initial average kinetic energy along the DPT normal mode vector. On a time scale of approximately 100 fs DPT has occurred and the excess energy is redistributed. For the DNA-embedded GC pair considerably more kinetic energy is required (30.0 kcal/mol) for DPT and the process is completed within one hydrogen vibration. The relevance of studies of isolated base pairs and base pair analogs in regard of reactions or properties involving DNA is discussed.     

Proton exchanges between phenols and ammonia or amines: a computational study

Density functional theory calculations at the B3LYP/6-31+G(d,p) level of theory have been performed to explore proton exchanges between phenols and ammonia or amines, which can be used to account for previous NMR experiments. For the parent phenol-NH(3) system, a transition state with a symmetric phenolate-NH(4)(+)-like structure, which lies about 35 kcal mol(-1) in energy above the hydrogen-bonded complex, has been successfully located. An intrinsic reaction coordinate (IRC) analysis indicates that the proton exchange is a concerted process, which can be roughly divided into four continuous subprocesses. A series of para-substituted phenol-NH(3) systems have been considered to investigate the substituent effect. Whereas introduction of an electron-withdrawing group on the phenol appreciably reduces the barrier, an opposite effect is observed for an electron-donating group. Moreover, it has been disclosed that there exists a good linear correlation between the activation barriers and the interaction energies between the phenols and NH(3), indicating the important role of proton transfer (or hydrogen bonding) in determining the proton exchange. Also considered are the proton exchanges between phenol and amines and those for some sterically hindered systems. The results show that the phenol tends to exchange hydrogen with the amines, preferably the secondary amines, and that the steric effect is favorable for the proton exchange, which imply that, as the IRC analysis suggested, besides the proton transfer, the flip of the ammonium-like moiety may play a significant role in the course of proton exchange. For all of these systems, we investigated the solvent effects and found that the barrier heights of proton exchange decrease remarkably as compared to those in a vacuum due to the ion pair feature of the transition state. Finally, we explored the phenol radical cation-NH(3) system; the barrierless proton transfer and remarkably low barrier (5.2 kcal mol(-1)) of proton exchange provide further evidence for the importance of proton transfer in the proton exchange.     

Aluminosilicate surfaces as promoters for peptide bond formation: an assessment of Bernal's hypothesis by ab initio methods

The role in prebiotic chemistry that Br?nsted and Lewis sites, both present at the surface of common aluminosilicates, may have played in favoring the peptide bond formation has been addressed by ab initio methods within a cluster approach. B3LYP/6-31+G(d,p) free energy potential energy surfaces have been fully characterized for the model reaction glycine + NH3 --> 2-NH2 acetamide (mimicking the true 2 Gly --> GlyGly one) occurring on (i) a Lewis site, (ii) a Br?nsted site, and (iii) a combined action of Lewis/Br?nsted sites. Compared to the gas-phase (gp) activation free energy of 50 kcal/mol, the Lewis site alone reduces the gp barrier to 41 kcal/mol, whereas the activation by the Br?nsted site dramatically reduces the barrier to about 18 kcal/mol. Nevertheless, formation of the prereactant complex in this latter case will rarely occur, since water will easily displace the glycine molecule interacting with the Br?nsted site. However, if a realistic feldspar surface with neighboring Br?nsted and Lewis sites is considered, the proper prereactant complex is highly stabilized by a simultaneous interaction with the Lewis and the Br?nsted sites, in such a way that the Lewis site strongly attaches the glycine molecule to the surface whereas the Br?nsted site efficiently catalyzes the condensation reaction, showing that the interplay between Lewis/Br?nsted sites is an important issue. The free energy barrier computed for the realistic feldspar surface model is 26 kcal/mol. The role of dispersive interactions on the free energy barrier and the stabilization of the final product, not accounted for by the B3LYP functional, have been estimated and shown to be substantial. Speculations about further elongation of the formed dipeptide have been put forward on the basis of the relatively strong interaction energy of the formed GlyGly dipeptide with the aluminosilicate surface.     

New model for a theoretical density functional theory investigation of the mechanism of the carbonic anhydrase: how does the internal bicarbonate rearrangement occur?

A theoretical density functional theory (DFT, B3LYP) investigation has been carried out on the catalytic cycle of the carbonic anhydrase. A model system including the Glu106 and Thr199 residues and the "deep" water molecule has been used. It has been found that the nucleophilic attack of the zinc-bound OH on the CO(2) molecule has a negligible barrier (only 1.2 kcal mol(-1)). This small value is due to a hydrogen-bond network involving Glu106, Thr199, and the deep water molecule. The two usually proposed mechanisms for the internal bicarbonate rearrangement have been carefully examined. In the presence of the two Glu106 and Thr199 residues, the direct proton transfer (Lipscomb mechanism) is a two-step process, which proceeds via a proton relay network characterized by two activation barriers of 4.4 and 9.0 kcal mol(-1). This pathway can effectively compete with a rotational mechanism (Lindskog mechanism), which has a barrier of 13.2 kcal mol(-1). The fast proton transfer found here is basically due to the effect of the Glu106 residue, which stabilizes an intermediate situation where the Glu106 fragment is protonated. In the absence of Glu106, the barrier for the proton transfer is much larger (32.3 kcal mol(-1)) and the Lindskog mechanism becomes favored.     

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.     

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.     

Different ESIPT Mechanisms for Angular‐Shaped Quinacridone in Toluene and Dimethyl Formamide (DMF) Solvents: A Theoretical Study

Adopting density functional theory (DFT) and time‐dependent density functional theory (TDDFT) methods, we investigat and present two different excited‐state intramolecular proton transfer (ESIPT) mechanisms of angular‐quinacridone (a‐QD) in both toluene and DMF,theoretically. Comparing the primary structural variations of a‐QD involved in the intramolecular hydrogen bond, we conclude that N1–H2⋯O3 should be strengthened in the S₁ state, which may facilitate the ESIPT process. Particularly, in toluene, the S₁‐state‐stable a‐QD enol* could not be located because of the non‐barrier ESIPT process. Concomitantly, infrared vibrational spectral analysis further verified the stability of the hydrogen bond. In addition, the role of charge–transfer interaction has been addressed under the frontier molecular orbitals (MOs), which depicts the nature of the electronic excited state and supports the ESIPT reaction. The potential energy curves according to variational N1–H2 coordinate demonstrates that the proton transfer process should occur spontaneously in toluene; however, in DMF, a low potential energy barrier of 0.493 kcal/mol is needed to complete the ESIPT reaction. Although this barrier of 0.493 kcal/mol is too low to make an important impact on the ESIPT reaction, just because of the existence of barrier, ESIPT mechanisms in toluene and DMF are different.     

Quantum and molecular mechanical study of the first proton transfer in the catalytic cycle of cytochrome P450cam and its mutant D251N

In the catalytic cycle of cytochrome P450cam, the hydroperoxo intermediate (Cpd 0) is formed by proton transfer from a reduced oxyheme complex (S5). This process is drastically slowed down when Asp251 is mutated to Asn (D251N). We report quantum mechanical/molecular mechanical (QM/MM) calculations that address this proton delivery in the doublet state through a hydrogen-bond network in the Asp251 channel, both for the wild-type enzyme and the D251N mutant, using four different active-site models. For the wild-type, we find a facile concerted mechanism for proton transfer from protonated Asp251 via Wat901 and Thr252 to the FeOO moiety, with a barrier of about 1 kcal/mol and a high exothermicity of more than 20 kcal/mol. In the D251N mutant with a neutral Asn251 residue, the proton transfer is almost thermoneutral or slightly exothermic in the three models considered. It is still very facile when the Asn251 residue adopts a conformation analogous to Asp251 in the wild-type enzyme, but the barrier increases significantly when the Asn251 side chain flips (as indicated by classical molecular dynamics simulations). This flip disrupts the hydrogen-bond network and hence the proton-transfer pathway, which causes a longer lifetime of S5 in the D251N mutant (consistent with experimental observations). The entry of an additional water molecule into the active site of D251N with flipped Asn251 regenerates the hydrogen-bond network and provides a viable mechanism for proton delivery in the mutant, with a moderate barrier of about 7 kcal/mol.     

Valence anions in complexes of adenine and 9-methyladenine with formic acid: stabilization by intermolecular proton transfer

Photoelectron spectra of adenine-formic acid (AFA(-)) and 9-methyladenine-formic acid (MAFA(-)) anionic complexes have been recorded with 2.540 eV photons. These spectra reveal broad features with maxima at 1.5-1.4 eV that indicate formation of stable valence anions in the gas phase. The neutral and anionic complexes of adenine/9-methyladenine and formic acid were also studied computationally at the B3LYP, second-order M?ller-Plesset, and coupled-cluster levels of theory with the 6-31++G** and aug-cc-pVDZ basis sets. The neutral complexes form cyclic hydrogen bonds, and the most stable dimers are bound by 17.7 and 16.0 kcal/mol for AFA and MAFA, respectively. The theoretical results indicate that the excess electron in both AFA(-) and MAFA(-) occupies a pi* orbital localized on adenine/9-methyladenine, and the adiabatic stability of the most stable anions amounts to 0.67 and 0.54 eV for AFA(-) and MAFA(-), respectively. The attachment of the excess electron to the complexes induces a barrier-free proton transfer (BFPT) from the carboxylic group of formic acid to a N atom of adenine or 9-methyladenine. As a result, the most stable structures of the anionic complexes can be characterized as neutral radicals of hydrogenated adenine (9-methyladenine) solvated by a deprotonated formic acid. The BFPT to the N atoms of adenine may be biologically relevant because some of these sites are not involved in the Watson-Crick pairing scheme and are easily accessible in the cellular environment. We suggest that valence anions of purines might be as important as those of pyrimidines in the process of DNA damage by low-energy electrons.     

Exploring a reaction mechanism for acetato ligand replacement in paddlewheel tetrakisacetatodirhodium (II,II) complex by ammonia: computational density functional theory study

This study focuses on the first step of interaction between DNA and the paddle-wheel dirhodium complex. The ammonia molecule was used to model the oligonucleotide sequence. The reaction was considered in neutral and acidic conditions, in gas phase, and in solvent, using the COSMO model. Molecular structures of the complexes were optimized in both models at the B3PW91/6-31G(d) level. The B3LYP functional and aug-cc-pvdz basis set were employed for single-point energy determination and electron distribution analyses. It was shown that in neutral solution the replacement of axial aqua ligand is mildly exoergic. The reaction is characterized by a relatively low activation barrier (10-12 kcal/mol), and, according to Eyring transition state theory, it proceeds very quickly. The breaking of the Rh-O(ac) bond in neutral solution is mildly endoergic (less than 1 kcal/mol) with an activation barrier of about 21 kcal/mol. However, this process can occur much more spontaneously (ΔG of -14 kcal/mol) when the dirhodium complex is protonated at the acetyl oxygen in remote position.     

Is the protein surrounding the active site critical for hydrogen peroxide reduction by selenoprotein glutathione peroxidase? An ONIOM study

In this ONIOM(QM:MM) study, we evaluate the role of the protein surroundings in the mechanism of H2O2 reduction catalyzed by the glutathione peroxidase enzyme, using the whole monomer (3113 atoms in 196 amino acid residues) as a model. A new optimization scheme that allows the full optimization of transition states for large systems has been utilized. It was found that in the presence of the surrounding protein the optimized active site structure bears a closer resemblance to the one in the X-ray structure than that without the surrounding protein. H2O2 reduction occurs through a two-step mechanism. In the first step, the selenolate anion (E-Se(-)) formation occurs with a barrier of 16.4 kcal/mol and is endothermic by 12.0 kcal/mol. The Gln83 residue plays the key role of the proton abstractor, which is in line with the experimental suggestion. In the second step, the O-O bond is cleaved, and selenenic acid (R-Se-OH) and a water molecule are formed. The calculated barrier for this process is 6.0 kcal/mol, and it is exothermic by 80.9 kcal/mol. The overall barrier of 18.0 kcal/mol for H2O2 reduction is in reasonable agreement with the experimentally measured barrier of 14.9 kcal/mol. The protein surroundings has been calculated to exert a net effect of only 0.70 kcal/mol (in comparison to the "active site only" model including solvent effects) on the overall barrier, which is most likely due to the active site being located at the enzyme surface.     

Does Silica Surface Catalyse Peptide Bond Formation? New Insights from First‐Principles Calculations

The role that silica surface could have played in prebiotic chemistry as a catalyst for peptide bond formation has been addressed at the B3LYP/6-31+G(d,p) level for a model reaction involving glycine and ammonia on a silica cluster mimicking an isolated terminal silanol group present at the silica surface. Hydrogen-bond complexation between glycine and the silanol is followed by the formation of the mixed surface anhydride Si(surf)-O-C(=O)-R, which has been suggested in the literature to activate the C=O bond towards nucleophilic attack by a second glycine molecule, here simulated by the simpler NH3 molecule. However, B3LYP/6-31+G(d,p) calculations show that formation of the surface mixed anhydride Si(surf)-O-C(=O)-R is disfavoured (delta(r)G298 approximately 6 kcal mol(-1)), and that the surface bond only moderately lowers the free-energy barrier of the nucleophilic attack responsible for peptide bond formation (deltaG298(double dagger) approximately 48 kcal mol(-1)) in comparison with the uncatalysed reaction (deltaG298(double dagger) approximately 52 kcal mol(-1)). A further decrease of the free-energy barrier of peptide bond formation (deltaG298(double dagger) approximately 41 kcal mol(-1)) is achieved by a single water molecule close to the reaction centre acting as a proton-transfer helper in the activated complex. A possible role of strained silica surface defects on the formation of the surface mixed anhydride Si(surf)-O-C(=O)-R has also been addressed.     

Mechanistic Study on the Monofunctional Binding of Hydrolysis Products of Non-planar Heterocyclic Complexes Trans-[PtCl_2NH_3（piperidine）] and Trans-[PtCl_2（piperidine）_2] to DNA Purine Bases

The monofunctional substitution reactions between trans-[PtCl(H2O)(NH3)(pip)]+,trans-[Pt(H2O)2(NH3)(pip)]2+,trans-[PtCl(H2O)(pip)2]+,trans-[Pt(H2O)2(pip)2]2+ (pip = piperidine) and adenine/guanine nucleotides are explored by using B3LYP hybrid functional and IEF-PCM salvation models. For the trans-[Pt(H2O)2(NH3)(pip)]2+ and trans-[PtCl(H2O)(NH3)(pip)]+ complexes,the computed barrier heights in aqueous solution are 13.5/13.5 and 11.6/11.6 kcal/mol from trans-Pt-chloroaqua complex to trans/cis-monoadduct for adenine and guanine,and the corresponding values are 20.7/20.7 and 18.8/18.8 kcal/mol from trans-Pt-diaqua complex to trans/cis-monoadduct for adenine and guanine,respectively. For trans-[PtCl(H2O)(pip)2]+ and trans-[Pt(H2O)2(pip)2]2+,the corresponding values are 21.5/21.3 and 19.4/19.4 kcal/mol,and 26.0/26.0 and 20.7/20.8 kal/mol for adenine and guanine,respectively. Our calculations demonstrate that the barrier heights of chloroaqua are lower than the corresponding values of diaqua for adenine and guanine. In addition,the free energies of activation for guanine in aqueous solution are all smaller than that for adenine,which predicts a preference of 1.9 kcal/mol when trans-[PtCl(H2O)(NH3)(pip)]+ and trans-[Pt(H2O)2(NH3)(pip)]2+ are the active agents and ～1.9 and ～ 5.3 kcal/mol when trans-[PtCl(H2O)(pip)2]+ and trans-[Pt(H2O)2(pip)2]2+ are the active agents,respectively. For the reaction of trans-Pt-chloroaqua (or diaqua) to cis-monoadduct,we obtain the same transition-state structure as from the reaction of trans-Pt-chloroaqua (or diaqua) to trans-monoadduct,which seems that the trans-Pt-chloroaqua (or diaqua) complex can generate trans-or cis-monoadduct via the same transition-state.     

Cu^2＋诱导甘氨酸质子迁移机理的理论研究

采用CCSD/6-31＋＋G＊＊//B3LYP/6-31＋＋G＊＊方法研究了Cu2＋诱导甘氨酸质子迁移的机理.优化得到了7个中性配合物和1个两性配合物,其中两性配合物最稳定,结合能为215.93 kcal mol-1.中性构型间通过分子内单键的旋转相互转化,C-N、C-C和C-O键旋转的能垒范围分别为1.62～2.49、0.27～7.80和2.27～16.97 kcal mol-1;中性构型N6经质子迁移变为两性构型,能垒为33.82 kcal mol-1.Cu2＋作用于甘氨酸,使甘氨酸N5原子负电荷减少超过0.5,降低了N5对H6原子的库仑吸引,钝化了共价键B（O3–H6）,动力学上不利于H6质子迁移;但是H6质子迁移后,形成的两性构型Z1却是热力学最稳定体系.     

Theoretical study of isomerization and dissociation of acetylene dication in the ground and excited electronic states

Ab initio calculations employing the configuration interaction method including Davidson's corrections for quadruple excitations have been carried out to unravel the dissociation mechanism of acetylene dication in various electronic states and to elucidate ultrafast acetylene-vinylidene isomerization recently observed experimentally. Both in the ground triplet and the lowest singlet electronic states of C2H2(2+) the proton migration barrier is shown to remain high, in the range of 50 kcal/mol. On the other hand, the barrier in the excited 2 3A" and 1 3A' states decreases to about 15 and 34 kcal/mol, respectively, indicating that the ultrafast proton migration is possible in these states, especially, in 2 3A", even at relatively low available vibrational energies. Rice-Ramsperger-Kassel-Marcus calculations of individual reaction-rate constants and product branching ratios indicate that if C2H(2)2+ dissociates from the ground triplet state, the major reaction products should be CCH+(3Sigma-)+H+ followed by CH+(3Pi)+CH+(1Sigma+) and with a minor contribution (approximately 1%) of C2H+(2A1)+C+(2P). In the lowest singlet state, C2H+(2A1)+C+(2P) are the major dissociation products at low available energies when the other channels are closed, whereas at Eint>5 eV, the CCH+(1A')+H+ products have the largest branching ratio, up to 70% and higher, that of CH+(1Sigma+)+CH+(1Sigma+) is in the range of 25%-27%, and the yield of C2H++C+ is only 2%-3%. The calculated product branching ratios at Eint approximately 17 eV are in qualitative agreement with the available experimental data. The appearance thresholds calculated for the CCH++H+, CH++CH+, and C2H++C+ products are 34.25, 35.12, and 34.55 eV. The results of calculations in the presence of strong electric field show that the field can make the vinylidene isomer unstable and the proton elimination spontaneous, but is unlikely to significantly reduce the barrier for the acetylene-vinylidene isomerization and to render the acetylene configuration unstable or metastable with respect to proton migration.     

The ground-state tunneling splitting of various carboxylic acid dimers

Carboxylic acid dimers in gas phase reveal ground-state tunneling splittings due to a double proton transfer between the two subunits. In this study we apply a recently developed accurate semiclassical method to determine the ground-state tunneling splittings of eight different carboxylic acid derivative dimers (formic acid, benzoic acid, carbamic acid, fluoro formic acid, carbonic acid, glyoxylic acid, acrylic acid, and N,N-dimethyl carbamic acid) and their fully deuterated analogs. The calculated splittings range from 5.3e-4 to 0.13 cm(-1) (for the deuterated species from 2.8e-7 to 3.3e-4 cm(-1)), thus indicating a strong substituent dependence of the splitting, which varies by more than two orders of magnitude. One reason for differences in the splittings could be addressed to different barriers heights, which vary from 6.3 to 8.8 kcal/mol, due to different mesomeric stabilization of the various transition states. The calculated splittings were compared to available experimental data and good agreement was found. A correlation could be found between the tunneling splitting and the energy barrier of the double proton transfer, as the splitting increases with increased strength of the hydrogen bonds. From this correlation an empirical formula was derived, which allows the prediction of the ground-state tunneling splitting of carboxylic acid dimers at a very low cost and the tunneling splittings for parahalogen substituted benzoic acid dimers is predicted.     

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.