Role of the catalytic triad and oxyanion hole in acetylcholinesterase catalysis: an ab initio QM/MM study

The initial step of the acylation reaction catalyzed by acetylcholinesterase (AChE) has been studied by a combined ab initio quantum mechanical/molecular mechanical (QM/MM) approach. The reaction proceeds through the nucleophilic addition of the Ser203 O to the carbonyl C of acetylcholine, and the reaction is facilitated by simultaneous proton transfer from Ser203 to His447. The calculated potential energy barrier at the MP2(6-31+G) QM/MM level is 10.5 kcal/mol, consistent with the experimental reaction rate. The third residue of the catalytic triad, Glu334, is found to be essential in stabilizing the transition state through electrostatic interactions. The oxyanion hole, formed by peptidic NH groups from Gly121, Gly122, and Ala204, is also found to play an important role in catalysis. Our calculations indicate that, in the AChE-ACh Michaelis complex, only two hydrogen bonds are formed between the carbonyl oxygen of ACh and the peptidic NH groups of Gly121 and Gly122. As the reaction proceeds, the distance between the carbonyl oxygen of ACh and NH group of Ala204 becomes smaller, and the third hydrogen bond is formed both in the transition state and in the tetrahedral intermediate.     

Mechanism of proteolysis in matrix metalloproteinase‐2 revealed by QM/MM modeling

The mechanism of enzymatic peptide hydrolysis in matrix metalloproteinase‐2 (MMP‐2) was studied at atomic resolution through quantum mechanics/molecular mechanics (QM/MM) simulations. An all‐atom three‐dimensional molecular model was constructed on the basis of a crystal structure from the Protein Data Bank (ID: 1QIB), and the oligopeptide Ace‐Gln‐Gly～Ile‐Ala‐Gly‐Nme was considered as the substrate. Two QM/MM software packages and several computational protocols were employed to calculate QM/MM energy profiles for a four‐step mechanism involving an initial nucleophilic attack followed by hydrogen bond rearrangement, proton transfer, and C? N bond cleavage. These QM/MM calculations consistently yield rather low overall barriers for the chemical steps, in the range of 5–10 kcal/mol, for diverse QM treatments (PBE0, B3LYP, and BB1K density functionals as well as local coupled cluster treatments) and two MM force fields (CHARMM and AMBER). It, thus, seems likely that product release is the rate‐limiting step in MMP‐2 catalysis. This is supported by an exploration of various release channels through QM/MM reaction path calculations and steered molecular dynamics simulations. © 2015 Wiley Periodicals, Inc.     

Theoretical perspectives on the reaction mechanism of serine proteases: the reaction free energy profiles of the acylation process

The reaction mechanism of serine proteases (trypsin), which catalyze peptide hydrolysis, is studied theoretically by ab initio QM/MM electronic structure calculations combined with Molecular Dynamics-Free Energy Perturbation calculations. We have calculated the entire reaction free energy profiles of the first reaction step of this enzyme (acylation process). The present calculations show that the rate-determining step of the acylation is the formation of the tetrahedral intermediate, and the breakdown of this intermediate has a small energy barrier. The calculated activation free energy for the acylation is approximately 17.8 kcal/mol at QM/MM MP2/(aug)-cc-pVDZ//HF/6-31(+)G/AMBER level, and this reaction is an exothermic process. MD simulations of the enzyme-substrate (ES) complex and the free enzyme in aqueous phase show that the substrate binding induces slight conformational changes around the active site, which favor the alignment of the reactive fragments (His57, Asp102, and Ser195) together in a reactive orientation. It is also shown that the proton transfer from Ser195 to His57 and the nucleophilic attack of Ser195 to the carbonyl carbon of the scissile bond of the substrate occur in a concerted manner. In this reaction, protein environment plays a crucial role to lowering the activation free energy by stabilizing the tetrahedral intermediate compared to the ES complex. The polarization energy calculations show that the enzyme active site is in a very polar environment because of the polar main chain contributions of protein. Also, the ground-state destabilization effect (steric strain) is not a major catalytic factor. The most important catalytic factor of stabilizing the tetrahedral intermediate is the electrostatic interaction between the active site and particular regions of protein: the main chain NH groups in Gly193 and Ser195 (so-called oxyanion hole region) stabilize negative charge generated on the carbonyl oxygen of the scissile bond, and the main chain carbonyl groups in Ile212 approximately Ser214 stabilize a positive charge generated on the imidazole ring of His57.     

