RNA-RNA interaction prediction based on multiple sequence alignments

Many computerized methods for RNA-RNA interaction structure prediction have been developed. Recently, $O(N^6)$ time and $O(N^4)$ space dynamic programming algorithms have become available that compute the partition function of RNA-RNA interaction complexes. However, few of these methods incorporate the knowledge concerning related sequences, thus relevant evolutionary information is often neglected from the structure determination. Therefore, it is of considerable practical interest to introduce a method taking into consideration both thermodynamic stability and sequence covariation. We present the \emph{a priori} folding algorithm \texttt{ripalign}, whose input consists of two (given) multiple sequence alignments (MSA). \texttt{ripalign} outputs (1) the partition function, (2) base-pairing probabilities, (3) hybrid probabilities and (4) a set of Boltzmann-sampled suboptimal structures consisting of canonical joint structures that are compatible to the alignments. Compared to the single sequence-pair folding algorithm \texttt{rip}, \texttt{ripalign} requires negligible additional memory resource. Furthermore, we incorporate possible structure constraints as input parameters into our algorithm. The algorithm described here is implemented in C as part of the \texttt{rip} package. The supplemental material, source code and input/output files can freely be downloaded from \url{http://www.combinatorics.cn/cbpc/ripalign.html}. \section{Contact} Christian Reidys \texttt{[email protected]}Comment: 8 pages, 9 figure

Comprehensive sequence-to-function mapping of cofactor-dependent RNA catalysis in the glmS ribozyme.

Massively parallel, quantitative measurements of biomolecular activity across sequence space can greatly expand our understanding of RNA sequence-function relationships. We report the development of an RNA-array assay to perform such measurements and its application to a model RNA: the core glmS ribozyme riboswitch, which performs a ligand-dependent self-cleavage reaction. We measure the cleavage rates for all possible single and double mutants of this ribozyme across a series of ligand concentrations, determining kcat and KM values for active variants. These systematic measurements suggest that evolutionary conservation in the consensus sequence is driven by maintenance of the cleavage rate. Analysis of double-mutant rates and associated mutational interactions produces a structural and functional mapping of the ribozyme sequence, revealing the catalytic consequences of specific tertiary interactions, and allowing us to infer structural rearrangements that permit certain sequence variants to maintain activity

Profile Context-Sensitive HMMs for Probabilistic Modeling of Sequences With Complex Correlations

The profile hidden Markov model is a specific type of HMM that is well suited for describing the common features of a set of related sequences. It has been extensively used in computational biology, where it is still one of the most popular tools. In this paper, we propose a new model called the profile context-sensitive HMM. Unlike traditional profile-HMMs, the proposed model is capable of describing complex long-range correlations between distant symbols in a consensus sequence. We also introduce a general algorithm that can be used for finding the optimal state-sequence of an observed symbol sequence based on the given profile-csHMM. The proposed model has an important application in RNA sequence analysis, especially in modeling and analyzing RNA pseudoknots

Reconstructing phylogeny from RNA secondary structure via simulated evolution

DNA sequences of genes encoding functional RNA molecules (e.g., ribosomal RNAs) are commonly used in phylogenetics (i.e. to infer evolutionary history). Trees derived from ribosomal RNA (rRNA) sequences, however, are inconsistent with other molecular data in investigations of deep branches in the tree of life. Since much of te functional constraints on the gene products (i.e. RNA molecules) relate to three-dimensional structure, rather than their actual sequences, accumulated mutations in the gene sequences may obscure phylogenetic signal over very large evolutionary time-scales. Variation in structure, however, may be suitable for phylogenetic inference even under extreme sequence divergence. To evaluate qualitatively the manner in which structural evolution relates to sequence change, we simulated the evolution of RNA sequences under various constraints on structural change

DNA editing in DNA/RNA hybrids by adenosine deaminases that act on RNA.

Adenosine deaminases that act on RNA (ADARs) carry out adenosine (A) to inosine (I) editing reactions with a known requirement for duplex RNA. Here, we show that ADARs also react with DNA/RNA hybrid duplexes. Hybrid substrates are deaminated efficiently by ADAR deaminase domains at dA-C mismatches and with E to Q mutations in the base flipping loop of the enzyme. For a long, perfectly matched hybrid, deamination is more efficient with full length ADAR2 than its isolated deaminase domain. Guide RNA strands for directed DNA editing by ADAR were used to target six different 2΄-deoxyadenosines in the M13 bacteriophage ssDNA genome. DNA editing efficiencies varied depending on the sequence context of the editing site consistent with known sequence preferences for ADARs. These observations suggest the reaction within DNA/RNA hybrids may be a natural function of human ADARs. In addition, this work sets the stage for development of a new class of genome editing tools based on directed deamination of 2΄-deoxyadenosines in DNA/RNA hybrids

Dimerisation of HIV-2 genomic RNA is linked to efficient RNA packaging, normal particle maturation and viral infectivity.

BACKGROUND: Retroviruses selectively encapsidate two copies of their genomic RNA, the Gag protein binding a specific RNA motif in the 5' UTR of the genome. In human immunodeficiency virus type 2 (HIV-2), the principal packaging signal (Psi) is upstream of the major splice donor and hence is present on all the viral RNA species. Cotranslational capture of the full length genome ensures specificity. HIV-2 RNA dimerisation is thought to occur at the dimer initiation site (DIS) located in stem-loop 1 (SL-1), downstream of the main packaging determinant. However, the HIV-2 packaging signal also contains a palindromic sequence (pal) involved in dimerisation. In this study, we analysed the role of the HIV-2 packaging signal in genomic RNA dimerisation in vivo and its implication in viral replication. RESULTS: Using a series of deletion and substitution mutants in SL-1 and the Psi region, we show that in fully infectious HIV-2, genomic RNA dimerisation is mediated by the palindrome pal. Mutation of the DIS had no effect on dimerisation or viral infectivity, while mutations in the packaging signal severely reduce both processes as well as RNA encapsidation. Electron micrographs of the Psi-deleted virions revealed a significant reduction in the proportion of mature particles and an increase in that of particles containing multiple cores. CONCLUSION: In addition to its role in RNA encapsidation, the HIV-2 packaging signal contains a palindromic sequence that is critical for genomic RNA dimerisation. Encapsidation of a dimeric genome seems required for the production of infectious mature particles, and provides a promising therapeutic target