RNA challenges for computational chemists

Some experimental results for the thermodynamics of RNA folding cannot be explained by simple pairwise hydrogen-bonding models. Such effects include the stabilities of isoguanosine-isocytidine (iG-iC) base pairs and of various 2 x 2 nucleotide internal loops. Presumably, these results can be explained by base stacking effects, which can be partitioned into Coulombic and overlap effects. We review experimental measurements that provide benchmarks for testing the approximations and theories used for modeling nucleic acids. Quantitative agreement between experiment and theory will indicate understanding of the interactions determining RNA stability and structure.     

On comparing two structured RNA multiple alignments

We present a method, called BlockMatch, for aligning two blocks, where a block is an RNA multiple sequence alignment with the consensus secondary structure of the alignment in Stockholm format. The method employs a quadratic-time dynamic programming algorithm for aligning columns and column pairs of the multiple alignments in the blocks. Unlike many other tools that can perform pairwise alignment of either single sequences or structures only, BlockMatch takes into account the characteristics of all the sequences in the blocks along with their consensus structures during the alignment process, thus being able to achieve a high-quality alignment result. We apply BlockMatch to phylogeny reconstruction on a set of 5S rRNA sequences taken from fifteen bacteria species. Experimental results showed that the phylogenetic tree generated by our method is more accurate than the tree constructed based on the widely used ClustalW tool. The BlockMatch algorithm is implemented into a web server, accessible at http://bioinformatics.njit.edu/blockmatch. A jar file of the program is also available for download from the web server.     

Evolution of RNA virus in spatially structured heterogeneous environments

A hallmark of the infectious cycle for many RNA viruses parasitizing multicellular hosts is the need to invade and successfully replicate in tissues that comprise a variety of cell types. Thus, multicellular hosts represent a heterogeneous environment to evolving viral populations. To understand viral adaptation to multicellular hosts, we took a double approach. First, we developed a mathematical model that served to make predictions concerning the dynamics of viral populations evolving in heterogeneous environments. Second, the predictions were tested by evolving vesicular stomatitis virus in vitro on a spatially structured environment formed by three different cell types. In the absence of gene flow, adaptation was tissue-specific, but fitness in all tissues decreased with migration rate. The performance in a given tissue was negatively correlated with its distance to the tissue hosting the population. This correlation decreased with migration rate.     

A structured retroviral RNA element that mediates nucleocytoplasmic export of intron-containing RNA.

A common feature of gene expression in all retroviruses is that unspliced, intron-containing RNA is exported to the cytoplasm despite the fact that cellular RNAs which contain introns are usually restricted to the nucleus. In complex retroviruses, the export of intron-containing RNA is mediated by specific viral regulatory proteins (e.g., human immunodeficiency virus type 1 [HIV-1] Rev) that bind to elements in the viral RNA. However, simpler retroviruses do not encode such regulatory proteins. Here we show that the genome of the simpler retrovirus Mason-Pfizer monkey virus (MPMV) contains an element that serves as an autonomous nuclear export signal for intron-containing RNA. This element is essential for MPMV replication; however, its function can be complemented by HIV-1 Rev and the Rev-responsive element. The element can also facilitate the export of cellular intron-containing RNA. These results suggest that the MPMV element mimics cellular RNA transport signals and mediates RNA export through interaction with endogenous cellular factors.     

RSEARCH: Finding homologs of single structured RNA sequences

Background

For many RNA molecules, secondary structure rather than primary sequence is the evolutionarily conserved feature. No programs have yet been published that allow searching a sequence database for homologs of a single RNA molecule on the basis of secondary structure.

Results

We have developed a program, RSEARCH, that takes a single RNA sequence with its secondary structure and utilizes a local alignment algorithm to search a database for homologous RNAs. For this purpose, we have developed a series of base pair and single nucleotide substitution matrices for RNA sequences called RIBOSUM matrices. RSEARCH reports the statistical confidence for each hit as well as the structural alignment of the hit. We show several examples in which RSEARCH outperforms the primary sequence search programs BLAST and SSEARCH. The primary drawback of the program is that it is slow. The C code for RSEARCH is freely available from our lab's website.

Conclusion

RSEARCH outperforms primary sequence programs in finding homologs of structured RNA sequences.

     

Novel role for RNase PH in the degradation of structured RNA