QM/MM study of mechanisms for compound I formation in the catalytic cycle of cytochrome P450cam

In the catalytic cycle of cytochrome P450cam, after molecular oxygen binds as a ligand to the heme iron atom to yield a ferrous dioxygen complex, there are fast proton transfers that lead to the formation of the active species, Compound I (Cpd I), which are not well understood because they occur so rapidly. In the present work, the conversion of the ferric hydroperoxo complex (Cpd 0) to Cpd I has been investigated by combined quantum-mechanical/molecular-mechanical (QM/MM) calculations. The residues Asp(251) and Glu(366) are considered as proton sources. In mechanism I, a proton is transported to the distal oxygen atom of the hydroperoxo group via a hydrogen bonding network to form protonated Cpd 0 (prot-Cpd0: FeOOH(2)), followed by heterolytic O-O bond cleavage that generates Cpd I and water. Although a local minimum is found for prot-Cpd0 in the Glu(366) channel, it is very high in energy (more than 20 kcal/mol above Cpd 0) and the barriers for its decay are only 3-4 kcal/mol (both toward Cpd 0 and Cpd I). In mechanism II, an initial O-O bond cleavage followed by a concomitant proton and electron transfer yields Cpd I and water. The rate-limiting step in mechanism II is O-O cleavage with a barrier of about 13-14 kcal/mol. According to the QM/MM calculations, the favored low-energy pathway to Cpd I is provided by mechanism II in the Asp(251) channel. Cpd 0 and Cpd I are of similar energies, with a slight preference for Cpd I.     

Mechanisms of antibiotic resistance: QM/MM modeling of the acylation reaction of a class A beta-lactamase with benzylpenicillin

