Benchmark ab initio proton affinity and gas-phase basicity of α-alanine based on coupled-cluster theory and statistical mechanics

We determine the proton affinity (PA) and gas-phase basicity (GB) of amino acid α-alanine at a chemically accurate level by performing explicitly-correlated CCSD(T)-F12b/aug-cc-pVDZ geometry optimizations and normal mode vibrational frequency calculations as well as CCSD(T)-F12b/aug-cc-pVTZ energy computations at the possible neutral and protonated geometries. Temperature effects at 298.15 K considering translational, rotational, and vibrational enthalpy and entropy corrections are obtained via standard statistical mechanics utilizing the molecular geometries and the harmonic vibrational energy levels. Both the amino nitrogen (N) and the carbonyl oxygen (O) atoms are proven to be potential protonation sites and a systematic conformational search reveals 3 N- and 9 O-protonated conformers in the 0.00–7.88 and 25.43–30.43 kcal/mol energy ranges at 0 K, respectively. The final computed PA and GB values at (0)298.15 K in case of N-protonation are (214.47)216.80 and 207.07 kcal/mol, respectively, whereas the corresponding values for O-protonation are (189.04)190.63 and 182.31 kcal/mol. The results of the benchmark high-level coupled-cluster computations are utilized to assess the accuracy of several lower-level cost-effective methods such as MP2 and density functional theory with various functionals (SOGGA11-X, M06-2X, PBE0, B3LYP, M06, TPSS).     

Quantifying the role of water in protein-carbohydrate interactions

Water-mediated interactions play a key role in carbohydrate-lectin binding, where the interactions involve a conserved water that is separated from the bulk solvent and present a bridge between the side chains of the protein and the carbohydrate ligand. To apply quantum mechanical methods to examine the role of conserved waters, we present an analysis in which the relevant carbohydrate atoms are modeled by methanol, and in which the protein is replaced by a limited number of amino acid side chains. Clusters containing a conserved water and a representative amino acid fragment were also examined to determine the influence of amino acid side chains on interaction energies. To quantify the differential binding energies of methanol versus water, quantum mechanical calculations were performed at the B3LYP/6-311++G(3df,3pd)//B3LYP/6-31+G(d) level in which either a methanol molecule was bound to the conserved water (liganded state) or in which a water molecule replaces the methanol (unliganded state). Not surprisingly, the binding of a water to clusters containing charged amino acid side chains was more favorable by 1.55 to 7.23 kcal/mol than that for the binding of a water to the corresponding pure water clusters. In contrast, the binding energy of water to clusters containing polar-uncharged amino acid side chains ranged from 4.35 kcal/mol less favorable to 4.72 kcal/mol more favorable than for binding to the analogous pure water clusters. The overall trend for the binding of methanol versus water, in any of the clusters, favored methanol by an average value of 1.05 kcal/mol. To extend these studies to a complex between a protein (Concanavalin A) and its carbohydrate ligand, a cluster was examined that contained the side chains of three key amino acids, namely asparagine, aspartate, and arginine, as well as a key water molecule, arranged as in the X-ray diffraction structure of Con A. Again, using methanol as a model for the endogenous carbohydrate ligand, energies of -5.94 kcal/mol and -5.70 kcal/mol were obtained for the binding of methanol and water, respectively, to the Con A-water cluster. The extent to which cooperativity enhanced the binding energies has been quantified in terms of nonadditive three-body contributions. In general, the binding of water or methanol to neutral dimers formed cooperative clusters; in contrast, the cooperativity in charged clusters depended on the overall geometry as well as the charge.     

Model studies of the chemisorption of hydrogen and oxygen on nickel surfaces

Atomic chemisorption of hydrogen and oxygen on the Ni(100) surface has been studied using an Effective Core Potential (ECP) approach described in a previous paper. Clusters of up to 50 nickel atoms have been used to model the surface. The computed chemisorption energies are 62 kcal/mol (exp. 63 kcal/mol) for hydrogen and 106 kcal/mol (exp. 115–130 kcal/mol) for oxygen. Correlating the adsorbate and the cluster-adsorbate bonds is extremely important for obtaining accceptable results, particularly for oxygen. Reasonable convergence of chemisorption energies is obtained with 40–50 cluster atoms for both hydrogen and oxygen. For hydrogen the addition of a third cluster layer stabilizes the results considerably. Both hydrogen and oxygen are adsorbed at (or close to) the four-fold hollow site. The calculated barriers for surface migration are also in good agreement with the experimental estimates. The calculated equilibrium heights above the surface are on the other hand too high compared with experiments. This disagreement is believed to be due to core-valence correlation effects, which are not incorporated in the present ECP. The cluster convergence for the height above the surface is much slower than for the chemisorption energy.     

