1
|
Weissman B, Ekesan Ş, Lin HC, Gardezi S, Li NS, Giese TJ, McCarthy E, Harris ME, York DM, Piccirilli JA. Dissociative Transition State in Hepatitis Delta Virus Ribozyme Catalysis. J Am Chem Soc 2023; 145:2830-2839. [PMID: 36706353 PMCID: PMC10112047 DOI: 10.1021/jacs.2c10079] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Ribonucleases and small nucleolytic ribozymes are both able to catalyze RNA strand cleavage through 2'-O-transphosphorylation, provoking the question of whether protein and RNA enzymes facilitate mechanisms that pass through the same or distinct transition states. Here, we report the primary and secondary 18O kinetic isotope effects for hepatitis delta virus ribozyme catalysis that reveal a dissociative, metaphosphate-like transition state in stark contrast to the late, associative transition states observed for reactions catalyzed by specific base, Zn2+ ions, or ribonuclease A. This new information provides evidence for a discrete ribozyme active site design that modulates the RNA cleavage pathway to pass through an altered transition state.
Collapse
Affiliation(s)
- Benjamin Weissman
- Department of Chemistry and Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois 60637, United States
| | - Şölen Ekesan
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Hsuan-Chun Lin
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Shahbaz Gardezi
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Nan-Sheng Li
- Department of Chemistry and Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois 60637, United States
| | - Timothy J Giese
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Erika McCarthy
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Michael E Harris
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Joseph A Piccirilli
- Department of Chemistry and Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
2
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 327] [Impact Index Per Article: 54.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|
3
|
Šponer J, Krepl M, Banáš P, Kührová P, Zgarbová M, Jurečka P, Havrila M, Otyepka M. How to understand atomistic molecular dynamics simulations of RNA and protein-RNA complexes? WILEY INTERDISCIPLINARY REVIEWS-RNA 2016; 8. [PMID: 27863061 DOI: 10.1002/wrna.1405] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2016] [Revised: 09/13/2016] [Accepted: 10/10/2016] [Indexed: 01/01/2023]
Abstract
We provide a critical assessment of explicit-solvent atomistic molecular dynamics (MD) simulations of RNA and protein/RNA complexes, written primarily for non-specialists with an emphasis to explain the limitations of MD. MD simulations can be likened to hypothetical single-molecule experiments starting from single atomistic conformations and investigating genuine thermal sampling of the biomolecules. The main advantage of MD is the unlimited temporal and spatial resolution of positions of all atoms in the simulated systems. Fundamental limitations are the short physical time-scale of simulations, which can be partially alleviated by enhanced-sampling techniques, and the highly approximate atomistic force fields describing the simulated molecules. The applicability and present limitations of MD are demonstrated on studies of tetranucleotides, tetraloops, ribozymes, riboswitches and protein/RNA complexes. Wisely applied simulations respecting the approximations of the model can successfully complement structural and biochemical experiments. WIREs RNA 2017, 8:e1405. doi: 10.1002/wrna.1405 For further resources related to this article, please visit the WIREs website.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| |
Collapse
|
4
|
Casalino L, Magistrato A. Structural, dynamical and catalytic interplay between Mg2+ ions and RNA. Vices and virtues of atomistic simulations. Inorganica Chim Acta 2016. [DOI: 10.1016/j.ica.2016.02.011] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
5
|
Lee TS, Radak BK, Harris ME, York DM. A Two-Metal-Ion-Mediated Conformational Switching Pathway for HDV Ribozyme Activation. ACS Catal 2016; 6:1853-1869. [PMID: 27774349 DOI: 10.1021/acscatal.5b02158] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
RNA enzymes serve as a potentially powerful platform from which to design catalysts and engineer new biotechnology. A fundamental understanding of these systems provides insight to guide design. The hepatitis delta virus ribozyme (HDVr) is a small, self-cleaving RNA motif widely distributed in nature, that has served as a paradigm for understanding basic principles of RNA catalysis. Nevertheless, questions remain regarding the precise roles of divalent metal ions and key nucleotides in catalysis. In an effort to establish a reaction mechanism model consistent with available experimental data, we utilize molecular dynamics simulations to explore different conformations and metal ion binding modes along the HDVr reaction path. Building upon recent crystallographic data, our results provide a dynamic model of the HDVr reaction mechanism involving a conformational switch between multiple non-canonical G25:U20 base pair conformations in the active site. These local nucleobase dynamics play an important role in catalysis by modulating the metal binding environments of two Mg2+ ions that support catalysis at different steps of the reaction pathway. The first ion plays a structural role by inducing a base pair flip necessary to obtain the catalytic fold in which C75 moves towards to the scissile phosphate in the active site. Ejection of this ion then permits a second ion to bind elsewhere in the active site and facilitate nucleophile activation. The simulations collectively describe a mechanistic scenario that is consistent with currently available experimental data from crystallography, phosphorothioate substitutions, and chemical probing studies. Avenues for further experimental verification are suggested.
Collapse
Affiliation(s)
- Tai-Sung Lee
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Brian K. Radak
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
- Argonne National Laboratory, Argonne, Illinois 60439, United State
| | - Michael E. Harris
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, Ohio 44106, United States
| | - Darrin M. York
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| |
Collapse
|
6
|
Panteva MT, Giambaşu GM, York DM. Force Field for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) Ions That Have Balanced Interactions with Nucleic Acids. J Phys Chem B 2015; 119:15460-70. [PMID: 26583536 DOI: 10.1021/acs.jpcb.5b10423] [Citation(s) in RCA: 85] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Divalent metal ions are of fundamental importance to the function and folding of nucleic acids. Divalent metal ion-nucleic acid interactions are complex in nature and include both territorial and site specific binding. Commonly employed nonbonded divalent ion models, however, are often parametrized against bulk ion properties and are subsequently utilized in biomolecular simulations without considering any data related to interactions at specific nucleic acid sites. Previously, we assessed the ability of 17 different nonbonded Mg(2+) ion models to reproduce different properties of Mg(2+) in aqueous solution including radial distribution functions, solvation free energies, water exchange rates, and translational diffusion coefficients. In the present work, we depart from the recently developed 12-6-4 potential models for divalent metal ions developed by Li and Merz and tune the pairwise parameters for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) binding dimethyl phosphate, adenosine, and guanosine in order to reproduce experimental site specific binding free energies derived from potentiometric pH titration data. We further apply these parameters to investigate a metal ion migration previously proposed to occur during the catalytic reaction of the hammerhead ribozyme. The new parameters are shown to be accurate and balanced for nucleic acid binding in comparison with available experimental data and provide an important tool for molecular dynamics and free energy simulations of nucleic acids where these ions may exhibit different binding modes.
Collapse
Affiliation(s)
- Maria T Panteva
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University , 174 Frelinghuysen Road, Piscataway, New Jersey 08854-8076, United States
| | - George M Giambaşu
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University , 174 Frelinghuysen Road, Piscataway, New Jersey 08854-8076, United States
| | - Darrin M York
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University , 174 Frelinghuysen Road, Piscataway, New Jersey 08854-8076, United States
| |
Collapse
|
7
|
Radak BK, Lee TS, Harris ME, York DM. Assessment of metal-assisted nucleophile activation in the hepatitis delta virus ribozyme from molecular simulation and 3D-RISM. RNA (NEW YORK, N.Y.) 2015; 21:1566-1577. [PMID: 26170378 PMCID: PMC4536318 DOI: 10.1261/rna.051466.115] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/17/2015] [Accepted: 05/26/2015] [Indexed: 06/04/2023]
Abstract
The hepatitis delta virus ribozyme is an efficient catalyst of RNA 2'-O-transphosphorylation and has emerged as a key experimental system for identifying and characterizing fundamental features of RNA catalysis. Recent structural and biochemical data have led to a proposed mechanistic model whereby an active site Mg(2+) ion facilitates deprotonation of the O2' nucleophile, and a protonated cytosine residue (C75) acts as an acid to donate a proton to the O5' leaving group as noted in a previous study. This model assumes that the active site Mg(2+) ion forms an inner-sphere coordination with the O2' nucleophile and a nonbridging oxygen of the scissile phosphate. These contacts, however, are not fully resolved in the crystal structure, and biochemical data are not able to unambiguously exclude other mechanistic models. In order to explore the feasibility of this model, we exhaustively mapped the free energy surfaces with different active site ion occupancies via quantum mechanical/molecular mechanical (QM/MM) simulations. We further incorporate a three-dimensional reference interaction site model for the solvated ion atmosphere that allows these calculations to consider not only the rate associated with the chemical steps, but also the probability of observing the system in the presumed active state with the Mg(2+) ion bound. The QM/MM results predict that a pathway involving metal-assisted nucleophile activation is feasible based on the rate-controlling transition state barrier departing from the presumed metal-bound active state. However, QM/MM results for a similar pathway in the absence of Mg(2+) are not consistent with experimental data, suggesting that a structural model in which the crystallographically determined Mg(2+) is simply replaced with Na(+) is likely incorrect. It should be emphasized, however, that these results hinge upon the assumption of the validity of the presumed Mg(2+)-bound starting state, which has not yet been definitively verified experimentally, nor explored in depth computationally. Thus, further experimental and theoretical study is needed such that a consensus view of the catalytic mechanism emerges.
Collapse
Affiliation(s)
- Brian K Radak
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8076, USA Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Tai-Sung Lee
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8076, USA
| | - Michael E Harris
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, Ohio 44106, USA
| | - Darrin M York
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8076, USA
| |
Collapse
|
8
|
Kellerman DL, Simmons KS, Pedraza M, Piccirilli JA, York DM, Harris ME. Determination of hepatitis delta virus ribozyme N(-1) nucleobase and functional group specificity using internal competition kinetics. Anal Biochem 2015; 483:12-20. [PMID: 25937290 DOI: 10.1016/j.ab.2015.04.024] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2015] [Revised: 04/17/2015] [Accepted: 04/23/2015] [Indexed: 12/16/2022]
Abstract
Biological catalysis involves interactions distant from the site of chemistry that can position the substrate for reaction. Catalysis of RNA 2'-O-transphosphorylation by the hepatitis delta virus (HDV) ribozyme is sensitive to the identity of the N(-1) nucleotide flanking the reactive phosphoryl group. However, the interactions that affect the conformation of this position, and in turn the 2'O nucleophile, are unclear. Here, we describe the application of multiple substrate internal competition kinetic analyses to understand how the N(-1) nucleobase contributes to HDV catalysis and test the utility of this approach for RNA structure-function studies. Internal competition reactions containing all four substrate sequence variants at the N(-1) position in reactions using ribozyme active site mutations at A77 and A78 were used to test a proposed base-pairing interaction. Mutants A78U, A78G, and A79G retain significant catalytic activity but do not alter the specificity for the N(-1) nucleobase. Effects of nucleobase analog substitutions at N(-1) indicate that U is preferred due to the ability to donate an H-bond in the Watson-Crick face and avoid minor groove steric clash. The results provide information essential for evaluating models of the HDV active site and illustrate multiple substrate kinetic analyses as a practical approach for characterizing structure-function relationships in RNA reactions.
Collapse
Affiliation(s)
- Daniel L Kellerman
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA
| | - Kandice S Simmons
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA
| | - Mayra Pedraza
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA
| | - Joseph A Piccirilli
- Department of Chemistry and Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, IL 60637, USA
| | - Darrin M York
- Center for Integrative Proteomics Research, BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Michael E Harris
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA.
| |
Collapse
|
9
|
Review of computer simulations of isotope effects on biochemical reactions: From the Bigeleisen equation to Feynman's path integral. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2015; 1854:1782-94. [PMID: 25936775 DOI: 10.1016/j.bbapap.2015.04.021] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/01/2015] [Revised: 04/13/2015] [Accepted: 04/22/2015] [Indexed: 12/18/2022]
Abstract
Enzymatic reactions are integral components in many biological functions and malfunctions. The iconic structure of each reaction path for elucidating the reaction mechanism in details is the molecular structure of the rate-limiting transition state (RLTS). But RLTS is very hard to get caught or to get visualized by experimentalists. In spite of the lack of explicit molecular structure of the RLTS in experiment, we still can trace out the RLTS unique "fingerprints" by measuring the isotope effects on the reaction rate. This set of "fingerprints" is considered as a most direct probe of RLTS. By contrast, for computer simulations, oftentimes molecular structures of a number of TS can be precisely visualized on computer screen, however, theoreticians are not sure which TS is the actual rate-limiting one. As a result, this is an excellent stage setting for a perfect "marriage" between experiment and theory for determining the structure of RLTS, along with the reaction mechanism, i.e., experimentalists are responsible for "fingerprinting", whereas theoreticians are responsible for providing candidates that match the "fingerprints". In this Review, the origin of isotope effects on a chemical reaction is discussed from the perspectives of classical and quantum worlds, respectively (e.g., the origins of the inverse kinetic isotope effects and all the equilibrium isotope effects are purely from quantum). The conventional Bigeleisen equation for isotope effect calculations, as well as its refined version in the framework of Feynman's path integral and Kleinert's variational perturbation (KP) theory for systematically incorporating anharmonicity and (non-parabolic) quantum tunneling, are also presented. In addition, the outstanding interplay between theory and experiment for successfully deducing the RLTS structures and the reaction mechanisms is demonstrated by applications on biochemical reactions, namely models of bacterial squalene-to-hopene polycyclization and RNA 2'-O-transphosphorylation. For all these applications, we used our recently-developed path-integral method based on the KP theory, called automated integration-free path-integral (AIF-PI) method, to perform ab initio path-integral calculations of isotope effects. As opposed to the conventional path-integral molecular dynamics (PIMD) and Monte Carlo (PIMC) simulations, values calculated from our AIF-PI path-integral method can be as precise as (not as accurate as) the numerical precision of the computing machine. Lastly, comments are made on the general challenges in theoretical modeling of candidates matching the experimental "fingerprints" of RLTS. This article is part of a Special Issue entitled: Enzyme Transition States from Theory and Experiment.
Collapse
|
10
|
Thaplyal P, Ganguly A, Hammes-Schiffer S, Bevilacqua PC. Inverse thio effects in the hepatitis delta virus ribozyme reveal that the reaction pathway is controlled by metal ion charge density. Biochemistry 2015; 54:2160-75. [PMID: 25799319 PMCID: PMC4824481 DOI: 10.1021/acs.biochem.5b00190] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
![]()
The
hepatitis delta virus (HDV) ribozyme self-cleaves in the presence
of a wide range of monovalent and divalent ions. Prior theoretical
studies provided evidence that self-cleavage proceeds via a concerted
or stepwise pathway, with the outcome dictated by the valency of the
metal ion. In the present study, we measure stereospecific thio effects
at the nonbridging oxygens of the scissile phosphate under a wide
range of experimental conditions, including varying concentrations
of diverse monovalent and divalent ions, and combine these with quantum
mechanical/molecular mechanical (QM/MM) free energy simulations on
the stereospecific thio substrates. The RP substrate gives large normal thio effects in the presence of all
monovalent ions. The SP substrate also
gives normal or no thio effects, but only for smaller monovalent and
divalent cations, such as Li+, Mg2+, Ca2+, and Sr2+; in contrast, sizable inverse thio
effects are found for larger monovalent and divalent cations, including
Na+, K+, NH4+, and Ba2+. Proton inventories are found to be unity in the presence
of the larger monovalent and divalent ions, but two in the presence
of Mg2+. Additionally, rate–pH profiles are inverted
for the low charge density ions, and only imidazole plus ammonium
ions rescue an inactive C75Δ variant in the absence of Mg2+. Results from the thio effect experiments, rate–pH
profiles, proton inventories, and ammonium/imidazole rescue experiments,
combined with QM/MM free energy simulations, support a change in the
mechanism of HDV ribozyme self-cleavage from concerted and metal ion-stabilized
to stepwise and proton transfer-stabilized as the charge density of
the metal ion decreases.
Collapse
Affiliation(s)
- Pallavi Thaplyal
- †Department of Chemistry and Center for RNA Molecular Biology, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Abir Ganguly
- ‡Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Sharon Hammes-Schiffer
- ‡Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Philip C Bevilacqua
- †Department of Chemistry and Center for RNA Molecular Biology, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| |
Collapse
|
11
|
Panteva MT, Dissanayake T, Chen H, Radak BK, Kuechler ER, Giambaşu GM, Lee TS, York DM. Multiscale methods for computational RNA enzymology. Methods Enzymol 2015; 553:335-74. [PMID: 25726472 PMCID: PMC4739856 DOI: 10.1016/bs.mie.2014.10.064] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
RNA catalysis is of fundamental importance to biology and yet remains ill-understood due to its complex nature. The multidimensional "problem space" of RNA catalysis includes both local and global conformational rearrangements, changes in the ion atmosphere around nucleic acids and metal ion binding, dependence on potentially correlated protonation states of key residues, and bond breaking/forming in the chemical steps of the reaction. The goal of this chapter is to summarize and apply multiscale modeling methods in an effort to target the different parts of the RNA catalysis problem space while also addressing the limitations and pitfalls of these methods. Classical molecular dynamics simulations, reference interaction site model calculations, constant pH molecular dynamics (CpHMD) simulations, Hamiltonian replica exchange molecular dynamics, and quantum mechanical/molecular mechanical simulations will be discussed in the context of the study of RNA backbone cleavage transesterification. This reaction is catalyzed by both RNA and protein enzymes, and here we examine the different mechanistic strategies taken by the hepatitis delta virus ribozyme and RNase A.
Collapse
Affiliation(s)
- Maria T Panteva
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Thakshila Dissanayake
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Haoyuan Chen
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Brian K Radak
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Erich R Kuechler
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - George M Giambaşu
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Tai-Sung Lee
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Darrin M York
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA.
| |
Collapse
|
12
|
Estarellas C, Otyepka M, Koča J, Banáš P, Krepl M, Šponer J. Molecular dynamic simulations of protein/RNA complexes: CRISPR/Csy4 endoribonuclease. Biochim Biophys Acta Gen Subj 2014; 1850:1072-1090. [PMID: 25450173 DOI: 10.1016/j.bbagen.2014.10.021] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2014] [Revised: 10/15/2014] [Accepted: 10/20/2014] [Indexed: 12/20/2022]
Abstract
BACKGROUND Many prokaryotic genomes comprise Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) offering defense against foreign nucleic acids. These immune systems are conditioned by the production of small CRISPR-derived RNAs matured from long RNA precursors. This often requires a Csy4 endoribonuclease cleaving the RNA 3'-end. METHODS We report extended explicit solvent molecular dynamic (MD) simulations of Csy4/RNA complex in precursor and product states, based on X-ray structures of product and inactivated precursor (55 simulations; ~3.7μs in total). RESULTS The simulations identify double-protonated His29 and deprotonated terminal phosphate as the likely dominant protonation states consistent with the product structure. We revealed potential substates consistent with Ser148 and His29 acting as the general base and acid, respectively. The Ser148 could be straightforwardly deprotonated through solvent and could without further structural rearrangements deprotonate the nucleophile, contrasting similar studies investigating the general base role of nucleobases in ribozymes. We could not locate geometries consistent with His29 acting as general base. However, we caution that the X-ray structures do not always capture the catalytically active geometries and then the reactive structures may be unreachable by the simulation technique. CONCLUSIONS We identified potential catalytic arrangement of the Csy4/RNA complex but we also report limitations of the simulation technique. Even for the dominant protonation state we could not achieve full agreement between the simulations and the structural data. GENERAL SIGNIFICANCE Potential catalytic arrangement of the Csy4/RNA complex is found. Further, we provide unique insights into limitations of simulations of protein/RNA complexes, namely, the influence of the starting experimental structures and force field limitations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Carolina Estarellas
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jaroslav Koča
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic; Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic.
| |
Collapse
|
13
|
Kapral GJ, Jain S, Noeske J, Doudna JA, Richardson DC, Richardson JS. New tools provide a second look at HDV ribozyme structure, dynamics and cleavage. Nucleic Acids Res 2014; 42:12833-46. [PMID: 25326328 PMCID: PMC4227795 DOI: 10.1093/nar/gku992] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme essential for processing viral transcripts during rolling circle viral replication. The first crystal structure of the cleaved ribozyme was solved in 1998, followed by structures of uncleaved, mutant-inhibited and ion-complexed forms. Recently, methods have been developed that make the task of modeling RNA structure and dynamics significantly easier and more reliable. We have used ERRASER and PHENIX to rebuild and re-refine the cleaved and cis-acting C75U-inhibited structures of the HDV ribozyme. The results correct local conformations and identify alternates for RNA residues, many in functionally important regions, leading to improved R values and model validation statistics for both structures. We compare the rebuilt structures to a higher resolution, trans-acting deoxy-inhibited structure of the ribozyme, and conclude that although both inhibited structures are consistent with the currently accepted hammerhead-like mechanism of cleavage, they do not add direct structural evidence to the biochemical and modeling data. However, the rebuilt structures (PDBs: 4PR6, 4PRF) provide a more robust starting point for research on the dynamics and catalytic mechanism of the HDV ribozyme and demonstrate the power of new techniques to make significant improvements in RNA structures that impact biologically relevant conclusions.
Collapse
Affiliation(s)
- Gary J Kapral
- Department of Biochemistry, Duke University, Durham, NC 27710, USA
| | - Swati Jain
- Department of Biochemistry, Duke University, Durham, NC 27710, USA Program in Computational Biology and Bioinformatics, Duke University, Durham, NC 27710, USA
| | - Jonas Noeske
- Department of Molecular and Cell Biology, University of California Berkeley, Berkeley, CA 94720, USA
| | - Jennifer A Doudna
- Department of Molecular and Cell Biology, University of California Berkeley, Berkeley, CA 94720, USA
| | | | | |
Collapse
|
14
|
Sripathi KN, Tay WW, Banáš P, Otyepka M, Šponer J, Walter NG. Disparate HDV ribozyme crystal structures represent intermediates on a rugged free-energy landscape. RNA (NEW YORK, N.Y.) 2014; 20:1112-28. [PMID: 24854621 PMCID: PMC4114689 DOI: 10.1261/rna.044982.114] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/23/2014] [Accepted: 04/16/2014] [Indexed: 06/03/2023]
Abstract
The hepatitis delta virus (HDV) ribozyme is a member of the class of small, self-cleaving catalytic RNAs found in a wide range of genomes from HDV to human. Both pre- and post-catalysis (precursor and product) crystal structures of the cis-acting genomic HDV ribozyme have been determined. These structures, together with extensive solution probing, have suggested that a significant conformational change accompanies catalysis. A recent crystal structure of a trans-acting precursor, obtained at low pH and by molecular replacement from the previous product conformation, conforms to the product, raising the possibility that it represents an activated conformer past the conformational change. Here, using fluorescence resonance energy transfer (FRET), we discovered that cleavage of this ribozyme at physiological pH is accompanied by a structural lengthening in magnitude comparable to previous trans-acting HDV ribozymes. Conformational heterogeneity observed by FRET in solution appears to have been removed upon crystallization. Analysis of a total of 1.8 µsec of molecular dynamics (MD) simulations showed that the crystallographically unresolved cleavage site conformation is likely correctly modeled after the hammerhead ribozyme, but that crystal contacts and the removal of several 2'-oxygens near the scissile phosphate compromise catalytic in-line fitness. A cis-acting version of the ribozyme exhibits a more dynamic active site, while a G-1 residue upstream of the scissile phosphate favors poor fitness, allowing us to rationalize corresponding changes in catalytic activity. Based on these data, we propose that the available crystal structures of the HDV ribozyme represent intermediates on an overall rugged RNA folding free-energy landscape.
Collapse
Affiliation(s)
- Kamali N. Sripathi
- Department of Medicinal Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1065, USA
| | - Wendy W. Tay
- Program in Chemical Biology, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
| | - Pavel Banáš
- Regional Centre of Advance Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, 771 46 Olomouc, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advance Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, 612 65 Brno, Czech Republic
- Masaryk University, Campus Bohunice, 625 00 Brno, Czech Republic
| | - Nils G. Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
| |
Collapse
|
15
|
Wong KY, Xu Y, York DM. Ab initio path-integral calculations of kinetic and equilibrium isotope effects on base-catalyzed RNA transphosphorylation models. J Comput Chem 2014; 35:1302-16. [PMID: 24841935 PMCID: PMC4096342 DOI: 10.1002/jcc.23628] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2014] [Accepted: 04/06/2014] [Indexed: 12/16/2022]
Abstract
Detailed understandings of the reaction mechanisms of RNA catalysis in various environments can have profound importance for many applications, ranging from the design of new biotechnologies to the unraveling of the evolutionary origin of life. An integral step in the nucleolytic RNA catalysis is self-cleavage of RNA strands by 2'-O-transphosphorylation. Key to elucidating a reaction mechanism is determining the molecular structure and bonding characteristics of transition state. A direct and powerful probe of transition state is measuring isotope effects on biochemical reactions, particularly if we can reproduce isotope effect values from quantum calculations. This article significantly extends the scope of our previous joint experimental and theoretical work in examining isotope effects on enzymatic and nonenzymatic 2'-O-transphosphorylation reaction models that mimic reactions catalyzed by RNA enzymes (ribozymes), and protein enzymes such as ribonuclease A (RNase A). Native reactions are studied, as well as reactions with thio substitutions representing chemical modifications often used in experiments to probe mechanism. Here, we report and compare results from eight levels of electronic-structure calculations for constructing the potential energy surfaces in kinetic and equilibrium isotope effects (KIE and EIE) computations, including a "gold-standard" coupled-cluster level of theory [CCSD(T)]. In addition to the widely used Bigeleisen equation for estimating KIE and EIE values, internuclear anharmonicity and quantum tunneling effects were also computed using our recently developed ab initio path-integral method, that is, automated integration-free path-integral method. The results of this work establish an important set of benchmarks that serve to guide calculations of KIE and EIE for RNA catalysis.
Collapse
Affiliation(s)
- Kin-Yiu Wong
- Department of Physics, High Performance Cluster Computing Centre, Institute of Computational and Theoretical Studies, Hong Kong Baptist University, 224 Waterloo Road, Kowloon Tong, Hong Kong; Institute of Research and Continuing Education, Hong Kong Baptist University (Shenzhen), Shenzhen, China
| | | | | |
Collapse
|
16
|
Šponer J, Banáš P, Jurečka P, Zgarbová M, Kührová P, Havrila M, Krepl M, Stadlbauer P, Otyepka M. Molecular Dynamics Simulations of Nucleic Acids. From Tetranucleotides to the Ribosome. J Phys Chem Lett 2014; 5:1771-82. [PMID: 26270382 DOI: 10.1021/jz500557y] [Citation(s) in RCA: 113] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
We present a brief overview of explicit solvent molecular dynamics (MD) simulations of nucleic acids. We explain physical chemistry limitations of the simulations, namely, the molecular mechanics (MM) force field (FF) approximation and limited time scale. Further, we discuss relations and differences between simulations and experiments, compare standard and enhanced sampling simulations, discuss the role of starting structures, comment on different versions of nucleic acid FFs, and relate MM computations with contemporary quantum chemistry. Despite its limitations, we show that MD is a powerful technique for studying the structural dynamics of nucleic acids with a fast growing potential that substantially complements experimental results and aids their interpretation.
Collapse
Affiliation(s)
- Jiří Šponer
- †Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- ‡CEITEC - Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Pavel Banáš
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petr Jurečka
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Marie Zgarbová
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petra Kührová
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Marek Havrila
- †Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- ‡CEITEC - Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Miroslav Krepl
- †Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petr Stadlbauer
- †Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
17
|
Abstract
RNA dynamics play a fundamental role in many cellular functions. However, there is no general framework to describe these complex processes, which typically consist of many structural maneuvers that occur over timescales ranging from picoseconds to seconds. Here, we classify RNA dynamics into distinct modes representing transitions between basins on a hierarchical free-energy landscape. These transitions include large-scale secondary-structural transitions at >0.1-s timescales, base-pair/tertiary dynamics at microsecond-to-millisecond timescales, stacking dynamics at timescales ranging from nanoseconds to microseconds, and other "jittering" motions at timescales ranging from picoseconds to nanoseconds. We review various modes within these three different tiers, the different mechanisms by which they are used to regulate function, and how they can be coupled together to achieve greater functional complexity.
Collapse
|
18
|
Ganguly A, Thaplyal P, Rosta E, Bevilacqua PC, Hammes-Schiffer S. Quantum mechanical/molecular mechanical free energy simulations of the self-cleavage reaction in the hepatitis delta virus ribozyme. J Am Chem Soc 2014; 136:1483-96. [PMID: 24383543 PMCID: PMC3954522 DOI: 10.1021/ja4104217] [Citation(s) in RCA: 62] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
The
hepatitis delta virus (HDV) ribozyme catalyzes a self-cleavage
reaction using a combination of nucleobase and metal ion catalysis.
Both divalent and monovalent ions can catalyze this reaction, although
the rate is slower with monovalent ions alone. Herein, we use quantum
mechanical/molecular mechanical (QM/MM) free energy simulations to
investigate the mechanism of this ribozyme and to elucidate the roles
of the catalytic metal ion. With Mg2+ at the catalytic
site, the self-cleavage mechanism is observed to be concerted with
a phosphorane-like transition state and a free energy barrier of ∼13
kcal/mol, consistent with free energy barrier values extrapolated
from experimental studies. With Na+ at the catalytic site,
the mechanism is observed to be sequential, passing through a phosphorane
intermediate, with free energy barriers of 2–4 kcal/mol for
both steps; moreover, proton transfer from the exocyclic amine of
protonated C75 to the nonbridging oxygen of the scissile phosphate
occurs to stabilize the phosphorane intermediate in the sequential
mechanism. To explain the slower rate observed experimentally with
monovalent ions, we hypothesize that the activation of the O2′
nucleophile by deprotonation and orientation is less favorable with
Na+ ions than with Mg2+ ions. To explore this
hypothesis, we experimentally measure the pKa of O2′ by kinetic and NMR methods and find it to be
lower in the presence of divalent ions rather than only monovalent
ions. The combined theoretical and experimental results indicate that
the catalytic Mg2+ ion may play three key roles: assisting
in the activation of the O2′ nucleophile, acidifying the general
acid C75, and stabilizing the nonbridging oxygen to prevent proton
transfer to it.
Collapse
Affiliation(s)
- Abir Ganguly
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | | | | | | | | |
Collapse
|
19
|
Havrila M, Réblová K, Zirbel CL, Leontis NB, Šponer J. Isosteric and nonisosteric base pairs in RNA motifs: molecular dynamics and bioinformatics study of the sarcin-ricin internal loop. J Phys Chem B 2013; 117:14302-19. [PMID: 24144333 PMCID: PMC3946555 DOI: 10.1021/jp408530w] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The sarcin-ricin RNA motif (SR motif) is one of the most prominent recurrent RNA building blocks that occurs in many different RNA contexts and folds autonomously, that is, in a context-independent manner. In this study, we combined bioinformatics analysis with explicit-solvent molecular dynamics (MD) simulations to better understand the relation between the RNA sequence and the evolutionary patterns of the SR motif. A SHAPE probing experiment was also performed to confirm the fidelity of the MD simulations. We identified 57 instances of the SR motif in a nonredundant subset of the RNA X-ray structure database and analyzed their base pairing, base-phosphate, and backbone-backbone interactions. We extracted sequences aligned to these instances from large rRNA alignments to determine the frequency of occurrence for different sequence variants. We then used a simple scoring scheme based on isostericity to suggest 10 sequence variants with a highly variable expected degree of compatibility with the SR motif 3D structure. We carried out MD simulations of SR motifs with these base substitutions. Nonisosteric base substitutions led to unstable structures, but so did isosteric substitutions which were unable to make key base-phosphate interactions. The MD technique explains why some potentially isosteric SR motifs are not realized during evolution. We also found that the inability to form stable cWW geometry is an important factor in the case of the first base pair of the flexible region of the SR motif. A comparison of structural, bioinformatics, SHAPE probing, and MD simulation data reveals that explicit solvent MD simulations neatly reflect the viability of different sequence variants of the SR motif. Thus, MD simulations can efficiently complement bioinformatics tools in studies of conservation patterns of RNA motifs and provide atomistic insight into the role of their different signature interactions.
Collapse
Affiliation(s)
- Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Kamila Réblová
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Craig L. Zirbel
- Department of Chemistry, Bowling Green State University, Bowling Green, OH 43403, USA
- Department of Mathematics and Statistics, Bowling Green State University, Bowling Green, OH 43403, USA
| | - Neocles B. Leontis
- Department of Chemistry, Bowling Green State University, Bowling Green, OH 43403, USA
- Department of Mathematics and Statistics, Bowling Green State University, Bowling Green, OH 43403, USA
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
20
|
Experimental and computational analysis of the transition state for ribonuclease A-catalyzed RNA 2'-O-transphosphorylation. Proc Natl Acad Sci U S A 2013; 110:13002-7. [PMID: 23878223 DOI: 10.1073/pnas.1215086110] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Enzymes function by stabilizing reaction transition states; therefore, comparison of the transition states of enzymatic and nonenzymatic model reactions can provide insight into biological catalysis. Catalysis of RNA 2'-O-transphosphorylation by ribonuclease A is proposed to involve electrostatic stabilization and acid/base catalysis, although the structure of the rate-limiting transition state is uncertain. Here, we describe coordinated kinetic isotope effect (KIE) analyses, molecular dynamics simulations, and quantum mechanical calculations to model the transition state and mechanism of RNase A. Comparison of the (18)O KIEs on the 2'O nucleophile, 5'O leaving group, and nonbridging phosphoryl oxygens for RNase A to values observed for hydronium- or hydroxide-catalyzed reactions indicate a late anionic transition state. Molecular dynamics simulations using an anionic phosphorane transition state mimic suggest that H-bonding by protonated His12 and Lys41 stabilizes the transition state by neutralizing the negative charge on the nonbridging phosphoryl oxygens. Quantum mechanical calculations consistent with the experimental KIEs indicate that expulsion of the 5'O remains an integral feature of the rate-limiting step both on and off the enzyme. Electrostatic interactions with positively charged amino acid site chains (His12/Lys41), together with proton transfer from His119, render departure of the 5'O less advanced compared with the solution reaction and stabilize charge buildup in the transition state. The ability to obtain a chemically detailed description of 2'-O-transphosphorylation transition states provides an opportunity to advance our understanding of biological catalysis significantly by determining how the catalytic modes and active site environments of phosphoryl transferases influence transition state structure.
Collapse
|
21
|
Wrzesinski J, Błaszczyk L, Wrońska M, Kasprowicz A, Stokowa-Sołtys K, Nagaj J, Szafraniec M, Kulinski T, Jeżowska-Bojczuk M, Ciesiołka J. Mapping the interactions of selected antibiotics and their Cu2+complexes with the antigenomic delta ribozyme. FEBS J 2013; 280:2652-64. [DOI: 10.1111/febs.12257] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2013] [Revised: 03/15/2013] [Accepted: 03/18/2013] [Indexed: 11/30/2022]
Affiliation(s)
- Jan Wrzesinski
- Institute of Bioorganic Chemistry; Polish Academy of Sciences; Poznań; Poland
| | - Leszek Błaszczyk
- Institute of Bioorganic Chemistry; Polish Academy of Sciences; Poznań; Poland
| | | | | | | | | | - Milena Szafraniec
- Institute of Bioorganic Chemistry; Polish Academy of Sciences; Poznań; Poland
| | - Tadeusz Kulinski
- Institute of Bioorganic Chemistry; Polish Academy of Sciences; Poznań; Poland
| | | | - Jerzy Ciesiołka
- Institute of Bioorganic Chemistry; Polish Academy of Sciences; Poznań; Poland
| |
Collapse
|
22
|
Guo Y, Zhang W. Molecular dynamics simulation of RNA pseudoknot unfolding pathway. WUHAN UNIVERSITY JOURNAL OF NATURAL SCIENCES 2013. [PMCID: PMC7149040 DOI: 10.1007/s11859-013-0905-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Many biological functions of RNA molecules are related to their pseudoknot structures. It is significant for predicting the structure and function of RNA that learning about the stability and the process of RNA pseudoknot folding and unfolding. The structural features of mouse mammary tumor virus (MMTV) RNA pseudoknot in different ion concentration, the unfolding process of the RNA pseudoknot, and the two hairpin helices that constitute the RNA pseudoknot were studied with all atom molecule dynamics simulation method in this paper. We found that the higher cation concentration can cause structure of the RNA molecules more stable, and ions played an indispensable role in keeping the structure of RNA molecules stable; the unfolding process of hairpin structure was corresponding to the antiprocess of its folding process. The main pathway of pseudoknot unfolding was that the inner base pair opened first, and then, the two helices, which formed the RNA pseudoknot opened decussately, while the folding pathway of the RNA pseudoknot was a helix folding after formation of the other helix. Therefore, the unfolding process of RNA pseudoknot is different from the antiprocess of its folding process, and the unfolding process of each helix in the RNA pseudoknot is similar to the hairpin structure’s unfolding process, which means that both are the unzipping process.
Collapse
|
23
|
Abstract
The role of pH-dependent protonation equilibrium in modulating RNA dynamics and function is one of the key unanswered questions in RNA biology. Molecular dynamics (MD) simulations can provide insight into the mechanistic roles of protonated nucleotides, but it is only capable of modeling fixed protonation states and requires prior knowledge of the key residue's protonation state. Recently, we developed a framework for constant pH molecular dynamics simulations (CPHMDMSλD) of nucleic acids, where the nucleotides' protonation states are modeled as dynamic variables that are coupled to the structural dynamics of the RNA. In the present study, we demonstrate the application of CPHMDMSλD to the lead-dependent ribozyme; establishing the validity of this approach for modeling complex RNA structures. We show that CPHMDMSλD accurately predicts the direction of the pKa shifts and reproduces experimentally-measured microscopic pKa values with an average unsigned error of 1.3 pKa units. The effects of coupled titration states in RNA structures are modeled, and the importance of conformation sampling is highlighted. The general accuracy of CPHMDMSλD simulations in reproducing pH-dependent observables reported in this work demonstrates that constant pH simulations provides a powerful tool to investigate pH-dependent processes in nucleic acids.
Collapse
Affiliation(s)
- Garrett B Goh
- Department of Chemistry, University of Michigan, 930 N. University, Ann Arbor, Michigan 48109, United States
| | | | | |
Collapse
|
24
|
Rouleau SG, Jodoin R, Bisaillon M, Perreault JP. Programming a highly structured ribozyme into complex allostery using RNA oligonucleotides. ACS Chem Biol 2012; 7:1802-6. [PMID: 22950875 DOI: 10.1021/cb300319m] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
RNA possesses great potential for expanding the toolbox currently available to synthetic biologists. Here, the modulation of the Hepatitis Delta Virus ribozyme's activity with a series of rationally designed aptamers and effector RNA oligonucleotides is described. The ribozyme was initially fused with an 18-nucleotide hairpin structure that abolished its self-cleaving activity. The binding of a 14-mer oligonucleotide to the hairpin rescued the self-cleavage in a concentration-dependent manner. This modified ribozyme was inserted into the 5' UTR of a reporter gene, and the resulting construct was used to demonstrate that it is possible to modulate the ribozyme activity in cellulo with the oligonucleotide. Subsequently, ribozymes possessing specialized aptamers respecting other logic gates were also successfully designed and found to be functional in vitro. To our knowledge, this is the first example of HDV ribozyme regulation by oligonucleotides, as well as the first allosteric regulation of HDV ribozyme in mammalian cells.
Collapse
Affiliation(s)
- Samuel G. Rouleau
- RNA group/Groupe
ARN, Département de biochimie,
Faculté de médecine et des sciences de la santé,
Pavillon de Recherche Appliquée sur le Cancer, Université de Sherbrooke, 3201 Jean-Mignault,
Sherbrooke, Québec J1E 4K8, Canada
| | - Rachel Jodoin
- RNA group/Groupe
ARN, Département de biochimie,
Faculté de médecine et des sciences de la santé,
Pavillon de Recherche Appliquée sur le Cancer, Université de Sherbrooke, 3201 Jean-Mignault,
Sherbrooke, Québec J1E 4K8, Canada
| | - Martin Bisaillon
- RNA group/Groupe
ARN, Département de biochimie,
Faculté de médecine et des sciences de la santé,
Pavillon de Recherche Appliquée sur le Cancer, Université de Sherbrooke, 3201 Jean-Mignault,
Sherbrooke, Québec J1E 4K8, Canada
| | - Jean-Pierre Perreault
- RNA group/Groupe
ARN, Département de biochimie,
Faculté de médecine et des sciences de la santé,
Pavillon de Recherche Appliquée sur le Cancer, Université de Sherbrooke, 3201 Jean-Mignault,
Sherbrooke, Québec J1E 4K8, Canada
| |
Collapse
|
25
|
Roca M, Aranda J, Moliner V, Tuñón I. Modeling methods for studying post-translational and transcriptional modifying enzymes. Curr Opin Chem Biol 2012; 16:465-71. [PMID: 23127358 DOI: 10.1016/j.cbpa.2012.10.014] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2012] [Revised: 09/20/2012] [Accepted: 10/11/2012] [Indexed: 10/27/2022]
Abstract
Biological catalysis is a complex chemical process that involves not only electronic reorganization in the substrate but also the reorganization of the catalyst. This complexity is even larger in the case of post-transcriptional and post-translational modifications which may involve the interaction between two biomacromolecules. However, the development over the past decades of new computational methods and strategies is offering a detailed molecular picture of the catalytic event and a deep understanding of the mechanisms of chemical reactions in biological environments. Here we review the efforts made in the last years to model catalysis in post-transcriptional and post-translational processes. We stress on the advantages and problems of the different computational strategies and their applicability in different cases.
Collapse
Affiliation(s)
- Maite Roca
- Departamento de Química Física, Universitat de València, 46100 Burjassot, Spain
| | | | | | | |
Collapse
|
26
|
Banáš P, Sklenovský P, Wedekind JE, Šponer J, Otyepka M. Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study. J Phys Chem B 2012; 116:12721-34. [PMID: 22998634 PMCID: PMC3505677 DOI: 10.1021/jp309230v] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Riboswitches often occur in the 5'-untranslated regions of bacterial mRNA where they regulate gene expression. The preQ(1) riboswitch controls the biosynthesis of a hypermodified nucleoside queuosine in response to binding the queuosine metabolic intermediate. Structures of the ligand-bound and ligand-free states of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis were determined recently by X-ray crystallography. We used multiple, microsecond-long molecular dynamics simulations (29 μs in total) to characterize the structural dynamics of preQ(1) riboswitches in both states. We observed different stabilities of the stem in the bound and free states, resulting in different accessibilities of the ribosome-binding site. These differences are related to different stacking interactions between nucleotides of the stem and the associated loop, which itself adopts different conformations in the bound and free states. We suggest that the loop not only serves to bind preQ(1) but also transmits information about ligand binding from the ligand-binding pocket to the stem, which has implications for mRNA accessibility to the ribosome. We explain functional results obscured by a high salt crystallization medium and help to refine regions of disordered electron density, which demonstrates the predictive power of our approach. Besides investigating the functional dynamics of the riboswitch, we have also utilized this unique small folded RNA system for analysis of performance of the RNA force field on the μs time scale. The latest AMBER parmbsc0χ(OL3) RNA force field is capable of providing stable trajectories of the folded molecule on the μs time scale. On the other hand, force fields that are not properly balanced lead to significant structural perturbations on the sub-μs time scale, which could easily lead to inappropriate interpretation of the simulation data.
Collapse
Affiliation(s)
- Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Petr Sklenovský
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Joseph E. Wedekind
- Department of Biochemistry and Biophysics, University of Rochester School of Medicine and Dentistry, 601 Elmwood Ave, Box 712, Rochester, NY 14620, USA
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
- CEITEC – Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
27
|
Bowman JC, Lenz TK, Hud NV, Williams LD. Cations in charge: magnesium ions in RNA folding and catalysis. Curr Opin Struct Biol 2012; 22:262-72. [PMID: 22595008 DOI: 10.1016/j.sbi.2012.04.006] [Citation(s) in RCA: 139] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2012] [Revised: 04/24/2012] [Accepted: 04/24/2012] [Indexed: 12/22/2022]
Affiliation(s)
- Jessica C Bowman
- School of Chemistry and Biochemistry, Parker H. Petit Institute for Bioengineering and Bioscience, Center for Ribosomal Origins and Evolution, Georgia Institute of Technology, Atlanta, GA 30332-0400, United States
| | | | | | | |
Collapse
|
28
|
Ganguly A, Bevilacqua PC, Hammes-Schiffer S. Quantum Mechanical/Molecular Mechanical Study of the HDV Ribozyme: Impact of the Catalytic Metal Ion on the Mechanism. J Phys Chem Lett 2011; 2:2906-2911. [PMID: 22163069 PMCID: PMC3233192 DOI: 10.1021/jz2013215] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
A recent crystal structure of the precleaved HDV ribozyme along with biochemical data support a mechanism for phosphodiester bond self-cleavage in which C75 acts as a general acid and bound Mg(2+) ion acts as a Lewis acid. Herein this precleaved crystal structure is used as the basis for quantum mechanical/molecular mechanical calculations. These calculations indicate that the self-cleavage reaction is concerted with a phosphorane-like transition state when a divalent ion, Mg(2+) or Ca(2+), is bound at the catalytic site but is sequential with a phosphorane intermediate when a monovalent ion, such as Na(+), is at this site. Electrostatic potential calculations suggest that the divalent metal ion at the catalytic site lowers the pK(a) of C75, leading to the concerted mechanism in which the proton is partially transferred to the leaving group in the phosphorane-like transition state. These observations are consistent with experimental data, including pK(a) measurements, reaction kinetics, and proton inventories with divalent and monovalent ions.
Collapse
|