1
|
Luo B, Zhang C, Ling X, Mukherjee S, Jia G, Xie J, Jia X, Liu L, Baulin EF, Luo Y, Jiang L, Dong H, Wei X, Bujnicki JM, Su Z. Cryo-EM reveals dynamics of Tetrahymena group I intron self-splicing. Nat Catal 2023. [DOI: 10.1038/s41929-023-00934-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2023]
|
2
|
Rees HC, Gogacz W, Li NS, Koirala D, Piccirilli JA. Structural Basis for Fluorescence Activation by Pepper RNA. ACS Chem Biol 2022; 17:1866-1875. [PMID: 35759696 PMCID: PMC9969808 DOI: 10.1021/acschembio.2c00290] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Pepper is a fluorogenic RNA aptamer tag that binds to a variety of benzylidene-cyanophenyl (HBC) derivatives with tight affinity and activates their fluorescence. To investigate how Pepper RNA folds to create a binding site for HBC, we used antibody-assisted crystallography to determine the structures of Pepper bound to HBC530 and HBC599 to 2.3 and 2.7 Å resolutions, respectively. The structural data show that Pepper folds into an elongated structure and organizes nucleotides within an internal bulge to create the ligand binding site, assisted by an out-of-plane platform created by tertiary interactions with an adjacent bulge. As predicted from a lack of K+ dependence, Pepper does not use a G-quadruplex to form a binding pocket for HBC. Instead, Pepper uses a unique base-quadruple·base-triple stack to sandwich the ligand with a U·G wobble pair. Site-bound Mg2+ ions support ligand binding structurally and energetically. This research provides insight into the structural features that allow the Pepper aptamer to bind HBC and show how Pepper's function may expand to allow the in vivo detection of other small molecules and metals.
Collapse
Affiliation(s)
- Huw C. Rees
- Department of Chemistry, University of Chicago, Chicago, Illinois, 60637, United States
| | - Wojciech Gogacz
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois, 60637, United States
| | - Nan-Sheng Li
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois, 60637, United States
| | - Deepak Koirala
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois, 60637, United States
| | - Joseph A. Piccirilli
- Department of Chemistry, University of Chicago, Chicago, Illinois, 60637, United States,Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois, 60637, United States,corresponding author
| |
Collapse
|
3
|
Cryo-EM structures of full-length Tetrahymena ribozyme at 3.1 Å resolution. Nature 2021; 596:603-607. [PMID: 34381213 PMCID: PMC8405103 DOI: 10.1038/s41586-021-03803-w] [Citation(s) in RCA: 55] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Accepted: 07/06/2021] [Indexed: 02/07/2023]
Abstract
Single-particle cryogenic electron microscopy (cryo-EM) has become a standard technique for determining protein structures at atomic resolution1-3. However, cryo-EM studies of protein-free RNA are in their early days. The Tetrahymena thermophila group I self-splicing intron was the first ribozyme to be discovered and has been a prominent model system for the study of RNA catalysis and structure-function relationships4, but its full structure remains unknown. Here we report cryo-EM structures of the full-length Tetrahymena ribozyme in substrate-free and bound states at a resolution of 3.1 Å. Newly resolved peripheral regions form two coaxially stacked helices; these are interconnected by two kissing loop pseudoknots that wrap around the catalytic core and include two previously unforeseen (to our knowledge) tertiary interactions. The global architecture is nearly identical in both states; only the internal guide sequence and guanosine binding site undergo a large conformational change and a localized shift, respectively, upon binding of RNA substrates. These results provide a long-sought structural view of a paradigmatic RNA enzyme and signal a new era for the cryo-EM-based study of structure-function relationships in ribozymes.
Collapse
|
4
|
Grotz KK, Cruz-León S, Schwierz N. Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties. J Chem Theory Comput 2021; 17:2530-2540. [PMID: 33720710 PMCID: PMC8047801 DOI: 10.1021/acs.jctc.0c01281] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2020] [Indexed: 12/31/2022]
Abstract
Magnesium ions play an essential role in many vital processes. To correctly describe their interactions in molecular dynamics simulations, an accurate parametrization is crucial. Despite the importance and considerable scientific effort, current force fields based on the commonly used 12-6 Lennard-Jones interaction potential fail to reproduce a variety of experimental solution properties. In particular, no parametrization exists so far that simultaneously reproduces the solvation free energy and the distance to the water oxygens in the first hydration shell. Moreover, current Mg2+ force fields significantly underestimate the rate of water exchange leading to unrealistically slow exchange kinetics. In order to make progress in the development of improved models, we systematically optimize the Mg2+ parameters in combination with the TIP3P water model in a much larger parameter space than previously done. The results show that a long-ranged interaction potential and modified Lorentz-Berthelot combination rules allow us to accurately reproduce multiple experimental properties including the solvation free energy, the distances to the oxygens of the first hydration shell, the hydration number, the activity coefficient derivative in MgCl2 solutions, the self-diffusion coefficient, and the binding affinity to the phosphate oxygen of RNA. Matching this broad range of thermodynamic properties, we present two sets of optimal parameters: MicroMg yields water exchange on the microsecond timescale in agreement with experiments. NanoMg yields water exchange on the nanosecond timescale facilitating the direct observation of ion-binding events. As shown for the example of the add A-riboswitch, the optimized parameters correctly reproduce the structure of specifically bound ions and permit the de novo prediction of Mg2+-binding sites in biomolecular simulations.
Collapse
Affiliation(s)
- Kara K. Grotz
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| | - Sergio Cruz-León
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| |
Collapse
|
5
|
Abdullin D, Schiemann O. Localization of metal ions in biomolecules by means of pulsed dipolar EPR spectroscopy. Dalton Trans 2021; 50:808-815. [PMID: 33416053 DOI: 10.1039/d0dt03596c] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Metal ions are important for the folding, structure, and function of biomolecules. Thus, knowing where their binding sites are located in proteins or oligonucleotides is a critical objective. X-ray crystallography and nuclear magnetic resonance are powerful methods in this respect, but both have their limitations. Here, a complementary method is highlighted in which paramagnetic metal ions are localized by means of trilateration using a combination of site-directed spin labeling and pulsed dipolar electron paramagnetic resonance spectroscopy. The working principle, the requirements, and the limitations of the method are critically discussed. Several applications of the method are outlined and compared with each other.
Collapse
Affiliation(s)
- Dinar Abdullin
- Institute of Physical and Theoretical Chemistry, University of Bonn, Wegelerstr. 12, 53115 Bonn, Germany.
| | | |
Collapse
|
6
|
Pallan PS, Lybrand TP, Schlegel MK, Harp JM, Jahns H, Manoharan M, Egli M. Incorporating a Thiophosphate Modification into a Common RNA Tetraloop Motif Causes an Unanticipated Stability Boost. Biochemistry 2020; 59:4627-4637. [PMID: 33275419 DOI: 10.1021/acs.biochem.0c00685] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
GNRA (N = A, C, G, or U; R = A or G) tetraloops are common RNA secondary structural motifs and feature a phosphate stacked atop a nucleobase. The rRNA sarcin/ricin loop (SRL) is capped by GApGA, and the phosphate p stacks on G. We recently found that regiospecific incorporation of a single dithiophosphate (PS2) but not a monothiophosphate (PSO) instead of phosphate in the backbone of RNA aptamers dramatically increases the binding affinity for their targets. In the RNA:thrombin complex, the key contribution to the 1000-fold tighter binding stems from an edge-on contact between PS2 and a phenylalanine ring. Here we investigated the consequences of replacing the SRL phosphate engaged in a face-on interaction with guanine with either PS2 or PSO for stability. We found that PS2···G and Rp-PSO···G contacts stabilize modified SRLs compared to the parent loop to unexpected levels: up to 6.3 °C in melting temperature Tm and -4.7 kcal/mol in ΔΔG°. Crystal structures demonstrate that the vertical distance to guanine for the closest sulfur is just 0.05 Å longer on average compared to that of oxygen despite the larger van der Waals radius of the former (1.80 Å for S vs 1.52 Å for O). The higher stability is enthalpy-based, and the negative charge as assessed by a neutral methylphosphonate modification plays only a minor role. Quantum mechanical/molecular mechanical calculations are supportive of favorable dispersion attraction interactions by sulfur making the dominant contribution. A stacking interaction between phosphate and guanine (SRL) or uracil (U-turn) is also found in newly classified RNA tetraloop families besides GNRA.
Collapse
Affiliation(s)
| | | | - Mark K Schlegel
- Alnylam Pharmaceuticals, 300 Third Street, Cambridge, Massachusetts 02142, United States
| | | | - Hartmut Jahns
- Alnylam Pharmaceuticals, 300 Third Street, Cambridge, Massachusetts 02142, United States
| | - Muthiah Manoharan
- Alnylam Pharmaceuticals, 300 Third Street, Cambridge, Massachusetts 02142, United States
| | | |
Collapse
|
7
|
Perlinska AP, Kalek M, Christian T, Hou YM, Sulkowska JI. Mg 2+-Dependent Methyl Transfer by a Knotted Protein: A Molecular Dynamics Simulation and Quantum Mechanics Study. ACS Catal 2020; 10:8058-8068. [PMID: 32904895 PMCID: PMC7462349 DOI: 10.1021/acscatal.0c00059] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2020] [Revised: 06/18/2020] [Indexed: 11/27/2022]
Abstract
![]()
Mg2+ is required for the catalytic activity of TrmD,
a bacteria-specific methyltransferase that is made up of a protein
topological knot-fold, to synthesize methylated m1G37-tRNA
to support life. However, neither the location of Mg2+ in
the structure of TrmD nor its role in the catalytic mechanism is known.
Using molecular dynamics (MD) simulations, we identify a plausible
Mg2+ binding pocket within the active site of the enzyme,
wherein the ion is coordinated by two aspartates and a glutamate.
In this position, Mg2+ additionally interacts with the
carboxylate of a methyl donor cofactor S-adenosylmethionine (SAM).
The computational results are validated by experimental mutation studies,
which demonstrate the importance of the Mg2+-binding residues
for the catalytic activity. The presence of Mg2+ in the
binding pocket induces SAM to adopt a unique bent shape required for
the methyl transfer activity and causes a structural reorganization
of the active site. Quantum mechanical calculations show that the
methyl transfer is energetically feasible only when Mg2+ is bound in the position revealed by the MD simulations, demonstrating
that its function is to align the active site residues within the
topological knot-fold in a geometry optimal for catalysis. The obtained
insights provide the opportunity for developing a strategy of antibacterial
drug discovery based on targeting of Mg2+-binding to TrmD.
Collapse
Affiliation(s)
- Agata P. Perlinska
- Centre of New Technologies, University of Warsaw, Warsaw 02-097, Poland
- College of Inter-Faculty Individual Studies in Mathematics and Natural Sciences, University of Warsaw, Warsaw 02-097, Poland
| | - Marcin Kalek
- Centre of New Technologies, University of Warsaw, Warsaw 02-097, Poland
| | - Thomas Christian
- Department of Biochemistry and Molecular Biology, Thomas Jefferson University, Philadelphia, Pennsylvania 19107, United States
| | - Ya-Ming Hou
- Department of Biochemistry and Molecular Biology, Thomas Jefferson University, Philadelphia, Pennsylvania 19107, United States
| | - Joanna I. Sulkowska
- Centre of New Technologies, University of Warsaw, Warsaw 02-097, Poland
- Faculty of Chemistry, University of Warsaw, Warsaw 02-097, Poland
| |
Collapse
|
8
|
Zou A, Lee S, Li J, Zhou R. Retained Stability of the RNA Structure in DNA Packaging Motor with a Single Mg2+ Ion Bound at the Double Mg-Clamp Structure. J Phys Chem B 2020; 124:701-707. [DOI: 10.1021/acs.jpcb.9b06428] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Aodong Zou
- Department of Engineering Mechanics, Zhejiang University, Hangzhou 310027, China
- Institute of Quantitative Biology and Department of Physics, Zhejiang University, Hangzhou 310027, China
| | - Sangyun Lee
- Computational Biological Center, IBM Thomas J. Watson Research Center, Yorktown Heights, New York 10598, United States
| | - Jingyuan Li
- Institute of Quantitative Biology and Department of Physics, Zhejiang University, Hangzhou 310027, China
| | - Ruhong Zhou
- Institute of Quantitative Biology and Department of Physics, Zhejiang University, Hangzhou 310027, China
- Computational Biological Center, IBM Thomas J. Watson Research Center, Yorktown Heights, New York 10598, United States
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
9
|
Gracia B, Al-Hashimi HM, Bisaria N, Das R, Herschlag D, Russell R. Hidden Structural Modules in a Cooperative RNA Folding Transition. Cell Rep 2019; 22:3240-3250. [PMID: 29562180 PMCID: PMC5894117 DOI: 10.1016/j.celrep.2018.02.101] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2017] [Revised: 02/05/2018] [Accepted: 02/26/2018] [Indexed: 02/01/2023] Open
Abstract
Large-scale, cooperative rearrangements underlie the functions of RNA in RNA-protein machines and gene regulation. To understand how such rearrangements are orchestrated, we used high-throughput chemical footprinting to dissect a seemingly concerted rearrangement in P5abc RNA, a paradigm of RNA folding studies. With mutations that systematically disrupt or restore putative structural elements, we found that this transition reflects local folding of structural modules, with modest and incremental cooperativity that results in concerted behavior. First, two distant secondary structure changes are coupled through a bridging three-way junction and Mg2+-dependent tertiary structure. Second, long-range contacts are formed between modules, resulting in additional cooperativity. Given the sparseness of RNA tertiary contacts after secondary structure formation, we expect that modular folding and incremental cooperativity are generally important for specifying functional structures while also providing productive kinetic paths to these structures. Additionally, we expect our approach to be useful for uncovering modularity in other complex RNAs.
Collapse
Affiliation(s)
- Brant Gracia
- Department of Molecular Biosciences and the Institute for Cellular and Molecular Biology, University of Texas at Austin, Austin, TX 78712, USA
| | - Hashim M Al-Hashimi
- Department of Biochemistry, Duke University Medical Center, Durham, NC 27710, USA; Department of Chemistry, Duke University Medical Center, Durham, NC 27710, USA
| | - Namita Bisaria
- Department of Biochemistry, Stanford University, Stanford, CA 94305, USA
| | - Rhiju Das
- Department of Biochemistry, Stanford University, Stanford, CA 94305, USA
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, CA 94305, USA
| | - Rick Russell
- Department of Molecular Biosciences and the Institute for Cellular and Molecular Biology, University of Texas at Austin, Austin, TX 78712, USA.
| |
Collapse
|
10
|
Yesselman JD, Eiler D, Carlson ED, Gotrik MR, d'Aquino AE, Ooms AN, Kladwang W, Carlson PD, Shi X, Costantino DA, Herschlag D, Lucks JB, Jewett MC, Kieft JS, Das R. Computational design of three-dimensional RNA structure and function. NATURE NANOTECHNOLOGY 2019; 14:866-873. [PMID: 31427748 PMCID: PMC7324284 DOI: 10.1038/s41565-019-0517-8] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/26/2017] [Accepted: 06/24/2019] [Indexed: 05/30/2023]
Abstract
RNA nanotechnology seeks to create nanoscale machines by repurposing natural RNA modules. The field is slowed by the current need for human intuition during three-dimensional structural design. Here, we demonstrate that three distinct problems in RNA nanotechnology can be reduced to a pathfinding problem and automatically solved through an algorithm called RNAMake. First, RNAMake discovers highly stable single-chain solutions to the classic problem of aligning a tetraloop and its sequence-distal receptor, with experimental validation from chemical mapping, gel electrophoresis, solution X-ray scattering and crystallography with 2.55 Å resolution. Second, RNAMake automatically generates structured tethers that integrate 16S and 23S ribosomal RNAs into single-chain ribosomal RNAs that remain uncleaved by ribonucleases and assemble onto messenger RNA. Third, RNAMake enables the automated stabilization of small-molecule binding RNAs, with designed tertiary contacts that improve the binding affinity of the ATP aptamer and improve the fluorescence and stability of the Spinach RNA in cell extracts and in living Escherichia coli cells.
Collapse
Affiliation(s)
- Joseph D Yesselman
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA
| | - Daniel Eiler
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, CO, USA
| | - Erik D Carlson
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, IL, USA
- Chemistry of Life Processes Institute, Northwestern University, Evanston, IL, USA
- Center for Synthetic Biology, Northwestern University, Evanston, IL, USA
| | - Michael R Gotrik
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA
| | - Anne E d'Aquino
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, IL, USA
- Center for Synthetic Biology, Northwestern University, Evanston, IL, USA
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, IL, USA
| | - Alexandra N Ooms
- Department of Cancer Genetics & Genomics, Stanford University School of Medicine, Stanford, CA, USA
| | - Wipapat Kladwang
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA
| | - Paul D Carlson
- Robert F. Smith School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, NY, USA
| | - Xuesong Shi
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA
| | - David A Costantino
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, CO, USA
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA
- Department of Chemistry, Stanford University School of Medicine, Stanford, CA, USA
- Stanford ChEM-H (Chemistry, Engineering, and Medicine for Human Health), Stanford University, Stanford, CA, USA
| | - Julius B Lucks
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, IL, USA
- Chemistry of Life Processes Institute, Northwestern University, Evanston, IL, USA
- Center for Synthetic Biology, Northwestern University, Evanston, IL, USA
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, IL, USA
| | - Michael C Jewett
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, IL, USA
- Chemistry of Life Processes Institute, Northwestern University, Evanston, IL, USA
- Center for Synthetic Biology, Northwestern University, Evanston, IL, USA
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, IL, USA
| | - Jeffrey S Kieft
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, CO, USA
| | - Rhiju Das
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA.
- Department of Physics, Stanford University, Stanford, CA, USA.
| |
Collapse
|
11
|
Bevilacqua PC, Harris ME, Piccirilli JA, Gaines C, Ganguly A, Kostenbader K, Ekesan Ş, York DM. An Ontology for Facilitating Discussion of Catalytic Strategies of RNA-Cleaving Enzymes. ACS Chem Biol 2019; 14:1068-1076. [PMID: 31095369 PMCID: PMC6661149 DOI: 10.1021/acschembio.9b00202] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
A predictive understanding of the mechanisms of RNA cleavage is important for the design of emerging technology built from biological and synthetic molecules that have promise for new biochemical and medicinal applications. Over the past 15 years, RNA cleavage reactions involving 2'-O-transphosphorylation have been discussed using a simplified framework introduced by Breaker that consists of four fundamental catalytic strategies (designated α, β, γ, and δ) that contribute to rate enhancement. As more detailed mechanistic data emerge, there is need for the framework to evolve and keep pace. We develop an ontology for discussion of strategies of enzymes that catalyze RNA cleavage via 2'-O-transphosphorylation that stratifies Breaker's framework into primary (1°), secondary (2°), and tertiary (3°) contributions to enable more precise interpretation of mechanism in the context of structure and bonding. Further, we point out instances where atomic-level changes give rise to changes in more than one catalytic contribution, a phenomenon we refer to as "functional blurring". We hope that this ontology will help clarify our conversations and pave the path forward toward a consensus view of these fundamental and fascinating mechanisms. The insight gained will deepen our understanding of RNA cleavage reactions catalyzed by natural protein and RNA enzymes, as well as aid in the design of new engineered DNA and synthetic enzymes.
Collapse
Affiliation(s)
- Philip C. Bevilacqua
- Department of Chemistry, Center for RNA Molecular Biology, and Department of Biochemistry, Microbiology, and Molecular Biology, The Pennsylvania State University, University Park, PA 16802, USA
| | - Michael E. Harris
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, USA
| | - Joseph A. Piccirilli
- Department of Chemistry, and Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, IL 60637, USA
| | - Colin Gaines
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087, USA
| | - Abir Ganguly
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087, USA
| | - Ken Kostenbader
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087, USA
| | - Şölen Ekesan
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087, USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087, USA
| |
Collapse
|
12
|
Abstract
The past decades have witnessed tremendous developments in our understanding of RNA biology. At the core of these advances have been studies aimed at discerning RNA structure and at understanding the forces that influence the RNA folding process. It is easy to take the present state of understanding for granted, but there is much to be learned by considering the path to our current understanding, which has been tortuous, with the birth and death of models, the adaptation of experimental tools originally developed for characterization of protein structure and catalysis, and the development of novel tools for probing RNA. In this review we tour the stages of RNA folding studies, considering them as "epochs" that can be generalized across scientific disciplines. These epochs span from the discovery of catalytic RNA, through biophysical insights into the putative primordial RNA World, to characterization of structured RNAs, the building and testing of models, and, finally, to the development of models with the potential to yield generalizable predictive and quantitative models for RNA conformational, thermodynamic, and kinetic behavior. We hope that this accounting will aid others as they navigate the many fascinating questions about RNA and its roles in biology, in the past, present, and future.
Collapse
Affiliation(s)
- Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, California 94305
- Department of Chemical Engineering, Stanford University, Stanford, California 94305
- Department of Chemistry, Stanford University, Stanford, California 94305
- Stanford ChEM-H (Chemistry, Engineering, and Medicine for Human Health), Stanford, California 94305
| | - Steve Bonilla
- Department of Biochemistry, Stanford University, Stanford, California 94305
- Department of Chemical Engineering, Stanford University, Stanford, California 94305
| | - Namita Bisaria
- Whitehead Institute for Biomedical Research, Cambridge, Massachusetts 02142
| |
Collapse
|
13
|
Woods CT, Laederach A. Classification of RNA structure change by 'gazing' at experimental data. Bioinformatics 2018; 33:1647-1655. [PMID: 28130241 PMCID: PMC5447233 DOI: 10.1093/bioinformatics/btx041] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2016] [Accepted: 01/20/2017] [Indexed: 11/12/2022] Open
Abstract
Motivation Mutations (or Single Nucleotide Variants) in folded RiboNucleic Acid structures that cause local or global conformational change are riboSNitches. Predicting riboSNitches is challenging, as it requires making two, albeit related, structure predictions. The data most often used to experimentally validate riboSNitch predictions is Selective 2' Hydroxyl Acylation by Primer Extension, or SHAPE. Experimentally establishing a riboSNitch requires the quantitative comparison of two SHAPE traces: wild-type (WT) and mutant. Historically, SHAPE data was collected on electropherograms and change in structure was evaluated by 'gel gazing.' SHAPE data is now routinely collected with next generation sequencing and/or capillary sequencers. We aim to establish a classifier capable of simulating human 'gazing' by identifying features of the SHAPE profile that human experts agree 'looks' like a riboSNitch. Results We find strong quantitative agreement between experts when RNA scientists 'gaze' at SHAPE data and identify riboSNitches. We identify dynamic time warping and seven other features predictive of the human consensus. The classSNitch classifier reported here accurately reproduces human consensus for 167 mutant/WT comparisons with an Area Under the Curve (AUC) above 0.8. When we analyze 2019 mutant traces for 17 different RNAs, we find that features of the WT SHAPE reactivity allow us to improve thermodynamic structure predictions of riboSNitches. This is significant, as accurate RNA structural analysis and prediction is likely to become an important aspect of precision medicine. Availability and Implementation The classSNitch R package is freely available at http://classsnitch.r-forge.r-project.org . Contact alain@email.unc.edu. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Chanin Tolson Woods
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA.,Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Alain Laederach
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA.,Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| |
Collapse
|
14
|
Bonilla S, Limouse C, Bisaria N, Gebala M, Mabuchi H, Herschlag D. Single-Molecule Fluorescence Reveals Commonalities and Distinctions among Natural and in Vitro-Selected RNA Tertiary Motifs in a Multistep Folding Pathway. J Am Chem Soc 2017; 139:18576-18589. [PMID: 29185740 PMCID: PMC5748328 DOI: 10.1021/jacs.7b08870] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
![]()
Decades
of study of the RNA folding problem have revealed that
diverse and complex structured RNAs are built from a common set of
recurring structural motifs, leading to the perspective that a generalizable
model of RNA folding may be developed from understanding of the folding
properties of individual structural motifs. We used single-molecule
fluorescence to dissect the kinetic and thermodynamic properties of
a set of variants of a common tertiary structural motif, the tetraloop/tetraloop-receptor
(TL/TLR). Our results revealed a multistep TL/TLR folding pathway
in which preorganization of the ubiquitous AA-platform submotif precedes
the formation of the docking transition state and tertiary A-minor
hydrogen bond interactions form after the docking transition state.
Differences in ion dependences between TL/TLR variants indicated the
occurrence of sequence-dependent conformational rearrangements prior
to and after the formation of the docking transition state. Nevertheless,
varying the junction connecting the TL/TLR produced a common kinetic
and ionic effect for all variants, suggesting that the global conformational
search and compaction electrostatics are energetically independent
from the formation of the tertiary motif contacts. We also found that in vitro-selected variants, despite their similar stability
at high Mg2+ concentrations, are considerably less stable
than natural variants under near-physiological ionic conditions, and
the occurrence of the TL/TLR sequence variants in Nature correlates
with their thermodynamic stability in isolation. Overall, our findings
are consistent with modular but complex energetic properties of RNA
structural motifs and will aid in the eventual quantitative description
of RNA folding from its secondary and tertiary structural elements.
Collapse
Affiliation(s)
- Steve Bonilla
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| | - Charles Limouse
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| | - Namita Bisaria
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| | - Magdalena Gebala
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| | - Hideo Mabuchi
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| | - Daniel Herschlag
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| |
Collapse
|
15
|
Yeldell SB, Ruble BK, Dmochowski IJ. Oligonucleotide modifications enhance probe stability for single cell transcriptome in vivo analysis (TIVA). Org Biomol Chem 2017; 15:10001-10009. [PMID: 29052679 PMCID: PMC5718921 DOI: 10.1039/c7ob02353g] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Single cell transcriptomics provides a powerful discovery tool for identifying new cell types and functions as well as a means to probe molecular features of the etiology and treatment of human diseases, including cancer. However, such analyses are limited by the difficulty of isolating mRNA from single cells within biological samples. We recently introduced a photochemical method for isolating mRNA from single living cells, Transcriptome In Vivo Analysis (TIVA). The TIVA probe is a "caged" polyU : polyA oligonucleotide hairpin designed to enter live tissue, where site-specific activation with 405 nm laser reveals the polyU-biotin strand to bind mRNA in a target cell, enabling subsequent mRNA isolation and sequencing. The TIVA method is well suited for analysis of living cells in resected tissue, but has not yet been applied to living cells in whole organisms. Adapting TIVA to this more challenging environment requires a probe with higher thermal stability, more robust caging, and greater nuclease resistance. In this paper we present modifications to the original TIVA probe with multiple aspects of enhanced stability. These newer probes utilize an extended 22mer polyU capture strand with two 9mer polyA blocking strands ("22/9/9") for higher thermal stability pre-photolysis and improved mRNA capture affinity post-photolysis. The "22/9/9 GC" probe features a terminal GC pair to reduce pre-photolysis interactions with mRNA by more than half. The "PS-22/9/9" probe features a phosphorothioated backbone, which extends serum stability from <1 h to at least 48 h, and also mediates uptake into cultured human fibroblasts.
Collapse
Affiliation(s)
- S B Yeldell
- Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, PA 19104-6323, USA.
| | | | | |
Collapse
|
16
|
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
|
17
|
Bingaman JL, Zhang S, Stevens DR, Yennawar NH, Hammes-Schiffer S, Bevilacqua PC. The GlcN6P cofactor plays multiple catalytic roles in the glmS ribozyme. Nat Chem Biol 2017; 13:439-445. [PMID: 28192411 PMCID: PMC5362308 DOI: 10.1038/nchembio.2300] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2016] [Accepted: 12/09/2016] [Indexed: 01/06/2023]
Abstract
RNA enzymes (ribozymes) have remarkably diverse biological roles despite having limited chemical diversity. Protein enzymes enhance their reactivity through recruitment of cofactors; likewise, the naturally occurring glmS ribozyme uses the glucosamine-6-phosphate (GlcN6P) organic cofactor for phosphodiester bond cleavage. Prior structural and biochemical studies have implicated GlcN6P as the general acid. Here we describe new catalytic roles of GlcN6P through experiments and calculations. Large stereospecific normal thio effects and a lack of metal-ion rescue in the holoribozyme indicate that nucleobases and the cofactor play direct chemical roles and align the active site for self-cleavage. Large stereospecific inverse thio effects in the aporibozyme suggest that the GlcN6P cofactor disrupts an inhibitory interaction of the nucleophile. Strong metal-ion rescue in the aporibozyme reveals that this cofactor also provides electrostatic stabilization. Ribozyme organic cofactors thus perform myriad catalytic roles, thereby allowing RNA to compensate for its limited functional diversity.
Collapse
Affiliation(s)
- Jamie L. Bingaman
- Department of Chemistry and Center for RNA Molecular
Biology, The Pennsylvania State University, University Park, Pennsylvania 16802,
United States
| | - Sixue Zhang
- Department of Chemistry, University of Illinois at
Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United
States
| | - David R. Stevens
- Department of Chemistry, University of Illinois at
Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United
States
| | - Neela H. Yennawar
- X-ray Crystallography Facility, Huck Institutes of the Life
Sciences, The Pennsylvania State University, 8 Althouse Laboratory, 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, United
States
| | - Philip C. Bevilacqua
- Department of Chemistry and Center for RNA Molecular
Biology, The Pennsylvania State University, University Park, Pennsylvania 16802,
United States
- Department of Biochemistry and Molecular Biology, The
Pennsylvania State University, University Park, Pennsylvania 16802, United
States
| |
Collapse
|
18
|
Leonarski F, D'Ascenzo L, Auffinger P. Mg2+ ions: do they bind to nucleobase nitrogens? Nucleic Acids Res 2017; 45:987-1004. [PMID: 27923930 PMCID: PMC5314772 DOI: 10.1093/nar/gkw1175] [Citation(s) in RCA: 64] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2016] [Revised: 11/10/2016] [Accepted: 11/14/2016] [Indexed: 01/28/2023] Open
Abstract
Given the many roles proposed for Mg2+ in nucleic acids, it is essential to accurately determine their binding modes. Here, we surveyed the PDB to classify Mg2+ inner-sphere binding patterns to nucleobase imine N1/N3/N7 atoms. Among those, purine N7 atoms are considered to be the best nucleobase binding sites for divalent metals. Further, Mg2+ coordination to N7 has been implied in several ribozyme catalytic mechanisms. We report that Mg2+ assigned near imine nitrogens derive mostly from poor interpretations of electron density patterns and are most often misidentified Na+, K+, NH4+ ions, water molecules or spurious density peaks. Consequently, apart from few documented exceptions, Mg2+ ions do not bind to N7 atoms. Without much of a surprise, Mn2+, Zn2+ and Cd2+, which have a higher affinity for nitrogens, may contact N7 atoms when present in crystallization buffers. In this respect, we describe for the first time a potential Zn2+ ribosomal binding site involving two purine N7 atoms. Further, we provide a set of guidelines to help in the assignment of Mg2+ in crystallographic, cryo-EM, NMR and model building practices and discuss implications of our findings related to ion substitution experiments.
Collapse
Affiliation(s)
- Filip Leonarski
- Université de Strasbourg, CNRS, Architecture et Réactivité de l'ARN, UPR9002, F-67000 Strasbourg, France
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Luigi D'Ascenzo
- Université de Strasbourg, CNRS, Architecture et Réactivité de l'ARN, UPR9002, F-67000 Strasbourg, France
| | - Pascal Auffinger
- Université de Strasbourg, CNRS, Architecture et Réactivité de l'ARN, UPR9002, F-67000 Strasbourg, France
| |
Collapse
|
19
|
Lenz TK, Norris AM, Hud NV, Williams LD. Protein-free ribosomal RNA folds to a near-native state in the presence of Mg2+. RSC Adv 2017. [DOI: 10.1039/c7ra08696b] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The assembled bacterial ribosome contains around 50 proteins and many counterions.
Collapse
Affiliation(s)
- Timothy K. Lenz
- Department of Chemistry & Biochemistry
- Georgia Institute of Technology
- Atlanta
- USA
| | - Ashlyn M. Norris
- Department of Chemistry & Biochemistry
- Georgia Institute of Technology
- Atlanta
- USA
| | - Nicholas V. Hud
- Department of Chemistry & Biochemistry
- Georgia Institute of Technology
- Atlanta
- USA
| | - Loren Dean Williams
- Department of Chemistry & Biochemistry
- Georgia Institute of Technology
- Atlanta
- USA
| |
Collapse
|
20
|
Gracia B, Xue Y, Bisaria N, Herschlag D, Al-Hashimi HM, Russell R. RNA Structural Modules Control the Rate and Pathway of RNA Folding and Assembly. J Mol Biol 2016; 428:3972-3985. [PMID: 27452365 PMCID: PMC5048535 DOI: 10.1016/j.jmb.2016.07.013] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2016] [Revised: 07/12/2016] [Accepted: 07/14/2016] [Indexed: 11/27/2022]
Abstract
Structured RNAs fold through multiple pathways, but we have little understanding of the molecular features that dictate folding pathways and determine rates along a given pathway. Here, we asked whether folding of a complex RNA can be understood from its structural modules. In a two-piece version of the Tetrahymena group I ribozyme, the separated P5abc subdomain folds to local native secondary and tertiary structure in a linked transition and assembles with the ribozyme core via three tertiary contacts: a kissing loop (P14), a metal core-receptor interaction, and a tetraloop-receptor interaction, the first two of which are expected to depend on native P5abc structure from the local transition. Native gel, NMR, and chemical footprinting experiments showed that mutations that destabilize the native P5abc structure slowed assembly up to 100-fold, indicating that P5abc folds first and then assembles with the core by conformational selection. However, rate decreases beyond 100-fold were not observed because an alternative pathway becomes dominant, with nonnative P5abc binding the core and then undergoing an induced-fit rearrangement. P14 is formed in the rate-limiting step along the conformational selection pathway but after the rate-limiting step along the induced-fit pathway. Strikingly, the assembly rate along the conformational selection pathway resembles that of an isolated kissing loop similar to P14, and the rate along the induced-fit pathway resembles that of an isolated tetraloop-receptor interaction. Our results indicate substantial modularity in RNA folding and assembly and suggest that these processes can be understood in terms of underlying structural modules.
Collapse
Affiliation(s)
- Brant Gracia
- Department of Molecular Biosciences and the Institute for Cellular and Molecular Biology, University of Texas at Austin, Austin, TX 78712, USA
| | - Yi Xue
- Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
| | - Namita Bisaria
- Department of Biochemistry, Stanford University, Stanford, CA 94305, USA
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, CA 94305, USA
| | - Hashim M Al-Hashimi
- Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
| | - Rick Russell
- Department of Molecular Biosciences and the Institute for Cellular and Molecular Biology, University of Texas at Austin, Austin, TX 78712, USA.
| |
Collapse
|
21
|
Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway. Proc Natl Acad Sci U S A 2016; 113:E4956-65. [PMID: 27493222 DOI: 10.1073/pnas.1525082113] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
The past decade has seen a wealth of 3D structural information about complex structured RNAs and identification of functional intermediates. Nevertheless, developing a complete and predictive understanding of the folding and function of these RNAs in biology will require connection of individual rate and equilibrium constants to structural changes that occur in individual folding steps and further relating these steps to the properties and behavior of isolated, simplified systems. To accomplish these goals we used the considerable structural knowledge of the folded, unfolded, and intermediate states of P4-P6 RNA. We enumerated structural states and possible folding transitions and determined rate and equilibrium constants for the transitions between these states using single-molecule FRET with a series of mutant P4-P6 variants. Comparisons with simplified constructs containing an isolated tertiary contact suggest that a given tertiary interaction has a stereotyped rate for breaking that may help identify structural transitions within complex RNAs and simplify the prediction of folding kinetics and thermodynamics for structured RNAs from their parts. The preferred folding pathway involves initial formation of the proximal tertiary contact. However, this preference was only ∼10 fold and could be reversed by a single point mutation, indicating that a model akin to a protein-folding contact order model will not suffice to describe RNA folding. Instead, our results suggest a strong analogy with a modified RNA diffusion-collision model in which tertiary elements within preformed secondary structures collide, with the success of these collisions dependent on whether the tertiary elements are in their rare binding-competent conformations.
Collapse
|
22
|
Gulen B, Petrov AS, Okafor CD, Vander Wood D, O'Neill EB, Hud NV, Williams LD. Ribosomal small subunit domains radiate from a central core. Sci Rep 2016; 6:20885. [PMID: 26876483 PMCID: PMC4753503 DOI: 10.1038/srep20885] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2015] [Accepted: 01/05/2016] [Indexed: 12/26/2022] Open
Abstract
The domain architecture of a large RNA can help explain and/or predict folding, function, biogenesis and evolution. We offer a formal and general definition of an RNA domain and use that definition to experimentally characterize the rRNA of the ribosomal small subunit. Here the rRNA comprising a domain is compact, with a self-contained system of molecular interactions. A given rRNA helix or stem-loop must be allocated uniquely to a single domain. Local changes such as mutations can give domain-wide effects. Helices within a domain have interdependent orientations, stabilities and interactions. With these criteria we identify a core domain (domain A) of small subunit rRNA. Domain A acts as a hub, linking the four peripheral domains and imposing orientational and positional restraints on the other domains. Experimental characterization of isolated domain A, and mutations and truncations of it, by methods including selective 2'OH acylation analyzed by primer extension and circular dichroism spectroscopy are consistent with our architectural model. The results support the utility of the concept of an RNA domain. Domain A, which exhibits structural similarity to tRNA, appears to be an essential core of the small ribosomal subunit.
Collapse
Affiliation(s)
- Burak Gulen
- School of Chemistry and Biochemistry, Georgia institute of Technology, Atlanta, Georgia 30332, United States of America
| | - Anton S Petrov
- School of Chemistry and Biochemistry, Georgia institute of Technology, Atlanta, Georgia 30332, United States of America
| | - C Denise Okafor
- School of Chemistry and Biochemistry, Georgia institute of Technology, Atlanta, Georgia 30332, United States of America
| | - Drew Vander Wood
- School of Chemistry and Biochemistry, Georgia institute of Technology, Atlanta, Georgia 30332, United States of America
| | - Eric B O'Neill
- School of Chemistry and Biochemistry, Georgia institute of Technology, Atlanta, Georgia 30332, United States of America
| | - Nicholas V Hud
- School of Chemistry and Biochemistry, Georgia institute of Technology, Atlanta, Georgia 30332, United States of America
| | - Loren Dean Williams
- School of Chemistry and Biochemistry, Georgia institute of Technology, Atlanta, Georgia 30332, United States of America
| |
Collapse
|
23
|
|
24
|
Probing the kinetic and thermodynamic consequences of the tetraloop/tetraloop receptor monovalent ion-binding site in P4-P6 RNA by smFRET. Biochem Soc Trans 2016; 43:172-8. [PMID: 25849913 DOI: 10.1042/bst20140268] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Structured RNA molecules play roles in central biological processes and understanding the basic forces and features that govern RNA folding kinetics and thermodynamics can help elucidate principles that underlie biological function. Here we investigate one such feature, the specific interaction of monovalent cations with a structured RNA, the P4-P6 domain of the Tetrahymena ribozyme. We employ single molecule FRET (smFRET) approaches as these allow determination of folding equilibrium and rate constants over a wide range of stabilities and thus allow direct comparisons without the need for extrapolation. These experiments provide additional evidence for specific binding of monovalent cations, Na+ and K+, to the RNA tetraloop-tetraloop receptor (TL-TLR) tertiary motif. These ions facilitate both folding and unfolding, consistent with an ability to help order the TLR for binding and further stabilize the tertiary contact subsequent to attainment of the folding transition state.
Collapse
|
25
|
Abstract
Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.
Collapse
Affiliation(s)
- Pascal Auffinger
- Architecture et Réactivité de l'ARN, Université de Strasbourg, IBMC, CNRS, 15 rue René Descartes, F-67084, Strasbourg, France.
| | - Luigi D'Ascenzo
- Architecture et Réactivité de l'ARN, Université de Strasbourg, IBMC, CNRS, 15 rue René Descartes, F-67084, Strasbourg, France.
| | - Eric Ennifar
- Architecture et Réactivité de l'ARN, Université de Strasbourg, IBMC, CNRS, 15 rue René Descartes, F-67084, Strasbourg, France.
| |
Collapse
|
26
|
Panteva MT, Giambaşu GM, York DM. Force Field for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) Ions That Have Balanced Interactions with Nucleic Acids. J Phys Chem B 2015; 119:15460-70. [PMID: 26583536 DOI: 10.1021/acs.jpcb.5b10423] [Citation(s) in RCA: 88] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Divalent metal ions are of fundamental importance to the function and folding of nucleic acids. Divalent metal ion-nucleic acid interactions are complex in nature and include both territorial and site specific binding. Commonly employed nonbonded divalent ion models, however, are often parametrized against bulk ion properties and are subsequently utilized in biomolecular simulations without considering any data related to interactions at specific nucleic acid sites. Previously, we assessed the ability of 17 different nonbonded Mg(2+) ion models to reproduce different properties of Mg(2+) in aqueous solution including radial distribution functions, solvation free energies, water exchange rates, and translational diffusion coefficients. In the present work, we depart from the recently developed 12-6-4 potential models for divalent metal ions developed by Li and Merz and tune the pairwise parameters for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) binding dimethyl phosphate, adenosine, and guanosine in order to reproduce experimental site specific binding free energies derived from potentiometric pH titration data. We further apply these parameters to investigate a metal ion migration previously proposed to occur during the catalytic reaction of the hammerhead ribozyme. The new parameters are shown to be accurate and balanced for nucleic acid binding in comparison with available experimental data and provide an important tool for molecular dynamics and free energy simulations of nucleic acids where these ions may exhibit different binding modes.
Collapse
Affiliation(s)
- Maria T Panteva
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University , 174 Frelinghuysen Road, Piscataway, New Jersey 08854-8076, United States
| | - George M Giambaşu
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University , 174 Frelinghuysen Road, Piscataway, New Jersey 08854-8076, United States
| | - Darrin M York
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University , 174 Frelinghuysen Road, Piscataway, New Jersey 08854-8076, United States
| |
Collapse
|
27
|
Panteva MT, GiambaȈsu GM, York DM. Comparison of structural, thermodynamic, kinetic and mass transport properties of Mg(2+) ion models commonly used in biomolecular simulations. J Comput Chem 2015; 36:970-82. [PMID: 25736394 PMCID: PMC4409555 DOI: 10.1002/jcc.23881] [Citation(s) in RCA: 66] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2014] [Revised: 01/17/2015] [Accepted: 02/08/2015] [Indexed: 01/09/2023]
Abstract
The prevalence of Mg(2+) ions in biology and their essential role in nucleic acid structure and function has motivated the development of various Mg(2+) ion models for use in molecular simulations. Currently, the most widely used models in biomolecular simulations represent a nonbonded metal ion as an ion-centered point charge surrounded by a nonelectrostatic pairwise potential that takes into account dispersion interactions and exchange effects that give rise to the ion's excluded volume. One strategy toward developing improved models for biomolecular simulations is to first identify a Mg(2+) model that is consistent with the simulation force fields that closely reproduces a range of properties in aqueous solution, and then, in a second step, balance the ion-water and ion-solute interactions by tuning parameters in a pairwise fashion where necessary. The present work addresses the first step in which we compare 17 different nonbonded single-site Mg(2+) ion models with respect to their ability to simultaneously reproduce structural, thermodynamic, kinetic and mass transport properties in aqueous solution. None of the models based on a 12-6 nonelectrostatic nonbonded potential was able to reproduce the experimental radial distribution function, solvation free energy, exchange barrier and diffusion constant. The models based on a 12-6-4 potential offered improvement, and one model in particular, in conjunction with the SPC/E water model, performed exceptionally well for all properties. The results reported here establish useful benchmark calculations for Mg(2+) ion models that provide insight into the origin of the behavior in aqueous solution, and may aid in the development of next-generation models that target specific binding sites in biomolecules.
Collapse
Affiliation(s)
- Maria T. Panteva
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University, 174 Frelinghuysen Road, Piscataway, NJ 08854-8076, USA
| | - George M. GiambaȈsu
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University, 174 Frelinghuysen Road, Piscataway, NJ 08854-8076, USA
| | - Darrin M. York
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University, 174 Frelinghuysen Road, Piscataway, NJ 08854-8076, USA
| |
Collapse
|
28
|
Eckstein F. Phosphorothioates, Essential Components of Therapeutic Oligonucleotides. Nucleic Acid Ther 2014; 24:374-87. [DOI: 10.1089/nat.2014.0506] [Citation(s) in RCA: 335] [Impact Index Per Article: 33.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Affiliation(s)
- Fritz Eckstein
- Max-Planck-Institut für Experimentelle Medizin, Göttingen, Germany
| |
Collapse
|
29
|
A divalent metal ion-dependent N(1)-methyl transfer to G37-tRNA. ACTA ACUST UNITED AC 2014; 21:1351-1360. [PMID: 25219964 DOI: 10.1016/j.chembiol.2014.07.023] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2014] [Revised: 07/20/2014] [Accepted: 07/22/2014] [Indexed: 01/09/2023]
Abstract
The catalytic mechanism of the majority of S-adenosyl methionine (AdoMet)-dependent methyl transferases requires no divalent metal ions. Here we report that methyl transfer from AdoMet to N(1) of G37-tRNA, catalyzed by the bacterial TrmD enzyme, is strongly dependent on divalent metal ions and that Mg(2+) is the most physiologically relevant. Kinetic isotope analysis, metal rescue, and spectroscopic measurements indicate that Mg(2+) is not involved in substrate binding, but in promoting methyl transfer. On the basis of the pH-activity profile indicating one proton transfer during the TrmD reaction, we propose a catalytic mechanism in which the role of Mg(2+) is to help to increase the nucleophilicity of N(1) of G37 and stabilize the negative developing charge on O(6) during attack on the methyl sulfonium of AdoMet. This work demonstrates how Mg(2+) contributes to the catalysis of AdoMet-dependent methyl transfer in one of the most crucial posttranscriptional modifications to tRNA.
Collapse
|
30
|
Abstract
Ions surround nucleic acids in what is referred to as an ion atmosphere. As a result, the folding and dynamics of RNA and DNA and their complexes with proteins and with each other cannot be understood without a reasonably sophisticated appreciation of these ions' electrostatic interactions. However, the underlying behavior of the ion atmosphere follows physical rules that are distinct from the rules of site binding that biochemists are most familiar and comfortable with. The main goal of this review is to familiarize nucleic acid experimentalists with the physical concepts that underlie nucleic acid-ion interactions. Throughout, we provide practical strategies for interpreting and analyzing nucleic acid experiments that avoid pitfalls from oversimplified or incorrect models. We briefly review the status of theories that predict or simulate nucleic acid-ion interactions and experiments that test these theories. Finally, we describe opportunities for going beyond phenomenological fits to a next-generation, truly predictive understanding of nucleic acid-ion interactions.
Collapse
Affiliation(s)
- Jan Lipfert
- Department of Bionanoscience, Kavli Institute of Nanoscience, Delft University of Technology, 2628 CJ Delft, Netherlands;
| | | | | | | |
Collapse
|
31
|
Tseng CK, Cheng SC. The spliceosome catalyzes debranching in competition with reverse of the first chemical reaction. RNA (NEW YORK, N.Y.) 2013; 19:971-81. [PMID: 23681507 PMCID: PMC3683931 DOI: 10.1261/rna.038638.113] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
Splicing of nuclear pre-mRNA occurs via two steps of the transesterification reaction, forming a lariat intermediate and product. The reactions are catalyzed by the spliceosome, a large ribonucleoprotein complex composed of five small nuclear RNAs and numerous protein factors. The spliceosome shares a similar catalytic core structure with that of fungal group II introns, which can self-splice using the same chemical mechanism. Like group II introns, both catalytic steps of pre-mRNA splicing can efficiently reverse on the affinity-purified spliceosome. The spliceosome also catalyzes a hydrolytic spliced-exon reopening reaction as observed in group II introns, indicating a strong link in their evolutionary relationship. We show here that, by arresting splicing after the first catalytic step, the purified spliceosome can catalyze debranching of lariat-intron-exon 2. The debranching reaction, although not observed in group II introns, has similar monovalent cation preferences as those for splicing catalysis of group II introns. The debranching reaction is in competition with the reverse Step 1 reaction influenced by the ionic environment and the structure of components binding near the catalytic center, suggesting that the catalytic center of the spliceosome can switch between different conformations to direct different chemical reactions.
Collapse
|
32
|
Li NS, Frederiksen JK, Piccirilli JA. Automated solid-phase synthesis of RNA oligonucleotides containing a nonbridging phosphorodithioate linkage via phosphorothioamidites. J Org Chem 2012; 77:9889-92. [PMID: 23050987 DOI: 10.1021/jo301834p] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
This work describes a general method for the synthesis of oligoribonucleotides containing a site-specific nonbridging phosphorodithioate linkage via automated solid-phase synthesis using 5'-O-DMTr-2'-O-TBS-ribonucleoside 3'-N,N-dimethyl-S-(2,4-dichlorobenzyl) phosphorothioamidites (2a-2d). The 3'-phosphorothioamidites (2a-2d) can be conveniently prepared in good yields (86-99%) via a one-pot reaction from the corresponding 5'-O-DMTr-2'-O-TBS-ribonucleosides (1a-1d).
Collapse
Affiliation(s)
- Nan-Sheng Li
- Department of Biochemistry & Molecular Biology, The University of Chicago, 929 East 57th Street, Chicago, Illinois 60637, United States.
| | | | | |
Collapse
|