A novel quantum mechanical/molecular mechanical approach to the free energy calculation for isomerization of glycine in aqueous solution

The free energy change associated with the isomerization reaction of glycine in water solution has been studied by a hybrid quantum mechanical/molecular mechanical (QM/MM) approach combined with the theory of energy representation (QM/MM-ER) recently developed. The solvation free energies for both neutral and zwitterionic form of glycine have been determined by means of the QM/MM-ER simulation. The contributions of the electronic polarization and the fluctuation of the QM solute to the solvation free energy have been investigated. It has been found that the contribution of the density fluctuation of the zwitterionic solute is estimated as -4.2 kcal/mol in the total solvation free energy of -46.1 kcal/mol, while that of the neutral form is computed as -3.0 kcal/mol in the solvation free energy of -15.6 kcal/mol. The resultant free energy change associated with the isomerization of glycine in water has been obtained as -7.8 kcal/mol, in excellent agreement with the experimental data of -7.3 or -7.7 kcal/mol, implying the accuracy of the QM/MM-ER approach. The results have also been compared with those computed by other methodologies such as the polarizable continuum model and the classical molecular simulation. The efficiency and advantage of the QM/MM-ER method has been discussed.     

Comparison of 6-31G*-based MST/SCRF and FEP evaluations of the free energies of hydration for small neutral molecules

A comparison between Miertus–Scrocco–Tomasi (MST) SCRF and free energy perturbation (FEP) estimates of the free energy of hydration of eight small neutral molecules is presented. In both cases, the 6-31G* molecular electrostatic potential is used to describe the electrostatic properties of the molecules. The results demonstrate the ability of both methodologies to provide useful theoretical estimates of the total free energy of hydration; the average errors are only 1.5 kcal/mol (FEP) and 0.8 kcal/mol (MST/SCRF). The largest errors in the FEP and MST/SCRF results are less than 1.5 kcal/mol for all molecules except acetic acid, where the FEP method overestimates the free energy of hydration by 3.3 kcal/mol. © John Wiley & Sons, Inc.     

A theoretical investigation of the relative stability of hydrated glycine and methylcarbamic acid--from water clusters to interstellar ices

We have theoretically investigated how the low-energy conformers of the neutral and the zwitterionic forms of glycine as well as methylcarbamic acid are stabilized by the presence water. The MP2/6-311++G(d,p) method was utilized to conduct calculations on glycine and methylcarbamic acid in both isolated clusters and in clusters embedded in the conductor-like polarizable continuum model (C-PCM), where the clusters explicitly contain between one and ten water molecules. The neutral forms of glycine and methylcarbamic acid were found to have similar hydration energies, whereas the neutral methylcarbamic acid was determined to be approximately 32 kJ mol(-1) more stable than the neutral glycine in the isolated clusters and 30 kJ mol(-1) more stable in the C-PCM embedded clusters. Both the number and strength of the hydrogen bonding interactions between water and the zwitterions drive the stability. This lowers the relative energy of the glycine zwitterion from 50 kJ mol(-1) above neutral glycine, when there are two water molecules in the clusters to 11 kJ mol(-1) below for the clusters containing ten water molecules. For the methylcarbamic acid clusters with two water molecules, the zwitterion is 51 kJ mol(-1) higher in energy than the neutral form, but it remains 13 kJ mol(-1) above the neutral methylcarbamic acid in the clusters containing ten water molecules. When the bulk water environment is simulated by the C-PCM calculations, we find both the methylcarbamic acid and glycine zwitterionic forms have similar energies at 20 kJ mol(-1) above the neutral methylcarbamic acid energy and 10 kJ mol(-1) lower than the neutral glycine energy. Although neither methylcarbamic acid nor glycine have been detected in the interstellar medium yet, our findings indicate that methylcarbamic acid is the more stable product from methylamine and carbon dioxide reactions in a water ice. This suggests that methylcarbamic acid likely plays a role in the intermediate steps if glycine is formed in the interstellar medium.     

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.     

Low-energy tautomers and conformers of neutral and protonated arginine.

The relative stabilities of zwitterionic and canonical forms of neutral arginine and of its protonated derivative were studied by using ab initio electronic structure methods. Trial structures were first identified at the PM3 level of theory with use of a genetic algorithm to systematically vary geometrical parameters. Further geometry optimizations of these structures were performed at the MP2 and B3LYP levels of theory with basis sets of the 6-31++G** quality. The final energies were determined at the CCSD/6-31++G** level and corrected for thermal effects determined at the B3LYP level. Two new nonzwitterionic structures of the neutral were identified, and one of them is the lowest energy structure found so far. The five lowest energy structures of neutral arginine are all nonzwitterionic in nature and are clustered within a narrow energy range of 2.3 kcal/mol. The lowest energy zwitterion structure is less stable than the lowest nonzwitterion structure by 4.0 kcal/mol. For no level of theory is a zwitterion structure suggested to be the global minimum. The calculated proton affinity of 256.3 kcal/mol and gas-phase basicity of 247.8 kcal/mol of arginine are in reasonable agreement with the measured values of 251.2 and 240.6 kcal/mol, respectively. The calculated vibrational characteristics of the low-energy structures of neutral arginine provide an alternative interpretation of the IR-CRLAS spectrum (Chapo et al. J. Am. Chem. Soc. 1998, 120, 12956-12957).     

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.     

