Universal and strain specific structure features of segment 8 genomic RNA of influenza A virus—application of 4-thiouridine photocrosslinking

RNA structure in the influenza A virus (IAV) has been the focus of several studies that have shown connections between conserved secondary structure motifs and their biological function in the virus replication cycle. Questions have arisen on how to best recognize and understand the pandemic properties of IAV strains from an RNA perspective, but determination of the RNA secondary structure has been challenging. Herein, we used chemical mapping to determine the secondary structure of segment 8 viral RNA (vRNA) of the pandemic A/California/04/2009 (H1N1) strain of IAV. Additionally, this long, naturally occurring RNA served as a model to evaluate RNA mapping with 4-thiouridine (4sU) crosslinking. We explored 4-thiouridine as a probe of nucleotides in close proximity, through its incorporation into newly transcribed RNA and subsequent photoactivation. RNA secondary structural features both universal to type A strains and unique to the A/California/04/2009 (H1N1) strain were recognized. 4sU mapping confirmed and facilitated RNA structure prediction, according to several rules: 4sU photocross-linking forms efficiently in the double-stranded region of RNA with some flexibility, in the ends of helices, and across bulges and loops when their structural mobility is permitted. This method highlighted three-dimensional properties of segment 8 vRNA secondary structure motifs and allowed to propose several long-range three-dimensional interactions. 4sU mapping combined with chemical mapping and bioinformatic analysis could be used to enhance the RNA structure determination as well as recognition of target regions for antisense strategies or viral RNA detection.     

Changes in 7SL RNA conformation during the signal recognition particle cycle.

The structure of 7SL RNA has been probed by chemical modification followed by primer extension, using four substrates: (i) naked 7SL RNA; (ii) free signal recognition particle (SRP); (iii) polysome bound SRP; and (iv) membrane bound SRP. Decreasing sensitivity to chemical modification between these different substrates suggests regions on 7SL RNA that: bind proteins associated with SRP might interact with ribosomes; and are protected by binding to membranes. Other areas increase in chemical sensitivity, exemplified by a tertiary interaction present in naked 7SL RNA but not in free SRP. Such changes suggest that 7SL RNA changes its conformation during the SRP cycle. These conformational changes could be a necessary component to move through the SRP cycle from one stage to the next.     

A microcomupter program for the computation of free energy of the secondary structure of RNA oligomers by Ninio's rules

An interactive microcomputer program for the rapid computation of the free (ΔG) of the secondary structure of RNA molecules is presented. The program assigns free energies (in kcal/mol) to helices; bulge, internal, and multi-branch loops; hairpins; unparied, and G:C termini according to Ninio's rules, and displays a running total during computation. It is written in ‘Microsoft Basic’ and is apllicable to virtually any Basic system with no modification. The program enables rapid bench-top determinations of free energies attributable to secondary structure features of RNA oligomers, and eliminates both the tedium and risk of error associated with manual calculations.     

Detection of periodic patterns in RNA sequences: the first encapsidated region of the TMV RNA

A method is developed to study the periodic properties of nucleotide sequences allowing the favoured pattern of the repeating unit, as well as the length and localization of this periodic segment to be determined simultaneously. The degree of periodicity is evaluated calculating the probabilities for random occurrence of the maximal deviations of the nucleotide composition in each phase, making use of the binomial formula.The nucleotide sequence of the tobacco mosaic virus (TMV) RNA responsible for recognition of the homologous protein (“assembly origin”, AO) (Zimmern & Butler, 1977) was investigated in order to find periodic regions of primary structure which might be essential in the recognition process. As a result the most periodic segments of the AO consisting of 31 and 17 nucleotides corresponding to the schemes GAU or GA1 have been found. However, the periodicities in these regions do not exceed that expected for random sequences. It can be considered as an evidence that in addition to peculiarities of primary structure, some other features such as RNA secondary or tertiary structure are essential in this interaction.For comparison the nucleotide sequences of the other fragments of TMV RNA as well as MS2 RNA, TYMV RNA, 16S rRNA and phage fd DNA were investigated by the same method.     

The structure of a rigorously conserved RNA element within the SARS virus genome

