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

Effects of the Y432S Cancer-Associated Variant on the Reaction Mechanism of Human DNA Polymerase κ

Human DNA polymerases are vital for genetic information management. Their function involves catalyzing the synthesis of DNA strands with unparalleled accuracy, which ensures the fidelity and stability of the human genomic blueprint. Several disease-associated mutations and their functional impact on DNA polymerases have been reported. One particular polymerase, human DNA polymerase kappa (Pol κ), has been reported to be susceptible to several cancer-associated mutations. The Y432S mutation in Pol κ, associated with various cancers, is of interest due to its impact on polymerization activity and markedly reduced thermal stability. Here, we have used computational simulations to investigate the functional consequences of the Y432S using classical molecular dynamics (MD) and coupled quantum mechanics/molecular mechanics (QM/MM) methods. Our findings suggest that Y432S induces structural alterations in domains responsible for nucleotide addition and ternary complex stabilization while retaining structural features consistent with possible catalysis in the active site. Calculations of the minimum energy path associated with the reaction mechanism of the wild type (WT) and Y432S Pol κ indicate that, while both enzymes are catalytically competent (in terms of energetics and the active site's geometries), the cancer mutation results in an endoergic reaction and an increase in the catalytic barrier. Interactions with a third magnesium ion and environmental effects on nonbonded interactions, particularly involving key residues, contribute to the kinetic and thermodynamic distinctions between the WT and mutant during the catalytic reaction. The energetics and electronic findings suggest that active site residues favor the catalytic reaction with dCTP3- over dCTP4-.

Mechanistic Insights into the Halophilic Xylosidase Xylo-1 and Its Role in Xylose Production

The Xylo-1 xylosidase, which belongs to the GH43 family, exhibits a high salt tolerance. The present study demonstrated that the catalytic activity of Xylo-1 increased by 195% in the presence of 5 M NaCl. Additionally, the half-life of Xylo-1 increased 25.9-fold in the presence of 1 M NaCl. Through comprehensive analysis including circular dichroism, fluorescence spectroscopy, and molecular dynamics simulations, we elucidated that the presence of Na+ ions increased the contact frequency between the surface acidic amino acids and the surrounding water molecules. This resulted in the stabilization of the surrounding hydration layer of Xylo-1. Additionally, Na+ ions also stabilized the substrate-binding conformation and the fluctuation of water molecules within the active site, which enhanced the catalytic activity of Xylo-1 by increasing the nucleophilic attack by the water molecules. Ultimately, the optimal reaction conditions for the production of xylose by synergistic catalysis with Xylo-1 and xylanase were determined. The results demonstrated that the conversion yield of the method was high for various sources of xylan, indicating the method could have potential industrial applications. This study explored the structure-activity relationship of catalysis in Xylo-1 under high-salt conditions, provides novel insights into the mechanism of halophilic enzymes, and serves as a reference for the industrial application of Xylo-1.

Multiple deprotonation paths of the nucleophile 3'-OH in the DNA synthesis reaction

DNA synthesis by polymerases is essential for life. Deprotonation of the nucleophile 3'-OH is thought to be the obligatory first step in the DNA synthesis reaction. We have examined each entity surrounding the nucleophile 3'-OH in the reaction catalyzed by human DNA polymerase (Pol) η and delineated the deprotonation process by combining mutagenesis with steady-state kinetics, high-resolution structures of in crystallo reactions, and molecular dynamics simulations. The conserved S113 residue, which forms a hydrogen bond with the primer 3'-OH in the ground state, stabilizes the primer end in the active site. Mutation of S113 to alanine destabilizes primer binding and reduces the catalytic efficiency. Displacement of a water molecule that is hydrogen bonded to the 3'-OH using the 2'-OH of a ribonucleotide or 2'-F has little effect on catalysis. Moreover, combining the S113A mutation with 2'-F replacement, which removes two potential hydrogen acceptors of the 3'-OH, does not reduce the catalytic efficiency. We conclude that the proton can leave the O3' via alternative paths, supporting the hypothesis that binding of the third Mg²⁺ initiates the reaction by breaking the α-β phosphodiester bond of an incoming deoxyribonucleoside triphosphate (dNTP).

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Uncovering a unique approach for damaged DNA replication: A computational investigation of a mutagenic tobacco-derived thymine lesion

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone is a potent nicotine carcinogen that leads to many DNA lesions, the most persistent being the O2-[4-oxo-4-(3-pyridyl)butyl]thymine adduct (POB-T). Although the experimental mutagenic profile for the minor groove POB-T lesion has been previously reported, the findings are puzzling in terms of the human polymerases involved. Specifically, while pol κ typically replicates minor groove adducts, in vivo studies indicate pol η replicates POB-T despite being known for processing major groove adducts. Our multiscale modeling approach reveals that the canonical (anti) glycosidic orientation of POB-T can fit in the pol κ active site, but only a unique (syn) POB-T conformation is accommodated by pol η. These distinct binding orientations rationalize the differential in vitro mutagenic spectra based on the preferential stabilization of dGTP and dTTP opposite the lesion for pol κ and η, respectively. Overall, by uncovering the first evidence for the replication of a damaged pyrimidine in the syn glycosidic orientation, the current work provides the insight necessary to clarify a discrepancy in the DNA replication literature, expand the biological role of the critical human pol η, and understand the mutational signature in human cancers associated with tobacco exposure.

Computational insights into the mutagenicity of two tobacco-derived carcinogenic DNA lesions

4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone is a potent carcinogen found in all tobacco products that leads to a variety of DNA lesions in cells, including O6-[4-oxo-4-(3-pyridyl)butyl]guanine (POB-G) and O6-[4-hydroxy-4-(3-pyridyl)butyl]guanine (PHB-G), which differ by only a single substituent in the bulky moiety. This work uses a multiscale computational approach to shed light on the intrinsic conformational and base-pairing preferences of POB-G and PHB-G, and the corresponding properties in DNA and the polymerase η active site. Our calculations reveal that both lesions form stable pairs with C and T, with the T pairs being the least distorted relative to canonical DNA. This rationalizes the experimentally reported mutational profile for POB-G and validates our computational model. The same approach predicts that PHB-G is more mutagenic than POB-G due to a difference in the bulky moiety hydrogen-bonding pattern, which increases the stability of the PHB-G:T pair. The mutagenicity of PHB-G is likely further increased by stabilization of an intercalated DNA conformation that is associated with deletion mutations. This work thereby uncovers structural explanations for the reported mutagenicity of POB-G, provides the first clues regarding the mutagenicity of PHB-G and complements a growing body of literature highlighting that subtle chemical changes can affect the biological outcomes of DNA adducts.

Computational Simulations of DNA Polymerases: Detailed Insights on Structure/Function/Mechanism from Native Proteins to Cancer Variants

Genetic information is vital in the cell cycle of DNA-based organisms. DNA polymerases (DNA Pols) are crucial players in transactions dealing with these processes. Therefore, the detailed understanding of the structure, function, and mechanism of these proteins has been the focus of significant effort. Computational simulations have been applied to investigate various facets of DNA polymerase structure and function. These simulations have provided significant insights over the years. This perspective presents the results of various computational studies that have been employed to research different aspects of DNA polymerases including detailed reaction mechanism investigation, mutagenicity of different metal cations, possible factors for fidelity synthesis, and discovery/functional characterization of cancer-related mutations on DNA polymerases.

Molecular Insights into the Translesion Synthesis of Benzyl-Guanine from Molecular Dynamics Simulations: Structural Evidence of Mutagenic and Nonmutagenic Replication

DNA can be damaged by many compounds in our environment, and the resulting damaged DNA is commonly replicated by translesion synthesis (TLS) polymerases. Because the mechanism and efficiency of TLS are affected by the type of DNA damage, obtaining information for a variety of DNA adducts is critical. However, there is no structural information for the insertion of a dNTP opposite an O6-dG adduct, which is a particularly harmful class of DNA lesions. We used molecular dynamics (MD) simulations to investigate structural and energetic parameters that dictate preferred dNTP insertion opposite O6-benzyl-guanine (Bz-dG) by DNA polymerase IV, a prototypical TLS polymerase. Specifically, MD simulations were completed on all possible ternary insertion complexes and ternary -1 base deletion complexes with different Bz-dG conformations. Our data suggests that the purines are unlikely to be inserted opposite anti- or syn-Bz-dG, and dTTP is unlikely to be inserted opposite syn-Bz-dG, because of changes in the active site conformation, including critical hydrogen-bonding interactions and/or reaction-ready parameters compared to natural dG replication. In contrast, a preserved active site conformation suggests that dCTP can be inserted opposite either anti- or syn-Bz-dG and dTTP can be inserted opposite anti-Bz-dG. This is the first structural explanation for the experimentally observed preferential insertion of dCTP and misincorporation of dTTP opposite Bz-dG. Furthermore, we provide atomic level insight into why Bz-dG replication does not lead to deletion mutations, which is in contrast with the replication outcomes of other adducts. These findings provide a basis for understanding the replication of related O6-dG adducts.

A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases

The enzymatic polymerization of DNA and RNA is the basis for genetic inheritance for all living organisms. It is catalyzed by the DNA/RNA polymerase (Pol) superfamily. Here, bioinformatics analysis reveals that the incoming nucleotide substrate always forms an H-bond between its 3'-OH and β-phosphate moieties upon formation of the Michaelis complex. This previously unrecognized H-bond implies a novel self-activated mechanism (SAM), which synergistically connects the in situ nucleophile formation with subsequent nucleotide addition and, importantly, nucleic acid translocation. Thus, SAM allows an elegant and efficient closed-loop sequence of chemical and physical steps for Pol catalysis. This is markedly different from previous mechanistic hypotheses. Our proposed mechanism is corroborated via ab initio QM/MM simulations on a specific Pol, the human DNA polymerase-η, an enzyme involved in repairing damaged DNA. The structural conservation of DNA and RNA Pols supports the possible extension of SAM to Pol enzymes from the three domains of life.

Exploring the mechanism of DNA polymerases by analyzing the effect of mutations of active site acidic groups in Polymerase β

Elucidating the catalytic mechanism of DNA polymerase is crucial for a progress in the understanding of the control of replication fidelity. This work tries to advance the mechanistic understanding by analyzing the observed effect of mutations of the acidic groups in the active site of Polymerase β as well as the pH effect on the rate constant. The analysis involves both empirical valence bond (EVB) free energy calculations and considerations of the observed pH dependence of the reaction. The combined analysis indicates that the proton transfer (PT) from the nucleophilic O3' has two possible pathways, one to D256 and the second to the bulk. We concluded based on calculations and the experimental pH profile that the most likely path for the wild-type (WT) and the D256E and D256A mutants is a PT to the bulk, although the WT may also use a PT to Asp 256. Our analysis highlights the need for very extensive sampling in the calculations of the activation barrier and also clearly shows that ab initio QM/MM calculations that do not involve extensive sampling are unlikely to give a clear quantitative picture of the reaction mechanism. Proteins 2016; 84:1644-1657. © 2016 Wiley Periodicals, Inc.

Computational Evaluation of Nucleotide Insertion Opposite Expanded and Widened DNA by the Translesion Synthesis Polymerase Dpo4

Expanded (x) and widened (y) deoxyribose nucleic acids (DNA) have an extra benzene ring incorporated either horizontally (xDNA) or vertically (yDNA) between a natural pyrimidine base and the deoxyribose, or between the 5- and 6-membered rings of a natural purine. Far-reaching applications for (x,y)DNA include nucleic acid probes and extending the natural genetic code. Since modified nucleobases must encode information that can be passed to the next generation in order to be a useful extension of the genetic code, the ability of translesion (bypass) polymerases to replicate modified bases is an active area of research. The common model bypass polymerase DNA polymerase IV (Dpo4) has been previously shown to successfully replicate and extend past a single modified nucleobase on a template DNA strand. In the current study, molecular dynamics (MD) simulations are used to evaluate the accommodation of expanded/widened nucleobases in the Dpo4 active site, providing the first structural information on the replication of (x,y)DNA. Our results indicate that the Dpo4 catalytic (palm) domain is not significantly impacted by the (x,y)DNA bases. Instead, the template strand is displaced to accommodate the increased C1’–C1’ base-pair distance. The structural insights unveiled in the present work not only increase our fundamental understanding of Dpo4 replication, but also reveal the process by which Dpo4 replicates (x,y)DNA, and thereby will contribute to the optimization of high fidelity and efficient polymerases for the replication of modified nucleobases.

The Closing Mechanism of DNA Polymerase I at Atomic Resolution

DNA polymerases must quickly and accurately distinguish between similar nucleic acids to form Watson-Crick base pairs and avoid DNA replication errors. Deoxynucleoside triphosphate (dNTP) binding to the DNA polymerase active site induces a large conformational change that is difficult to characterize experimentally on an atomic level. Here, we report an X-ray crystal structure of DNA polymerase I bound to DNA in the open conformation with a dNTP present in the active site. We use this structure to computationally simulate the open to closed transition of DNA polymerase in the presence of a Watson-Crick base pair. Our microsecond simulations allowed us to characterize the key steps involved in active site assembly, and propose the sequence of events involved in the prechemistry steps of DNA polymerase catalysis. They also reveal new features of the polymerase mechanism, such as a conserved histidine as a potential proton acceptor from the primer 3'-hydroxyl.

Applications of quantum mechanical/molecular mechanical methods to the chemical insertion step of DNA and RNA polymerization

We review theoretical attempts to model the chemical insertion reactions of nucleoside triphosphates catalyzed by the nucleic acid polymerases using combined quantum mechanical/molecular mechanical methodology. Due to an existing excellent database of high-resolution X-ray crystal structures, the DNA polymerase β system serves as a useful template for discussion and comparison. The convergence of structures of high-quality complexes and continued developments of theoretical techniques suggest a bright future for understanding the global features of nucleic acid polymerization.

Ultrafast water dynamics at the interface of the polymerase-DNA binding complex

DNA polymerases slide on DNA during replication, and the interface must be mobile for various conformational changes. The role of lubricant interfacial water is not understood. In this report, we systematically characterize the water dynamics at the interface and in the active site of a tight binding polymerase (pol β) in its binary complex and ternary state using tryptophan as a local optical probe. Using femtosecond spectroscopy, we observed that upon DNA recognition the surface hydration water is significantly confined and becomes bound water at the interface, but the dynamics are still ultrafast and occur on the picosecond time scale. These interfacial water molecules are not trapped but are mobile in the heterogeneous binding nanospace. Combining our findings with our previous observation of ultrafast water motions at the interface of a loose binding polymerase (Dpo4), we conclude that the binding interface is dynamic and the water molecules in various binding clefts, channels, and caves are mobile and even fluid with different levels of mobility for loose or tight binding polymerases. Such a dynamic interface should be general to all DNA polymerase complexes to ensure the biological function of DNA synthesis.

Direct probing of solvent accessibility and mobility at the binding interface of polymerase (Dpo4)-DNA complex

Water plays essential structural and dynamical roles in protein-DNA recognition through contributing to enthalpic or entropic stabilization of binding complex and by mediating intermolecular interactions and fluctuations for biological function. These interfacial water molecules are confined by the binding partners in nanospace, but in many cases they are highly mobile and exchange with outside bulk solution. Here, we report our studies of the interfacial water dynamics in the binary and ternary complexes of a polymerase (Dpo4) with DNA and an incoming nucleotide using a site-specific tryptophan probe with femtosecond resolution. By systematic comparison of the interfacial water motions and local side chain fluctuations in the apo, binary, and ternary states of Dpo4, we observed that the DNA binding interface and active site are dynamically solvent accessible and the interfacial water dynamics are similar to the surface hydration water fluctuations on picosecond time scales. Our molecular dynamics simulations also show the binding interface full of water molecules and nonspecific weak interactions. Such a fluid binding interface facilitates the polymerase sliding on DNA for fast translocation whereas the spacious and mobile hydrated active site contributes to the low fidelity of the lesion-bypass Y-family DNA polymerase.

Substrate-mediated proton relay mechanism for the religation reaction in topoisomerase II

The DNA religation reaction of yeast type II topoisomerase (topo II) was investigated to elucidate its metal-dependent general acid/base catalysis. Quantum mechanical/molecular mechanical calculations were performed for the topo II religation reaction, and the proton transfer pathway was examined. We found a substrate-mediated proton transfer of the topo II religation reaction, which involves the 3' OH nucleophile, the reactive phosphate, water, Arg781, and Tyr782. Metal A stabilizes the transition states, which is consistent with a two-metal mechanism in topo II. This pathway may be required for the cleavage/religation reaction of topo IA and II and will provide a general explanation for the catalytic mechanism in the topo IA and II.

Structures of intermediates along the catalytic cycle of terminal deoxynucleotidyltransferase: dynamical aspects of the two-metal ion mechanism

Terminal deoxynucleotidyltransferase (Tdt) is a non-templated eukaryotic DNA polymerase of the polX family that is responsible for the random addition of nucleotides at the V(D)J junctions of immunoglobulins and T-cell receptors. Here we describe a series of high-resolution X-ray structures that mimic the pre-catalytic state, the post-catalytic state and a competent state that can be transformed into the two other ones in crystallo via the addition of dAMPcPP and Zn(2+), respectively. We examined the effect of Mn(2+), Co(2+) and Zn(2+) because they all have a marked influence on the kinetics of the reaction. We demonstrate a dynamic role of divalent transition metal ions bound to site A: (i) Zn(2+) (or Co(2+)) in Metal A site changes coordination from octahedral to tetrahedral after the chemical step, which explains the known higher affinity of Tdt for the primer strand when these ions are present, and (ii) metal A has to leave to allow the translocation of the primer strand and to clear the active site, a typical feature for a ratchet-like mechanism. Except for Zn(2+), the sugar puckering of the primer strand 3' terminus changes from C2'-endo to C3'-endo during catalysis. In addition, our data are compatible with a scheme where metal A is the last component that binds to the active site to complete its productive assembly, as already inferred in human pol beta. The new structures have potential implications for modeling pol mu, a closely related polX implicated in the repair of DNA double-strand breaks, in a complex with a DNA synapsis.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

Addressing open questions about phosphate hydrolysis pathways by careful free energy mapping

The nature and mechanism of phosphate hydrolysis reactions are of great interest in view of the crucial role of these reactions in key biological processes. Although it is becoming clearer that the ultimate way of resolving mechanistic controversies must involve reliable theoretical studies, it is not widely realized that such studies cannot be performed at present by using most existing automated ways and that only careful systematic studies can lead to meaningful conclusions. The present work clarifies the above point by considering the hydrolysis of phosphate monoesters. The clarification starts by defining the actual issues that should be addressed in careful studies and by highlighting the problems with studies that ignore the need for unique mechanistic definitions (e.g., works that confuse associative and dissociative pathways). We then focus on the analysis of the proton transfer (PT) pathways in phosphate hydrolysis and on recent suggestions that PT involves more than one water molecule. Here we point out that most of the studies that found a proton transfer through several water molecules have not involved a sufficient systematic search of the relevant reaction coordinates. This includes both energy minimization approaches as well as a recent metadynamics (MTD) simulation study. To illustrate the crucial need of exploring the potential surfaces reliably, rather than relying on automated approaches, we present here a very careful study of the free energy landscape along a 3D reaction coordinate (RC) exploring both the standard 2D RC, comprised of the attacking and leaving group reaction coordinates, as well as of the proton transfer (PT) coordinate. Our study points out that QM/MM minimization or MTD studies that concluded that the hydrolysis of phosphate monoesters involves a PT through several water molecules, have not explored carefully the single water (1W) path (that involves a direct PT form the attacking water molecule to the phosphate oxygen). Furthermore, we identified the most likely reason for the difficulty in finding the 1W path by QM/MM minimization methods, as well as by the current MTD simulations. We also discuss the problems with current studies that challenge the phosphate as a base mechanism and emphasize that all recent studies found associative/concerted paths (although many have not realized the meaning of their results). Finally, although we clearly do not have the last word about the 1W versus 2W paths we believe that we illustrated that the crucial mechanistic problems with alternative pathways should not be resolved by just running black box search approaches.

Perspective: pre-chemistry conformational changes in DNA polymerase mechanisms

In recent papers, there has been a lively exchange concerning theories for enzyme catalysis, especially the role of protein dynamics/pre-chemistry conformational changes in the catalytic cycle of enzymes. Of particular interest is the notion that substrate-induced conformational changes that assemble the polymerase active site prior to chemistry are required for DNA synthesis and impact fidelity (i.e., substrate specificity). High-resolution crystal structures of DNA polymerase β representing intermediates of substrate complexes prior to the chemical step are available. These structures indicate that conformational adjustments in both the protein and substrates must occur to achieve the requisite geometry of the reactive participants for catalysis. We discuss computational and kinetic methods to examine possible conformational change pathways that lead from the observed crystal structure intermediates to the final structures poised for chemistry. The results, as well as kinetic data from site-directed mutagenesis studies, are consistent with models requiring pre-chemistry conformational adjustments in order to achieve high fidelity DNA synthesis. Thus, substrate-induced conformational changes that assemble the polymerase active site prior to chemistry contribute to DNA synthesis even when they do not represent actual rate-determining steps for chemistry.

Preferred WMSA catalytic mechanism of the nucleotidyl transfer reaction in human DNA polymerase κ elucidates error-free bypass of a bulky DNA lesion

Human DNA Pol κ is a polymerase enzyme, specialized for near error-free bypass of certain bulky chemical lesions to DNA that are derived from environmental carcinogens present in tobacco smoke, automobile exhaust and cooked food. By employing ab initio QM/MM-MD (Quantum Mechanics/Molecular Mechanics-Molecular Dynamics) simulations with umbrella sampling, we have determined the entire free energy profile of the nucleotidyl transfer reaction catalyzed by Pol κ and provided detailed mechanistic insights. Our results show that a variant of the Water Mediated and Substrate Assisted (WMSA) mechanism that we previously deduced for Dpo4 and T7 DNA polymerases is preferred for Pol κ as well, suggesting its broad applicability. The hydrogen on the 3'-OH primer terminus is transferred through crystal and solvent waters to the γ-phosphate of the dNTP, followed by the associative nucleotidyl transfer reaction; this is facilitated by a proton transfer from the γ-phosphate to the α,β-bridging oxygen as pyrophosphate leaves, to neutralize the evolving negative charge. MD simulations show that the near error-free incorporation of dCTP opposite the major benzo[a]pyrene-derived dG lesion is compatible with the WMSA mechanism, allowing for an essentially undisturbed pentacovalent phosphorane transition state, and explaining the bypass of this lesion with little mutation by Pol κ.

The structural role of Mg2+ ions in a class I RNA polymerase ribozyme: a molecular simulation study

According to the RNA world hypothesis, self-replicating ribozymes, storing the genetic information and being able to perform catalysis, were the constituents of the first living organisms. In particular, RNA polymerase ribozymes, similar to current proteinaceous enzymatic polymerases, may have been able to promote the synthesis of RNA strands in a primitive world. Polymerase catalysis is usually assisted by Mg(2+) ions, but it is not always trivial to find out experimentally the number of Mg(2+) ions placed in the active site as well as the identity and the number of their coordination ligands. Here, we addressed this issue in an artificial class I ligase ribozyme. On the basis of a recently solved crystal structure, we constructed computational models of reactant and product states of this ribozyme, considering monometallic and bimetallic species. Our models were relaxed by force field based molecular dynamics (MD) simulations and mixed quantum-classical (QM/MM) MD. The structural and dynamical properties of these models were consistent with experimental data and were validated by a comparison with the catalytic sites of proteinaceous DNA and RNA polymerases. Consistently with enzymatic polymerases, our results suggest that class I RNA ligases most probably contain two magnesium ions in the active site and they may, therefore, catalyze the junction of two RNA strands via "a two Mg(2+) ions" mechanism.

Prechemistry versus preorganization in DNA replication fidelity

The molecular origin of nucleotide insertion catalysis and fidelity of DNA polymerases is explored by means of computational simulations. Special attention is paid to the examination of the validity of proposals that invoke prechemistry effects, checkpoints concepts, and dynamical effects. The simulations reproduce the observed fidelity in Pol β, starting with the relevant observed X-ray structures of the complex with the right (R) and wrong (W) nucleotides. The generation of free energy surfaces for the R and W systems also allowed us to analyze different proposals about the origin of the fidelity and to reach several important conclusions. It is found that the potential of mean force (PMF) obtained by proper sampling does not support QM/MM-based proposals of a large barrier before the prechemistry state. Furthermore, examination of dynamical proposals by the renormalization approach indicates that the motions from open to close configurations do not contribute to catalysis or fidelity. Finally we discuss and analyze the induced fit concept and show that, despite its importance, it does not explain fidelity. That is, the fidelity is apparently due to the change in the preorganization of the chemical site, as a result of the relaxation of the binding site upon binding of the incorrect nucleotide. Finally and importantly, since the issue is the barrier associated with the enzyme-substrate (ES)/DNA complex at the chemical transition state and not the path to this complex formation (unless this path involves rate determining steps), it is also not useful to invoke checkpoints while discussing fidelity.

The glmS ribozyme: use of a small molecule coenzyme by a gene-regulatory RNA

The glmS ribozyme is the first known example of a natural ribozyme that has evolved to require binding of an exogenous small molecule for activity. In Gram-positive bacteria, this RNA domain is part of the messenger RNA (mRNA) encoding the essential enzyme that synthesizes glucosamine-6-phosphate (GlcN6P). When present at physiologic concentration, this small molecule binds to the glmS ribozyme and uncovers a latent self-cleavage activity that ultimately leads to degradation of the mRNA. Biochemical and structural studies reveal that the RNA adopts a rigid fold stabilized by three pseudoknots and the packing of a peripheral domain against the ribozyme core. GlcN6P binding to this pre-organized RNA does not induce conformational changes; rather, the small molecule functions as a coenzyme, providing a catalytically essential amine group to the active site. The ribozyme is not a passive player, however. Active site functional groups are essential for catalysis, even in the presence of GlcN6P. In addition to being a superb experimental system with which to analyze how RNA catalysts can exploit small molecule coenzymes to broaden their chemical versatility, the presence of the glmS ribozyme in numerous pathogenic bacteria make this RNA an attractive target for the development of new antibiotics and antibacterial strategies.

Computational delineation of the catalytic step of a high-fidelity DNA polymerase

The Bacillus fragment, belonging to a class of high-fidelity polymerases, demonstrates high processivity (adding approximately 115 bases per DNA binding event) and exceptional accuracy (1 error in 10(6) nucleotide incorporations) during DNA replication. We present analysis of structural rearrangements and energetics just before and during the chemical step (phosphodiester bond formation) using a combination of classical molecular dynamics, mixed quantum mechanics molecular mechanics simulations, and free energy computations. We find that the reaction is associative, proceeding via the two-metal-ion mechanism, and requiring the proton on the terminal primer O3' to transfer to the pyrophosphate tail of the incoming nucleotide before the formation of the pentacovalent transition state. Different protonation states for key active site residues direct the system to alternative pathways of catalysis and we estimate a free energy barrier of approximately 12 kcal/mol for the chemical step. We propose that the protonation of a highly conserved catalytic aspartic acid residue is essential for the high processivity demonstrated by the enzyme and suggest that global motions could be part of the reaction free energy landscape.

Conformational coupling, bridge helix dynamics and active site dehydration in catalysis by RNA polymerase

Molecular dynamics simulation of Thermus thermophilus (Tt) RNA polymerase (RNAP) in a catalytic conformation demonstrates that the active site dNMP-NTP base pair must be substantially dehydrated to support full active site closing and optimum conditions for phosphodiester bond synthesis. In silico mutant beta R428A RNAP, which was designed based on substitutions at the homologous position (Rpb2 R512) of Saccharomyces cerevisiae (Sc) RNAP II, was used as a reference structure to compare to Tt RNAP in simulations. Long range conformational coupling linking a dynamic segment of the bridge alpha-helix, the extended fork loop, the active site, and the trigger loop-trigger helix is apparent and adversely affected in beta R428A RNAP. Furthermore, bridge helix bending is detected in the catalytic structure, indicating that bridge helix dynamics may regulate phosphodiester bond synthesis as well as translocation. An active site "latch" assembly that includes a key trigger helix residue Tt beta' H1242 and highly conserved active site residues beta E445 and R557 appears to help regulate active site hydration/dehydration. The potential relevance of these observations in understanding RNAP and DNAP induced fit and fidelity is discussed.

At the dawn of the 21st century: Is dynamics the missing link for understanding enzyme catalysis?

Enzymes play a key role in almost all biological processes, accelerating a variety of metabolic reactions as well as controlling energy transduction, the transcription, and translation of genetic information, and signaling. They possess the remarkable capacity to accelerate reactions by many orders of magnitude compared to their uncatalyzed counterparts, making feasible crucial processes that would otherwise not occur on biologically relevant timescales. Thus, there is broad interest in understanding the catalytic power of enzymes on a molecular level. Several proposals have been put forward to try to explain this phenomenon, and one that has rapidly gained momentum in recent years is the idea that enzyme dynamics somehow contributes to catalysis. This review examines the dynamical proposal in a critical way, considering basically all reasonable definitions, including (but not limited to) such proposed effects as "coupling between conformational and chemical motions," "landscape searches" and "entropy funnels." It is shown that none of these proposed effects have been experimentally demonstrated to contribute to catalysis, nor are they supported by consistent theoretical studies. On the other hand, it is clarified that careful simulation studies have excluded most (if not all) dynamical proposals. This review places significant emphasis on clarifying the role of logical definitions of different catalytic proposals, and on the need for a clear formulation in terms of the assumed potential surface and reaction coordinate. Finally, it is pointed out that electrostatic preorganization actually accounts for the observed catalytic effects of enzymes, through the corresponding changes in the activation free energies.

The EVB as a quantitative tool for formulating simulations and analyzing biological and chemical reactions

Recent years have seen dramatic improvements in computer power, allowing ever more challenging problems to be approached. In light of this, it is imperative to have a quantitative model for examining chemical reactivity, both in the condensed phase and in solution, as well as to accurately quantify physical organic chemistry (particularly as experimental approaches can often be inconclusive). Similarly, computational approaches allow for great progress in studying enzyme catalysis, as they allow for the separation of the relevant energy contributions to catalysis. Due to the complexity of the problems that need addressing, there is a need for an approach that can combine reliability with an ability to capture complex systems in order to resolve long-standing controversies in a unique way. Herein, we will demonstrate that the empirical valence bond (EVB) approach provides a powerful way to connect the classical concepts of physical organic chemistry to the actual energies of enzymatic reactions by means of computation. Additionally, we will discuss the proliferation of this approach, as well as attempts to capture its basic chemistry and repackage it under different names. We believe that the EVB approach is the most powerful tool that is currently available for studies of chemical processes in the condensed phase in general and enzymes in particular, particularly when trying to explore the different proposals about the origin of the catalytic power of enzymes.

Polymerase-tailored variations in the water-mediated and substrate-assisted mechanism for nucleotidyl transfer: insights from a study of T7 DNA polymerase

The nucleotidyl transfer reaction catalyzed by DNA polymerases is the critical step governing the accurate transfer of genetic information during DNA replication, and its malfunctioning can cause mutations leading to human diseases, including cancer. Here, utilizing ab initio quantum mechanical/molecular mechanical calculations with free-energy perturbation, we carried out an extensive investigation of the nucleotidyl transfer reaction mechanism in the well-characterized high-fidelity replicative DNA polymerase from phage T7. Our defined mechanism entails an initial concerted deprotonation of a conserved crystal water molecule with protonation of the gamma-phosphate of the deoxynucleotide triphosphate(dNTP) via a solvent water molecule, and then the proton on the primer 3'-terminus is transferred to the resulting hydroxide ion. Subsequently, the nucleophilic attack takes place, with the formation of a metastable pentacovalent phosphorane intermediate. Finally, the pyrophosphate leaves, facilitated by the relay of the proton on the gamma-phosphate to the alpha-beta bridging oxygen via solvent water. The computed activation free-energy barrier is consistent with kinetic data for the chemistry step with correct nucleotide incorporation in T7 DNA polymerase. This variant of the water-mediated and substrate-assisted mechanism has features tailored to the structure of the T7 DNA polymerase. However, a unifying theme in the water-mediated and substrate-assisted mechanism is the cycling through crystal and solvent water molecules of the proton originating from the primer 3'-terminus to the alpha-beta bridging oxygen of the deoxynucleotide triphosphate; this neutralizes the evolving negative charge as pyrophosphate leaves and restores the polymerase to its pre-chemistry state. These unifying features are likely requisite elements for nucleotidyl transfer reactions.