Synergy of intramolecular hydrogen-bonding network in myo-inositol 2-monophosphate: theoretical investigations into the electronic structure, proton transfer, and pKa

This work demonstrates the pivotal role that an intramolecular hydrogen-bonding network (intra-HBN) plays in the determination of the conformation of myo-inositol 2-monophosphate (Ins(2)P1), a member of the inositol phosphate family of compounds, which are important participants in the role that phosphates play in biological and environmental chemistry. For biologically significant compounds that contain phosphate and hydroxyl groups, Ins(2)P1 is a model system for studying both the primary forces that determine their conformations and their chemical properties from the effect of phosphate group addition. We performed ab initio calculations to determine the intra-HBN within important thermally accessible conformations for neutral Ins(2)P1 and its anions, Ins(2)P1(1-) and Ins(2)P1(2-). The results show that the global minima prefer 1a/5e structures where the phosphate group is in the axial position with all -OH groups in the equatorial positions. The calculations of transition state structures for ring inversion at each ionization state predict an activation energy of 18.16 kcal/mol for the neutral species in water, while the activation energy is lower for the charged compounds, 15.62 kcal/mol for Ins(2)P1(1-) and 12.48 kcal/mol for Ins(2)P1(2-). The pK(a) values of Ins(2)P1 were calculated by modeling the solvent as a polarizable continuum medium (PCM) and as explicit solvent molecules. These values are in good agreement with experimental data. A novel four-center pattern of hydrogen bonding was found to stabilize the system. The intramolecular proton transfer across a low barrier hydrogen bond between the charged phosphate and hydroxyl groups was found to occur under standard conditions with an activation energy that is less than 0.5 kcal/mol.     

Vinylogous Mannich reactions: some theoretical studies on the origins of diastereoselectivity

[reaction: see text] Ab initio calculations at the RHF/3-21G level were used to investigate the limiting transition states in the addition of 2-methoxyfuran to a pyrrolinium ion. Four stationary points were found on the potential energy surface with relative energies of threo Diels-Alder, 0.0 kcal/mol; erythro open, 0.9 kcal/mol; erythro Diels-Alder, 1.3 kcal/mol; and threo open, 1.8 kcal/mol.     

Binding energies of protonated betaine complexes: a probe of zwitterion structure in the gas phase

The dissociation kinetics of proton-bound dimers of betaine with molecules of comparable gas-phase basicity were investigated using blackbody infrared radiative dissociation (BIRD). Threshold dissociation energies were obtained from these data using master equation modeling. For bases that have comparable or higher gas-phase basicity, the binding energy of the protonated base.betaine complex is ~1.4 eV. For molecules that are ~2 kcal/mol or more less basic, the dissociation energy of the complexes is ~1.2 eV. The higher binding energy of the former is attributed to an ion-zwitterion structure which has a much larger ion-dipole interaction. The lower binding energy for molecules that are ~2 kcal/mol or more less basic indicates that an ion-molecule structure is more favored. Semiempirical calculations at both the AM1 and PM3 levels indicate the most stable ion-molecule structure is one in which the base interacts with the charged quaternary ammonium end of betaine. These results indicate that the measurement of binding energies of neutral molecules to biological ions could provide a useful probe for the presence of zwitterions and salt bridges in the gas phase. From the BIRD data, the gas-phase basicity of betaine obtained from the kinetic method is found to be 239.2 +/- 1.0 kcal/mol. This value is in excellent agreement with the value of 239.3 kcal/mol (298 K) from ab initio calculations at the MP2/6-31+g** level. The measured value is slightly higher than those reported previously. This difference is attributed to entropy effects. The lower ion internal energy and longer time frame of BIRD experiments should provide values closer to those at standard temperature.     

Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains

Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.     

A new alligator-clip compound for molecular electronics

