On the Reaction of Glycerol Dehydratase with But‐3‐ene‐1,2‐diol

Intriguing inactivation : Calculations suggest that the ability of relatively high‐energy radical intermediates to inactivate glycerol dehydratase (GDH) may reflect a general and hitherto unidentified inactivation mechanism in the reaction of coenzyme B₁₂‐dependent enzymes and 3‐unsaturated 1,2‐diols (see scheme; AdoCbl: adenosylcobalamin or coenzyme B₁₂).

     

Understanding the Role of Water in Aqueous Ruthenium‐Catalyzed Transfer Hydrogenation of Ketones.

We report an accurate computational study of the role of water in transfer hydrogenation of formaldehyde with a ruthenium‐based catalyst using a water‐specific model. Our results suggest that the reaction mechanism in aqueous solution is significantly different from that in the gas phase or in methanol solution. Previous theoretical studies have shown a concerted hydride and proton transfer in the gas phase (M. Yamakawa, H. Ito, R. Noyori, J. Am. Chem. Soc. 2000 , 122, 1466–1478;J.‐W. Handgraaf, J. N. H. Reek, E. J. Meijer, Organometallics 2003 , 22, 3150–3157; D. A. Alonso, P. Brandt, S. J. M. Nordin, P. G. Andersson, J. Am. Chem. Soc. 1999 , 121, 9580–9588; D. G. I. Petra, J. N. H. Reek, J.‐W. Handgraaf, E. J. Meijer, P. Dierkes, P. C. J. Kamer, J. Brussee, H. E. Schoemaker, P. W. N. M. van Leeuwen, Chem. Eur. J. 2000 , 6, 2818–2829), whereas a delayed, solvent‐mediated proton transfer has been observed in methanol solution (J.‐W. Handgraaf, E. J. Meijer, J. Am. Chem. Soc. 2007 , 129, 3099–3103). In aqueous solution, a concerted transition state is observed, as in the previous studies. However, only the hydride is transferred at that point, whereas the proton is transferred later by a water molecule instead of the catalyst.     

Structural Characteristics of Graphane‐Type C and BN Nanostructures by Periodic Local MP2 Approach

The structural characteristics of fully‐hydrogenated carbon and boron nitride mono‐ and multilayer slabs, together with nanotubes derived from the slabs, are investigated mainly by means of periodic local second‐order Møller–Plesset perturbation (LMP2) calculations and the results are compared with Hartree–Fock (HF), density functional theory (DFT), and dispersion function‐augmented DFT (DFT‐D) obtained ones. The investigated systems are structurally analogous to (111) and (110) slabs of diamond, where the hydrogenated (111) slab of diamond corresponds to the experimentally known graphane. Multilayering of monolayers and nanotubes is energetically favorable at the LMP2 level for both C and BN, while HF and DFT are not able to reproduce this behavior for CH systems. The work highlights the importance of utilizing methods capable of properly describing weak interactions in the investigation of dispersively‐bound systems such as the multilayered graphanes and the corresponding nanotubes.     

Theoretical Study on the Complexes of Benzene with Isoelectronic Nitrogen‐Containing Heterocycles

The π–π interactions between benzene and the aromatic nitrogen heterocycles pyridine, pyrimidine, 1,3,5‐triazine, 1,2,3‐triazine, 1,2,4,5‐tetrazine, and 1,2,3,4,5‐pentazine are systematically investigated. The T‐shaped structures of all complexes studied exhibit a contraction of the C? H bond accompanied by a rather large blue shift (40–52 cm^?1) of its stretching frequency, and they are almost isoenergetic with the corresponding displaced‐parallel structures at reliable levels of theory. With increasing number of nitrogen atoms in the heterocycle, the geometries, frequencies, energies, percentage of s character at C, and the electron density in the C? H σ antibonding orbital of the complexes all increase or decrease systematically. Decomposition analysis of the total binding energy showed that for all the complexes, the dispersion energy is the dominant attractive contribution, and a rather large attraction originating from electrostatic contribution is compensated by its exchange counterpart.     

Static and Frequency‐Dependent Dipole–Dipole Polarizabilities of All Closed‐Shell Atoms up to Radium: A Four‐Component Relativistic DFT Study

