A theoretical study of the mechanism of the nucleotidyl transfer reaction catalyzed by yeast RNA polymerase II

The mechanism of the nucleotidyl transfer reaction catalyzed by yeast RNA polymerase II has been investigated using molecular mechanics and quantum mechanics methods.Molecular dynamics(MD) simulations were carried out using the TIP3 water model and generalized solvent boundary potential(GSBP) by CHARMM based on the X-ray crystal structure.Two models of the ternary elongation complex were constructed based on CHARMM MD calculations.All the species including reactants,transition states,intermediates,and products were optimized using the DFT-PBE method coupled with the basis set DZVP and the auxiliary basis set GEN-A2.Three pathways were explored using the DFT method.The most favorable reaction pathway involves indirect proton migration from the RNA primer 3’-OH to the oxygen atom of-phosphate via a solvent water molecule,proton rotation from the oxygen atom of-phosphate to the-phosphate side,the RNA primer 3’-O nucleophilic attack on the-phosphorus atom,and P-O bond breakage.The corresponding reaction potential profile was obtained.The rate limiting step,with a barrier height of 21.5 kcal/mol,is the RNA primer 3’-O nucleophilic attack,rather than the commonly considered proton transfer process.A high-resolution crystal structure including crystallographic water molecules is required for further studies.     

Theoretical investigation on the mechanisms of transfer hydrogenation of ketones catalyzed by iridium complexes

The detailed catalytic mechanisms on transfer hydrogenation of ketones are explored by employing the representative reaction of 3-pentanone and 2-propanol catalyzed by the model complex IrH₃[(Me₂PC₂H₄)₂NH], derived from the catalyst IrH₃[(ⁱPr₂PC₂H₄)₂NH], with the aid of the density functional theory calculations. The geometrical transformation from an octahedron to a Y-type involved in the catalytic cycle is also elucidated in terms of molecular theory of transition metal complexes. The trend for the variation of Ir-N bond distance is also analyzed.     

Direct hydride transfer in the reaction mechanism of quinoprotein alcohol dehydrogenases: a quantum mechanical investigation

Oxidation of alcohols by direct hydride transfer to the pyrroloquinoline quinone (PQQ) cofactor of quinoprotein alcohol dehydrogenases has been studied using ab initio quantum mechanical methods. Energies and geometries were calculated at the 6-31G(d,p) level of theory. Comparison of the results obtained for PQQ and several derivatives with available structural and spectroscopic data served to judge the feasibility of the calculations. The role of calcium in the enzymatic reaction mechanism has been investigated. Transition state searches have been conducted at the semiempirical and STO-3G(d) level of theory. It is concluded that hydride transfer from the Calpha-position of the substrate alcohol (or aldehyde) directly to the C(5) carbon of PQQ is energetically feasible. Copyright 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1732-1749, 2001     

Hydride transfer reaction catalyzed by hyperthermophilic dihydrofolate reductase is dominated by quantum mechanical tunneling and is promoted by both inter- and intramonomeric correlated motions

Simulations of hydride and deuteride transfer catalyzed by dihydrofolate reductase from the hyperthermophile Thermotoga maritima (TmDHFR) are presented. TmDHFR was modeled with its active homodimeric quaternary structure, where each monomer has three subdomains. The potential energy function was a combined quantum mechanical and molecular mechanical potential (69 atoms were treated quantum mechanically, and 35 287, by molecular mechanics). The calculations of the rate constants by ensemble-averaged variational transition state theory with multidimensional tunneling predicted that hydride and deuteride transfer at 278 K proceeded with 81 and 80% by tunneling. These percentages decreased to 50 and 49% at 338 K. The kinetic isotope effect was dominated by contributions of bound vibrations and decreased from 3.0 to 2.2 over the temperature range. The calculated rates for hydride and deuteride transfer catalyzed by the hypothetical monomer were smaller by approximately 2 orders of magnitude. At 298 K tunneling contributed 73 and 66% to hydride and deuteride transfer in the monomer. The decreased catalytic efficiency of the monomer was therefore not the result of a decrease of the tunneling contribution but an increase in the quasi-classical activation free energy. The catalytic effect was associated in the dimer with correlated motions between domains as well as within and between subunits. The intrasubunit correlated motions were decreased in the monomer when compared to both native dimeric TmDHFR and monomeric E. coli enzyme. TmDHFR and its E. coli homologue involve similar patterns of correlated interactions that affect the free energy barrier of hydride transfer despite only 27% sequence identity and different quaternary structures.     

Benchmark calculations on models of the phosphoryl transfer reaction catalyzed by protein kinase A

