Efficient pairwise RNA structure prediction and alignment using sequence alignment constraints

BACKGROUND: We are interested in the problem of predicting secondary structure for small sets of homologous RNAs, by incorporating limited comparative sequence information into an RNA folding model. The Sankoff algorithm for simultaneous RNA folding and alignment is a basis for approaches to this problem. There are two open problems in applying a Sankoff algorithm: development of a good unified scoring system for alignment and folding and development of practical heuristics for dealing with the computational complexity of the algorithm. RESULTS: We use probabilistic models (pair stochastic context-free grammars, pairSCFGs) as a unifying framework for scoring pairwise alignment and folding. A constrained version of the pairSCFG structural alignment algorithm was developed which assumes knowledge of a few confidently aligned positions (pins). These pins are selected based on the posterior probabilities of a probabilistic pairwise sequence alignment. CONCLUSION: Pairwise RNA structural alignment improves on structure prediction accuracy relative to single sequence folding. Constraining on alignment is a straightforward method of reducing the runtime and memory requirements of the algorithm. Five practical implementations of the pairwise Sankoff algorithm – this work (Consan), David Mathews' Dynalign, Ian Holmes' Stemloc, Ivo Hofacker's PMcomp, and Jan Gorodkin's FOLDALIGN – have comparable overall performance with different strengths and weaknesses

Can Clustal-style progressive pairwise alignment of multiple sequences be used in RNA secondary structure prediction?

<p>Abstract</p> <p>Background</p> <p>In ribonucleic acid (RNA) molecules whose function depends on their final, folded three-dimensional shape (such as those in ribosomes or spliceosome complexes), the secondary structure, defined by the set of internal basepair interactions, is more consistently conserved than the primary structure, defined by the sequence of nucleotides.</p> <p>Results</p> <p>The research presented here investigates the possibility of applying a progressive, pairwise approach to the alignment of multiple RNA sequences by simultaneously predicting an energy-optimized consensus secondary structure. We take an existing algorithm for finding the secondary structure common to two RNA sequences, Dynalign, and alter it to align profiles of multiple sequences. We then explore the relative successes of different approaches to designing the tree that will guide progressive alignments of sequence profiles to create a multiple alignment and prediction of conserved structure.</p> <p>Conclusion</p> <p>We have found that applying a progressive, pairwise approach to the alignment of multiple ribonucleic acid sequences produces highly reliable predictions of conserved basepairs, and we have shown how these predictions can be used as constraints to improve the results of a single-sequence structure prediction algorithm. However, we have also discovered that the amount of detail included in a consensus structure prediction is highly dependent on the order in which sequences are added to the alignment (the guide tree), and that if a consensus structure does not have sufficient detail, it is less likely to provide useful constraints for the single-sequence method.</p

RNA secondary structure prediction from multi-aligned sequences

