Which is more likely in trichodiene biosynthesis: hydride or proton transfer?

The mechanisms proposed for enzyme-catalyzed formation of the sesquiterpene natural product trichodiene consistently include a step involving a 1,4-hydride transfer. Using quantum chemical methods (B3LYP/6-31+G(d,p) and mPW1PW91/6-31+G(d,p)), we discovered two alternative pathways for transformation of the intermediate bisabolyl cation to the cuprenyl cation, one of which--a proton-transfer pathway--appears to be much more energetically favorable (by more than 10 kcal/mol) than the hydride transfer pathways usually proposed.     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.     

Comparison of Some Computational Methods for Geometry Optimisation of Phosphorus Acid Derivatives

Several economical methods for geometry optimization, that should be applicable to larger molecules, have been evaluated for 19 phosphorus acid derivatives. MP2/cc-pVDZ geometry optimizations are used as reference points and the geometries obtained from the other methods are evaluated with respect to deviations in bond lengths and angles, from the reference geometries. The geometry optimization methods are also compared to the much used B3LYP/6-31G(d) method. Single point energies obtained by subsequent EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) calculations on the respective equilibrium geometries are also reported relative to the energies obtained from the reference geometries. The geometries from HF/MIDI! optimizations were closer to those of the references than the geometries of the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! optimizations. The EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) single point energies obtained from the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! geometries gave a mean absolute deviation (MAD) from that of the reference geometries of 1.4-3.9 kcal mol m 1 . The HF/MIDI! geometries, however, gave EDF1/6-31+G(d) and B3LYP/6-31+G(d,p) energies with a MAD of only about 0.5 and 0.55 kcal mol m 1 respectively from the energies obtained with the reference geometries. Thus, use of HF/MIDI! for geometry optimization of phosphorus acids is a method that gives geometries of near-MP2 quality, resulting in a fair accuracy of energies in subsequent single point calculations, at a much lower computational cost other methods that give similar accuracies.     

Computational study on the mechanisms and energetics of trimethylindium reactions with H2O and H2S

The reactions of trimethylindium (TMIn) with H2O and H2S are relevant to the chemical vapor deposition of indium oxide and indium sulfide thin films. The mechanisms and energetics of these reactions in the gas phase have been investigated by density functional theory and ab initio calculations using the CCSD(T)/[6-31G(d,p)+Lanl2dz]//B3LYP/[6-31G(d,p)+Lanl2dz] and CCSD(T)/[6-31G(d,p)+Lanl2dz] //MP2/[6-31G(d,p)+Lanl2dz] methods. The results of both methods are in good agreement for the optimized geometries and relative energies. When TMIn reacts with H2O and H2S, initial molecular complexes [(CH3)3In:OH2 (R1)] and [(CH3)3In:SH2 (R2)] are formed with 12.6 and 3.9 kcal/mol binding energies. Elimination of a CH4 molecule from each complex occurs with a similar energy barrier at TS1 (19.9 kcal/mol) and at TS3 (22.1 kcal/mol), respectively, giving stable intermediates (CH3)2InOH and (CH3)2InSH. The elimination of the second CH4 molecule from these intermediate products, however, has to overcome very high and much different barriers of 66.1 and 53.2 kcal/mol, respectively. In the case of DMIn with H2O and H2S reactions, formation of both InO and InS is exothermic by 3.1 and 30.8 kcal/mol respectively. On the basis of the predicted heats of formation of R1 and R2 at 0 K and -20.1 and 43.6 kcal/mol, the heats of formation of (CH3)2InOH, (CH3)2InSH, CH3InO, CH3InS, InO, and InS are estimated to be -20.6, 31.8, and 29.0 and 48.4, 35.5, and 58.5 kcal/mol, respectively. The values for InO and InS are in good agreement with available experimental data. A similar study on the reactions of (CH3)2In with H2O and H2S has been carried out; in these reactions CH3InOH and CH3InSH were found to be the key intermediate products.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

Investigation of ethylene dimerization over faujasite zeolite by the ONIOM method.

Ethylene dimerization was investigated by using an 84T cluster of faujasite zeolite modeled by the ONIOM3(MP2/6-311++G(d,p):HF/6-31G(d):UFF) method. Concerted and stepwise mechanisms were evaluated. In the stepwise mechanism, the reaction proceeds by protonation of ethylene to form the surface ethoxide and then C--C bond formation between the ethoxide and the second ethylene molecule to give the butoxide product. The first step is rate-determining and has an activation barrier of 30.06 kcal mol(-1). The ethoxide intermediate is rather reactive and readily reacts with another ethylene molecule with a smaller activation energy of 28.87 kcal mol(-1). In the concerted mechanism, the reaction occurs in one step of simultaneous protonation and C--C bond formation. The activation barrier is calculated to be 38.08 kcal mol(-1). Therefore, the stepwise mechanism should dominate in ethylene dimerization.     

