Quantum dynamics of hydride transfer catalyzed by bimetallic electrophilic catalysis: synchronous motion of Mg(2+) and H(-) in xylose isomerase

Xylose isomerase exhibits a bridged-bimetallic active-site motif in which the substrate is bound to two metals connected by a glutamate bridge, and X-ray crystallographic studies suggest that metal movement is involved in the hydride transfer rate-controlling catalytic step. Here we report classical/quantal dynamical simulations of this step that provide new insight into the metal motion. The potential energy surface is calculated by treating xylose with semiempirical molecular orbital theory augmented by a simple valence bond potential and the rest of the system by molecular mechanics. The rate constant for the hydride-transfer step was calculated by ensemble-averaged dynamical simulations including both variational transition-state theory for determination of the statistically averaged dynamical bottleneck and optimized multidimensional tunneling calculations. The dynamics calculations include 25 317 atoms, with quantized vibrational free energy in 89 active-site degrees of freedom, and with 32 atoms moving through static secondary zone transition-state configurations in the quantum tunneling simulation. Our simulations show that the average Mg-Mg distance R increases monotonically as a function of the hydride-transfer progress variable z. The range of the average R along the reaction path is consistent with the X-ray structure, thus providing a dynamical demonstration of the postulated role of Mg in catalysis. We also predicted the primary deuterium kinetic isotope effect (KIE) for the chemical step. We calculated a KIE of 3.8 for xylose at 298 K, which is consistent with somewhat smaller experimentally observed KIEs for glucose substrate at higher temperatures. More than half of our KIE is due to tunneling; neglecting quantum effects on the reaction coordinate reduces the calculated KIE to 1.8.     

Ab initio path‐integral calculations of kinetic and equilibrium isotope effects on base‐catalyzed RNA transphosphorylation models

Detailed understandings of the reaction mechanisms of RNA catalysis in various environments can have profound importance for many applications, ranging from the design of new biotechnologies to the unraveling of the evolutionary origin of life. An integral step in the nucleolytic RNA catalysis is self‐cleavage of RNA strands by 2′‐O‐transphosphorylation. Key to elucidating a reaction mechanism is determining the molecular structure and bonding characteristics of transition state. A direct and powerful probe of transition state is measuring isotope effects on biochemical reactions, particularly if we can reproduce isotope effect values from quantum calculations. This article significantly extends the scope of our previous joint experimental and theoretical work in examining isotope effects on enzymatic and nonenzymatic 2′‐O‐transphosphorylation reaction models that mimic reactions catalyzed by RNA enzymes (ribozymes), and protein enzymes such as ribonuclease A (RNase A). Native reactions are studied, as well as reactions with thio substitutions representing chemical modifications often used in experiments to probe mechanism. Here, we report and compare results from eight levels of electronic‐structure calculations for constructing the potential energy surfaces in kinetic and equilibrium isotope effects (KIE and EIE) computations, including a “gold‐standard” coupled‐cluster level of theory [CCSD(T)]. In addition to the widely used Bigeleisen equation for estimating KIE and EIE values, internuclear anharmonicity and quantum tunneling effects were also computed using our recently developed ab initio path‐integral method, that is, automated integration‐free path‐integral method. The results of this work establish an important set of benchmarks that serve to guide calculations of KIE and EIE for RNA catalysis. © 2014 Wiley Periodicals, Inc.     

Elucidating the mechanism for the reduction of nitrite by copper nitrite reductase--a contribution from quantum chemical studies

Density functional methods have been applied to investigate the properties of the active site of copper-containing nitrite reductases and possible reaction mechanisms for the enzyme catalysis. The results for a model of the active site indicate that a hydroxyl intermediate is not formed during the catalytic cycle, but rather a state with a protonated nitrite bound to the reduced copper. Electron affinity calculations indicate that reduction of the T2 copper site does not occur immediately after nitrite binding. Proton affinity calculations are indicative of substantial pK(a) differences between different states of the T2 site. The calculations further suggest that the reaction does not proceed until uptake of a second proton from the bulk solution. They also indicate that Asp-92 may play both a key role as a proton donor to the substrate, and a structural role in promoting catalysis. In the D92N mutant another base, presumably a nearby histidine (His-249) may take the role as the proton donor. On the basis of these model calculations and available experimental evidence, an ordered reaction mechanism for the reduction of nitrite is suggested. An investigation of the binding modes of the nitric oxide product and the nitrite substrate to the model site has also been made, indicating that nitric oxide prefers to bind in an end-on fashion to the reduced T2 site.     