Understanding the mechanisms by which beta-lactamases destroy beta-lactam antibiotics is potentially vital in developing effective therapies to overcome bacterial antibiotic resistance. Class A beta-lactamases are the most important and common type of these enzymes. A key process in the reaction mechanism of class A beta-lactamases is the acylation of the active site serine by the antibiotic. We have modeled the complete mechanism of acylation with benzylpenicillin, using a combined quantum mechanical and molecular mechanical (QM/MM) method (B3LYP/6-31G+(d)//AM1-CHARMM22). All active site residues directly involved in the reaction, and the substrate, were treated at the QM level, with reaction energies calculated at the hybrid density functional (B3LYP/6-31+Gd) level. Structures and interactions with the protein were modeled by the AM1-CHARMM22 QM/MM approach. Alternative reaction coordinates and mechanisms have been tested by calculating a number of potential energy surfaces for each step of the acylation mechanism. The results support a mechanism in which Glu166 acts as the general base. Glu166 deprotonates an intervening conserved water molecule, which in turn activates Ser70 for nucleophilic attack on the antibiotic. This formation of the tetrahedral intermediate is calculated to have the highest barrier of the chemical steps in acylation. Subsequently, the acylenzyme is formed with Ser130 as the proton donor to the antibiotic thiazolidine ring, and Lys73 as a proton shuttle residue. The presented mechanism is both structurally and energetically consistent with experimental data. The QM/MM energy barrier (B3LYP/ 6-31G+(d)//AM1-CHARMM22) for the enzymatic reaction of 9 kcal mol(-1) is consistent with the experimental activation energy of about 12 kcal mol(-1). The effects of essential catalytic residues have been investigated by decomposition analysis. The results demonstrate the importance of the "oxyanion hole" in stabilizing the transition state and the tetrahedral intermediate. In addition, Asn132 and a number of charged residues in the active site have been identified as being central to the stabilizing effect of the enzyme. These results will be potentially useful in the development of stable beta-lactam antibiotics and for the design of new inhibitors.     

Quantum mechanics/molecular mechanics investigation of the chemical reaction in Dpo4 reveals water-dependent pathways and requirements for active site reorganization

The nucleotidyl-transfer reaction coupled with the conformational transitions in DNA polymerases is critical for maintaining the fidelity and efficiency of DNA synthesis. We examine here the possible reaction pathways of a Y-family DNA polymerase, Sulfolobus solfataricus DNA polymerase IV (Dpo4), for the correct insertion of dCTP opposite 8-oxoguanine using the quantum mechanics/molecular mechanics (QM/MM) approach, both from a chemistry-competent state and a crystal closed state. The latter examination is important for understanding pre-chemistry barriers to interpret the entire enzyme mechanism, since the crystal closed state is not an ideal state for initiating the chemical reaction. The most favorable reaction path involves initial deprotonation of O3'H via two bridging water molecules to O1A, overcoming an overall potential energy barrier of approximately 20.0 kcal/mol. The proton on O1A-P(alpha) then migrates to the gamma-phosphate oxygen of the incoming nucleotide as O3' attacks P(alpha), and the P(alpha)-O3A bond breaks. The other possible pathway in which the O3'H proton is transferred directly to O1A on P(alpha) has an overall energy barrier of 25.0 kcal/mol. In both reaction paths, the rate-limiting step is the initial deprotonation, and the trigonal-bipyramidal configuration for P(alpha) occurs during the concerted bond formation (O3'-P(alpha)) and breaking (P(alpha)-O3A), indicating the associative nature of the chemical reaction. In contrast, the Dpo4/DNA complex with an imperfect active-site geometry corresponding to the crystal state must overcome a much higher activation energy barrier (29.0 kcal/mol) to achieve a tightly organized site due to hindered O3'H deprotonation stemming from larger distances and distorted conformation of the proton acceptors. This significant difference demonstrates that the pre-chemistry reorganization in Dpo4 costs approximately 4.0 to 9.0 kcal/mol depending on the primer terminus environment. Compared to the higher fidelity DNA polymerase beta from the X-family, Dpo4 has a higher chemical reaction barrier (20.0 vs 15.0 kcal/mol) due to the more solvent-exposed active site.     

Computational mutation analysis of hydrogen abstraction and radical rearrangement steps in the catalysis of coenzyme B12-dependent diol dehydratase

A mutation analysis of the catalytic functions of active-site residues of coenzyme B(12)-dependent diol dehydratase in the conversion of 1,2-propanediol to 1,1-propanediol has been carried out by using QM/MM computations. Mutants His143Ala, Glu170Gln, Glu170Ala, and Glu170Ala/Glu221Ala were considered to estimate the impact of the mutations of His143 and Glu170. In the His143Ala mutant the activation energy for OH migration increased to 16.4 from 11.5 kcal mol(-1) in the wild-type enzyme. The highest activation energy, 19.6 kcal mol(-1), was measured for hydrogen back-abstraction in this reaction. The transition state for OH migration is not sufficiently stabilized by the hydrogen-bonding interaction formed between the spectator OH group and Gln170 in the Glu170Gln mutant, which demonstrates that a strong proton acceptor is required to promote OH migration. In the Glu170Ala mutant, a new strong hydrogen bond is formed between the spectator OH group and Glu221. A computed activation energy of 13.6 kcal mol(-1) for OH migration in the Glu170Ala mutant is only 2.1 kcal mol(-1) higher than the corresponding barrier in the wild-type enzyme. Despite the low activation barrier, the Glu170Ala mutant is inactive because the subsequent hydrogen back-abstraction is energetically demanding in this mutant. OH migration is not feasible in the Glu170Ala/Glu221Ala mutant because the activation barrier for OH migration is greatly increased by the loss of COO(-) groups near the spectator OH group. This result indicates that the effect of partial deprotonation of the spectator OH group is the most important factor in reducing the activation barrier for OH migration in the conversion of 1,2-propanediol to 1,1-propanediol catalyzed by diol dehydratase.     