We used the B3LYP flavor of density functional calculations to study new alligator-clip compounds for molecular electronics with platinum electrodes. First, with commonly used S-based linkage molecule 3-methyl-1,2-dithiolane (MDTL) we find that after chemisorption on Pt(1 1 1) the most stable structure is ring-opened with a binding energy of 32.44 kcal/mol. Among several alternative alligator-clip compounds we find that P-based molecules lead to much higher binding energies. For the ring-closed structure of 3-methyl-1,2-diphospholane (MDPL) a binding energy of 47.72 kcal/mol was calculated and even 54.88 kcal/mol for the ring-opened molecule. Thus, MDPL provides a more stable link to the metal surface and might increase the conductance.     

Probing phenylalanine/adenine pi-stacking interactions in protein complexes with explicitly correlated and CCSD(T) computations

To examine the effects of pi-stacking interactions between aromatic amino acid side chains and adenine bearing ligands in crystalline protein structures, 26 toluene/(N9-methyl)adenine model configurations have been constructed from protein/ligand crystal structures. Full geometry optimizations with the MP2 method cause the 26 crystal structures to collapse to six unique structures. The complete basis set (CBS) limit of the CCSD(T) interaction energies has been determined for all 32 structures by combining explicitly correlated MP2-R12 computations with a correction for higher-order correlation effects from CCSD(T) calculations. The CCSD(T) CBS limit interaction energies of the 26 crystal structures range from -3.19 to -6.77 kcal mol (-1) and average -5.01 kcal mol (-1). The CCSD(T) CBS limit interaction energies of the optimized complexes increase by roughly 1.5 kcal mol (-1) on average to -6.54 kcal mol (-1) (ranging from -5.93 to -7.05 kcal mol (-1)). Corrections for higher-order correlation effects are extremely important for both sets of structures and are responsible for the modest increase in the interaction energy after optimization. The MP2 method overbinds the crystal structures by 2.31 kcal mol (-1) on average compared to 4.50 kcal mol (-1) for the optimized structures.     

Mechanism of C-terminal intein cleavage in protein splicing from QM/MM molecular dynamics simulations

Protein splicing is a post-translational process in which a biologically inactive protein is activated by the release of a segment denoted as an intein. The process involves four steps. In the third, the scission of the intein takes place after the cyclization of the last amino acid of the segment, an asparagine. Little is known about the chemical reaction necessary for this cyclization. Experiments demonstrate that two histidines (the penultimate amino acid of the intein, and a histidine located 10 amino acids upstream) are relevant in the cyclization of the asparagine. We have investigated the mechanism and determinants of reaction in the GyrA intein focusing on the requirements for asparagine activation for its cyclization. First, the influence that the protonation states of these two histidines have on the orientation of the asparagine side chain is investigated by means of molecular dynamics simulation. Molecular dynamics simulations using the CHARMM27 force field were carried out on the three possible protonation states for each of these two histidines. The results indicate that the only protonation state in which the conformation of the system is suitable for cyclization is when the penultimate histidine is fully protonated (positively charged), and the upstream histidine is in the His(ε) neutral tautomeric form. The free energy profile for the reaction in which the asparagine is activated by a proton transfer to the upstream histidine is presented, computed by hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics at the SCCDFTB/CHARMM27 level of theory. The calculated free energy barrier for the reaction is 19.0 kcal mol(-1). B3LYP/6-31+G(d) QM/MM single-point calculations give a qualitatively a similar energy profile, although with somewhat higher energy barriers, in good agreement with the value derived from experiment of 25 kcal mol(-1) at 60 °C. QM/MM molecular dynamics simulations of the reactant, activated reactant and intermediate states highlight the importance of the Arg181-Val182-Asp183 segment in catalysing the reaction. Overall, the results indicate that nucleophilic activation of the asparagine for its cyclization by the upstream histidine acting as the base is a plausible mechanism for the C-terminal cleavage in protein splicing.     

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.     

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

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

Magnitude of the CH/pi interaction in the gas phase: experimental and theoretical determination of the accurate interaction energy in benzene-methane

The accurate CH/pi interaction energy of the benzene-methane model system was experimentally and theoretically determined. In the experiment, mass analyzed threshold ionization spectroscopy was applied to the benzene-methane cluster in the gas phase, prepared in a supersonic molecular beam. The binding energy in the neutral ground state of the cluster, which is regarded as the CH/pi interaction energy for this model system, was evaluated from the dissociation threshold measurements of the cluster cation. The experimentally determined binding energy (D(0)) was 1.03-1.13 kcal/mol. The interaction energy of the model system was calculated by ab initio molecular orbital methods. The estimated CCSD(T) interaction energy at the basis set limit (D(e)) was -1.43 kcal/mol. The calculated binding energy (D(0)) after the vibrational zero-point energy correction (1.13 kcal/mol) agrees well with the experimental value. The effects of basis set and electron correlation correction procedure on the calculated CH/pi interaction energy were evaluated. Accuracy of the calculated interaction energies by DFT methods using BLYP, B3LYP, PW91 and PBE functionals was also discussed.