Escherichia coli contains multiple 3' to 5' RNases, of which two, RNase PH and polynucleotide phosphorylase (PNPase), use inorganic phosphate as a nucleophile to catalyze RNA cleavage. It is known that an absence of these two enzymes causes growth defects, but the basis for these defects has remained undefined. To further an understanding of the function of these enzymes, the degradation pattern of different cellular RNAs was analyzed. It was observed that an absence of both enzymes results in the appearance of novel mRNA degradation fragments. Such fragments were also observed in strains containing mutations in RNase R and PNPase, enzymes whose collective absence is known to cause an accumulation of structured RNA fragments. Additional experiments indicated that the growth defects of strains containing RNase R and PNPase mutations were exacerbated upon RNase PH removal. Taken together, these observations suggested that RNase PH could play a role in structured RNA degradation. Biochemical experiments with RNase PH demonstrated that this enzyme digests through RNA duplexes of moderate stability. In addition, mapping and sequence analysis of an mRNA degradation fragment that accumulates in the absence of the phosphorolytic enzymes revealed the presence of an extended stem-loop motif at the 3' end. Overall, these results indicate that RNase PH plays a novel role in the degradation of structured RNAs and provides a potential explanation for the growth defects caused by an absence of the phosphorolytic RNases.     

Informatics in neuroscience

The application of informatics to neuroscience goes far beyond 'traditional' bioinformatics modalities such as DNA sequences. In this review, we describe how informatics is being used to study the nervous system at multiple levels, spanning scales from molecules to behavior. The continuing development of standards for data exchange and interoperability, together with increasing awareness and acceptance of the importance of data sharing, are among the key efforts required to advance the field.     

The paradoxical behavior of a highly structured misfolded intermediate in RNA folding

Like many structured RNAs, the Tetrahymena group I ribozyme is prone to misfolding. Here we probe a long-lived misfolded species, referred to as M, and uncover paradoxical aspects of its structure and folding. Previous work indicated that a non-native local secondary structure, termed alt P3, led to formation of M during folding in vitro. Surprisingly, hydroxyl radical footprinting, fluorescence measurements with site-specifically incorporated 2-aminopurine, and functional assays indicate that the native P3, not alt P3, is present in the M state. The paradoxical behavior of alt P3 presumably arises because alt P3 biases folding toward M, but, after commitment to this folding pathway and before formation of M, alt P3 is replaced by P3. Further, structural and functional probes demonstrate that the misfolded ribozyme contains extensive native structure, with only local differences between the two states, and the misfolded structure even possesses partial catalytic activity. Despite the similarity of these structures, re-folding of M to the native state is very slow and is strongly accelerated by urea, Na+, and increased temperature and strongly impeded by Mg2+ and the presence of native peripheral contacts. The paradoxical observations of extensive native structure within the misfolded species but slow conversion of this species to the native state are readily reconciled by a model in which the misfolded state is a topological isomer of the native state, and computational results support the feasibility of this model. We speculate that the complex topology of RNA secondary structures and the inherent rigidity of RNA helices render kinetic traps due to topological isomers considerably more common for RNA than for proteins.     

Evidence for ligandable sites in structured RNA throughout the Protein Data Bank

RNA has attracted considerable attention as a target for small molecules. However, methods to identify, study, and characterize suitable RNA targets have lagged behind strategies for protein targets. One approach that has received considerable attention for protein targets has been to utilize computational analysis to investigate ligandable “pockets” on proteins that are amenable to small molecule binding. These studies have shown that selected physical properties of pockets are important parameters that govern the ability of a structure to bind to small molecules. This work describes a similar analysis to study pockets on all RNAs in the Protein Data Bank (PDB). Using parameters such as buriedness, hydrophobicity, volume, and other properties, the set of all RNAs is analyzed and compared to all proteins. Considerable overlap is observed between the properties of pockets on RNAs and proteins. Thus, many RNAs are capable of populating conformations with pockets that are likely suitable for small molecule binding. Further, principal moment of inertia (PMI) calculations reveal that liganded RNAs exist in diverse structural space, much of which overlaps with protein structural space. Taken together, these results suggest that complex folded RNAs adopt unique structures with pockets that may represent viable opportunities for small molecule targeting.     

The highly structured encapsidation signal of MuLV RNA is involved in the nuclear export of its unspliced RNA