We test the performance of four‐component relativistic density functional theory by calculating the static and frequency‐dependent electric dipole–dipole polarizabilities of all (ground‐state) closed‐shell atoms up to Ra. We consider 12 nonrelativistic functionals, including three asymptotically shape‐corrected functionals, by using two smooth interpolation schemes introduced by the Baerends group: the gradient‐regulated asymptotic connection (GRAC) procedure and the statistical averaging of (model) orbital potentials (SAOP). Basis sets of doubly augmented triple‐zeta quality are used. The results are compared to experimental data or to accurate ab initio results. The reference static electric dipole polarizability of palladium has been obtained by finite‐field calculations using the coupled‐cluster singles, doubles, and perturbative triples method within this work. The best overall performance is obtained using hybrid functionals and their GRAC shape‐corrected versions. The performance of SAOP is among the best for nonhybrid functionals for Group 18 atoms but its precision degrades when considering the full set of atoms. In general, we find that conclusions based on results obtained for the rare‐gas atoms are not necessarily representative of the complete set of atoms. GRAC cannot be used with effective core potentials since the asymptotic correction is switched on in the core region.     

Model studies of the optical rotation,and theoretical determination of its sign for β‐pinene and trans‐pinane

We have carried out extensive studies on the basis set dependence of the calculated specific optical rotation (OR) in molecules at the level of the time–dependent Hartree–Fock and density functional approximations. To reach the limits of the basis set saturation, we have devised an artificial model, the asymmetrically deformed (chiral) methane (CM) molecule. This small system permits to use basis sets which are prohibitively large for real chiral molecules and yet shows all the important features of the basis set dependence of the OR values. The convergence of the OR has been studied with n‐aug‐cc‐pVXZ basis sets of Dunning up to the 6–ζ. In a parallel series of calculations, we have used the recently developed large polarized (LPolX) basis sets. The relatively small LPolX sets have been shown to be competitive to very large n‐aug‐cc‐pVXZ basis sets. The conclusions reached in calculations of OR in CM concerning the usefulness of LPolX basis sets have been further tested on (S)‐methyloxirane and (S)‐fluoro‐oxirane. The smallest set of the LPolX family (LPol–ds) has been found to yield OR values of similar quality as those obtained with much larger Dunning's aug‐cc‐pVQZ basis set. These results have encouraged us to carry out the OR calculations with LPol–ds basis sets for systems as large as β‐pinene and trans‐pinane. In both cases, our calculations have lead to the correct sign of the OR value in these molecules. This makes the relatively small LPol–ds basis sets likely to be useful in OR calculations for large molecules. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Nucleophile‐Dependent Regio‐ and Stereoselective Ring Opening of 1‐Azoniabicyclo[3.1.0]hexane Tosylate

1‐[(1R)‐(1‐Phenylethyl)]‐1‐azoniabicyclo[3.1.0]hexane tosylate was generated as a stable bicyclic aziridinium salt from the corresponding 2‐(3‐hydroxypropyl)aziridine upon reaction with p‐toluenesulfonyl anhydride. This bicyclic aziridinium ion was then treated with various nucleophiles including halides, azide, acetate, and cyanide in CH₃CN to afford either piperidines or pyrrolidines through regio‐ and stereoselective ring opening, mediated by the characteristics of the applied nucleophile. On the basis of DFT calculations, ring‐opening reactions under thermodynamic control yield piperidines, whereas reactions under kinetic control can yield both piperidines and pyrrolidines depending on the activation energies for both pathways.     

Enthalpy and Entropy Barriers Explain the Effects of Topology on the Kinetics of Zeolite‐Catalyzed Reactions

