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伪结预测是RNA研究的一个难点问题。文中提出一种基于堆积协变信息与最小自由能的RNA伪结预测方法。该方法使用已知结构的RNA比对序列(ClustalW比对和结构比对)测试此方法, 侧重考虑相邻碱基对之间相互作用形成的堆积协变信息, 并结合最小自由能方法对碱基配对综合评分, 通过逐步迭代求得含伪结的RNA二级结构。结果表明, 此方法能正确预测伪结, 其平均敏感性和特异性优于参考算法, 并且结构比对的预测性能比ClustalW比对的预测性能更加稳定。文中同时讨论了不同协变信息权重因子对预测性能的影响, 发现权重因子比值在l1: l2=5:1时, 预测性能达到最优。     

基于堆积协变信息与最小自由能预测含伪结的RNA二级结构

RNA伪结预测是RNA研究的一个难点问题。文中提出一种基于堆积协变信息与最小自由能的RNA伪结预测方法。该方法使用已知结构的RNA比对序列(ClustalW比对和结构比对)测试此方法, 侧重考虑相邻碱基对之间相互作用形成的堆积协变信息, 并结合最小自由能方法对碱基配对综合评分, 通过逐步迭代求得含伪结的RNA二级结构。结果表明, 此方法能正确预测伪结, 其平均敏感性和特异性优于参考算法, 并且结构比对的预测性能比ClustalW比对的预测性能更加稳定。文中同时讨论了不同协变信息权重因子对预测性能的影响, 发现权重因子比值在l1: l2=5:1时, 预测性能达到最优。     

Structural Basis for the Discriminative Recognition of N6-Methyladenosine RNA by the Human YT521-B Homology Domain Family of Proteins

N⁶-Methyladenosine (m⁶A) is the most abundant internal modification in RNA and is specifically recognized by YT521-B homology (YTH) domain-containing proteins. Recently we reported that YTHDC1 prefers guanosine and disfavors adenosine at the position preceding the m⁶A nucleotide in RNA and preferentially binds to the GG(m⁶A)C sequence. Now we systematically characterized the binding affinities of the YTH domains of three other human proteins and yeast YTH domain protein Pho92 and determined the crystal structures of the YTH domains of human YTHDF1 and yeast Pho92 in complex with a 5-mer m⁶A RNA, respectively. Our binding and structural data revealed that the YTH domain used a conserved aromatic cage to recognize m⁶A. Nevertheless, none of these YTH domains, except YTHDC1, display sequence selectivity at the position preceding the m⁶A modification. Structural comparison of these different YTH domains revealed that among those, only YTHDC1 harbors a distinctly selective binding pocket for the nucleotide preceding the m⁶A nucleotide.     

Simulating the Folding Pathway of RNA Secondary Structure Using the Modified Ant Colony Algorithm

<正> A new method for simulating the folding pathway of RNA secondary structure using the modified ant colony algorithmis proposed.For a given RNA sequence,the set of all possible stems is obtained and the energy of each stem iscalculated and stored at the initial stage.Furthermore,a more realistic formula is used to compute the energy ofmulti-branch loop in the following iteration.Then a folding pathway is simulated,including such processes as constructionof the heuristic information,the rule of initializing the pheromone,the mechanism of choosing the initial andnext stem and the strategy of updating the pheromone between two different stems.Finally by testing RNA sequences withknown secondary structures from the public databases,we analyze the experimental data to select appropriate values forparameters.The measure indexes show that our procedure is more consistent with phylogenetically proven structures thansoftware RNAstructure sometimes and more effective than the standard Genetic Algorithm.     

Specificity of mRNA Folding and Its Association with Evolutionarily Adaptive mRNA Secondary Structures

The secondary structure is a fundamental feature of both non-coding RNAs (ncRNAs) and messenger RNAs (mRNAs). However, our understanding of the secondary structures of mRNAs, especially those of the coding regions, remains elusive, likely due to translation and the lack of RNA-binding proteins that sustain the consensus structure like those binding to ncRNAs. Indeed, mRNAs have recently been found to adopt diverse alternative structures, but the overall functional significance remains untested. We hereby approach this problem by estimating the folding specificity, i.e., the probability that a fragment of an mRNA folds back to the same partner once refolded. We show that the folding specificity of mRNAs is lower than that of ncRNAs and exhibits moderate evolutionary conservation. Notably, we find that specific rather than alternative folding is likely evolutionarily adaptive since specific folding is frequently associated with functionally important genes or sites within a gene. Additional analysis in combination with ribosome density suggests the ability to modulate ribosome movement as one potential functional advantage provided by specific folding. Our findings reveal a novel facet of the RNA structurome with important functional and evolutionary implications and indicate a potential method for distinguishing the mRNA secondary structures maintained by natural selection from molecular noise.     

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     