We have solved the three-dimensional crystal structure of the stem-loop II motif (s2m) RNA element of the SARS virus genome to 2.7-Å resolution. SARS and related coronaviruses and astroviruses all possess a motif at the 3′ end of their RNA genomes, called the s2m, whose pathogenic importance is inferred from its rigorous sequence conservation in an otherwise rapidly mutable RNA genome. We find that this extreme conservation is clearly explained by the requirement to form a highly structured RNA whose unique tertiary structure includes a sharp 90° kink of the helix axis and several novel longer-range tertiary interactions. The tertiary base interactions create a tunnel that runs perpendicular to the main helical axis whose interior is negatively charged and binds two magnesium ions. These unusual features likely form interaction surfaces with conserved host cell components or other reactive sites required for virus function. Based on its conservation in viral pathogen genomes and its absence in the human genome, we suggest that these unusual structural features in the s2m RNA element are attractive targets for the design of anti-viral therapeutic agents. Structural genomics has sought to deduce protein function based on three-dimensional homology. Here we have extended this approach to RNA by proposing potential functions for a rigorously conserved set of RNA tertiary structural interactions that occur within the SARS RNA genome itself. Based on tertiary structural comparisons, we propose the s2m RNA binds one or more proteins possessing an oligomer-binding-like fold, and we suggest a possible mechanism for SARS viral RNA hijacking of host protein synthesis, both based upon observed s2m RNA macromolecular mimicry of a relevant ribosomal RNA fold.     

SHAMS: Combining chemical modification of RNA with mass spectrometry to examine polypurine tract-containing RNA/DNA hybrids

Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) has gained popularity as a facile method of examining RNA structure both in vitro and in vivo, exploiting accessibility of the ribose 2′-OH to acylation by N-methylisatoic anhydride (NMIA) in unpaired or flexible configurations. Subsequent primer extension terminates at the site of chemical modification, and these products are fractionated by high-resolution gel electrophoresis. When applying SHAPE to investigate structural features associated with the wild-type and analog-substituted polypurine tract (PPT)–containing RNA/DNA hybrids, their size (20–25 base pairs) rendered primer extension impractical. As an alternative method of detection, we reasoned that chemical modification could be combined with tandem mass spectrometry, relying on the mass increment of RNA fragments containing the NMIA adduct (M_r = 133 Da). Using this approach, we demonstrate both specific modification of the HIV-1 PPT RNA primer and variations in its acylation pattern induced by replacing template nucleotides with a non-hydrogen-bonding thymine isostere. Our selective 2′-hydroxyl acylation analyzed by mass spectrometry strategy (SHAMS) should find utility when examining the structure of small RNA fragments or RNA/DNA hybrids where primer extension cannot be performed.     

An RNA Phage Lab: MS2 in Walter Fiers’ Laboratory of Molecular Biology in Ghent, from Genetic Code to Gene and Genome, 1963–1976

The importance of viruses as model organisms is well-established in molecular biology and Max Delbrück’s phage group set standards in the DNA phage field. In this paper, I argue that RNA phages, discovered in the 1960s, were also instrumental in the making of molecular biology. As part of experimental systems, RNA phages stood for messenger RNA (mRNA), genes and genome. RNA was thought to mediate information transfers between DNA and proteins. Furthermore, RNA was more manageable at the bench than DNA due to the availability of specific RNases, enzymes used as chemical tools to analyse RNA. Finally, RNA phages provided scientists with a pure source of mRNA to investigate the genetic code, genes and even a genome sequence. This paper focuses on Walter Fiers’ laboratory at Ghent University (Belgium) and their work on the RNA phage MS2. When setting up his Laboratory of Molecular Biology, Fiers planned a comprehensive study of the virus with a strong emphasis on the issue of structure. In his lab, RNA sequencing, now a little-known technique, evolved gradually from a means to solve the genetic code, to a tool for completing the first genome sequence. Thus, I follow the research pathway of Fiers and his ‘RNA phage lab’ with their evolving experimental system from 1960 to the late 1970s. This study illuminates two decisive shifts in post-war biology: the emergence of molecular biology as a discipline in the 1960s in Europe and of genomics in the 1990s.     

Fingerprinting the junctions of RNA structure by an open-paddlewheel diruthenium compound

RNA function is determined by its structural organization. The RNA structure consists of the combination of distinct secondary structure motifs connected by junctions that play an essential role in RNA folding. Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) probing is an established methodology to analyze the secondary structure of long RNA molecules in solution, which provides accurate data about unpaired nucleotides. However, the residues located at the junctions of RNA structures usually remain undetected. Here we report an RNA probing method based on the use of a novel open-paddlewheel diruthenium (OPW-Ru) compound [Ru₂Cl₂(µ-DPhF)₃(DMSO)] (DPhF = N,N′-diphenylformamidinate). This compound has four potential coordination sites in a singular disposition to establish covalent bonds with substrates. As a proof of concept, we have analyzed the reactivity of OPW-Ru toward RNA using two viral internal ribosome entry site (IRES) elements whose function depends on the structural organization of the molecule. Our study suggests that the compound OPW-Ru preferentially attacks at positions located one or two nucleotides away from junctions or bulges of the RNA structure. The OPW-Ru fingerprinting data differ from that obtained by other chemical reagents and provides new information about RNA structure features.     