p27Kip1 (p27) influences cell division by regulating nuclear cyclin-dependent kinases. Before binding, p27 is at least partially disordered and folds upon binding its Cdk/cyclin targets. 30-40% of human proteins, including p27, are predicted to contain disordered segments, and have been termed intrinsically unstructured proteins (IUPs). Unfortunately, the inherent dynamics of IUPs hamper detailed analysis of their structure/function relationships. Here, we describe the use of molecular dynamics (MD) computations and solution NMR spectroscopy to reveal that several segments of the p27 kinase inhibitory domain (p27-KID), in addition to the previously characterized helical segment, exist as highly populated, intrinsically folded structural units (IFSUs). Several IFSUs resemble structural features of bound p27-KID, while another exhibits alternative conformations. Interestingly, the highly conserved, specificity determining segment of p27 is shown to be highly disordered. Elucidation of IFSUs within p27-KID allows consideration of their influences on the thermodynamics and kinetics of Cdk/cyclin binding. The degree to which IFSUs are populated within p27-KID is surprising and suggests that other putative IUPs contain IFSUs that may be studied using similar techniques.     

生物多样性信息学研究进展

生物多样性信息学是一门蓬勃发展的新学科。它将现代的信息技术带入生物多样性及其相关学科的研究领域。它在生物多样性基础数据的数字化、模型工具和各种工具软件的开发、数据整合, 以及全球、地区和国家尺度生物多样性信息网络等多个方面的发展, 向我们展示了未来在全球范围内自由、免费共享生物多样性数据和信息, 以及人们行动起来共同关注、调查与监测野外生物多样性的前景。目前, 已有大量数字化的物种编目、标本馆标本、多媒体影像、研究文献等生物多样性基础信息可以通过互联网检索和利用。其中, 最值得关注的是一些成功的国际性研究项目, 如物种2000、全球生物多样性信息网络、生命条形码以及网络生命大百科全书。这些项目的成功不仅体现在对大量基础信息和数据的发布, 而且它们通过与生物多样性信息标准TDWG(Biodiversity Information Standards: TDWG)的合作, 推动了达尔文核心标准(Darwin Core)等一些重要的生物多样性信息标准的应用, 以及地区和国家性生物多样性信息节点的建立, 这些都为将来全球范围生物多样性信息的共享和数据交换奠定了重要基础。在数字化信息的基础上, 研究人员也开发了一些在特定研究领域应用的数据挖掘和模型工具, 例如基于数字化标本的地理分布预测工具MAXENT, 分类学专家知识管理的LifeDesk。公民科学理念的发展则向我们展示了公众和科学爱好者广泛参与以互联网为基础的生物多样性信息学研究活动。因此, 生物多样性信息学的发展前景广阔, 它将为我们实现全球保护战略目标, 应对生物多样性危机, 解决全球气候变化条件下生物多样性资源管理和利用建立坚实的信息基础。     

Programmed ribosomal frameshifting in SIV is induced by a highly structured RNA stem-loop

Simian immunodeficiency virus (SIV), like its human homologues (HIV-1, HIV-2), requires a -1 translational frameshift event to properly synthesize all of the proteins required for viral replication. The frameshift mechanism is dependent upon a seven-nucleotide slippery sequence and a downstream RNA structure. In SIV, the downstream RNA structure has been proposed to be either a stem-loop or a pseudoknot. Here, we report the functional, structural and thermodynamic characterization of the SIV frameshift site RNA. Translational frameshift assays indicate that a stem-loop structure is sufficient to promote efficient frameshifting in vitro. NMR and thermodynamic studies of SIV RNA constructs of varying length further support the absence of any pseudoknot interaction and indicate the presence of a stable stem-loop structure. We determined the structure of the SIV frameshift-inducing RNA by NMR. The structure reveals a highly ordered 12 nucleotide loop containing a sheared G-A pair, cross-strand adenine stacking, two G-C base-pairs, and a novel CCC triloop turn. The loop structure and its high thermostability preclude pseudoknot formation. Sequence conservation and modeling studies suggest that HIV-2 RNA forms the same structure. We conclude that, like the main sub-groups of HIV-1, SIV and HIV-2 utilize stable stem-loop structures to function as a thermodynamic barrier to translation, thereby inducing ribosomal pausing and frameshifting.