Towards the Accurate Thermodynamic Characterization of Enzyme Reaction Mechanisms

We employed QM/MM molecular dynamics (MD) simulations to characterize the rate-limiting step of the glycosylation reaction of pancreatic α-amylase with combined DFT/molecular dynamics methods (PBE/def2-SVP : AMBER). Upon careful choice of four starting active site conformations based on thorough reactivity criteria, Gibbs energy profiles were calculated with umbrella sampling simulations within a statistical convergence of 1–2 kcal ⋅ mol⁻¹. Nevertheless, Gibbs activation barriers and reaction energies still varied from 11.0 to 16.8 kcal ⋅ mol⁻¹ and −6.3 to +3.8 kcal ⋅ mol⁻¹ depending on the starting conformations, showing that despite significant state-of-the-art QM/MM MD sampling (0.5 ns/profile) the result still depends on the starting structure. The results supported the one step dissociative mechanism of Asp197 glycosylation preceded by an acid-base reaction by the Glu233, which are qualitatively similar to those from multi-PES QM/MM studies, and thus support the use of the latter to determine enzyme reaction mechanisms.     

Computational determination of fundamental pathway and activation barriers for acetohydroxyacid synthase‐catalyzed condensation reactions of α‐keto acids

Acetohydroxyacid synthase (AHAS) is the first common enzyme in the biosynthetic pathway leading to the production of various branched‐chain amino acids. AHAS is recognized as a promising target for new antituberculosis drugs, antibacterial drugs, and herbicides. Extensive first‐principles quantum mechanical (QM) and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations have enabled us, in this study, to uncover the fundamental reaction pathway, determine the activation barriers, and obtain valuable insights concerning the specific roles of key amino acid residues for the common steps of AHAS‐catalyzed condensation reactions of α‐keto acids. The computational results reveal that the rate‐determining step of the AHAS‐catalyzed reactions is the second reaction step and that the most important amino acid residues involved in the catalysis include Glu144′, Gln207′, Gly121′, and Gly511 that form favorable hydrogen bonds with the reaction center (consisting of atoms from the substrate and cofactor) during the reaction process. In addition, Glu144′ also accepts a proton from cofactor thiamin diphosphate (ThDP) through hydrogen bonding during the catalytic reaction. The favorable interactions between the reaction center and protein environment remarkably stabilize the transition state and, thus, lower the activation barrier for the rate‐determining reaction step by ～20 kcal/mol. The activation barrier calculated for the rate‐determining step is in good agreement with the experimental activation barrier. The detailed structural and mechanistic insights should be valuable for rational design of novel, potent AHAS inhibitors that may be used as promising new anti‐tuberculosis drugs, antibacterial drugs, and/or herbicides to overcome drug resistance problem. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods for enzymic reactions. II. An energy decomposition analysis

