Showing papers on "Nucleic acid secondary structure published in 2016"
TL;DR: New genome-wide RNA structure-probing methods are described and the information emerging from their application is synthesized.
Abstract: Single-stranded RNA molecules fold into extraordinarily complicated secondary and tertiary structures as a result of intramolecular base pairing. In vivo, these RNA structures are not static. Instead, they are remodeled in response to changes in the prevailing physicochemical environment of the cell and as a result of intermolecular base pairing and interactions with RNA-binding proteins. Remarkable technical advances now allow us to probe RNA secondary structure at single-nucleotide resolution and genome-wide, both in vitro and in vivo. These data sets provide new glimpses into the RNA universe. Analyses of RNA structuromes in HIV, yeast, Arabidopsis, and mammalian cells and tissues have revealed regulatory effects of RNA structure on messenger RNA (mRNA) polyadenylation, splicing, translation, and turnover. Application of new methods for genome-wide identification of mRNA modifications, particularly methylation and pseudouridylation, has shown that the RNA "epitranscriptome" both influences and is influenced by RNA structure. In this review, we describe newly developed genome-wide RNA structure-probing methods and synthesize the information emerging from their application.
183 citations
TL;DR: The observed imino proton nuclear magnetic resonance resonances and Förster resonance energy transfer efficiencies suggest that m(6)A selectively destabilizes the portion of the hairpin stem where the U5-tract is located, increasing the solvent accessibility of the neighboring bases while maintaining the overall hairpin structure.
160 citations
TL;DR: In vivo click selective 2-hydroxyl acylation and profiling experiment captures RNA secondary structure at a transcriptome-wide level by measuring nucleotide flexibility at base resolution and can reveal in vivo RNA-binding protein imprints or facilitate the dissection of RNA post-transcriptional modifications.
Abstract: icSHAPE (in vivo click selective 2-hydroxyl acylation and profiling experiment) captures RNA secondary structure at a transcriptome-wide level by measuring nucleotide flexibility at base resolution. Living cells are treated with the icSHAPE chemical NAI-N3 followed by selective chemical enrichment of NAI-N3-modified RNA, which provides an improved signal-to-noise ratio compared with similar methods leveraging deep sequencing. Purified RNA is then reverse-transcribed to produce cDNA, with SHAPE-modified bases leading to truncated cDNA. After deep sequencing of cDNA, computational analysis yields flexibility scores for every base across the starting RNA population. The entire experimental procedure can be completed in ∼5 d, and the sequencing and bioinformatics data analysis take an additional 4-5 d with no extensive computational skills required. Comparing in vivo and in vitro icSHAPE measurements can reveal in vivo RNA-binding protein imprints or facilitate the dissection of RNA post-transcriptional modifications. icSHAPE reactivities can additionally be used to constrain and improve RNA secondary structure prediction models.
129 citations
TL;DR: Large-scale combinatorial mutagenesis of unfavorable sequence contexts revealed a consistent pattern whereby mutations that increased motif accessibility improved protein binding and regulatory activity and suggested that mutations that alter sequence context commonly affect RBP binding and regulation.
108 citations
TL;DR: One such method, selective 2'-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq), is described, which measures nucleotide resolution flexibility information for RNAs in vitro and in vivo.
75 citations
TL;DR: The origins, form, and function of RNA secondary structure are reviewed, focusing on plants, and an overview of methods for probing secondary structure is provided, from physical methods such as X-ray crystallography and nuclear magnetic resonance (NMR) imaging to chemical and nuclease probing methods.
Abstract: RNA transcripts fold into secondary structures via intricate patterns of base pairing. These secondary structures impart catalytic, ligand binding, and scaffolding functions to a wide array of RNAs, forming a critical node of biological regulation. Among their many functions, RNA structural elements modulate epigenetic marks, alter mRNA stability and translation, regulate alternative splicing, transduce signals, and scaffold large macromolecular complexes. Thus, the study of RNA secondary structure is critical to understanding the function and regulation of RNA transcripts. Here, we review the origins, form, and function of RNA secondary structure, focusing on plants. We then provide an overview of methods for probing secondary structure, from physical methods such as X-ray crystallography and nuclear magnetic resonance (NMR) imaging to chemical and nuclease probing methods. Combining these latter methods with high-throughput sequencing has enabled them to scale across whole transcriptomes, yielding tremendous new insights into the form and function of RNA secondary structure.
73 citations
TL;DR: In insights into specific structural features that increase the difficulty of finding sequences that fold into a target RNA secondary structure, the Eterna100 benchmark of 100 secondary structure design challenges that span a large range in design difficulty is compiled to help test future efforts.
71 citations
TL;DR: In this article, a computational sequence database mining approach was used to identify two Class 2 CRISPR-Cas systems (subtype VI-B) that lack Cas1 and Cas2 and encompass a single large effector protein, Cas13b, along with one of two previously uncharacterized associated proteins, Csx27 or Csax28.
Abstract: CRISPR-Cas adaptive immune systems defend microbes against foreign nucleic acids via RNA- guided endonucleases. Using a computational sequence database mining approach, we identify two Class 2 CRISPR-Cas systems (subtype VI-B) that lack Cas1 and Cas2 and encompass a single large effector protein, Cas13b, along with one of two previously uncharacterized associated proteins, Csx27 or Csx28. We establish that these CRISPR-Cas systems can achieve RNA interference when heterologously expressed. Through a combination of biochemical and genetic experiments, we show that Cas13b processes its own CRISPR array with short and long direct repeats, cleaves target RNA, and exhibits collateral RNase activity. Using an E. coli essential gene screen, we demonstrate that Cas13b has a double-sided protospacer-flanking sequence and elucidate RNA secondary structure requirements for targeting. We also find that Csx27 represses, whereas Csx28 enhances, Cas13b-mediated RNA interference. Characterization of these CRISPR systems creates opportunities to develop tools to manipulate and monitor cellular transcripts.
71 citations
TL;DR: It is demonstrated that an RNA partition function can be used to calculate the exact probability of formation of hairpin loops, internal loops, bulge loops, or multibranch loops at a given position and that highly probable loops are more likely to be present in the known structure than less probable loops.
Abstract: RNA secondary structure prediction is widely used to analyze RNA sequences. In an RNA partition function calculation, free energy nearest neighbor parameters are used in a dynamic programming algorithm to estimate statistical properties of the secondary structure ensemble. Previously, partition functions have largely been used to estimate the probability that a given pair of nucleotides form a base pair, the conditional stacking probability, the accessibility to binding of a continuous stretch of nucleotides, or a representative sample of RNA structures. Here it is demonstrated that an RNA partition function can also be used to calculate the exact probability of formation of hairpin loops, internal loops, bulge loops, or multibranch loops at a given position. This calculation can also be used to estimate the probability of formation of specific helices. Benchmarking on a set of RNA sequences with known secondary structures indicated that loops that were calculated to be more probable were more likely to be present in the known structure than less probable loops. Furthermore, highly probable loops are more likely to be in the known structure than the set of loops predicted in the lowest free energy structures.
62 citations
TL;DR: A computational approach that leverages the experimental design of the HT-SELEX protocol, RNA secondary structure, and the potential presence of many secondary motifs to identify sequence-structure motifs that show a signature of selection, AptaTRACE can identify low-abundance motifs and could lower HT- SELEX cost and time by reducing the number of selection cycles required.
Abstract: Aptamers, short RNA or DNA molecules that bind distinct targets with high affinity and specificity, can be identified using high-throughput systematic evolution of ligands by exponential enrichment (HT-SELEX), but scalable analytic tools for understanding sequence-function relationships from diverse HT-SELEX data are not available. Here we present AptaTRACE, a computational approach that leverages the experimental design of the HT-SELEX protocol, RNA secondary structure, and the potential presence of many secondary motifs to identify sequence-structure motifs that show a signature of selection. We apply AptaTRACE to identify nine motifs in C-C chemokine receptor type 7 targeted by aptamers in an in vitro cell-SELEX experiment. We experimentally validate two aptamers whose binding required both sequence and structural features. AptaTRACE can identify low-abundance motifs, and we show through simulations that, because of this, it could lower HT-SELEX cost and time by reducing the number of selection cycles required.
53 citations
TL;DR: HyperFold, a multistrand RNA/DNA structure prediction approach for predicting nucleic acid complexes that can contain pseudoknots, is presented and it is shown that HyperFold also performs competitively compared to other published folding algorithms.
Abstract: RNA is an attractive material for the creation of molecular logic gates that release programmed functionalities only in the presence of specific molecular interaction partners. Here we present HyperFold, a multistrand RNA/DNA structure prediction approach for predicting nucleic acid complexes that can contain pseudoknots. We show that HyperFold also performs competitively compared to other published folding algorithms. We performed a large variety of RNA/DNA hybrid reassociation experiments for different concentrations, DNA toehold lengths, and G+C content and find that the observed tendencies for reassociation correspond well to computational predictions. Importantly, we apply this method to the design and experimental verification of a two-stranded RNA molecular switch that upon binding to a single-stranded RNA toehold disease-marker trigger mRNA changes its conformation releasing an shRNA-like Dicer substrate structure. To demonstrate the concept, connective tissue growth factor (CTGF) mRNA and enhance...
TL;DR: The first comprehensive comparison comprising 14 methods that predict general intermolecular basepairs is provided, finding that maintaining high accuracy across biologically different data sets and increasing input lengths remains a huge challenge, causing implications for de novo transcriptome-wide searches.
Abstract: RNA-RNA interactions are fast emerging as a major functional component in many newly discovered non-coding RNAs. Basepairing is believed to be a major contributor to the stability of these intermolecular interactions, much like intramolecular basepairs formed in RNA secondary structure. As such, using algorithms similar to those for predicting RNA secondary structure, computational methods have been recently developed for the prediction of RNA-RNA interactions. We provide the first comprehensive comparison comprising 14 methods that predict general intermolecular basepairs. To evaluate these, we compile an extensive data set of 54 experimentally confirmed fungal snoRNA-rRNA interactions and 102 bacterial sRNA-mRNA interactions. We test the performance accuracy of all methods, evaluating the effects of tool settings, sequence length, and multiple sequence alignment usage and quality. Our results show that-unlike for RNA secondary structure prediction--the overall best performing tools are non-comparative energy-based tools utilizing accessibility information that predict short interactions on this data set. Furthermore, we find that maintaining high accuracy across biologically different data sets and increasing input lengths remains a huge challenge, causing implications for de novo transcriptome-wide searches. Finally, we make our interaction data set publicly available for future development and benchmarking efforts.
TL;DR: RBP-Var provides an easy-to-use web interface that allows users to rapidly find whether SNVs of interest can transform the secondary structure of RNA and identify RBPs whose binding may be subsequently disrupted, and is useful in selecting candidate SNVs for further functional studies and exploring causal SNVs underlying human diseases.
Abstract: Transcription factors bind to the genome by forming specific contacts with the primary DNA sequence; however, RNA-binding proteins (RBPs) have greater scope to achieve binding specificity through the RNA secondary structure. It has been revealed that single nucleotide variants (SNVs) that alter RNA structure, also known as RiboSNitches, exhibit 3-fold greater local structure changes than replicates of the same DNA sequence, demonstrated by the fact that depletion of RiboSNitches could result in the alteration of specific RNA shapes at thousands of sites, including 3' UTRs, binding sites of microRNAs and RBPs. However, the network between SNVs and post-transcriptional regulation remains unclear. Here, we developed RBP-Var, a database freely available at http://www.rbp-var.biols.ac.cn/, which provides annotation of functional variants involved in post-transcriptional interaction and regulation. RBP-Var provides an easy-to-use web interface that allows users to rapidly find whether SNVs of interest can transform the secondary structure of RNA and identify RBPs whose binding may be subsequently disrupted. RBP-Var integrates DNA and RNA biology to understand how various genetic variants and post-transcriptional mechanisms cooperate to orchestrate gene expression. In summary, RBP-Var is useful in selecting candidate SNVs for further functional studies and exploring causal SNVs underlying human diseases.
TL;DR: A global analysis of diverse virus types for the impact of nucleotide preference on the predicted structure of the RNA genome that is packaged in virion particles and on the codon usage in the viral open reading frames finds the nucleotide bias was the major determinant of the virus-specific codon usages, thus limiting a role for codon selection and translational control.
Abstract: We are interested in the influence of nucleotide composition on the fundamental characteristics of the virus RNA genome. Most RNA viruses have genomes with a distinct nucleotide composition, e.g. ranging from minimally 12.9 % to maximally 40.3 % (C- and U-count, respectively, in coronavirus HKU). We present a global analysis of diverse virus types, including plus-strand, minus-strand and double-strand RNA viruses, for the impact of this nucleotide preference on the predicted structure of the RNA genome that is packaged in virion particles and on the codon usage in the viral open reading frames. Several virus-specific features will be described, but also some general conclusions were drawn. Without exception, the virus-specific nucleotide bias was enriched in the unpaired, single-stranded regions of the RNA genome, thus creating an even more striking virus-specific signature. We present a simple mechanism that is based on elementary aspects of RNA structure folding to explain this general trend. In general, the nucleotide bias was the major determinant of the virus-specific codon usages, thus limiting a role for codon selection and translational control. We will discuss molecular and evolutionary scenarios that may be responsible for the diverse nucleotide biases of RNA viruses.
TL;DR: The results suggest that the conserved secondary structures predicted in the M segment are involved in the production of infectious viral particles during IAV replication.
Abstract: As well as encoding viral proteins, genomes of RNA viruses harbor secondary and tertiary RNA structures that have been associated with functions essential for successful replication and propagation. Here, we identified stem-loop structures that are extremely conserved among 1,884 M segment sequences of influenza A virus (IAV) strains from various subtypes and host species using computational and evolutionary methods. These structures were predicted within the 3' and 5' ends of the coding regions of M1 and M2, respectively, where packaging signals have been previously proposed to exist. These signals are thought to be required for the incorporation of a single copy of 8 different negative-strand RNA segments (vRNAs) into an IAV particle. To directly test the functionality of conserved stem-loop structures, we undertook reverse genetic experiments to introduce synonymous mutations designed to disrupt secondary structures predicted at 3 locations and found them to attenuate infectivity of recombinant virus. In one mutant, predicted to disrupt stem loop structure at nucleotide positions 219-240, attenuation was more evident at increased temperature and was accompanied by an increase in the production of defective virus particles. Our results suggest that the conserved secondary structures predicted in the M segment are involved in the production of infectious viral particles during IAV replication.
TL;DR: RNAComposer is a user-friendly and freely available server for 3D structure prediction of RNA up to 500 nucleotide residues and this method allows incorporation of distance restraints derived from the experimental data to strengthen the 3D predictions.
TL;DR: This work presents the first experimentally-derived exploration of potential secondary structure in an influenza A naked (protein-free) genomic segment, suggesting that vRNA self-folding may generate helixes and loops that are important at one or more stages of the influenza life cycle.
Abstract: Influenza A is a negative sense RNA virus that kills hundreds of thousands of humans each year. Base pairing in RNA is very favorable, but possibilities for RNA secondary structure of the influenza genomic RNA have not been investigated. This work presents the first experimentally-derived exploration of potential secondary structure in an influenza A naked (protein-free) genomic segment. Favorable folding regions are revealed by in vitro chemical structure mapping, thermodynamics, bioinformatics, and binding to isoenergetic microarrays of an entire natural sequence of the 875 nt segment 8 vRNA and of a smaller fragment. Segment 8 has thermodynamically stable and evolutionarily conserved RNA structure and encodes essential viral proteins NEP and NS1. This suggests that vRNA self-folding may generate helixes and loops that are important at one or more stages of the influenza life cycle.
TL;DR: The novel software, ParasoR, is designed to run on a computer cluster and enables the exact computation of various structural features of long RNA sequences under the constraint of maximal base-pairing distance and is shown to be very useful for analyzing the structural properties of longRNA sequences.
Abstract: RNA secondary structure around splice sites is known to assist normal splicing by promoting spliceosome recognition. However, analyzing the structural properties of entire intronic regions or pre-mRNA sequences has been difficult hitherto, owing to serious experimental and computational limitations, such as low read coverage and numerical problems. Our novel software, “ParasoR”, is designed to run on a computer cluster and enables the exact computation of various structural features of long RNA sequences under the constraint of maximal base-pairing distance. ParasoR divides dynamic programming (DP) matrices into smaller pieces, such that each piece can be computed by a separate computer node without losing the connectivity information between the pieces. ParasoR directly computes the ratios of DP variables to avoid the reduction of numerical precision caused by the cancellation of a large number of Boltzmann factors. The structural preferences of mRNAs computed by ParasoR shows a high concordance with those determined by high-throughput sequencing analyses. Using ParasoR, we investigated the global structural preferences of transcribed regions in the human genome. A genome-wide folding simulation indicated that transcribed regions are significantly more structural than intergenic regions after removing repeat sequences and k-mer frequency bias. In particular, we observed a highly significant preference for base pairing over entire intronic regions as compared to their antisense sequences, as well as to intergenic regions. A comparison between pre-mRNAs and mRNAs showed that coding regions become more accessible after splicing, indicating constraints for translational efficiency. Such changes are correlated with gene expression levels, as well as GC content, and are enriched among genes associated with cytoskeleton and kinase functions. We have shown that ParasoR is very useful for analyzing the structural properties of long RNA sequences such as mRNAs, pre-mRNAs, and long non-coding RNAs whose lengths can be more than a million bases in the human genome. In our analyses, transcribed regions including introns are indicated to be subject to various types of structural constraints that cannot be explained from simple sequence composition biases. ParasoR is freely available at https://github.com/carushi/ParasoR
.
TL;DR: The results indicate that RNA editing is a quantitative trait subject to genetic variation and that evolutionary constraints have given rise to distinct genetic architectures in the two canonical types of RNA editing.
Abstract: RNA editing refers to post-transcriptional processes that alter the base sequence of RNA. Recently, hundreds of new RNA editing targets have been reported. However, the mechanisms that determine the specificity and degree of editing are not well understood. We examined quantitative variation of site-specific editing in a genetically diverse multiparent population, Diversity Outbred mice, and mapped polymorphic loci that alter editing ratios globally for C-to-U editing and at specific sites for A-to-I editing. An allelic series in the C-to-U editing enzyme Apobec1 influences the editing efficiency of Apob and 58 additional C-to-U editing targets. We identified 49 A-to-I editing sites with polymorphisms in the edited transcript that alter editing efficiency. In contrast to the shared genetic control of C-to-U editing, most of the variable A-to-I editing sites were determined by local nucleotide polymorphisms in proximity to the editing site in the RNA secondary structure. Our results indicate that RNA editing is a quantitative trait subject to genetic variation and that evolutionary constraints have given rise to distinct genetic architectures in the two canonical types of RNA editing.
TL;DR: It is found that the 5-formyl group could increase duplex thermal stability and enhance base pairing specificity and increase base stacking between the f5C and the neighboring bases while not causing significant global and local structure perturbations.
Abstract: 5-Formylcytidine (f(5)C), a previously discovered natural nucleotide in the mitochondrial tRNA of many species including human, has been recently detected as the oxidative product of 5-methylcytidine (m(5)C) through 5-hydroxymethylcytidine (hm(5)C) in total RNA of mammalian cells. The discovery indicated that these cytosine derivatives in RNA might also play important epigenetic roles similar as in DNA, which has been intensively investigated in the past few years. In this paper, we studied the base pairing specificity of f(5)C in different RNA duplex contexts. We found that the 5-formyl group could increase duplex thermal stability and enhance base pairing specificity. We present three high-resolution crystal structures of an octamer RNA duplex [5'-GUA(f(5)C)GUAC-3']2 that have been solved under three crystallization conditions with different buffers and pH values. Our results showed that the 5-formyl group is located in the same plane as the cytosine base and forms an intra-residue hydrogen bond with the amino group in the N4 position. In addition, this modification increases the base stacking between the f(5)C and the neighboring bases while not causing significant global and local structure perturbations. This work provides insights into the effects of 5-formylcytosine on RNA duplex.
TL;DR: Structural Surfer provides an easy-to-use database and visualization interface for allowing users to interrogate the currently available transcriptome-wide RNA secondary structure information for mammals.
Abstract: Background
RNA molecules fold into complex three-dimensional shapes, guided by the pattern of hydrogen bonding between nucleotides. This pattern of base pairing, known as RNA secondary structure, is critical to their cellular function. Recently several diverse methods have been developed to assay RNA secondary structure on a transcriptome-wide scale using high-throughput sequencing. Each approach has its own strengths and caveats, however there is no widely available tool for visualizing and comparing the results from these varied methods.
TL;DR: BEAM (BEAr Motif finder), a novel method for structural motif discovery from a set of unaligned RNAs, taking advantage of a recently developed encoding for RNA secondary structure named BEAR (Brand nEw Alphabet for RNAs) and of evolutionary substitution rates of secondary structure elements.
Abstract: Functional RNA regions are often related to recurrent secondary structure patterns (or motifs), which can exert their role in several different ways, particularly in dictating the interaction with RNA-binding proteins, and acting in the regulation of a large number of cellular processes. Among the available motif-finding tools, the majority focuses on sequence patterns, sometimes including secondary structure as additional constraints to improve their performance. Nonetheless, secondary structures motifs may be concurrent to their sequence counterparts or even encode a stronger functional signal. Current methods for searching structural motifs generally require long pipelines and/or high computational efforts or previously aligned sequences. Here, we present BEAM (BEAr Motif finder), a novel method for structural motif discovery from a set of unaligned RNAs, taking advantage of a recently developed encoding for RNA secondary structure named BEAR (Brand nEw Alphabet for RNAs) and of evolutionary substitution rates of secondary structure elements. Tested in a varied set of scenarios, from small- to large-scale, BEAM is successful in retrieving structural motifs even in highly noisy data sets, such as those that can arise in CLIP-Seq or other high-throughput experiments.
TL;DR: Surface plasmon resonance (SPR) binding analysis suggests that the accessibility to G4 recognition elements is different among the three structures and supports the idea that G4 ligands might be shaped to achieve structure selectivity in a biological context.
Abstract: G-rich RNA and DNA oligonucleotides derived from the human telomeric sequence were assembled onto addressable cyclopeptide platforms through oxime ligations and copper-catalyzed azide-alkyne cycloaddition (CuAAc) reactions. The resulting conjugates were able to fold into highly stable RNA and DNA:RNA hybrid G-quadruplex (G4) architectures as demonstrated by UV, circular dichroism (CD), and NMR spectroscopic analysis. Whereas rationally designed parallel RNA and DNA:RNA hybrid G4 topologies could be obtained, we could not force the formation of an antiparallel RNA G4 structure, thus supporting the idea that this topology is strongly disfavored. The binding affinities of four representative G4 ligands toward the discrete RNA and DNA:RNA hybrid G4 topologies were compared to the one obtained with the corresponding DNA G4 structure. Surface plasmon resonance (SPR) binding analysis suggests that the accessibility to G4 recognition elements is different among the three structures and supports the idea that G4 ligands might be shaped to achieve structure selectivity in a biological context.
TL;DR: The RAG resource for large tree graphs up to 13 vertices is expanded, thereby cataloging more than 10 times as many secondary structures and applying clustering algorithms based on features of RNA secondary structures translated from known tertiary structures to suggest which hypothetical large RNA motifs can be considered "RNA-like".
TL;DR: NMR was used to investigate destabilization of the full-length TAR (trans activating response element) RNA by NC, which is involved in the first strand transfer step of reverse transcription and significantly accelerated the formation of the initial complex at a ratio of 1:59.
Abstract: The mature HIV-1 nucleocapsid protein NCp7 (NC) plays a key role in reverse transcription facilitating the two obligatory strand transfers. Several properties contribute to its efficient chaperon activity: preferential binding to single-stranded regions, nucleic acid aggregation, helix destabilization, and rapid dissociation from nucleic acids. However, little is known about the relationships between these different properties, which are complicated by the ability of the protein to recognize particular HIV-1 stem-loops, such as SL1, SL2, and SL3, with high affinity and without destabilizing them. These latter properties are important in the context of genome packaging, during which NC is part of the Gag precursor. We used NMR to investigate destabilization of the full-length TAR (trans activating response element) RNA by NC, which is involved in the first strand transfer step of reverse transcription. NC was used at a low protein:nucleotide (nt) ratio of 1:59 in these experiments. NMR data for the imino protons of TAR identified most of the base pairs destabilized by NC. These base pairs were adjacent to the loops in the upper part of the TAR hairpin rather than randomly distributed. Gel retardation assays showed that conversion from the initial TAR-cTAR complex to the fully annealed form occurred much more slowly at the 1:59 ratio than at the higher ratios classically used. Nevertheless, NC significantly accelerated the formation of the initial complex at a ratio of 1:59.
TL;DR: A probabilistic model is learned from bona fide binding sites determined by CLIP and applied genome wide to generate high specificity binding site predictions and leverages motif conservation, RNA secondary structure, and the tendency of RBP binding sites to cluster together.
Abstract: RNA-binding proteins (RBPs) are critical components of post-transcriptional gene expression regulation. However, their binding sites have until recently been difficult to determine due to the apparent low specificity of RBPs for their target transcripts and the lack of high-throughput assays for analyzing binding sites genome wide. Here we present a bioinformatics method for predicting RBP binding motif sites on a genome-wide scale that leverages motif conservation, RNA secondary structure, and the tendency of RBP binding sites to cluster together. A probabilistic model is learned from bona fide binding sites determined by CLIP and applied genome wide to generate high specificity binding site predictions.
TL;DR: In this paper, the conformation of the ELVd−RNAs was determined in silico and in vitro using three different softwares, and the results were essentially consistent among themselves.
Abstract: Eggplant latent viroid (ELVd), like other members of family Avsunviroidae, replicates in plastids through a symmetric rolling-circle mechanism in which elongation of RNA strands is most likely catalyzed by a nuclear-encoded polymerase (NEP) translocated to plastids. Here we have addressed where NEP initiates transcription of viroid strands. Because this step is presumably directed by sequence/structural motifs, we have previously determined the conformation of the monomeric linear (+) and (-) RNAs of ELVd resulting from hammerhead-mediated self-cleavage. In silico predictions with 3 softwares led to similar bifurcated conformations for both ELVd strands. In vitro examination by non-denaturing PAGE showed that they migrate as prominent single bands, with the ELVd (+) RNA displaying a more compact conformation as revealed by its faster electrophoretic mobility. In vitro SHAPE analysis corroborated the ELVd conformations derived from thermodynamics-based predictions in silico. Moreover, sequence analysis of 94 full-length natural ELVd variants disclosed co-variations, and mutations converting canonical into wobble pairs or vice versa, which confirmed in vivo most of the stems predicted in silico and in vitro, and additionally helped to introduce minor structural refinements. Therefore, results from the 3 experimental approaches were essentially consistent among themselves. Application to RNA preparations from ELVd-infected tissue of RNA ligase-mediated rapid amplification of cDNA ends, combined with pretreatments to modify the 5' ends of viroid strands, mapped the transcription initiation sites of ELVd (+) and (-) strands in vivo at different sequence/structural motifs, in contrast with the situation previously observed in 2 other members of the family Avsunviroidae.
TL;DR: A set of papers within this issue explore a diversity of new strategies for solving, modeling and designing RNA structures and many of them expand upon a set of working themes that were generated by the Telluride meeting group.
TL;DR: DNA intercalation has been very useful for engineering DNA-based functional materials and its application in materials engineering is still being investigated.
Abstract: DNA intercalation has been very useful for engineering DNA-based functional materials. It is generally expected that the intercalation phenomenon in RNA would be similar to that in DNA. Here we note that the neighbor-exclusion principle is violated in RNA by naphthalene-based cationic probes, in contrast to the fact that it is usually valid in DNA. All the intercalation structures are responsible for the fluorescence, where small naphthalene moieties are intercalated in between bases via π-π interactions. The structure is aided by hydrogen bonds between the cationic moieties and the ribose-phosphate backbone, which results in specific selectivity for RNA over DNA. This experimentally observed mechanism is supported by computationally reproducing the fluorescence and CD data. MD simulations confirm the unfolding of RNA due to the intercalation of probes. Elucidation of the mechanism of selective sensing for RNA over DNA would be highly beneficial for dynamical observation of RNA which is essential for studying its biological roles.
TL;DR: Comparison of the models with other RNA structural elements, as well as similarity searches of human transcriptome and viral genomes demonstrate that SLα and SLβ are specific to HBV transcripts, however, they are well conserved among the hepadnaviruses of non-human primates.
Abstract: Many viruses contain RNA elements that modulate splicing and/or promote nuclear export of their RNAs. The RNAs of the major human pathogen, hepatitis B virus (HBV) contain a large (~600 bases) composite cis-acting 'post-transcriptional regulatory element' (PRE). This element promotes expression from these naturally intronless transcripts. Indeed, the related woodchuck hepadnavirus PRE (WPRE) is used to enhance expression in gene therapy and other expression vectors. These PRE are likely to act through a combination of mechanisms, including promotion of RNA nuclear export. Functional components of both the HBV PRE and WPRE are 2 conserved RNA cis-acting stem-loop (SL) structures, SLα and SLβ. They are within the coding regions of polymerase (P) gene, and both P and X genes, respectively. Based on previous studies using mutagenesis and/or nuclear magnetic resonance (NMR), here we propose 2 covariance models for SLα and SLβ. The model for the 30-nucleotide SLα contains a G-bulge and a CNGG(U) apical loop of which the first and the fourth loop residues form a CG pair and the fifth loop residue is bulged out, as observed in the NMR structure. The model for the 23-nucleotide SLβ contains a 7-base-pair stem and a 9-nucleotide loop. Comparison of the models with other RNA structural elements, as well as similarity searches of human transcriptome and viral genomes demonstrate that SLα and SLβ are specific to HBV transcripts. However, they are well conserved among the hepadnaviruses of non-human primates, the woodchuck and ground squirrel.