The methylation of ethene, propene, and trans‐2‐butene on zeolites H‐ZSM‐58 (DDR), H‐ZSM‐22 (TON), and H‐ZSM‐5 (MFI) is studied to elucidate the particular influence of topology on the kinetics of zeolite‐catalyzed reactions. H‐ZSM‐58 and H‐ZSM‐22 are found to display overall lower methylation rates compared to H‐ZSM‐5 and also different trends in methylation rates with increasing alkene size. These variations may be rationalized based on a decomposition of the free‐energy barriers into enthalpic and entropic contributions, which reveals that the lower methylation rates on H‐ZSM‐58 and H‐ZSM‐22 have virtually opposite reasons. On H‐ZSM‐58, the lower methylation rates are caused by higher enthalpy barriers, owing to inefficient stabilization of the reaction intermediates in the large cage‐like pores. On the other hand, on H‐ZSM‐22, the methylation rates mostly suffer from higher entropy barriers, because excessive entropy losses are incurred inside the narrow‐channel structure. These results show that the kinetics of crucial elementary steps hinge on the balance between proper stabilization of the reaction intermediates inside the zeolite pores and the resulting entropy losses. These fundamental insights into their inner workings are indispensable for ultimately selecting or designing better zeolite catalysts.     

Theoretical Study of the Stability of the DNA Duplexes Modified by a Series of Hydrophobic Base Analogues

The geometries of a 13 mer of a DNA double helix (5′‐GCGTAC A CATGCG‐3′) were determined by molecular dynamics simulations using a Cornell et al. empirical force field. The bases in the central base pair (shown in bold) were replaced (one or both) by a series of hydrophobic base analogues (phenyl, biphenyl, phenylnaphathalene, phenylanthracene and phenylphenanthrene). Due to the large fluctuations of the systems, an average geometry could not be determined. The interaction energies of the Model A, which consisted of three central steps of a duplex without a sugar phosphate backbone, taken from molecular dynamics simulations (geometry sampled every 1 ps), were calculated by the self‐consistent charge density functional based tight‐binding (SCC‐DFTB‐D) method and were subsequently averaged. The higher the stability of the systems the higher the aromaticity of the base analogues. To estimate the desolvation energy of the duplex, the COSMO continuum solvent model was used and the calculations were provided on a larger model, Model B (the three central steps of the duplex with a sugar phosphate backbone neutralised by H atoms), taken from molecular dynamics simulations (geometry sampled every 200 ps) and subsequently averaged. The selectivity of the base analogue pairs was ascertained (Model B) by including the desolvation energy and the interaction energy of both strands, as determined by the SCC‐DFTB‐D method. The highest selectivity was found for a phenylphenanthrene. Replacing the nucleic acid bases with a base analogue leads to structural changes of the central pair. Only with the smallest base analogues (phenyl) does the central base pair stay planar. When passing to larger base analogues the central base pair is usually stacked.     

Quantum chemical modeling of 1,1‐proton transfer reaction catalyzed by a cofactor‐independent α‐methylacyl‐CoA racemase

α‐Methylacyl‐CoA racemases (AMACR) are essential enzymes for branched‐chain lipids and drugs metabolism. AMACR catalyzes the chiral inversion of (2R) and (2S)‐methylacyl‐CoA esters in both directions. In this study, we investigated the catalytic mechanism of Mycobacterium tuberculosis (MCR) α‐methylacyl‐CoA racemase by using the density functional theory with the hybrid functional B3LYP. Our calculations elucidate and support the mechanism proposed by Prasenjit Bhaumik. His126 and Asp156 serve as the acid/base‐pair residues in the 1,1‐proton transfer catalytic reaction. From the optimized structures, it can be seen that an enolate intermediate is formed and the possibility of forming a ketene or a carbanion intermediate is excluded. By comparing the energy barriers, we could consider that the deprotonation step is the rate‐determined step in the invert direction from (S)‐ to (R)‐enantiomer. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Computational investigation of hydrogen abstraction from 2-aminoethanol by the 1,5-dideoxyribose-5-yl radical: a model study of a reaction occurring in the active site of ethanolamine ammonia lyase