QM/MM methods have been developed as a computationally feasible solution to QM simulation of chemical processes, such as enzyme-catalyzed reactions, within a more approximate MM representation of the condensed-phase environment. However, there has been no independent method for checking the quality of this representation, especially for highly nonisotropic protein environments such as those surrounding enzyme active sites. Hence, the validity of QM/MM methods is largely untested. Here we use the possibility of performing all-QM calculations at the semiempirical PM3 level with a linear-scaling method (MOZYME) to assess the performance of a QM/MM method (PM3/AMBER94 force field). Using two model pathways for the hydride-ion transfer reaction of the enzyme dihydrofolate reductase studied previously (Titmuss et al., Chem Phys Lett 2000, 320, 169-176), we have analyzed the reaction energy contributions (QM, QM/MM, and MM) from the QM/MM results and compared them with analogous-region components calculated via an energy partitioning scheme implemented into MOZYME. This analysis further divided the MOZYME components into Coulomb, resonance and exchange energy terms. For the model in which the MM coordinates are kept fixed during the reaction, we find that the MOZYME and QM/MM total energy profiles agree very well, but that there are significant differences in the energy components. Most significantly there is a large change (approximately 16 kcal/mol) in the MOZYME MM component due to polarization of the MM region surrounding the active site, and which arises mostly from MM atoms close to (<10 A) the active-site QM region, which is not modelled explicitly by our QM/MM method. However, for the model where the MM coordinates are allowed to vary during the reaction, we find large differences in the MOZYME and QM/MM total energy profiles, with a discrepancy of 52 kcal/mol between the relative reaction (product-reactant) energies. This is largely due to a difference in the MM energies of 58 kcal/mol, of which we can attribute approximately 40 kcal/mol to geometry effects in the MM region and the remainder, as before, to MM region polarization. Contrary to the fixed-geometry model, there is no correlation of the MM energy changes with distance from the QM region, nor are they contributed by only a few residues. Overall, the results suggest that merely extending the size of the QM region in the QM/MM calculation is not a universal solution to the MOZYME- and QM/MM-method differences. They also suggest that attaching physical significance to MOZYME Coulomb, resonance and exchange components is problematic. Although we conclude that it would be possible to reparameterize the QM/MM force field to reproduce MOZYME energies, a better way to account for both the effects of the protein environment and known deficiencies in semiempirical methods would be to parameterize the force field based on data from DFT or ab initio QM linear-scaling calculations. Such a force field could be used efficiently in MD simulations to calculate free energies.     

Molecular dynamics and free energy perturbation study of hydride-ion transfer step in dihydrofolate reductase using combined quantum and molecular mechanical model

We used molecular dynamics simulation and free energy perturbation (FEP) methods to investigate the hydride-ion transfer step in the mechanism for the nicotinamide adenine dinucleotide phosphate (NADPH)-dependent reduction of a novel substrate by the enzyme dihydrofolate reductase (DHFR). The system is represented by a coupled quantum mechanical and molecular mechanical (QM/MM) model based on the AM1 semiempirical molecular orbital method for the reacting substrate and NADPH cofactor fragments, the AMBER force field for DHFR, and the TIP3P model for solvent water. The FEP calculations were performed for a number of choices for the QM system. The substrate, 8-methylpterin, was treated quantum mechanically in all the calculations, while the larger cofactor molecule was partitioned into various QM and MM regions with the addition of “link” atoms (F, CH₃, and H). Calculations were also carried out with the entire NADPH molecule treated by QM. The free energies of reaction and the net charges on the NADPH fragments were used to determine the most appropriate QM/MM model. The hydride-ion transfer was also carried out over several FEP pathways, and the QM and QM/MM component free energies thus calculated were found to be state functions (i.e., independent of pathway). A ca. 10 kcal/mol increase in free energy for the hydride-ion transfer with an activation barrier of ca. 30 kcal/mol was calculated. The increase in free energy on the hydride-ion transfer arose largely from the QM/MM component. Analysis of the QM/MM energy components suggests that, although a number of charged residues may contribute to the free energy change through long-range electrostatic interactions, the only interaction that can account for the 10 kcal/mol increase in free energy is the hydrogen bond between the carboxylate side chain of Glu30 (avian DHFR) and the activated (protonated) substrate. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 977–988, 1998     

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.     

Ab initio QM/MM study of class A beta-lactamase acylation: dual participation of Glu166 and Lys73 in a concerted base promotion of Ser70

Beta-lactamase acquisition is the most prevalent basis for Gram-negative bacteria resistance to the beta-lactam antibiotics. The mechanism used by the most common class A Gram-negative beta-lactamases is serine acylation followed by hydrolytic deacylation, destroying the beta-lactam. The ab initio quantum mechanical/molecular mechanical (QM/MM) calculations, augmented by extensive molecular dynamics simulations reported herein, describe the serine acylation mechanism for the class A TEM-1 beta-lactamase with penicillanic acid as substrate. Potential energy surfaces (based on approximately 350 MP2/6-31+G calculations) reveal the proton movements that govern Ser70 tetrahedral formation and then collapse to the acyl-enzyme. A remarkable duality of mechanism for tetrahedral formation is implicated. Following substrate binding, the pathway initiates by a low energy barrier (5 kcal mol(-1)) and an energetically favorable transfer of a proton from Lys73 to Glu166, through the catalytic water molecule and Ser70. This gives unprotonated Lys73 and protonated Glu166. Tetrahedral formation ensues in a concerted general base process, with Lys73 promoting Ser70 addition to the beta-lactam carbonyl. Moreover, the three-dimensional potential energy surface also shows that the previously proposed pathway, involving Glu166 as the general base promoting Ser70 through a conserved water molecule, exists in competition with the Lys73 process. The existence of two routes to the tetrahedral species is fully consistent with experimental data for mutant variants of the TEM beta-lactamase.     