YTH Domain: A Family of N6-methyladenosine (m6A) Readers

Like protein and DNA, different types of RNA molecules undergo various modifications. Accumulating evidence suggests that these RNA modifications serve as sophisticated codes to mediate RNA behaviors and many important biological functions. N⁶-methyladenosine (m⁶A) is the most abundant internal RNA modification found in a variety of eukaryotic RNAs, including but not limited to mRNAs, tRNAs, rRNAs, and long non-coding RNAs (lncRNAs). In mammalian cells, m⁶A can be incorporated by a methyltransferase complex and removed by demethylases, which ensures that the m⁶A modification is reversible and dynamic. Moreover, m⁶A is recognized by the YT521-B homology (YTH) domain-containing proteins, which subsequently direct different complexes to regulate RNA signaling pathways, such as RNA metabolism, RNA splicing, RNA folding, and protein translation. Herein, we summarize the recent progresses made in understanding the molecular mechanisms underlying the m⁶A recognition by YTH domain-containing proteins, which would shed new light on m⁶A-specific recognition and provide clues to the future identification of reader proteins of many other RNA modifications.     

RNA Folding Argues Against a Hot-Start Origin of Life

Opinion is strongly divided on whether life arose on earth under hot or cold conditions, the hot-start and cold-start scenarios, respectively. The origin of life close to deep thermal vents appears as the majority opinion among biologists, but there is considerable biochemical evidence that high temperatures are incompatible with an RNA world. To be functional, RNA has to fold into a three-dimensional structure. We report both theoretical and experimental results on RNA folding and show that (as expected) hot conditions strongly reduce RNA folding. The theoretical results come from energy-minimization calculations of the average extent of folding of RNA, mainly from 0–90°C, for both random sequences and tRNA sequences. The experimental results are from circular-dichroism measurements of tRNA over a similar range of temperatures. The quantitative agreement between calculations and experiment is remarkable, even to the shape of the curves indicating the cooperative nature of RNA folding and unfolding. These results provide additional evidence for a lower temperature stage being necessary in the origin of life. Received: 1 March 2000 / Accepted: 14 June 2000     

Nanomanipulation of Single RNA Molecules by Optical Tweezers

A large portion of the human genome is transcribed but not translated. In this post genomic era, regulatory functions of RNA have been shown to be increasingly important. As RNA function often depends on its ability to adopt alternative structures, it is difficult to predict RNA three-dimensional structures directly from sequence. Single-molecule approaches show potentials to solve the problem of RNA structural polymorphism by monitoring molecular structures one molecule at a time. This work presents a method to precisely manipulate the folding and structure of single RNA molecules using optical tweezers. First, methods to synthesize molecules suitable for single-molecule mechanical work are described. Next, various calibration procedures to ensure the proper operations of the optical tweezers are discussed. Next, various experiments are explained. To demonstrate the utility of the technique, results of mechanically unfolding RNA hairpins and a single RNA kissing complex are used as evidence. In these examples, the nanomanipulation technique was used to study folding of each structural domain, including secondary and tertiary, independently. Lastly, the limitations and future applications of the method are discussed.     

Effects of N2,N2-dimethylguanosine on RNA structure and stability: crystal structure of an RNA duplex with tandem m2 2G:A pairs

Methylation of the exocyclic amino group of guanine is a relatively common modification in rRNA and tRNA. Single methylation (N²-methylguanosine, m²G) is the second most frequently encountered nucleoside analog in Escherichia coli rRNAs. The most prominent case of dual methylation (N²,N²-dimethylguanosine, m²₂G) is found in the majority of eukaryotic tRNAs at base pair m²₂G26:A44. The latter modification eliminates the ability of the N² function to donate in hydrogen bonds and alters its pairing behavior, notably vis-à-vis C. Perhaps a less obvious consequence of the N²,N²-dimethyl modification is its role in controlling the pairing modes between G and A. We have determined the crystal structure of a 13-mer RNA duplex with central tandem m²₂G:A pairs. In the structure both pairs adopt an imino-hydrogen bonded, pseudo-Watson–Crick conformation. Thus, the sheared conformation frequently seen in tandem G:A pairs is avoided due to a potential steric clash between an N²-methyl group and the major groove edge of A. Additionally, for a series of G:A containing self-complementary RNAs we investigated how methylation affects competitive hairpin versus duplex formation based on UV melting profile analysis.     

Structures of Human ALKBH5 Demethylase Reveal a Unique Binding Mode for Specific Single-stranded N6-Methyladenosine RNA Demethylation

