A base-pairing model of duplex formation. I. Watson-Crick pairing geometries

We present a base-pairing model of oligonucleotide duplex formation and show in detail its equivalence to the nearest-neighbor dimer methods from fits to free energy of duplex formation data for short DNA-DNA and DNA-RNA hybrids containing only Watson-Crick pairs. For completeness, the corresponding RNA-RNA parameters are included. In this approach, the connection between rank-deficient polymer and rank-determinant oligonucleotide parameter sets for DNA duplexes is transparent. The method is generalized to include RNA-DNA hybrids where the rank-deficient model with 11 dimer parameters in fact provides slightly improved predictions relative to the standard method with 16 independent dimer parameters (DeltaG mean errors of 4.5 and 5.4%, respectively).     

Revolutions in RNA secondary structure prediction

RNA structure formation is hierarchical and, therefore, secondary structure, the sum of canonical base-pairs, can generally be predicted without knowledge of the three-dimensional structure. Secondary structure prediction algorithms evolved from predicting a single, lowest free energy structure to their current state where statistics can be determined from the thermodynamic ensemble. This article reviews the free energy minimization technique and the salient revolutions in the dynamic programming algorithm methods for secondary structure prediction. Emphasis is placed on highlighting the recently developed method, which statistically samples structures from the complete Boltzmann ensemble.     

RNA二级结构数据库

RNA分子众多、结构复杂、功能重要,已经成为当前重要的研究热点之一。RNA的功能与结构密切相关,伴随RNA分子及功能的发现,建立了有关RNA二级结构的数据库,一方面有助于理解RNA功能的结构基础,一方面有助于开发各种有关RNA结构的预测模型。本文对近年常见的RNA二级结构数据库作一概述,希望有助于相关工作者更好地了解与应用相关数据。     

RNA二级子结构的计数

多核苷酸的二级结构可视为一类顶点标号平面图,通常通过枚举每类RNA二级结构图的各种子图来计算其递推公式．本文给出了若干限制端环长度的RNA二级子结构的递推公式及渐近值。     

Measuring the thermodynamics of RNA secondary structure formation

The thermodynamics of RNA secondary structure formation in small model systems provides a database for predicting RNA structure from sequence. Methods for making these measurements are reviewed with emphasis on optical methods and treatment of experimental errors. Analysis of experimental results in terms of simple nearest-neighbor models is presented. Some measured sequence dependences of non-Watson-Crick motifs are discussed. © 1998 John Wiley & Sons, Inc. Biopoly 44: 309–319, 1997     

基于动态权重匹配的RNA折叠算法

本文在最大权重匹配（Maximum Weighted Matching,MWM）算法的基顾础上引入与茎区长度相关的动态权重,采用一种递归算法逐步寻找具有最大权重和的茎区。从而最终确定RNA的二级结构．该算法避开了繁杂的自由能计算,同样也能达到较高的预测精确度并且还能预测到大多数类型的潜在假结（pseudoknots）．     

Semiautomated improvement of RNA alignments