Protein phosphorylation has been proved to be of great importance in many stages of cell life. In the last few years, its reaction mechanism has been extensively studied. In this work we present the analysis of the performances of several computational methods with different computational costs (from multilevel to semiempirical) to point out the best method to be used at each level in the study of phosphoryl transfer. Finally, we center on the semiempirical methods, and mainly on the AM1/d Hamiltonian with different sets of parameters, which will permit hybrid quantum mechanics/molecular mechanics (QM/MM) free energy calculations on big models at an acceptable computational cost. We have used quite a large set of molecules and model reactions to test the computational methods, reproducing all the chemical steps involved in the mainly accepted reaction pathways for the protein phosphorylation. In the end, we also present the results for an enlarged model, cut out from an entire biological model: we compare the 2-D PES at the B3LYP and AM1/d levels with the purpose of obtaining a correction for the semiempirical method. The AM1/d-PhoT semiempirical parameterization corrected using single-point energy calculations at the B3LYP/MG3S level seems to be suitable to carry out reliable QM/MM calculations of the complete biological system.     

A triad interaction in the fingers subdomain of DNA polymerase beta controls polymerase activity

DNA polymerase beta (pol beta) is the main polymerase involved in the base excision repair pathway responsible for repairing damaged bases in the DNA. Previous studies on the H285D mutant of pol beta suggested that the C-terminal region of the polymerase is important for polymerase function. In this study, the C-terminal region of pol beta was mutated to assess its role in polymerization. Kinetic experiments showed that the C-terminal region is required for wild-type polymerase activity. Additionally, an interaction between the fingers and palm subdomain revealed itself to be required for polymerase activity. The E316R mutant of pol beta was shown to have a 29,000-fold reduction in polymerization rate with no reduction in nucleotide binding, suggesting that there exists a noncovalent mechanistic step between nucleotide binding and nucleophilic attack of the primer 3'-hydroxyl group on the α-PO(4) of the nucleotide. Molecular modeling studies of the E316R mutant demonstrate that disrupting the interaction between Arg182 and Glu316 disrupts the packing of side chains in the hydrophobic hinge region and may be hampering the conformational change during polymerization. Taken together, these data demonstrate that the triad interaction of Arg182, Glu316, and Arg333 is crucial for polymerase function.     

DFT investigation on the mechanism of the deacetylation reaction catalyzed by LpxC

LpxC is a key enzyme in the biochemical synthesis of Lipid A, an important outer cell-membrane component found in a number of pathogenic bacteria. Using DFT, we have investigated the binding of the substrate within its active site as well as the deacetylation mechanism it catalyzes. The substrate is found to preferentially coordinate to the active site Zn2+ via its carbonyl oxygen between a Zn2+-bound H2O and an adjacent threonine (Thr191). Furthermore, upon substrate binding a nearby Glu78 residue is found to readily deprotonate the remaining Zn2+-bound H2O. Unlike several related metallopeptidases, the mechanism of LpxC is found to proceed via four steps; (i) initial hydroxylation of the substrates' carbonyl carbon to give a gem-diolate intermediate, (ii) protonation of the amide nitrogen by the histidine His265-H+, (iii) a barrier-less change in the active site-intermediate hydrogen-bond network and finally, (iv) C-N bond cleavage. Notably, the rate-determining step of the mechanism of LpxC is found to be the initial hydroxylation while the final C-N bond cleavage occurs with an overall barrier of 23.6 kJ mol-1. Furthermore, LpxC uses a general acid/base pair mechanism as indicated by the fact that both His265-H+ and Glu78 are accordingly involved.     

A systematic treatment of three-dimensional quantum mechanical reaction coordinates

Summary The rigorous, collinear, canonical point transformation method with hyper-hyperbolic coordinates is extended to the infinite central mass problem in three dimensions. The initial transformation performed is (x_A, y _A, z _A, x _C, y _C, z _C) (, , , r, R, ), where (, , ) are the Euler angles; r and R are the AB and BC interatomic distances, respectively, and is the angle between r and R. A second transformation is then performed to (, , , , , ), where is the reaction coordinate mimicking the reaction path, and is the vibrational coordinate of the diatom. The transformed spaces are all one-to-one mappings from the original spaces, and thus do not have any three-to-one regions. The transformed momenta and Hamiltonians are derived, and are Hermitian in their respective transformed spaces.     

Inverse-electron-demand hetero-Diels-Alder reaction of beta,gamma-unsaturated alpha-ketophosphonates catalyzed by prolinal dithioacetals

Some novel prolinal dithioacetal derivatives were studied as catalysts for the inverse-electron-demand hetero-Diels-Alder reaction of enolizable aldehydes and beta,gamma-unsaturated alpha-ketophosphonates. The corresponding 5,6-dihydro-4H-pyran-2-ylphosphonates were obtained in good ee values (up to 94% ee).     

Finding the possible mechanisms for a given type of overall reaction