Structural insight into antisense gapmer-RNA oligomer duplexes through molecular dynamics simulations

There is an extensive research carrying out on antisense technology and the molecules entering into clinical trials are increasing rapidly. Phosphorothioate (PS) is a chemical modification in which nonbridged oxygen is replaced with a sulfur, consequently providing resistance against nuclease activity. The 2'-4' conformationally restricted nucleoside has the structural features of both 2'-O-methoxy ethyl RNA (MOE), which shows good toxicity profile, and locked nucleic acid (LNA), which shows good binding affinity towards the target RNA. These modifications have been studied and suggested that they can be a potential therapeutic agents in antisense therapy. Mipomersen (ISIS 301012), which contains the novel nucleoside modification has been used to target to apolipoprotein (Apo B), which reduces LDL cholesterol by 6–41%. In this study, classical molecular dynamics (MD) simulations were performed on six different antisense gapmer/target-RNA oligomer duplexes (LNA-PS-LNA/RNA, RcMOE-PS-RcMOE/RNA, ScMOE-PS-ScMOE/RNA, MOE-PS-MOE/RNA, PS-DNA/RNA and DNA/RNA) to investigate the structural dynamics, stability and solvation properties. The LNA, MOE nucleotides present in respective duplexes are showing the structure of A-form and the PS-DNA nucleotides resemble the structure of B-form helix with respect to some of the helical parameters. Free energy calculations suggest that the oligomer, which contains LNA binds to the RNA strongly than other modifications as shown in experimental results. The MOE modified nucleotide, which although had a lower binding affinity but higher solvent accessible surface area (SASA) compared to the other modifications, may be influencing the toxicity and hence may be used it in Mipomersen, the second antisense molecule which is approved by FDA.

Communicated by Ramaswamy H. Sarma     

A model system for investigating lineshape/structure correlations in RNA site-directed spin labeling

In RNA site-directed spin labeling (SDSL) studies, structural and dynamic information at the individual RNA nucleotide level is derived from the observed electron paramagnetic resonance spectrum of a covalently attached nitroxide. A systematic approach for RNA SDSL is to establish a library that categorizes observed spectral lineshapes based on known RNA structures, thus enabling lineshape-based structure identification at any RNA site. To establish the first RNA SDSL library, selective secondary structure elements have been systematically engineered into a model RNA. Nitroxide lineshapes reporting features specific to each element were obtained utilizing a new avidin-tethering scheme for suppressing spectral effects due to uniform RNA tumbling. The data demonstrated two key features required for a SDSL library with a predicting power: (i) spectral divergence--distinctive lineshape for different elements; and (ii) spectral convergence--similar lineshape for the same element in different contexts. This sets the foundation for further RNA SDSL library development.     

Energy landscape of k-point mutants of an RNA molecule

MOTIVATION: A k-point mutant of a given RNA sequence s = s(1), ..., s(n) is an RNA sequence s' = s'(1),..., s'(n) obtained by mutating exactly k-positions in s; i.e. Hamming distance between s and s' equals k. To understand the effect of pointwise mutation in RNA, we consider the distribution of energies of all secondary structures of k-point mutants of a given RNA sequence. RESULTS: Here we describe a novel algorithm to compute the mean and standard deviation of energies of all secondary structures of k-point mutants of a given RNA sequence. We then focus on the tail of the energy distribution and compute, using the algorithm AMSAG, the k-superoptimal structure; i.e. the secondary structure of a < or =k-point mutant having least free energy over all secondary structures of all k'-point mutants of a given RNA sequence, for k' < or = k. Evidence is presented that the k-superoptimal secondary structure is often closer, as measured by base pair distance and two additional distance measures, to the secondary structure derived by comparative sequence analysis than that derived by the Zuker minimum free energy structure of the original (wild type or unmutated) RNA.     

A knowledge base of the chemical compounds of intermediary metabolism