We have developed a semiautomated RNA sequence editor (SARSE) that integrates tools for analyzing RNA alignments. The editor highlights different properties of the alignment by color, and its integrated analysis tools prevent the introduction of errors when doing alignment editing. SARSE readily connects to external tools to provide a flexible semiautomatic editing environment. A new method, Pcluster, is introduced for dividing the sequences of an RNA alignment into subgroups with secondary structure differences. Pcluster was used to evaluate 574 seed alignments obtained from the Rfam database and we identified 71 alignments with significant prediction of inconsistent base pairs and 102 alignments with significant prediction of novel base pairs. Four RNA families were used to illustrate how SARSE can be used to manually or automatically correct the inconsistent base pairs detected by Pcluster: the mir-399 RNA, vertebrate telomase RNA (vert-TR), bacterial transfer-messenger RNA (tmRNA), and the signal recognition particle (SRP) RNA. The general use of the method is illustrated by the ability to accommodate pseudoknots and handle even large and divergent RNA families. The open architecture of the SARSE editor makes it a flexible tool to improve all RNA alignments with relatively little human intervention. Online documentation and software are available at (http://sarse.ku.dk).     

Incorporating global features of RNA motifs in predictions for an ensemble of secondary structures for encapsidated MS2 bacteriophage RNA

The secondary structure of encapsidated MS2 genomic RNA poses an interesting RNA folding challenge. Cryoelectron microscopy has demonstrated that encapsidated MS2 RNA is well-ordered. Models of MS2 assembly suggest that the RNA hairpin-protein interactions and the appropriate placement of hairpins in the MS2 RNA secondary structure can guide the formation of the correct icosahedral particle. The RNA hairpin motif that is recognized by the MS2 capsid protein dimers, however, is energetically unfavorable, and thus free energy predictions are biased against this motif. Computer programs called Crumple, Sliding Windows, and Assembly provide useful tools for prediction of viral RNA secondary structures when the traditional assumptions of RNA structure prediction by free energy minimization may not apply. These methods allow incorporation of global features of the RNA fold and motifs that are difficult to include directly in minimum free energy predictions. For example, with MS2 RNA the experimental data from SELEX experiments, crystallography, and theoretical calculations of the path for the series of hairpins can be incorporated in the RNA structure prediction, and thus the influence of free energy considerations can be modulated. This approach thoroughly explores conformational space and generates an ensemble of secondary structures. The predictions from this new approach can test hypotheses and models of viral assembly and guide construction of complete three-dimensional models of virus particles.     

正在崛起的RNA纳米技术

RNA纳米技术得益于纽约大学西曼(Nadrian C.Seeman)教授开创的DNA纳米技术,RNA是由腺嘌呤(A)、尿嘧啶(U)、鸟嘌呤(G)和胞嘧啶(C)构成的一种核糖核酸高分子,与DNA的Watson-Crick碱基配对(A-T,G-C)的双螺旋链的结构不完全一样,RNA的二级结构里经常出现一些非传统的碱基配对如环环相互作用,这些非传统配对促使RNA分子折叠成刚性结构。本文综述了正在崛起的RNA纳米技术,列举了一些著名的实验,如郭培宣(Peixuan Guo)等从自然界的phi29噬菌体中发现的pRNA纳米马达是由六个小RNA分子构成的六环结构,Jaeger等发展了RNA构造术(RNA-tectonics),根据已知的RNA分子的碱基和非传统配对,他们设计利用小RNA分子构造二聚体、一维线性多聚体、和二维网状的七巧板迷宫(jigsaw puzzle)等图案,用tRNA分子或设计用几条RNA分子来构建多面体如立方体和八面体等立体结构等。RNA纳米技术正在崛起,它将在医学、生物技术、合成生物学和纳米技术领域扮演重要的角色。     

Statistics of RNA melting kinetics

We present and study the behavior of a simple kinetic model for the melting of RNA secondary structures, given that those structures are known. The model is then used as a map that. assigns structure dependent overall rate constants of melting (or refolding) to a sequence. This induces a landscape of reaction rates, or activation energies, over the space of sequences with fixed length. We study the distribution and the correlation structure of these activation energies. Correspondence to: P. Schuster     

RNA二级结构预测系统构建

运用下列RNA二级结构预测算法:碱基最大配对方法、Zuker极小化自由能方法、螺旋区最优堆积、螺旋区随机堆积和所有可能组合方法与基于一级螺旋区的RNA二级结构绘图技术, 构建了RNA二级结构预测系统Rnafold. 另外, 通过随机选取20个tRNA序列, 从自由能和三叶草结构两个方面比较了前4种二级结构预测算法, 并运用t检验方法分析了自由能的统计学差别. 从三叶草结构来看, 以随机堆积方法最好, 其次是螺旋区最优堆积方法和Zuker算法, 以碱基最大配对方法最差. 最后, 分析了两种极小化自由能方法之间的差别.     

RNA模式分析进展

研究表明,R N A模式在基因表达调控方面起着重要作用.由于RN A模式不仅与初级序列有关,更多的表现为高级结构(一般为二级结构)的保守性,所以R N A模式的识别比D N A模式的识别要复杂的多.近十几年里,对R N A模式分析作了大量的计算方面的研究,包括:R N A结构的预测、识别和已知的类型相似的R N A模式、在一组功能或进化相关的基因中找出共同的R N A模式.这里对上述3个方面的计算方法的发展和研究进行了综述.     

Thermodynamic characterization of naturally occurring RNA tetraloops

Although tetraloops are one of the most frequently occurring secondary structure motifs in RNA, less than one-third of the 30 most frequently occurring RNA tetraloops have been thermodynamically characterized. Therefore, 24 stem–loop sequences containing common tetraloops were optically melted, and the thermodynamic parameters ΔH°, ΔS°, ΔG°_37, and T_M for each stem–loop were determined. These new experimental values, on average, are 0.7 kcal/mol different from the values predicted for these tetraloops using the model proposed by Vecenie CJ, Morrow CV, Zyra A, Serra MJ. 2006. Biochemistry 45: 1400–1407. The data for the 24 tetraloops reported here were then combined with the data for 28 tetraloops that were published previously. A new model, independent of terminal mismatch data, was derived to predict the free energy contribution of previously unmeasured tetraloops. The average absolute difference between the measured values and the values predicted using this proposed model is 0.4 kcal/mol. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA stem–loops containing tetraloops and, furthermore, should allow for improved prediction of secondary structure from sequence. It was also shown that tetraloops within the sequence 5′-GCCNNNNGGC-3′ are, on average, 0.6 kcal/mol more stable than the same tetraloop within the sequence 5′-GGCNNNNGCC-3′. More systemic studies are required to determine the full extent of non-nearest-neighbor effects on tetraloop stability.     

The tedious task of finding homologous noncoding RNA genes

User-driven in silico RNA homology search is still a nontrivial task. In part, this is the consequence of a limited precision of the computational tools in spite of recent exciting progress in this area, and to a certain extent, computational costs are still problematic in practice. An important, and as we argue here, dominating issue is the dependence on good curated (secondary) structural alignments of the RNAs. These are often hard to obtain, not so much because of an inherent limitation in the available data, but because they require substantial manual curation, an effort that is rarely acknowledged. Here, we qualitatively describe a realistic scenario for what a “regular user” (i.e., a nonexpert in a particular RNA family) can do in practice, and what kind of results are likely to be achieved. Despite the indisputable advances in computational RNA biology, the conclusion is discouraging: BLAST still works better or equally good as other methods unless extensive expert knowledge on the RNA family is included. However, when good curated data are available the recent development yields further improvements in finding remote homologs. Homology search beyond the reach of BLAST hence is not at all a routine task.     

RNA2D3D: A program for Generating,Viewing, and Comparing 3-Dimensional Models of RNA

Abstract

Using primary and secondary structure information of an RNA molecule, the program RNA2D3D automatically and rapidly produces a first-order approximation of a 3-dimensional conformation consistent with this information. Applicable to structures of arbitrary branching complexity and pseudoknot content, it features efficient interactive graphical editing for the removal of any overlaps introduced by the initial generating procedure and for making conformational changes favorable to targeted features and subsequent refinement. With emphasis on fast exploration of alternative 3D conformations, one may interactively add or delete base-pairs, adjacent stems can be coaxially stacked or unstacked, single strands can be shaped to accommodate special constraints, and arbitrary subsets can be defined and manipulated as rigid bodies. Compaction, whereby base stacking within stems is optimally extended into connecting single strands, is also available as a means of strategically making the structures more compact and revealing folding motifs. Subsequent refinement of the first-order approximation, of modifications, and for the imposing of tertiary constraints is assisted with standard energy refinement techniques. Previously determined coordinates for any part of the molecule are readily incorporated, and any part of the modeled structure can be output as a PDB or XYZ file. Illustrative applications in the areas of ribozymes, viral kissing loops, viral internal ribosome entry sites, and nanobiology are presented.     

Ligand-inducible formation of RNA pseudoknot

Here, we demonstrate that a series of naphthyridine derivatives, naphthyridine carbamate tetramer (NCTn), can bind to the RNA CGG/CGG triad comprised of two single-stranded regions of a hairpin loop and a tail. Complete suppression of the binding by a single mutation indicated simultaneous binding of NCTn between hairpin loop and single stranded tail, leading to the formation of NCTn-induced pseudoknot.     

General combinatorics of RNA secondary structure

The total number of RNA secondary structures of a given length with minimal hairpin loop length m(m>0) and with minimal stack length l(l>0) is computed, under the assumption that all base pairs can occur. Asymptotics are derived from the determination of recurrence relations of decomposition properties.