1
|
Deng P, Tan SQ, Yang QY, Fu L, Wu Y, Zhu HZ, Sun L, Bao Z, Lin Y, Zhang QC, Wang H, Wang J, Liu JJG. Structural RNA components supervise the sequential DNA cleavage in R2 retrotransposon. Cell 2023; 186:2865-2879.e20. [PMID: 37301196 DOI: 10.1016/j.cell.2023.05.032] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2023] [Revised: 03/14/2023] [Accepted: 05/20/2023] [Indexed: 06/12/2023]
Abstract
Retroelements are the widespread jumping elements considered as major drivers for genome evolution, which can also be repurposed as gene-editing tools. Here, we determine the cryo-EM structures of eukaryotic R2 retrotransposon with ribosomal DNA target and regulatory RNAs. Combined with biochemical and sequencing analysis, we reveal two essential DNA regions, Drr and Dcr, required for recognition and cleavage. The association of 3' regulatory RNA with R2 protein accelerates the first-strand cleavage, blocks the second-strand cleavage, and initiates the reverse transcription starting from the 3'-tail. Removing 3' regulatory RNA by reverse transcription allows the association of 5' regulatory RNA and initiates the second-strand cleavage. Taken together, our work explains the DNA recognition and RNA supervised sequential retrotransposition mechanisms by R2 machinery, providing insights into the retrotransposon and application reprogramming.
Collapse
Affiliation(s)
- Pujuan Deng
- Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, School of Life Sciences, Tsinghua University, Beijing 100084, China
| | - Shun-Qing Tan
- Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, School of Life Sciences, Tsinghua University, Beijing 100084, China
| | - Qi-Yu Yang
- Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, School of Life Sciences, Tsinghua University, Beijing 100084, China
| | - Liangzheng Fu
- State Key Laboratory of Stem Cell and Reproductive Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yachao Wu
- State Key Laboratory of Stem Cell and Reproductive Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Han-Zhou Zhu
- Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, School of Life Sciences, Tsinghua University, Beijing 100084, China
| | - Lei Sun
- Shandong Provincial Key Laboratory of Animal Cell and Developmental Biology, School of Life Sciences, Shandong University, Qingdao 266237, China
| | - Zhangbin Bao
- Tsinghua-Peking Center for Life Sciences, School of Life Sciences, Tsinghua University, Beijing 100084, China; IDG/McGovern Institute for Brain Research, Tsinghua University, Beijing 100084, China
| | - Yi Lin
- Tsinghua-Peking Center for Life Sciences, School of Life Sciences, Tsinghua University, Beijing 100084, China; IDG/McGovern Institute for Brain Research, Tsinghua University, Beijing 100084, China
| | - Qiangfeng Cliff Zhang
- Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, School of Life Sciences, Tsinghua University, Beijing 100084, China
| | - Haoyi Wang
- State Key Laboratory of Stem Cell and Reproductive Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Jia Wang
- Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, School of Life Sciences, Tsinghua University, Beijing 100084, China.
| | - Jun-Jie Gogo Liu
- Beijing Advanced Innovation Center for Structural Biology & Frontier Research Center for Biological Structure, School of Life Sciences, Tsinghua University, Beijing 100084, China; Tsinghua-Peking Center for Life Sciences, School of Life Sciences, Tsinghua University, Beijing 100084, China.
| |
Collapse
|
2
|
Moudgal N, Arhin G, Frank AT. Using Unassigned NMR Chemical Shifts to Model RNA Secondary Structure. J Phys Chem A 2022; 126:2739-2745. [PMID: 35470661 DOI: 10.1021/acs.jpca.2c00456] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
NMR-derived chemical shifts are sensitive probes of RNA structure. However, the need to assign NMR spectra hampers their utility as a direct source of structural information. In this report, we describe a simple method that uses unassigned 2D NMR spectra to model the secondary structure of RNAs. As in the case of assigned chemical shifts, we could use unassigned chemical shift data to reweight conformational libraries such that the highest weighted structure closely resembles their reference NMR structure. Furthermore, the application of our approach to the 3'- and 5'-UTR of the SARS-CoV-2 genome yields structures that are, for the most part, consistent with the secondary structure models derived from chemical probing data. Therefore, we expect the framework we describe here will be useful as a general strategy for rapidly generating preliminary structural RNA models directly from unassigned 2D NMR spectra. As we demonstrated for the 337-nt and 472-nt UTRs of SARS-CoV-2, our approach could be especially valuable for modeling the secondary structures of large RNA.
Collapse
Affiliation(s)
- Neel Moudgal
- Saline High School, 1300 Campus Pkwy, Saline, Michigan 48176, United States
| | - Grace Arhin
- Biophysics Program, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109, United States
| | - Aaron T Frank
- Biophysics Program, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109, United States.,Chemistry Department, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
3
|
Mandal T, Dhara K, Parui N, Dash J. Domino Relay Olefin Metathesis of Triallyl Oxindole and Indole Precursors to Access Cyclic Indoxyls and Carbazoles. ChemCatChem 2020. [DOI: 10.1002/cctc.202000813] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Tirtha Mandal
- School of Chemical Sciences Indian Association for the Cultivation of Science Jadavpur, Kolkata 700032 India
| | - Kalyan Dhara
- School of Chemical Sciences Indian Association for the Cultivation of Science Jadavpur, Kolkata 700032 India
| | - Nabin Parui
- School of Chemical Sciences Indian Association for the Cultivation of Science Jadavpur, Kolkata 700032 India
| | - Jyotirmayee Dash
- School of Chemical Sciences Indian Association for the Cultivation of Science Jadavpur, Kolkata 700032 India
| |
Collapse
|
4
|
Antunes D, Jorge NAN, Caffarena ER, Passetti F. Using RNA Sequence and Structure for the Prediction of Riboswitch Aptamer: A Comprehensive Review of Available Software and Tools. Front Genet 2018; 8:231. [PMID: 29403526 PMCID: PMC5780412 DOI: 10.3389/fgene.2017.00231] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2017] [Accepted: 12/21/2017] [Indexed: 12/14/2022] Open
Abstract
RNA molecules are essential players in many fundamental biological processes. Prokaryotes and eukaryotes have distinct RNA classes with specific structural features and functional roles. Computational prediction of protein structures is a research field in which high confidence three-dimensional protein models can be proposed based on the sequence alignment between target and templates. However, to date, only a few approaches have been developed for the computational prediction of RNA structures. Similar to proteins, RNA structures may be altered due to the interaction with various ligands, including proteins, other RNAs, and metabolites. A riboswitch is a molecular mechanism, found in the three kingdoms of life, in which the RNA structure is modified by the binding of a metabolite. It can regulate multiple gene expression mechanisms, such as transcription, translation initiation, and mRNA splicing and processing. Due to their nature, these entities also act on the regulation of gene expression and detection of small metabolites and have the potential to helping in the discovery of new classes of antimicrobial agents. In this review, we describe software and web servers currently available for riboswitch aptamer identification and secondary and tertiary structure prediction, including applications.
Collapse
Affiliation(s)
- Deborah Antunes
- Scientific Computing Program (PROCC), Computational Biophysics and Molecular Modeling Group, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
| | - Natasha A N Jorge
- Laboratory of Functional Genomics and Bioinformatics, Oswaldo Cruz Institute, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil.,Laboratory of Gene Expression Regulation, Carlos Chagas Institute, Fundação Oswaldo Cruz, Curitiba, Brazil
| | - Ernesto R Caffarena
- Scientific Computing Program (PROCC), Computational Biophysics and Molecular Modeling Group, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
| | - Fabio Passetti
- Laboratory of Functional Genomics and Bioinformatics, Oswaldo Cruz Institute, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil.,Laboratory of Gene Expression Regulation, Carlos Chagas Institute, Fundação Oswaldo Cruz, Curitiba, Brazil
| |
Collapse
|
5
|
Abstract
Dynamic programming methods for predicting RNA secondary structure often use thermodynamics and experimental restraints and/or constraints to limit folding space. Chemical mapping results typically restrain certain nucleotides not to be in AU or GC pairs. Two-dimensional nuclear magnetic resonance (NMR) spectra can reveal the order of AU, GC, and GU pairs in double helixes. This chapter describes a program, NMR-assisted prediction of secondary structure and chemical shifts (NAPSS-CS), that constrains possible secondary structures on the basis of the NMR determined order and 5'-3' direction of AU, GC, and GU pairs in helixes. NAPSS-CS minimally requires input of the order of base pairs as determined from nuclear Overhauser effect spectroscopy (NOESY) of imino protons. The program deduces the 5'-3' direction of the base pairs if certain chemical shifts are also input. Secondary structures predicted by the program provide assignments of input chemical shifts to particular nucleotides in the sequence, thus facilitating an important step for determination of the three dimensional structure by NMR. The method is particularly useful for revealing pseudoknots and an example is provided. The method may also allow determination of secondary structures when a sequence folds into two structures that exchange slowly.
Collapse
|
6
|
Kauffmann AD, Kennedy SD, Zhao J, Turner DH. Nuclear Magnetic Resonance Structure of an 8 × 8 Nucleotide RNA Internal Loop Flanked on Each Side by Three Watson-Crick Pairs and Comparison to Three-Dimensional Predictions. Biochemistry 2017; 56:3733-3744. [PMID: 28700212 DOI: 10.1021/acs.biochem.7b00201] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The prediction of RNA three-dimensional structure from sequence alone has been a long-standing goal. High-resolution, experimentally determined structures of simple noncanonical pairings and motifs are critical to the development of prediction programs. Here, we present the nuclear magnetic resonance structure of the (5'CCAGAAACGGAUGGA)2 duplex, which contains an 8 × 8 nucleotide internal loop flanked by three Watson-Crick pairs on each side. The loop is comprised of a central 5'AC/3'CA nearest neighbor flanked by two 3RRs motifs, a known stable motif consisting of three consecutive sheared GA pairs. Hydrogen bonding patterns between base pairs in the loop, the all-atom root-mean-square deviation for the loop, and the deformation index were used to compare the structure to automated predictions by MC-sym, RNA FARFAR, and RNAComposer.
Collapse
Affiliation(s)
- Andrew D Kauffmann
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States.,Center for RNA Biology, University of Rochester , Rochester, New York 14627, United States
| | - Scott D Kennedy
- Department of Biochemistry and Biophysics, School of Medicine & Dentistry, University of Rochester , Rochester, New York 14642, United States.,Center for RNA Biology, University of Rochester , Rochester, New York 14627, United States
| | - Jianbo Zhao
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States.,Center for RNA Biology, University of Rochester , Rochester, New York 14627, United States
| | - Douglas H Turner
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States.,Center for RNA Biology, University of Rochester , Rochester, New York 14627, United States
| |
Collapse
|
7
|
Abstract
Experimental probing data can be used to improve the accuracy of RNA secondary structure prediction. The software package RNAstructure can take advantage of enzymatic cleavage data, FMN cleavage data, traditional chemical modification reactivity data, and SHAPE reactivity data for secondary structure modeling. This chapter provides protocols for using experimental probing data with RNAstructure to restrain or constrain RNA secondary structure prediction.
Collapse
Affiliation(s)
- Zhenjiang Zech Xu
- Department of Biochemistry & Biophysics, University of Rochester Medical Center, 601 Elmwood Avenue, Box 712, Rochester, NY, 14642, USA
- Center for RNA Biology, University of Rochester Medical Center, 601 Elmwood Avenue, Box 712, Rochester, NY, 14642, USA
| | - David H Mathews
- Department of Biochemistry & Biophysics, University of Rochester Medical Center, 601 Elmwood Avenue, Box 712, Rochester, NY, 14642, USA.
- Center for RNA Biology, University of Rochester Medical Center, 601 Elmwood Avenue, Box 712, Rochester, NY, 14642, USA.
- Department of Biostatistics & Computational Biology, University of Rochester Medical Center, 601 Elmwood Avenue, Box 712, Rochester, NY, 14642, USA.
| |
Collapse
|
8
|
Abstract
Knowledge of RNA secondary structure is often sufficient to identify relationships between the structure of RNA and processing pathways, and the design of therapeutics. Nuclear magnetic resonance (NMR) can identify types of nucleotide base pairs and the sequence, thus limiting possible secondary structures. Because NMR experiments, like chemical mapping, are performed in solution, not in single crystals, experiments can be initiated as soon as the biomolecule is expressed and purified. This chapter summarizes NMR methods that permit rapid identification of RNA secondary structure, information that can be used as supplements to chemical mapping, and/or as preliminary steps required for 3D structure determination. The primary aim is to provide guidelines to enable a researcher with minimal knowledge of NMR to quickly extract secondary structure information from basic datasets. Instrumental and sample considerations that can maximize data quality are discussed along with some details for optimal data acquisition and processing parameters. Approaches for identifying base pair types in both unlabeled and isotopically labeled RNA are covered. Common problems, such as missing signals and overlaps, and approaches to address them are considered. Programs under development for merging NMR data with structure prediction algorithms are briefly discussed.
Collapse
Affiliation(s)
- Scott D Kennedy
- Department of Biochemistry and Biophysics, School of Medicine and Dentistry, University of Rochester, Rochester, NY, 14642, USA.
| |
Collapse
|
9
|
Chen JL, Bellaousov S, Tubbs JD, Kennedy SD, Lopez MJ, Mathews DH, Turner DH. Nuclear Magnetic Resonance-Assisted Prediction of Secondary Structure for RNA: Incorporation of Direction-Dependent Chemical Shift Constraints. Biochemistry 2015; 54:6769-82. [PMID: 26451676 PMCID: PMC4666457 DOI: 10.1021/acs.biochem.5b00833] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Knowledge
of RNA
structure is necessary to determine structure–function relationships
and to facilitate design of potential therapeutics.
RNA secondary structure prediction can be improved by applying constraints
from nuclear magnetic resonance (NMR) experiments to a dynamic programming
algorithm. Imino proton walks from NOESY spectra reveal double-stranded
regions. Chemical shifts of protons in GH1, UH3, and UH5 of GU pairs,
UH3, UH5, and AH2 of AU pairs, and GH1 of GC pairs were analyzed to
identify constraints for the 5′ to 3′ directionality
of base pairs in helices. The 5′ to 3′ directionality
constraints were incorporated into an NMR-assisted prediction of secondary
structure (NAPSS-CS) program. When it was tested on 18 structures,
including nine pseudoknots, the sensitivity and positive predictive
value were improved relative to those of three unrestrained programs.
The prediction accuracy for the pseudoknots improved the most. The
program also facilitates assignment of chemical shifts to individual
nucleotides, a necessary step for determining three-dimensional structure.
Collapse
Affiliation(s)
- Jonathan L Chen
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States
| | - Stanislav Bellaousov
- Department of Biochemistry and Biophysics, University of Rochester School of Medicine and Dentistry , Rochester, New York 14642, United States
| | - Jason D Tubbs
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States
| | - Scott D Kennedy
- Department of Biochemistry and Biophysics, University of Rochester School of Medicine and Dentistry , Rochester, New York 14642, United States
| | - Michael J Lopez
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States
| | - David H Mathews
- Department of Biochemistry and Biophysics, University of Rochester School of Medicine and Dentistry , Rochester, New York 14642, United States.,Center for RNA Biology, University of Rochester , Rochester, New York 14642, United States
| | - Douglas H Turner
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States.,Center for RNA Biology, University of Rochester , Rochester, New York 14642, United States
| |
Collapse
|
10
|
Kierzek R, Turner DH, Kierzek E. Microarrays for identifying binding sites and probing structure of RNAs. Nucleic Acids Res 2014; 43:1-12. [PMID: 25505162 PMCID: PMC4288193 DOI: 10.1093/nar/gku1303] [Citation(s) in RCA: 75] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Oligonucleotide microarrays are widely used in various biological studies. In this review, application of oligonucleotide microarrays for identifying binding sites and probing structure of RNAs is described. Deep sequencing allows fast determination of DNA and RNA sequence. High-throughput methods for determination of secondary structures of RNAs have also been developed. Those methods, however, do not reveal binding sites for oligonucleotides. In contrast, microarrays directly determine binding sites while also providing structural insights. Microarray mapping can be used over a wide range of experimental conditions, including temperature, pH, various cations at different concentrations and the presence of other molecules. Moreover, it is possible to make universal microarrays suitable for investigations of many different RNAs, and readout of results is rapid. Thus, microarrays are used to provide insight into oligonucleotide sequences potentially able to interfere with biological function. Better understanding of structure-function relationships of RNA can be facilitated by using microarrays to find RNA regions capable to bind oligonucleotides. That information is extremely important to design optimal sequences for antisense oligonucleotides and siRNA because both bind to single-stranded regions of target RNAs.
Collapse
Affiliation(s)
- Ryszard Kierzek
- Institute of Bioorganic Chemistry Polish Academy of Sciences, 61-704 Poznan, Noskowskiego 12/14, Poland
| | - Douglas H Turner
- Department of Chemistry, University of Rochester, Rochester, NY 14627, USA
| | - Elzbieta Kierzek
- Institute of Bioorganic Chemistry Polish Academy of Sciences, 61-704 Poznan, Noskowskiego 12/14, Poland
| |
Collapse
|
11
|
Strulson CA, Boyer JA, Whitman EE, Bevilacqua PC. Molecular crowders and cosolutes promote folding cooperativity of RNA under physiological ionic conditions. RNA (NEW YORK, N.Y.) 2014; 20:331-47. [PMID: 24442612 PMCID: PMC3923128 DOI: 10.1261/rna.042747.113] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/28/2013] [Accepted: 11/22/2013] [Indexed: 05/21/2023]
Abstract
Folding mechanisms of functional RNAs under idealized in vitro conditions of dilute solution and high ionic strength have been well studied. Comparatively little is known, however, about mechanisms for folding of RNA in vivo where Mg(2+) ion concentrations are low, K(+) concentrations are modest, and concentrations of macromolecular crowders and low-molecular-weight cosolutes are high. Herein, we apply a combination of biophysical and structure mapping techniques to tRNA to elucidate thermodynamic and functional principles that govern RNA folding under in vivo-like conditions. We show by thermal denaturation and SHAPE studies that tRNA folding cooperativity increases in physiologically low concentrations of Mg(2+) (0.5-2 mM) and K(+) (140 mM) if the solution is supplemented with physiological amounts (∼ 20%) of a water-soluble neutral macromolecular crowding agent such as PEG or dextran. Low-molecular-weight cosolutes show varying effects on tRNA folding cooperativity, increasing or decreasing it based on the identity of the cosolute. For those additives that increase folding cooperativity, the gain is manifested in sharpened two-state-like folding transitions for full-length tRNA over its secondary structural elements. Temperature-dependent SHAPE experiments in the absence and presence of crowders and cosolutes reveal extent of cooperative folding of tRNA on a nucleotide basis and are consistent with the melting studies. Mechanistically, crowding agents appear to promote cooperativity by stabilizing tertiary structure, while those low molecular cosolutes that promote cooperativity stabilize tertiary structure and/or destabilize secondary structure. Cooperative folding of functional RNA under physiological-like conditions parallels the behavior of many proteins and has implications for cellular RNA folding kinetics and evolution.
Collapse
Affiliation(s)
- Christopher A. Strulson
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Center for RNA Molecular Biology, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - Joshua A. Boyer
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Center for RNA Molecular Biology, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - Elisabeth E. Whitman
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Center for RNA Molecular Biology, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - Philip C. Bevilacqua
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Center for RNA Molecular Biology, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Corresponding authorE-mail
| |
Collapse
|
12
|
Tubbs J, Condon DE, Kennedy S, Hauser M, Bevilacqua P, Turner DH. The nuclear magnetic resonance of CCCC RNA reveals a right-handed helix, and revised parameters for AMBER force field torsions improve structural predictions from molecular dynamics. Biochemistry 2013; 52:996-1010. [PMID: 23286901 PMCID: PMC3571207 DOI: 10.1021/bi3010347] [Citation(s) in RCA: 64] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2012] [Revised: 12/27/2012] [Indexed: 12/29/2022]
Abstract
The sequence dependence of RNA energetics is important for predicting RNA structure. Hairpins with C(n) loops are consistently less stable than hairpins with other loops, which suggests the structure of C(n) regions could be unusual in the "unfolded" state. For example, previous nuclear magnetic resonance (NMR) evidence suggested that polycytidylic acid forms a left-handed helix. In this study, UV melting experiments show that the hairpin formed by r(5'GGACCCCCGUCC) is less stable than r(5'GGACUUUUGUCC). NMR spectra for single-stranded C(4) oligonucleotide, mimicking the unfolded hairpin loop, are consistent with a right-handed A-form-like helix. Comparisons between NMR spectra and molecular dynamics (MD) simulations suggest that recent reparametrizations, parm99χ_YIL and parm99TOR, of the AMBER parm99 force field improve the agreement between structural features for C(4) determined by NMR and predicted by MD. Evidently, the force field revisions to parm99 improve the modeling of RNA energetics and therefore structure.
Collapse
Affiliation(s)
- Jason
D. Tubbs
- Department
of Chemistry, University of Rochester,
Rochester, New York 14627,
United States
- Center for
RNA Biology, University of Rochester, Rochester,
New York 14642,
United States
| | - David E. Condon
- Department
of Chemistry, University of Rochester,
Rochester, New York 14627,
United States
- Center for
RNA Biology, University of Rochester, Rochester,
New York 14642,
United States
| | - Scott
D. Kennedy
- Center for
RNA Biology, University of Rochester, Rochester,
New York 14642,
United States
- Department
of Biochemistry and
Biophysics, University of Rochester, Rochester,
New York 14642, United States
| | - Melanie Hauser
- Department of Chemistry, Buena Vista University, Storm Lake, Iowa 50588, United
States
- Department of Chemistry
and
Center for RNA Molecular Biology, The Pennsylvania State
University, University Park, Pennsylvania 16802, United
States
| | - Philip
C. Bevilacqua
- Department of Chemistry
and
Center for RNA Molecular Biology, The Pennsylvania State
University, University Park, Pennsylvania 16802, United
States
| | - Douglas H. Turner
- Department
of Chemistry, University of Rochester,
Rochester, New York 14627,
United States
- Center for
RNA Biology, University of Rochester, Rochester,
New York 14642,
United States
| |
Collapse
|
13
|
Abstract
The purpose of this section is to detail methods for the computational prediction of RNA secondary structure. This protocol is intended to provide an easy entry into the field of RNA structure prediction for those wishing to utilize it in their research and to suggest 'best practices' for going from sequence to secondary structure depending on the available data.
Collapse
Affiliation(s)
- Walter N Moss
- Department of Chemistry, University of Rochester, Rochester, NY, USA.
| |
Collapse
|
14
|
Chen J, Dishler AL, Kennedy SD, Yildirim I, Liu B, Turner DH, Serra MJ. Testing the nearest neighbor model for canonical RNA base pairs: revision of GU parameters. Biochemistry 2012; 51:3508-22. [PMID: 22490167 PMCID: PMC3335265 DOI: 10.1021/bi3002709] [Citation(s) in RCA: 69] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2012] [Indexed: 11/30/2022]
Abstract
Thermodynamic parameters for GU pairs are important for predicting the secondary structures of RNA and for finding genomic sequences that code for structured RNA. Optical melting curves were measured for 29 RNA duplexes with GU pairs to improve nearest neighbor parameters for predicting stabilities of helixes. The updated model eliminates a prior penalty assumed for terminal GU pairs. Six additional duplexes with the 5'GG/3'UU motif were added to the single representation in the previous database. This revises the ΔG°(37) for the 5'GG/3'UU motif from an unfavorable 0.5 kcal/mol to a favorable -0.2 kcal/mol. Similarly, the ΔG°(37) for the 5'UG/3'GU motif changes from 0.3 to -0.6 kcal/mol. The correlation coefficients between predicted and experimental ΔG°(37), ΔH°, and ΔS° for the expanded database are 0.95, 0.89, and 0.87, respectively. The results should improve predictions of RNA secondary structure.
Collapse
Affiliation(s)
- Jonathan
L. Chen
- Department
of Chemistry, University of Rochester,
Rochester, New York 14627, United States
| | - Abigael L. Dishler
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, United States
| | - Scott D. Kennedy
- Department of Biochemistry and Biophysics, University of Rochester School of Medicine and Dentistry, Rochester, New York 14642, United States
| | - Ilyas Yildirim
- Department
of Chemistry, University of Rochester,
Rochester, New York 14627, United States
| | - Biao Liu
- Department
of Chemistry, University of Rochester,
Rochester, New York 14627, United States
| | - Douglas H. Turner
- Department
of Chemistry, University of Rochester,
Rochester, New York 14627, United States
- Center for RNA Biology, University of Rochester, Rochester, New York 14627, United States
| | - Martin J. Serra
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, United States
| |
Collapse
|
15
|
Bahrami A, Clos LJ, Markley JL, Butcher SE, Eghbalnia HR. RNA-PAIRS: RNA probabilistic assignment of imino resonance shifts. JOURNAL OF BIOMOLECULAR NMR 2012; 52:289-302. [PMID: 22359049 PMCID: PMC3480180 DOI: 10.1007/s10858-012-9603-z] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/06/2011] [Accepted: 01/08/2012] [Indexed: 05/13/2023]
Abstract
The significant biological role of RNA has further highlighted the need for improving the accuracy, efficiency and the reach of methods for investigating RNA structure and function. Nuclear magnetic resonance (NMR) spectroscopy is vital to furthering the goals of RNA structural biology because of its distinctive capabilities. However, the dispersion pattern in the NMR spectra of RNA makes automated resonance assignment, a key step in NMR investigation of biomolecules, remarkably challenging. Herein we present RNA Probabilistic Assignment of Imino Resonance Shifts (RNA-PAIRS), a method for the automated assignment of RNA imino resonances with synchronized verification and correction of predicted secondary structure. RNA-PAIRS represents an advance in modeling the assignment paradigm because it seeds the probabilistic network for assignment with experimental NMR data, and predicted RNA secondary structure, simultaneously and from the start. Subsequently, RNA-PAIRS sets in motion a dynamic network that reverberates between predictions and experimental evidence in order to reconcile and rectify resonance assignments and secondary structure information. The procedure is halted when assignments and base-parings are deemed to be most consistent with observed crosspeaks. The current implementation of RNA-PAIRS uses an initial peak list derived from proton-nitrogen heteronuclear multiple quantum correlation ((1)H-(15)N 2D HMQC) and proton-proton nuclear Overhauser enhancement spectroscopy ((1)H-(1)H 2D NOESY) experiments. We have evaluated the performance of RNA-PAIRS by using it to analyze NMR datasets from 26 previously studied RNAs, including a 111-nucleotide complex. For moderately sized RNA molecules, and over a range of comparatively complex structural motifs, the average assignment accuracy exceeds 90%, while the average base pair prediction accuracy exceeded 93%. RNA-PAIRS yielded accurate assignments and base pairings consistent with imino resonances for a majority of the NMR resonances, even when the initial predictions are only modestly accurate. RNA-PAIRS is available as a public web-server at http://pine.nmrfam.wisc.edu/RNA/.
Collapse
Affiliation(s)
- Arash Bahrami
- National Magnetic Resonance Facility at Madison, Madison, WI, USA
| | - Lawrence J. Clos
- National Magnetic Resonance Facility at Madison, Madison, WI, USA
| | - John L. Markley
- National Magnetic Resonance Facility at Madison, Madison, WI, USA. Biochemistry Department, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Samuel E. Butcher
- National Magnetic Resonance Facility at Madison, Madison, WI, USA. Biochemistry Department, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Hamid R. Eghbalnia
- Department of Molecular and Cellular Physiology, University of Cincinnati, P.O. Box 670576, Cincinnati, OH 45267-0576, USA
| |
Collapse
|
16
|
Puglisi EV, Puglisi JD. Secondary structure of the HIV reverse transcription initiation complex by NMR. J Mol Biol 2011; 410:863-74. [PMID: 21763492 DOI: 10.1016/j.jmb.2011.04.024] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2011] [Revised: 04/06/2011] [Accepted: 04/11/2011] [Indexed: 01/11/2023]
Abstract
Initiation of reverse transcription of genomic RNA is a key early step in replication of the human immunodeficiency virus (HIV) upon infection of a host cell. Viral reverse transcriptase initiates from a specific RNA-RNA complex formed between a host transfer RNA (tRNA(Lys)(3)) and a region at the 5' end of genomic RNA; the 3' end of the tRNA acts as a primer for reverse transcription of genomic RNA. We report here the secondary structure of the HIV genomic RNA-human tRNA(Lys)(3) initiation complex using heteronuclear nuclear magnetic resonance methods. We show that both RNAs undergo large-scale conformational changes upon complex formation. Formation of the 18-bp primer helix with the 3' end of tRNA(Lys)(3) drives large conformational rearrangements of the tRNA at the 5' end while maintaining the anticodon loop for potential loop-loop interactions. HIV RNA forms an intramolecular helix adjacent to the intermolecular primer helix. This helix, which must be broken by reverse transcription, likely acts as a kinetic block to reverse transcription.
Collapse
Affiliation(s)
- Elisabetta Viani Puglisi
- Department of Structural Biology, Stanford University School of Medicine, Stanford, CA 94305-5126, USA.
| | | |
Collapse
|
17
|
Lerman YV, Kennedy SD, Shankar N, Parisien M, Major F, Turner DH. NMR structure of a 4 x 4 nucleotide RNA internal loop from an R2 retrotransposon: identification of a three purine-purine sheared pair motif and comparison to MC-SYM predictions. RNA (NEW YORK, N.Y.) 2011; 17:1664-77. [PMID: 21778280 PMCID: PMC3162332 DOI: 10.1261/rna.2641911] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/24/2011] [Accepted: 05/08/2011] [Indexed: 05/31/2023]
Abstract
The NMR solution structure is reported of a duplex, 5'GUGAAGCCCGU/3'UCACAGGAGGC, containing a 4 × 4 nucleotide internal loop from an R2 retrotransposon RNA. The loop contains three sheared purine-purine pairs and reveals a structural element found in other RNAs, which we refer to as the 3RRs motif. Optical melting measurements of the thermodynamics of the duplex indicate that the internal loop is 1.6 kcal/mol more stable at 37°C than predicted. The results identify the 3RRs motif as a common structural element that can facilitate prediction of 3D structure. Known examples include internal loops having the pairings: 5'GAA/3'AGG, 5'GAG/3'AGG, 5'GAA/3'AAG, and 5'AAG/3'AGG. The structural information is compared with predictions made with the MC-Sym program.
Collapse
Affiliation(s)
- Yelena V. Lerman
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | - Scott D. Kennedy
- Department of Biochemistry and Biophysics, University of Rochester School of Medicine and Dentistry, Rochester, New York 14642, USA
| | - Neelaabh Shankar
- Department of Biochemistry and Biophysics, University of Rochester School of Medicine and Dentistry, Rochester, New York 14642, USA
| | - Marc Parisien
- Department of Computer Science and Operations Research, University of Montreal, Montreal, Quebec H3C CJ7, Canada
| | - Francois Major
- Department of Computer Science and Operations Research, University of Montreal, Montreal, Quebec H3C CJ7, Canada
| | - Douglas H. Turner
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
18
|
Duchardt-Ferner E, Ferner J, Wöhnert J. Direkte Identifizierung von nichtkanonischen RNA-Strukturelementen durch den Nachweis von OH⋅⋅⋅OP-, NH⋅⋅⋅OP- und NH2⋅⋅⋅OP-Wasserstoffbrücken mit NMR-Spektroskopie in Lösung. Angew Chem Int Ed Engl 2011. [DOI: 10.1002/ange.201101743] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
19
|
Duchardt-Ferner E, Ferner J, Wöhnert J. Rapid identification of noncanonical RNA structure elements by direct detection of OH···O=P, NH···O=P, and NH2···O=P hydrogen bonds in solution NMR spectroscopy. Angew Chem Int Ed Engl 2011; 50:7927-30. [PMID: 21837618 DOI: 10.1002/anie.201101743] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2011] [Revised: 05/31/2011] [Indexed: 12/24/2022]
Affiliation(s)
- Elke Duchardt-Ferner
- Institute of Molecular Biosciences, Center for Biomolecular Magnetic Resonance (BMRZ), Johann-Wolfgang-Goethe-Universität Frankfurt, Max-von-Laue-Strasse 9, 60438 Frankfurt, Germany
| | | | | |
Collapse
|
20
|
Moss WN, Eickbush DG, Lopez MJ, Eickbush TH, Turner DH. The R2 retrotransposon RNA families. RNA Biol 2011; 8:714-8. [PMID: 21734471 DOI: 10.4161/rna.8.5.16033] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Analysis of the R2 retrotransposons from multiple silkmoth and fruitfly species have revealed three segments that contain conserved RNA secondary structures. These conserved structures play important roles in the propagation of the R2 element, including R2 RNA processing and transposon integration into the host genome as well as a likely role in translation. Two of the structured regions comprise protein binding sites: one is located in the 3' UTR and the other is in the 5' UTR close to the putative start of the R2 open reading frame (ORF). The 3' structure was deduced from chemical mapping and sequence comparison. The 5' structure was determined using a combination of chemical mapping, oligonucleotide binding, NMR and sequence analysis and contains an unusual pseudoknot structure. The third structure occurs at the 5' end of the R2 RNA and is responsible for self-cleavage of the 5' end of the element from a 28S ribosomal RNA co-transcript. A structure for this fragment was proposed based on motif searching and sequence comparison. There is remarkable similarity in sequence and structure to the hepatitis delta virus (HDV) ribozyme. Seed alignments for the 5' structure and the R2 ribozyme, containing representative sequences and consensus structures, have been submitted to the Rfam database.
Collapse
Affiliation(s)
- Walter N Moss
- Department of Chemistry, University of Rochester, Rochester, NY, USA.
| | | | | | | | | |
Collapse
|
21
|
Andronescu M, Condon A, Hoos HH, Mathews DH, Murphy KP. Computational approaches for RNA energy parameter estimation. RNA (NEW YORK, N.Y.) 2010; 16:2304-18. [PMID: 20940338 PMCID: PMC2995392 DOI: 10.1261/rna.1950510] [Citation(s) in RCA: 65] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
Methods for efficient and accurate prediction of RNA structure are increasingly valuable, given the current rapid advances in understanding the diverse functions of RNA molecules in the cell. To enhance the accuracy of secondary structure predictions, we developed and refined optimization techniques for the estimation of energy parameters. We build on two previous approaches to RNA free-energy parameter estimation: (1) the Constraint Generation (CG) method, which iteratively generates constraints that enforce known structures to have energies lower than other structures for the same molecule; and (2) the Boltzmann Likelihood (BL) method, which infers a set of RNA free-energy parameters that maximize the conditional likelihood of a set of reference RNA structures. Here, we extend these approaches in two main ways: We propose (1) a max-margin extension of CG, and (2) a novel linear Gaussian Bayesian network that models feature relationships, which effectively makes use of sparse data by sharing statistical strength between parameters. We obtain significant improvements in the accuracy of RNA minimum free-energy pseudoknot-free secondary structure prediction when measured on a comprehensive set of 2518 RNA molecules with reference structures. Our parameters can be used in conjunction with software that predicts RNA secondary structures, RNA hybridization, or ensembles of structures. Our data, software, results, and parameter sets in various formats are freely available at http://www.cs.ubc.ca/labs/beta/Projects/RNA-Params.
Collapse
Affiliation(s)
- Mirela Andronescu
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195, USA.
| | | | | | | | | |
Collapse
|
22
|
Liang R, Kierzek E, Kierzek R, Turner DH. Comparisons between chemical mapping and binding to isoenergetic oligonucleotide microarrays reveal unexpected patterns of binding to the Bacillus subtilis RNase P RNA specificity domain. Biochemistry 2010; 49:8155-68. [PMID: 20557101 PMCID: PMC2938832 DOI: 10.1021/bi100286n] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2010] [Revised: 06/16/2010] [Indexed: 11/29/2022]
Abstract
Microarrays with isoenergetic pentamer and hexamer 2'-O-methyl oligonucleotide probes with LNA (locked nucleic acid) and 2,6-diaminopurine substitutions were used to probe the binding sites on the RNase P RNA specificity domain of Bacillus subtilis. Unexpected binding patterns were revealed. Because of their enhanced binding free energies, isoenergetic probes can break short duplexes, merge adjacent loops, and/or induce refolding. This suggests new approaches to the rational design of short oligonucleotide therapeutics but limits the utility of microarrays for providing constraints for RNA structure determination. The microarray results are compared to results from chemical mapping experiments, which do provide constraints. Results from both types of experiments indicate that the RNase P RNA folds similarly in 1 M Na(+) and 10 mM Mg(2+).
Collapse
Affiliation(s)
- Ruiting Liang
- Department of Chemistry, University of Rochester, Rochester, New York 14627
| | - Elzbieta Kierzek
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 60-714 Poznan, Noskowskiego 12/14, Poland
| | - Ryszard Kierzek
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 60-714 Poznan, Noskowskiego 12/14, Poland
| | - Douglas H. Turner
- Department of Chemistry, University of Rochester, Rochester, New York 14627
| |
Collapse
|
23
|
Mathews DH, Moss WN, Turner DH. Folding and finding RNA secondary structure. Cold Spring Harb Perspect Biol 2010; 2:a003665. [PMID: 20685845 DOI: 10.1101/cshperspect.a003665] [Citation(s) in RCA: 106] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Optimal exploitation of the expanding database of sequences requires rapid finding and folding of RNAs. Methods are reviewed that automate folding and discovery of RNAs with algorithms that couple thermodynamics with chemical mapping, NMR, and/or sequence comparison. New functional noncoding RNAs in genome sequences can be found by combining sequence comparison with the assumption that functional noncoding RNAs will have more favorable folding free energies than other RNAs. When a new RNA is discovered, experiments and sequence comparison can restrict folding space so that secondary structure can be rapidly determined with the help of predicted free energies. In turn, secondary structure restricts folding in three dimensions, which allows modeling of three-dimensional structure. An example from a domain of a retrotransposon is described. Discovery of new RNAs and their structures will provide insights into evolution, biology, and design of therapeutics. Applications to studies of evolution are also reviewed.
Collapse
Affiliation(s)
- David H Mathews
- Department of Biochemistry and Biophysics and Center for RNA Biology, University of Rochester School of Medicine and Dentistry, Rochester, New York 14642, USA
| | | | | |
Collapse
|
24
|
Hammond NB, Tolbert BS, Kierzek R, Turner DH, Kennedy SD. RNA internal loops with tandem AG pairs: the structure of the 5'GAGU/3'UGAG loop can be dramatically different from others, including 5'AAGU/3'UGAA. Biochemistry 2010; 49:5817-27. [PMID: 20481618 PMCID: PMC2900907 DOI: 10.1021/bi100332r] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
![]()
Thermodynamic stabilities of 2 × 2 nucleotide tandem AG internal loops in RNA range from −1.3 to +3.4 kcal/mol at 37 °C and are not predicted well with a hydrogen-bonding model. To provide structural information to facilitate development of more sophisticated models for the sequence dependence of stability, we report the NMR solution structures of five RNA duplexes: (rGACGAGCGUCA)2, (rGACUAGAGUCA)2, (rGACAAGUGUCA)2, (rGGUAGGCCA)2, and (rGACGAGUGUCA)2. The structures of these duplexes are compared to that of the previously solved (rGGCAGGCC)2 (Wu, M., SantaLucia, J., Jr., and Turner, D. H. (1997) Biochemistry 36, 4449−4460). For loops bounded by Watson−Crick pairs, the AG and Watson−Crick pairs are all head-to-head imino-paired (cis Watson−Crick/Watson−Crick). The structures suggest that the sequence-dependent stability may reflect non-hydrogen-bonding interactions. Of the two loops bounded by G-U pairs, only the 5′UAGG/3′GGAU loop adopts canonical UG wobble pairing (cis Watson−Crick/Watson−Crick), with AG pairs that are only weakly imino-paired. Strikingly, the 5′GAGU/3′UGAG loop has two distinct duplex conformations, the major of which has both guanosine residues (G4 and G6 in (rGACGAGUGUCA)2) in a syn glycosidic bond conformation and forming a sheared GG pair (G4-G6*, GG trans Watson−Crick/Hoogsteen), both uracils (U7 and U7*) flipped out of the helix, and an AA pair (A5-A5*) in a dynamic or stacked conformation. These structures provide benchmarks for computational investigations into interactions responsible for the unexpected differences in loop free energies and structure.
Collapse
Affiliation(s)
- Nicholas B Hammond
- Department of Chemistry, University of Rochester, Rochester, New York 14627-0216, USA
| | | | | | | | | |
Collapse
|
25
|
Wachowius F, Höbartner C. Chemical RNA modifications for studies of RNA structure and dynamics. Chembiochem 2010; 11:469-80. [PMID: 20135663 DOI: 10.1002/cbic.200900697] [Citation(s) in RCA: 74] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Falk Wachowius
- Research Group Nucleic Acid Chemistry, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Göttingen, Germany
| | | |
Collapse
|
26
|
Quarrier S, Martin JS, Davis-Neulander L, Beauregard A, Laederach A. Evaluation of the information content of RNA structure mapping data for secondary structure prediction. RNA (NEW YORK, N.Y.) 2010; 16:1108-17. [PMID: 20413617 PMCID: PMC2874162 DOI: 10.1261/rna.1988510] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/05/2009] [Accepted: 02/20/2010] [Indexed: 05/21/2023]
Abstract
Structure mapping experiments (using probes such as dimethyl sulfate [DMS], kethoxal, and T1 and V1 RNases) are used to determine the secondary structures of RNA molecules. The process is iterative, combining the results of several probes with constrained minimum free-energy calculations to produce a model of the structure. We aim to evaluate whether particular probes provide more structural information, and specifically, how noise in the data affects the predictions. Our approach involves generating "decoy" RNA structures (using the sFold Boltzmann sampling procedure) and evaluating whether we are able to identify the correct structure from this ensemble of structures. We show that with perfect information, we are always able to identify the optimal structure for five RNAs of known structure. We then collected orthogonal structure mapping data (DMS and RNase T1 digest) under several solution conditions using our high-throughput capillary automated footprinting analysis (CAFA) technique on two group I introns of known structure. Analysis of these data reveals the error rates in the data under optimal (low salt) and suboptimal solution conditions (high MgCl(2)). We show that despite these errors, our computational approach is less sensitive to experimental noise than traditional constraint-based structure prediction algorithms. Finally, we propose a novel approach for visualizing the interaction of chemical and enzymatic mapping data with RNA structure. We project the data onto the first two dimensions of a multidimensional scaling of the sFold-generated decoy structures. We are able to directly visualize the structural information content of structure mapping data and reconcile multiple data sets.
Collapse
Affiliation(s)
- Scott Quarrier
- Biomedical Sciences Program, University at Albany, Albany, New York 12208, USA
| | | | | | | | | |
Collapse
|
27
|
Liu B, Shankar N, Turner DH. Fluorescence competition assay measurements of free energy changes for RNA pseudoknots. Biochemistry 2010; 49:623-34. [PMID: 19921809 PMCID: PMC2808147 DOI: 10.1021/bi901541j] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
![]()
RNA pseudoknots have important functions, and thermodynamic stability is a key to predicting pseudoknots in RNA sequences and to understanding their functions. Traditional methods, such as UV melting and differential scanning calorimetry, for measuring RNA thermodynamics are restricted to temperature ranges around the melting temperature for a pseudoknot. Here, we report RNA pseudoknot free energy changes at 37 °C measured by fluorescence competition assays. Sequence-dependent studies for the loop 1−stem 2 region reveal (1) the individual nearest-neighbor hydrogen bonding (INN-HB) model provides a reasonable estimate for the free energy change when a Watson−Crick base pair in stem 2 is changed, (2) the loop entropy can be estimated by a statistical polymer model, although some penalty for certain loop sequences is necessary, and (3) tertiary interactions can significantly stabilize pseudoknots and extending the length of stem 2 may alter tertiary interactions such that the INN-HB model does not predict the net effect of adding a base pair. The results can inform writing of algorithms for predicting and/or designing RNA secondary structures.
Collapse
Affiliation(s)
- Biao Liu
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | | | | |
Collapse
|
28
|
Reuter JS, Mathews DH. RNAstructure: software for RNA secondary structure prediction and analysis. BMC Bioinformatics 2010; 11:129. [PMID: 20230624 PMCID: PMC2984261 DOI: 10.1186/1471-2105-11-129] [Citation(s) in RCA: 1271] [Impact Index Per Article: 90.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2009] [Accepted: 03/15/2010] [Indexed: 11/16/2022] Open
Abstract
Background To understand an RNA sequence's mechanism of action, the structure must be known. Furthermore, target RNA structure is an important consideration in the design of small interfering RNAs and antisense DNA oligonucleotides. RNA secondary structure prediction, using thermodynamics, can be used to develop hypotheses about the structure of an RNA sequence. Results RNAstructure is a software package for RNA secondary structure prediction and analysis. It uses thermodynamics and utilizes the most recent set of nearest neighbor parameters from the Turner group. It includes methods for secondary structure prediction (using several algorithms), prediction of base pair probabilities, bimolecular structure prediction, and prediction of a structure common to two sequences. This contribution describes new extensions to the package, including a library of C++ classes for incorporation into other programs, a user-friendly graphical user interface written in JAVA, and new Unix-style text interfaces. The original graphical user interface for Microsoft Windows is still maintained. Conclusion The extensions to RNAstructure serve to make RNA secondary structure prediction user-friendly. The package is available for download from the Mathews lab homepage at http://rna.urmc.rochester.edu/RNAstructure.html.
Collapse
Affiliation(s)
- Jessica S Reuter
- Department of Biochemistry & Biophysics and Center for RNA Biology, University of Rochester Medical Center, 601 Elmwood Avenue, Box 712, Rochester, NY 14642, USA
| | | |
Collapse
|
29
|
Puffer B, Kreutz C, Rieder U, Ebert MO, Konrat R, Micura R. 5-Fluoro pyrimidines: labels to probe DNA and RNA secondary structures by 1D 19F NMR spectroscopy. Nucleic Acids Res 2010; 37:7728-40. [PMID: 19843610 PMCID: PMC2794194 DOI: 10.1093/nar/gkp862] [Citation(s) in RCA: 72] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
19F NMR spectroscopy has proved to be a valuable tool to monitor functionally important conformational transitions of nucleic acids. Here, we present a systematic investigation on the application of 5-fluoro pyrimidines to probe DNA and RNA secondary structures. Oligonucleotides with the propensity to adapt secondary structure equilibria were chosen as model systems and analyzed by 1D 19F and 1H NMR spectroscopy. A comparison with the unmodified analogs revealed that the equilibrium characteristics of the bistable DNA and RNA oligonucleotides were hardly affected upon fluorine substitution at C5 of pyrimidines. This observation was in accordance with UV spectroscopic melting experiments which demonstrated that single 5-fluoro substitutions in double helices lead to comparable thermodynamic stabilities. Thus, 5-fluoro pyrimidine labeling of DNA and RNA can be reliably applied for NMR based nucleic acid secondary structure evaluation. Furthermore, we developed a facile synthetic route towards 5-fluoro cytidine phosphoramidites that enables their convenient site-specific incorporation into oligonucleotides by solid-phase synthesis.
Collapse
Affiliation(s)
- Barbara Puffer
- Institute of Organic Chemistry, Center for Molecular Biosciences (CMBI), University of Innsbruck, 6020 Innsbruck, Austria, Laboratory of Organic Chemistry, ETH Zürich, 8093 Zürich, Switzerland and Max Perutz Laboratories, Vienna Biocenter, University of Vienna, 1030 Vienna, Austria
| | | | | | | | | | | |
Collapse
|
30
|
Abstract
RNA pseudoknots are important for function. Three-dimensional structural information is available, insights into factors affecting pseudoknot stability are being reported, and computer programs are available for predicting pseudoknots.
Collapse
Affiliation(s)
- Biao Liu
- Department of Chemistry120 Trustee RoadUniversity of RochesterRochester, NY 14627USA
| | - David H Mathews
- Department of Biochemistry and Biophysics, School of Medicine and Dentistry601 Elmwood AveUniversity of RochesterRochester, NY 14642USA
- Center for RNA Biology, School of Medicine and Dentistry, University of RochesterRochester, NY 14642USA
| | - Douglas H Turner
- Department of Chemistry120 Trustee RoadUniversity of RochesterRochester, NY 14627USA
- Center for RNA Biology, School of Medicine and Dentistry, University of RochesterRochester, NY 14642USA
| |
Collapse
|
31
|
Kierzek E. Binding of short oligonucleotides to RNA: studies of the binding of common RNA structural motifs to isoenergetic microarrays. Biochemistry 2009; 48:11344-56. [PMID: 19835418 DOI: 10.1021/bi901264v] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Binding of short oligonucleotides to RNA is important for many biological processes. On the basis of RNAi phenomena, antisense, and ribozyme approaches, it is useful in the inhibition of biological functions. To be considered as potential therapeutics, oligonucleotides must bind strongly and selectively to a complementary fragment of target RNA. Microarray technologies also involve the binding of oligonucleotide probes to DNA or RNA. Herein, the hybridization of common structural motifs of RNA, i.e., hairpins, internal loops, bulges, 3'- and 5'-dangling ends, and pseudoknots to isoenergetic microarray probes is presented. The analysis demonstrates that microarray probes bind to bulges, internal loops, and dangling ends as expected. Probes may also bind to terminal helixes, however, possibly due to the rearrangement of base pairs. These results suggest that isoenergetic microarray mapping can provide data to facilitate and improve RNA secondary structure prediction. However, optimal results require combination with chemical and/or enzymatic mapping.
Collapse
Affiliation(s)
- Elzbieta Kierzek
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 60-704 Poznan, Noskowskiego 12/14, Poland.
| |
Collapse
|
32
|
Rieder U, Lang K, Kreutz C, Polacek N, Micura R. Evidence for pseudoknot formation of class I preQ1 riboswitch aptamers. Chembiochem 2009; 10:1141-4. [PMID: 19382115 DOI: 10.1002/cbic.200900155] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
All in a knot. The smallest riboswitch forms a pseudoknot in solution. This is demonstrated for preQ(1) class I aptamers by mutational analysis in combination with (1)H NMR-based structure probing. How pseudoknot formation mediates the mRNA response through its expression platform is now open for investigation.
Collapse
Affiliation(s)
- Ulrike Rieder
- Institute of Organic Chemistry, Center for Molecular Biosciences CMBI, 6020 Innsbruck. Austria
| | | | | | | | | |
Collapse
|
33
|
Harmanci AO, Sharma G, Mathews DH. Stochastic sampling of the RNA structural alignment space. Nucleic Acids Res 2009; 37:4063-75. [PMID: 19429694 PMCID: PMC2709569 DOI: 10.1093/nar/gkp276] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
A novel method is presented for predicting the common secondary structures and alignment of two homologous RNA sequences by sampling the 'structural alignment' space, i.e. the joint space of their alignments and common secondary structures. The structural alignment space is sampled according to a pseudo-Boltzmann distribution based on a pseudo-free energy change that combines base pairing probabilities from a thermodynamic model and alignment probabilities from a hidden Markov model. By virtue of the implicit comparative analysis between the two sequences, the method offers an improvement over single sequence sampling of the Boltzmann ensemble. A cluster analysis shows that the samples obtained from joint sampling of the structural alignment space cluster more closely than samples generated by the single sequence method. On average, the representative (centroid) structure and alignment of the most populated cluster in the sample of structures and alignments generated by joint sampling are more accurate than single sequence sampling and alignment based on sequence alone, respectively. The 'best' centroid structure that is closest to the known structure among all the centroids is, on average, more accurate than structure predictions of other methods. Additionally, cluster analysis identifies, on average, a few clusters, whose centroids can be presented as alternative candidates. The source code for the proposed method can be downloaded at http://rna.urmc.rochester.edu.
Collapse
Affiliation(s)
- Arif Ozgun Harmanci
- Department of Electrical and Computer Engineering, University of Rochester, Hopeman 204, RC Box 270126, Rochester, NY 14627, USA
| | | | | |
Collapse
|
34
|
Kierzek E, Christensen SM, Eickbush TH, Kierzek R, Turner DH, Moss WN. Secondary structures for 5' regions of R2 retrotransposon RNAs reveal a novel conserved pseudoknot and regions that evolve under different constraints. J Mol Biol 2009; 390:428-42. [PMID: 19397915 DOI: 10.1016/j.jmb.2009.04.048] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2009] [Revised: 04/17/2009] [Accepted: 04/21/2009] [Indexed: 12/24/2022]
Abstract
Sequences from the 5' region of R2 retrotransposons of four species of silk moth are reported. In Bombyx mori, this region of the R2 messenger RNA contains a binding site for R2 protein and mediates interactions critical to R2 element insertion into the host genome. A model of secondary structure for a segment of this RNA is proposed on the basis of binding to oligonucleotide microarrays, chemical mapping, and comparative sequence analysis. Five conserved secondary structures are identified, including a novel pseudoknot. There is an apparent transition from an entirely RNA structure coding function in most of the 5' segment to a protein coding function near the 3' end. This suggests that local regions evolved under separate functional constraints (structural, coding, or both).
Collapse
Affiliation(s)
- Elzbieta Kierzek
- Department of Chemistry, University of Rochester, NY 14627-0216, USA
| | | | | | | | | | | |
Collapse
|
35
|
Cao S, Chen SJ. Predicting structures and stabilities for H-type pseudoknots with interhelix loops. RNA (NEW YORK, N.Y.) 2009; 15:696-706. [PMID: 19237463 PMCID: PMC2661829 DOI: 10.1261/rna.1429009] [Citation(s) in RCA: 76] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/21/2008] [Accepted: 01/10/2009] [Indexed: 05/20/2023]
Abstract
RNA pseudoknots play a critical role in RNA-related biology from the assembly of ribosome to the regulation of viral gene expression. A predictive model for pseudoknot structure and stability is essential for understanding and designing RNA structure and function. A previous statistical mechanical theory allows us to treat canonical H-type RNA pseudoknots that contain no intervening loop between the helices (see S. Cao and S.J. Chen [2006] in Nucleic Acids Research, Vol. 34; pp. 2634-2652). Biologically significant RNA pseudoknots often contain interhelix loops. Predicting the structure and stability for such more-general pseudoknots remains an unsolved problem. In the present study, we develop a predictive model for pseudoknots with interhelix loops. The model gives conformational entropy, stability, and the free-energy landscape from RNA sequences. The main features of this new model are the computation of the conformational entropy and folding free-energy base on the complete conformational ensemble and rigorous treatment for the excluded volume effects. Extensive tests for the structural predictions show overall good accuracy with average sensitivity and specificity equal to 0.91 and 0.91, respectively. The theory developed here may be a solid starting point for first-principles modeling of more complex, larger RNAs.
Collapse
Affiliation(s)
- Song Cao
- Department of Physics, University of Missouri, Columbia, 65211, USA
| | | |
Collapse
|
36
|
|