Trapping of an intermediate in the reaction catalyzed by flavin-dependent thymidylate synthase

Thymidylate is a DNA nucleotide that is essential to all organisms and is synthesized by the enzyme thymidylate synthase (TSase). Several human pathogens rely on an alternative flavin-dependent thymidylate synthase (FDTS), which differs from the human TSase both in structure and molecular mechanism. It has recently been shown that FDTS catalysis does not rely on an enzymatic nucleophile and that the proposed reaction intermediates are not covalently bound to the enzyme during catalysis, an important distinction from the human TSase. Here we report the chemical trapping, isolation, and identification of a derivative of such an intermediate in the FDTS-catalyzed reaction. The chemically modified reaction intermediate is consistent with currently proposed FDTS mechanisms that do not involve an enzymatic nucleophile, and it has never been observed during any other TSase reaction. These findings establish the timing of the methylene transfer during FDTS catalysis. The presented methodology provides an important experimental tool for further studies of FDTS, which may assist efforts directed toward the rational design of inhibitors as leads for future antibiotics.     

Kinetic isotope effects on the rate-limiting step of heme oxygenase catalysis indicate concerted proton transfer/heme hydroxylation

Heme oxygenase (HO) catalyzes the O2 and NADPH/cytochrome P450 reductase-dependent conversion of heme to biliverdin, free iron ion, and CO through a process in which the heme participates as both dioxygen-activating prosthetic group and substrate. We earlier confirmed that the first step of HO catalysis is a monooxygenation in which the addition of one electron and two protons to the HO oxy-ferroheme produces ferric-alpha-meso-hydroxyheme (h). Cryoreduction/EPR and ENDOR measurements further showed that hydroperoxo-ferri-HO converts directly to h in a single kinetic step without formation of a Compound I. We here report details of that rate-limiting step. One-electron 77 K cryoreduction of human oxy-HO and annealing at 200 K generates a structurally relaxed hydroperoxo-ferri-HO species, denoted R. We here report the cryoreduction/annealing experiments that directly measure solvent and secondary kinetic isotope effects (KIEs) of the rate-limiting R --> h conversion, using enzyme prepared with meso-deuterated heme and in H2O/D2O buffers to measure the solvent KIE (solv-KIE), and the secondary KIE (sec-KIE) associated with the conversion. This approach is unique in that KIEs measured by monitoring the rate-limiting step are not susceptible to masking by KIEs of other processes, and these results represent the first direct measurement of the KIEs of product formation by a kinetically competent reaction intermediate in any dioxygen-activating heme enzyme.The observation of both solv-KIE(298) = 1.8 and sec-KIE(298) = 0.8 (inverse) indicates that the rate-limiting step for formation of h by HO is a concerted process: proton transfer to the hydroperoxo-ferri-heme through the distal-pocket H-bond network, likely from a carboxyl group acting as a general acid catalyst, occurring in synchrony with bond formation between the terminal hydroperoxo-oxygen atom and the alpha-meso carbon to form a tetrahedral hydroxylated-heme intermediate. Subsequent rearrangement and loss of H2O then generates h.     

Aromatic hydroxylation in a copper bis(imine) complex mediated by a micro-eta2:eta2 peroxo dicopper core: a mechanistic scenario

Detailed mechanistic studies on the ligand hydroxylation reaction mediated by a copper bis(imine) complex are presented. Starting from a structural analysis of the CuI complex and the CuII product with a hydroxylated ligand, the optical absorption and vibrational spectra of starting material and product are analyzed. Kinetic analysis of the ligand hydroxylation reaction shows that O2 binding is the rate-limiting step. The reaction proceeds much faster in methanol than in acetonitrile. Moreover, an inverse kinetic isotope effect (KIE) is evidenced for the reaction in acetonitrile, which is attributed to a sterically congested transition state leading to the peroxo adduct. In methanol, however, no KIE is observed. A DFT analysis of the oxygenation reaction mediated by the micro-eta2:eta2 peroxo core demonstrates that the major barrier after O2 binding corresponds to electrophilic attack on the arene ring. The relevant orbital interaction occurs between the sigma* orbital of the Cu2O2 unit and the HOMO of the ligand. On the basis of the activation energy for the rate-limiting step (18.3 kcal mol(-1)) this reaction is thermally allowed, in agreement with the experimental observation. The calculations also predict the presence of a stable dienone intermediate which, however, escaped experimental detection so far. Reasons for these findings are considered. The implications of the results for the mechanism of tyrosinase are discussed.     