Hydrogen abstraction from 2-aminoethanol by the 5'-deoxyadenosyl radical, which is formed upon Co--C bond homolysis in coenzyme B(12), was investigated by theoretical means with employment of the DFT (B3LYP) and ab initio (MP2) approaches. As a model system for the 5'-deoxyadenosyl moiety the computationally less demanding 1,5-dideoxyribose was employed; two conformers, which differ in ring conformation (C2- and C3-endo), were considered. If hydrogen is abstracted from "free" substrate by the C2-endo conformer of the 1,5-dideoxyribose-5-yl radical, the activation enthalpy is 16.7 kcal mol(-1); with the C3-endo counterpart, the value is 17.3 kcal mol(-1). These energetic requirements are slightly above the activation enthalpy limit (15 kcal mol(-1)) determined experimentally for the rate-determining step of the sequence, that is, hydrogen delivery from 5'-deoxyadenosine to the product radical. The activation enthalpy is lower when the substrate interacts with at least one amino acid from the active site. According to the computations, when a His model system partially protonates the substrate the activation enthalpy is 4.5 kcal mol(-1) for the C3-endo conformer and 5.8 kcal mol(-1) for the C2-endo counterpart. As hydrogen abstraction from the fully as well as the partially protonated substrate is preceded by the formation of quite stable encounter complexes, the actual activation barriers are around 13-15 kcal mol(-1). A synergistic interaction of 2-aminoethanol with two amino acids where His partially protonates the NH(2) group and Asp partially deprotonates the OH group of the substrate results in an activation enthalpy of 12.4 kcal mol(-1) for the C3-endo conformer and 13.2 kcal mol(-1) for the C2-endo counterpart. However, if encounter complexes exist in the active site, the actual activation barriers are much higher (>25 kcal mol(-1)) than that reported for the rate-determining step. These findings together with previous computations suggest that the energetics of the initial hydrogen abstraction decrease with an interaction of the substrate with only a protonating auxiliary, but for the rearrangement of the radical the synergistic effects of two auxiliaries are essential to pull the barrier below the limit of 15 kcal mol(-1).     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.     

Accurate pKa Calculation of the Conjugate Acids of Alkanolamines,Alkaloids and Nucleotide Bases by Quantum Chemical Methods

The pK_a of the conjugate acids of alkanolamines, neurotransmitters, alkaloid drugs and nucleotide bases are calculated with density functional methods (B3LYP, M08‐HX and M11‐L) and ab initio methods (SCS‐MP2, G3). Implicit solvent effects are included with a conductor‐like polarizable continuum model (CPCM) and universal solvation models (SMD, SM8). G3, SCS‐MP2 and M11‐L methods coupled with SMD and SM8 solvation models perform well for alkanolamines with mean unsigned errors below 0.20 pK_a units, in all cases. Extending this method to the pK_a calculation of 35 nitrogen‐containing compounds spanning 12 pK_a units showed an excellent correlation between experimental and computational pK_a values of these 35 amines with the computationally low‐cost SM8/M11‐L density functional approach.     

Quantum chemical study on the abstraction reaction of alkylidenegermylene with oxirane and thiirane

The mechanisms of theion reaction of alkylidenegermylene with oxirane and thiirane have been characterized detail in using density functional theory, as well as ab initio method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies were calculated by CCSD(T)/6-311G(d)//B3LYP/6-311G(d,p) method for the involved conformations. The results show that the reaction pathways for both reactions consist of two ways: (1) the reactants can yield bent products (P1; P4) by syn-isomers; (2) the reactants can also yield three-membered products (P2; P5) by anti-isomers, which then further react with oxirane and thiirane to form the ultimate products (P3-1, P3-2; P6). Furthermore, a comparison with alkylidenecarbene, oxirane, and thiirane was done.     

First detection of a selenenyl fluoride ArSe-F by NMR spectroscopy: the nature of Ar2Se2/XeF2 and ArSe-SiMe3/XeF2 reagents

Arylselenenyl fluorides ArSeF are obtained from diselenides Ar2Se2 or arylselenotrimethylsilanes ArSe-SiMe3, and XeF2. They are detected by low-temperature 19F and 77Se NMR spectroscopy. Substitution in the ortho position of the aromatic ring to provide electronic or steric protection is a requirement for their formation. ArSe--F compounds decompose according to 3 ArSe-F-->[ArSe-SeF2Ar]+ArSe-F-->ArSeF3+Ar2Se2. Reaction energies for this disproportionation as well as that of the sulfur and tellurium homologues have been calculated with MP2, CCSD(T,) and B3 LYP methods. They were found to be increasingly exothermic in the sequence S相似文献