A systematic way to derive all thea priori possible mechanisms classified according to the numbers, ϱ, of elementary or molecular reaction steps was presented recently. The method is now applied to overall reactions of the type A +B ⇒ C+D. The “laminar” (elementary steps' stoichiometric coefficients unity) mechanisms that turn out possible are derived and listed for the important ϱ = 2 case. The ϱ = 3 three step mechanisms were reported in another paper. Many examples from actual chemical, biochemical mechanisms, such as S_E1, S_N1, S_E2, S_N2, activated complex theory of abstraction reactions, free radical chain propagation and biochemical electron transport chains are given. By far most mechanisms usually encountered are of the laminar type, although the method also gives the “turbulent” ones, i.e. with some stoichiometric coefficients in elementary steps larger than one or some species occurring in many of the steps. Alexander von Humboldt Senior Scientist Awardee, West Germany. Predoctoral Research Assistant, Yale University 1973–74. Then at the Roswell Park Memorial Institute, Department of Experimental Therapeutics, Buffalo, New York 14203.     

A new mechanism of acyl group transfer in the reaction catalyzed by D-glyceraldehyde-3-phosphate dehydrogenase

D-Glyceraldehyde-3-phosphate dehydrogenase (GAPDH, EC 1.2.1.12) catalyzes the oxidative phosphorylation of its substrate in a two-step reaction. Asa result of the first, oxidativestep, the covalent intermediate where in 3-phosphoglyceroyl moiety is bound to Cys149 of the active center is subjected to nucleophilic attack by inorganic phosphate, but remains resistant to hydrolytic decomposition. This ensures tight coupling of oxidation with phosphorylation in glycolysis. In this article, we present the experimental evidence for the conversion of GAPDH into a form capable of performing the reaction in the absence of inorganic phosphate. The structural basis for this conversion is the oxidation of a cysteine residue (probably Cys 153) into a sulfenic acid derivative under mild conditions to affect the integrity of the essential Cys 149. As a result, an intram olecular transfer of 3-phosphoglyceroyl group from the active center Cys 149 to Cys 153 becomes possible with subsequent hydrolysis of the sulfenyl carboxylate intermediate.     

Processive replication of single DNA molecules in a nanopore catalyzed by phi29 DNA polymerase

Coupling nucleic acid processing enzymes to nanoscale pores allows controlled movement of individual DNA or RNA strands that is reported as an ionic current/time series. Hundreds of individual enzyme complexes can be examined in single-file order at high bandwidth and spatial resolution. The bacteriophage phi29 DNA polymerase (phi29 DNAP) is an attractive candidate for this technology, due to its remarkable processivity and high affinity for DNA substrates. Here we show that phi29 DNAP-DNA complexes are stable when captured in an electric field across the α-hemolysin nanopore. DNA substrates were activated for replication at the nanopore orifice by exploiting the 3'-5' exonuclease activity of wild-type phi29 DNAP to excise a 3'-H terminal residue, yielding a primer strand 3'-OH. In the presence of deoxynucleoside triphosphates, DNA synthesis was initiated, allowing real-time detection of numerous sequential nucleotide additions that was limited only by DNA template length. Translocation of phi29 DNAP along DNA substrates was observed in real time at ?ngstrom-scale precision as the template strand was drawn through the nanopore lumen during replication.     

An investigation into the mechanism of the imidazolidinone catalyzed cascade reaction

The cascade reaction of α,β-unsaturated butyric aldehydes with 2-methyl furan and chlorinated quinone catalyzed by a (2S,5S)-5-benzyl-2-tert-butyl-3-methylimidazolidin-4-one·TFA was investigated by using density functional theory (DFT) calculations at the PCM(EtOAc)/B3LYP/6-311++G(d,p)//B3LYP/6-31G(d) level to (a) confirm the detailed reaction mechanism and key factors controlling the enantioselectivity; and (b) check the models of the iminium ion formation and hydrolysis process that were carried out in another reaction. Two favorable reaction channels, corresponding to the enantioselectivity of the (2R,3S)-product and (2S,3S)-product, have been characterized. The enantioselectivity is controlled by the steps involved in the formation of the C–C bond and the C–Cl bond in the iminium catalysis and the enamine catalysis, respectively. The calculated results explain the reaction mechanism and the enantioselectivity, which are in agreement with experimental observations, and may be helpful for understanding the reaction mechanism of similar cascade reactions.     

Am1 and pm3 transition structure for the hydride transfer. A model of reaction catalyzed by dihydrofolate reductase

The transition state structure for the hydride transfer in dihydrofolate reductase, DHFR, enzyme has been calculated with analytical gradients at semiempirical levels: AM1 and PM3. The geometry, electronic structure and transition vector components are qualitatively semiempirical level independent. Comparing the transition structures for the hydride transfer step in models of liver alcohol dehydrogenase, formate dehydrogenase, lactate dehydrogenase, and glutathione reductase, the geometries of these stationary points are transferable and invariant. The topology of the transition structures in these enzymes resembles the one calculated in this paper.     

A quantum mechanics study on the reaction mechanism of chalcone formation from p-coumaroyl-CoA and malonyl-CoA catalyzed by chalcone synthase

Chalcone synthase catalyzes formation of phenylpropanoid chalcone from one p-coumaroyl-coenzyme A (CoA) and three malonyl-CoA molecules. In order to elucidate structural and energetic features of the reaction mechanism, we performed the quantum mechanics calculations and obtained the following results. In loading step, only a tetrahedral intermediate is located without transition state (TS). Our results indicate that His303 acts as a H₃₁ donor, but not a hydrogen bond donor, to stabilize the intermediate formation. In decarboxylation step, the reaction proceeds via a TS and is sensitive to the environment. In elongation step, a tetrahedral TS is located. All of the results above support the reaction mechanism and further complement the proposal of Noel JP et al.     

Optimizing hierarchical equations of motion for quantum dissipation and quantifying quantum bath effects on quantum transfer mechanisms

We present an optimized hierarchical equations of motion theory for quantum dissipation in multiple Brownian oscillators bath environment, followed by a mechanistic study on a model donor-bridge-acceptor system. We show that the optimal hierarchy construction, via the memory-frequency decomposition for any specified Brownian oscillators bath, is generally achievable through a universal pre-screening search. The algorithm goes by identifying the candidates for the best be just some selected Padé spectrum decomposition based schemes, together with a priori accuracy control criterions on the sole approximation, the white-noise residue ansatz, involved in the hierarchical construction. Beside the universal screening search, we also analytically identify the best for the case of Drude dissipation and that for the Brownian oscillators environment without strongly underdamped bath vibrations. For the mechanistic study, we quantify the quantum nature of bath influence and further address the issue of localization versus delocalization. Proposed are a reduced system entropy measure and a state-resolved constructive versus destructive interference measure. Their performances on quantifying the correlated system-environment coherence are exemplified in conjunction with the optimized hierarchical equations of motion evaluation of the model system dynamics, at some representing bath parameters and temperatures. Analysis also reveals the localization to delocalization transition as temperature decreases.     

水催化HBrO与OH自由基反应的微观机理

采用密度泛函理论,分别在B3LYP/6-311++g(d,p)和B3LYP/aug-cc-PVTZ理论水平下,系统研究了无水和水催化的OH自由基与HBrO反应,即HBrO+OH和HBrO+OH+H_2O 2个反应的微观反应机理,给出了所有可能发生的反应路径,并指出能量最低的反应通道.对于没有水参与的反应,由于OH自由基进攻HBrO方式不同,存在顺式方向和反式方向2种进攻方式的反应路径;当有一分子水参与反应时,考虑HBrO H_2O复合物与OH自由基的反应和HBrO与H_2O OH复合物2种反应情况,共发现4条不同的反应路径.这2种反应的所有路径均是在OH自由基提取氢之前以氢键复合物形式存在,反应过程均为无势垒加合过程,总反应为放热反应.水对目标反应起催化作用,有效地降低了反应的势垒,可以加快OH自由基和HBrO的消耗速度.     

A density functional theory investigation on the formation mechanisms of DNA interstrand crosslinks induced by chloroethylnitrosoureas

Chloroethylnitrosoureas (CENUs) are an important family of alkylating agents used in the clinical treatment of cancer. Their anticancer mechanism primarily involves the formation of DNA interstrand crosslinks (ICLs) induced by the chloroethyldiazonium ion derived from the decomposition of CENUs. In this work, the mechanism for the formation of ICLs was investigated by density functional theory (DFT) with B3LYP, wB97XD, and M062X functinoals using conductor‐like polarizable continuum model solvent model. Three pathways leading to the formation of three types of G–C crosslinks were compared. G(N1)–C(N3) crosslink is predicted to be the dominant crosslinking product other than G(O⁶)–C(N⁴) and G(N²)–C(O²) crosslinks, which is consistent with the previous results obtained from QM/MM computations. The results indicate that the formation of the G(N1)–C(N3) crosslink via pathway A is the most favorable mechanism from both kinetic and thermodynamic standpoints. In this pathway, the chloroethyldiazonium ion alkylates guanine on the O⁶ site followed by intramolecular cyclization to form O⁶,N1‐ethanoguanine ( 4 ). The cytosine then reacts with intermediate 4 on the C^α atom to yield the G(N1)–C(N3) crosslink. This work provides reasonable explanations for the supposed mechanism of CENUs‐induced ICLs formation obtained from experimental investigations. © 2012 Wiley Periodicals, Inc.