Catalytic mechanism of glycosyltransferases: hybrid quantum mechanical/molecular mechanical study of the inverting N-acetylglucosaminyltransferase I

The Golgi glycosyltransferase, N-acetylglucosaminyltransferase I (GnT-I), catalyzes the transfer of a GlcNAc residue from the donor UDP-GlcNAc to the C2-hydroxyl group of a mannose residue in the trimannosyl core of the Man5GlcNAc2-Asn-X oligosaccharide. The catalytic mechanism of GnT-I was investigated using a hybrid quantum mechanical/molecular mechanical (QM/MM) method with a QM part containing 88 atoms treated with density functional theory (DFT) at the BP/TZP level. The remaining parts of a GnT-I complex, altogether 5633 atoms, were modeled using the AMBER molecular force field. A theoretical model of a Michaelis complex was built using the X-ray structure of GnT-I in complex with the donor having geometrical features consistent with kinetic studies. The QM(DFT)/MM model identified a concerted SN2-type of transition state with D291 as the catalytic base for the reaction in the enzyme active site. The TS model features nearly simultaneous nucleophilic addition and dissociation steps accompanied by the transfer of the nucleophile proton Hb2 to the catalytic base D291. The structure of the TS model is characterized by the Ob2-C1 and C1-O1 bond distances of 1.912 and 2.542 A, respectively. The activation energy for the proposed reaction mechanism was estimated to be approximately 19 kcal mol-1. The calculated alpha-deuterium kinetic isotope effect of 1.060 is consistent with the proposed reaction mechanism. Theoretical results also identified interactions between the Hb6 and beta-phosphate oxygen of the UDP and a low-barrier hydrogen bond between the nucleophile and the catalytic base D291. It is proposed that these interactions contribute to a stabilization of TS. This modeling study provided detailed insight into the mechanism of the GlcNAc transfer catalyzed by GnT-I, which is the first step in the conversion of high mannose oligosaccharides to complex and hybrid N-glycan structures.     

Convergence in the QM‐only and QM/MM modeling of enzymatic reactions: A case study for acetylene hydratase

We report systematic quantum mechanics‐only (QM‐only) and QM/molecular mechanics (MM) calculations on an enzyme‐catalyzed reaction to assess the convergence behavior of QM‐only and QM/MM energies with respect to the size of the chosen QM region. The QM and MM parts are described by density functional theory (typically B3LYP/def2‐SVP) and the CHARMM force field, respectively. Extending our previous work on acetylene hydratase with QM regions up to 157 atoms (Liao and Thiel, J. Chem. Theory Comput. 2012, 8, 3793), we performed QM/MM geometry optimizations with a QM region M4 composed of 408 atoms, as well as further QM/MM single‐point calculations with even larger QM regions up to 657 atoms. A charge deletion analysis was conducted for the previously used QM/MM model ( M3a , with a QM region of 157 atoms) to identify all MM residues with strong electrostatic contributions to the reaction energetics (typically more than 2 kcal/mol), which were then included in M4 . QM/MM calculations with this large QM region M4 lead to the same overall mechanism as the previous QM/MM calculations with M3a , but there are some variations in the relative energies of the stationary points, with a mean absolute deviation (MAD) of 2.7 kcal/mol. The energies of the two relevant transition states are close to each other at all levels applied (typically within 2 kcal/mol), with the first (second) one being rate‐limiting in the QM/MM calculations with M3a ( M4 ). QM‐only gas‐phase calculations give a very similar energy profile for QM region M4 (MAD of 1.7 kcal/mol), contrary to the situation for M3a where we had previously found significant discrepancies between the QM‐only and QM/MM results (MAD of 7.9 kcal/mol). Extension of the QM region beyond M4 up to M7 (657 atoms) leads to only rather small variations in the relative energies from single‐point QM‐only and QM/MM calculations (MAD typically about 1–2 kcal/mol). In the case of acetylene hydratase, a model with 408 QM atoms thus seems sufficient to achieve convergence in the computed relative energies to within 1–2 kcal/mol.Copyright © 2013 Wiley Periodicals, Inc.     