Structures, internal rotor potentials, and thermochemical properties for a series of nitrocarbonyls, nitroolefins, corresponding nitrites, and their carbon centered radicals

Structures, enthalpy (Δ(f)H°(298)), entropy (S°(T)), and heat capacity (C(p)(T)) are determined for a series of nitrocarbonyls, nitroolefins, corresponding nitrites, and their carbon centered radicals using the density functional B3LYP and composite CBS-QB3 calculations. Enthalpies of formation (Δ(f)H°(298)) are determined at the B3LYP/6-31G(d,p), B3LYP/6-31+G(2d,2p), and composite CBS-QB3 levels using several work reactions for each species. Entropy (S) and heat capacity (C(p)(T)) values from vibration, translational, and external rotational contributions are calculated using the rigid-rotor-harmonic-oscillator approximation based on the vibration frequencies and structures obtained from the density functional studies. Contribution to Δ(f)H(T), S, and C(p)(T) from the analysis on the internal rotors is included. Recommended values for enthalpies of formation of the most stable conformers of nitroacetone cc(═o)cno2, acetonitrite cc(═o)ono, nitroacetate cc(═o)no2, and acetyl nitrite cc(═o)ono are -51.6 kcal mol(-1), -51.3 kcal mol(-1), -45.4 kcal mol(-1), and -58.2 kcal mol(-1), respectively. The calculated Δ(f)H°(298) for nitroethylene c═cno2 is 7.6 kcal mol(-1) and for vinyl nitrite c═cono is 7.2 kcal mol(-1). We also found an unusual phenomena: an intramolecular transfer reaction (isomerization) with a low barrier (3.6 kcal mol(-1)) in the acetyl nitrite. The NO of the nitrite (R-ONO) in CH(3)C(═O')ONO moves to the C═O' oxygen in a motion of a stretching frequency and then a shift to the carbonyl oxygen (marked as O' for illustration purposes).     

Theoretical study of the addition and abstraction reactions of hydroxyl radical with uracil

The addition as well as abstraction reactions of hydroxyl radical (OH) with the nucleic acid base, uracil (U), in the gas phase has been explored at the B3LYP/6-31+G(d,p) level of density functional theory (DFT). The energy barrier of the OH addition to both the C5 and C6 positions of the uracil is less than 1 kcal/mol while the hydrogen abstractions (H-abstractions) from either the N1 or the N3 positions are ∼9.5 kcal/mol. Further the energetics of these reactions are assessed by applying the effect of aqueous medium through the polarizable continuum model (PCM). Both the gas and the solution phase data established that the thermodynamic and kinetic factors are more favorable for the OH addition to either C5 or C6 positions of the uracil than the H-abstraction reactions. Moreover, calculations at the MPW1K/6-31+G(d,p), CCSD(T)/6-31+G(d,p)//B3LYP/6-31+G(d,p) and CCSD(T)/6-31+G(d,p)//MPW1K/6-31+G(d,p) levels of theoretical methods qualitatively supported the B3LYP/6-31+G(d,p) results.     

A theoretical study of the mechanisms and regiochemistry of the reactions of 5-alkoxyoxazole with thioaldehydes,nitroso compounds,and aldehydes

A theoretical study based on B3LYP/6-31G calculations has been applied to the mechanisms and regiochemistry of reactions of 5-alkoxyoxazole with thioaldehydes, nitroso compounds, and aldehydes. All three reactions adopt similar mechanisms, which start with Diels-Alder (DA) reactions, followed by either a novel, concerted ring-opening-ring-closing (RORC) step to transfer the DA adduct to 2-alkoxycarbonyl-3-thiazoline and 2-alkoxycarbonyl-3-oxazoline for thioaldehydes and aldehydes, respectively, or stepwise ring-opening and ring-closing steps to generate 1,2,4-oxadiazoline for nitroso compounds. The reactions of 5-alkoxyoxazole with thioaldehydes and nitroso compounds can be conducted under thermal reaction conditions due to the 10 kcal/mol activation barriers for their rate-determining DA reactions. By contrast, the reaction of 5-alkoxyoxazole with aldehydes cannot take place under thermal conditions, since this bimolecular reaction has the rate-determining RORC transition state higher than the reactants by 30.5 kcal/mol.     

Formation of 2-hexene by cationic dimerization of propene: an ab initio and density functional theory study

The formation of a 2-hexene radical cation from a propene radical cation and a neutral propene molecule is investigated by means of ab initio UHF and spin projected MP2 calculations, as well as the SVWN and B3LYP levels of density functional theory. A stable addition complex, with loose CC bonds, is found. To proceed from the addition complex to the product, a locally planar transition state must be passed, with a migrating hydrogen located half-way between the donating and the accepting carbon atoms. At the highest computational levels considered, PMP2/6-31G(d,p)//MP2/3-21G and B3LYP/6-31G(d,p), this transition state lies approximately 11 and 13 kcal/mol, respectively, above the addition complex. The high barrier is believed to be one reason why radical cation oligomerization of propene has not been detected experimentally, in contrast to the case of ethene, where the corresponding barrier is only a few kcal/mol. Received: 17 December 1996 / Accepted: 12 May 1997     

A theoretical study of the effect of a tetraalkylammonium counterion on the hydrogen bond strength in Z-hydrogen maleate.

High-level ab initio calculations (B3LYP/6-31+G and QCISD(T)/6-311+G**) were carried out to resolve the disagreement between recent experimental and computational estimates of the relative strength of the intramolecular hydrogen bond in Z-hydrogen maleate anion with respect to the normal hydrogen bond in maleic acid. The computational estimates for the strength of the intramolecular hydrogen bond in the gas-phase maleate anion are in a range of 14-28 kcal/mol depending on the choice of the reference structure. Computational data suggest that the electrostatic influence of a counterion such as a tetraalkylammonium cation can considerably weaken the hydrogen bonding interaction (by 1.5-2 times) in the complexed hydrogen maleate anion relative to that in the naked anion. The estimated internal H-bonding energies for a series of Z-maleate/R4N+ salts (R = CH3, C2H5, CH3CH2CH2CH2) range from 8 to 13 kcal/mol. The calculated energy differences between the E- and Z-hydrogen maleates complexed to Me4N+, Et4N+, and Bu4N+ cation are 4.9 (B3LYP/6-31+G(d,p)) and 5.7 and 5.8 kcal/mol (B3LYP/6-31G(d)). It is also demonstrated that the sodium cation exerts a similar electrostatic influence on the hydrogen bond strength in bifluoride anion (FHF-). The present study shows that while low-barrier short hydrogen bonds can exist in the gas phase (the barrier for the hydrogen transfer in maleate anion is only 0.2 kcal/mol at the QCISD(T)/6-311+G//QCISD/6-31+G level), whether they can also be strong in condensed media or not depends on how their interactions with their immediate environment affect their strength.     

Quadricyclane radical cation rearrangements: a computational study of the transformations to 1,3,5-cycloheptatriene and norbornadiene

An alternative skeletal rearrangement of the quadricyclane radical cation (Q*+) explains the side products formed in the one-electron oxidation to norbornadiene. First, the bicyclo[2.2.1]hepta-2-ene-5-yl-7-ylium radical cation, with an activation energy of 14.9 kcal mol(-1), is formed. Second, this species can further rearrange to 1,3,5-cycloheptatriene through two plausible paths, that is, a multistep mechanism with two shallow intermediates and a stepwise path in which the bicyclo[3.2.0]hepta-2,6-diene radical cation is an intermediate. The multistep rearrangement has a rate-limiting step with an estimated activation energy of 16.5 kcal mol(-1), which is 2.8 kcal mol(-1) lower in energy than the stepwise mechanism. However, the lowest activation energy is found for the Q*+ cycloreversion to norbornadiene that has a transition structure, in close correspondence with earlier studies, and an activation energy of 10.1 kcal mol(-1), which agrees well with the experimental estimate of 9.3 kcal mol(-1). The computational estimates of activation energies were done using the CCSD(T)/6-311+G(d,p) method with geometries optimized on the B3LYP/6-311+G(d,p) level, combined with B3LYP/6-311+G(d,p) frequencies.     

Quantum mechanical study of the ring-closing reaction of the hexatriene radical cation

The ring-closing reaction of hexatriene radical cation 1(*)(+) to 1,3-cyclohexadiene radical cation 2(*)(+) was studied computationally at the B3LYP/6-31G* and QCISD(T)/6-311G*//QCISD/6-31G* levels of theory. Both, concerted and stepwise mechanisms were initially considered for this reaction. Upon evaluation at the B3LYP level of theory, three of the possible pathways-a concerted C(2)-symmetric via transition structure 3(*)(+) and stepwise C(1)-symmetric pathways involving three-membered ring intermediate 5(*)(+) and four-membered ring intermediate 6(*)(+)-were rejected due to high-energy stationary points along the reaction pathway. The two remaining pathways were found to be of competing energy. The first proceeds through the asymmetric, concerted transition structure 4(*)(+) with an activation barrier E(a) = 16.2 kcal/mol and an overall exothermicity of -23.8 kcal/mol. The second pathway, beginning from the cis,cis,trans rotamer of 1(*)(+), proceeds by a stepwise pathway to the cyclohexadiene product with an overall exothermicity of -18.6 kcal/mol. The activation energy for the rate-determining step in this process, the formation of the intermediate bicyclo[3.1.0]hex-2-ene via transition structure 9(*)(+), was found to be 20.4 kcal/mol. More rigorous calculations of a smaller subsection of the potential energy hypersurface at the QCISD(T)//QCISD level confirmed these findings and emphasized the importance of conformational control of the reactant.     

Reactions of 1,3-cyclohexadiene with singlet oxygen. A theoretical study.

A thorough study of the reaction of singlet oxygen with 1,3-cyclohexadiene has been made at the B3LYP/6-31G(d) and CASPT2(12e,10o) levels. The initial addition reaction follows a stepwise diradical pathway to form cyclohexadiene endoperoxide with an activation barrier of 6.5 kcal/mol (standard level = CASPT2(12e,10o)/6-31G(d); geometries and zero-point corrections at B3LYP/6-31G(d)), which is consistent with an experimental value of 5.5 kcal/mol. However, as the enthalpy of the transition structure for the second step is lower than the diradical intermediate, the reaction might also be viewed as a nonsynchronous concerted reaction. In fact, the concertedness of the reaction is temperature dependent since entropy differences create a free energy barrier for the second step of 1.8 kcal/mol at 298 K. There are two ene reactions; one is a concerted mechanism (DeltaH(double dagger) = 8.8 kcal/mol) to 1-hydroperoxy-2,5-cyclohexadiene (5), while the other, which forms 1-hydroperoxy-2,4-cyclohexadiene (18), passes through the same diradical intermediate (9) as found on the pathway to endoperoxide. The major pathway from the endoperoxide is O-O bond cleavage (22.0 kcal/mol barrier) to form a 1,4-diradical (25), which is 13.9 kcal/mol less stable than the endoperoxide. From the diradical, two low-energy pathways exist, one to epoxyketone (29) and the other to the diepoxide (27), where both products are known to be formed experimentally with a product ratio sensitive to the nature of substitutents. A significantly higher activation barrier leads to C-C bond cleavage and direct formation of maleic aldehyde plus ethylene.     

Homoaromaticity in carbene intermediates

Analyses of isodesmic reactions for singlet 7-carbenanorbornene (4S), 8-carbena-endo-tricyclo[3.2.1.0(2,4)]octane (5S), 3-carbenabicyclo[3.1.0]hexane (3S), 2-carbenanorbornene 6S, and 2-carbenabicyclooctadiene 7S at the B3LYP/6-311+G(3df, 2p)//B3LYP/6-31G level provide stabilization energies of 13.83, 13. 50, 3.00, -2.22, and -3.01 kcal/mol, respectively. The C7 carbene in 4S and the C8 carbene center in 5S are strongly bent toward the double bond and cyclopropane ring, respectively, in contrast to their related triplets, 4T and 5T and parent hydrocarbons. The geometric change for 3S compared to 3T or to parent bicyclo[3.1. 0]hexane is minimal. Comparison of the stability of 6S with 2-carbenanorbornane and the geometry of 6S with that of 6T and also with the singlet and triplet 2-carbenanorbornane suggests very modest bridging. The stabilization energy and geometry of 2-carbenabicyclooctadiene 7 resemble antihomoaromatic bicyclooctadienyl cation 9 rather than the related homoaromatic bicyclooctadienyl anion 8. The diamagnetic susceptibility exhaltations (Deltachi) of 3, 4, 5, and 7, calculated at the B3LYP/6-311+G(2d,p)//B3LYP/6-31G level, are -0.7, 22.7, 26.0, and -10.3 cgs-ppm, respectively. The singlet-triplet energy differences, DeltaE(TS), for carbenacyclohexane, carbenacyclopentane, 3, 4, 5, 6, 7, and 2-carbenabicyclo[3.2.1]oct-3-ene at the B3LYP/6-311+G(3df, 2p)//B3LYP/6-31G level are 3.4, 10.3, 8.8, 27.1, 25.9, 12.7, -4.0, and -0.9 kcal/mol, respectively.     

Hydride affinity scale of various substituted arylcarbeniums in acetonitrile

Combined with the integral equation formalism polarized continuum model (IEFPCM), the hydride affinities of 96 various acylcarbenium ions in the gas phase and CH(3)CN were estimated by using the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-311++G(2df,2p)//B3LYP/6-31+G(d), and BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) methods for the first time. The results show that the combination of the BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) method and IEFPCM could successfully predict the hydride affinities of arylcarbeniums in MeCN with a precision of about 3 kcal/mol. On the basis of the calculated results from the BLYP method, it can be found that the hydride affinity scale of the 96 arylcarbeniums in MeCN ranges from -130.76 kcal/mol for NO(2)-PhCH(+)-CN to -63.02 kcal/mol for p-(Me)(2)N-PhCH(+)-N(Me)(2), suggesting most of the arylcarbeniums are good hydride acceptors. Examination of the effect of the number of phenyl rings attached to the carbeniums on the hydride affinities shows that the increase of the hydride affinities takes place linearly with increasing number of benzene rings in the arylcarbeniums. Analyzing the effect of the substituents on the hydride affinities of arylcarbeniums indicates that electron-donating groups decrease the hydride affinities and electron-withdrawing groups show the opposite effect. The hydride affinities of arylcarbeniums are linearly dependent on the sum of the Hammett substituent parameters σ(p)(+). Inspection of the correlation of the solution-phase hydride affinities with gas-phase hydride affinities and aqueous-phase pK(R)(+) values reveals a remarkably good correspondence of ΔG(H(-)A)(R(+)) with both the gas-phase relative hydride affinities only if the α substituents X have no large electron-donating or -withdrawing properties and the pK(R)(+) values even though the media are dramatically different. The solution-phase hydride affinities also have a linear relationship with the electrophilicity parameter E, and this dependence can certainly serve as one of the most effective ways to estimate the new E values from ΔG(H(-)A)(R(+)) or vice versa. Combining the hydride affinities and the reduction potentials of the arylcarbeniums, we obtained the bond homolytic dissociation Gibbs free energy changes of the C-H bonds in the corresponding hydride adducts in acetonitrile, ΔG(HD)(RH), and found that the effects of the substituent on ΔG(HD)(RH) are very small. Simple thermodynamic analytic platforms for the three C-H cleavage modes were constructed. It is evident that the present work would be helpful in understanding the nature of the stabilities of the carbeniums and mechanisms of the hydride transfers between carbeniums and other hydride donors.     

Ab initio and experimental studies on the hetero-Diels-Alder and cheletropic additions of sulfur dioxide to (E)-1-methoxybutadiene: a mechanism involving three molecules of SO(2)

Kinetics on the cheletropic addition of sulfur dioxide to (E)-1-methoxybutadiene (1) to give the corresponding sulfolene 2 (2-methoxy-2,5-dihydrothiophene-1,1-dioxide) gave the rate law d[2]/dt = k[1][SO(2)](x)() with x = 2.6 +/- 0.2 at 198 K. Under these conditions, no sultine 3 [(2RS,6RS)-6-methoxy-3,6-dihydro-1,2-oxathiin-2-oxide] resulting from a hetero-Diels-Alder addition was observed, and the cheletropic elimination 2 --> 1 + SO(2) did not occur. Ab initio and DFT quantum calculations confirmed that the cheletropic addition 1 + SO(2) --> 2 follows two parallel mechanisms, one involving two molecules of SO(2) and the transition structure with DeltaG(++) = 18.2 +/- 0.2 kcal/mol at 198 K (exptl); 22.5-22.7 kcal/mol [B3LYP/6-31G(d,p)], the other one involving three molecules of SO(2) with DeltaG(++) = 18.9 +/- 0.1 kcal/mol at 198 K (exptl); 19.7 kcal/mol [B3LYP/6-31G(d,p)]. The mechanism involving only one molecule of SO(2) in the transition structure requires a higher activation energy, DeltaG(++) = 25.2 kcal/mol [B3LYP/6-31G(d,p)]. Comparison of the geometries and energetics of the structures involved into the 1 + SO(2) --> 2, 3 and 1 + 2SO(2) --> 2, 3 + SO(2) reactions obtained by ab initio and DFT methods suggest that the latter calculation techniques can be used to study the cycloadditions of sulfur dioxide. The calculations predict that the hetero-Diels-Alder addition 1 + SO(2) --> 3 also prefers a mechanism in which three molecules of SO(2) are involved in the cycloaddition transition structure. At 198 K and in SO(2) solutions, the entropy cost (TDeltaS(++)) is overcompensated by the specific solvation by SO(2) in the transition structures of both the cheletropic and hetero-Diels-Alder reactions of (E)-1-methoxybutadiene with SO(2).     

Mechanism of formation of the internal aldimine in pyridoxal 5'-phosphate-dependent enzymes

In this paper we studied the mechanism of formation of the internal aldimine, a common intermediate to most pyridoxal 5'-phosphate (PLP)-dependent enzymes. A large model based on the crystal structure from the human ornithine decarboxylase (ODC) enzyme was constructed and in total accounts for 504 atoms. The reaction mechanism was investigated using the ONIOM methodology (B3LYP/6-31G(d)//AM1), and the final energies were calculated with the M06/6-311++G(2d,2p)//B3LYP/6-31G(d) level of theory. It was demonstrated that the reaction is accomplished in three sequential steps: (i) the nucleophilic attack of Lysine69 to PLP, (ii) the carbinolamine formation, and (iii) a final dehydration step. For the carbinolamine formation, several mechanistic hypotheses were explored, and the preferred pathway assigns a key role for the conserved active site Cys360. The overall reaction is exergonic in -9.1 kcal/mol, and the rate-limiting step is the dehydration step (E(a) = 13.5 kcal/mol). For the first time, we provide an atomistic portrait of this mechanism in an enzymatic environment. Moreover, we were able to assign a novel role to Cys360 in the ODC reaction mechanism that was never proposed.     

o-Quinone methide as alkylating agent of nitrogen, oxygen, and sulfur nucleophiles. The role of H-bonding and solvent effects on the reactivity through a DFT computational study

The reactivity of the alkylating agent o-quinone methide (o-QM) toward NH(3), H(2)O, and H(2)S, prototypes of nitrogen-, oxygen-, and sulfur-centered nucleophiles, has been studied by quantum chemical methods in the frame of DF theory (B3LYP) in reactions modeling its reactivity in water with biological nucleophiles. The computational analysis explores the reaction of NH(3), H(2)O, and H(2)S with o-QM, both free and H-bonded to a discrete water molecule, with the aim to rationalize the specific and general effect of the solvent on o-QM reactivity. Optimizations of stationary points were done at the B3LYP level using several basis sets [6-31G(d), 6-311+G(d,p), adding d and f functions to the S atom, 6-311+G(d,p),S(2df), and AUG-cc-pVTZ]. The activation energies calculated for the addition reactions were found to be reduced by the assistance of a water molecule, which makes easier the proton-transfer process in these alkylation reactions by at least 12.9, 10.5, and 6.0 kcal mol(-1) [at the B3LYP/AUG-cc-pVTZ//B3LYP/6-311+G(d,p) level], for ammonia, water, and hydrogen sulfide, respectively. A proper comparison of an uncatalyzed with a water-catalyzed reaction mechanism has been made on the basis of activation Gibbs free energies. In gas-phase alkylation of ammonia and water by o-QM, reactions assisted by an additional water molecule H-bonded to o-QM (water-catalyzed mechanism) are favored over their uncatalyzed counterparts by 5.6 and 4.0 kcal mol(-1) [at the B3LYP/6-311+G(d,p) level], respectively. In contrast, the hydrogen sulfide alkylation reaction in the gas phase shows a slight preference for a direct alkylation without water assistance, even though the free energy difference (DeltaDeltaG(#)) between the two reaction mechanisms is very small (by 1.0 kcal mol(-1) at the B3LYP/6-311+G(d,p),S(2df) level of theory). The bulk solvent effect, evaluated by the C-PCM model, significantly modifies the relative importance of the uncatalyzed and water-assisted alkylation mechanism by o-QM in comparison to the case in the gas phase. Unexpectedly, the uncatalyzed mechanism becomes highly favored over the catalyzed one in the alkylation reaction of ammonia (by 7.0 kcal mol(-1)) and hydrogen sulfide (by 4.0 kcal mol(-1)). In contrast, activation induced by water complexation still plays an important role in the o-QM hydration reaction in water as solvent.