1
|
Hou Q, Chatterjee S, Lund PE, Suddala KC, Walter NG. Single-molecule FRET observes opposing effects of urea and TMAO on structurally similar meso- and thermophilic riboswitch RNAs. Nucleic Acids Res 2023; 51:11345-11357. [PMID: 37855661 PMCID: PMC10639078 DOI: 10.1093/nar/gkad866] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2023] [Revised: 09/17/2023] [Accepted: 09/27/2023] [Indexed: 10/20/2023] Open
Abstract
Bacteria live in a broad range of environmental temperatures that require adaptations of their RNA sequences to maintain function. Riboswitches are regulatory RNAs that change conformation upon typically binding metabolite ligands to control bacterial gene expression. The paradigmatic small class-I preQ1 riboswitches from the mesophile Bacillus subtilis (Bsu) and the thermophile Thermoanaerobacter tengcongensis (Tte) adopt similar pseudoknot structures when bound to preQ1. Here, we use UV-melting analysis combined with single-molecule detected chemical denaturation by urea to compare the thermodynamic and kinetic folding properties of the two riboswitches, and the urea-countering effects of trimethylamine N-oxide (TMAO). Our results show that, first, the Tte riboswitch is more thermotolerant than the Bsu riboswitch, despite only subtle sequence differences. Second, using single-molecule FRET, we find that urea destabilizes the folded pseudoknot structure of both riboswitches, yet has a lower impact on the unfolding kinetics of the thermodynamically less stable Bsu riboswitch. Third, our analysis shows that TMAO counteracts urea denaturation and promotes folding of both the riboswitches, albeit with a smaller effect on the more stable Tte riboswitch. Together, these findings elucidate how subtle sequence adaptations in a thermophilic bacterium can stabilize a common RNA structure when a new ecological niche is conquered.
Collapse
Affiliation(s)
- Qian Hou
- Single Molecule Analysis Group, Department of Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
- Tri-Institutional PhD Program in Chemical Biology, Weill Cornell Medicine, The Rockefeller University, Memorial Sloan Kettering Cancer Center, NY, NY 10021, USA
| | - Surajit Chatterjee
- Single Molecule Analysis Group, Department of Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
- Department of Physics, Case Western Reserve University, Cleveland, Ohio 44106, USA
| | - Paul E Lund
- Single Molecule Analysis Group, Department of Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
| | - Krishna C Suddala
- Single Molecule Analysis Group, Department of Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
- Laboratory of Molecular Biology, National Institute of Diabetes and Digestive and Kidney Diseases, NIH, Bethesda, MD 20892, USA
| | - Nils G Walter
- Single Molecule Analysis Group, Department of Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
| |
Collapse
|
2
|
Martinez-Monge A, Pastor I, Bustamante C, Manosas M, Ritort F. Measurement of the specific and non-specific binding energies of Mg 2+ to RNA. Biophys J 2022; 121:3010-3022. [PMID: 35864738 PMCID: PMC9463699 DOI: 10.1016/j.bpj.2022.07.020] [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: 11/02/2022] Open
Abstract
Determining the non-specific and specific electrostatic contributions of magnesium binding to RNA is a challenging problem. We introduce a single-molecule method based on measuring the folding energy of a native RNA in magnesium and at its equivalent sodium concentration. The latter is defined so that the folding energy in sodium equals the non-specific electrostatic contribution in magnesium. The sodium equivalent can be estimated according to the empirical 100/1 rule (1 M NaCl is equivalent to 10 mM MgCl2), which is a good approximation for most RNAs. The method is applied to an RNA three-way junction (3WJ) that contains specific Mg2+ binding sites and misfolds into a double hairpin structure without binding sites. We mechanically pull the RNA with optical tweezers and use fluctuation theorems to determine the folding energies of the native and misfolded structures in magnesium (10 mM MgCl2) and at the equivalent sodium condition (1 M NaCl). While the free energies of the misfolded structure are equal in magnesium and sodium, they are not for the native structure, the difference being due to the specific binding energy of magnesium to the 3WJ, which equals ΔG≃ 10 kcal/mol. Besides stabilizing the 3WJ, Mg2+ also kinetically rescues it from the misfolded structure over timescales of tens of seconds in a force-dependent manner. The method should generally be applicable to determine the specific binding energies of divalent cations to other tertiary RNAs.
Collapse
Affiliation(s)
- A Martinez-Monge
- Small Biosystems Lab, Departament de Física de la Matèria Condensada, Facultat de Física, Universitat de Barcelona, Carrer de Martí i Franquès, 1, 08028 Barcelona, Spain
| | - Isabel Pastor
- Small Biosystems Lab, Departament de Física de la Matèria Condensada, Facultat de Física, Universitat de Barcelona, Carrer de Martí i Franquès, 1, 08028 Barcelona, Spain; Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, 08028 Barcelona, Spain
| | - Carlos Bustamante
- Departments of Chemistry, Physics and Molecular and Cell Biology, University of California Berkeley, Berkeley, California; Howard Hughes Medical Institute University of California Berkeley, Berkeley, California; Kavli Energy Nanosciences Institute, University of California Berkeley, Berkeley, California
| | - Maria Manosas
- Small Biosystems Lab, Departament de Física de la Matèria Condensada, Facultat de Física, Universitat de Barcelona, Carrer de Martí i Franquès, 1, 08028 Barcelona, Spain; Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, 08028 Barcelona, Spain.
| | - Felix Ritort
- Small Biosystems Lab, Departament de Física de la Matèria Condensada, Facultat de Física, Universitat de Barcelona, Carrer de Martí i Franquès, 1, 08028 Barcelona, Spain; Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, 08028 Barcelona, Spain.
| |
Collapse
|
3
|
Henning-Knechtel A, Thirumalai D, Kirmizialtin S. Differences in ion-RNA binding modes due to charge density variations explain the stability of RNA in monovalent salts. SCIENCE ADVANCES 2022; 8:eabo1190. [PMID: 35857829 PMCID: PMC9299541 DOI: 10.1126/sciadv.abo1190] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
The stability of RNA increases as the charge density of the alkali metal cations increases. The molecular mechanism for this phenomenon remains elusive. To fill this gap, we performed all-atom molecular dynamics pulling simulations of HIV-1 trans-activation response RNA. We first established that the free energy landscape obtained in the simulations is in excellent agreement with the single-molecule optical tweezer experiments. The origin of the stronger stability in sodium compared to potassium is found to be due to the differences in the charge density-related binding modes. The smaller hydrated sodium ion preferentially binds to the highly charged phosphates that have high surface area. In contrast, the larger potassium ions interact with the major grooves. As a result, more cations condense around phosphate groups in the case of sodium ions, leading to the reduction of electrostatic repulsion. Because the proposed mechanism is generic, we predict that the same conclusions are valid for divalent alkaline earth metal cations.
Collapse
Affiliation(s)
| | - D. Thirumalai
- Department of Chemistry, University of Texas, Austin, TX 78712, USA
- Corresponding author. (D.T.); (S.K.)
| | - Serdal Kirmizialtin
- Chemistry Program, Math and Sciences, New York University Abu Dhabi, Abu Dhabi, UAE
- Corresponding author. (D.T.); (S.K.)
| |
Collapse
|
4
|
Progress toward SHAPE Constrained Computational Prediction of Tertiary Interactions in RNA Structure. Noncoding RNA 2021; 7:ncrna7040071. [PMID: 34842779 PMCID: PMC8628965 DOI: 10.3390/ncrna7040071] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Revised: 10/29/2021] [Accepted: 11/02/2021] [Indexed: 01/04/2023] Open
Abstract
As more sequencing data accumulate and novel puzzling genetic regulations are discovered, the need for accurate automated modeling of RNA structure increases. RNA structure modeling from chemical probing experiments has made tremendous progress, however accurately predicting large RNA structures is still challenging for several reasons: RNA are inherently flexible and often adopt many energetically similar structures, which are not reliably distinguished by the available, incomplete thermodynamic model. Moreover, computationally, the problem is aggravated by the relevance of pseudoknots and non-canonical base pairs, which are hardly predicted efficiently. To identify nucleotides involved in pseudoknots and non-canonical interactions, we scrutinized the SHAPE reactivity of each nucleotide of the 188 nt long lariat-capping ribozyme under multiple conditions. Reactivities analyzed in the light of the X-ray structure were shown to report accurately the nucleotide status. Those that seemed paradoxical were rationalized by the nucleotide behavior along molecular dynamic simulations. We show that valuable information on intricate interactions can be deduced from probing with different reagents, and in the presence or absence of Mg2+. Furthermore, probing at increasing temperature was remarkably efficient at pointing to non-canonical interactions and pseudoknot pairings. The possibilities of following such strategies to inform structure modeling software are discussed.
Collapse
|
5
|
Thompson MK, Kiourlappou M, Davis I. Ribo-Pop: simple, cost-effective, and widely applicable ribosomal RNA depletion. RNA (NEW YORK, N.Y.) 2020; 26:1731-1742. [PMID: 32759389 PMCID: PMC7566562 DOI: 10.1261/rna.076562.120] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Accepted: 07/31/2020] [Indexed: 06/11/2023]
Abstract
The measurement of RNA abundance derived from massively parallel sequencing experiments is an essential technique. Methods that reduce ribosomal RNA levels are usually required prior to sequencing library construction because ribosomal RNA typically comprises the vast majority of a total RNA sample. For some experiments, ribosomal RNA depletion is favored over poly(A) selection because it offers a more inclusive representation of the transcriptome. However, methods to deplete ribosomal RNA are generally proprietary, complex, inefficient, applicable to only specific species, or compatible with only a narrow range of RNA input levels. Here, we describe Ribo-Pop (ribosomal RNA depletion for popular use), a simple workflow and antisense oligo design strategy that we demonstrate works over a wide input range and can be easily adapted to any organism with a sequenced genome. We provide a computational pipeline for probe selection, a streamlined 20-min protocol, and ready-to-use oligo sequences for several organisms. We anticipate that our simple and generalizable "open source" design strategy would enable virtually any laboratory to pursue full transcriptome sequencing in their organism of interest with minimal time and resource investment.
Collapse
Affiliation(s)
- Mary Kay Thompson
- Department of Biochemistry, University of Oxford, Oxford, OX1 3QU, United Kingdom
| | - Maria Kiourlappou
- Department of Biochemistry, University of Oxford, Oxford, OX1 3QU, United Kingdom
| | - Ilan Davis
- Department of Biochemistry, University of Oxford, Oxford, OX1 3QU, United Kingdom
| |
Collapse
|
6
|
Kasprzak WK, Ahmed NA, Shapiro BA. Modeling ligand docking to RNA in the design of RNA-based nanostructures. Curr Opin Biotechnol 2020; 63:16-25. [DOI: 10.1016/j.copbio.2019.10.010] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2019] [Accepted: 10/30/2019] [Indexed: 12/30/2022]
|
7
|
Meyer M, Walbott H, Oliéric V, Kondo J, Costa M, Masquida B. Conformational adaptation of UNCG loops upon crowding. RNA (NEW YORK, N.Y.) 2019; 25:1522-1531. [PMID: 31427457 PMCID: PMC6795143 DOI: 10.1261/rna.072694.119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/22/2019] [Accepted: 08/01/2019] [Indexed: 06/10/2023]
Abstract
If the A-form helix is the major structural motif found in RNA, the loops that cap them constitute the second most important family of motifs. Among those, two are overrepresented, GNRA and UNCG tetraloops. Recent surveys of RNA structures deposited in the PDB show that GNRA and UNCG tetraloops can adopt tertiary folds that are very different from their canonical conformations, characterized by the presence of a U-turn of a Z-turn, respectively. Crystallographic data from both a lariat-capping (LC) ribozyme and a group II intron ribozyme reveal that a given UUCG tetraloop can adopt a distinct fold depending on its structural environment. Specifically, when the crystal packing applies relaxed constraints on the loop, the canonical Z-turn conformation is observed. In contrast, a highly packed environment induces "squashing" of the tetraloop by distorting its sugar-phosphate backbone in a specific way that expels the first and fourth nucleobases out of the loop, and falls in van der Waals distance of the last base pair of the helix, taking the place of the pair formed between the first and fourth residues in Z-turn loops. The biological relevance of our observations is supported by the presence of similarly deformed loops in the highly packed environment of the ribosome and in a complex between a dsRNA and a RNase III. The finding that Z-turn loops change conformation under higher molecular packing suggests that, in addition to their demonstrated role in stabilizing RNA folding, they may contribute to the three-dimensional structure of RNA by mediating tertiary interactions with distal residues.
Collapse
Affiliation(s)
| | - Hélène Walbott
- Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Université Paris-Sud, Université Paris-Saclay, 91198, Gif-sur-Yvette cedex, France
| | - Vincent Oliéric
- Paul Scherrer Institute, Swiss Light Source, 5232 Villigen PSI, Switzerland
| | - Jiro Kondo
- Department of Materials and Life Sciences, Sophia University, 7-1 Kioi-cho, Chiyoda-ku, 102-8554 Tokyo, Japan
| | - Maria Costa
- Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Université Paris-Sud, Université Paris-Saclay, 91198, Gif-sur-Yvette cedex, France
| | - Benoît Masquida
- UMR7156 GMGM Université de Strasbourg - CNRS, 67084 Strasbourg, France
| |
Collapse
|
8
|
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
|
9
|
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
|
10
|
Plumridge A, Katz AM, Calvey GD, Elber R, Kirmizialtin S, Pollack L. Revealing the distinct folding phases of an RNA three-helix junction. Nucleic Acids Res 2018; 46:7354-7365. [PMID: 29762712 PMCID: PMC6101490 DOI: 10.1093/nar/gky363] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2018] [Revised: 04/17/2018] [Accepted: 04/24/2018] [Indexed: 01/08/2023] Open
Abstract
Remarkable new insight has emerged into the biological role of RNA in cells. RNA folding and dynamics enable many of these newly discovered functions, calling for an understanding of RNA self-assembly and conformational dynamics. Because RNAs pass through multiple structures as they fold, an ensemble perspective is required to visualize the flow through fleetingly populated sets of states. Here, we combine microfluidic mixing technology and small angle X-ray scattering (SAXS) to measure the Mg-induced folding of a small RNA domain, the tP5abc three helix junction. Our measurements are interpreted using ensemble optimization to select atomically detailed structures that recapitulate each experimental curve. Structural ensembles, derived at key stages in both time-resolved studies and equilibrium titrations, reproduce the features of known intermediates, and more importantly, offer a powerful new structural perspective on the time-progression of folding. Distinct collapse phases along the pathway appear to be orchestrated by specific interactions with Mg ions. These key interactions subsequently direct motions of the backbone that position the partners of tertiary contacts for later bonding, and demonstrate a remarkable synergy between Mg and RNA across numerous time-scales.
Collapse
Affiliation(s)
- Alex Plumridge
- School of Applied and Engineering Physics, Cornell University, Ithaca, NY, USA
| | - Andrea M Katz
- School of Applied and Engineering Physics, Cornell University, Ithaca, NY, USA
| | - George D Calvey
- School of Applied and Engineering Physics, Cornell University, Ithaca, NY, USA
| | - Ron Elber
- Department of Chemistry and Institute for Computational Engineering and Sciences (ICES) University of Texas at Austin, Austin, TX, USA
| | - Serdal Kirmizialtin
- Chemistry Program, Science Division, New York University Abu Dhabi, Abu Dhabi, United Arab Emirates
| | - Lois Pollack
- School of Applied and Engineering Physics, Cornell University, Ithaca, NY, USA
| |
Collapse
|
11
|
Š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
|
12
|
Liu L, Hyeon C. Contact Statistics Highlight Distinct Organizing Principles of Proteins and RNA. Biophys J 2017; 110:2320-2327. [PMID: 27276250 DOI: 10.1016/j.bpj.2016.04.020] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2016] [Revised: 03/15/2016] [Accepted: 04/01/2016] [Indexed: 11/18/2022] Open
Abstract
Although both RNA and proteins have densely packed native structures, chain organizations of these two biopolymers are fundamentally different. Motivated by the recent discoveries in chromatin folding that interphase chromosomes have territorial organization with signatures pointing to metastability, we analyzed the biomolecular structures deposited in the Protein Data Bank and found that the intrachain contact probabilities, P(s) as a function of the arc length s, decay in power-law ∼s(-γ) over the intermediate range of s, 10 ≲ s ≲ 110. We found that the contact probability scaling exponent is γ ≈ 1.11 for large RNA (N > 110), γ ≈ 1.41 for small-sized RNA (N < 110), and γ ≈ 1.65 for proteins. Given that Gaussian statistics is expected for a fully equilibrated chain in polymer melts, the deviation of γ-value from γ = 1.5 for the subchains of large RNA in the native state suggests that the chain configuration of RNA is not fully equilibrated. It is visually clear that folded structures of large-sized RNA (N ≳ 110) adopt crumpled structures, partitioned into modular multidomains assembled by proximal sequences along the chain, whereas the polypeptide chain of folded proteins looks better mixed with the rest of the structure. Our finding of γ ≈ 1 for large RNA might be an ineluctable consequence of the hierarchical ordering of the secondary to tertiary elements in the folding process.
Collapse
Affiliation(s)
- Lei Liu
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Republic of Korea
| | - Changbong Hyeon
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Republic of Korea.
| |
Collapse
|
13
|
Xu H, Plaut B, Zhu X, Chen M, Mavinkurve U, Maiti A, Song G, Murari K, Mandal M. Direct Observation of Folding Energy Landscape of RNA Hairpin at Mechanical Loading Rates. J Phys Chem B 2017; 121:2220-2229. [PMID: 28248503 DOI: 10.1021/acs.jpcb.6b10362] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
By applying a controlled mechanical load using optical tweezers, we measured the diffusive barrier crossing in a 49 nt long P5ab RNA hairpin. We find that in the free-energy landscape the barrier height (G‡) and transition distance (x‡) are dependent on the loading rate (r) along the pulling direction, x, as predicted by Bell. The barrier shifted toward the initial state, whereas ΔG‡ reduced significantly from 50 to 5 kT, as r increased from 0 to 32 pN/s. However, the equilibrium work (ΔG) during strand separation, as estimated by Crook's fluctuation theorem, remained unchanged at different rates. Previously, helix formation and denaturation have been described as two-state (F ↔ U) transitions for P5ab. Herein, we report three intermediate states I1, I, and I2 located at 4, 11, and 16 nm respectively, from the folded conformation. The intermediates were observed only when the hairpin was subjected to an optimal r, 7.6 pN/s. The results indicate that the complementary strands in P5ab can zip and unzip through complex routes, whereby mismatches act as checkpoints and often impose barriers. The study highlights the significance of loading rates in force-spectroscopy experiments that are increasingly being used to measure the folding properties of biomolecules.
Collapse
Affiliation(s)
- Huizhong Xu
- Department of Physics, ‡Department of Mathematical Sciences, §Department of Computer Science, and ∥Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - Benjamin Plaut
- Department of Physics, ‡Department of Mathematical Sciences, §Department of Computer Science, and ∥Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - Xiran Zhu
- Department of Physics, ‡Department of Mathematical Sciences, §Department of Computer Science, and ∥Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - Maverick Chen
- Department of Physics, ‡Department of Mathematical Sciences, §Department of Computer Science, and ∥Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - Udit Mavinkurve
- Department of Physics, ‡Department of Mathematical Sciences, §Department of Computer Science, and ∥Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - Anindita Maiti
- Department of Physics, ‡Department of Mathematical Sciences, §Department of Computer Science, and ∥Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - Guangtao Song
- Department of Physics, ‡Department of Mathematical Sciences, §Department of Computer Science, and ∥Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - Krishna Murari
- Department of Physics, ‡Department of Mathematical Sciences, §Department of Computer Science, and ∥Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - Maumita Mandal
- Department of Physics, ‡Department of Mathematical Sciences, §Department of Computer Science, and ∥Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
14
|
Welty R, Hall KB. Nucleobases Undergo Dynamic Rearrangements during RNA Tertiary Folding. J Mol Biol 2016; 428:4490-4502. [PMID: 27693721 DOI: 10.1016/j.jmb.2016.09.015] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2016] [Revised: 08/31/2016] [Accepted: 09/15/2016] [Indexed: 02/07/2023]
Abstract
The tertiary structure of the GTPase center (GAC) of 23S ribosomal RNA (rRNA) as seen in cocrystals is extremely compact. It is stabilized by long-range hydrogen bonds and nucleobase stacking and by a triloop that forms within its three-way junction. Its folding pathway from secondary structure to tertiary structure has not been previously observed, but it was shown to require Mg2+ ions in equilibrium experiments. The fluorescent nucleotide 2-aminopurine was substituted at selected sites within the 60-nt GAC. Fluorescence intensity changes upon addition of MgCl2 were monitored over a time-course from 1ms to 100s as the RNA folds. The folding pathway is revealed here to be hierarchical through several intermediates. Observation of the nucleobases during folding provides a new perspective on the process and the pathway, revealing the dynamics of nucleobase conformational exchange during the folding transitions.
Collapse
Affiliation(s)
- Robb Welty
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO 63110, USA.
| | - Kathleen B Hall
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO 63110, USA.
| |
Collapse
|
15
|
Xue Y, Gracia B, Herschlag D, Russell R, Al-Hashimi HM. Visualizing the formation of an RNA folding intermediate through a fast highly modular secondary structure switch. Nat Commun 2016; 7:ncomms11768. [PMID: 27292179 PMCID: PMC4909931 DOI: 10.1038/ncomms11768] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2015] [Accepted: 04/26/2016] [Indexed: 12/28/2022] Open
Abstract
Intermediates play important roles in RNA folding but can be difficult to characterize when short-lived or not significantly populated. By combining (15)N relaxation dispersion NMR with chemical probing, we visualized a fast (kex=k1+k-1≈423 s(-1)) secondary structural switch directed towards a low-populated (∼3%) partially folded intermediate in tertiary folding of the P5abc subdomain of the 'Tetrahymena' group I intron ribozyme. The secondary structure switch changes the base-pairing register across the P5c hairpin, creating a native-like structure, and occurs at rates of more than two orders of magnitude faster than tertiary folding. The switch occurs robustly in the absence of tertiary interactions, Mg(2+) or even when the hairpin is excised from the three-way junction. Fast, highly modular secondary structural switches may be quite common during RNA tertiary folding where they may help smoothen the folding landscape by allowing folding to proceed efficiently via additional pathways.
Collapse
Affiliation(s)
- Yi Xue
- Department of Biochemistry, Duke Center for RNA Biology, Duke University Medical Center, Durham, North Carolina 27710, USA
| | - Brant Gracia
- Department of Molecular Biosciences, Institute for Cellular and Molecular Biology, University of Texas at Austin, Austin, Texas 78712, USA
| | - Daniel Herschlag
- Department of Biochemistry, Beckman Center, Stanford University, Stanford, California 94305, USA.,Department of Chemistry, Stanford University, Stanford, California 94305, USA.,Department of Chemical Engineering, Stanford University, Stanford, California 94305, USA.,Chemistry, Engineering, and Medicine for Human Health (ChEM-H) Institute, Stanford University, Stanford, California 94305, USA
| | - Rick Russell
- Department of Molecular Biosciences, Institute for Cellular and Molecular Biology, University of Texas at Austin, Austin, Texas 78712, USA
| | - Hashim M Al-Hashimi
- Department of Biochemistry, Duke Center for RNA Biology, Duke University Medical Center, Durham, North Carolina 27710, USA.,Department of Chemistry, Duke University, Durham, Stanford, North Carolina 27710, USA
| |
Collapse
|
16
|
Xu X, Yu T, Chen SJ. Understanding the kinetic mechanism of RNA single base pair formation. Proc Natl Acad Sci U S A 2016; 113:116-21. [PMID: 26699466 PMCID: PMC4711849 DOI: 10.1073/pnas.1517511113] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.
Collapse
Affiliation(s)
- Xiaojun Xu
- Department of Physics, University of Missouri, Columbia, MO 65211; Department of Biochemistry, University of Missouri, Columbia, MO 65211; Informatics Institute, University of Missouri, Columbia, MO 65211
| | - Tao Yu
- Department of Physics, University of Missouri, Columbia, MO 65211; Department of Biochemistry, University of Missouri, Columbia, MO 65211; Informatics Institute, University of Missouri, Columbia, MO 65211; Department of Physics, Jianghan University, Wuhan, Hubei 430056, China
| | - Shi-Jie Chen
- Department of Physics, University of Missouri, Columbia, MO 65211; Department of Biochemistry, University of Missouri, Columbia, MO 65211; Informatics Institute, University of Missouri, Columbia, MO 65211;
| |
Collapse
|
17
|
Angulo J, Ulryck N, Deforges J, Chamond N, Lopez-Lastra M, Masquida B, Sargueil B. LOOP IIId of the HCV IRES is essential for the structural rearrangement of the 40S-HCV IRES complex. Nucleic Acids Res 2015; 44:1309-25. [PMID: 26626152 PMCID: PMC4756818 DOI: 10.1093/nar/gkv1325] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2015] [Accepted: 11/11/2015] [Indexed: 12/14/2022] Open
Abstract
As obligatory intracellular parasites, viruses rely on cellular machines to complete their life cycle, and most importantly they recruit the host ribosomes to translate their mRNA. The Hepatitis C viral mRNA initiates translation by directly binding the 40S ribosomal subunit in such a way that the initiation codon is correctly positioned in the P site of the ribosome. Such a property is likely to be central for many viruses, therefore the description of host-pathogen interaction at the molecular level is instrumental to provide new therapeutic targets. In this study, we monitored the 40S ribosomal subunit and the viral RNA structural rearrangement induced upon the formation of the binary complex. We further took advantage of an IRES viral mutant mRNA deficient for translation to identify the interactions necessary to promote translation. Using a combination of structure probing in solution and molecular modeling we establish a whole atom model which appears to be very similar to the one obtained recently by cryoEM. Our model brings new information on the complex, and most importantly reveals some structural rearrangement within the ribosome. This study suggests that the formation of a ‘kissing complex’ between the viral RNA and the 18S ribosomal RNA locks the 40S ribosomal subunit in a conformation proficient for translation.
Collapse
Affiliation(s)
- Jenniffer Angulo
- CNRS UMR 8015, Laboratoire de cristallographie et RMN Biologiques, Université Paris Descartes, 4 avenue de l'Observatoire, 75270 Paris Cedex 06, France Laboratorio de Virología Molecular, Instituto Milenio de Inmunología e Inmunoterapia, Centro de Investigaciones Médicas, Escuela de Medicina, Pontificia Universidad Católica de Chile, Marcoleta 391, Santiago, Chile
| | - Nathalie Ulryck
- CNRS UMR 8015, Laboratoire de cristallographie et RMN Biologiques, Université Paris Descartes, 4 avenue de l'Observatoire, 75270 Paris Cedex 06, France
| | - Jules Deforges
- CNRS UMR 8015, Laboratoire de cristallographie et RMN Biologiques, Université Paris Descartes, 4 avenue de l'Observatoire, 75270 Paris Cedex 06, France
| | - Nathalie Chamond
- CNRS UMR 8015, Laboratoire de cristallographie et RMN Biologiques, Université Paris Descartes, 4 avenue de l'Observatoire, 75270 Paris Cedex 06, France
| | - Marcelo Lopez-Lastra
- Laboratorio de Virología Molecular, Instituto Milenio de Inmunología e Inmunoterapia, Centro de Investigaciones Médicas, Escuela de Medicina, Pontificia Universidad Católica de Chile, Marcoleta 391, Santiago, Chile
| | - Benoît Masquida
- UMR 7156 Génétique Moléculaire Génomique Microbiologie, CNRS - Université de Strasbourg, Strasbourg, France
| | - Bruno Sargueil
- CNRS UMR 8015, Laboratoire de cristallographie et RMN Biologiques, Université Paris Descartes, 4 avenue de l'Observatoire, 75270 Paris Cedex 06, France
| |
Collapse
|
18
|
Mustoe AM, Al-Hashimi HM, Brooks CL. Secondary structure encodes a cooperative tertiary folding funnel in the Azoarcus ribozyme. Nucleic Acids Res 2015; 44:402-12. [PMID: 26481360 PMCID: PMC4705646 DOI: 10.1093/nar/gkv1055] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2015] [Accepted: 10/03/2015] [Indexed: 12/20/2022] Open
Abstract
A requirement for specific RNA folding is that the free-energy landscape discriminate against non-native folds. While tertiary interactions are critical for stabilizing the native fold, they are relatively non-specific, suggesting additional mechanisms contribute to tertiary folding specificity. In this study, we use coarse-grained molecular dynamics simulations to explore how secondary structure shapes the tertiary free-energy landscape of the Azoarcus ribozyme. We show that steric and connectivity constraints posed by secondary structure strongly limit the accessible conformational space of the ribozyme, and that these so-called topological constraints in turn pose strong free-energy penalties on forming different tertiary contacts. Notably, native A-minor and base-triple interactions form with low conformational free energy, while non-native tetraloop/tetraloop–receptor interactions are penalized by high conformational free energies. Topological constraints also give rise to strong cooperativity between distal tertiary interactions, quantitatively matching prior experimental measurements. The specificity of the folding landscape is further enhanced as tertiary contacts place additional constraints on the conformational space, progressively funneling the molecule to the native state. These results indicate that secondary structure assists the ribozyme in navigating the otherwise rugged tertiary folding landscape, and further emphasize topological constraints as a key force in RNA folding.
Collapse
Affiliation(s)
- Anthony M Mustoe
- Department of Biophysics, University of Michigan, Ann Arbor, MI 48109, USA
| | - Hashim M Al-Hashimi
- Department of Biochemistry and Chemistry, Duke University School of Medicine, Durham, NC 27710, USA
| | - Charles L Brooks
- Department of Biophysics, University of Michigan, Ann Arbor, MI 48109, USA Department of Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
| |
Collapse
|
19
|
Shaw E, St-Pierre P, McCluskey K, Lafontaine DA, Penedo JC. Using sm-FRET and denaturants to reveal folding landscapes. Methods Enzymol 2015; 549:313-41. [PMID: 25432755 DOI: 10.1016/b978-0-12-801122-5.00014-3] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
RNA folding studies aim to clarify the relationship among sequence, tridimensional structure, and biological function. In the last decade, the application of single-molecule fluorescence resonance energy transfer (sm-FRET) techniques to investigate RNA structure and folding has revealed the details of conformational changes and timescale of the process leading to the formation of biologically active RNA structures with subnanometer resolution on millisecond timescales. In this review, we initially summarize the first wave of single-molecule FRET-based RNA techniques that focused on analyzing the influence of mono- and divalent metal ions on RNA function, and how these studies have provided very valuable information about folding pathways and the presence of intermediate and low-populated states. Next, we describe a second generation of single-molecule techniques that combine sm-FRET with the use of chemical denaturants as an emerging powerful approach to reveal information about the dynamics and energetics of RNA folding that remains hidden using conventional sm-FRET approaches. The main advantages of using the competing interplay between folding agents such as metal ions and denaturants to observe and manipulate the dynamics of RNA folding and RNA-ligand interactions is discussed in the context of the adenine riboswitch aptamer.
Collapse
Affiliation(s)
- Euan Shaw
- SUPA School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife, United Kingdom
| | - Patrick St-Pierre
- RNA Group, Department of Biology, Faculty of Science, Université de Sherbrooke, Sherbrooke, Québec, Canada
| | - Kaley McCluskey
- SUPA School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife, United Kingdom
| | - Daniel A Lafontaine
- RNA Group, Department of Biology, Faculty of Science, Université de Sherbrooke, Sherbrooke, Québec, Canada.
| | - J Carlos Penedo
- SUPA School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife, United Kingdom; Biomedical Sciences Research Complex, University of St. Andrews, St. Andrews, Fife, United Kingdom.
| |
Collapse
|
20
|
MD Simulations of tRNA and Aminoacyl-tRNA Synthetases: Dynamics, Folding, Binding, and Allostery. Int J Mol Sci 2015; 16:15872-902. [PMID: 26184179 PMCID: PMC4519929 DOI: 10.3390/ijms160715872] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2015] [Revised: 07/07/2015] [Accepted: 07/08/2015] [Indexed: 12/21/2022] Open
Abstract
While tRNA and aminoacyl-tRNA synthetases are classes of biomolecules that have been extensively studied for decades, the finer details of how they carry out their fundamental biological functions in protein synthesis remain a challenge. Recent molecular dynamics (MD) simulations are verifying experimental observations and providing new insight that cannot be addressed from experiments alone. Throughout the review, we briefly discuss important historical events to provide a context for how far the field has progressed over the past few decades. We then review the background of tRNA molecules, aminoacyl-tRNA synthetases, and current state of the art MD simulation techniques for those who may be unfamiliar with any of those fields. Recent MD simulations of tRNA dynamics and folding and of aminoacyl-tRNA synthetase dynamics and mechanistic characterizations are discussed. We highlight the recent successes and discuss how important questions can be addressed using current MD simulations techniques. We also outline several natural next steps for computational studies of AARS:tRNA complexes.
Collapse
|
21
|
Mustoe AM, Brooks CL, Al-Hashimi HM. Topological constraints are major determinants of tRNA tertiary structure and dynamics and provide basis for tertiary folding cooperativity. Nucleic Acids Res 2014; 42:11792-804. [PMID: 25217593 PMCID: PMC4191394 DOI: 10.1093/nar/gku807] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Recent studies have shown that basic steric and connectivity constraints encoded at the secondary structure level are key determinants of 3D structure and dynamics in simple two-way RNA junctions. However, the role of these topological constraints in higher order RNA junctions remains poorly understood. Here, we use a specialized coarse-grained molecular dynamics model to directly probe the thermodynamic contributions of topological constraints in defining the 3D architecture and dynamics of transfer RNA (tRNA). Topological constraints alone restrict tRNA's allowed conformational space by over an order of magnitude and strongly discriminate against formation of non-native tertiary contacts, providing a sequence independent source of folding specificity. Topological constraints also give rise to long-range correlations between the relative orientation of tRNA's helices, which in turn provides a mechanism for encoding thermodynamic cooperativity between distinct tertiary interactions. These aspects of topological constraints make it such that only several tertiary interactions are needed to confine tRNA to its native global structure and specify functionally important 3D dynamics. We further show that topological constraints are conserved across tRNA's different naturally occurring secondary structures. Taken together, our results emphasize the central role of secondary-structure-encoded topological constraints in defining RNA 3D structure, dynamics and folding.
Collapse
Affiliation(s)
- Anthony M Mustoe
- Department of Biophysics, University of Michigan, Ann Arbor, MI 48109, USA
| | - Charles L Brooks
- Department of Biophysics, University of Michigan, Ann Arbor, MI 48109, USA Department of Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
| | - Hashim M Al-Hashimi
- Department of Biochemistry and Chemistry, Duke University School of Medicine, Durham, NC 27710, USA
| |
Collapse
|
22
|
Tanaka T, Furuta H, Ikawa Y. Installation of orthogonality to the interface that assembles two modular domains in the Tetrahymena group I ribozyme. J Biosci Bioeng 2014; 117:407-12. [DOI: 10.1016/j.jbiosc.2013.10.008] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2013] [Revised: 10/04/2013] [Accepted: 10/07/2013] [Indexed: 01/08/2023]
|
23
|
Abstract
RNA dynamics play a fundamental role in many cellular functions. However, there is no general framework to describe these complex processes, which typically consist of many structural maneuvers that occur over timescales ranging from picoseconds to seconds. Here, we classify RNA dynamics into distinct modes representing transitions between basins on a hierarchical free-energy landscape. These transitions include large-scale secondary-structural transitions at >0.1-s timescales, base-pair/tertiary dynamics at microsecond-to-millisecond timescales, stacking dynamics at timescales ranging from nanoseconds to microseconds, and other "jittering" motions at timescales ranging from picoseconds to nanoseconds. We review various modes within these three different tiers, the different mechanisms by which they are used to regulate function, and how they can be coupled together to achieve greater functional complexity.
Collapse
|
24
|
Mustoe AM, Al-Hashimi HM, Brooks CL. Coarse grained models reveal essential contributions of topological constraints to the conformational free energy of RNA bulges. J Phys Chem B 2014; 118:2615-27. [PMID: 24547945 PMCID: PMC3983386 DOI: 10.1021/jp411478x] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
![]()
Recent studies have shown that simple
stereochemical constraints
encoded at the RNA secondary structure level significantly restrict
the orientation of RNA helices across two-way junctions and yield
physically reasonable distributions of RNA 3D conformations. Here
we develop a new coarse-grain model, TOPRNA, that is optimized for
exploring detailed aspects of these topological constraints in complex
RNA systems. Unlike prior models, TOPRNA effectively treats RNAs as
collections of semirigid helices linked by freely rotatable single
strands, allowing us to isolate the effects of secondary structure
connectivity and sterics on 3D structure. Simulations of bulge junctions
show that TOPRNA captures new aspects of topological constraints,
including variations arising from deviations in local A-form structure,
translational displacements of the helices, and stereochemical constraints
imposed by bulge-linker nucleotides. Notably, these aspects of topological
constraints define free energy landscapes that coincide with the distribution
of bulge conformations in the PDB. Our simulations also quantitatively
reproduce NMR RDC measurements made on HIV-1 TAR at low salt concentrations,
although not for different TAR mutants or at high salt concentrations.
Our results confirm that topological constraints are an important
determinant of bulge conformation and dynamics and demonstrate the
utility of TOPRNA for studying the topological constraints of complex
RNAs.
Collapse
Affiliation(s)
- Anthony M Mustoe
- Departments of Biophysics and ‡Chemistry, University of Michigan , 930 North University Avenue, Ann Arbor, Michigan 48109, United States
| | | | | |
Collapse
|
25
|
Li R, Ge HW, Cho SS. Sequence-dependent base-stacking stabilities guide tRNA folding energy landscapes. J Phys Chem B 2013; 117:12943-52. [PMID: 23841777 DOI: 10.1021/jp402114p] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The folding of bacterial tRNAs with disparate sequences has been observed to proceed in distinct folding mechanisms despite their structural similarity. To explore the folding landscapes of tRNA, we performed ion concentration-dependent coarse-grained TIS model MD simulations of several E. coli tRNAs to compare their thermodynamic melting profiles to the classical absorbance spectra of Crothers and co-workers. To independently validate our findings, we also performed atomistic empirical force field MD simulations of tRNAs, and we compared the base-to-base distances from coarse-grained and atomistic MD simulations to empirical base-stacking free energies. We then projected the free energies to the secondary structural elements of tRNA, and we observe distinct, parallel folding mechanisms whose differences can be inferred on the basis of their sequence-dependent base-stacking stabilities. In some cases, a premature, nonproductive folding intermediate corresponding to the Ψ hairpin loop must backtrack to the unfolded state before proceeding to the folded state. This observation suggests a possible explanation for the fast and slow phases observed in tRNA folding kinetics.
Collapse
Affiliation(s)
- Rongzhong Li
- Department of Physics, Wake Forest University , Winston-Salem, North Carolina 27106, United States
| | | | | |
Collapse
|
26
|
Pabit SA, Sutton JL, Chen H, Pollack L. Role of ion valence in the submillisecond collapse and folding of a small RNA domain. Biochemistry 2013; 52:1539-46. [PMID: 23398396 DOI: 10.1021/bi3016636] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Following the addition of ions to trigger folding, RNA molecules undergo a transition from rigid, extended states to a compact ensemble. Determining the time scale for this collapse provides important insights into electrostatic contributions to RNA folding; however, it can be challenging to isolate the effects of purely nonspecific collapse, e.g., relaxation due to backbone charge compensation, from the concurrent formation of some tertiary contacts. To solve this problem, we decoupled nonspecific collapse from tertiary folding using a single-point mutation to eliminate tertiary contacts in the small RNA subdomain known as tP5abc. Microfluidic mixing with microsecond time resolution and Förster resonance energy transfer detection provides insight into the ionic strength-dependent transition from extended to compact ensembles. Differences in reaction rates are detected when folding is initiated by monovalent or divalent ions, consistent with equilibrium measurements illustrating the enhanced screening of divalent ions relative to monovalent ions at the same ionic strength. Ion-driven collapse is fast, and a comparison of the collapse time of the wild-type and mutant tP5abc suggests that site binding of Mg(2+) occurs on submillisecond time scales.
Collapse
Affiliation(s)
- Suzette A Pabit
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, USA
| | | | | | | |
Collapse
|
27
|
Ray SS, Pal SK. RNA secondary structure prediction using soft computing. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2013; 10:2-17. [PMID: 23702539 DOI: 10.1109/tcbb.2012.159] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
Prediction of RNA structure is invaluable in creating new drugs and understanding genetic diseases. Several deterministic algorithms and soft computing-based techniques have been developed for more than a decade to determine the structure from a known RNA sequence. Soft computing gained importance with the need to get approximate solutions for RNA sequences by considering the issues related with kinetic effects, cotranscriptional folding, and estimation of certain energy parameters. A brief description of some of the soft computing-based techniques, developed for RNA secondary structure prediction, is presented along with their relevance. The basic concepts of RNA and its different structural elements like helix, bulge, hairpin loop, internal loop, and multiloop are described. These are followed by different methodologies, employing genetic algorithms, artificial neural networks, and fuzzy logic. The role of various metaheuristics, like simulated annealing, particle swarm optimization, ant colony optimization, and tabu search is also discussed. A relative comparison among different techniques, in predicting 12 known RNA secondary structures, is presented, as an example. Future challenging issues are then mentioned.
Collapse
|