1
|
Rissone P, Severino A, Pastor I, Ritort F. Universal cold RNA phase transitions. Proc Natl Acad Sci U S A 2024; 121:e2408313121. [PMID: 39150781 PMCID: PMC11348302 DOI: 10.1073/pnas.2408313121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2024] [Accepted: 07/12/2024] [Indexed: 08/18/2024] Open
Abstract
RNA's diversity of structures and functions impacts all life forms since primordia. We use calorimetric force spectroscopy to investigate RNA folding landscapes in previously unexplored low-temperature conditions. We find that Watson-Crick RNA hairpins, the most basic secondary structure elements, undergo a glass-like transition below [Formula: see text]C where the heat capacity abruptly changes and the RNA folds into a diversity of misfolded structures. We hypothesize that an altered RNA biochemistry, determined by sequence-independent ribose-water interactions, outweighs sequence-dependent base pairing. The ubiquitous ribose-water interactions lead to universal RNA phase transitions below TG, such as maximum stability at [Formula: see text]C where water density is maximum, and cold denaturation at [Formula: see text]C. RNA cold biochemistry may have a profound impact on RNA function and evolution.
Collapse
Affiliation(s)
- Paolo Rissone
- Small Biosystems Lab, Condensed Matter Physics Department, Universitat de Barcelona, Barcelona 08028, Spain
| | - Aurélien Severino
- Small Biosystems Lab, Condensed Matter Physics Department, Universitat de Barcelona, Barcelona 08028, Spain
| | - Isabel Pastor
- Small Biosystems Lab, Condensed Matter Physics Department, Universitat de Barcelona, Barcelona 08028, Spain
| | - Felix Ritort
- Small Biosystems Lab, Condensed Matter Physics Department, Universitat de Barcelona, Barcelona 08028, Spain
- Institut de Nanociència i Nanotecnologia, Universitat de Barcelona, Barcelona 08028, Spain
| |
Collapse
|
2
|
Bukina V, Božič A. Context-dependent structure formation of hairpin motifs in bacteriophage MS2 genomic RNA. Biophys J 2024:S0006-3495(24)00526-5. [PMID: 39118324 DOI: 10.1016/j.bpj.2024.08.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2024] [Revised: 07/17/2024] [Accepted: 08/05/2024] [Indexed: 08/10/2024] Open
Abstract
Many functions of ribonucleic acid (RNA) rely on its ability to assume specific sequence-structure motifs. Packaging signals found in certain RNA viruses are one such prominent example of functional RNA motifs. These signals are short hairpin loops that interact with coat proteins and drive viral self-assembly. As they are found in different positions along the much longer genomic RNA, the formation of their correct structure occurs as a part of a larger context. Any changes to this context can consequently lead to changes in the structure of the motifs themselves. In fact, previous studies have shown that structure and function of RNA motifs can be highly context sensitive to the flanking sequence surrounding them. However, in what ways different flanking sequences influence the structure of an RNA motif they surround has yet to be studied in detail. We focus on a hairpin-rich region of the RNA genome of bacteriophage MS2-a well-studied RNA virus with a wide potential for use in biotechnology-and systematically examine context-dependent structural stability of 14 previously identified hairpin motifs, which include putative and confirmed packaging signals. Combining secondary and tertiary RNA structure prediction of the hairpin motifs placed in different contexts, ranging from the native genomic sequence to random RNA sequences and unstructured poly-U sequences, we determine different measures of motif structural stability. In this way, we show that while some motif structures can be stable in any context, others require specific context provided by the genome. Our results demonstrate the importance of context in RNA structure formation and how changes in the flanking sequence of an RNA motif sometimes lead to drastic changes in its structure. Structural stability of a motif in different contexts could provide additional insights into its functionality as well as assist in determining whether it remains functional when intentionally placed in other contexts.
Collapse
Affiliation(s)
- Veronika Bukina
- Department of Theoretical Physics, Jožef Stefan Institute, Ljubljana, Slovenia; Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana, Slovenia
| | - Anže Božič
- Department of Theoretical Physics, Jožef Stefan Institute, Ljubljana, Slovenia.
| |
Collapse
|
3
|
Tarafder S, Bhattacharya D. lociPARSE: a locality-aware invariant point attention model for scoring RNA 3D structures. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.11.04.565599. [PMID: 37961488 PMCID: PMC10635153 DOI: 10.1101/2023.11.04.565599] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/15/2023]
Abstract
A scoring function that can reliably assess the accuracy of a 3D RNA structural model in the absence of experimental structure is not only important for model evaluation and selection but also useful for scoring-guided conformational sampling. However, high-fidelity RNA scoring has proven to be difficult using conventional knowledge-based statistical potentials and currently-available machine learning-based approaches. Here we present lociPARSE, a locality-aware invariant point attention architecture for scoring RNA 3D structures. Unlike existing machine learning methods that estimate superposition-based root mean square deviation (RMSD), lociPARSE estimates Local Distance Difference Test (lDDT) scores capturing the accuracy of each nucleotide and its surrounding local atomic environment in a superposition-free manner, before aggregating information to predict global structural accuracy. Tested on multiple datasets including CASP15, lociPARSE significantly outperforms existing statistical potentials (rsRNASP, cgRNASP, DFIRE-RNA, and RASP) and machine learning methods (ARES and RNA3DCNN) across complementary assessment metrics. lociPARSE is freely available at https://github.com/Bhattacharya-Lab/lociPARSE.
Collapse
Affiliation(s)
- Sumit Tarafder
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia, 24061, USA
| | | |
Collapse
|
4
|
Jang SS, Ray KK, Lynall DG, Shepard KL, Nuckolls C, Gonzalez RL. RNA adapts its flexibility to efficiently fold and resist unfolding. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.27.595525. [PMID: 38853856 PMCID: PMC11160689 DOI: 10.1101/2024.05.27.595525] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2024]
Abstract
Recent studies have demonstrated that the mechanisms through which biopolymers like RNA interconvert between multiple folded structures are critical for their cellular functions. A major obstacle to elucidating these mechanisms is the lack of experimental approaches that can resolve these interconversions between functionally relevant biomolecular structures. Here, using a nano-electronic device with microsecond time resolution, we dissect the complete set of structural rearrangements executed by an ultra-stable RNA, the UUCG stem-loop, at the single-molecule level. We show that the stem-loop samples at least four conformations along two folding pathways leading to two distinct folded structures, only one of which has been previously observed. By modulating its flexibility, the stem-loop can adaptively select between these pathways, enabling it to both fold rapidly and resist unfolding. This paradigm of stabilization through compensatory changes in flexibility broadens our understanding of stable RNA structures and is expected to serve as a general strategy employed by all biopolymers.
Collapse
Affiliation(s)
- Sukjin S. Jang
- Department of Chemistry, Columbia University, New York, NY 10027 USA
| | - Korak Kumar Ray
- Department of Chemistry, Columbia University, New York, NY 10027 USA
| | - David G. Lynall
- Department of Electrical Engineering, Columbia University, New York, NY 10027 USA
| | - Kenneth L. Shepard
- Department of Electrical Engineering, Columbia University, New York, NY 10027 USA
| | - Colin Nuckolls
- Department of Chemistry, Columbia University, New York, NY 10027 USA
| | - Ruben L. Gonzalez
- Department of Chemistry, Columbia University, New York, NY 10027 USA
| |
Collapse
|
5
|
Rinaldi S, Moroni E, Rozza R, Magistrato A. Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation. J Chem Theory Comput 2024; 20:993-1018. [PMID: 38287883 DOI: 10.1021/acs.jctc.3c01239] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2024]
Abstract
Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.
Collapse
Affiliation(s)
- Silvia Rinaldi
- National Research Council of Italy (CNR) - Institute of Chemistry of OrganoMetallic Compounds (ICCOM), c/o Area di Ricerca CNR di Firenze Via Madonna del Piano 10, 50019 Sesto Fiorentino, Florence, Italy
| | - Elisabetta Moroni
- National Research Council of Italy (CNR) - Institute of Chemical Sciences and Technologies (SCITEC), via Mario Bianco 9, 20131 Milano, Italy
| | - Riccardo Rozza
- National Research Council of Italy (CNR) - Institute of Material Foundry (IOM) c/o International School for Advanced Studies (SISSA), Via Bonomea, 265, 34136 Trieste, Italy
| | - Alessandra Magistrato
- National Research Council of Italy (CNR) - Institute of Material Foundry (IOM) c/o International School for Advanced Studies (SISSA), Via Bonomea, 265, 34136 Trieste, Italy
| |
Collapse
|
6
|
Bose R, Saleem I, Mustoe AM. Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states. Cell Chem Biol 2024; 31:17-35. [PMID: 38199037 PMCID: PMC10842484 DOI: 10.1016/j.chembiol.2023.12.010] [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] [Received: 09/11/2023] [Revised: 11/21/2023] [Accepted: 12/12/2023] [Indexed: 01/12/2024]
Abstract
RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.
Collapse
Affiliation(s)
- Ritwika Bose
- Therapeutic Innovation Center (THINC), Department of Biochemistry and Molecular Pharmacology, Baylor College of Medicine, Houston, TX, USA
| | - Irfana Saleem
- Therapeutic Innovation Center (THINC), Department of Biochemistry and Molecular Pharmacology, Baylor College of Medicine, Houston, TX, USA
| | - Anthony M Mustoe
- Therapeutic Innovation Center (THINC), Department of Biochemistry and Molecular Pharmacology, Baylor College of Medicine, Houston, TX, USA; Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA.
| |
Collapse
|
7
|
Sarrazin-Gendron R, Waldispühl J, Reinharz V. Classification and Identification of Non-canonical Base Pairs and Structural Motifs. Methods Mol Biol 2024; 2726:143-168. [PMID: 38780731 DOI: 10.1007/978-1-0716-3519-3_7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
The 3D structures of many ribonucleic acid (RNA) loops are characterized by highly organized networks of non-canonical interactions. Multiple computational methods have been developed to annotate structures with those interactions or automatically identify recurrent interaction networks. By contrast, the reverse problem that aims to retrieve the geometry of a look from its sequence or ensemble of interactions remains much less explored. In this chapter, we will describe how to retrieve and build families of conserved structural motifs using their underlying network of non-canonical interactions. Then, we will show how to assign sequence alignments to those families and use the software BayesPairing to build statistical models of structural motifs with their associated sequence alignments. From this model, we will apply BayesPairing to identify in new sequences regions where those loop geometries can occur.
Collapse
Affiliation(s)
| | | | - Vladimir Reinharz
- Department of Computer Science, Université du Québec à Montréal, Montreal, QC, Canada.
| |
Collapse
|
8
|
Hori N, Thirumalai D. Watching ion-driven kinetics of ribozyme folding and misfolding caused by energetic and topological frustration one molecule at a time. Nucleic Acids Res 2023; 51:10737-10751. [PMID: 37758176 PMCID: PMC10602927 DOI: 10.1093/nar/gkad755] [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] [Received: 04/17/2023] [Revised: 08/23/2023] [Accepted: 09/05/2023] [Indexed: 10/03/2023] Open
Abstract
Folding of ribozymes into well-defined tertiary structures usually requires divalent cations. How Mg2+ ions direct the folding kinetics has been a long-standing unsolved problem because experiments cannot detect the positions and dynamics of ions. To address this problem, we used molecular simulations to dissect the folding kinetics of the Azoarcus ribozyme by monitoring the path each molecule takes to reach the folded state. We quantitatively establish that Mg2+ binding to specific sites, coupled with counter-ion release of monovalent cations, stimulate the formation of secondary and tertiary structures, leading to diverse pathways that include direct rapid folding and trapping in misfolded structures. In some molecules, key tertiary structural elements form when Mg2+ ions bind to specific RNA sites at the earliest stages of the folding, leading to specific collapse and rapid folding. In others, the formation of non-native base pairs, whose rearrangement is needed to reach the folded state, is the rate-limiting step. Escape from energetic traps, driven by thermal fluctuations, occurs readily. In contrast, the transition to the native state from long-lived topologically trapped native-like metastable states is extremely slow. Specific collapse and formation of energetically or topologically frustrated states occur early in the assembly process.
Collapse
Affiliation(s)
- Naoto Hori
- Department of Chemistry, University of Texas, Austin, TX 78712, USA
- School of Pharmacy, University of Nottingham, Nottingham, UK
| | - D Thirumalai
- Department of Chemistry, University of Texas, Austin, TX 78712, USA
- Department of Physics, University of Texas, Austin, TX 78712, USA
| |
Collapse
|
9
|
Kaiser C, Vogel M, Appel B, Weigand J, Müller S, Suess B, Wachtveitl J. Magnesium Ion-Driven Folding and Conformational Switching Kinetics of Tetracycline Binding Aptamer: Implications for in vivo Riboswitch Engineering. J Mol Biol 2023; 435:168253. [PMID: 37640152 DOI: 10.1016/j.jmb.2023.168253] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Revised: 08/22/2023] [Accepted: 08/22/2023] [Indexed: 08/31/2023]
Abstract
Engineering in vitro selected RNA aptamers into in vivo functional riboswitches represents a long-standing challenge in molecular biology. The highly specific aptamer domain of the riboswitch undergoes a conformational adjustment in response to ligand sensing, which in turn exerts the regulatory function. Besides essential factors like structural complexity and ligand binding kinetics, the active role of magnesium ions in stabilizing RNA tertiary structures and assisting in ligand binding can be a vital criterion. We present spectroscopic studies on the magnesium ion-driven folding of the Tetracycline binding aptamer. Using fluorescent labels, the aptamer pre-folding and subsequent ligand binding is monitored by magnesium titration experiments and time-resolved stopped-flow measurements. A minimum concentration of 0.5 mM magnesium is required to fold into a magnesium ion-stabilized binding-competent state with a preformed binding pocket. Tetracycline binding causes a pronounced conformational change that results in the establishment of the triple helix core motif, and that further propagates towards the closing stem. By a dynamic acquisition of magnesium ions, a kink motif is formed at the intersection of the triple helix and closing stem regions. This ultimately entails a stabilization of the closing stem which is discussed as a key element in the regulatory function of the Tetracycline aptamer.
Collapse
Affiliation(s)
- Christoph Kaiser
- Institute for Physical and Theoretical Chemistry, Goethe University, Frankfurt/Main, Max-von-Laue Str. 9, D-60438, Germany.
| | - Marc Vogel
- Department of Biology, Technical University Darmstadt, Darmstadt, Schnittspahnstraße 10, D-64287 Darmstadt, Germany
| | - Bettina Appel
- Institute of Biochemistry, University of Greifswald, Felix-Hausdorff-Str. 4, D-17489 Greifswald, Germany
| | - Julia Weigand
- Department of Biology, Technical University Darmstadt, Darmstadt, Schnittspahnstraße 10, D-64287 Darmstadt, Germany; Institute of Pharmaceutical Chemistry, Philipps-University Marburg, Marburg, Marbacher Weg 6, D-35037, Germany. https://twitter.com/WachtveitlLab
| | - Sabine Müller
- Institute of Biochemistry, University of Greifswald, Felix-Hausdorff-Str. 4, D-17489 Greifswald, Germany
| | - Beatrix Suess
- Department of Biology, Technical University Darmstadt, Darmstadt, Schnittspahnstraße 10, D-64287 Darmstadt, Germany; Centre for Synthetic Biology, Technical University Darmstadt, Darmstadt, Germany.
| | - Josef Wachtveitl
- Institute for Physical and Theoretical Chemistry, Goethe University, Frankfurt/Main, Max-von-Laue Str. 9, D-60438, Germany.
| |
Collapse
|
10
|
Siddika MA, Oi H, Hidaka K, Sugiyama H, Endo M, Matsumura S, Ikawa Y. Structural Expansion of Catalytic RNA Nanostructures through Oligomerization of a Cyclic Trimer of Engineered Ribozymes. Molecules 2023; 28:6465. [PMID: 37764241 PMCID: PMC10535472 DOI: 10.3390/molecules28186465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2023] [Revised: 08/29/2023] [Accepted: 08/30/2023] [Indexed: 09/29/2023] Open
Abstract
The multimolecular assembly of three-dimensionally structured proteins forms their quaternary structures, some of which have high geometric symmetry. The size and complexity of protein quaternary structures often increase in a hierarchical manner, with simpler, smaller structures serving as units for larger quaternary structures. In this study, we exploited oligomerization of a ribozyme cyclic trimer to achieve larger ribozyme-based RNA assembly. By installing kissing loop (KL) interacting units to one-, two-, or three-unit RNA molecules in the ribozyme trimer, we constructed dimers, open-chain oligomers, and branched oligomers of ribozyme trimer units. One type of open-chain oligomer preferentially formed a closed tetramer containing 12 component RNAs to provide 12 ribozyme units. We also observed large assembly of ribozyme trimers, which reached 1000 nm in size.
Collapse
Affiliation(s)
- Mst. Ayesha Siddika
- Graduate School of Innovative Life Science, University of Toyama, Toyama 930-8555, Toyama, Japan (S.M.)
| | - Hiroki Oi
- Department of Chemistry, Graduate School of Science and Engineering, University of Toyama, Toyama 930-8555, Toyama, Japan
| | - Kumi Hidaka
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8501, Kyoto, Japan
| | - Hiroshi Sugiyama
- Institute for Integrated Cell-Material Sciences, Kyoto University, Kyoto 606-8501, Kyoto, Japan; (H.S.); (M.E.)
| | - Masayuki Endo
- Institute for Integrated Cell-Material Sciences, Kyoto University, Kyoto 606-8501, Kyoto, Japan; (H.S.); (M.E.)
- Organization for Research and Development of Innovative Science and Technology, Kansai University, Suita 564-8680, Osaka, Japan
| | - Shigeyoshi Matsumura
- Graduate School of Innovative Life Science, University of Toyama, Toyama 930-8555, Toyama, Japan (S.M.)
- Department of Chemistry, Graduate School of Science and Engineering, University of Toyama, Toyama 930-8555, Toyama, Japan
| | - Yoshiya Ikawa
- Graduate School of Innovative Life Science, University of Toyama, Toyama 930-8555, Toyama, Japan (S.M.)
- Department of Chemistry, Graduate School of Science and Engineering, University of Toyama, Toyama 930-8555, Toyama, Japan
| |
Collapse
|
11
|
Kulkarni M, Thangappan J, Deb I, Wu S. Comparative analysis of RNA secondary structure accuracy on predicted RNA 3D models. PLoS One 2023; 18:e0290907. [PMID: 37656749 PMCID: PMC10473517 DOI: 10.1371/journal.pone.0290907] [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] [Received: 03/10/2023] [Accepted: 08/18/2023] [Indexed: 09/03/2023] Open
Abstract
RNA structure is conformationally dynamic, and accurate all-atom tertiary (3D) structure modeling of RNA remains challenging with the prevailing tools. Secondary structure (2D) information is the standard prerequisite for most RNA 3D modeling. Despite several 2D and 3D structure prediction tools proposed in recent years, one of the challenges is to choose the best combination for accurate RNA 3D structure prediction. Here, we benchmarked seven small RNA PDB structures (40 to 90 nucleotides) with different topologies to understand the effects of different 2D structure predictions on the accuracy of 3D modeling. The current study explores the blind challenge of 2D to 3D conversions and highlights the performances of de novo RNA 3D modeling from their predicted 2D structure constraints. Our results show that conformational sampling-based methods such as SimRNA and IsRNA1 depend less on 2D accuracy, whereas motif-based methods account for 2D evidence. Our observations illustrate the disparities in available 3D and 2D prediction methods and may further offer insights into developing topology-specific or family-specific RNA structure prediction pipelines.
Collapse
Affiliation(s)
- Mandar Kulkarni
- R&D Center, PharmCADD Co. Ltd., Dong-gu, Busan, Republic of Korea
| | | | - Indrajit Deb
- R&D Center, PharmCADD Co. Ltd., Dong-gu, Busan, Republic of Korea
| | - Sangwook Wu
- R&D Center, PharmCADD Co. Ltd., Dong-gu, Busan, Republic of Korea
- Department of Physics, Pukyong National University, Busan, Republic of Korea
| |
Collapse
|
12
|
Dupont MJ, Major F. D-ORB: A Web Server to Extract Structural Features of Related But Unaligned RNA Sequences. J Mol Biol 2023; 435:168181. [PMID: 37468182 DOI: 10.1016/j.jmb.2023.168181] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2022] [Revised: 06/02/2023] [Accepted: 06/06/2023] [Indexed: 07/21/2023]
Abstract
Identifying the common structural elements of functionally related RNA sequences (family) is usually based on an alignment of the sequences, which is often subject to human bias and may not be accurate. The resulting covariance model (CM) provides probabilities for each base to covary with another, which allows to support evolutionarily the formation of double helical regions and possibly pseudoknots. The coexistence of alternative folds in RNA, resulting from its dynamic nature, may lead to the potential omission of motifs by CM. To overcome this limitation, we present D-ORB, a system of algorithms that identifies overrepresented motifs in the secondary conformational landscapes of a family when compared to those of unrelated sequences. The algorithms are bundled into an easy-to-use website allowing users to submit a family, and optionally provide unrelated sequences. D-ORB produces a non-pseudoknotted secondary structure based on the overrepresented motifs, a deep neural network classifier and two decision trees. When used to model an Rfam family, D-ORB fits overrepresented motifs in the corresponding Rfam structure; more than a hundred Rfam families have been modeled. The statistical approach behind D-ORB derives the structural composition of an RNA family, making it a valuable tool for analyzing and modeling it. Its easy-to-use interface and advanced algorithms make it an essential resource for researchers studying RNA structure. D-ORB is available at https://d-orb.major.iric.ca/.
Collapse
Affiliation(s)
- Mathieu J Dupont
- Department of Computer Science and Operations Research, and the Institute for Research in Immunology and Cancer, Université de Montréal, Montreal, Quebec H3C 3J7, Canada
| | - François Major
- Department of Computer Science and Operations Research, and the Institute for Research in Immunology and Cancer, Université de Montréal, Montreal, Quebec H3C 3J7, Canada. https://twitter.com/francois_major
| |
Collapse
|
13
|
Yao HT, Lorenz R, Hofacker IL, Stadler PF. Mono-valent salt corrections for RNA secondary structures in the ViennaRNA package. Algorithms Mol Biol 2023; 18:8. [PMID: 37516881 PMCID: PMC10386259 DOI: 10.1186/s13015-023-00236-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2023] [Accepted: 07/10/2023] [Indexed: 07/31/2023] Open
Abstract
BACKGROUND RNA features a highly negatively charged phosphate backbone that attracts a cloud of counter-ions that reduce the electrostatic repulsion in a concentration dependent manner. Ion concentrations thus have a large influence on folding and stability of RNA structures. Despite their well-documented effects, salt effects are not handled consistently by currently available secondary structure prediction algorithms. Combining Debye-Hückel potentials for line charges and Manning's counter-ion condensation theory, Einert et al. (Biophys J 100: 2745-2753, 2011) modeled the energetic contributions of monovalent cations on loops and helices. RESULTS The model of Einert et al. is adapted to match the structure of the dynamic programming recursion of RNA secondary structure prediction algorithms. An empirical term describing the salt dependence of the duplex initiation energy is added to improve co-folding predictions for two or more RNA strands. The slightly modified model is implemented in the ViennaRNA package in such way that only the energy parameters but not the algorithmic structure is affected. A comparison with data from the literature show that predicted free energies and melting temperatures are in reasonable agreement with experiments. CONCLUSION The new feature in the ViennaRNA package makes it possible to study effects of salt concentrations on RNA folding in a systematic manner. Strictly speaking, the model pertains only to mono-valent cations, and thus covers the most important parameter, i.e., the NaCl concentration. It remains a question for future research to what extent unspecific effects of bi- and tri-valent cations can be approximated in a similar manner. AVAILABILITY Corrections for the concentration of monovalent cations are available in the ViennaRNA package starting from version 2.6.0.
Collapse
Affiliation(s)
- Hua-Ting Yao
- Department of Theoretical Chemistry, University of Vienna, Währinger Straße 17, 1090, Vienna, Austria.
| | - Ronny Lorenz
- Department of Theoretical Chemistry, University of Vienna, Währinger Straße 17, 1090, Vienna, Austria
| | - Ivo L Hofacker
- Department of Theoretical Chemistry, University of Vienna, Währinger Straße 17, 1090, Vienna, Austria
- Research Group Bioinformatics and Computational Biology, Faculty of Computer Science, University of Vienna, Währingerstraße 29, 1090, Vienna, Austria
| | - Peter F Stadler
- Department of Theoretical Chemistry, University of Vienna, Währinger Straße 17, 1090, Vienna, Austria.
- Bioinformatics Group, Department of Computer Science, and Interdisciplinary Center for Bioinformatics, Universität Leipzig, Härtelstraße 16-18, 04107, Leipzig, Germany.
- Competence Center for Scalable Data Services and Solutions Dresden/Leipzig, Interdisciplinary Center for Bioinformatics, German Centre for Integrative Biodiversity Research (iDiv), and Leipzig Research Center for Civilization Diseases, Universität Leipzig, Augustusplatz 12, 04107, Leipzig, Germany.
- Max Planck Institute for Mathematics in the Sciences, Inselstraße 22, 04109, Leipzig, Germany.
- Facultad de Ciencias, Universidad National de Colombia, Sede Bogotá, Ciudad Universitaria, 111321, Bogotá, D.C., Colombia.
- Santa Fe Institute, 1399 Hyde Park Rd., NM87501, Santa Fe, USA.
| |
Collapse
|
14
|
Mustoe AM, Weidmann CA, Weeks KM. Single-Molecule Correlated Chemical Probing: A Revolution in RNA Structure Analysis. Acc Chem Res 2023; 56:763-775. [PMID: 36917683 PMCID: PMC10078950 DOI: 10.1021/acs.accounts.2c00782] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
RNA molecules convey biological information both in their linear sequence and in their base-paired secondary and tertiary structures. Chemical probing experiments, which involve treating an RNA with a reagent that modifies conformationally dynamic nucleotides, have broadly enabled examination of short- and long-range RNA structure in diverse contexts, including in living cells. For decades, chemical probing experiments have been interpreted in a per-nucleotide way, such that the reactivity measured at each nucleotide reports the average structure at a position over all RNA molecules within a sample. However, there are numerous important cases where per-nucleotide chemical probing falls short, including for RNAs that are bound by proteins, RNAs that form complex higher order structures, and RNAs that sample multiple conformations.Recent experimental and computational innovations have started a revolution in RNA structure analysis by transforming chemical probing into a massively parallel, single-molecule experiment. Enabled by a specialized reverse transcription strategy called mutational profiling (MaP), multiple chemical modification events can be measured within individual RNA molecules. Nucleotides that communicate structurally through direct base pairing or large-scale folding-unfolding transitions will react with chemical probes in a correlated manner, thereby revealing structural complexity hidden to conventional approaches. These single-molecule correlated chemical probing (smCCP) experiments can be interpreted to directly identify nucleotides that base pair (the PAIR-MaP strategy) and to reveal long-range, through-space structural communication (RING-MaP). Correlated probing can also define the thermodynamic populations of complex RNA ensembles (DANCE-MaP). Complex RNA-protein networks can be interrogated by cross-linking proteins to RNA and measuring correlations between cross-linked positions (RNP-MaP).smCCP thus visualizes RNA secondary and higher-order structure with unprecedented accuracy, defining novel structures, RNA-protein interaction networks, time-resolved dynamics, and allosteric structural switches. These strategies are not mutually exclusive; in favorable cases, multiple levels of RNA structure ─ base pairing, through-space structural communication, and equilibrium ensembles ─ can be resolved concurrently. The physical experimentation required for smCCP is profoundly simple, and experiments are readily performed in cells on RNAs of any size, including large noncoding RNAs and mRNAs. Single-molecule correlated chemical probing is paving the way for a new generation of biophysical studies on RNA in living systems.
Collapse
Affiliation(s)
- Anthony M. Mustoe
- Verna and Marrs McClean Department of Biochemistry and Molecular Biology, Department of Molecular and Human Genetics, and Therapeutic Innovation Center (THINC), One Baylor Plaza, Baylor College of Medicine, Houston, TX 77030
| | - Chase A. Weidmann
- Department of Biological Chemistry, Center for RNA Biomedicine, 1150 W. Medical Center Drive, University of Michigan Medical School, Ann Arbor, MI 48109
| | - Kevin M. Weeks
- Department of Chemistry, University of North Carolina, Chapel Hill NC 27599-3290
| |
Collapse
|
15
|
Hollar A, Bursey H, Jabbari H. Pseudoknots in RNA Structure Prediction. Curr Protoc 2023; 3:e661. [PMID: 36779804 DOI: 10.1002/cpz1.661] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/14/2023]
Abstract
RNA molecules play active roles in the cell and are important for numerous applications in biotechnology and medicine. The function of an RNA molecule stems from its structure. RNA structure determination is time consuming, challenging, and expensive using experimental methods. Thus, much research has been directed at RNA structure prediction through computational means. Many of these methods focus primarily on the secondary structure of the molecule, ignoring the possibility of pseudoknotted structures. However, pseudoknots are known to play functional roles in many RNA molecules or in their method of interaction with other molecules. Improving the accuracy and efficiency of computational methods that predict pseudoknots is an ongoing challenge for single RNA molecules, RNA-RNA interactions, and RNA-protein interactions. To improve the accuracy of prediction, many methods focus on specific applications while restricting the length and the class of the pseudoknotted structures they can identify. In recent years, computational methods for structure prediction have begun to catch up with the impressive developments seen in biotechnology. Here, we provide a non-comprehensive overview of available pseudoknot prediction methods and their best-use cases. © 2023 Wiley Periodicals LLC.
Collapse
Affiliation(s)
- Andrew Hollar
- Department of Computer Science, University of Victoria, Victoria, Canada
| | - Hunter Bursey
- Department of Computer Science, University of Victoria, Victoria, Canada
| | - Hosna Jabbari
- Department of Computer Science, University of Victoria, Victoria, Canada
| |
Collapse
|
16
|
Jang SS, Dubnik S, Hon J, Hellenkamp B, Lynall DG, Shepard KL, Nuckolls C, Gonzalez RL. Characterizing the Conformational Free-Energy Landscape of RNA Stem-Loops Using Single-Molecule Field-Effect Transistors. J Am Chem Soc 2023; 145:402-412. [PMID: 36547391 PMCID: PMC10025942 DOI: 10.1021/jacs.2c10218] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We have developed and used single-molecule field-effect transistors (smFETs) to characterize the conformational free-energy landscape of RNA stem-loops. Stem-loops are one of the most common RNA structural motifs and serve as building blocks for the formation of complex RNA structures. Given their prevalence and integral role in RNA folding, the kinetics of stem-loop (un)folding has been extensively characterized using both experimental and computational approaches. Interestingly, these studies have reported vastly disparate timescales of (un)folding, which has been interpreted as evidence that (un)folding of even simple stem-loops occurs on a highly rugged conformational energy landscape. Because smFETs do not rely on fluorophore reporters of conformation or mechanical (un)folding forces, they provide a unique approach that has allowed us to directly monitor tens of thousands of (un)folding events of individual stem-loops at a 200 μs time resolution. Our results show that under our experimental conditions, stem-loops (un)fold over a 1-200 ms timescale during which they transition between ensembles of unfolded and folded conformations, the latter of which is composed of at least two sub-populations. The 1-200 ms timescale of (un)folding we observe here indicates that smFETs report on complete (un)folding trajectories in which unfolded conformations of the RNA spend long periods of time wandering the free-energy landscape before sampling one of several misfolded conformations or the natively folded conformation. Our findings highlight the extremely rugged landscape on which even the simplest RNA structural elements fold and demonstrate that smFETs are a unique and powerful approach for characterizing the conformational free-energy of RNA.
Collapse
Affiliation(s)
- Sukjin S. Jang
- Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 10027, USA
| | - Sarah Dubnik
- Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 10027, USA
| | - Jason Hon
- Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 10027, USA
| | - Björn Hellenkamp
- Department of Electrical Engineering, Columbia University, 3000 Broadway, New York, 10027, USA
| | - David G. Lynall
- Department of Electrical Engineering, Columbia University, 3000 Broadway, New York, 10027, USA
| | - Kenneth L. Shepard
- Department of Electrical Engineering, Columbia University, 3000 Broadway, New York, 10027, USA
| | - Colin Nuckolls
- Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 10027, USA
| | - Ruben L. Gonzalez
- Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 10027, USA
| |
Collapse
|
17
|
Mollica L, Cupaioli FA, Rossetti G, Chiappori F. An overview of structural approaches to study therapeutic RNAs. Front Mol Biosci 2022; 9:1044126. [PMID: 36387283 PMCID: PMC9649582 DOI: 10.3389/fmolb.2022.1044126] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Accepted: 10/18/2022] [Indexed: 11/07/2023] Open
Abstract
RNAs provide considerable opportunities as therapeutic agent to expand the plethora of classical therapeutic targets, from extracellular and surface proteins to intracellular nucleic acids and its regulators, in a wide range of diseases. RNA versatility can be exploited to recognize cell types, perform cell therapy, and develop new vaccine classes. Therapeutic RNAs (aptamers, antisense nucleotides, siRNA, miRNA, mRNA and CRISPR-Cas9) can modulate or induce protein expression, inhibit molecular interactions, achieve genome editing as well as exon-skipping. A common RNA thread, which makes it very promising for therapeutic applications, is its structure, flexibility, and binding specificity. Moreover, RNA displays peculiar structural plasticity compared to proteins as well as to DNA. Here we summarize the recent advances and applications of therapeutic RNAs, and the experimental and computational methods to analyze their structure, by biophysical techniques (liquid-state NMR, scattering, reactivity, and computational simulations), with a focus on dynamic and flexibility aspects and to binding analysis. This will provide insights on the currently available RNA therapeutic applications and on the best techniques to evaluate its dynamics and reactivity.
Collapse
Affiliation(s)
- Luca Mollica
- Department of Medical Biotechnologies and Translational Medicine, L.I.T.A/University of Milan, Milan, Italy
| | | | | | - Federica Chiappori
- National Research Council—Institute for Biomedical Technologies, Milan, Italy
| |
Collapse
|
18
|
Box-shaped ribozyme octamer formed by face-to-face dimerization of a pair of square-shaped ribozyme tetramers. J Biosci Bioeng 2022; 134:195-202. [PMID: 35810135 DOI: 10.1016/j.jbiosc.2022.06.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2022] [Revised: 06/14/2022] [Accepted: 06/14/2022] [Indexed: 11/21/2022]
Abstract
Naturally occurring ribozymes with defined three-dimensional (3D) structures serve as promising platforms for the design and construction of artificial RNA nanostructures. We constructed a hexameric ribozyme nanostructure by face-to-face dimerization of a pair of triangular ribozyme trimers, unit RNAs of which were derived from the Tetrahymena group I ribozyme. In this study, we have expanded the dimerization strategy to a square-shaped ribozyme tetramer by introducing four pillar units. The resulting box-shaped nanostructures, which contained eight ribozyme units, can be assembled from either four or two components of their unit RNAs.
Collapse
|
19
|
Bushhouse DZ, Choi EK, Hertz LM, Lucks JB. How does RNA fold dynamically? J Mol Biol 2022; 434:167665. [PMID: 35659535 DOI: 10.1016/j.jmb.2022.167665] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Revised: 05/26/2022] [Accepted: 05/27/2022] [Indexed: 10/18/2022]
Abstract
Recent advances in interrogating RNA folding dynamics have shown the classical model of RNA folding to be incomplete. Here, we pose three prominent questions for the field that are at the forefront of our understanding of the importance of RNA folding dynamics for RNA function. The first centers on the most appropriate biophysical framework to describe changes to the RNA folding energy landscape that a growing RNA chain encounters during transcriptional elongation. The second focuses on the potential ubiquity of strand displacement - a process by which RNA can rapidly change conformations - and how this process may be generally present in broad classes of seemingly different RNAs. The third raises questions about the potential importance and roles of cellular protein factors in RNA conformational switching. Answers to these questions will greatly improve our fundamental knowledge of RNA folding and function, drive biotechnological advances that utilize engineered RNAs, and potentially point to new areas of biology yet to be discovered.
Collapse
Affiliation(s)
- David Z Bushhouse
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, Illinois 60208, USA; Center for Synthetic Biology, Northwestern University, Evanston, Illinois 60208, USA
| | - Edric K Choi
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, Illinois 60208, USA; Center for Synthetic Biology, Northwestern University, Evanston, Illinois 60208, USA
| | - Laura M Hertz
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, Illinois 60208, USA; Center for Synthetic Biology, Northwestern University, Evanston, Illinois 60208, USA
| | - Julius B Lucks
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, Illinois 60208, USA; Center for Synthetic Biology, Northwestern University, Evanston, Illinois 60208, USA; Department of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, USA; Center for Water Research, Northwestern University, Evanston, Illinois 60208, USA; Center for Engineering Sustainability and Resilience, Northwestern University, Evanston, Illinois 60208, USA.
| |
Collapse
|
20
|
Mlýnský V, Janeček M, Kührová P, Fröhlking T, Otyepka M, Bussi G, Banáš P, Šponer J. Toward Convergence in Folding Simulations of RNA Tetraloops: Comparison of Enhanced Sampling Techniques and Effects of Force Field Modifications. J Chem Theory Comput 2022; 18:2642-2656. [PMID: 35363478 DOI: 10.1021/acs.jctc.1c01222] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Atomistic molecular dynamics simulations represent an established technique for investigation of RNA structural dynamics. Despite continuous development, contemporary RNA simulations still suffer from suboptimal accuracy of empirical potentials (force fields, ffs) and sampling limitations. Development of efficient enhanced sampling techniques is important for two reasons. First, they allow us to overcome the sampling limitations, and second, they can be used to quantify ff imbalances provided they reach a sufficient convergence. Here, we study two RNA tetraloops (TLs), namely the GAGA and UUCG motifs. We perform extensive folding simulations and calculate folding free energies (ΔGfold°) with the aim to compare different enhanced sampling techniques and to test several modifications of the nonbonded terms extending the AMBER OL3 RNA ff. We demonstrate that replica-exchange solute tempering (REST2) simulations with 12-16 replicas do not show any sign of convergence even when extended to a timescale of 120 μs per replica. However, the combination of REST2 with well-tempered metadynamics (ST-MetaD) achieves good convergence on a timescale of 5-10 μs per replica, improving the sampling efficiency by at least 2 orders of magnitude. Effects of ff modifications on ΔGfold° energies were initially explored by the reweighting approach and then validated by new simulations. We tested several manually prepared variants of the gHBfix potential which improve stability of the native state of both TLs by ∼2 kcal/mol. This is sufficient to conveniently stabilize the folded GAGA TL while the UUCG TL still remains under-stabilized. Appropriate adjustment of van der Waals parameters for C-H···O5' base-phosphate interaction may further stabilize the native states of both TLs by ∼0.6 kcal/mol.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Janeček
- Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petra Kührová
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| | - Thorben Fröhlking
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic.,IT4Innovations, VSB─Technical University of Ostrava, 17. listopadu 2172/15, 708 00 Ostrava-Poruba, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| |
Collapse
|
21
|
Cheng L, White EN, Brandt NL, Yu AM, Chen AA, Lucks J. Cotranscriptional RNA strand exchange underlies the gene regulation mechanism in a purine-sensing transcriptional riboswitch. Nucleic Acids Res 2022; 50:12001-12018. [PMID: 35348734 PMCID: PMC9756952 DOI: 10.1093/nar/gkac102] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2021] [Revised: 01/17/2022] [Accepted: 02/02/2022] [Indexed: 12/24/2022] Open
Abstract
RNA folds cotranscriptionally to traverse out-of-equilibrium intermediate structures that are important for RNA function in the context of gene regulation. To investigate this process, here we study the structure and function of the Bacillus subtilis yxjA purine riboswitch, a transcriptional riboswitch that downregulates a nucleoside transporter in response to binding guanine. Although the aptamer and expression platform domain sequences of the yxjA riboswitch do not completely overlap, we hypothesized that a strand exchange process triggers its structural switching in response to ligand binding. In vivo fluorescence assays, structural chemical probing data and experimentally informed secondary structure modeling suggest the presence of a nascent intermediate central helix. The formation of this central helix in the absence of ligand appears to compete with both the aptamer's P1 helix and the expression platform's transcriptional terminator. All-atom molecular dynamics simulations support the hypothesis that ligand binding stabilizes the aptamer P1 helix against central helix strand invasion, thus allowing the terminator to form. These results present a potential model mechanism to explain how ligand binding can induce downstream conformational changes by influencing local strand displacement processes of intermediate folds that could be at play in multiple riboswitch classes.
Collapse
Affiliation(s)
- Luyi Cheng
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, IL 60208, USA
| | - Elise N White
- Department of Biological Sciences, Rensselaer Polytechnic Institute, Troy, NY 12180, USA
| | - Naomi L Brandt
- Department of Chemistry and the RNA Institute, University at Albany, Albany, NY 12222, USA
| | - Angela M Yu
- Department of Electrical and Computer Engineering, University of Washington, Seattle, WA 98195, USA
| | - Alan A Chen
- Correspondence may also be addressed to Alan A. Chen. Tel: +1 518 437 4420;
| | - Julius B Lucks
- To whom correspondence should be addressed. Tel: +1 847 467 2943;
| |
Collapse
|
22
|
Lai LB, Lai SM, Szymanski ES, Kapur M, Choi EK, Al-Hashimi HM, Ackerman SL, Gopalan V. Structural basis for impaired 5' processing of a mutant tRNA associated with defects in neuronal homeostasis. Proc Natl Acad Sci U S A 2022; 119:e2119529119. [PMID: 35238631 PMCID: PMC8915964 DOI: 10.1073/pnas.2119529119] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2021] [Accepted: 01/09/2022] [Indexed: 02/06/2023] Open
Abstract
SignificanceUnderstanding and treating neurological disorders are global priorities. Some of these diseases are engendered by mutations that cause defects in the cellular synthesis of transfer RNAs (tRNAs), which function as adapter molecules that translate messenger RNAs into proteins. During tRNA biogenesis, ribonuclease P catalyzes removal of the transcribed sequence upstream of the mature tRNA. Here, we focus on a cytoplasmic tRNAArgUCU that is expressed specifically in neurons and, when harboring a particular point mutation, contributes to neurodegeneration in mice. Our results suggest that this mutation favors stable alternative structures that are not cleaved by mouse ribonuclease P and motivate a paradigm that may help to understand the molecular basis for disease-associated mutations in other tRNAs.
Collapse
Affiliation(s)
- Lien B. Lai
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, OH 43210
- Center for RNA Biology, The Ohio State University, Columbus, OH 43210
| | - Stella M. Lai
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, OH 43210
- Center for RNA Biology, The Ohio State University, Columbus, OH 43210
| | - Eric S. Szymanski
- Department of Biochemistry, School of Medicine, Duke University, Durham, NC 27710
| | - Mridu Kapur
- Department of Cellular and Molecular Medicine, Section of Neurobiology, University of California San Diego, La Jolla, CA 92093
- Howard Hughes Medical Institute, University of California San Diego, La Jolla, CA 92093
| | - Edric K. Choi
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, OH 43210
- Center for RNA Biology, The Ohio State University, Columbus, OH 43210
| | - Hashim M. Al-Hashimi
- Department of Biochemistry, School of Medicine, Duke University, Durham, NC 27710
| | - Susan L. Ackerman
- Department of Cellular and Molecular Medicine, Section of Neurobiology, University of California San Diego, La Jolla, CA 92093
- Howard Hughes Medical Institute, University of California San Diego, La Jolla, CA 92093
| | - Venkat Gopalan
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, OH 43210
- Center for RNA Biology, The Ohio State University, Columbus, OH 43210
| |
Collapse
|
23
|
Guo ZH, Yuan L, Tan YL, Zhang BG, Shi YZ. RNAStat: An Integrated Tool for Statistical Analysis of RNA 3D Structures. FRONTIERS IN BIOINFORMATICS 2022; 1:809082. [PMID: 36303785 PMCID: PMC9580920 DOI: 10.3389/fbinf.2021.809082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Accepted: 12/17/2021] [Indexed: 11/13/2022] Open
Abstract
The 3D architectures of RNAs are essential for understanding their cellular functions. While an accurate scoring function based on the statistics of known RNA structures is a key component for successful RNA structure prediction or evaluation, there are few tools or web servers that can be directly used to make comprehensive statistical analysis for RNA 3D structures. In this work, we developed RNAStat, an integrated tool for making statistics on RNA 3D structures. For given RNA structures, RNAStat automatically calculates RNA structural properties such as size and shape, and shows their distributions. Based on the RNA structure annotation from DSSR, RNAStat provides statistical information of RNA secondary structure motifs including canonical/non-canonical base pairs, stems, and various loops. In particular, the geometry of base-pairing/stacking can be calculated in RNAStat by constructing a local coordinate system for each base. In addition, RNAStat also supplies the distribution of distance between any atoms to the users to help build distance-based RNA statistical potentials. To test the usability of the tool, we established a non-redundant RNA 3D structure dataset, and based on the dataset, we made a comprehensive statistical analysis on RNA structures, which could have the guiding significance for RNA structure modeling. The python code of RNAStat, the dataset used in this work, and corresponding statistical data files are freely available at GitHub (https://github.com/RNA-folding-lab/RNAStat).
Collapse
Affiliation(s)
- Zhi-Hao Guo
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan, China
| | - Li Yuan
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Ben-Gong Zhang
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
- *Correspondence: Ya-Zhou Shi,
| |
Collapse
|
24
|
Wienecke A, Laederach A. A novel algorithm for ranking RNA structure candidates. Biophys J 2022; 121:7-10. [PMID: 34896370 PMCID: PMC8758412 DOI: 10.1016/j.bpj.2021.12.004] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Revised: 12/02/2021] [Accepted: 12/06/2021] [Indexed: 01/07/2023] Open
Abstract
RNA research is advancing at an ever increasing pace. The newest and most state-of-the-art instruments and techniques have made possible the discoveries of new RNAs, and they have carried the field to new frontiers of disease research, vaccine development, therapeutics, and architectonics. Like proteins, RNAs show a marked relationship between structure and function. A deeper grasp of RNAs requires a finer understanding of their elaborate structures. In pursuit of this, cutting-edge experimental and computational structure-probing techniques output several candidate geometries for a given RNA, each of which is perfectly aligned with experimentally determined parameters. Identifying which structure is the most accurate, however, remains a major obstacle. In recent years, several algorithms have been developed for ranking candidate RNA structures in order from most to least probable, though their levels of accuracy and transparency leave room for improvement. Most recently, advances in both areas are demonstrated by rsRNASP, a novel algorithm proposed by Tan et al. rsRNASP is a residue-separation-based statistical potential for three-dimensional structure evaluation, and it outperforms the leading algorithms in the field.
Collapse
Affiliation(s)
- Anastacia Wienecke
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina,Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina
| | - Alain Laederach
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina,Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina,Corresponding author
| |
Collapse
|
25
|
rsRNASP: A residue-separation-based statistical potential for RNA 3D structure evaluation. Biophys J 2022; 121:142-156. [PMID: 34798137 PMCID: PMC8758408 DOI: 10.1016/j.bpj.2021.11.016] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 10/23/2021] [Accepted: 11/10/2021] [Indexed: 01/07/2023] Open
Abstract
Knowledge-based statistical potentials have been shown to be rather effective in protein 3-dimensional (3D) structure evaluation and prediction. Recently, several statistical potentials have been developed for RNA 3D structure evaluation, while their performances are either still at a low level for the test datasets from structure prediction models or dependent on the "black-box" process through neural networks. In this work, we have developed an all-atom distance-dependent statistical potential based on residue separation for RNA 3D structure evaluation, namely rsRNASP, which is composed of short- and long-ranged potentials distinguished by residue separation. The extensive examinations against available RNA test datasets show that rsRNASP has apparently higher performance than the existing statistical potentials for the realistic test datasets with large RNAs from structure prediction models, including the newly released RNA-Puzzles dataset, and is comparable to the existing top statistical potentials for the test datasets with small RNAs or near-native decoys. In addition, rsRNASP is superior to RNA3DCNN, a recently developed scoring function through 3D convolutional neural networks. rsRNASP and the relevant databases are available to the public.
Collapse
|
26
|
Fang X, Gallego J, Wang YX. Deriving RNA topological structure from SAXS. Methods Enzymol 2022; 677:479-529. [DOI: 10.1016/bs.mie.2022.08.037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
27
|
Zerze GH, Piaggi PM, Debenedetti PG. A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States. J Phys Chem B 2021; 125:13685-13695. [PMID: 34890201 DOI: 10.1021/acs.jpcb.1c08038] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. We isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.
Collapse
Affiliation(s)
- Gül H Zerze
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Pablo M Piaggi
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
28
|
Poblete S, Božič A, Kanduč M, Podgornik R, Guzman HV. RNA Secondary Structures Regulate Adsorption of Fragments onto Flat Substrates. ACS OMEGA 2021; 6:32823-32831. [PMID: 34901632 PMCID: PMC8655909 DOI: 10.1021/acsomega.1c04774] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Accepted: 11/04/2021] [Indexed: 06/14/2023]
Abstract
RNA is a functionally rich molecule with multilevel, hierarchical structures whose role in the adsorption to molecular substrates is only beginning to be elucidated. Here, we introduce a multiscale simulation approach that combines a tractable coarse-grained RNA structural model with an interaction potential of a structureless flat adsorbing substrate. Within this approach, we study the specific role of stem-hairpin and multibranch RNA secondary structure motifs on its adsorption phenomenology. Our findings identify a dual regime of adsorption for short RNA fragments with and without the secondary structure and underline the adsorption efficiency in both cases as a function of the surface interaction strength. The observed behavior results from an interplay between the number of contacts formed at the surface and the conformational entropy of the RNA molecule. The adsorption phenomenology of RNA seems to persist also for much longer RNAs as qualitatively observed by comparing the trends of our simulations with a theoretical approach based on an ideal semiflexible polymer chain.
Collapse
Affiliation(s)
- Simón Poblete
- Instituto
de Ciencias Físicas y Matemáticas, Universidad Austral de Chile, Valdivia 5091000, Chile
- Computational
Biology Lab, Fundación Ciencia &
Vida, Santiago 7780272, Chile
| | - Anže Božič
- Department
of Theoretical Physics, Jožef Stefan
Institute, SI-1000 Ljubljana, Slovenia
| | - Matej Kanduč
- Department
of Theoretical Physics, Jožef Stefan
Institute, SI-1000 Ljubljana, Slovenia
| | - Rudolf Podgornik
- School
of Physical Sciences and Kavli Institute for Theoretical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
- Institute
of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Wenzhou
Institute of the University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
- Department
of Physics, Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia
| | - Horacio V. Guzman
- Department
of Theoretical Physics, Jožef Stefan
Institute, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
29
|
Abstract
To exert their functions, RNAs adopt diverse structures, ranging from simple secondary to complex tertiary and quaternary folds. In vivo, RNA folding starts with RNA transcription, and a wide variety of processes are coupled to co-transcriptional RNA folding events, including the regulation of fundamental transcription dynamics, gene regulation by mechanisms like attenuation, RNA processing or ribonucleoprotein particle formation. While co-transcriptional RNA folding and associated co-transcriptional processes are by now well accepted as pervasive regulatory principles in all organisms, investigations into the role of the transcription machinery in co-transcriptional folding processes have so far largely focused on effects of the order in which RNA regions are produced and of transcription kinetics. Recent structural and structure-guided functional analyses of bacterial transcription complexes increasingly point to an additional role of RNA polymerase and associated transcription factors in supporting co-transcriptional RNA folding by fostering or preventing strategic contacts to the nascent transcripts. In general, the results support the view that transcription complexes can act as RNA chaperones, a function that has been suggested over 30 years ago. Here, we discuss transcription complexes as RNA chaperones based on recent examples from bacterial transcription.
Collapse
Affiliation(s)
- Nelly Said
- Freie Universität Berlin, Department Biology, Chemistry, Pharmacy, Institute of Chemistry and Biochemistry, Laboratory of Structural Biochemistry, Berlin, Germany
| | - Markus C Wahl
- Freie Universität Berlin, Department Biology, Chemistry, Pharmacy, Institute of Chemistry and Biochemistry, Laboratory of Structural Biochemistry, Berlin, Germany.,Helmholtz-Zentrum Berlin Für Materialien Und Energie, Macromolecular Crystallography, Berlin, Germany
| |
Collapse
|
30
|
Manz C, Kobitski AY, Samanta A, Nienhaus K, Jäschke A, Nienhaus GU. Exploring the energy landscape of a SAM-I riboswitch. J Biol Phys 2021; 47:371-386. [PMID: 34698957 PMCID: PMC8603990 DOI: 10.1007/s10867-021-09584-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2021] [Accepted: 08/25/2021] [Indexed: 11/30/2022] Open
Abstract
SAM-I riboswitches regulate gene expression through transcription termination upon binding a S-adenosyl-L-methionine (SAM) ligand. In previous work, we characterized the conformational energy landscape of the full-length Bacillus subtilis yitJ SAM-I riboswitch as a function of Mg2+ and SAM ligand concentrations. Here, we have extended this work with measurements on a structurally similar ligand, S-adenosyl-l-homocysteine (SAH), which has, however, a much lower binding affinity. Using single-molecule Förster resonance energy transfer (smFRET) microscopy and hidden Markov modeling (HMM) analysis, we identified major conformations and determined their fractional populations and dynamics. At high Mg2+ concentration, FRET analysis yielded four distinct conformations, which we assigned to two terminator and two antiterminator states. In the same solvent, but with SAM added at saturating concentrations, four states persisted, although their populations, lifetimes and interconversion dynamics changed. In the presence of SAH instead of SAM, HMM revealed again four well-populated states and, in addition, a weakly populated ‘hub’ state that appears to mediate conformational transitions between three of the other states. Our data show pronounced and specific effects of the SAM and SAH ligands on the RNA conformational energy landscape. Interestingly, both SAM and SAH shifted the fractional populations toward terminator folds, but only gradually, so the effect cannot explain the switching action. Instead, we propose that the noticeably accelerated dynamics of interconversion between terminator and antiterminator states upon SAM binding may be essential for control of transcription.
Collapse
Affiliation(s)
- Christoph Manz
- Institute of Applied Physics, Karlsruhe Institute of Technology, Wolfgang-Gaede-Str. 1, 76131, Karlsruhe, Germany
| | - Andrei Yu Kobitski
- Institute of Applied Physics, Karlsruhe Institute of Technology, Wolfgang-Gaede-Str. 1, 76131, Karlsruhe, Germany
| | - Ayan Samanta
- Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Im Neuenheimer Feld 364, 69120, Heidelberg, Germany.,Department of Chemistry, Uppsala University, Box 538, 751 21, Uppsala, Sweden
| | - Karin Nienhaus
- Institute of Applied Physics, Karlsruhe Institute of Technology, Wolfgang-Gaede-Str. 1, 76131, Karlsruhe, Germany
| | - Andres Jäschke
- Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Im Neuenheimer Feld 364, 69120, Heidelberg, Germany
| | - Gerd Ulrich Nienhaus
- Institute of Applied Physics, Karlsruhe Institute of Technology, Wolfgang-Gaede-Str. 1, 76131, Karlsruhe, Germany. .,Institute of Nanotechnology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany. .,Institute of Biological and Chemical Systems, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany. .,Department of Physics, University of Illinois at Urbana-Champaign, 1110 W. Green St, Urbana, IL, 61801, USA.
| |
Collapse
|
31
|
Predicting RNA Secondary Structure Using In Vitro and In Vivo Data. Methods Mol Biol 2021. [PMID: 34694602 DOI: 10.1007/978-1-0716-1851-6_2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
Abstract
The new flow of high-throughput RNA secondary structure data coming from different techniques allowed the further development of machine learning approaches. We developed CROSS and CROSSalive, two algorithms trained on experimental data able to predict the RNA secondary structure propensity both in vitro and in vivo. Since the in vivo folding of RNA molecules depends on multiple factors due to the cellular crowded environment, prediction is a complex problem that needs additional calculations for the interaction with proteins and other molecules. In the following chapter, we will describe the differences in predicting RNA secondary structure propensity using experimental data as input for an Artificial Neural Network (ANN) in vitro and in vivo.
Collapse
|
32
|
Gunawardhana SM, Holmstrom ED. Apolar chemical environments compact unfolded RNAs and can promote folding. BIOPHYSICAL REPORTS 2021; 1. [PMID: 35382036 PMCID: PMC8978554 DOI: 10.1016/j.bpr.2021.100004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
It is well documented that the structure, and thus function, of nucleic acids depends on the chemical environment surrounding them, which often includes potential proteinaceous binding partners. The nonpolar amino acid side chains of these proteins will invariably alter the polarity of the local chemical environment around the nucleic acid. However, we are only beginning to understand how environmental polarity generally influences the structural and energetic properties of RNA folding. Here, we use a series of aqueous-organic cosolvent mixtures to systematically modulate the solvent polarity around two different RNA folding constructs that can form either secondary or tertiary structural elements. Using single-molecule Förster resonance energy transfer spectroscopy to simultaneously monitor the structural and energetic properties of these RNAs, we show that the unfolded conformations of both model RNAs become more compact in apolar environments characterized by dielectric constants less than that of pure water. In the case of tertiary structure formation, this compaction also gives rise to more energetically favorable folding. We propose that these physical changes arise from an enhanced accumulation of counterions in the low dielectric environment surrounding the unfolded RNA.
Collapse
Affiliation(s)
| | - Erik D Holmstrom
- Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas.,Department of Chemistry, University of Kansas, Lawrence, Kansas
| |
Collapse
|
33
|
High-throughput dissection of the thermodynamic and conformational properties of a ubiquitous class of RNA tertiary contact motifs. Proc Natl Acad Sci U S A 2021; 118:2109085118. [PMID: 34373334 DOI: 10.1073/pnas.2109085118] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Despite RNA's diverse secondary and tertiary structures and its complex conformational changes, nature utilizes a limited set of structural "motifs"-helices, junctions, and tertiary contact modules-to build diverse functional RNAs. Thus, in-depth descriptions of a relatively small universe of RNA motifs may lead to predictive models of RNA tertiary conformational landscapes. Motifs may have different properties depending on sequence and secondary structure, giving rise to subclasses that expand the universe of RNA building blocks. Yet we know very little about motif subclasses, given the challenges in mapping conformational properties in high throughput. Previously, we used "RNA on a massively parallel array" (RNA-MaP), a quantitative, high-throughput technique, to study thousands of helices and two-way junctions. Here, we adapt RNA-MaP to study the thermodynamic and conformational properties of tetraloop/tetraloop receptor (TL/TLR) tertiary contact motifs, analyzing 1,493 TLR sequences from different classes. Clustering analyses revealed variability in TL specificity, stability, and conformational behavior. Nevertheless, natural GAAA/11ntR TL/TLRs, while varying in tertiary stability by ∼2.5 kcal/mol, exhibited conserved TL specificity and conformational properties. Thus, RNAs may tune stability without altering the overall structure of these TL/TLRs. Furthermore, their stability correlated with natural frequency, suggesting thermodynamics as the dominant selection pressure. In contrast, other TL/TLRs displayed heterogenous conformational behavior and appear to not be under strong thermodynamic selection. Our results build toward a generalizable model of RNA-folding thermodynamics based on the properties of isolated motifs, and our characterized TL/TLR library can be used to engineer RNAs with predictable thermodynamic and conformational behavior.
Collapse
|
34
|
Poblete S, Guzman HV. Structural 3D Domain Reconstruction of the RNA Genome from Viruses with Secondary Structure Models. Viruses 2021; 13:1555. [PMID: 34452420 PMCID: PMC8402887 DOI: 10.3390/v13081555] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2021] [Revised: 07/28/2021] [Accepted: 07/30/2021] [Indexed: 02/07/2023] Open
Abstract
Three-dimensional RNA domain reconstruction is important for the assembly, disassembly and delivery functionalities of a packed proteinaceus capsid. However, to date, the self-association of RNA molecules is still an open problem. Recent chemical probing reports provide, with high reliability, the secondary structure of diverse RNA ensembles, such as those of viral genomes. Here, we present a method for reconstructing the complete 3D structure of RNA genomes, which combines a coarse-grained model with a subdomain composition scheme to obtain the entire genome inside proteinaceus capsids based on secondary structures from experimental techniques. Despite the amount of sampling involved in the folded and also unfolded RNA molecules, advanced microscope techniques can provide points of anchoring, which enhance our model to include interactions between capsid pentamers and RNA subdomains. To test our method, we tackle the satellite tobacco mosaic virus (STMV) genome, which has been widely studied by both experimental and computational communities. We provide not only a methodology to structurally analyze the tertiary conformations of the RNA genome inside capsids, but a flexible platform that allows the easy implementation of features/descriptors coming from both theoretical and experimental approaches.
Collapse
Affiliation(s)
- Simón Poblete
- Instituto de Ciencias Físicas y Matemáticas, Universidad Austral de Chile, Valdivia 5091000, Chile
- Chile and Computational Biology Lab, Fundación Ciencia & Vida, Santiago 7780272, Chile
| | - Horacio V. Guzman
- Department of Theoretical Physics, Jožef Stefan Institute, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
35
|
Hsu CF, Chang KC, Chen YL, Hsieh PS, Lee AI, Tu JY, Chen YT, Wen JD. Formation of frameshift-stimulating RNA pseudoknots is facilitated by remodeling of their folding intermediates. Nucleic Acids Res 2021; 49:6941-6957. [PMID: 34161580 PMCID: PMC8266650 DOI: 10.1093/nar/gkab512] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Revised: 05/27/2021] [Accepted: 06/04/2021] [Indexed: 12/15/2022] Open
Abstract
Programmed –1 ribosomal frameshifting is an essential regulation mechanism of translation in viruses and bacteria. It is stimulated by mRNA structures inside the coding region. As the structure is unfolded repeatedly by consecutive translating ribosomes, whether it can refold properly each time is important in performing its function. By using single-molecule approaches and molecular dynamics simulations, we found that a frameshift-stimulating RNA pseudoknot folds sequentially through its upstream stem S1 and downstream stem S2. In this pathway, S2 folds from the downstream side and tends to be trapped in intermediates. By masking the last few nucleotides to mimic their gradual emergence from translating ribosomes, S2 can be directed to fold from the upstream region. The results show that the intermediates are greatly suppressed, suggesting that mRNA refolding may be modulated by ribosomes. Moreover, masking the first few nucleotides of S1 favors the folding from S2 and yields native pseudoknots, which are stable enough to retrieve the masked nucleotides. We hypothesize that translating ribosomes can remodel an intermediate mRNA structure into a stable conformation, which may in turn stimulate backward slippage of the ribosome. This supports an interactive model of ribosomal frameshifting and gives an insightful account addressing previous experimental observations.
Collapse
Affiliation(s)
- Chiung-Fang Hsu
- Institute of Molecular and Cellular Biology, National Taiwan University, Taipei 10617, Taiwan
| | - Kai-Chun Chang
- Institute of Molecular and Cellular Biology, National Taiwan University, Taipei 10617, Taiwan
| | - Yi-Lan Chen
- Genome and Systems Biology Degree Program, National Taiwan University and Academia Sinica, Taipei 10617, Taiwan
| | - Po-Szu Hsieh
- Institute of Molecular and Cellular Biology, National Taiwan University, Taipei 10617, Taiwan
| | - An-I Lee
- Institute of Molecular and Cellular Biology, National Taiwan University, Taipei 10617, Taiwan
| | - Jui-Yun Tu
- Institute of Molecular and Cellular Biology, National Taiwan University, Taipei 10617, Taiwan
| | - Yu-Ting Chen
- Institute of Molecular and Cellular Biology, National Taiwan University, Taipei 10617, Taiwan
| | - Jin-Der Wen
- Institute of Molecular and Cellular Biology, National Taiwan University, Taipei 10617, Taiwan.,Genome and Systems Biology Degree Program, National Taiwan University and Academia Sinica, Taipei 10617, Taiwan
| |
Collapse
|
36
|
Christy TW, Giannetti CA, Houlihan G, Smola MJ, Rice GM, Wang J, Dokholyan NV, Laederach A, Holliger P, Weeks KM. Direct Mapping of Higher-Order RNA Interactions by SHAPE-JuMP. Biochemistry 2021; 60:1971-1982. [PMID: 34121404 PMCID: PMC8256721 DOI: 10.1021/acs.biochem.1c00270] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Higher-order structure governs function for many RNAs. However, discerning this structure for large RNA molecules in solution is an unresolved challenge. Here, we present SHAPE-JuMP (selective 2'-hydroxyl acylation analyzed by primer extension and juxtaposed merged pairs) to interrogate through-space RNA tertiary interactions. A bifunctional small molecule is used to chemically link proximal nucleotides in an RNA structure. The RNA cross-link site is then encoded into complementary DNA (cDNA) in a single, direct step using an engineered reverse transcriptase that "jumps" across cross-linked nucleotides. The resulting cDNAs contain a deletion relative to the native RNA sequence, which can be detected by sequencing, that indicates the sites of cross-linked nucleotides. SHAPE-JuMP measures RNA tertiary structure proximity concisely across large RNA molecules at nanometer resolution. SHAPE-JuMP is especially effective at measuring interactions in multihelix junctions and loop-to-helix packing, enables modeling of the global fold for RNAs up to several hundred nucleotides in length, facilitates ranking of structural models by consistency with through-space restraints, and is poised to enable solution-phase structural interrogation and modeling of complex RNAs.
Collapse
Affiliation(s)
- Thomas W. Christy
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina, Chapel Hill, North Carolina 27599
| | - Catherine A. Giannetti
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290
| | - Gillian Houlihan
- MRC Laboratory of Molecular Biology, Francis Crick Avenue, Cambridge CB2 0QH, UK
| | - Matthew J. Smola
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290
| | - Greggory M. Rice
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290
| | - Jian Wang
- Departments of Pharmacology, and Biochemistry and Molecular Biology, Penn State University College of Medicine, Hershey, PA 17033, USA
| | - Nikolay V. Dokholyan
- Departments of Pharmacology, and Biochemistry and Molecular Biology, Penn State University College of Medicine, Hershey, PA 17033, USA
- Departments of Chemistry, and Biomedical Engineering, Pennsylvania State University, University Park, PA 16802
| | - Alain Laederach
- Department of Biology, University of North Carolina, Chapel Hill, North Carolina 27599
| | - Philipp Holliger
- MRC Laboratory of Molecular Biology, Francis Crick Avenue, Cambridge CB2 0QH, UK
| | - Kevin M. Weeks
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290
| |
Collapse
|
37
|
Chen J, Chen X, Sun LZ, Xu XJ, Luo MB. Translocation of a looped polymer threading through a nanopore. SOFT MATTER 2021; 17:4342-4351. [PMID: 33908563 DOI: 10.1039/d1sm00007a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Recent experiments reported that the complicated translocation dynamics of a looped DNA chain through a nanopore can be detected by ionic current blockade profiles. Inspired by the experimental results, we systematically study the translocation dynamics of a looped polymer, formed by three building blocks of a loop in the middle and two tails of the same length connected with the loop, by using Langevin dynamics simulations. Based on two entering modes (tail-leading and loop-leading) and three translocation orders (loop-tail-tail, tail-loop-tail, and tail-tail-loop), the translocation of the looped polymer is classified into six translocation pathways, corresponding to different current blockade profiles. The probabilities of the six translocation pathways are dependent on the loop length, polymer length, and pore radius. Moreover, the translocation times of the entire polymer and the loop are investigated. We find that the two translocation times show different dependencies on the translocation pathways and on the lengths of the loop and the entire polymer.
Collapse
Affiliation(s)
- Jia Chen
- Department of Physics, Zhejiang University, Hangzhou 310027, China.
| | - Xian Chen
- Department of Physics, Zhejiang University, Hangzhou 310027, China.
| | - Li-Zhen Sun
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China.
| | - Xiao-Jun Xu
- Institute of Bioinformatics and Medical Engineering, Jiangsu University of Technology, Changzhou, 213001, China
| | - Meng-Bo Luo
- Department of Physics, Zhejiang University, Hangzhou 310027, China.
| |
Collapse
|
38
|
Sung HL, Sengupta A, Nesbitt D. Smaller molecules crowd better: Crowder size dependence revealed by single-molecule FRET studies and depletion force modeling analysis. J Chem Phys 2021; 154:155101. [PMID: 33887926 DOI: 10.1063/5.0045492] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The cell is an extremely crowded environment, which is known to have a profound impact on the thermodynamics, functionality, and conformational stability of biomolecules. Speculations from recent theoretical molecular dynamics studies suggest an intriguing size dependence to such purely entropic crowding effects, whereby small molecular weight crowders under constant enthalpy conditions are more effective than larger crowders on a per volume basis. If experimentally confirmed, this would be profoundly significant, as the cellular cytoplasm is also quite concentrated in smaller molecular weight solutes such as inorganic ions, amino acids, and various metabolites. The challenge is to perform such studies isolating entropic effects under isoenthalpic conditions. In this work, we first present results from single-molecule FRET spectroscopy (smFRET) on the molecular size-dependent crowding stabilization of a simple RNA tertiary motif (the GAAA tetraloop-tetraloop receptor), indeed providing evidence in support of the surprising notion in the crowding literature that "smaller is better." Specifically, systematic smFRET studies as a function of crowder solute size reveal that smaller molecules both significantly increase the RNA tertiary folding rate and, yet, simultaneously decrease the unfolding rate, predicting strongly size-dependent stabilization of RNA tertiary structures under crowded cellular conditions. The size dependence of these effects has been explored via systematic variation of crowder size over a broad range of molecular weights (90-3000 amu). Furthermore, corresponding temperature dependent studies indicate the systematic changes in the folding equilibrium to be predominantly entropic in origin, i.e., consistent with a fundamental picture of entropic molecular crowding without additional enthalpic interactions. Most importantly, all trends in the single-molecule crowding data can be quantitatively recapitulated by a simple analytic depletion force model, whereby excluded volume interactions represent the major thermodynamic driving force toward folding. Our study, thus, not only provides experimental evidence and theoretical support for small molecule crowding but also predicts further enhancement of crowding effects for even smaller molecules on a per volume basis.
Collapse
Affiliation(s)
- Hsuan-Lei Sung
- JILA, National Institute of Standards and Technology and University of Colorado, Boulder, Colorado 80309, USA
| | - Abhigyan Sengupta
- Biophysics Department, Technical University of Munich, Garching, Germany
| | - David Nesbitt
- JILA, National Institute of Standards and Technology and University of Colorado, Boulder, Colorado 80309, USA
| |
Collapse
|
39
|
Sun L, Li P, Ju X, Rao J, Huang W, Ren L, Zhang S, Xiong T, Xu K, Zhou X, Gong M, Miska E, Ding Q, Wang J, Zhang QC. In vivo structural characterization of the SARS-CoV-2 RNA genome identifies host proteins vulnerable to repurposed drugs. Cell 2021; 184:1865-1883.e20. [PMID: 33636127 PMCID: PMC7871767 DOI: 10.1016/j.cell.2021.02.008] [Citation(s) in RCA: 129] [Impact Index Per Article: 43.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2020] [Revised: 10/01/2020] [Accepted: 02/02/2021] [Indexed: 01/10/2023]
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the ongoing coronavirus disease 2019 (COVID-19) pandemic. Understanding of the RNA virus and its interactions with host proteins could improve therapeutic interventions for COVID-19. By using icSHAPE, we determined the structural landscape of SARS-CoV-2 RNA in infected human cells and from refolded RNAs, as well as the regulatory untranslated regions of SARS-CoV-2 and six other coronaviruses. We validated several structural elements predicted in silico and discovered structural features that affect the translation and abundance of subgenomic viral RNAs in cells. The structural data informed a deep-learning tool to predict 42 host proteins that bind to SARS-CoV-2 RNA. Strikingly, antisense oligonucleotides targeting the structural elements and FDA-approved drugs inhibiting the SARS-CoV-2 RNA binding proteins dramatically reduced SARS-CoV-2 infection in cells derived from human liver and lung tumors. Our findings thus shed light on coronavirus and reveal multiple candidate therapeutics for COVID-19 treatment.
Collapse
Affiliation(s)
- Lei Sun
- MOE Key Laboratory of Bioinformatics, Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, Center for Synthetic and Systems Biology, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, Beijing 100084, China
| | - Pan Li
- MOE Key Laboratory of Bioinformatics, Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, Center for Synthetic and Systems Biology, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, Beijing 100084, China
| | - Xiaohui Ju
- Center for Infectious Disease Research, School of Medicine, Tsinghua University, Beijing 100084, China
| | - Jian Rao
- NHC Key Laboratory of Systems Biology of Pathogens and Christophe Mérieux Laboratory, Institute of Pathogen Biology, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100730, China
| | - Wenze Huang
- MOE Key Laboratory of Bioinformatics, Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, Center for Synthetic and Systems Biology, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, Beijing 100084, China
| | - Lili Ren
- NHC Key Laboratory of Systems Biology of Pathogens and Christophe Mérieux Laboratory, Institute of Pathogen Biology, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100730, China; Key Laboratory of Respiratory Disease Pathogenomics, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing 100730, China
| | - Shaojun Zhang
- MOE Key Laboratory of Bioinformatics, Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, Center for Synthetic and Systems Biology, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, Beijing 100084, China
| | - Tuanlin Xiong
- MOE Key Laboratory of Bioinformatics, Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, Center for Synthetic and Systems Biology, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, Beijing 100084, China
| | - Kui Xu
- MOE Key Laboratory of Bioinformatics, Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, Center for Synthetic and Systems Biology, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, Beijing 100084, China
| | - Xiaolin Zhou
- MOE Key Laboratory of Bioinformatics, Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, Center for Synthetic and Systems Biology, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, Beijing 100084, China
| | - Mingli Gong
- Center for Infectious Disease Research, School of Medicine, Tsinghua University, Beijing 100084, China
| | - Eric Miska
- Wellcome Trust/Cancer Research UK Gurdon Institute, Department of Genetics, University of Cambridge, Cambridge CB2 1QN, UK; Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge, CB10 1SA, UK
| | - Qiang Ding
- Center for Infectious Disease Research, School of Medicine, Tsinghua University, Beijing 100084, China.
| | - Jianwei Wang
- NHC Key Laboratory of Systems Biology of Pathogens and Christophe Mérieux Laboratory, Institute of Pathogen Biology, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100730, China; Key Laboratory of Respiratory Disease Pathogenomics, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing 100730, China.
| | - Qiangfeng Cliff Zhang
- MOE Key Laboratory of Bioinformatics, Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, Center for Synthetic and Systems Biology, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, Beijing 100084, China.
| |
Collapse
|
40
|
Zhang D, Li J, Chen SJ. IsRNA1: De Novo Prediction and Blind Screening of RNA 3D Structures. J Chem Theory Comput 2021; 17:1842-1857. [PMID: 33560836 DOI: 10.1021/acs.jctc.0c01148] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Modeling structures and functions of large ribonucleic acid (RNAs) especially with complicated topologies is highly challenging due to the inefficiency of large conformational sampling and the presence of complicated tertiary interactions. To address this problem, one highly promising approach is coarse-grained modeling. Here, following an iterative simulated reference state approach to decipher the correlations between different structural parameters, we developed a potent coarse-grained RNA model named as IsRNA1 for RNA studies. Molecular dynamics simulations in the IsRNA1 can predict the native structures of small RNAs from a sequence and fold medium-sized RNAs into near-native tertiary structures with the assistance of secondary structure constraints. A large-scale benchmark test on RNA 3D structure prediction shows that IsRNA1 exhibits improved performance for relatively large RNAs of complicated topologies, such as large stem-loop structures and structures containing long-range tertiary interactions. The advantages of IsRNA1 include the consideration of the correlations between the different structural variables, the appropriate characterization of canonical base-pairing and base-stacking interactions, and the better sampling for the backbone conformations. Moreover, a blind screening protocol was developed based on IsRNA1 to identify good structural models from a pool of candidates without prior knowledge of the native structures.
Collapse
Affiliation(s)
- Dong Zhang
- Department of Physics, Department of Biochemistry, and Institute of Data Science and Informatics, University of Missouri, Columbia, Missouri 65211, United States
| | - Jun Li
- Department of Physics, Department of Biochemistry, and Institute of Data Science and Informatics, University of Missouri, Columbia, Missouri 65211, United States
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Institute of Data Science and Informatics, University of Missouri, Columbia, Missouri 65211, United States
| |
Collapse
|
41
|
Despons L, Martin F. How Many Messenger RNAs Can Be Translated by the START Mechanism? Int J Mol Sci 2020; 21:ijms21218373. [PMID: 33171614 PMCID: PMC7664666 DOI: 10.3390/ijms21218373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2020] [Revised: 10/29/2020] [Accepted: 11/05/2020] [Indexed: 11/22/2022] Open
Abstract
Translation initiation is a key step in the protein synthesis stage of the gene expression pathway of all living cells. In this important process, ribosomes have to accurately find the AUG start codon in order to ensure the integrity of the proteome. “Structure Assisted RNA Translation”, or “START”, has been proposed to use stable secondary structures located in the coding sequence to augment start site selection by steric hindrance of the progression of pre-initiation complex on messenger RNA. This implies that such structures have to be located downstream and at on optimal distance from the AUG start codon (i.e., downstream nucleotide +16). In order to assess the importance of the START mechanism in the overall mRNA translation process, we developed a bioinformatic tool to screen coding sequences for such stable structures in a 50 nucleotide-long window spanning the nucleotides from +16 to +65. We screened eight bacterial genomes and six eukaryotic genomes. We found stable structures in 0.6–2.5% of eukaryotic coding sequences. Among these, approximately half of them were structures predicted to form G-quadruplex structures. In humans, we selected 747 structures. In bacteria, the coding sequences from Gram-positive bacteria contained 2.6–4.2% stable structures, whereas the structures were less abundant in Gram-negative bacteria (0.2–2.7%). In contrast to eukaryotes, putative G-quadruplex structures are very rare in the coding sequence of bacteria. Altogether, our study reveals that the START mechanism seems to be an ancient strategy to facilitate the start codon recognition that is used in different kingdoms of life.
Collapse
|
42
|
Sun LZ, Jiang Y, Zhou Y, Chen SJ. RLDOCK: A New Method for Predicting RNA-Ligand Interactions. J Chem Theory Comput 2020; 16:7173-7183. [PMID: 33095555 DOI: 10.1021/acs.jctc.0c00798] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
The ability to accurately predict the binding site, binding pose, and binding affinity for ligand-RNA binding is important for RNA-targeted drug design. Here, we describe a new computational method, RLDOCK, for predicting the binding site and binding pose for ligand-RNA binding. By developing an energy-based scoring function, we sample exhaustively all of the possible binding sites with flexible ligand conformations for a ligand-RNA pair based on the geometric and energetic scores. The model distinguishes from other approaches in three notable features. First, the model enables exhaustive scanning of all of the possible binding sites, including multiple alternative or coexisting binding sites, for a given ligand-RNA pair. Second, the model is based on a new energy-based scoring function developed here. Third, the model employs a novel multistep screening algorithm to improve computational efficiency. Specifically, first, for each binding site, we use a gird-based energy map to rank the binding sites according to the minimum Lennard-Jones potential energy for the different ligand poses. Second, for a given selected binding site, we predict the ligand pose using a two-step algorithm. In the first step, we quickly identify the probable ligand poses using a coarse-grained simplified energy function. In the second step, for each of the probable ligand poses, we predict the ligand poses using a refined energy function. Tests of the RLDOCK for a set of 230 RNA-ligand-bound structures indicate that RLDOCK can successfully predict ligand poses for 27.8, 58.3, and 69.6% of all of the test cases with the root-mean-square deviation within 1.0, 2.0, and 3.0 Å, respectively, for the top three predicted docking poses. The computational method presented here may enable the development of a new, more comprehensive framework for the prediction of ligand-RNA binding with an ensemble of RNA conformations and the metal-ion effects.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China.,Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, United States
| | - Yangwei Jiang
- Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, United States
| | - Yuanzhe Zhou
- Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, United States
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, United States
| |
Collapse
|
43
|
Ehrhardt JE, Weeks KM. Time-Resolved, Single-Molecule, Correlated Chemical Probing of RNA. J Am Chem Soc 2020; 142:18735-18740. [PMID: 33095984 DOI: 10.1021/jacs.0c06221] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Capturing the folding dynamics of large, functionally important RNAs has relied primarily on global measurements of structure or on per-nucleotide chemical probing. These approaches infer, but do not directly measure, through-space structural interactions. Here we introduce trimethyloxonium (TMO) as a chemical probe for RNA. TMO alkylates RNA at high levels in seconds, and thereby enables time-resolved, single-molecule, through-space probing of RNA folding using the RING-MaP correlated chemical probing framework. Time-resolved correlations in the RNase P RNA-a functional RNA with a complex structure stabilized by multiple noncanonical interactions-revealed that a long-range tertiary interaction guides native RNA folding for both secondary and tertiary structure. This unanticipated nonhierarchical folding mechanism was directly validated by examining the consequences of concise disruption of the through-space interaction. Single-molecule, time-resolved RNA structure probing with TMO is poised to reveal a wide range of dynamic RNA folding processes and principles of RNA folding.
Collapse
Affiliation(s)
- Jeffrey E Ehrhardt
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| | - Kevin M Weeks
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| |
Collapse
|
44
|
Zhao C, Zhang D, Jiang Y, Chen SJ. Modeling Loop Composition and Ion Concentration Effects in RNA Hairpin Folding Stability. Biophys J 2020; 119:1439-1455. [PMID: 32949490 PMCID: PMC7568001 DOI: 10.1016/j.bpj.2020.07.042] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Revised: 06/12/2020] [Accepted: 07/08/2020] [Indexed: 12/21/2022] Open
Abstract
The ability to accurately predict RNA hairpin structure and stability for different loop sequences and salt conditions is important for understanding, modeling, and designing larger RNA folds. However, traditional RNA secondary structure models cannot treat loop-sequence and ionic effects on RNA hairpin folding. Here, we describe a general, three-dimensional (3D) conformation-based computational method for modeling salt concentration-dependent conformational distributions and the detailed 3D structures for a set of three RNA hairpins that contain a variable, 15-nucleotide loop sequence. For a given RNA sequence, the new, to our knowledge, method integrates a Vfold2D two-dimensional structure folding model with IsRNA coarse-grained molecular dynamics 3D folding simulations and Monte Carlo tightly bound ion estimations of ion-mediated electrostatic interactions. The model predicts free-energy landscapes for the different RNA hairpin-forming sequences with variable salt conditions. The theoretically predicted results agree with the experimental fluorescence measurements, validating the strategy. Furthermore, the theoretical model goes beyond the experimental results by enabling in-depth 3D structural analysis, revealing energetic mechanisms for the sequence- and salt-dependent folding stability. Although the computational framework presented here is developed for RNA hairpin systems, the general method may be applied to investigate other RNA systems, such as multiway junctions or pseudoknots in mixed metal ion solutions.
Collapse
Affiliation(s)
- Chenhan Zhao
- Department of Physics, Department of Biochemistry, and Institute for Data Science and Informatics, University of Missouri, Columbia, Missouri
| | - Dong Zhang
- Department of Physics, Department of Biochemistry, and Institute for Data Science and Informatics, University of Missouri, Columbia, Missouri
| | - Yangwei Jiang
- Department of Physics, Department of Biochemistry, and Institute for Data Science and Informatics, University of Missouri, Columbia, Missouri
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Institute for Data Science and Informatics, University of Missouri, Columbia, Missouri.
| |
Collapse
|
45
|
Osmer PS, Singh G, Boris-Lawrie K. A New Approach to 3D Modeling of Inhomogeneous Populations of Viral Regulatory RNA. Viruses 2020; 12:v12101108. [PMID: 33003639 PMCID: PMC7650772 DOI: 10.3390/v12101108] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Revised: 09/24/2020] [Accepted: 09/27/2020] [Indexed: 12/17/2022] Open
Abstract
Tertiary structure (3D) is the physical context of RNA regulatory activity. Retroviruses are RNA viruses that replicate through the proviral DNA intermediate transcribed by hosts. Proviral transcripts form inhomogeneous populations due to variable structural ensembles of overlapping regulatory RNA motifs in the 5′-untranslated region (UTR), which drive RNAs to be spliced or translated, and/or dimerized and packaged into virions. Genetic studies and structural techniques have provided fundamental input constraints to begin predicting HIV 3D conformations in silico. Using SimRNA and sets of experimentally-determined input constraints of HIVNL4-3 trans-activation responsive sequence (TAR) and pairings of unique-5′ (U5) with dimerization (DIS) or AUG motifs, we calculated a series of 3D models that differ in proximity of 5′-Cap and the junction of TAR and PolyA helices; configuration of primer binding site (PBS)-segment; and two host cofactors binding sites. Input constraints on U5-AUG pairings were most compatible with intramolecular folding of 5′-UTR motifs in energetic minima. Introducing theoretical constraints predicted metastable PolyA region drives orientation of 5′-Cap with TAR, U5 and PBS-segment helices. SimRNA and the workflow developed herein provides viable options to predict 3D conformations of inhomogeneous populations of large RNAs that have been intractable to conventional ensemble methods.
Collapse
Affiliation(s)
- Patrick S. Osmer
- Department of Astronomy, The Ohio State University, Columbus, OH 43210, USA;
| | - Gatikrushna Singh
- Department of Veterinary and Biomedical Sciences, University of Minnesota, Saint Paul, MN 55108, USA;
| | - Kathleen Boris-Lawrie
- Department of Veterinary and Biomedical Sciences, University of Minnesota, Saint Paul, MN 55108, USA;
- Correspondence: ; Tel.: +1-612-625-2100
| |
Collapse
|
46
|
Himmelstoß M, Erharter K, Renard E, Ennifar E, Kreutz C, Micura R. 2'- O-Trifluoromethylated RNA - a powerful modification for RNA chemistry and NMR spectroscopy. Chem Sci 2020; 11:11322-11330. [PMID: 34094374 PMCID: PMC8162808 DOI: 10.1039/d0sc04520a] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2020] [Accepted: 09/18/2020] [Indexed: 11/21/2022] Open
Abstract
New RNA modifications are needed to advance our toolbox for targeted manipulation of RNA. In particular, the development of high-performance reporter groups facilitating spectroscopic analysis of RNA structure and dynamics, and of RNA-ligand interactions has attracted considerable interest. To this end, fluorine labeling in conjunction with 19F-NMR spectroscopy has emerged as a powerful strategy. Appropriate probes for RNA previously focused on single fluorine atoms attached to the 5-position of pyrimidine nucleobases or at the ribose 2'-position. To increase NMR sensitivity, trifluoromethyl labeling approaches have been developed, with the ribose 2'-SCF3 modification being the most prominent one. A major drawback of the 2'-SCF3 group, however, is its strong impact on RNA base pairing stability. Interestingly, RNA containing the structurally related 2'-OCF3 modification has not yet been reported. Therefore, we set out to overcome the synthetic challenges toward 2'-OCF3 labeled RNA and to investigate the impact of this modification. We present the syntheses of 2'-OCF3 adenosine and cytidine phosphoramidites and their incorporation into oligoribonucleotides by solid-phase synthesis. Importantly, it turns out that the 2'-OCF3 group has only a slight destabilizing effect when located in double helical regions which is consistent with the preferential C3'-endo conformation of the 2'-OCF3 ribose as reflected in the 3 J (H1'-H2') coupling constants. Furthermore, we demonstrate the exceptionally high sensitivity of the new label in 19F-NMR analysis of RNA structure equilibria and of RNA-small molecule interactions. The study is complemented by a crystal structure at 0.9 Å resolution of a 27 nt hairpin RNA containing a single 2'-OCF3 group that well integrates into the minor groove. The new label carries high potential to outcompete currently applied fluorine labels for nucleic acid NMR spectroscopy because of its significantly advanced performance.
Collapse
Affiliation(s)
- Maximilian Himmelstoß
- University of Innsbruck, Institute of Organic Chemistry, Center for Molecular Biosciences (CMBI) Innrain 80-82 6020 Innsbruck Austria
| | - Kevin Erharter
- University of Innsbruck, Institute of Organic Chemistry, Center for Molecular Biosciences (CMBI) Innrain 80-82 6020 Innsbruck Austria
| | - Eva Renard
- Université de Strasbourg, Architecture et Réactivité de l'ARN-CNRS UPR 9002, Institut de Biologie Moléculaire et Cellulaire 67000 Strasbourg France
| | - Eric Ennifar
- Université de Strasbourg, Architecture et Réactivité de l'ARN-CNRS UPR 9002, Institut de Biologie Moléculaire et Cellulaire 67000 Strasbourg France
| | - Christoph Kreutz
- University of Innsbruck, Institute of Organic Chemistry, Center for Molecular Biosciences (CMBI) Innrain 80-82 6020 Innsbruck Austria
| | - Ronald Micura
- University of Innsbruck, Institute of Organic Chemistry, Center for Molecular Biosciences (CMBI) Innrain 80-82 6020 Innsbruck Austria
| |
Collapse
|
47
|
Xu X, Chen SJ. Topological constraints of RNA pseudoknotted and loop-kissing motifs: applications to three-dimensional structure prediction. Nucleic Acids Res 2020; 48:6503-6512. [PMID: 32491164 PMCID: PMC7337929 DOI: 10.1093/nar/gkaa463] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2020] [Accepted: 05/19/2020] [Indexed: 01/23/2023] Open
Abstract
An RNA global fold can be described at the level of helix orientations and relatively flexible loop conformations that connect the helices. The linkage between the helices plays an essential role in determining the structural topology, which restricts RNA local and global folds, especially for RNA tertiary structures involving cross-linked base pairs. We quantitatively analyze the topological constraints on RNA 3D conformational space, in particular, on the distribution of helix orientations, for pseudoknots and loop-loop kissing structures. The result shows that a viable conformational space is predominantly determined by the motif type, helix size, and loop size, indicating a strong topological coupling between helices and loops in RNA tertiary motifs. Moreover, the analysis indicates that (cross-linked) tertiary contacts can cause much stronger topological constraints on RNA global fold than non-cross-linked base pairs. Furthermore, based on the topological constraints encoded in the 2D structure and the 3D templates, we develop a 3D structure prediction approach. This approach can be further combined with structure probing methods to expand the capability of computational prediction for large RNA folds.
Collapse
Affiliation(s)
- Xiaojun Xu
- Institute of Bioinformatics and Medical Engineering, Jiangsu University of Technology, Changzhou, Jiangsu 213001, China
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, MO 65211, USA
| |
Collapse
|
48
|
Mishra B, Balaji A, Beesetti H, Swaminathan S, Aduri R. The RNA secondary structural variation in the cyclization elements of the dengue genome and the possible implications in pathogenicity. Virusdisease 2020; 31:299-307. [PMID: 32904896 PMCID: PMC7458965 DOI: 10.1007/s13337-020-00615-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2020] [Accepted: 07/20/2020] [Indexed: 10/23/2022] Open
Abstract
Dengue virus (DENV), the causative agent of dengue fever and severe dengue, exists as four antigenically different serotypes. These serotypes are further classified into genotypes and have varying degrees of pathogenicity. The 5' and 3' ends of the genomic RNA play a critical role in the viral life cycle. A global scale study of the RNA structural variation among the sero- and genotypes was carried out to correlate RNA structure with pathogenicity. We found that the GC rich stem and rigid loop structure of the 5' end of the genomic RNA of DENV 2 differs significantly from the others. The observed variation in base composition and base pairing may confer structural and functional advantage in highly virulent strains. This variation in the structure may influence the ease of cyclization and recruitment of viral RNA polymerase, NS5 RdRp, thereby affecting the pathogenicity of these strains.
Collapse
Affiliation(s)
- Bibhudutta Mishra
- Department of Biological Sciences, Birla Institute of Technology and Science, Pilani, K K Birla Goa Campus, Zuarinagar, South Goa, Goa 403 726 India
| | - Advait Balaji
- Department of Biological Sciences, Birla Institute of Technology and Science, Pilani, K K Birla Goa Campus, Zuarinagar, South Goa, Goa 403 726 India
| | - Hemalatha Beesetti
- Department of Biological Sciences, Birla Institute of Technology and Science, Pilani, Hyderabad Campus, Jawahar Nagar, Shameerpet Mandal, Hyderabad, Telangana 500 078 India
- Present Address: Molecular Medicine Division, Recombinant Gene Products Group, International Centre for Genetic Engineering and Biotechnology, New Delhi, 110067 India
| | - Sathyamangalam Swaminathan
- Department of Biological Sciences, Birla Institute of Technology and Science, Pilani, Hyderabad Campus, Jawahar Nagar, Shameerpet Mandal, Hyderabad, Telangana 500 078 India
- Present Address: Molecular Medicine Division, Recombinant Gene Products Group, International Centre for Genetic Engineering and Biotechnology, New Delhi, 110067 India
| | - Raviprasad Aduri
- Department of Biological Sciences, Birla Institute of Technology and Science, Pilani, K K Birla Goa Campus, Zuarinagar, South Goa, Goa 403 726 India
| |
Collapse
|
49
|
Fard EM, Moradi S, Salekdeh NN, Bakhshi B, Ghaffari MR, Zeinalabedini M, Salekdeh GH. Plant isomiRs: origins, biogenesis, and biological functions. Genomics 2020; 112:3382-3395. [DOI: 10.1016/j.ygeno.2020.06.019] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2020] [Revised: 05/22/2020] [Accepted: 06/10/2020] [Indexed: 12/14/2022]
|
50
|
Prosdocimi F, Zamudio GS, Palacios-Pérez M, Torres de Farias S, V. José M. The Ancient History of Peptidyl Transferase Center Formation as Told by Conservation and Information Analyses. Life (Basel) 2020; 10:life10080134. [PMID: 32764248 PMCID: PMC7459865 DOI: 10.3390/life10080134] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Revised: 07/24/2020] [Accepted: 07/31/2020] [Indexed: 12/19/2022] Open
Abstract
The peptidyl transferase center (PTC) is the catalytic center of the ribosome and forms part of the 23S ribosomal RNA. The PTC has been recognized as the earliest ribosomal part and its origins embodied the First Universal Common Ancestor (FUCA). The PTC is frequently assumed to be highly conserved along all living beings. In this work, we posed the following questions: (i) How many 100% conserved bases can be found in the PTC? (ii) Is it possible to identify clusters of informationally linked nucleotides along its sequence? (iii) Can we propose how the PTC was formed? (iv) How does sequence conservation reflect on the secondary and tertiary structures of the PTC? Aiming to answer these questions, all available complete sequences of 23S ribosomal RNA from Bacteria and Archaea deposited on GenBank database were downloaded. Using a sequence bait of 179 bp from the PTC of Thermus termophilus, we performed an optimum pairwise alignment to retrieve the PTC region from 1424 filtered 23S rRNA sequences. These PTC sequences were multiply aligned, and the conserved regions were assigned and observed along the primary, secondary, and tertiary structures. The PTC structure was observed to be more highly conserved close to the adenine located at the catalytical site. Clusters of interrelated, co-evolving nucleotides reinforce previous assumptions that the PTC was formed by the concatenation of proto-tRNAs and important residues responsible for its assembly were identified. The observed sequence variation does not seem to significantly affect the 3D structure of the PTC ribozyme.
Collapse
Affiliation(s)
- Francisco Prosdocimi
- Laboratório de Biologia Teórica e de Sistemas, Instituto de Bioquímica Médica Leopoldo de Meis, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21.941-902, Brazil
- Theoretical Biology Group, Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México, Ciudad Universitaria, CDMX 04510, Mexico; (G.S.Z.); (M.P.-P.)
- Correspondence: (F.P.); (M.V.J.)
| | - Gabriel S. Zamudio
- Theoretical Biology Group, Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México, Ciudad Universitaria, CDMX 04510, Mexico; (G.S.Z.); (M.P.-P.)
| | - Miryam Palacios-Pérez
- Theoretical Biology Group, Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México, Ciudad Universitaria, CDMX 04510, Mexico; (G.S.Z.); (M.P.-P.)
| | - Sávio Torres de Farias
- Laboratório de Genética Evolutiva Paulo Leminsk, Departamento de Biologia Molecular, Universidade Federal da Paraíba, João Pessoa, Paraíba 58051-900, Brazil;
| | - Marco V. José
- Theoretical Biology Group, Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México, Ciudad Universitaria, CDMX 04510, Mexico; (G.S.Z.); (M.P.-P.)
- Correspondence: (F.P.); (M.V.J.)
| |
Collapse
|