This paper describes a publicly available knowledge base ofthe chemical compounds involved in intermediary metabolism.We consider the motivations for constructing a knowledge baseof metabolic compounds, the methodology by which it was constructed,and the information that it currently contains. Currently theknowledge base describes 981 compounds, listing for each: synonymsfor its name, a systematic name, CAS registry number, chemicalformula, molecular weight, chemical structure and two–dimensionaldisplay coordinates for the structure. The Compound KnowledgeBase (CompoundKB) illustrates several methodological principlesthat should guide the development of biological knowledge bases.I argue that biological datasets should be made available inmultiple representations to increase their accessibility toend users, and I present multiple representations of the CompoundKB(knowledge base, relational data base and ASN. 1 representations).I also analyze the general characteristics of these representationsto provide an understanding of their relative advantages anddisadvantages. Another principle is that the error rate of biologicaldata bases should be estimated and documented—this analysisis performed for the CompoundKB.     

Structure,properties, and biological relevance of the DNA and RNA G-quadruplexes: Overview 50 years after their discovery

G-quadruplexes (G4s), which are known to have important roles in regulation of key biological processes in both normal and pathological cells, are the most actively studied non-canonical structures of nucleic acids. In this review, we summarize the results of studies published in recent years that change significantly scientific views on various aspects of our understanding of quadruplexes. Modern notions on the polymorphism of DNA quadruplexes, on factors affecting thermodynamics and kinetics of G4 folding–unfolding, on structural organization of multiquadruplex systems, and on conformational features of RNA G4s and hybrid DNA–RNA G4s are discussed. Here we report the data on location of G4 sequence motifs in the genomes of eukaryotes, bacteria, and viruses, characterize G4-specific small-molecule ligands and proteins, as well as the mechanisms of their interactions with quadruplexes. New information on the structure and stability of G4s in telomeric DNA and oncogene promoters is discussed as well as proof being provided on the occurrence of G-quadruplexes in cells. Prominence is given to novel experimental techniques (single molecule manipulations, optical and magnetic tweezers, original chemical approaches, G4 detection in situ, in-cell NMR spectroscopy) that facilitate breakthroughs in the investigation of the structure and functions of G-quadruplexes.     

Sampled ensemble neutrality as a feature to classify potential structured RNAs

Background

Structured RNAs have many biological functions ranging from catalysis of chemical reactions to gene regulation. Yet, many homologous structured RNAs display most of their conservation at the secondary or tertiary structure level. As a result, strategies for structured RNA discovery rely heavily on identification of sequences sharing a common stable secondary structure. However, correctly distinguishing structured RNAs from surrounding genomic sequence remains challenging, especially during de novo discovery. RNA also has a long history as a computational model for evolution due to the direct link between genotype (sequence) and phenotype (structure). From these studies it is clear that evolved RNA structures, like protein structures, can be considered robust to point mutations. In this context, an RNA sequence is considered robust if its neutrality (extent to which single mutant neighbors maintain the same secondary structure) is greater than that expected for an artificial sequence with the same minimum free energy structure.

Results

In this work, we bring concepts from evolutionary biology to bear on the structured RNA de novo discovery process. We hypothesize that alignments corresponding to structured RNAs should consist of neutral sequences. We evaluate several measures of neutrality for their ability to distinguish between alignments of structured RNA sequences drawn from Rfam and various decoy alignments. We also introduce a new measure of RNA structural neutrality, the structure ensemble neutrality (SEN). SEN seeks to increase the biological relevance of existing neutrality measures in two ways. First, it uses information from an alignment of homologous sequences to identify a conserved biologically relevant structure for comparison. Second, it only counts base-pairs of the original structure that are absent in the comparison structure and does not penalize the formation of additional base-pairs.

Conclusion

We find that several measures of neutrality are effective at separating structured RNAs from decoy sequences, including both shuffled alignments and flanking genomic sequence. Furthermore, as an independent feature classifier to identify structured RNAs, SEN yields comparable performance to current approaches that consider a variety of features including stability and sequence identity. Finally, SEN outperforms other measures of neutrality at detecting mutational robustness in bacterial regulatory RNA structures.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-014-1203-8) contains supplementary material, which is available to authorized users.     

The contributions of dsRNA structure to Dicer specificity and efficiency

