1
|
Zhang Y, Zhang J, Wan H, Wu Z, Xu H, Zhang Z, Wang Y, Wang J. New Insights into the Dependence of CPEB3 Ribozyme Cleavage on Mn 2+ and Mg 2. J Phys Chem Lett 2024; 15:2708-2714. [PMID: 38427973 DOI: 10.1021/acs.jpclett.3c03221] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2024]
Abstract
CPEB3 ribozyme is a self-cleaving RNA that occurs naturally in mammals and requires divalent metal ions for efficient activity. Ribozymes exhibit preferences for specific metal ions, but the exact differences in the catalytic mechanisms of various metal ions on the CPEB3 ribozyme remain unclear. Our findings reveal that Mn2+ functions as a more effective cofactor for CPEB3 ribozyme catalysis compared to Mg2+, as confirmed by its stronger binding affinity to CPEB3 by EPR. Cleavage assays of CPEB3 mutants and molecular docking analyses further showed that excessive Mn2+ ions can bind to a second binding site near the catalytic site, hindering CPEB3 catalytic efficiency and contributing to the Mn2+ bell-shaped curve. These results implicate a pivotal role for the local nucleobase-Mn2+ interactions in facilitating RNA folding and modulating the directed attack of nucleophilic reagents. Our study provides new insights and experimental evidence for exploring the divalent cation dependent cleavage mechanism of the CPEB3 ribozyme.
Collapse
Affiliation(s)
- Yaoyao Zhang
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jing Zhang
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Hengjia Wan
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Ziwei Wu
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Huangtao Xu
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
| | - Zhe Zhang
- Institute for Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Yujuan Wang
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- International Magnetobiology Frontier Research Center (iMFRC), Science Island, Hefei 230031, China
| | - Junfeng Wang
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- University of Science and Technology of China, Hefei, Anhui 230026, China
- School of Basic Medical Sciences, Anhui Medical University, Hefei, Anhui 230032, China
- International Magnetobiology Frontier Research Center (iMFRC), Science Island, Hefei 230031, China
| |
Collapse
|
2
|
Chen CC, Han J, Chinn CA, Rounds JS, Li X, Nikan M, Myszka M, Tong L, Passalacqua LFM, Bredy T, Wood MA, Luptak A. Inhibition of Cpeb3 ribozyme elevates CPEB3 protein expression and polyadenylation of its target mRNAs and enhances object location memory. eLife 2024; 13:e90116. [PMID: 38319152 PMCID: PMC10919898 DOI: 10.7554/elife.90116] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Accepted: 02/05/2024] [Indexed: 02/07/2024] Open
Abstract
A self-cleaving ribozyme that maps to an intron of the cytoplasmic polyadenylation element-binding protein 3 (Cpeb3) gene is thought to play a role in human episodic memory, but the underlying mechanisms mediating this effect are not known. We tested the activity of the murine sequence and found that the ribozyme's self-scission half-life matches the time it takes an RNA polymerase to reach the immediate downstream exon, suggesting that the ribozyme-dependent intron cleavage is tuned to co-transcriptional splicing of the Cpeb3 mRNA. Our studies also reveal that the murine ribozyme modulates maturation of its harboring mRNA in both cultured cortical neurons and the hippocampus: inhibition of the ribozyme using an antisense oligonucleotide leads to increased CPEB3 protein expression, which enhances polyadenylation and translation of localized plasticity-related target mRNAs, and subsequently strengthens hippocampal-dependent long-term memory. These findings reveal a previously unknown role for self-cleaving ribozyme activity in regulating experience-induced co-transcriptional and local translational processes required for learning and memory.
Collapse
Affiliation(s)
- Claire C Chen
- Department of Pharmaceutical Sciences, University of California, IrvineIrvineUnited States
| | - Joseph Han
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California, IrvineIrvineUnited States
| | - Carlene A Chinn
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California, IrvineIrvineUnited States
| | - Jacob S Rounds
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California, IrvineIrvineUnited States
| | - Xiang Li
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California, IrvineIrvineUnited States
| | | | - Marie Myszka
- Department of Chemistry, University of California, IrvineIrvineUnited States
| | - Liqi Tong
- Institute for Memory Impairments and Neurological Disorders, University of California, IrvineIrvineUnited States
| | - Luiz FM Passalacqua
- Department of Pharmaceutical Sciences, University of California, IrvineIrvineUnited States
| | - Timothy Bredy
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California, IrvineIrvineUnited States
| | - Marcelo A Wood
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California, IrvineIrvineUnited States
| | - Andrej Luptak
- Department of Pharmaceutical Sciences, University of California, IrvineIrvineUnited States
- Department of Chemistry, University of California, IrvineIrvineUnited States
- Department of Molecular Biology and Biochemistry, University of California, IrvineIrvineUnited States
| |
Collapse
|
3
|
Das R, Kretsch RC, Simpkin AJ, Mulvaney T, Pham P, Rangan R, Bu F, Keegan RM, Topf M, Rigden DJ, Miao Z, Westhof E. Assessment of three-dimensional RNA structure prediction in CASP15. Proteins 2023; 91:1747-1770. [PMID: 37876231 PMCID: PMC10841292 DOI: 10.1002/prot.26602] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Revised: 08/21/2023] [Accepted: 09/07/2023] [Indexed: 10/26/2023]
Abstract
The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty-two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and x-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as noncanonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.
Collapse
Affiliation(s)
- Rhiju Das
- Department of Biochemistry, Stanford University School of Medicine, CA USA
- Biophysics Program, Stanford University School of Medicine, CA USA
- Howard Hughes Medical Institute, Stanford University, CA USA
| | | | - Adam J. Simpkin
- Institute of Systems, Molecular & Integrative Biology, The University of Liverpool, UK
| | - Thomas Mulvaney
- Centre for Structural Systems Biology (CSSB), Leibniz-Institut für Virologie (LIV), Hamburg, Germany
- University Medical Center Hamburg-Eppendorf (UKE), Hamburg, Germany
| | - Phillip Pham
- Department of Biochemistry, Stanford University School of Medicine, CA USA
| | - Ramya Rangan
- Biophysics Program, Stanford University School of Medicine, CA USA
| | - Fan Bu
- Guangzhou Laboratory, Guangzhou International Bio Island, Guangzhou 510005, China
- Division of Life Sciences and Medicine, University of Science and Technology of China, Hefei 230036, Anhui, China
| | - Ronan M. Keegan
- Institute of Systems, Molecular & Integrative Biology, The University of Liverpool, UK
- Life Science, Diamond Light Source, Harwell Science, UK
| | - Maya Topf
- Centre for Structural Systems Biology (CSSB), Leibniz-Institut für Virologie (LIV), Hamburg, Germany
- University Medical Center Hamburg-Eppendorf (UKE), Hamburg, Germany
| | - Daniel J. Rigden
- Institute of Systems, Molecular & Integrative Biology, The University of Liverpool, UK
| | - Zhichao Miao
- GMU-GIBH Joint School of Life Sciences, The Guangdong-Hong Kong-Macau Joint Laboratory for Cell Fate Regulation and Diseases, Guangzhou National Laboratory, Guangzhou Medical University
- Shanghai Key Laboratory of Anesthesiology and Brain Functional Modulation, Clinical Research Center for Anesthesiology and Perioperative Medicine, Translational Research Institute of Brain and Brain-Like Intelligence, Shanghai Fourth People's Hospital, School of Medicine, Tongji University, Shanghai 200434, China
| | - Eric Westhof
- Architecture et Réactivité de l’ARN, Institut de Biologie Moléculaire et Cellulaire du CNRS, Université de Strasbourg, F-67084, Strasbourg, France
| |
Collapse
|
4
|
Chen CC, Han J, Chinn CA, Rounds JS, Li X, Nikan M, Myszka M, Tong L, Passalacqua LFM, Bredy TW, Wood MA, Lupták A. Inhibition of CPEB3 ribozyme elevates CPEB3 protein expression and polyadenylation of its target mRNAs, and enhances object location memory. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.06.07.543953. [PMID: 37333407 PMCID: PMC10274809 DOI: 10.1101/2023.06.07.543953] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/20/2023]
Abstract
A self-cleaving ribozyme that maps to an intron of the cytoplasmic polyadenylation element binding protein 3 (CPEB3) gene is thought to play a role in human episodic memory, but the underlying mechanisms mediating this effect are not known. We tested the activity of the murine sequence and found that the ribozyme's self-scission half-life matches the time it takes an RNA polymerase to reach the immediate downstream exon, suggesting that the ribozyme-dependent intron cleavage is tuned to co-transcriptional splicing of the CPEB3 mRNA. Our studies also reveal that the murine ribozyme modulates maturation of its harboring mRNA in both cultured cortical neurons and the hippocampus: inhibition of the ribozyme using an antisense oligonucleotide leads to increased CPEB3 protein expression, which enhances polyadenylation and translation of localized plasticity-related target mRNAs, and subsequently strengthens hippocampal-dependent long-term memory. These findings reveal a previously unknown role for self-cleaving ribozyme activity in regulating experience-induced co-transcriptional and local translational processes required for learning and memory.
Collapse
Affiliation(s)
- Claire C. Chen
- Department of Pharmaceutical Sciences, University of California–Irvine, Irvine, California 92697, United States
| | - Joseph Han
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California–Irvine, Irvine, California 92697, United States
| | - Carlene A. Chinn
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California–Irvine, Irvine, California 92697, United States
| | - Jacob S. Rounds
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California–Irvine, Irvine, California 92697, United States
| | - Xiang Li
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California–Irvine, Irvine, California 92697, United States
| | - Mehran Nikan
- Ionis Pharmaceuticals, 2855 Gazelle Court, Carlsbad, CA 92010, USA
| | - Marie Myszka
- Department of Chemistry, University of California–Irvine, Irvine, California 92697, United States
| | - Liqi Tong
- Institute for Memory Impairments and Neurological Disorders, University of California–Irvine, Irvine, California 92697, United States
| | - Luiz F. M. Passalacqua
- Department of Pharmaceutical Sciences, University of California–Irvine, Irvine, California 92697, United States
| | - Timothy W. Bredy
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California–Irvine, Irvine, California 92697, United States
| | - Marcelo A. Wood
- Department of Neurobiology and Behavior, Center for the Neurobiology of Learning and Memory, University of California–Irvine, Irvine, California 92697, United States
| | - Andrej Lupták
- Department of Pharmaceutical Sciences, University of California–Irvine, Irvine, California 92697, United States
- Department of Chemistry, University of California–Irvine, Irvine, California 92697, United States
- Department of Molecular Biology and Biochemistry, University of California–Irvine, Irvine, California 92697, United States
| |
Collapse
|
5
|
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
|
6
|
Sieg JP, McKinley LN, Huot MJ, Yennawar NH, Bevilacqua PC. The Metabolome Weakens RNA Thermodynamic Stability and Strengthens RNA Chemical Stability. Biochemistry 2022; 61:2579-2591. [PMID: 36306436 PMCID: PMC9669196 DOI: 10.1021/acs.biochem.2c00488] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We examined the complex network of interactions among RNA, the metabolome, and divalent Mg2+ under conditions that mimic the Escherichia coli cytoplasm. We determined Mg2+ binding constants for the top 15 E. coli metabolites, comprising 80% of the metabolome by concentration at physiological pH and monovalent ion concentrations. These data were used to inform the development of an artificial cytoplasm that mimics in vivo E. coli conditions, which we term "Eco80". We empirically determined that the mixture of E. coli metabolites in Eco80 approximated single-site binding behavior toward Mg2+ in the biologically relevant free Mg2+ range of ∼0.5 to 3 mM Mg2+, using a Mg2+-sensitive fluorescent dye. Effects of Eco80 conditions on the thermodynamic stability, chemical stability, structure, and catalysis of RNA were examined. We found that Eco80 conditions lead to opposing effects on the thermodynamic and chemical stabilities of RNA. In particular, the thermodynamic stability of RNA helices was weakened by 0.69 ± 0.12 kcal/mol, while the chemical stability was enhanced ∼2-fold, which can be understood using the speciation of Mg2+ between weak and strong Mg2+-metabolite complexes in Eco80. Overall, the use of Eco80 reflects RNA function in vivo and enhances the biological relevance of mechanistic studies of RNA.
Collapse
Affiliation(s)
- Jacob P. Sieg
- Department of Chemistry, Pennsylvania State University, University Park, PA 16802
- Center for RNA Molecular Biology, Pennsylvania State University, University Park, PA 16802
| | - Lauren N. McKinley
- Department of Chemistry, Pennsylvania State University, University Park, PA 16802
- Center for RNA Molecular Biology, Pennsylvania State University, University Park, PA 16802
| | - Melanie J. Huot
- Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, PA 16802
- Department of Biology, Pennsylvania State University, University Park, PA 16802
| | - Neela H. Yennawar
- The Huck Institutes of the Life Sciences, The Pennsylvania State University, University Park, PA 16802
| | - Philip C. Bevilacqua
- Department of Chemistry, Pennsylvania State University, University Park, PA 16802
- Center for RNA Molecular Biology, Pennsylvania State University, University Park, PA 16802
- Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, PA 16802
| |
Collapse
|
7
|
Abstract
In recent years, it has become clear that RNA molecules are involved in almost all vital cellular processes and pathogenesis of human disorders. The functional diversity of RNA comes from its structural richness. Although composed of only four nucleotides, RNA molecules present a plethora of secondary and tertiary structures critical for intra and intermolecular contacts with other RNAs and ligands (proteins, small metabolites, etc.). In order to fully understand RNA function it is necessary to define its spatial structure. Crystallography, nuclear magnetic resonance and cryogenic electron microscopy have demonstrated considerable success in determining the structures of biologically important RNA molecules. However, these powerful methods require large amounts of sample. Despite their limitations, chemical synthesis and in vitro transcription are usually employed to obtain milligram quantities of RNA for structural studies, delivering simple and effective methods for large-scale production of homogenous samples. The aim of this paper is to provide an overview of methods for large-scale RNA synthesis with emphasis on chemical synthesis and in vitro transcription. We also present our own results of testing the efficiency of these approaches in order to adapt the material acquisition strategy depending on the desired RNA construct.
Collapse
|
8
|
Ruszkowska A, Zheng YY, Mao S, Ruszkowski M, Sheng J. Structural Insights Into the 5′UG/3′GU Wobble Tandem in Complex With Ba2+ Cation. Front Mol Biosci 2022; 8:762786. [PMID: 35096964 PMCID: PMC8793689 DOI: 10.3389/fmolb.2021.762786] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2021] [Accepted: 11/24/2021] [Indexed: 11/25/2022] Open
Abstract
G•U wobble base pair frequently occurs in RNA structures. The unique chemical, thermodynamic, and structural properties of the G•U pair are widely exploited in RNA biology. In several RNA molecules, the G•U pair plays key roles in folding, ribozyme catalysis, and interactions with proteins. G•U may occur as a single pair or in tandem motifs with different geometries, electrostatics, and thermodynamics, further extending its biological functions. The metal binding affinity, which is essential for RNA folding, catalysis, and other interactions, differs with respect to the tandem motif type due to the different electrostatic potentials of the major grooves. In this work, we present the crystal structure of an RNA 8-mer duplex r[UCGUGCGA]2, providing detailed structural insights into the tandem motif I (5′UG/3′GU) complexed with Ba2+ cation. We compare the electrostatic potential of the presented motif I major groove with previously published structures of tandem motifs I, II (5′GU/3′UG), and III (5′GG/3′UU). A local patch of a strongly negative electrostatic potential in the major groove of the presented structure forms the metal binding site with the contributions of three oxygen atoms from the tandem. These results give us a better understanding of the G•U tandem motif I as a divalent metal binder, a feature essential for RNA functions.
Collapse
Affiliation(s)
| | - Ya Ying Zheng
- Department of Chemistry, The RNA Institute, University at Albany, State University of New York, Albany, NY, United States
| | - Song Mao
- Department of Chemistry, The RNA Institute, University at Albany, State University of New York, Albany, NY, United States
| | - Milosz Ruszkowski
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, Poznan, Poland
| | - Jia Sheng
- Department of Chemistry, The RNA Institute, University at Albany, State University of New York, Albany, NY, United States
| |
Collapse
|
9
|
Passalacqua LFM, Dingilian AI, Lupták A. Single-pass transcription by T7 RNA polymerase. RNA (NEW YORK, N.Y.) 2020; 26:2062-2071. [PMID: 32958559 PMCID: PMC7668259 DOI: 10.1261/rna.076778.120] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Accepted: 09/04/2020] [Indexed: 06/11/2023]
Abstract
RNA molecules can be conveniently synthesized in vitro by the T7 RNA polymerase (T7 RNAP). In some experiments, such as cotranscriptional biochemical analyses, continuous synthesis of RNA is not desired. Here, we propose a method for a single-pass transcription that yields a single transcript per template DNA molecule using the T7 RNAP system. We hypothesized that stalling the polymerase downstream from the promoter region and subsequent cleavage of the promoter by a restriction enzyme (to prevent promoter binding by another polymerase) would allow synchronized production of a single transcript per template. The single-pass transcription was verified in two different scenarios: a short self-cleaving ribozyme and a long mRNA. The results show that a controlled single-pass transcription using T7 RNAP allows precise measurement of cotranscriptional ribozyme activity, and this approach will facilitate the study of other kinetic events.
Collapse
Affiliation(s)
- Luiz F M Passalacqua
- Department of Pharmaceutical Sciences, University of California, Irvine, California 92697, USA
| | - Armine I Dingilian
- Department of Pharmaceutical Sciences, University of California, Irvine, California 92697, USA
| | - Andrej Lupták
- Department of Pharmaceutical Sciences, University of California, Irvine, California 92697, USA
- Department of Chemistry, University of California, Irvine, California 92697, USA
- Department of Molecular Biology and Biochemistry, University of California, Irvine, California 92697, USA
| |
Collapse
|
10
|
Micura R, Höbartner C. Fundamental studies of functional nucleic acids: aptamers, riboswitches, ribozymes and DNAzymes. Chem Soc Rev 2020; 49:7331-7353. [PMID: 32944725 DOI: 10.1039/d0cs00617c] [Citation(s) in RCA: 102] [Impact Index Per Article: 25.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
This review aims at juxtaposing common versus distinct structural and functional strategies that are applied by aptamers, riboswitches, and ribozymes/DNAzymes. Focusing on recently discovered systems, we begin our analysis with small-molecule binding aptamers, with emphasis on in vitro-selected fluorogenic RNA aptamers and their different modes of ligand binding and fluorescence activation. Fundamental insights are much needed to advance RNA imaging probes for detection of exo- and endogenous RNA and for RNA process tracking. Secondly, we discuss the latest gene expression-regulating mRNA riboswitches that respond to the alarmone ppGpp, to PRPP, to NAD+, to adenosine and cytidine diphosphates, and to precursors of thiamine biosynthesis (HMP-PP), and we outline new subclasses of SAM and tetrahydrofolate-binding RNA regulators. Many riboswitches bind protein enzyme cofactors that, in principle, can catalyse a chemical reaction. For RNA, however, only one system (glmS ribozyme) has been identified in Nature thus far that utilizes a small molecule - glucosamine-6-phosphate - to participate directly in reaction catalysis (phosphodiester cleavage). We wonder why that is the case and what is to be done to reveal such likely existing cellular activities that could be more diverse than currently imagined. Thirdly, this brings us to the four latest small nucleolytic ribozymes termed twister, twister-sister, pistol, and hatchet as well as to in vitro selected DNA and RNA enzymes that promote new chemistry, mainly by exploiting their ability for RNA labelling and nucleoside modification recognition. Enormous progress in understanding the strategies of nucleic acids catalysts has been made by providing thorough structural fundaments (e.g. first structure of a DNAzyme, structures of ribozyme transition state mimics) in combination with functional assays and atomic mutagenesis.
Collapse
Affiliation(s)
- Ronald Micura
- Institute of Organic Chemistry and Center for Molecular Biosciences Innsbruck CMBI, Leopold-Franzens University Innsbruck, Innsbruck, Austria.
| | | |
Collapse
|
11
|
Kalra K, Gorle S, Cavallo L, Oliva R, Chawla M. Occurrence and stability of lone pair-π and OH-π interactions between water and nucleobases in functional RNAs. Nucleic Acids Res 2020; 48:5825-5838. [PMID: 32392301 PMCID: PMC7293021 DOI: 10.1093/nar/gkaa345] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2019] [Revised: 04/21/2020] [Accepted: 05/07/2020] [Indexed: 01/15/2023] Open
Abstract
We identified over 1000 instances of water-nucleobase stacking contacts in a variety of RNA molecules from a non-redundant set of crystal structures with resolution ≤3.0 Å. Such contacts may be of either the lone pair-π (lp-π) or the OH-π type, in nature. The distribution of the distances of the water oxygen from the nucleobase plane peaks at 3.5 Å for A, G and C, and approximately at 3.1-3.2 Å for U. Quantum mechanics (QM) calculations confirm, as expected, that the optimal energy is reached at a shorter distance for the lp-π interaction as compared to the OH-π one (3.0 versus 3.5 Å). The preference of each nucleobase for either type of interaction closely correlates with its electrostatic potential map. Furthermore, QM calculations show that for all the nucleobases a favorable interaction, of either the lp-π or the OH-π type, can be established at virtually any position of the water molecule above the nucleobase skeleton, which is consistent with the uniform projection of the OW atoms over the nucleobases ring we observed in the experimental occurrences. Finally, molecular dynamics simulations of a model system for the characterization of water-nucleobase stacking contacts confirm the stability of these interactions also under dynamic conditions.
Collapse
Affiliation(s)
- Kanav Kalra
- King Abdullah University of Science and Technology (KAUST), Physical Sciences and Engineering Division, Kaust Catalysis Center, Thuwal 23955-6900, Saudi Arabia
| | - Suresh Gorle
- Department of Biochemistry and Molecular Biology, University of Texas Medical Branch, Galveston, TX 77555, USA
| | - Luigi Cavallo
- King Abdullah University of Science and Technology (KAUST), Physical Sciences and Engineering Division, Kaust Catalysis Center, Thuwal 23955-6900, Saudi Arabia
| | - Romina Oliva
- Department of Sciences and Technologies, University Parthenope of Naples, Centro Direzionale Isola C4, I-80143 Naples, Italy
| | - Mohit Chawla
- King Abdullah University of Science and Technology (KAUST), Physical Sciences and Engineering Division, Kaust Catalysis Center, Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
12
|
Ma L, Kartik S, Liu B, Liu J. From general base to general acid catalysis in a sodium-specific DNAzyme by a guanine-to-adenine mutation. Nucleic Acids Res 2019; 47:8154-8162. [PMID: 31276580 PMCID: PMC6736077 DOI: 10.1093/nar/gkz578] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2019] [Revised: 06/20/2019] [Accepted: 06/26/2019] [Indexed: 02/07/2023] Open
Abstract
Recently, a few Na+-specific RNA-cleaving DNAzymes were reported, where nucleobases are likely to play critical roles in catalysis. The NaA43 and NaH1 DNAzymes share the same 16-nt Na+-binding motif, but differ in one or two nucleotides in a small catalytic loop. Nevertheless, they display an opposite pH-dependency, implicating distinct catalytic mechanisms. In this work, rational mutation studies locate a catalytic adenine residue, A22, in NaH1, while previous studies found a guanine (G23) to be important for the catalysis of NaA43. Mutation with pKa-perturbed analogs, such as 2-aminopurine (∼3.8), 2,6-diaminopurine (∼5.6) and hypoxanthine (∼8.7) affected the overall reaction rate. Therefore, we propose that the N1 position of G23 (pKa ∼6.6) in NaA43 functions as a general base, while that of A22 (pKa ∼6.3) in NaH1 as a general acid. Further experiments with base analogs and a phosphorothioate-modified substrate suggest that the exocyclic amine in A22 and both of the non-bridging oxygens at the scissile phosphate are important for catalysis for NaH1. This is an interesting example where single point mutations can change the mechanism of cleavage from general base to general acid, and it can also explain this Na+-dependent DNAzyme scaffold being sensitive to a broad range of metal ions and molecules.
Collapse
Affiliation(s)
- Lingzi Ma
- Department of Chemistry, Waterloo Institute for Nanotechnology, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Sanjana Kartik
- Department of Chemistry, Waterloo Institute for Nanotechnology, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Biwu Liu
- Department of Chemistry, Waterloo Institute for Nanotechnology, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Juewen Liu
- Department of Chemistry, Waterloo Institute for Nanotechnology, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
13
|
Messina KJ, Kierzek R, Tracey MA, Bevilacqua PC. Small Molecule Rescue and Glycosidic Conformational Analysis of the Twister Ribozyme. Biochemistry 2019; 58:4857-4868. [PMID: 31742390 PMCID: PMC6901379 DOI: 10.1021/acs.biochem.9b00742] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The number of self-cleaving ribozymes has increased sharply in recent years, giving rise to elaborations of the four known ribozyme catalytic strategies, α, β, γ, and δ. One such extension is utilized by the twister ribozyme, which is hypothesized to conduct δ, or general acid catalysis, via N3 of the syn adenine +1 nucleobase indirectly via buffer catalysis at biological pH and directly at lower pH. Herein, we test the δ catalysis role of A1 via chemical rescue and the catalytic relevance of the syn orientation of the nucleobase by conformational analysis. Using inhibited twister ribozyme variants with A1(N3) deaza or A1 abasic modifications, we observe >100-fold chemical rescue effects in the presence of protonatable biological small molecules such as imidazole and histidine, similar to observed rescue values previously reported for C75U/C76Δ in the HDV ribozyme. Brønsted plots for the twister variants support a model in which small molecules rescue catalytic activity via a proton transfer mechanism, suggesting that A1 in the wild type is involved in proton transfer, most likely general acid catalysis. Additionally, through glycosidic conformational analysis in an appropriate background that accommodates the bromine atom, we observe that an 8BrA1-modified twister ribozyme is up to 10-fold faster than a nonmodified A1 ribozyme, supporting crystallographic data that show that A1 is syn when conducting proton transfer. Overall, this study provides functional evidence that the nucleotide immediately downstream of the cleavage site participates directly or indirectly in general acid-base catalysis in the twister ribozyme while occupying the syn conformation.
Collapse
Affiliation(s)
- Kyle J. Messina
- Department of Chemistry, Pennsylvania State University, University Park, PA 16802
- Center for RNA Molecular Biology, Pennsylvania State University, University Park, PA 16802
| | - Ryszard Kierzek
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61-704, Poznan, Poland
| | - Matthew A. Tracey
- Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, PA 16802
- Current Address: Department of Biochemistry and Biophysics, University of Pennsylvania Perelman School of Medicine, Philadelphia, PA 19104
| | - Philip C. Bevilacqua
- Department of Chemistry, Pennsylvania State University, University Park, PA 16802
- Center for RNA Molecular Biology, Pennsylvania State University, University Park, PA 16802
- Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, PA 16802
| |
Collapse
|
14
|
Yamagami R, Huang R, Bevilacqua PC. Cellular Concentrations of Nucleotide Diphosphate-Chelated Magnesium Ions Accelerate Catalysis by RNA and DNA Enzymes. Biochemistry 2019; 58:3971-3979. [PMID: 31512860 DOI: 10.1021/acs.biochem.9b00578] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
RNAs are involved in myriad cellular events. In general, RNA function is affected by cellular conditions. For instance, molecular crowding promotes RNA folding through compaction of the RNA. Metabolites generally destabilize RNA secondary structure, which improves RNA folding cooperativity and increases the fraction of functional RNA. Our recent studies demonstrate that cellular concentrations of amino acid-chelated magnesium (aaCM) stimulate RNA folding and catalysis while protecting RNAs from magnesium ion-induced degradation. However, effects of other cellular magnesium ion chelators on RNA function have not been tested. Herein, we report that nucleotide diphosphate-chelated magnesium, which is of intermediate strength, promotes RNA catalysis much like aaCM. Nucleotides are some of the major metabolites in cells and have one to three phosphates, which have increasingly tight binding of magnesium. On the basis of binding calculations, ∼85% ATP, ∼40% ADP, and only 5% AMP are estimated to possess a magnesium ion under cellular conditions of 0.50 mM Mg2+free. We tested the self-cleaving activity of the hammerhead ribozyme in the presence of these chelated magnesium species. Our results indicate that NTP-chelated magnesium and NMP-chelated magnesium do not appreciably stimulate RNA catalysis, whereas NDP-chelated magnesium promotes RNA catalysis up to 6.5-fold. Inspired by NDP, we observed similar stimulatory effects for several other Mg2+ diphosphate-containing metabolites, including UDP-GlcNAC and UDP-Glc; in addition, we found similar effects for a DNAzyme. Thus, rate stimulatory effects are general with respect to the diphosphate and nucleic acid enzyme. These results implicate magnesium-chelated diphosphate metabolites as general facilitators of RNA function inside cells.
Collapse
Affiliation(s)
- Ryota Yamagami
- Department of Chemistry , Pennsylvania State University , University Park , Pennsylvania 16802 , United States.,Center for RNA Molecular Biology , Pennsylvania State University , University Park , Pennsylvania 16802 , United States
| | - Ruochuan Huang
- Department of Chemistry , Pennsylvania State University , University Park , Pennsylvania 16802 , United States
| | - Philip C Bevilacqua
- Department of Chemistry , Pennsylvania State University , University Park , Pennsylvania 16802 , United States.,Center for RNA Molecular Biology , Pennsylvania State University , University Park , Pennsylvania 16802 , United States.,Department of Biochemistry and Molecular Biology , Pennsylvania State University , University Park , Pennsylvania 16802 , United States
| |
Collapse
|
15
|
Yamagami R, Kayedkhordeh M, Mathews DH, Bevilacqua PC. Design of highly active double-pseudoknotted ribozymes: a combined computational and experimental study. Nucleic Acids Res 2019; 47:29-42. [PMID: 30462314 PMCID: PMC6326823 DOI: 10.1093/nar/gky1118] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2018] [Accepted: 10/24/2018] [Indexed: 01/02/2023] Open
Abstract
Design of RNA sequences that adopt functional folds establishes principles of RNA folding and applications in biotechnology. Inverse folding for RNAs, which allows computational design of sequences that adopt specific structures, can be utilized for unveiling RNA functions and developing genetic tools in synthetic biology. Although many algorithms for inverse RNA folding have been developed, the pseudoknot, which plays a key role in folding of ribozymes and riboswitches, is not addressed in most algorithms. For the few algorithms that attempt to predict pseudoknot-containing ribozymes, self-cleavage activity has not been tested. Herein, we design double-pseudoknot HDV ribozymes using an inverse RNA folding algorithm and test their kinetic mechanisms experimentally. More than 90% of the positively designed ribozymes possess self-cleaving activity, whereas more than 70% of negative control ribozymes, which are predicted to fold to the necessary structure but with low fidelity, do not possess it. Kinetic and mutation analyses reveal that these RNAs cleave site-specifically and with the same mechanism as the WT ribozyme. Most ribozymes react just 50- to 80-fold slower than the WT ribozyme, and this rate can be improved to near WT by modification of a junction. Thus, fast-cleaving functional ribozymes with multiple pseudoknots can be designed computationally.
Collapse
Affiliation(s)
- Ryota Yamagami
- Department of Chemistry, Pennsylvania State University, University Park, PA 16802, USA.,Center for RNA Molecular Biology, Pennsylvania State University, University Park, PA 16802, USA
| | - Mohammad Kayedkhordeh
- Department of Biochemistry & Biophysics and Center for RNA Biology, University of Rochester Medical Center, Rochester, New York, NY 14642, USA
| | - David H Mathews
- Department of Biochemistry & Biophysics and Center for RNA Biology, University of Rochester Medical Center, Rochester, New York, NY 14642, USA.,Department of Biostatistics & Computational Biology, University of Rochester Medical Center, Rochester, New York, NY 14642, USA
| | - Philip C Bevilacqua
- Department of Chemistry, Pennsylvania State University, University Park, PA 16802, USA.,Center for RNA Molecular Biology, Pennsylvania State University, University Park, PA 16802, USA.,Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, PA 16802, USA
| |
Collapse
|
16
|
Abstract
Self-cleaving ribozymes are RNAs that catalyze position-specific cleavage of their phosphodiester backbone. The cleavage site of the newly discovered hatchet ribozyme is located at the very 5′ end of its consensus secondary structure motif. Here we report on the 2.1-Å crystal structure of the hatchet ribozyme in the product state, which defines its intricate tertiary fold and identifies key residues lining the catalytic pocket. This in turn has allowed us to propose a model of the precatalytic state structure and a role in catalysis for a conserved guanine. This study therefore provides a structure-based platform toward an improved understanding of the catalytic mechanism of hatchet ribozymes. Small self-cleaving ribozymes catalyze site-specific cleavage of their own phosphodiester backbone with implications for viral genome replication, pre-mRNA processing, and alternative splicing. We report on the 2.1-Å crystal structure of the hatchet ribozyme product, which adopts a compact pseudosymmetric dimeric scaffold, with each monomer stabilized by long-range interactions involving highly conserved nucleotides brought into close proximity of the scissile phosphate. Strikingly, the catalytic pocket contains a cavity capable of accommodating both the modeled scissile phosphate and its flanking 5′ nucleoside. The resulting modeled precatalytic conformation incorporates a splayed-apart alignment at the scissile phosphate, thereby providing structure-based insights into the in-line cleavage mechanism. We identify a guanine lining the catalytic pocket positioned to contribute to cleavage chemistry. The functional relevance of structure-based insights into hatchet ribozyme catalysis is strongly supported by cleavage assays monitoring the impact of selected nucleobase and atom-specific mutations on ribozyme activity.
Collapse
|
17
|
Cellular conditions of weakly chelated magnesium ions strongly promote RNA stability and catalysis. Nat Commun 2018; 9:2149. [PMID: 29858572 PMCID: PMC5984629 DOI: 10.1038/s41467-018-04415-1] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2018] [Accepted: 04/23/2018] [Indexed: 01/06/2023] Open
Abstract
Most RNA folding studies have been performed under non-physiological conditions of high concentrations (≥10 mM) of Mg2+free, while actual cellular concentrations of Mg2+free are only ~1 mM in a background of greater than 50 mM Mg2+total. To uncover cellular behavior of RNA, we devised cytoplasm mimic systems that include biological concentrations of amino acids, which weakly chelate Mg2+. Amino acid-chelated Mg2+ (aaCM) of ~15 mM dramatically increases RNA folding and prevents RNA degradation. Furthermore, aaCM enhance self-cleavage of several different ribozymes, up to 100,000-fold at Mg2+free of just 0.5 mM, indirectly through RNA compaction. Other metabolites that weakly chelate magnesium offer similar beneficial effects, which implies chelated magnesium may enhance RNA function in the cell in the same way. Overall, these results indicate that the states of Mg2+ should not be limited to free and bound only, as weakly bound Mg2+ strongly promotes RNA function under cellular conditions.
Collapse
|
18
|
Lu J, Koo SC, Weissman BP, Harris ME, Li NS, Piccirilli JA. Evidence That Nucleophile Deprotonation Exceeds Bond Formation in the HDV Ribozyme Transition State. Biochemistry 2018; 57:3465-3472. [PMID: 29733591 DOI: 10.1021/acs.biochem.8b00031] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Steric constraints imposed by the active sites of protein and RNA enzymes pose major challenges to the investigation of structure-function relationships within these systems. As a strategy to circumvent such constraints in the HDV ribozyme, we have synthesized phosphoramidites from propanediol derivatives and incorporated them at the 5'-termini of RNA and DNA oligonucleotides to generate a series of novel substrates with nucleophiles perturbed electronically through geminal fluorination. In nonenzymatic, hydroxide-catalyzed intramolecular transphosphorylation of the DNA substrates, pH-rate profiles revealed that fluorine substitution reduces the maximal rate and the kinetic p Ka, consistent with the expected electron-withdrawing effect. In HDV ribozyme reactions, we observed that the RNA substrates undergo transphosphorylation relatively efficiently, suggesting that the conformational constraints imposed by a ribofuranose ring are not strictly required for ribozyme catalysis. In contrast to the nonenzymatic reactions, however, substrate fluorination modestly increases the ribozyme reaction rate, consistent with a mechanism in which (1) the 2'-hydroxyl nucleophile exists predominantly in its neutral, protonated form in the ground state and (2) the 2'-hydroxyl bears some negative charge in the rate-determining step, consistent with a transition state in which the extent of 2'-OH deprotonation exceeds the extent of P-O bond formation.
Collapse
Affiliation(s)
- Jun Lu
- Department of Biochemistry and Molecular Biology and Department of Chemistry , University of Chicago , 929 East 57th Street , Chicago , Illinois 60637 , United States
| | - Selene C Koo
- Department of Biochemistry and Molecular Biology and Department of Chemistry , University of Chicago , 929 East 57th Street , Chicago , Illinois 60637 , United States
| | - Benjamin P Weissman
- Department of Biochemistry and Molecular Biology and Department of Chemistry , University of Chicago , 929 East 57th Street , Chicago , Illinois 60637 , United States
| | - Michael E Harris
- Department of Chemistry , University of Florida , 214 Leigh Hall , Gainesville , Florida 32611 , United States
| | - Nan-Sheng Li
- Department of Biochemistry and Molecular Biology and Department of Chemistry , University of Chicago , 929 East 57th Street , Chicago , Illinois 60637 , United States
| | - Joseph A Piccirilli
- Department of Biochemistry and Molecular Biology and Department of Chemistry , University of Chicago , 929 East 57th Street , Chicago , Illinois 60637 , United States
| |
Collapse
|
19
|
Š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: 336] [Impact Index Per Article: 56.0] [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
|
20
|
Webb CHT, Lupták A. Kinetic Parameters of trans Scission by Extended HDV-like Ribozymes and the Prospect for the Discovery of Genomic trans-Cleaving RNAs. Biochemistry 2018; 57:1440-1450. [PMID: 29388767 DOI: 10.1021/acs.biochem.7b00789] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Hepatitis delta virus (HDV)-like ribozymes are self-cleaving catalytic RNAs with a widespread distribution in nature and biological roles ranging from self-scission during rolling-circle replication in viroids to co-transcriptional processing of eukaryotic retrotransposons, among others. The ribozymes fold into a double pseudoknot with a common catalytic core motif and highly variable peripheral domains. Like other self-cleaving ribozymes, HDV-like ribozymes can be converted into trans-acting catalytic RNAs by bisecting the self-cleaving variants at non-essential loops. Here we explore the trans-cleaving activity of ribozymes derived from the largest examples of the ribozymes (drz-Agam-2 family), which contain an extended domain between the substrate strand and the rest of the RNA. When this peripheral domain is bisected at its distal end, the substrate strand is recognized through two helices, rather than just one 7 bp helix common among the HDV ribozymes, resulting in stronger binding and increased sequence specificity. Kinetic characterization of the extended trans-cleaving ribozyme revealed an efficient trans-cleaving system with a surprisingly high KM', supporting a model that includes a recently proposed activation barrier related to the assembly of the catalytically competent ribozyme. The ribozymes also exhibit a very long koff for the products (∼2 weeks), resulting in a trade-off between sequence specificity and turnover. Finally, structure-based searches for the catalytic cores of these ribozymes in the genome of the mosquito Anopheles gambiae, combined with sequence searches for their putative substrates, revealed two potential ribozyme-substrate pairs that may represent examples of natural trans-cleaving ribozymes.
Collapse
Affiliation(s)
- Chiu-Ho T Webb
- Department of Molecular Biology and Biochemistry , University of California-Irvine , Irvine , California 92697 , United States
| | - Andrej Lupták
- Department of Molecular Biology and Biochemistry , University of California-Irvine , Irvine , California 92697 , United States.,Department of Pharmaceutical Sciences , University of California-Irvine , Irvine , California 92697 , United States.,Department of Chemistry , University of California-Irvine , Irvine , California 92697 , United States
| |
Collapse
|
21
|
Frankel EA, Bevilacqua PC. Complexity in pH-Dependent Ribozyme Kinetics: Dark pK a Shifts and Wavy Rate-pH Profiles. Biochemistry 2017; 57:483-488. [PMID: 29271644 DOI: 10.1021/acs.biochem.7b00784] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Charged bases occur in RNA enzymes, or ribozymes, where they play key roles in catalysis. Cationic bases donate protons and perform electrostatic catalysis, while anionic bases accept protons. We previously published simulations of rate-pH profiles for ribozymes in terms of species plots for the general acid and general base that have been useful for understanding how ribozymes respond to pH. In that study, we did not consider interaction between the general acid and general base or interaction with other species on the RNA. Since that report, diverse small ribozyme classes have been discovered, many of which have charged nucleobases or metal ions in the active site that can either directly interact and participate in catalysis or indirectly interact as "influencers". Herein, we simulate experimental rate-pH profiles in terms of species plots in which reverse protonated charged nucleobases interact. These analyses uncover two surprising features of pH-dependent enzyme kinetics. (1) Cooperativity between the general acid and general base enhances population of the functional forms of a ribozyme and manifests itself as hidden or "dark" pKa shifts, real pKa shifts that accelerate the reaction but are not readily observed by standard experimental approaches, and (2) influencers favorably shift the pKas of proton-transferring nucleobases and manifest themselves as "wavy" rate-pH profiles. We identify parallels with the protein enzyme literature, including reverse protonation and wavelike behavior, while pointing out that RNA is more prone to reverse protonation. The complexities uncovered, which arise from simple pairwise interactions, should aid deconvolution of complex rate-pH profiles for RNA and protein enzymes and suggest veiled catalytic devices for promoting catalysis that can be tested by experiment and calculation.
Collapse
Affiliation(s)
- Erica A Frankel
- Department of Chemistry, Pennsylvania State University , University Park, Pennsylvania 16802, United States.,Center for RNA Molecular Biology, Pennsylvania State University , University Park, Pennsylvania 16802, United States
| | - Philip C Bevilacqua
- Department of Chemistry, Pennsylvania State University , University Park, Pennsylvania 16802, United States.,Center for RNA Molecular Biology, Pennsylvania State University , University Park, Pennsylvania 16802, United States.,Department of Biochemistry and Molecular Biology, Pennsylvania State University , University Park, Pennsylvania 16802, United States
| |
Collapse
|
22
|
Liu X, Chen Y, Fierke CA. Inner-Sphere Coordination of Divalent Metal Ion with Nucleobase in Catalytic RNA. J Am Chem Soc 2017; 139:17457-17463. [PMID: 29116782 PMCID: PMC6020041 DOI: 10.1021/jacs.7b08755] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Identification of the function of metal ions and the RNA moieties, particularly nucleobases, that bind metal ions is important in RNA catalysis. Here we combine single-atom and abasic substitutions to probe functions of conserved nucleobases in ribonuclease P (RNase P). Structural and biophysical studies of bacterial RNase P propose direct coordination of metal ions by the nucleobases of conserved uridine and guanosine in helix P4 of the RNA subunit (P RNA). To biochemically probe the function of metal ion interactions, we substituted the universally conserved bulged uridine (U51) in the P4 helix of circularly permuted Bacillus subtilis P RNA with 4-thiouridine, 4-deoxyuridine, and abasic modifications and G378/379 with 2-aminopurine, N7-deazaguanosine, and 6-thioguanosine. The functional group modifications of U51 decrease RNase P-catalyzed phosphodiester bond cleavage 16- to 23-fold, as measured by the single-turnover cleavage rate constant. The activity of the 4-thiouridine RNase P is partially rescued by addition of Cd(II) or Mn(II) ions. This is the first time a metal-rescue experiment provides evidence for inner-sphere divalent metal ion coordination with a nucleobase. Modifications of G379 modestly decrease the cleavage activity of RNase P, suggesting outer-sphere coordination of O6 on G379 to a metal ion. These data provide biochemical evidence for catalytically important interactions of the P4 helix of P RNA with metal ions, demonstrating that the bulged uridine coordinates at least one catalytic metal ion through an inner-sphere interaction. The combination of single-atom and abasic nucleotide substitutions provides a powerful strategy to probe functions of conserved nucleobases in large RNAs.
Collapse
Affiliation(s)
- Xin Liu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Yu Chen
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Carol A. Fierke
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Biological Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
23
|
Ren A, Micura R, Patel DJ. Structure-based mechanistic insights into catalysis by small self-cleaving ribozymes. Curr Opin Chem Biol 2017; 41:71-83. [PMID: 29107885 PMCID: PMC7955703 DOI: 10.1016/j.cbpa.2017.09.017] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2017] [Revised: 09/22/2017] [Accepted: 09/29/2017] [Indexed: 01/01/2023]
Abstract
Small self-cleaving ribozymes are widely distributed in nature and are essential for rolling-circle-based replication of satellite and pathogenic RNAs. Earlier structure-function studies on the hammerhead, hairpin, glmS, hepatitis delta virus and Varkud satellite ribozymes have provided insights into their overall architecture, their catalytic active site organization, and the role of nearby nucleobases and hydrated divalent cations in facilitating general acid-base and electrostatic-mediated catalysis. This review focuses on recent structure-function research on active site alignments and catalytic mechanisms of the Rzb hammerhead ribozyme, as well as newly-identified pistol, twister and twister-sister ribozymes. In contrast to an agreed upon mechanistic understanding of self-cleavage by Rzb hammerhead and pistol ribozymes, there exists a divergence of views as to the cleavage site alignments and catalytic mechanisms adopted by twister and twister-sister ribozymes. One approach to resolving this conundrum would be to extend the structural studies from currently available pre-catalytic conformations to their transition state mimic vanadate counterparts for both ribozymes.
Collapse
Affiliation(s)
- Aiming Ren
- Life Sciences Institute, Zhejiang University, Hangzhou 310058, China
| | - Ronald Micura
- Institute of Organic Chemistry, Leopold Franzens University, Innsbruck A6020, Austria
| | - Dinshaw J Patel
- Structural Biology Program, Memorial Sloan-Kettering Cancer Center, New York, NY 10065, USA.
| |
Collapse
|
24
|
Mlýnský V, Kührová P, Jurečka P, Šponer J, Otyepka M, Banáš P. Mapping the Chemical Space of the RNA Cleavage and Its Implications for Ribozyme Catalysis. J Phys Chem B 2017; 121:10828-10840. [DOI: 10.1021/acs.jpcb.7b09129] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Vojtěch Mlýnský
- 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
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via
Bonomea 265, 34136 Trieste, Italy
| | - 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
| | - 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
| | - Jiří Šponer
- 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
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 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
| | - 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
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
25
|
Passalacqua LFM, Jimenez RM, Fong JY, Lupták A. Allosteric Modulation of the Faecalibacterium prausnitzii Hepatitis Delta Virus-like Ribozyme by Glucosamine 6-Phosphate: The Substrate of the Adjacent Gene Product. Biochemistry 2017; 56:6006-6014. [PMID: 29045794 DOI: 10.1021/acs.biochem.7b00879] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Self-cleaving ribozymes were discovered 30 years ago and have been found throughout nature, from bacteria to animals, but little is known about their biological functions and regulation, particularly how cofactors and metabolites alter their activity. A hepatitis delta virus-like self-cleaving ribozyme maps upstream of a phosphoglucosamine mutase (glmM) open reading frame in the genome of the human gut bacterium Faecalibacterium prausnitzii. The presence of a ribozyme in the untranslated region of glmM suggests a regulation mechanism of gene expression. In the bacterial hexosamine biosynthesis pathway, the enzyme glmM catalyzes the isomerization of glucosamine 6-phosphate into glucosamine 1-phosphate. In this study, we investigated the effect of these metabolites on the co-transcriptional self-cleavage rate of the ribozyme. Our results suggest that glucosamine 6-phosphate, but not glucosamine 1-phosphate, is an allosteric ligand that increases the self-cleavage rate of drz-Fpra-1, providing the first known example of allosteric modulation of a self-cleaving ribozyme by the substrate of the adjacent gene product. Given that the ribozyme is activated by the glmM substrate, but not the product, this allosteric modulation may represent a potential feed-forward mechanism of gene expression regulation in bacteria.
Collapse
Affiliation(s)
- Luiz F M Passalacqua
- Department of Pharmaceutical Sciences, University of California , Irvine, California 92697, United States
| | - Randi M Jimenez
- Department of Molecular Biology and Biochemistry, University of California , Irvine, California 92697, United States
| | - Jennifer Y Fong
- Department of Pharmaceutical Sciences, University of California , Irvine, California 92697, United States
| | - Andrej Lupták
- Department of Pharmaceutical Sciences, University of California , Irvine, California 92697, United States.,Department of Molecular Biology and Biochemistry, University of California , Irvine, California 92697, United States.,Department of Chemistry, University of California , Irvine, California 92697, United States
| |
Collapse
|
26
|
Gaines CS, York DM. Model for the Functional Active State of the TS Ribozyme from Molecular Simulation. Angew Chem Int Ed Engl 2017. [DOI: 10.1002/ange.201705608] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Colin S. Gaines
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry & Chemical Biology; Rutgers University; 174 Frelinghuysen Road Piscataway NJ 08854-8076 USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry & Chemical Biology; Rutgers University; 174 Frelinghuysen Road Piscataway NJ 08854-8076 USA
| |
Collapse
|
27
|
Gaines CS, York DM. Model for the Functional Active State of the TS Ribozyme from Molecular Simulation. Angew Chem Int Ed Engl 2017; 56:13392-13395. [PMID: 28763583 DOI: 10.1002/anie.201705608] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2017] [Indexed: 12/23/2022]
Abstract
Recently, a crystal structure has been reported of a new catalytic RNA, the TS ribozyme, that has been identified through comparative genomics and is believed to be a metalloribozyme having novel mechanistic features. Although this data provides invaluable structural information, analysis suggests a conformational change is required to arrive at a catalytically relevant state. We report results of molecular simulations that predict a spontaneous local rearrangement of the active site, leading to solution structures consistent with available functional data and providing competing mechanistic hypotheses that can be experimentally tested. The two competing hypotheses differ in the proposed identity of the catalytic general acid: either a water molecule coordinating a Mg2+ ion bound at the Watson-Crick edge of residue C7, or the N3 position of residue C7 itself.
Collapse
Affiliation(s)
- Colin S Gaines
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry & Chemical Biology, Rutgers University, 174 Frelinghuysen Road, Piscataway, NJ, 08854-8076, USA
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry & Chemical Biology, Rutgers University, 174 Frelinghuysen Road, Piscataway, NJ, 08854-8076, USA
| |
Collapse
|
28
|
Liu Y, Wilson TJ, Lilley DM. The structure of a nucleolytic ribozyme that employs a catalytic metal ion. Nat Chem Biol 2017; 13:508-513. [PMID: 28263963 PMCID: PMC5392355 DOI: 10.1038/nchembio.2333] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2016] [Accepted: 12/21/2016] [Indexed: 01/05/2023]
Abstract
The TS ribozyme (originally called "twister sister") is a catalytic RNA. We present a crystal structure of the ribozyme in a pre-reactive conformation. Two co-axial helical stacks are organized by a three-way junction and two tertiary contacts. Five divalent metal ions are directly coordinated to RNA ligands, making important contributions to the RNA architecture. The scissile phosphate lies in a quasihelical loop region that is organized by a network of hydrogen bonding. A divalent metal ion is directly bound to the nucleobase 5' to the scissile phosphate, with an inner-sphere water molecule positioned to interact with the O2' nucleophile. The rate of ribozyme cleavage correlated in a log-linear manner with divalent metal ion pKa, consistent with proton transfer in the transition state, and we propose that the bound metal ion is a likely general base for the cleavage reaction. Our data indicate that the TS ribozyme functions predominantly as a metalloenzyme.
Collapse
Affiliation(s)
- Yijin Liu
- Cancer Research UK Nucleic Acid Structure Research Group, MSI/WTB Complex, The University of Dundee, Dow Street, Dundee DD1 5EH, U.K
| | - Timothy J. Wilson
- Cancer Research UK Nucleic Acid Structure Research Group, MSI/WTB Complex, The University of Dundee, Dow Street, Dundee DD1 5EH, U.K
| | - David M.J. Lilley
- Cancer Research UK Nucleic Acid Structure Research Group, MSI/WTB Complex, The University of Dundee, Dow Street, Dundee DD1 5EH, U.K
| |
Collapse
|
29
|
Structural and Biochemical Properties of Novel Self-Cleaving Ribozymes. Molecules 2017; 22:molecules22040678. [PMID: 28441772 PMCID: PMC6154101 DOI: 10.3390/molecules22040678] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2016] [Revised: 04/07/2017] [Accepted: 04/13/2017] [Indexed: 01/20/2023] Open
Abstract
Fourteen well-defined ribozyme classes have been identified to date, among which nine are site-specific self-cleaving ribozymes. Very recently, small self-cleaving ribozymes have attracted renewed interest in their structure, biochemistry, and biological function since the discovery, during the last three years, of four novel ribozymes, termed twister, twister sister, pistol, and hatchet. In this review, we mainly address the structure, biochemistry, and catalytic mechanism of the novel ribozymes. They are characterized by distinct active site architectures and divergent, but similar, biochemical properties. The cleavage activities of the ribozymes are highly dependent upon divalent cations, pH, and base-specific mutations, which can cause changes in the nucleotide arrangement and/or electrostatic potential around the cleavage site. It is most likely that a guanine and adenine in close proximity of the cleavage site are involved in general acid-base catalysis. In addition, metal ions appear to play a structural rather than catalytic role although some of their crystal structures have shown a direct metal ion coordination to a non-bridging phosphate oxygen at the cleavage site. Collectively, the structural and biochemical data of the four newest ribozymes could contribute to advance our mechanistic understanding of how self-cleaving ribozymes accomplish their efficient site-specific RNA cleavages.
Collapse
|
30
|
Abstract
Twister RNAs represent a recently discovered class of natural ribozymes that promote rapid cleaving of RNA backbones. Although an abundance of theoretical, biochemical, and structural data exist for several members of the twister class, disagreements about the architecture and mechanism of its active site have emerged. Historically, such storms regarding mechanistic details typically occur soon after each new self-cleaving ribozyme class is reported, but paths forward exist to quickly reach calmer conditions.
Collapse
Affiliation(s)
- Ronald R. Breaker
- Department
of Molecular, Cellular and Developmental Biology, ‡Department of Molecular
Biophysics and Biochemistry, §Howard Hughes Medical Institute, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
31
|
Š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
|
32
|
Ucisik MN, Bevilacqua PC, Hammes-Schiffer S. Molecular Dynamics Study of Twister Ribozyme: Role of Mg(2+) Ions and the Hydrogen-Bonding Network in the Active Site. Biochemistry 2016; 55:3834-46. [PMID: 27295275 PMCID: PMC5127262 DOI: 10.1021/acs.biochem.6b00203] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The recently discovered twister ribozyme is thought to utilize general acid-base catalysis in its self-cleavage mechanism, but the roles of nucleobases and metal ions in the mechanism are unclear. Herein, molecular dynamics simulations of the env22 twister ribozyme are performed to elucidate the structural and equilibrium dynamical properties, as well as to examine the role of Mg(2+) ions and possible candidates for the general base and acid in the self-cleavage mechanism. The active site region and the ends of the pseudoknots were found to be less mobile than other regions of the ribozyme, most likely providing structural stability and possibly facilitating catalysis. A purported catalytic Mg(2+) ion and the closest neighboring Mg(2+) ion remained chelated and relatively immobile throughout the microsecond trajectories, although removal of these Mg(2+) ions did not lead to any significant changes in the structure or equilibrium motions of the ribozyme on the microsecond time scale. In addition, a third metal ion, a Na(+) ion remained close to A1(O5'), the leaving group atom, during the majority of the microsecond trajectories, suggesting that it might stabilize the negative charge on A1(O5') during self-cleavage. The locations of these cations and their interactions with key nucleotides in the active site suggest that they may be catalytically relevant. The P1 stem is partially melted at its top and bottom in the crystal structure and further unwinds in the trajectories. The simulations also revealed an interconnected network comprised of hydrogen-bonding and π-stacking interactions that create a relatively rigid network around the self-cleavage site. The nucleotides involved in this network are among the highly conserved nucleotides in twister ribozymes, suggesting that this interaction network may be important to structure and function.
Collapse
Affiliation(s)
- Melek N Ucisik
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801-3364, United States
| | - Philip C Bevilacqua
- Department of Chemistry, Department of Biochemistry & Molecular Biology, and Center for RNA Molecular Biology, The Pennsylvania State University , University Park, Pennsylvania 16802, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801-3364, United States
| |
Collapse
|
33
|
Topological constraints of structural elements in regulation of catalytic activity in HDV-like self-cleaving ribozymes. Sci Rep 2016; 6:28179. [PMID: 27302490 PMCID: PMC4908430 DOI: 10.1038/srep28179] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2016] [Accepted: 05/31/2016] [Indexed: 11/09/2022] Open
Abstract
Self-cleaving ribozymes fold into intricate structures, which orient active site groups into catalytically competent conformations. Most ribozyme families have distinct catalytic cores stabilized by tertiary interactions between domains peripheral to those cores. We show that large hepatitis delta virus (HDV)-like ribozymes are activated by peripheral domains that bring two helical segments, P1 and P2, into proximity – a “pinch” that results in rate acceleration by almost three orders of magnitude. Kinetic analysis of ribozymes with systematically altered length and stability of the peripheral domain revealed that about one third of its free energy of formation is used to lower an activation energy barrier, likely related to a rate-limiting conformational change leading to the pre-catalytic state. These findings provide a quantitative view of enzyme regulation by peripheral domains and may shed light on the energetics of allosteric regulation.
Collapse
|
34
|
D'Ascenzo L, Leonarski F, Vicens Q, Auffinger P. 'Z-DNA like' fragments in RNA: a recurring structural motif with implications for folding, RNA/protein recognition and immune response. Nucleic Acids Res 2016; 44:5944-56. [PMID: 27151194 PMCID: PMC4937326 DOI: 10.1093/nar/gkw388] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2016] [Accepted: 04/28/2016] [Indexed: 01/09/2023] Open
Abstract
Since the work of Alexander Rich, who solved the first Z-DNA crystal structure, we have known that d(CpG) steps can adopt a particular structure that leads to forming left-handed helices. However, it is still largely unrecognized that other sequences can adopt ‘left-handed’ conformations in DNA and RNA, in double as well as single stranded contexts. These ‘Z-like’ steps involve the coexistence of several rare structural features: a C2’-endo puckering, a syn nucleotide and a lone pair–π stacking between a ribose O4’ atom and a nucleobase. This particular arrangement induces a conformational stress in the RNA backbone, which limits the occurrence of Z-like steps to ≈0.1% of all dinucleotide steps in the PDB. Here, we report over 600 instances of Z-like steps, which are located within r(UNCG) tetraloops but also in small and large RNAs including riboswitches, ribozymes and ribosomes. Given their complexity, Z-like steps are probably associated with slow folding kinetics and once formed could lock a fold through the formation of unique long-range contacts. Proteins involved in immunologic response also specifically recognize/induce these peculiar folds. Thus, characterizing the conformational features of these motifs could be a key to understanding the immune response at a structural level.
Collapse
Affiliation(s)
- Luigi D'Ascenzo
- Architecture et Réactivité de l'ARN, Université de Strasbourg, Institut de Biologie Moléculaire et Cellulaire du CNRS, Strasbourg 67084, France
| | - Filip Leonarski
- Architecture et Réactivité de l'ARN, Université de Strasbourg, Institut de Biologie Moléculaire et Cellulaire du CNRS, Strasbourg 67084, France Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Quentin Vicens
- Architecture et Réactivité de l'ARN, Université de Strasbourg, Institut de Biologie Moléculaire et Cellulaire du CNRS, Strasbourg 67084, France
| | - Pascal Auffinger
- Architecture et Réactivité de l'ARN, Université de Strasbourg, Institut de Biologie Moléculaire et Cellulaire du CNRS, Strasbourg 67084, France
| |
Collapse
|
35
|
Skilandat M, Rowinska-Zyrek M, Sigel RKO. Secondary structure confirmation and localization of Mg2+ ions in the mammalian CPEB3 ribozyme. RNA (NEW YORK, N.Y.) 2016; 22:750-763. [PMID: 26966151 PMCID: PMC4836649 DOI: 10.1261/rna.053843.115] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/08/2015] [Accepted: 01/04/2016] [Indexed: 06/05/2023]
Abstract
Most of today's knowledge of the CPEB3 ribozyme, one of the few small self-cleaving ribozymes known to occur in humans, is based on comparative studies with the hepatitis delta virus (HDV) ribozyme, which is highly similar in cleavage mechanism and probably also in structure. Here we present detailed NMR studies of the CPEB3 ribozyme in order to verify the formation of the predicted nested double pseudoknot in solution. In particular, the influence of Mg(2+), the ribozyme's crucial cofactor, on the CPEB3 structure is investigated. NMR titrations, Tb(3+)-induced cleavage, as well as stoichiometry determination by hydroxyquinoline sulfonic acid fluorescence and equilibrium dialysis, are used to evaluate the number, location, and binding mode of Mg(2+)ions. Up to eight Mg(2+)ions interact site-specifically with the ribozyme, four of which are bound with high affinity. The global fold of the CPEB3 ribozyme, encompassing 80%-90% of the predicted base pairs, is formed in the presence of monovalent ions alone. Low millimolar concentrations of Mg(2+)promote a more compact fold and lead to the formation of additional structures in the core of the ribozyme, which contains the inner small pseudoknot and the active site. Several Mg(2+)binding sites, which are important for the functional fold, appear to be located in corresponding locations in the HDV and CPEB3 ribozyme, demonstrating the particular relevance of Mg(2+)for the nested double pseudoknot structure.
Collapse
Affiliation(s)
- Miriam Skilandat
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | | | - Roland K O Sigel
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| |
Collapse
|
36
|
|
37
|
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
|
38
|
Sripathi KN, Banáš P, Réblová K, Šponer J, Otyepka M, Walter NG. Wobble pairs of the HDV ribozyme play specific roles in stabilization of active site dynamics. Phys Chem Chem Phys 2015; 17:5887-900. [PMID: 25631765 DOI: 10.1039/c4cp05083e] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The hepatitis delta virus (HDV) is the only known human pathogen whose genome contains a catalytic RNA motif (ribozyme). The overall architecture of the HDV ribozyme is that of a double-nested pseudoknot, with two GU pairs flanking the active site. Although extensive studies have shown that mutation of either wobble results in decreased catalytic activity, little work has focused on linking these mutations to specific structural effects on catalytic fitness. Here we use molecular dynamics simulations based on an activated structure to probe the active site dynamics as a result of wobble pair mutations. In both wild-type and mutant ribozymes, the in-line fitness of the active site (as a measure of catalytic proficiency) strongly depends on the presence of a C75(N3H3+)N1(O5') hydrogen bond, which positions C75 as the general acid for the reaction. Our mutational analyses show that each GU wobble supports catalytically fit conformations in distinct ways; the reverse G25U20 wobble promotes high in-line fitness, high occupancy of the C75(N3H3+)G1(O5') general-acid hydrogen bond and stabilization of the G1U37 wobble, while the G1U37 wobble acts more locally by stabilizing high in-line fitness and the C75(N3H3+)G1(O5') hydrogen bond. We also find that stable type I A-minor and P1.1 hydrogen bonding above and below the active site, respectively, prevent local structural disorder from spreading and disrupting global conformation. Taken together, our results define specific, often redundant architectural roles for several structural motifs of the HDV ribozyme active site, expanding the known roles of these motifs within all HDV-like ribozymes and other structured RNAs.
Collapse
Affiliation(s)
- Kamali N Sripathi
- Department of Medicinal Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1065, USA
| | | | | | | | | | | |
Collapse
|
39
|
Mir A, Chen J, Robinson K, Lendy E, Goodman J, Neau D, Golden BL. Two Divalent Metal Ions and Conformational Changes Play Roles in the Hammerhead Ribozyme Cleavage Reaction. Biochemistry 2015; 54:6369-81. [PMID: 26398724 PMCID: PMC4710350 DOI: 10.1021/acs.biochem.5b00824] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The hammerhead ribozyme is a self-cleaving RNA broadly dispersed across all kingdoms of life. Although it was the first of the small, nucleolytic ribozymes discovered, the mechanism by which it catalyzes its reaction remains elusive. The nucleobase of G12 is well positioned to be a general base, but it is unclear if or how this guanine base becomes activated for proton transfer. Metal ions have been implicated in the chemical mechanism, but no interactions between divalent metal ions and the cleavage site have been observed crystallographically. To better understand how this ribozyme functions, we have solved crystal structures of wild-type and G12A mutant ribozymes. We observe a pH-dependent conformational change centered around G12, consistent with this nucleotide becoming deprotonated. Crystallographic and kinetic analysis of the G12A mutant reveals a Zn(2+) specificity switch suggesting a direct interaction between a divalent metal ion and the purine at position 12. The metal ion specificity switch and the pH-rate profile of the G12A mutant suggest that the minor imino tautomer of A12 serves as the general base in the mutant ribozyme. We propose a model in which the hammerhead ribozyme rearranges prior to the cleavage reaction, positioning two divalent metal ions in the process. The first metal ion, positioned near G12, becomes directly coordinated to the O6 keto oxygen, to lower the pKa of the general base and organize the active site. The second metal ion, positioned near G10.1, bridges the N7 of G10.1 and the scissile phosphate and may participate directly in the cleavage reaction.
Collapse
Affiliation(s)
- Aamir Mir
- Department of Biochemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Ji Chen
- Department of Biochemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Kyle Robinson
- Department of Biochemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Emma Lendy
- Department of Biochemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Jaclyn Goodman
- Department of Biochemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - David Neau
- Department of Chemistry and Chemical Biology, Cornell University, Northeastern Collaborative Access Team, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Barbara L. Golden
- Department of Biochemistry, Purdue University, West Lafayette, Indiana 47907, United States,Corresponding Author: Telephone: (765) 496-6165;
| |
Collapse
|
40
|
Jimenez RM, Polanco JA, Lupták A. Chemistry and Biology of Self-Cleaving Ribozymes. Trends Biochem Sci 2015; 40:648-661. [PMID: 26481500 DOI: 10.1016/j.tibs.2015.09.001] [Citation(s) in RCA: 127] [Impact Index Per Article: 14.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2015] [Revised: 08/28/2015] [Accepted: 09/01/2015] [Indexed: 11/26/2022]
Abstract
Self-cleaving ribozymes were discovered 30 years ago, but their biological distribution and catalytic mechanisms are only beginning to be defined. Each ribozyme family is defined by a distinct structure, with unique active sites accelerating the same transesterification reaction across the families. Biochemical studies show that general acid-base catalysis is the most common mechanism of self-cleavage, but metal ions and metabolites can be used as cofactors. Ribozymes have been discovered in highly diverse genomic contexts throughout nature, from viroids to vertebrates. Their biological roles include self-scission during rolling-circle replication of RNA genomes, co-transcriptional processing of retrotransposons, and metabolite-dependent gene expression regulation in bacteria. Other examples, including highly conserved mammalian ribozymes, suggest that many new biological roles are yet to be discovered.
Collapse
Affiliation(s)
- Randi M Jimenez
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, USA
| | - Julio A Polanco
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, USA
| | - Andrej Lupták
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, USA; Department of Pharmaceutical Sciences, University of California, Irvine, CA, USA; Department of Chemistry, University of California, Irvine, CA, USA.
| |
Collapse
|
41
|
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
|
42
|
Koo SC, Lu J, Li NS, Leung E, Das SR, Harris ME, Piccirilli JA. Transition State Features in the Hepatitis Delta Virus Ribozyme Reaction Revealed by Atomic Perturbations. J Am Chem Soc 2015; 137:8973-82. [PMID: 26125657 PMCID: PMC4758122 DOI: 10.1021/jacs.5b01189] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Endonucleolytic ribozymes constitute a class of non-coding RNAs that catalyze single-strand RNA scission. With crystal structures available for all of the known ribozymes, a major challenge involves relating functional data to the physically observed RNA architecture. In the case of the hepatitis delta virus (HDV) ribozyme, there are three high-resolution crystal structures, the product state of the reaction and two precursor variants, with distinct mechanistic implications. Here, we develop new strategies to probe the structure and catalytic mechanism of a ribozyme. First, we use double-mutant cycles to distinguish differences in functional group proximity implicated by the crystal structures. Second, we use a corrected form of the Brønsted equation to assess the functional significance of general acid catalysis in the system. Our results delineate the functional relevance of atomic interactions inferred from structure, and suggest that the HDV ribozyme transition state resembles the cleavage product in the degree of proton transfer to the leaving group.
Collapse
Affiliation(s)
- Selene C. Koo
- Departments of Biochemistry & Molecular Biology, and Chemistry, The University of Chicago, Chicago, IL 60637
| | - Jun Lu
- Departments of Biochemistry & Molecular Biology, and Chemistry, The University of Chicago, Chicago, IL 60637
| | - Nan-Sheng Li
- Departments of Biochemistry & Molecular Biology, and Chemistry, The University of Chicago, Chicago, IL 60637
| | - Edward Leung
- Departments of Biochemistry & Molecular Biology, and Chemistry, The University of Chicago, Chicago, IL 60637
| | - Subha R. Das
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, PA 15213
| | - Michael E. Harris
- Department of Biochemistry and Center for Proteomics, Case Western Reserve University School of Medicine, Cleveland, OH 44106
| | - Joseph A. Piccirilli
- Departments of Biochemistry & Molecular Biology, and Chemistry, The University of Chicago, Chicago, IL 60637
| |
Collapse
|
43
|
Huang M, York DM. Linear free energy relationships in RNA transesterification: theoretical models to aid experimental interpretations. Phys Chem Chem Phys 2015; 16:15846-55. [PMID: 24961771 DOI: 10.1039/c4cp01050g] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
RNA cleavage transesterification is of fundamental reaction in biology that is catalyzed by both protein and RNA enzymes. In this work, a series of RNA transesterification model reactions with a wide range of leaving groups are investigated with density-functional calculations in an aqueous solvation environment in order to study linear free energy relationships (LFERs) and their connection to transition state structure and bonding. Overall, results obtained from the polarizable continuum solvation model with UAKS radii produce the best linear correlations and closest overall agreement with experimental results. Reactions with a poor leaving group are predicted to proceed via a stepwise mechanism with a late transition state that is rate controlling. As leaving group becomes more acidic and labile, the barriers of both early and late transition states decrease. LFERs for each transition state are computed, with the late transition state barrier showing greater sensitivity to leaving group pKa. For sufficiently enhanced leaving groups, the reaction mechanism transits to a concerted mechanism characterized by a single early transition state. Further linear relationships were derived for bond lengths and bond orders as a function of leaving group pKa and rate constant values that can be used for prediction. This work provides important benchmark linear free energy data that allows a molecular-level characterization of the structure and bonding of the transition states for this important class of phosphoryl transfer reactions. The relations reported herein can be used to aid in the interpretation of data obtained from experimental studies of non-catalytic and catalytic mechanisms.
Collapse
Affiliation(s)
- Ming Huang
- Scientific Computation, University of Minnesota, 207 Pleasant St. SE, Minneapolis, MN 55455-0431, USA
| | | |
Collapse
|
44
|
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
|
45
|
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
|
46
|
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
|
47
|
Mlýnský V, Walter NG, Šponer J, Otyepka M, Banáš P. The role of an active site Mg(2+) in HDV ribozyme self-cleavage: insights from QM/MM calculations. Phys Chem Chem Phys 2015; 17:670-9. [PMID: 25412464 PMCID: PMC4256098 DOI: 10.1039/c4cp03857f] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The hepatitis delta virus (HDV) ribozyme is a catalytic RNA motif embedded in the human pathogenic HDV RNA. It catalyzes self-cleavage of its sugar-phosphate backbone with direct participation of the active site cytosine C75. Biochemical and structural data support a general acid role of C75. Here, we used hybrid quantum mechanical/molecular mechanical (QM/MM) calculations to probe the reaction mechanism and changes in Gibbs energy along the ribozyme's reaction pathway with an N3-protonated C75H(+) in the active site, which acts as the general acid, and a partially hydrated Mg(2+) ion with one deprotonated, inner-shell coordinated water molecule that acts as the general base. We followed eight reaction paths with a distinct position and coordination of the catalytically important active site Mg(2+) ion. For six of them, we observed feasible activation barriers ranging from 14.2 to 21.9 kcal mol(-1), indicating that the specific position of the Mg(2+) ion in the active site is predicted to strongly affect the kinetics of self-cleavage. The deprotonation of the U-1(2'-OH) nucleophile and the nucleophilic attack of the resulting U-1(2'-O(-)) on the scissile phosphodiester are found to be separate steps, as deprotonation precedes the nucleophilic attack. This sequential mechanism of the HDV ribozyme differs from the concerted nucleophilic activation and attack suggested for the hairpin ribozyme. We estimate the pKa of the U-1(2'-OH) group to range from 8.8 to 11.2, suggesting that it is lowered by several units from that of a free ribose, comparable to and most likely smaller than the pKa of the solvated active site Mg(2+) ion. Our results thus support the notion that the structure of the HDV ribozyme, and particularly the positioning of the active site Mg(2+) ion, facilitate deprotonation and activation of the 2'-OH nucleophile.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic.
| | | | | | | | | |
Collapse
|
48
|
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
|
49
|
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
|
50
|
Crystal structure and mechanistic investigation of the twister ribozyme. Nat Chem Biol 2014; 10:739-44. [DOI: 10.1038/nchembio.1587] [Citation(s) in RCA: 105] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2014] [Accepted: 06/16/2014] [Indexed: 02/07/2023]
|