Insight into enzymatic C-F bond formation from QM and QM/MM calculations

The C-F bond-forming step in the fluorinase, the only native fluorination enzyme characterized to date, has been studied. The enzyme catalyzes the reaction between S-adenosyl-L-methionine (SAM) and fluoride ions to form 5'-fluoro-5'-deoxyadenosine (5'-FDA) and L-methionine. To obtain an insight into the mechanism of this unusual enzymatic reaction and to elucidate the role of the enzyme in catalysis, we have explored the conformational energetics of SAM and the intrinsic reactivity patterns of SAM and fluoride with DFT (BP86) and continuum solvent methods, before investigating the full enzymatic system with combined DFT/CHARMM calculations. We find that the enzymatic reaction follows an S(N)2 reaction mechanism, concurring with the intrinsic reactivity preferences in solution. The formation of sulfur ylides is thermodynamically strongly disfavored, and an alternative elimination-addition mechanism involving the concerted anti-Markovnikov addition of HF to an enol ether is energetically viable, but kinetically prohibitive. The S(N)2 activation energy is 92 (112) kJ mol(-)(1) in solution, but only 53 (63) kJ mol(-1) in the enzyme, and the reaction energy in the enzyme is -25 (-34) kJ mol(-1) (values in parentheses are B3LYP single-point energies). The fluorinase thus lowers the barrier for C-F bond formation by 39 (49) kJ mol(-)(1). A decomposition analysis shows that the major role of the enzyme is in the preparation and positioning of the substrates.     

Dinuclear Zn(II) complex catalyzed phosphodiester cleavage proceeds via a concerted mechanism: a density functional theory study

Density functional theory (DFT) calculations were used to study the mechanism for the cleavage reaction of the RNA analogue HpPNP (HpPNP = 2-hydroxypropyl-4-nitrophenyl phosphate) catalyzed by the dinuclear Zn(II) complex of 1,3-bis(1,4,7-triazacyclonon-1-yl)-2-hydroxypropane (Zn(2)(L(2)O)). We present a binding mode in which each terminal phosphoryl oxygen atom binds to one zinc center, respectively, and the nucleophilic 2-hydroxypropyl group coordinates to one of the zinc ions, while the hydroxide from deprotonation of a water molecule coordinates to the other zinc ion. Our calculations found a concerted mechanism for the HpPNP cleavage with a 16.5 kcal/mol reaction barrier. An alternative proposed stepwise mechanism through a pentavalent oxyphosphorane dianion reaction intermediate for the HpPNP cleavage was found to be less feasible with a significantly higher energy barrier. In this stepwise mechanism, the deprotonation of the nucleophilic 2-hydroxypropyl group is accompanied with nucleophilic attack in the rate-determining step. Calculations of the nucleophile (18)O kinetic isotope effect (KIE) and leaving (18)O KIE for the concerted mechanism are in reasonably good agreement with the experimental values. Our results indicate a specific-base catalysis mechanism takes place in which the deprotonation of the nucleophilic 2-hydroxypropyl group occurs in a pre-equilibrium step followed by a nucleophilic attack on the phosphorus center. Detailed comparison of the geometric and electronic structure for the HpPNP cleavage reaction mechanisms in the presence/absence of catalyst revealed that the catalyst significantly altered the determining-step transition state to become far more associative or tight, that is, bond formation to the nucleophile was remarkably more advanced than leaving group bond fission in the catalyzed mechanism. Our results are consistent with and provide a reliable interpretation for the experimental observations that suggest the reaction occurs by a concerted mechanism (see Humphry, T.; Iyer, S.; Iranzo, O.; Morrow, J. R.; Richard, J. P.; Paneth, P.; Hengge, A. C. J. Am. Chem. Soc. 2008, 130, 17858-17866) and has a specific-base catalysis character (see Yang, M.-Y.; Iranzo, O.; Richard, J. P.; Morrow, J. R. J. Am. Chem. Soc. 2005, 127, 1064-1065).     

Theoretical studies on the deacylation step of serine protease catalysis in the gas phase, in solution, and in elastase

The deacylation step of serine protease catalysis is studied using DFT and ab initio QM/MM calculations combined with MD/umbrella sampling calculations. Free energies of the entire reaction are calculated in the gas phase, in a continuum solvent, and in the enzyme elastase. The calculations show that a concerted mechanism in the gas phase is replaced by a stepwise mechanism when solvent effects or an acetate ion are added to the reference system, with the tetrahedral intermediate being a shallow minimum on the free energy surface. In the enzyme, the tetrahedral intermediate is a relatively stable species ( approximately 7 kcal/mol lower in energy than the transition state), mainly due to the electrostatic effects of the oxyanion hole and Asp102. It is formed in the first step of the reaction, as a result of a proton transfer from the nucleophilic water to His57 and of an attack of the remaining hydroxyl on the ester carbonyl. This is the rate-determining step of the reaction, which requires approximately 22 kcal/mol for activation, approximately 5 kcal/mol less than the reference reaction in water. In the second stage of the reaction, only small energy barriers are detected to facilitate the proton transfer from His57 to Ser195 and the breakdown of the tetrahedral intermediate. Those are attributed mainly to a movement of Ser195 and to a rotation of the His57 side chain. During the rotation, the imidazolium ion is stabilized by a strong H-bond with Asp102, and the C(epsilon)(1)-H...O H-bond with Ser214 is replaced by one with Thr213, suggesting that a "ring-flip mechanism" is not necessary as a driving force for the reaction. The movements of His57 and Ser195 are highly correlated with rearrangements of the binding site, suggesting that product release may be implicated in the deacylation process.     

Organocatalytic Enantioselective Pictet–Spengler Reaction of α-Ketoesters: Development and Application to the Total Synthesis of (+)-Alstratine A

We report herein an asymmetric Pictet–Spengler reaction of α-ketoesters. In the presence of a catalytic amount of simple alanine-derived squaramide and p-nitrobenzoic acid, reaction of tryptamines with methyl 2-oxoalkanoates afforded the corresponding 1-alkyl-1-methoxycarbonyl tetrahydro-β-carbolines (THBCs) in high yields and ee values. A primary kinetic isotope effect (KIE=4.5) using C2-deteurium-labelled tryptamine indicates that rearomatization through deprotonation of the pentahydro-β-carbolinium ion could be the rate- and enantioselectivity-determining step. A concise enantioselective total synthesis of (+)-alstratine A, a hexacyclic cagelike monoterpene indole alkaloid, featuring this reaction as a key step, was subsequently accomplished. Remeasurement of the [a]_D value of the natural product indicates that natural alstratine A is dextrorotatory rather than levorotatory as it was initially reported in the isolation paper.     

Proposed mechanism for the reaction catalyzed by a diterpene cyclase,aphidicolan-16beta-ol synthase: experimental results on biomimetic cyclization and examination of the cyclization pathway by ab initio calculations

To examine the mechanism of the cyclization reaction catalyzed by aphidicolan-16beta-ol synthase (ACS), which is a key enzyme in the biosynthesis of diterpene aphidicolin, a specific inhibitor of DNA polymerase alpha, skeletal rearrangement of 2a and biomimetic cyclization of 4b were employed. The structures of the reaction products, which reflect penultimate cation intermediates, allowed us to propose a detailed reaction pathway for the Lewis acid-catalyzed cyclizations and rearrangements. Isolation of these products in an aphidicolin-producing fungus led us to speculate that the mechanism of the ACS-catalyzed cyclization reaction is the same as that of a nonenzymatic reaction. Ab initio calculations of the acid-catalyzed reaction intermediates and the transition states indicate that the overall reaction catalyzed by ACS is an exothermic process though the reaction proceeds via an energetically disfavored secondary cation-like transition state. In conjunction with the solvent effect in the acid-catalyzed reactions, this indicates that the actual role of ACS is to provide a template which enforces conformations of the intermediate cations leading to the productive cyclization although it has been believed that the cation-pi interaction between cation intermediates and aromatic amino acid residues in the active site is important for the enzymatic catalysis. This study provided important information on the role of various cationic species, especially secondary cation-like structures, in both nonenzymatic and enzymatic reactions.     

手性噁唑磷烷-硼烷催化还原苯乙酮对映选择性的理论研究

应用密度泛函理论BHandH/6-31G**计算方法研究新型手性非金属催化剂1,3,2-噁唑磷烷-硼烷催化还原苯乙酮的对映选择性反应机理, 确定了在反应途径上的反应物、络合物、过渡态, 中间体和对映体中间产物. 计算结果表明, 该对映选择性还原反应是两个平行的分步反应, 对映体产物主要是(R)构型. 采用同样的方法研究了对映选择性还原反应在甲苯中的溶剂化效应. SCRF计算显示, 甲苯溶剂不改变反应的机理, 但能降低各驻点的能量和(R)反应的位垒, 有利于催化还原反应朝生成(R)构型中间产物的方向进行.     

Multiple isotope effect study of the acid-catalyzed hydrolysis of formamide

Multiple isotope effects were measured at the reactive center of formamide during acid-catalyzed hydrolysis in water at 25 degrees C. The mechanism involves a rapid pre-equilibrium protonation of the carbonyl oxygen, followed by the formation of at least one tetrahedral intermediate, which does not appreciably exchange its carbonyl oxygen with the solvent (kh/kex = 55). The pKa for formamide was determined by 15N NMR and found to be about -2.0. The formyl-hydrogen kinetic isotope effect (KIE) is indicative of a transition state that is highly tetrahedral (Dkobs = 0.79); the carbonyl-carbon KIE (13kobs = 1.031) is in agreement with this conclusion. The small leaving-nitrogen KIE (15kobs = 1.0050) is consistent with some step prior to breaking the C-N bond as rate-determining. The carbonyl-oxygen KIE (18kobs = 0.996) points to attack of water as the rate-determining step. On the basis of these results, a mechanism is proposed in which attachment of the nucleophile to a protonated formamide molecule is rate determining.     

Intramolecular kinetic isotope effect in hydride transfer from dihydroacridine to a quinolinium ion. Rejection of a proposed two-step mechanism with a kinetically significant intermediate

The intramolecular kinetic isotope effect (KIE) for hydride transfer from 10-methyl-9,10-dihydroacridine to 1-benzyl-3-cyanoquinolinium ion has been found to be 5-6 by both (1)H NMR and mass spectrometry. This KIE is consistent with other hydride transfers. It is inconsistent with the high intermolecular KIEs derived by fitting to a two-step mechanism with a kinetically significant intermediate complex, and it is inconsistent with the strong temperature dependence of those KIEs. We therefore reject the two-step mechanism for this reaction, and we suggest that other cases proposed to follow this mechanism are in error.     

胶束催化作用下实现聚苯乙烯的氯甲基化

用紫外分光光度法测定了离子型表面活性剂胶束溶液对聚苯乙烯四氯化碳溶液的增溶性能; 用胶束催化法实现了聚苯乙烯的氯甲基化, 用红外光谱法和佛尔哈德法表征了氯甲基化聚苯乙烯的化学结构与组成; 通过比较阴、 阳离子表面活性剂及结构不同的阳离子表面活性剂的催化效果, 探索了胶束催化的作用机理, 考察了表面活性剂结构对催化作用的影响规律. 结果表明, 表面活性剂胶束溶液可增溶聚苯乙烯的四氯化碳溶液, 随着四氯化碳在胶束中的增溶, 聚苯乙烯可转移至表面活性剂的胶束中; 胶束催化是实现聚苯乙烯的氯甲基化的有效途经, 仅用3.35 g/L的十六甲基三甲基溴化铵(CTAB),  于65 ℃， 5 h内即可使聚苯乙烯大分子链中的苯环氯甲基化程度达到37%; 聚苯乙烯与甲醛、 氯化氢的反应过程由亲电取代和亲核取代串联而成, 阳离子表面活性剂比阴离子表面活性剂的催化作用更加有效, 说明亲核取代是慢步骤; 阳离子表面活性剂疏水链越长, 对聚苯乙烯的增溶效果越好, 催化作用越强.     

Computational study of phosphatase activity in soluble epoxide hydrolase: high efficiency through a water bridge mediated proton shuttle

Recently, a new branch of fatty acid metabolism has been opened by the novel phosphatase activity found in the N-terminal domain of the, hence bifunctional, soluble epoxide hydrolase (sEH). Importantly, this finding has also provided a new site for drug targeting in sEH's activity regulation. Classical MD and hybrid Car-Parrinello QM/MM calculations have been performed to investigate the reaction mechanism of the phosphoenzyme intermediate formation in the first step of the catalysis. The results support a concerted multi-event reaction mechanism: (1) a dissociative in-line nucleophilic substitution for the phosphoryl transfer reaction; (2) a double proton transfer involved in the formation of a good leaving group in the transition state. The presence of a water bridge in the substrate/enzyme complex allowed an efficient proton shuttle, showing its key role in speeding up the catalysis. The calculated free energy of the favored catalytic pathway is approximately 19 kcal/mol, in excellent agreement with experimental data.     

Quantum mechanical/molecular mechanical study on the mechanism of the enzymatic Baeyer-Villiger reaction

We report a combined quantum mechanical/molecular mechanical (QM/MM) study on the mechanism of the enzymatic Baeyer-Villiger reaction catalyzed by cyclohexanone monooxygenase (CHMO). In QM/MM geometry optimizations and reaction path calculations, density functional theory (B3LYP/TZVP) is used to describe the QM region consisting of the substrate (cyclohexanone), the isoalloxazine ring of C4a-peroxyflavin, the side chain of Arg-329, and the nicotinamide ring and the adjacent ribose of NADP(+), while the remainder of the enzyme is represented by the CHARMM force field. QM/MM molecular dynamics simulations and free energy calculations at the semiempirical OM3/CHARMM level employ the same QM/MM partitioning. According to the QM/MM calculations, the enzyme-reactant complex contains an anionic deprotonated C4a-peroxyflavin that is stabilized by strong hydrogen bonds with the Arg-329 residue and the NADP(+) cofactor. The CHMO-catalyzed reaction proceeds via a Criegee intermediate having pronounced anionic character. The initial addition reaction has to overcome an energy barrier of about 9 kcal/mol. The formed Criegee intermediate occupies a shallow minimum on the QM/MM potential energy surface and can undergo fragmentation to the lactone product by surmounting a second energy barrier of about 7 kcal/mol. The transition state for the latter migration step is the highest point on the QM/MM energy profile. Gas-phase reoptimizations of the QM region lead to higher barriers and confirm the crucial role of the Arg-329 residue and the NADP(+) cofactor for the catalytic efficiency of CHMO. QM/MM calculations for the CHMO-catalyzed oxidation of 4-methylcyclohexanone reproduce and rationalize the experimentally observed (S)-enantioselectivity for this substrate, which is governed by the conformational preferences of the corresponding Criegee intermediate and the subsequent transition state for the migration step.     

Selective oxidation of styrene on an oxygen-adsorbed Au(111): A density functional theory study

The reaction mechanism for the styrene selective oxidation on the oxygen preadsorbed Au(111) surface has been studied by the density functional theory calculation with the periodic slab model. The calculated results showed that the process of reaction includes two steps: forming the oxametallacycle intermediate (OMME) and then producing the products. It was found that the second step, from OMME to product is the rate-controlling step, which is similar to ethylene selective oxidation on Ag. Importantly, the present density-functional-theory calculation results suggested that the mechanism via the OMME (2) (i.e. the preadsorbed atomic oxygen bound to the CH2 group involved in C6H5--CH=CH2) to produce styrene epoxide is kinetically favored than that of OMME (1).     

Formation of yohimbanones via a novel rearrangement

An investigation of the reaction of tryptophan and its methyl ester with ninhydrin has been conducted. In the reaction of tryptophan with ninhydrin only one product, yohimbanone (1), is isolated. In contrast, an intermediate, Pictet-Spengler product 3, is isolated from the reaction of tryptophan methyl ester with ninhydrin. The isolation of this intermediate provides support for a proposed mechanism of this novel rearrangement, which is presented herein.