Dicer processes long double-stranded RNA (dsRNA) and pre-microRNAs to generate the functional intermediates (short interfering RNAs and microRNAs) of the RNA interference pathway. Here we identify features of RNA structure that affect Dicer specificity and efficiency. The data presented show that various attributes of the 3' end structure, including overhang length and sequence composition, play a primary role in determining the position of Dicer cleavage in both dsRNA and unimolecular, short hairpin RNA (shRNA). We also demonstrate that siRNA end structure affects overall silencing functionality. Awareness of these new features of Dicer cleavage specificity as it is related to siRNA functionality provides a more detailed understanding of the RNAi mechanism and can shape the development of hairpins with enhanced functionality.     

Global RNA Structure Analysis of Poliovirus Identifies a Conserved RNA Structure Involved in Viral Replication and Infectivity

The genomes of RNA viruses often contain RNA structures that are crucial for translation and RNA replication and may play additional, uncharacterized roles during the viral replication cycle. For the picornavirus family member poliovirus, a number of functional RNA structures have been identified, but much of its genome, especially the open reading frame, has remained uncharacterized. We have now generated a global RNA structure map of the poliovirus genome using a chemical probing approach that interrogates RNA structure with single-nucleotide resolution. In combination with orthogonal evolutionary analyses, we uncover several conserved RNA structures in the open reading frame of the viral genome. To validate the ability of our global analyses to identify functionally important RNA structures, we further characterized one of the newly identified structures, located in the region encoding the RNA-dependent RNA polymerase, 3D^pol, by site-directed mutagenesis. Our results reveal that the structure is required for viral replication and infectivity, since synonymous mutants are defective in these processes. Furthermore, these defects can be partially suppressed by mutations in the viral protein 3C^pro, which suggests the existence of a novel functional interaction between an RNA structure in the 3D^pol-coding region and the viral protein(s) 3C^pro and/or its precursor 3CD^pro.     

Base composition studies on mitochondrial 4 S RNA from rat liver and Morris hepatomas 5123D and 7777.

The major and modified base composition of mitochondrial 4 S RNA from rat liver and from Morris hepatomas 5123D and 7777 has been determined for 16 constituents using a chemical tritium-derivative method. The base composition of these mitochondrial 4 S RNA preparations was compared with the base composition of cytoplasmic and bacterial (Escherichia coli B and Bacillus subtilis) 4-S RNAs. The results of these studies are: 1. When compared with cytoplasmic 4 S RNA, the liver and hepatoma mitochondrial 4-S RNAs are characterized by high (A + U)/(G + C) ratios and low overall degrees of base methylation and modification. 2. The mammalian mitochondrial 4-S RNAs are qualitatively even more different from the bacterial 4-S RNAs than from their cytoplasmic counterparts. Thus, several modified constituents found in both cytoplasmic and mitochondrial 4 S RNA are absent from the bacterial 4-S RNAs. 3. Mitochondrial 4S RNA from both hepatomas was found to be under-methylated and undermodified when compared with normal liver mitochondrial 4S RNA. This trend is more pronounced for the rapidly growing hepatoma 7777 (i.e., 17% undermethylation) than for the more slowly growing hepatoma 5123D (i.e., 8% undermethylation). These findings are discussed in relationship to (1) results of other authors on composition of mitochondrial 4 S RNA, (2) special features of structure and biosynthesis of mitochondrial 4 S RNA, (3) the possible evolutionary origin of mitochondria and (4) the possible role played by aberrant mitochondrial 4 S RNA in altered mitochondrial protein synthesis in tumors.     

基于注意力机制的RNA碱基关联图预测方法

目的 长链非编码RNA在遗传、代谢和基因表达调控等方面发挥着重要作用。然而，传统的实验方法解析RNA的三级结构耗时长、费用高且操作要求高。此外，通过计算方法来预测RNA的三级结构在近十年来无突破性进展。因此，需要提出新的预测算法来准确的预测RNA的三级结构。所以，本文发展可以用于提高RNA三级结构预测准确性的碱基关联图预测方法。方法 为了利用RNA理化特征信息，本文应用多层全卷积神经网络和循环神经网络的深度学习算法来预测RNA碱基间的接触概率，并通过注意力机制处理RNA序列中碱基间相互依赖的特征。结果 通过多层神经网络与注意力机制结合，本文方法能够有效得到RNA特征值中局部和全局的信息，提高了模型的鲁棒性和泛化能力。检验计算表明，所提出模型对序列长度L的4种标准（L/10、L/5、L/2、L）碱基关联图的预测准确率分别达到0.84、0.82、0.82和0.75。结论 基于注意力机制的深度学习预测算法能够提高RNA碱基关联图预测的准确率，从而帮助RNA三级结构的预测。