N⁶-Methyladenosine (m⁶A) is the most prevalent internal RNA modification in eukaryotes. ALKBH5 belongs to the AlkB family of dioxygenases and has been shown to specifically demethylate m⁶A in single-stranded RNA. Here we report crystal structures of ALKBH5 in the presence of either its cofactors or the ALKBH5 inhibitor citrate. Catalytic assays demonstrate that the ALKBH5 catalytic domain can demethylate both single-stranded RNA and single-stranded DNA. We identify the TCA cycle intermediate citrate as a modest inhibitor of ALKHB5 (IC₅₀, ∼488 μm). The structural analysis reveals that a loop region of ALKBH5 is immobilized by a disulfide bond that apparently excludes the binding of dsDNA to ALKBH5. We identify the m⁶A binding pocket of ALKBH5 and the key residues involved in m⁶A recognition using mutagenesis and ITC binding experiments.     

Roles of N6‐methyladenosine (m6A) RNA modifications in urological cancers

Epigenetics has long been a hot topic in the field of scientific research. The scope of epigenetics usually includes chromatin remodelling, DNA methylation, histone modifications, non‐coding RNAs and RNA modifications. In recent years, RNA modifications have emerged as important regulators in a variety of physiological processes and in disease progression, especially in human cancers. Among the various RNA modifications, m⁶A is the most common. The function of m⁶A modifications is mainly regulated by 3 types of proteins: m⁶A methyltransferases (writers), m⁶A demethylases (erasers) and m⁶A‐binding proteins (readers). In this review, we focus on RNA m⁶A modification and its relationship with urological cancers, particularly focusing on its roles and potential clinical applications.     

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.     

Non-nearest-neighbor dependence of the stability for RNA group II single-nucleotide bulge loops

Thirty-one RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters ΔH°, ΔS°, ΔG°(37), and T(M) for each sequence were determined. The bulge loops were of the group II variety, where the bulged nucleotide is identical to one of its nearest neighbors, leading to ambiguity as to the exact position of the bulge. The data were used to develop a model to predict the free energy of an RNA duplex containing a single-nucleotide bulge. The destabilization of the duplex by the bulge was primarily related to the stability of the stems adjacent to the bulge. Specifically, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. Since there is an ambiguity of the bulge position for group II bulges, several different stem combinations are possible. The destabilization of group II bulge loops is similar to the destabilization of group I bulge loops, if the second least stable stem is used to predict the influence of the group II bulge. In-line structure probing of the group II bulge loop embedded in a hairpin indicates that the bulged nucleotide is the one positioned farther from the hairpin loop.     

A range of complex probabilistic models for RNA secondary structure prediction that includes the nearest-neighbor model and more

The standard approach for single-sequence RNA secondary structure prediction uses a nearest-neighbor thermodynamic model with several thousand experimentally determined energy parameters. An attractive alternative is to use statistical approaches with parameters estimated from growing databases of structural RNAs. Good results have been reported for discriminative statistical methods using complex nearest-neighbor models, including CONTRAfold, Simfold, and ContextFold. Little work has been reported on generative probabilistic models (stochastic context-free grammars [SCFGs]) of comparable complexity, although probabilistic models are generally easier to train and to use. To explore a range of probabilistic models of increasing complexity, and to directly compare probabilistic, thermodynamic, and discriminative approaches, we created TORNADO, a computational tool that can parse a wide spectrum of RNA grammar architectures (including the standard nearest-neighbor model and more) using a generalized super-grammar that can be parameterized with probabilities, energies, or arbitrary scores. By using TORNADO, we find that probabilistic nearest-neighbor models perform comparably to (but not significantly better than) discriminative methods. We find that complex statistical models are prone to overfitting RNA structure and that evaluations should use structurally nonhomologous training and test data sets. Overfitting has affected at least one published method (ContextFold). The most important barrier to improving statistical approaches for RNA secondary structure prediction is the lack of diversity of well-curated single-sequence RNA secondary structures in current RNA databases.     

Interaction of rice and human SRP19 polypeptides with signal recognition particle RNA

The signal recognition particle (SRP) controls the transport of secretory proteins into and across lipid bilayers. SRP-like ribonucleoprotein complexes exist in all organisms, including plants. We characterized the rice SRP RNA and its primary RNA binding protein, SRP19. The secondary structure of the rice SRP RNA was similar to that found in other eukaryotes; however, as in other plant SRP RNAs, a GUUUCA hexamer sequence replaced the highly conserved GNRA-tetranucleotide loop motif at the apex of helix 8. The small domain of the rice SRP RNA was reduced considerably. Structurally, rice SRP19 lacked two small regions that can be present in other SRP19 homologues. Conservative structure prediction and site-directed mutagenesis of rice and human SRP19 polypeptides indicated that binding to the SRP RNAs occurred via a loop that is present in the N-domain of both proteins. Rice SRP19 protein was able to form a stable complex with the rice SRP RNA in vitro. Furthermore, heterologous ribonucleoprotein complexes with components of the human SRP were assembled, thus confirming a high degree of structural and functional conservation between plant and mammalian SRP components.