1
|
Yu M, Shen Z, Zhang S, Zhang Y, Zhao H, Zhang L. The active components of Erzhi wan and their anti-Alzheimer's disease mechanisms determined by an integrative approach of network pharmacology, bioinformatics, molecular docking, and molecular dynamics simulation. Heliyon 2024; 10:e33761. [PMID: 39027618 PMCID: PMC11255520 DOI: 10.1016/j.heliyon.2024.e33761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 06/21/2024] [Accepted: 06/26/2024] [Indexed: 07/20/2024] Open
Abstract
Erzhi Wan (EZW), a classic Traditional Chinese Medicine formula, has shown promise as a potential therapeutic option for Alzheimer's disease (AD), yet its mechanism remains elusive. Herein, we employed an integrative in-silico approach to investigate the active components and their mechanisms against AD. We screened four active components with blood-brain barrier permeabilities from TCMSP, along with 307 corresponding targets predicted by SwissTargetPrediction, PharmMapper, and TCMbank websites. Then, we retrieved 2260 AD-related targets from Genecards, OMIM, and NCBI databases. Furthermore, we constructed the protein-protein interaction (PPI) network of the intersected targets via the STRING database and performed the GO and KEGG enrichment analyses using the "clusterProfiler" R package. The results showed that the intersected targets were intimately related to the p53/PI3K/Akt signaling pathway, serotonergic synapse, and response to oxygen level. Subsequently, 25 core targets were found differentially expressed in brain regions by bioinformatics analyses of GEO datasets of clinical samples from the Alzdata database. The binding sites and stabilities between the active components and the core targets were investigated by the molecular docking approach using Autodock 4.2.6 software, followed by pocket detection and druggability assessment via the DoGSiteScorer server. The results showed that acacetin, β-sitosterol, and 3-O-acetyldammarenediol-II strongly interacted with the druggable pockets of AR, CASP8, POLB, and PREP. Eventually, the docking results were further cross-referenced with the literature research and validated by 100 ns of molecular dynamics simulations using GROMACS software. Binding free energies were calculated via MM/PBSA strategy combined with interaction entropy. The simulation results indicated stable bindings between four docking pairs including acacetin-AR, acacetin-CASP8, β-sitosterol-POLB, and 3-O-acetyldammarenediol-II-PREP. Overall, our study demonstrated a theoretical basis for how three active components of EZW confer efficacy against AD. It provides a promising reference for subsequent research regarding drug discoveries and clinical applications.
Collapse
Affiliation(s)
- Meng Yu
- Innovative Institute of Chinese Medicine and Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, Shandong, 250355, China
| | - Zhongqi Shen
- Institute of Chinese Medical Literature and Culture, Shandong University of Traditional Chinese Medicine, Jinan, Shandong, 250355, China
| | - Shaozhi Zhang
- Innovative Institute of Chinese Medicine and Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, Shandong, 250355, China
| | - Yang Zhang
- Innovative Institute of Chinese Medicine and Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, Shandong, 250355, China
| | - Hongwei Zhao
- Innovative Institute of Chinese Medicine and Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, Shandong, 250355, China
| | - Longfei Zhang
- Experimental Center, Shandong University of Traditional Chinese Medicine, Jinan, Shandong, 250355, China
| |
Collapse
|
2
|
Maffeo C, Chou HY, Aksimentiev A. Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations. ADVANCED THEORY AND SIMULATIONS 2019; 2:1800191. [PMID: 31728433 PMCID: PMC6855400 DOI: 10.1002/adts.201800191] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2018] [Indexed: 12/15/2022]
Abstract
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.
Collapse
Affiliation(s)
- Christopher Maffeo
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| | - Han-Yi Chou
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| | - Aleksei Aksimentiev
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| |
Collapse
|
3
|
Affiliation(s)
- Vito Genna
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| | - Elisa Donati
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| | - Marco De Vivo
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| |
Collapse
|
4
|
Raper AT, Reed AJ, Suo Z. Kinetic Mechanism of DNA Polymerases: Contributions of Conformational Dynamics and a Third Divalent Metal Ion. Chem Rev 2018; 118:6000-6025. [DOI: 10.1021/acs.chemrev.7b00685] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Austin T. Raper
- Department of Chemistry and Biochemistry, Ohio State Biochemistry Program, The Ohio State University, Columbus, Ohio 43210, United States
| | - Andrew J. Reed
- Department of Chemistry and Biochemistry, Ohio State Biochemistry Program, The Ohio State University, Columbus, Ohio 43210, United States
| | - Zucai Suo
- Department of Chemistry and Biochemistry, Ohio State Biochemistry Program, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
5
|
Vyas R, Reed AJ, Raper AT, Zahurancik WJ, Wallenmeyer PC, Suo Z. Structural basis for the D-stereoselectivity of human DNA polymerase β. Nucleic Acids Res 2017; 45:6228-6237. [PMID: 28402499 PMCID: PMC5449621 DOI: 10.1093/nar/gkx252] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2016] [Revised: 03/24/2017] [Accepted: 04/03/2017] [Indexed: 12/20/2022] Open
Abstract
Nucleoside reverse transcriptase inhibitors (NRTIs) with L-stereochemistry have long been an effective treatment for viral infections because of the strong D-stereoselectivity exhibited by human DNA polymerases relative to viral reverse transcriptases. The D-stereoselectivity of DNA polymerases has only recently been explored structurally and all three DNA polymerases studied to date have demonstrated unique stereochemical selection mechanisms. Here, we have solved structures of human DNA polymerase β (hPolβ), in complex with single-nucleotide gapped DNA and L-nucleotides and performed pre-steady-state kinetic analysis to determine the D-stereoselectivity mechanism of hPolβ. Beyond a similar 180° rotation of the L-nucleotide ribose ring seen in other studies, the pre-catalytic ternary crystal structures of hPolβ, DNA and L-dCTP or the triphosphate forms of antiviral drugs lamivudine ((-)3TC-TP) and emtricitabine ((-)FTC-TP) provide little structural evidence to suggest that hPolβ follows the previously characterized mechanisms of D-stereoselectivity. Instead, hPolβ discriminates against L-stereochemistry through accumulation of several active site rearrangements that lead to a decreased nucleotide binding affinity and incorporation rate. The two NRTIs escape some of the active site selection through the base and sugar modifications but are selected against through the inability of hPolβ to complete thumb domain closure.
Collapse
Affiliation(s)
- Rajan Vyas
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH 43210, USA
| | - Andrew J. Reed
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH 43210, USA
- The Ohio State Biochemistry Program, The Ohio State University, Columbus, OH 43210, USA
| | - Austin T. Raper
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH 43210, USA
- The Ohio State Biochemistry Program, The Ohio State University, Columbus, OH 43210, USA
| | - Walter J. Zahurancik
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH 43210, USA
- The Ohio State Biochemistry Program, The Ohio State University, Columbus, OH 43210, USA
| | - Petra C. Wallenmeyer
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH 43210, USA
| | - Zucai Suo
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH 43210, USA
- The Ohio State Biochemistry Program, The Ohio State University, Columbus, OH 43210, USA
| |
Collapse
|
6
|
Warshel A, Bora RP. Perspective: Defining and quantifying the role of dynamics in enzyme catalysis. J Chem Phys 2017; 144:180901. [PMID: 27179464 DOI: 10.1063/1.4947037] [Citation(s) in RCA: 134] [Impact Index Per Article: 19.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Enzymes control chemical reactions that are key to life processes, and allow them to take place on the time scale needed for synchronization between the relevant reaction cycles. In addition to general interest in their biological roles, these proteins present a fundamental scientific puzzle, since the origin of their tremendous catalytic power is still unclear. While many different hypotheses have been put forward to rationalize this, one of the proposals that has become particularly popular in recent years is the idea that dynamical effects contribute to catalysis. Here, we present a critical review of the dynamical idea, considering all reasonable definitions of what does and does not qualify as a dynamical effect. We demonstrate that no dynamical effect (according to these definitions) has ever been experimentally shown to contribute to catalysis. Furthermore, the existence of non-negligible dynamical contributions to catalysis is not supported by consistent theoretical studies. Our review is aimed, in part, at readers with a background in chemical physics and biophysics, and illustrates that despite a substantial body of experimental effort, there has not yet been any study that consistently established a connection between an enzyme's conformational dynamics and a significant increase in the catalytic contribution of the chemical step. We also make the point that the dynamical proposal is not a semantic issue but a well-defined scientific hypothesis with well-defined conclusions.
Collapse
Affiliation(s)
- Arieh Warshel
- Department of Chemistry, University of Southern California, SGM 418, 3620 McClintock Avenue, Los Angeles, California 90089, USA
| | - Ram Prasad Bora
- Department of Chemistry, University of Southern California, SGM 418, 3620 McClintock Avenue, Los Angeles, California 90089, USA
| |
Collapse
|
7
|
Lehmann J. Induced fit of the peptidyl-transferase center of the ribosome and conformational freedom of the esterified amino acids. RNA (NEW YORK, N.Y.) 2017; 23:229-239. [PMID: 27879432 PMCID: PMC5238797 DOI: 10.1261/rna.057273.116] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/01/2015] [Accepted: 11/18/2016] [Indexed: 06/06/2023]
Abstract
The catalytic site of most enzymes can efficiently handle only one substrate. In contrast, the ribosome is capable of polymerizing at a similar rate at least 20 different kinds of amino acids from aminoacyl-tRNA carriers while using just one catalytic site, the peptidyl-transferase center (PTC). An induced-fit mechanism has been uncovered in the PTC, but a possible connection between this mechanism and the uniform handling of the substrates has not been investigated. We present an analysis of published ribosome structures supporting the hypothesis that the induced fit eliminates unreactive rotamers predominantly populated for some A-site aminoacyl esters before induction. We show that this hypothesis is fully consistent with the wealth of kinetic data obtained with these substrates. Our analysis reveals that induction constrains the amino acids into a reactive conformation in a side-chain independent manner. It allows us to highlight the rationale of the PTC structural organization, which confers to the ribosome the very unusual ability to handle large as well as small substrates.
Collapse
Affiliation(s)
- Jean Lehmann
- Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Université Paris-Sud, Campus Paris-Saclay, 91198 Gif-sur-Yvette, France
| |
Collapse
|
8
|
Klvaňa M, Bren U, Florián J. Uniform Free-Energy Profiles of the P-O Bond Formation and Cleavage Reactions Catalyzed by DNA Polymerases β and λ. J Phys Chem B 2016; 120:13017-13030. [PMID: 27992186 PMCID: PMC5217713 DOI: 10.1021/acs.jpcb.6b08581] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
![]()
Human
X-family DNA polymerases β (Polβ) and λ
(Polλ) catalyze the nucleotidyl-transfer reaction in the base
excision repair pathway of the cellular DNA damage response. Using
empirical valence bond and free-energy perturbation simulations, we
explore the feasibility of various mechanisms for the deprotonation
of the 3′-OH group of the primer DNA strand, and the subsequent
formation and cleavage of P–O bonds in four Polβ, two
truncated Polλ (tPolλ), and two tPolλ Loop1 mutant
(tPolλΔL1) systems differing in the initial X-ray crystal
structure and nascent base pair. The average calculated activation
free energies of 14, 18, and 22 kcal mol–1 for Polβ,
tPolλ, and tPolλΔL1, respectively, reproduce the
trend in the observed catalytic rate constants. The most feasible
reaction pathway consists of two successive steps: specific base (SB)
proton transfer followed by rate-limiting concerted formation and
cleavage of the P–O bonds. We identify linear free-energy relationships
(LFERs) which show that the differences in the overall activation
and reaction free energies among the eight studied systems are determined
by the reaction free energy of the SB proton transfer. We discuss
the implications of the LFERs and suggest pKa of the 3′-OH group as a predictor of the catalytic
rate of X-family DNA polymerases.
Collapse
Affiliation(s)
- Martin Klvaňa
- Laboratory of Physical Chemistry and Chemical Thermodynamics, Faculty of Chemistry and Chemical Technology, University of Maribor , Smetanova ulica 17, 2000 Maribor, Slovenia.,Department of Chemistry and Biochemistry, Loyola University Chicago , 1032 W. Sheridan Road, Chicago, Illinois 60660, United States
| | - Urban Bren
- Laboratory of Physical Chemistry and Chemical Thermodynamics, Faculty of Chemistry and Chemical Technology, University of Maribor , Smetanova ulica 17, 2000 Maribor, Slovenia.,Laboratory for Molecular Modeling, National Institute of Chemistry , Hajdrihova ulica 19, 1001 Ljubljana, Slovenia
| | - Jan Florián
- Department of Chemistry and Biochemistry, Loyola University Chicago , 1032 W. Sheridan Road, Chicago, Illinois 60660, United States
| |
Collapse
|
9
|
Genna V, Vidossich P, Ippoliti E, Carloni P, De Vivo M. A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases. J Am Chem Soc 2016; 138:14592-14598. [PMID: 27530537 DOI: 10.1021/jacs.6b05475] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
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.
Collapse
Affiliation(s)
- Vito Genna
- Laboratory of Molecular Modeling & Drug Discovery, Istituto Italiano di Tecnologia , Via Morego 30, 16163, Genoa, Italy.,IAS-5/INM-9 Computational Biomedicine and JARA-HPC, Forschungszentrum Jülich , Wilhelm-Johnen-Strasse, 52428 Jülich, Germany
| | - Pietro Vidossich
- IAS-5/INM-9 Computational Biomedicine and JARA-HPC, Forschungszentrum Jülich , Wilhelm-Johnen-Strasse, 52428 Jülich, Germany
| | - Emiliano Ippoliti
- IAS-5/INM-9 Computational Biomedicine and JARA-HPC, Forschungszentrum Jülich , Wilhelm-Johnen-Strasse, 52428 Jülich, Germany
| | - Paolo Carloni
- IAS-5/INM-9 Computational Biomedicine and JARA-HPC, Forschungszentrum Jülich , Wilhelm-Johnen-Strasse, 52428 Jülich, Germany
| | - Marco De Vivo
- Laboratory of Molecular Modeling & Drug Discovery, Istituto Italiano di Tecnologia , Via Morego 30, 16163, Genoa, Italy.,IAS-5/INM-9 Computational Biomedicine and JARA-HPC, Forschungszentrum Jülich , Wilhelm-Johnen-Strasse, 52428 Jülich, Germany
| |
Collapse
|
10
|
Gaspard P. Kinetics and thermodynamics of exonuclease-deficient DNA polymerases. Phys Rev E 2016; 93:042419. [PMID: 27176340 DOI: 10.1103/physreve.93.042419] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Indexed: 05/02/2023]
Abstract
A kinetic theory is developed for exonuclease-deficient DNA polymerases, based on the experimental observation that the rates depend not only on the newly incorporated nucleotide, but also on the previous one, leading to the growth of Markovian DNA sequences from a Bernoullian template. The dependencies on nucleotide concentrations and template sequence are explicitly taken into account. In this framework, the kinetic and thermodynamic properties of DNA replication, in particular, the mean growth velocity, the error probability, and the entropy production are calculated analytically in terms of the rate constants and the concentrations. Theory is compared with numerical simulations for the DNA polymerases of T7 viruses and human mitochondria.
Collapse
Affiliation(s)
- Pierre Gaspard
- Center for Nonlinear Phenomena and Complex Systems, Université Libre de Bruxelles, Code Postal 231, Campus Plaine, B-1050 Brussels, Belgium
| |
Collapse
|
11
|
Gaspard P. Kinetics and thermodynamics of DNA polymerases with exonuclease proofreading. Phys Rev E 2016; 93:042420. [PMID: 27176341 DOI: 10.1103/physreve.93.042420] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Indexed: 05/02/2023]
Abstract
Kinetic theory and thermodynamics are applied to DNA polymerases with exonuclease activity, taking into account the dependence of the rates on the previously incorporated nucleotide. The replication fidelity is shown to increase significantly thanks to this dependence at the basis of the mechanism of exonuclease proofreading. In particular, this dependence can provide up to a 100-fold lowering of the error probability under physiological conditions. Theory is compared with numerical simulations for the DNA polymerases of T7 viruses and human mitochondria.
Collapse
Affiliation(s)
- Pierre Gaspard
- Center for Nonlinear Phenomena and Complex Systems, Université Libre de Bruxelles, Code Postal 231, Campus Plaine, B-1050 Brussels, Belgium
| |
Collapse
|
12
|
Moscato B, Swain M, Loria JP. Induced Fit in the Selection of Correct versus Incorrect Nucleotides by DNA Polymerase β. Biochemistry 2015; 55:382-95. [PMID: 26678253 DOI: 10.1021/acs.biochem.5b01213] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
DNA polymerase β (Pol β) repairs single-nucleotide gapped DNA (sngDNA) by enzymatic incorporation of the Watson-Crick partner nucleotide at the gapped position opposite the templating nucleotide. The process by which the matching nucleotide is incorporated into a sngDNA sequence has been relatively well-characterized, but the process of discrimination from nucleotide misincorporation remains unclear. We report here NMR spectroscopic characterization of full-length, uniformly labeled Pol β in apo, sngDNA-bound binary, and ternary complexes containing matching and mismatching nucleotide. Our data indicate that, while binding of the correct nucleotide to the binary complex induces chemical shift changes consistent with the process of enzyme closure, the ternary Pol β complex containing a mismatching nucleotide exhibits no such changes and appears to remain in an open, unstable, binary-like conformation. Our findings support an induced-fit mechanism for polymerases in which a closed ternary complex can only be achieved in the presence of matching nucleotide.
Collapse
Affiliation(s)
- Beth Moscato
- Department of Chemistry, Yale University , 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Monalisa Swain
- Department of Chemistry, Yale University , 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - J Patrick Loria
- Department of Chemistry, Yale University , 225 Prospect Street, New Haven, Connecticut 06520, United States.,Department of Molecular Biophysics and Biochemistry, Yale University , 260 Whitney Avenue, New Haven, Connecticut 06520, United States
| |
Collapse
|
13
|
Bora RP, Mills MJL, Frushicheva MP, Warshel A. On the Challenge of Exploring the Evolutionary Trajectory from Phosphotriesterase to Arylesterase Using Computer Simulations. J Phys Chem B 2015; 119:3434-45. [DOI: 10.1021/jp5124025] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Ram Prasad Bora
- Department
of Chemistry, University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089, United States
| | - Matthew J. L. Mills
- Deconstruction Division, Joint BioEnergy Institute, Emeryville, California, United States & Biomass Science and Conversion Technology Department, Sandia National Laboratories, Livermore, California 94550, United States
| | - Maria P. Frushicheva
- Department
of Chemical Engineering, Massachusetts Institute of Technology, 25 Ames
Street, Cambridge, Massachusetts 02142, United States
| | - Arieh Warshel
- Department
of Chemistry, University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089, United States
| |
Collapse
|
14
|
Elias AA, Cisneros GA. Computational study of putative residues involved in DNA synthesis fidelity checking in Thermus aquaticus DNA polymerase I. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2014; 96:39-75. [PMID: 25443954 DOI: 10.1016/bs.apcsb.2014.06.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
A fidelity-checking site for DNA polymerase I has been proposed based on recent single-molecule Förster resonance energy transfer studies. The checking site is believed to ensure proper base pairing of the newly inserted nucleotide. Computational studies have been utilized to predict residues involved in this putative checking site on the Klenow and Bacillus fragments. Here, we employ energy decomposition analysis, electrostatic free energy response, and noncovalent interaction plots to identify the residues involved in the hypothesized checking site in the homologous Klenow fragment from Thermus aquaticus (Klentaq). Our results indicate multiple protein residues that show altered interactions for three mispairs compared to the correctly paired DNA dimer. Many of these residues are also conserved along A family polymerases.
Collapse
Affiliation(s)
- Angela A Elias
- Department of Chemistry, Wayne State University, Detroit, Michigan, USA
| | - G Andrés Cisneros
- Department of Chemistry, Wayne State University, Detroit, Michigan, USA.
| |
Collapse
|
15
|
Abstract
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.
Collapse
|
16
|
Ram Prasad B, Kamerlin SCL, Florián J, Warshel A. Prechemistry barriers and checkpoints do not contribute to fidelity and catalysis as long as they are not rate limiting. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1288-6] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
17
|
Klvaňa M, Murphy DL, Jeřábek P, Goodman MF, Warshel A, Sweasy JB, Florián J. Catalytic effects of mutations of distant protein residues in human DNA polymerase β: theory and experiment. Biochemistry 2012; 51:8829-43. [PMID: 23013478 DOI: 10.1021/bi300783t] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We carried out free-energy calculations and transient kinetic experiments for the insertion of the right (dC) and wrong (dA) nucleotides by wild-type (WT) and six mutant variants of human DNA polymerase β (Pol β). Since the mutated residues in the point mutants, I174S, I260Q, M282L, H285D, E288K, and K289M, were not located in the Pol β catalytic site, we assumed that the WT and its point mutants share the same dianionic phosphorane transition-state structure of the triphosphate moiety of deoxyribonucleotide 5'-triphosphate (dNTP) substrate. On the basis of this assumption, we have formulated a thermodynamic cycle for calculating relative dNTP insertion efficiencies, Ω = (k(pol)/K(D))(mut)/(k(pol)/K(D))(WT) using free-energy perturbation (FEP) and linear interaction energy (LIE) methods. Kinetic studies on five of the mutants have been published previously using different experimental conditions, e.g., primer-template sequences. We have performed a presteady kinetic analysis for the six mutants for comparison with wild-type Pol β using the same conditions, including the same primer/template DNA sequence proximal to the dNTP insertion site used for X-ray crystallographic studies. This consistent set of kinetic and structural data allowed us to eliminate the DNA sequence from the list of factors that can adversely affect calculated Ω values. The calculations using the FEP free energies scaled by 0.5 yielded 0.9 and 1.1 standard deviations from the experimental log Ω values for the insertion of the right and wrong dNTP, respectively. We examined a hybrid FEP/LIE method in which the FEP van der Waals term for the interaction of the mutated amino acid residue with its surrounding environment was replaced by the corresponding van der Waals term calculated using the LIE method, resulting in improved 0.4 and 1.0 standard deviations from the experimental log Ω values. These scaled FEP and FEP/LIE methods were also used to predict log Ω for R283A and R283L Pol β mutants.
Collapse
Affiliation(s)
- Martin Klvaňa
- Department of Chemistry, Loyola University, Chicago, Illinois 60626, United States
| | | | | | | | | | | | | |
Collapse
|
18
|
Graham SE, Syeda F, Cisneros GA. Computational prediction of residues involved in fidelity checking for DNA synthesis in DNA polymerase I. Biochemistry 2012; 51:2569-78. [PMID: 22397306 DOI: 10.1021/bi201856m] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Recent single-molecule Förster resonance energy transfer studies of DNA polymerase I have led to the proposal of a postinsertion fidelity-checking site. This site is hypothesized to ensure proper base pairing of the newly inserted nucleotide. To help test this hypothesis, we have used energy decomposition, electrostatic free energy response, and noncovalent interaction analysis analyses to identify residues involved in this putative checking site. We have used structures of DNA polymerase I from two different organisms, the Klenow fragment from Escherichia coli and the Bacillus fragment from Bacillus stearothermophilus. Our results point to several residues that show altered interactions for three mispairs compared to the correctly paired DNA dimer. Furthermore, many of these residues are conserved among A family polymerases. The identified residues provide potential targets for mutagenesis studies for investigation of the fidelity-checking site hypothesis.
Collapse
Affiliation(s)
- Sarah E Graham
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | | | | |
Collapse
|
19
|
Ram Prasad B, Warshel A. Prechemistry versus preorganization in DNA replication fidelity. Proteins 2011; 79:2900-19. [PMID: 21905114 DOI: 10.1002/prot.23128] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2011] [Revised: 06/30/2011] [Accepted: 07/05/2011] [Indexed: 01/30/2023]
Abstract
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.
Collapse
Affiliation(s)
- B Ram Prasad
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | |
Collapse
|
20
|
Klvaňa M, Jeřábek P, Goodman MF, Florián J. An abridged transition state model to derive structure, dynamics, and energy components of DNA polymerase β fidelity. Biochemistry 2011; 50:7023-32. [PMID: 21739967 DOI: 10.1021/bi200790s] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We show how a restricted reaction surface can be used to facilitate the calculation of biologically important contributions of active site geometries and dynamics to DNA polymerase fidelity. Our analysis, using human DNA polymerase beta (pol β), is performed within the framework of an electrostatic linear free energy response (EFER) model. The structure, dynamics, and energetics of pol β-DNA-dNTP interactions are computed between two points on the multidimensional reaction free energy surface. "Point 1" represents a ground state activation intermediate (GSA), which is obtained by deprotonating the terminal 3'OH group of the primer DNA strand. "Point 2" is the transition state (PTS) for the attack of the 3'O(-) (O(nuc)) on the P(α) atom of dNTP substrate, having the electron density of a dianionic phosphorane intermediate. Classical molecular dynamics simulations are used to compute the geometric and dynamic contributions to the formation of right and wrong O(nuc)-P chemical bonds. Matched dCTP·G and mismatched dATP·G base pairs are used to illustrate the analysis. Compared to the dCTP·G base pair, the dATP·G mismatch has fewer GSA configurations with short distances between O(nuc) and P(α) atoms and between the oxygen in the scissile P-O bond (O(lg)) and the nearest structural water. The thumb subdomain conformation of the GSA complex is more open for the mismatch, and the H-bonds in the mispair become more extended during the nucleophilic attack than in the correct pair. The electrostatic contributions of pol β and DNA residues to catalysis of the right and wrong P-O(nuc) bond formation are 5.3 and 3.1 kcal/mol, respectively, resulting in an 80-fold contribution to fidelity. The EFER calculations illustrate the considerable importance of Arg183 and an O(lg)-proximal water molecule to pol β fidelity.
Collapse
Affiliation(s)
- Martin Klvaňa
- Department of Chemistry, Loyola University Chicago, Chicago, Illinois 60626, United States
| | | | | | | |
Collapse
|
21
|
Abstract
Recent years have witnessed a tremendous explosion in computational power, which in turn has resulted in great progress in the complexity of the biological and chemical problems that can be addressed by means of all-atom simulations. Despite this, however, our computational time is not infinite, and in fact many of the key problems of the field were resolved long before the existence of the current levels of computational power. This review will start by presenting a brief historical overview of the use of multiscale simulations in biology, and then present some key developments in the field, highlighting several cases where the use of a physically sound simplification is clearly superior to a brute-force approach. Finally, some potential future directions will be discussed.
Collapse
|
22
|
Dryga A, Warshel A. Renormalizing SMD: the renormalization approach and its use in long time simulations and accelerated PMF calculations of macromolecules. J Phys Chem B 2010; 114:12720-8. [PMID: 20836533 PMCID: PMC2948080 DOI: 10.1021/jp1056122] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Simulations of long time process in condensed phases, in general, and in biomolecules, in particular, present a major challenge that cannot be overcome at present by brute force molecular dynamics (MD) approaches. This work takes the renormalization method, intruded by us sometime ago, and establishes its reliability and potential in extending the time scale of molecular simulations. The validation involves a truncated gramicidin system in the gas phase. This system is small enough to allow for very long explicit simulations and sufficiently complex to present the physics of realistic ion channels. The renormalization approach is found to be reliable and arguably presents the first approach that allows one to exploit the otherwise problematic steered molecular dynamics (SMD) treatments in quantitative and meaningful studies. It is established that we can reproduce the long time behavior of large systems by using Langevin dynamics (LD) simulations of a renormalized implicit model. This is done without spending the enormous time needed to obtain such trajectories in the explicit system. The present study also provides a promising advance in accelerated evaluation of free energy barriers. This is done by adjusting the effective potential in the implicit model to reproduce the same passage time as that obtained in the explicit model under the influence of an external force. Here having a reasonable effective friction provides a way to extract the potential of mean force (PMF) without investing the time needed for regular PMF calculations. The renormalization approach, which is illustrated here in realistic calculations, is expected to provide a major help in studies of complex landscapes and in exploring long time dynamics of biomolecules.
Collapse
Affiliation(s)
- Anatoly Dryga
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062
| | - Arieh Warshel
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062
| |
Collapse
|
23
|
Nashine VC, Hammes-Schiffer S, Benkovic SJ. Coupled motions in enzyme catalysis. Curr Opin Chem Biol 2010; 14:644-51. [PMID: 20729130 DOI: 10.1016/j.cbpa.2010.07.020] [Citation(s) in RCA: 137] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2010] [Revised: 07/14/2010] [Accepted: 07/24/2010] [Indexed: 10/19/2022]
Abstract
Recently, the hypothesis that protein motions are involved in enzymatic turnover has gained significant attention. We review cases where there is evidence that protein motions are rate-limiting in the overall catalytic cycle and examine experimental and theoretical evidence for how such motions enhance the probability of sampling the transition state configurations relative to the ground state. The impact of tunneling, the possible role of vibrational coupling and the value of conformational chemical landscapes are also scrutinized.
Collapse
Affiliation(s)
- Vishal C Nashine
- Department of Chemistry, The Pennsylvania State University, 414 Wartik Laboratory, University Park, PA 16802, USA
| | | | | |
Collapse
|
24
|
Bhattacharyya M, Vishveshwara S. Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis. BMC STRUCTURAL BIOLOGY 2010; 10:27. [PMID: 20704697 PMCID: PMC2929236 DOI: 10.1186/1472-6807-10-27] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/25/2010] [Accepted: 08/12/2010] [Indexed: 11/11/2022]
Abstract
Background One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results In this study, we have explored the ligand induced conformational changes in H.pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H.pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co-operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.
Collapse
|
25
|
Kamerlin SCL, Warshel A. At the dawn of the 21st century: Is dynamics the missing link for understanding enzyme catalysis? Proteins 2010; 78:1339-75. [PMID: 20099310 PMCID: PMC2841229 DOI: 10.1002/prot.22654] [Citation(s) in RCA: 345] [Impact Index Per Article: 24.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
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.
Collapse
Affiliation(s)
- Shina C. L. Kamerlin
- Department of Chemistry, University of Southern California, 3620 McClintock Ave., Los Angeles CA-90089, USA
| | - Arieh Warshel
- Department of Chemistry, University of Southern California, 3620 McClintock Ave., Los Angeles CA-90089, USA
| |
Collapse
|
26
|
Kamerlin SCL, Mavri J, Warshel A. Examining the case for the effect of barrier compression on tunneling, vibrationally enhanced catalysis, catalytic entropy and related issues. FEBS Lett 2010; 584:2759-66. [PMID: 20433839 DOI: 10.1016/j.febslet.2010.04.062] [Citation(s) in RCA: 66] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2010] [Revised: 04/21/2010] [Accepted: 04/22/2010] [Indexed: 10/19/2022]
Abstract
The idea that tunneling is enhanced by the compression of the donor-acceptor distance has attracted significant interest. In particular, recent studies argued that this proposal is consistent with pressure effects on enzymatic reactions, and that the observed pressure effects support the idea of vibrationally enhanced catalysis. However, a careful analysis of the current works reveals serious inconsistencies in the evidence presented to support these hypotheses. Apparently, tunneling decreases upon compression, and external pressure does not lead to the applicable compression of the free energy surface. Additionally, pressure experiments do not provide actual evidence for vibrationally enhanced catalysis. Finally, the temperature dependence of the entropy change in hydride transfer reactions is shown to reflect simple electrostatic effects.
Collapse
|
27
|
Messer BM, Roca M, Chu ZT, Vicatos S, Kilshtain AV, Warshel A. Multiscale simulations of protein landscapes: using coarse-grained models as reference potentials to full explicit models. Proteins 2010; 78:1212-27. [PMID: 20052756 DOI: 10.1002/prot.22640] [Citation(s) in RCA: 54] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Evaluating the free-energy landscape of proteins and the corresponding functional aspects presents a major challenge for computer simulation approaches. This challenge is due to the complexity of the landscape and the enormous computer time needed for converging simulations. The use of simplified coarse-grained (CG) folding models offers an effective way of sampling the landscape but such a treatment, however, may not give the correct description of the effect of the actual protein residues. A general way around this problem that has been put forward in our early work (Fan et al., Theor Chem Acc 1999;103:77-80) uses the CG model as a reference potential for free-energy calculations of different properties of the explicit model. This method is refined and extended here, focusing on improving the electrostatic treatment and on demonstrating key applications. These applications include: evaluation of changes of folding energy upon mutations, calculations of transition-states binding free energies (which are crucial for rational enzyme design), evaluations of catalytic landscape, and evaluations of the time-dependent responses to pH changes. Furthermore, the general potential of our approach in overcoming major challenges in studies of structure function correlation in proteins is discussed.
Collapse
Affiliation(s)
- Benjamin M Messer
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062, USA
| | | | | | | | | | | |
Collapse
|
28
|
Rucker R, Oelschlaeger P, Warshel A. A binding free energy decomposition approach for accurate calculations of the fidelity of DNA polymerases. Proteins 2010; 78:671-80. [PMID: 19842163 DOI: 10.1002/prot.22596] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
DNA polymerase beta (pol beta) is a small eukaryotic enzyme with the ability to repair short single-stranded DNA gaps that has found use as a model system for larger replicative DNA polymerases. For all DNA polymerases, the factors determining their catalytic power and fidelity are the interactions between the bases of the base pair, amino acids near the active site, and the two magnesium ions. In this report, we study effects of all three aspects on human pol beta transition state (TS) binding free energies by reproducing a consistent set of experimentally determined data for different structures. Our calculations comprise the combination of four different base pairs (incoming pyrimidine nucleotides incorporated opposite both matched and mismatched purines) with four different pol beta structures (wild type and three mutants). We generate three fragments of the incoming deoxynucleoside 5'-triphosphate-TS and run separate calculations for the neutral base part and the highly charged triphosphate part, using different dielectric constants in order to account for the specific dielectric response. This new approach improves our ability to predict the effect of matched and mismatched base pairing and of mutations in DNA polymerases on fidelity and may be a useful tool in studying the potential of DNA polymerase mutations in the development of cancer. It also supports our point of view with regards to the origin of the structural control of fidelity, allowing for a quantified description of the fidelity of DNA polymerases.
Collapse
Affiliation(s)
- Robert Rucker
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | | |
Collapse
|
29
|
Lee HR, Wang M, Konigsberg W. The reopening rate of the fingers domain is a determinant of base selectivity for RB69 DNA polymerase. Biochemistry 2009; 48:2087-98. [PMID: 19228036 DOI: 10.1021/bi8016284] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Two divalent metal ions are required for nucleotide incorporation by DNA polymerases. Here we use the bacteriophage RB69 DNA polymerase (RB69 pol) and the metal ion exchange-inert nucleotide analogue rhodium(III) deoxythymidine triphosphate (Rh.dTTP) to investigate the requirements of metal binding to the "A" site and to the "B" site, independently. We show that while binding of a metal ion to the A site is required for the nucleotidyl transfer reaction to occur, this metal binding is insufficient to initiate the prechemistry enzyme isomerization that has been observed with this polymerase. Moreover, we show that binding of a deoxynucleoside triphosphate (dNTP), in the absence of a catalytic metal ion, is sufficient to induce this conformational change. In this report, we also present several lines of evidence (from pulse-chase, rapid chemical quench-flow, and stopped-flow fluorescence experiments) for the reverse rate of the enzyme isomerization, closed to open, of a DNA polymerase complex. The implications of these data for the fidelity of DNA polymerization by RB69 pol are discussed.
Collapse
Affiliation(s)
- Harold R Lee
- Molecular Biophysics and Biochemistry Department, Yale University, 333 Cedar Street, New Haven, Connecticut 06520, USA
| | | | | |
Collapse
|
30
|
Effect of N2-guanyl modifications on early steps in catalysis of polymerization by Sulfolobus solfataricus P2 DNA polymerase Dpo4 T239W. J Mol Biol 2009; 395:1007-18. [PMID: 19969000 DOI: 10.1016/j.jmb.2009.11.071] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2009] [Revised: 11/06/2009] [Accepted: 11/30/2009] [Indexed: 10/20/2022]
Abstract
Translesion DNA polymerases are more efficient at bypass of many DNA adducts than replicative polymerases. Previous work with the translesion polymerase Sulfolobus solfataricus Dpo4 showed a decrease in catalytic efficiency during bypass of bulky N(2)-alkyl guanine (G) adducts with N(2)-isobutylG showing the largest effect, decreasing approximately 120-fold relative to unmodified deoxyguanosine (Zhang, H., Eoff, R. L., Egli, M., Guengerich, F. P. Versatility of Y-family Sulfolobus solfataricus DNA polymerase Dpo4 in translation synthesis past bulky N(2)-alkylguanine adducts. J. Biol. Chem. 2009; 284: 3563-3576). The effect of adduct size on individual catalytic steps has not been easy to decipher because of the difficulty of distinguishing early noncovalent steps from phosphodiester bond formation. We developed a mutant with a single Trp (T239W) to monitor fluorescence changes associated with a conformational change that occurs after binding a correct 2'-deoxyribonucleoside triphosphate (Beckman, J. W., Wang, Q., Guengerich, F. P. Kinetic analysis of nucleotide insertion by a Y-family DNA polymerase reveals conformational change both prior to and following phosphodiester bond formation as detected by tryptophan fluorescence. J. Biol. Chem. 2008; 283: 36711-36723) and, in the present work, utilized this approach to monitor insertion opposite N(2)-alkylG-modified oligonucleotides. We estimated maximal rates for the forward conformational step, which coupled with measured rates of product formation yielded rate constants for the conformational step (both directions) during insertion opposite several N(2)-alkylG adducts. With the smaller N(2)-alkylG adducts, the conformational rate constants were not changed dramatically (<3-fold), indicating that the more sensitive steps are phosphodiester bond formation and partitioning into inactive complexes. With the larger adducts (>or=(2-naphthyl)methyl), the absence of fluorescence changes suggests impaired ability to undergo an appropriate conformational change, consistent with previous structural work.
Collapse
|
31
|
Beard WA, Shock DD, Batra VK, Pedersen LC, Wilson SH. DNA polymerase beta substrate specificity: side chain modulation of the "A-rule". J Biol Chem 2009; 284:31680-9. [PMID: 19759017 DOI: 10.1074/jbc.m109.029843] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Apurinic/apyrimidinic (AP) sites are continuously generated in genomic DNA. Left unrepaired, AP sites represent noninstructional premutagenic lesions that are impediments to DNA synthesis. When DNA polymerases encounter an AP site, they generally insert dAMP. This preferential insertion is referred to as the A-rule. Crystallographic structures of DNA polymerase (pol) beta, a family X polymerase, with active site mismatched nascent base pairs indicate that the templating (i.e. coding) base is repositioned outside of the template binding pocket thereby diminishing interactions with the incorrect incoming nucleotide. This effectively produces an abasic site because the template pocket is devoid of an instructional base. However, the template pocket is not empty; an arginine residue (Arg-283) occupies the space vacated by the templating nucleotide. In this study, we analyze the kinetics of pol beta insertion opposite an AP site and show that the preferential incorporation of dAMP is lost with the R283A mutant. The crystallographic structures of pol beta bound to gapped DNA with an AP site analog (tertrahydrofuran) in the gap (binary complex) and with an incoming nonhydrolyzable dATP analog (ternary complex) were solved. These structures reveal that binding of the dATP analog induces a closed polymerase conformation, an unstable primer terminus, and an upstream shift of the templating residue even in the absence of a template base. Thus, dATP insertion opposite an abasic site and dATP misinsertions have common features.
Collapse
Affiliation(s)
- William A Beard
- Laboratory of Structural Biology, NIEHS, National Institutes of Health, Research Triangle Park, North Carolina 27709, USA
| | | | | | | | | |
Collapse
|
32
|
Kamerlin SCL, McKenna CE, Goodman MF, Goondman MF, Warshel A. A computational study of the hydrolysis of dGTP analogues with halomethylene-modified leaving groups in solution: implications for the mechanism of DNA polymerases. Biochemistry 2009; 48:5963-71. [PMID: 19391628 DOI: 10.1021/bi900140c] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
DNA polymerases make up a family of enzymes responsible for regulating DNA replication and repair, which in turn maintains the integrity of the genome. However, despite intensive kinetic, crystallographic, and computational studies, elucidation of the detailed enzymatic mechanism still presents a significant challenge. We recently developed an alternative strategy for exploring the fidelity and mechanism of DNA polymerases, by probing leaving group effects on nucleotidyl transfer using a series of dGTP bisphosphonate analogues in which the beta,gamma-bridging oxygen was replaced by a series of substituted methylene groups (X = CYZ, where Y and Z = H, halogen, or another substituent). Pre-steady state kinetic measurements of DNA polymerase-catalyzed incorporation of correctly base paired (R) and mispaired (W) analogues demonstrated a strong linear free energy relationship (LFER) between the polymerase rate constant (k(pol)) and the highest pK(a) of the free bisphosphonic acid corresponding to the leaving group. However, unexpectedly, the data segregated into two distinctly different linear correlations depending on the nature of the substituent. The discrepancy between the two lines was considerably greater when the dGTP analogue formed an incorrect (G.T) rather than a correct (G.C) base pair, although the reason for this phenomenon remains unexplained. Here, we have evaluated the complete free energy surfaces for bisphosphonate hydrolysis in aqueous solution and evaluated the corresponding LFER. Our study, which employs several alternative solvation models, finds a split of the calculated LFER for the mono- and dihalogen compounds into two parallel lines, reflecting their behavior in the polymerase-catalyzed condensation reaction. We suggest that the division into two linear subsets may be a generalized solvation phenomenon involving the overall electrostatic interaction between the substrates and their surroundings and would also be observed in polar solvents in the absence of the enzyme, if the reaction in solvent is in fact identical to that of the enzyme. However, the amplified differences between the LFER lines for the incorporation of matched and mismatched deoxynucleotides probably reflects the differences in the electrostatic interaction between the TS charges in the polymerase active site. An understanding of the mechanism of this reaction in solution could thereby provide a steppingstone for understanding the factors governing the fidelity of DNA polymerases.
Collapse
Affiliation(s)
- Shina C L Kamerlin
- Department of Chemistry, University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089, USA.
| | | | | | | | | |
Collapse
|
33
|
Bachorz RA, Lupica G, Gutowski M, Haranczyk M. Electrostatic potential maps of damaged DNA studied by image analysis tools. 8-Oxoguanine and abasic site lesions. J Mol Model 2009; 15:817-27. [PMID: 19130099 DOI: 10.1007/s00894-008-0434-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2008] [Accepted: 12/05/2008] [Indexed: 11/28/2022]
Abstract
Changes of electrostatic potential around the DNA molecule resulting from chemical modifications of nucleotides may play a role in enzymatic recognition of damaged sites. The electrostatic potential around the DNA fragments containing either the intact guanine-cytosine pair or 8-oxoguanine-cytosine or the guanine-abasic site was projected on a cylindrical surface around the double helix. The 2D maps of EP of intact and damaged DNA fragments were compared using image analysis methods. Occurrence of abasic site and 8-oxoguanine lesions were found to be reflected in the EP maps. In the case of the 8-oxoguanine lesion, the two phosphate groups and countercations of the damaged strand are moved away from the lesion in opposite directions, whereas they are moved in the same direction in the case of the abasic site lesion. The characteristic features of 8-oxoguanine lesion might be identified in the major groove, whereas the features of abasic site lesion the minor groove.
Collapse
Affiliation(s)
- Rafał A Bachorz
- Chemical Sciences Division, Pacific Northwest National Laboratory, Richland, WA 99352, USA
| | | | | | | |
Collapse
|
34
|
Sampoli Benítez BA, Arora K, Balistreri L, Schlick T. Mismatched base-pair simulations for ASFV Pol X/DNA complexes help interpret frequent G*G misincorporation. J Mol Biol 2008; 384:1086-97. [PMID: 18955064 DOI: 10.1016/j.jmb.2008.10.025] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2008] [Revised: 10/02/2008] [Accepted: 10/06/2008] [Indexed: 10/21/2022]
Abstract
DNA polymerase X (pol X) from the African swine fever virus is a 174-amino-acid repair polymerase that likely participates in a viral base excision repair mechanism, characterized by low fidelity. Surprisingly, pol X's insertion rate of the G*G mispair is comparable to that of the four Watson-Crick base pairs. This behavior is in contrast with another X-family polymerase, DNA polymerase beta (pol beta), which inserts G*G mismatches poorly, and has higher DNA repair fidelity. Using molecular dynamics simulations, we previously provided support for an induced-fit mechanism for pol X in the presence of the correct incoming nucleotide. Here, we perform molecular dynamics simulations of pol X/DNA complexes with different incoming incorrect nucleotides in various orientations [C*C, A*G, and G*G (anti) and A*G and G*G (syn)] and compare the results to available kinetic data and prior modeling. Intriguingly, the simulations reveal that the G*G mispair with the incoming nucleotide in the syn configuration undergoes large-scale conformational changes similar to that observed in the presence of correct base pair (G*C). The base pairing in the G*G mispair is achieved via Hoogsteen hydrogen bonding with an overall geometry that is well poised for catalysis. Simulations for other mismatched base pairs show that an intermediate closed state is achieved for the A*G and G*G mispair with the incoming dGTP in anti conformation, while the protein remains near the open conformation for the C*C and the A*G syn mismatches. In addition, catalytic site geometry and base pairing at the nascent template-incoming nucleotide interaction reveal distortions and misalignments that range from moderate for A*G anti to worst for the C*C complex. These results agree well with kinetic data for pol X and provide a structural/dynamic basis to explain, at atomic level, the fidelity of this polymerase compared with other members of the X family. In particular, the more open and pliant active site of pol X, compared to pol beta, allows pol X to accommodate bulkier mismatches such as guanine opposite guanine, while the more structured and organized pol beta active site imposes higher discrimination, which results in higher fidelity. The possibility of syn conformers resonates with other low-fidelity enzymes such as Dpo4 (from the Y family), which readily accommodate oxidative lesions.
Collapse
Affiliation(s)
- Benedetta A Sampoli Benítez
- Department of Natural Sciences and Mathematics, Marymount Manhattan College, 221 East 71st Street, New York, NY 10021, USA
| | | | | | | |
Collapse
|
35
|
Wang Y, Schlick T. Quantum mechanics/molecular mechanics investigation of the chemical reaction in Dpo4 reveals water-dependent pathways and requirements for active site reorganization. J Am Chem Soc 2008; 130:13240-50. [PMID: 18785738 PMCID: PMC3195406 DOI: 10.1021/ja802215c] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- Yanli Wang
- Department of Chemistry and Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York University, New York, NY 10012
| | - Tamar Schlick
- Department of Chemistry and Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York University, New York, NY 10012
| |
Collapse
|
36
|
Roca M, Messer B, Hilvert D, Warshel A. On the relationship between folding and chemical landscapes in enzyme catalysis. Proc Natl Acad Sci U S A 2008; 105:13877-82. [PMID: 18779576 PMCID: PMC2544547 DOI: 10.1073/pnas.0803405105] [Citation(s) in RCA: 75] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2008] [Indexed: 11/18/2022] Open
Abstract
Elucidating the relationship between the folding landscape of enzymes and their catalytic power has been one of the challenges of modern enzymology. The present work explores this issue by using a simplified folding model to generate the free-energy landscape of an enzyme and then to evaluate the activation barriers for the chemical step in different regions of the landscape. This approach is used to investigate the recent finding that an engineered monomeric chorismate mutase exhibits catalytic efficiency similar to the naturally occurring dimer even though it exhibits the properties of an intrinsically disordered molten globule. It is found that the monomer becomes more confined than its native-like counterpart upon ligand binding but still retains a wider catalytic region. Although the overall rate acceleration is still determined by reduction of the reorganization energy, the detailed contribution of different barriers yields a more complex picture for the chemical process than that of a single path. This work provides insight into the relationship between folding landscapes and catalysis. The computational approach used here may also provide a powerful strategy for modeling single-molecule experiments and designing enzymes.
Collapse
Affiliation(s)
- Maite Roca
- *Department of Chemistry, University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, CA, 90089-1062; and
| | - Benjamin Messer
- *Department of Chemistry, University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, CA, 90089-1062; and
| | - Donald Hilvert
- Laboratory of Organic Chemistry, Eidgenössiche Technische Hochschule Zürich, Hönggerberg HCI F 339, CH-8093 Zürich, Switzerland
| | - Arieh Warshel
- *Department of Chemistry, University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, CA, 90089-1062; and
| |
Collapse
|
37
|
The associative nature of adenylyl transfer catalyzed by T4 DNA ligase. Proc Natl Acad Sci U S A 2008; 105:8563-8. [PMID: 18562298 DOI: 10.1073/pnas.0709140105] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
DNA ligase seals nicks in dsDNA using chemical energy of the phosphoanhydride bond in ATP or NAD(+) and assistance of a divalent metal cofactor Mg(2+). Molecular details of ligase catalysis are essential for understanding the mechanism of metal-promoted phosphoryl transfer reactions in the living cell responsible for a wide range of processes, e.g., DNA replication and transcription, signaling and differentiation, energy coupling and metabolism. Here we report a single-turnover (31)P solid-state NMR study of adenylyl transfer catalyzed by DNA ligase from bacteriophage T4. Formation of a high-energy covalent ligase-nucleotide complex is triggered in situ by the photo release of caged Mg(2+), and sequentially formed intermediates are monitored by NMR. Analyses of reaction kinetics and chemical-shift changes indicate that the pentacoordinated phosphorane intermediate builds up to 35% of the total reacting species after 4-5 h of reaction. This is direct experimental evidence of the associative nature of adenylyl transfer catalyzed by DNA ligase. NMR spectroscopy in rotating solids is introduced as an analytical tool for recording molecular movies of reaction processes. Presented work pioneers a promising direction in structural studies of biochemical transformations.
Collapse
|
38
|
Johnson KA. Role of induced fit in enzyme specificity: a molecular forward/reverse switch. J Biol Chem 2008; 283:26297-301. [PMID: 18544537 DOI: 10.1074/jbc.r800034200] [Citation(s) in RCA: 140] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Enzyme structures solved with and without bound substrate often show that substrate-induced conformational changes bring catalytic residues into alignment, alter the local environment, and position the substrate for catalysis. Although the structural data are compelling, the role of conformational changes in enzyme specificity has been controversial in that specificity is a kinetic property that is not easy to predict based upon structure alone. Recent studies on DNA polymerization have illuminated the role of substrate-induced conformational changes in enzyme specificity by showing that the rate at which the enzyme opens to release the bound substrate is a key kinetic parameter. The slow release of a correct substrate commits it to the forward reaction so that specificity is determined solely by the rate of substrate binding, including the isomerization step, and not by the slower rate of the chemical reaction. In contrast, fast dissociation of an incorrect substrate favors release rather than reaction. Thus, the conformational change acts as a molecular switch to select the right substrate and to recognize and disfavor the reaction of an incorrect substrate. A conformational switch may also favor release rather than reverse reaction of the product.
Collapse
Affiliation(s)
- Kenneth A Johnson
- Department of Chemistry and Biochemistry, Institute of Cellular and Molecular Biology, University of Texas, Austin, Texas 78712, USA.
| |
Collapse
|
39
|
Xiang Y, Warshel A. Quantifying free energy profiles of proton transfer reactions in solution and proteins by using a diabatic FDFT mapping. J Phys Chem B 2008; 112:1007-15. [PMID: 18166038 DOI: 10.1021/jp076931f] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Reliable studies of proton transfer (PT) reactions in solution and in enzymes by combined quantum mechanical/molecular mechanics (QM/MM) approaches with an ab initio description of the quantum region present a major challenge to computational chemists. The main problem is the need for extensive computer time to evaluate the QM energy, which in turn makes it extremely challenging to perform proper configurational sampling. The present work presents a new effective way for performing such calculations by using the frozen density functional (FDFT) approach to generate diabatic surfaces that are used to generate a mapping potential that takes the system from the reactant to the product state. The resulting umbrella sampling/free energy perturbation (US/FEP) mapping is done in full analogy with the approach used in the empirical valence bond (EVB) treatment, moving from the diabatic mapping potential to the adiabatic ground state surface, while an ab initio Hamiltonian is used for the QM part. The present approach provides a particularly effective way for evaluating the free energy associated with both the substrate and the solvent motions. This allows us to obtain a free energy barrier that properly reflects the solute entropy. This advance allows one to obtain ab initio QM/MM (QM(ai)/MM) free energy surfaces for very challenging cases such as the autodissociation of water in water, proton transfer between methanol and water in water, and the effect of Mg2+ ion on such a reaction. We also consider as a benchmark the initial PT reaction in the catalytic cycle of triose phosphate isomerase and obtain excellent results without any adjustable parameters. Our results point out that the present implementation of the FDFT approach provides a very promising approach for evaluating QM(ai)/MM free energy surfaces.
Collapse
Affiliation(s)
- Yun Xiang
- Department of Chemistry, University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089-1062, USA
| | | |
Collapse
|