Combined QM(DFT)/MM molecular dynamics simulations of the deamination of cytosine by yeast cytosine deaminase (yCD)

Extensive combined quantum mechanical (B3LYP/6‐31G*) and molecular mechanical (QM/MM) molecular dynamics simulations have been performed to elucidate the hydrolytic deamination mechanism of cytosine to uracil catalyzed by the yeast cytosine deaminase (yCD). Though cytosine has no direct binding to the zinc center, it reacts with the water molecule coordinated to zinc, and the adjacent conserved Glu64 serves as a general acid/base to shuttle protons from water to cytosine. The overall reaction consists of several proton‐transfer processes and nucleophilic attacks. A tetrahedral intermediate adduct of cytosine and water binding to zinc is identified and similar to the crystal structure of yCD with the inhibitor 2‐pyrimidinone. The rate‐determining step with the barrier of 18.0 kcal/mol in the whole catalytic cycle occurs in the process of uracil departure where the proton transfer from water to Glu64 and nucleophilic attack of the resulting hydroxide anion to C2 of the uracil ring occurs synchronously. © 2016 Wiley Periodicals, Inc.     

A quantum mechanical investigation of possible mechanisms for the nucleotidyl transfer reaction catalyzed by DNA polymerase beta

Several quantum mechanical (QM) and hybrid quantum/molecular mechanical (QM/MM) studies have been employed recently to analyze the nucleotidyl transfer reaction in DNA polymerase beta (pol beta). Our examination reveals strong dependence of the reported mechanism on the initial molecular model. Thus, we explore here several model systems by QM methods to investigate pol beta's possible pathway variations. Although our most favorable pathway involves a direct proton transfer from O3'(primer) to O2alpha(Palpha), we also discuss other initial proton-transfer steps--to an adjacent water, to triphosphate, or to aspartic units--and the stabilizing effect of crystallographic water molecules in the active site. Our favored reaction route has an energetically undemanding initial step of less than 1.0 kcal/mol (at the B3LYP/6-31G(d,p) level), and involves a slight rearrangement in the geometry of the active site. This is followed by two major steps: (1) direct proton transfer from O3'(primer) to O2alpha(Palpha) leading to the formation of a pentavalent, trigonal bipyramidal Palpha center, via an associative mechanism, at a cost of about 28 kcal/mol, and (2) breakage of the triphosphate unit (exothermic process, approximately 22 kcal/mol) that results in the full transfer of the nucleotide to the DNA and the formation of pyrophosphate. These energy values are expected to be lower in the physical system when full protein effects are incorporated. We also discuss variations from this dominant pathway, and their impact on the overall repair process. Our calculated barrier for the chemical reaction clearly indicates that chemistry is rate-limiting overall for correct nucleotide insertion in pol beta, in accord with other studies. Protonation studies on relevant intermediates suggest that, although protonation at a single aspartic residue may occur, the addition of a second proton to the system significantly disturbs the active site. We conclude that the active site rearrangement step necessary to attain a reaction-competent geometry is essential and closely related to the "pre-chemistry" avenue described recently as a key step in the overall kinetic cycle of DNA polymerases. Thus, our work emphasizes the many possible ways for DNA polymerase beta's chemical reaction to occur, determined by the active site environment and initial models.     

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

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