It has been well accepted that the RNA secondary structures of most functional non-coding RNAs (ncRNAs) are closely related to their functions and are conserved during evolution. Hence, prediction of conserved secondary structures from evolutionarily related sequences is one important task in RNA bioinformatics; the methods are useful not only to further functional analyses of ncRNAs but also to improve the accuracy of secondary structure predictions and to find novel functional RNAs from the genome. In this review, I focus on common secondary structure prediction from a given aligned RNA sequence, in which one secondary structure whose length is equal to that of the input alignment is predicted. I systematically review and classify existing tools and algorithms for the problem, by utilizing the information employed in the tools and by adopting a unified viewpoint based on maximum expected gain (MEG) estimators. I believe that this classification will allow a deeper understanding of each tool and provide users with useful information for selecting tools for common secondary structure predictions.Comment: A preprint of an invited review manuscript that will be published in a chapter of the book `Methods in Molecular Biology'. Note that this version of the manuscript may differ from the published versio

RNAstructure: software for RNA secondary structure prediction and analysis

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Results</p> <p>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.</p> <p>Conclusion</p> <p>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 <url>http://rna.urmc.rochester.edu/RNAstructure.html</url>.</p

Clinical Manifestations and Case Management of Ebola Haemorrhagic Fever caused by a newly identified virus strain, Bundibugyo, Uganda, 2007-2008

A confirmed Ebola haemorrhagic fever (EHF) outbreak in Bundibugyo, Uganda, November 2007-February 2008, was caused by a putative new species (Bundibugyo ebolavirus). It included 93 putative cases, 56 laboratory-confirmed cases, and 37 deaths (CFR = 25%). Study objectives are to describe clinical manifestations and case management for 26 hospitalised laboratory-confirmed EHF patients. Clinical findings are congruous with previously reported EHF infections. The most frequently experienced symptoms were non-bloody diarrhoea (81%), severe headache (81%), and asthenia (77%). Seven patients reported or were observed with haemorrhagic symptoms, six of whom died. Ebola care remains difficult due to the resource-poor setting of outbreaks and the infection-control procedures required. However, quality data collection is essential to evaluate case definitions and therapeutic interventions, and needs improvement in future epidemics. Organizations usually involved in EHF case management have a particular responsibility in this respect

PicXAA-R: Efficient structural alignment of multiple RNA sequences using a greedy approach

<p>Abstract</p> <p>Background</p> <p>Accurate and efficient structural alignment of non-coding RNAs (ncRNAs) has grasped more and more attentions as recent studies unveiled the significance of ncRNAs in living organisms. While the Sankoff style structural alignment algorithms cannot efficiently serve for multiple sequences, mostly progressive schemes are used to reduce the complexity. However, this idea tends to propagate the early stage errors throughout the entire process, thereby degrading the quality of the final alignment. For multiple protein sequence alignment, we have recently proposed PicXAA which constructs an accurate alignment in a non-progressive fashion.</p> <p>Results</p> <p>Here, we propose PicXAA-R as an extension to PicXAA for greedy structural alignment of ncRNAs. PicXAA-R efficiently grasps both folding information within each sequence and local similarities between sequences. It uses a set of probabilistic consistency transformations to improve the posterior base-pairing and base alignment probabilities using the information of all sequences in the alignment. Using a graph-based scheme, we greedily build up the structural alignment from sequence regions with high base-pairing and base alignment probabilities.</p> <p>Conclusions</p> <p>Several experiments on datasets with different characteristics confirm that PicXAA-R is one of the fastest algorithms for structural alignment of multiple RNAs and it consistently yields accurate alignment results, especially for datasets with locally similar sequences. PicXAA-R source code is freely available at: <url>http://www.ece.tamu.edu/~bjyoon/picxaa/</url>.</p

Ebola and Marburg Hemorrhagic Fevers: Neglected Tropical Diseases?

Ebola hemorrhagic fever (EHF) and Marburg hemorrhagic fever (MHF) are rare viral diseases, endemic to central Africa. The overall burden of EHF and MHF is small in comparison to the more common protozoan, helminth, and bacterial diseases typically referred to as neglected tropical diseases (NTDs). However, EHF and MHF outbreaks typically occur in resource-limited settings, and many aspects of these outbreaks are a direct consequence of impoverished conditions. We will discuss aspects of EHF and MHF disease, in comparison to the “classic” NTDs, and examine potential ways forward in the prevention and control of EHF and MHF in sub-Saharan Africa, as well as examine the potential for application of novel vaccines or antiviral drugs for prevention or control of EHF and MHF among populations at highest risk for disease

Polyfunctional Hiv-Specific Antibody Responses Are Associated with Spontaneous Hiv Control

Elite controllers (ECs) represent a unique model of a functional cure for HIV-1 infection as these individuals develop HIV-specific immunity able to persistently suppress viremia. Because accumulating evidence suggests that HIV controllers generate antibodies with enhanced capacity to drive antibody-dependent cellular cytotoxicity (ADCC) that may contribute to viral containment, we profiled an array of extra-neutralizing antibody effector functions across HIV-infected populations with varying degrees of viral control to define the characteristics of antibodies associated with spontaneous control. While neither the overall magnitude of antibody titer nor individual effector functions were increased in ECs, a more functionally coordinated innate immune–recruiting response was observed. Specifically, ECs demonstrated polyfunctional humoral immune responses able to coordinately recruit ADCC, other NK functions, monocyte and neutrophil phagocytosis, and complement. This functionally coordinated response was associated with qualitatively superior IgG3/IgG1 responses, whereas HIV-specific IgG2/IgG4 responses, prevalent among viremic subjects, were associated with poorer overall antibody activity. Rather than linking viral control to any single activity, this study highlights the critical nature of functionally coordinated antibodies in HIV control and associates this polyfunctionality with preferential induction of potent antibody subclasses, supporting coordinated antibody activity as a goal in strategies directed at an HIV-1 functional cure

Recovery after spinal cord relapse in multiple sclerosis is predicted by radial diffusivity

Background: The aim of this study was to determine whether the diffusion tensor-derived radial diffusivity and axial diffusivity, measured in the cortico-spinal tract in the cervical cord, predict clinical recovery after a cord relapse in patients with multiple sclerosis, and change over time