The Role of Structural Elements in Dimerization of RNA Fragments in Avian Retroviruses

The slipped loop structure, earlier identified as an unusual DNA structure, was found to be a possible element of the RNA folding. In order to experimentally test this suggestion, model oligoribonucleotides capable of forming the SLS were synthesized. Treatment of the oligoribonucleotides with nuclease S1 and RNases specific for single- and double-stranded RNA demonstrated the steric possibility of SLS formation. To determine the possible functional role of SLS-RNA, various naturally occurring RNAs were screened in silico. Among the most interesting findings were dimerization initiation sites of avian retroviral genomic RNAs. Analysis of RNA from 31 viruses showed that formation of the intermolecular SLS during RNA dimerization is theoretically possible, competing with the formation of an alternative hairpin structure. Identification of the secondary structure of selected RNA dimers employing nuclease digestion techniques as well as covariance analysis of the retroviral RNA dimerization initiation site sequences were used to show that the alternative conformation (loop–loop interaction of two hairpins, or kissing hairpins) is the most preferred. Alternative structures and conformational transitions in RNA dimerization mechanisms in avian retroviruses are discussed.     

RNA剪接因子结构与功能研究进展

RNA剪接是一个多步骤、形成多种中间状态复合物的复杂过程.尽管在已经发现的一百多种pre-mRNA剪接相关因子中,仅研究了约8%相关蛋白质的空间结构,已充分显示对剪接相关因子三维晶体结构以及溶液结构的测定与研究,在理解RNA剪接的复杂机理以及生物学特性中具有不可替代的重要意义.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

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.     

From knotted to nested RNA structures: a variety of computational methods for pseudoknot removal

Pseudoknots are abundant in RNA structures. Many computational analyses require pseudoknot-free structures, which means that some of the base pairs in the knotted structure must be disregarded to obtain a nested structure. There is a surprising diversity of methods to perform this pseudoknot removal task, but these methods are often poorly described and studies can therefore be difficult to reproduce (in part, because different procedures may be intuitively obvious to different investigators). Here we provide a variety of algorithms for pseudoknot removal, some of which can incorporate sequence or alignment information in the removal process. We demonstrate that different methods lead to different results, which might affect structure-based analyses. This work thus provides a starting point for discussion of the extent to which these different methods recapture the underlying biological reality. We provide access to reference implementations through a web interface (at http://www.ibi.vu.nl/programs/k2nwww), and the source code is available in the PyCogent project.     

Artificial ribonucleases: from combinatorial libraries to efficient catalysts of RNA cleavage

Combinatorial libraries of small organic compounds capable of cleaving RNA were synthesized. The compounds contain benzene ring substituted with two residues of bis quaternary salt of diazabicyclo[2.2.2]octane (DABCO) bearing hydrophobic fragments of different length and structure, attached to DABCO at the bridge position. These compounds, lacking traditional functionalities involved in transesterification reaction, exhibit pronounced RNA cleavage activity. To identify the most active artificial ribonucleases, sublibraries and truncated libraries, containing compounds lacking one of substituents were synthesized. Analysis of ribonuclease activity of truncated libraries resulted in identification of the most active compounds, which are characterized by the presence of at least one long oligomethylene substituent.     

Structural Constraints and Mutational Bias in the Evolutionary Restoration of a Severe Deletion in RNA Phage MS2

A 4-nucleotide (nt) deletion was made in the 36-nt-long intercistronic region separating the coat and replicase genes of the single-stranded RNA phage MS2. This region is the focus of several RNA structures conferring high fitness. One such element is the operator hairpin, which, in the course of infection, will bind a coat-protein dimer, thereby precluding further replicase synthesis and initiating encapsidation. Another structure is a long-distance base pairing (MJ) controlling replicase expression. The 4-nt deletion does not directly affect the operator hairpin but it disrupts the MJ pairing. Its main effect, however, is a frame shift in the overlapping lysis gene. This gene starts in the upstream coat gene, runs through the 36-nt-long intercistronic region, and ends in the downstream replicase cistron. Here we report and interpret the spectrum of solutions that emerges when the crippled phage is evolved. Four different solutions were obtained by sequencing 40 plaques. Three had cured the frame shift in the lysis gene by inserting one nt in the loop of the operator hairpin causing its inactivation. Yet these low-fitness revertants could further improve themselves when evolved. The inactivated operator was replaced by a substitute and thereafter these revertants found several ways to restore control over the replicase gene. To allow for the evolutionary enrichment of low-probability but high-fitness revertants, we passaged lysate samples before plating. Revertants obtained in this way also restored the frame shift, but not at the expense of the operator. By taking larger and larger lysates samples for such bulk evolution, ever higher-fitness and lower-frequency revertants surfaced. Only one made it back to wild type. As a rule, however, revertants moved further and further away from the wild-type sequence because restorative mutations are, in the majority of cases, selected for their capacity to improve the phenotype by optimizing one of several potential alternative RNA foldings that emerge as a result of the initial deletion. This illustrates the role of structural constraints which limit the path of subsequent restorative mutations. [Reviewing Editor: Dr. John Hulsenbeck]     

The promise of cryo-EM to explore RNA structural dynamics

Conformational dynamics are essential to macromolecular function. This is certainly true of RNA, whose ability to undergo programmed conformational dynamics is essential to create and regulate complex biological processes. However, methods to easily and simultaneously interrogate both the structure and conformational dynamics of fully functional RNAs in isolation and in complex with proteins have not historically been available. Due to its ability to image and classify single particles, cryogenic electron microscopy (cryo-EM) has the potential to address this gap and may be particularly amenable to exploring structural dynamics within the three-dimensional folds of biologically active RNAs. We discuss the possibilities and current limitations of applying cryo-EM to simultaneously study RNA structure and conformational dynamics, and present one example that illustrates this (as of yet) not fully realized potential.     

ProbKnot: Fast prediction of RNA secondary structure including pseudoknots

It is a significant challenge to predict RNA secondary structures including pseudoknots. Here, a new algorithm capable of predicting pseudoknots of any topology, ProbKnot, is reported. ProbKnot assembles maximum expected accuracy structures from computed base-pairing probabilities in O(N²) time, where N is the length of the sequence. The performance of ProbKnot was measured by comparing predicted structures with known structures for a large database of RNA sequences with fewer than 700 nucleotides. The percentage of known pairs correctly predicted was 69.3%. Additionally, the percentage of predicted pairs in the known structure was 61.3%. This performance is the highest of four tested algorithms that are capable of pseudoknot prediction. The program is available for download at: http://rna.urmc.rochester.edu/RNAstructure.html.     

Evolution and Taxonomy of Positive-Strand RNA Viruses: Implications of Comparative Analysis of Amino Acid Sequences

Abstract

Despite the rapid mutational change that is typical of positive-strand RNA viruses, enzymes mediating the replication and expression of virus genomes contain arrays of conserved sequence motifs. Proteins with such motifs include RNA-dependent RNA polymerase, putative RNA helicase, chymotrypsin-like and papain-like proteases, and methyltransferases. The genes for these proteins form partially conserved modules in large subsets of viruses. A concept of the virus genome as a relatively evolutionarily stable “core” of housekeeping genes accompanied by a much more flexible “shell” consisting mostly of genes coding for virion components and various accessory proteins is discussed. Shuffling of the “shell” genes including genome reorganization and recombination between remote groups of viruses is considered to be one of the major factors of virus evolution.

Multiple alignments for the conserved viral proteins were constructed and used to generate the respective phylogenetic trees. Based primarily on the tentative phylogeny for the RNA-dependent RNA polymerase, which is the only universally conserved protein of positive-strand RNA viruses, three large classes of viruses, each consisting of distinct smaller divisions, were delineated. A strong correlation was observed between this grouping and the tentative phylogenies for the other conserved proteins as well as the arrangement of genes encoding these proteins in the virus genome. A comparable correlation with the polymerase phylogeny was not found for genes encoding virion components or for genome expression strategies. It is surmised that several types of arrangement of the “shell” genes as well as basic mechanisms of expression could have evolved independently in different evolutionary lineages.

The grouping revealed by phylogenetic analysis may provide the basis for revision of virus classification, and phylogenetic taxonomy of positive-strand RNA viruses is outlined. Some of the phylogenetically derived divisions of positive-strand RNA viruses also include double-stranded RNA viruses, indicating that in certain cases the type of genome nucleic acid may not be a reliable taxonomic criterion for viruses.

Hypothetical evolutionary scenarios for positive-strand RNA viruses are proposed. It is hypothesized that all positive-strand RNA viruses and some related double-stranded RNA viruses could have evolved from a common ancestor virus that contained genes for RNA-dependent RNA polymerase, a chymotrypsin-related protease that also functioned as the capsid protein, and possibly an RNA helicase.     

A DEAD-box RNA helicase promotes thermodynamic equilibration of kinetically trapped RNA structures in vivo

RNAs must assemble into specific structures in order to carry out their biological functions, but in vitro RNA folding reactions produce multiple misfolded structures that fail to exchange with functional structures on biological time scales. We used carefully designed self-cleaving mRNAs that assemble through well-defined folding pathways to identify factors that differentiate intracellular and in vitro folding reactions. Our previous work showed that simple base-paired RNA helices form and dissociate with the same rate and equilibrium constants in vivo and in vitro. However, exchange between adjacent secondary structures occurs much faster in vivo, enabling RNAs to quickly adopt structures with the lowest free energy. We have now used this approach to probe the effects of an extensively characterized DEAD-box RNA helicase, Mss116p, on a series of well-defined RNA folding steps in yeast. Mss116p overexpression had no detectable effect on helix formation or dissociation kinetics or on the stability of interdomain tertiary interactions, consistent with previous evidence that intracellular factors do not affect these folding parameters. However, Mss116p overexpression did accelerate exchange between adjacent helices. The nonprocessive nature of RNA duplex unwinding by DEAD-box RNA helicases is consistent with a branch migration mechanism in which Mss116p lowers barriers to exchange between otherwise stable helices by the melting and annealing of one or two base pairs at interhelical junctions. These results suggest that the helicase activity of DEAD-box proteins like Mss116p distinguish intracellular RNA folding pathways from nonproductive RNA folding reactions in vitro and allow RNA structures to overcome kinetic barriers to thermodynamic equilibration in vivo.     

A low-density cDNA microarray with a unique reference RNA: pattern recognition analysis for IFN efficacy prediction to HCV as a model

We have designed and established a low-density (295 genes) cDNA microarray for the prediction of IFN efficacy in hepatitis C patients. To obtain a precise and consistent microarray data, we collected a data set from three spots for each gene (mRNA) and using three different scanning conditions. We also established an artificial reference RNA representing pseudo-inflammatory conditions from established hepatocyte cell lines supplemented with synthetic RNAs to 48 inflammatory genes. We also developed a novel algorithm that replaces the standard hierarchical-clustering method and allows handling of the large data set with ease. This algorithm utilizes a standard space database (SSDB) as a key scale to calculate the Mahalanobis distance (MD) from the center of gravity in the SSDB. We further utilized sMD (divided by parameter k: MD/k) to reduce MD number as a predictive value. The efficacy prediction of conventional IFN mono-therapy was 100% for non-responder (NR) vs. transient responder (TR)/sustained responder (SR) (P < 0.0005). Finally, we show that this method is acceptable for clinical application.     

不同体重肩斑银鳞蛛和大腹园蛛圆网的结构特征

对在自然条件下不同体重肩斑银鳞蛛Leucauge blanda和大腹园蛛Araneus ventricosus所织圆网的结构特征进行测量研究.结果表明:体重小于35.1 mg的肩斑银鳞蛛所织圆网的捕丝长度、捕食面面积、捕丝间距和半径丝根数均与个体体重呈显著正相关,而体重大于35.1 mg的个体中,这种关系并不显著,但其圆网的上、下部捕丝长度比与体重呈显著负相关,即圆网随个体的体重增加而表现出更强的不对称性;体重小于144.9 mg的大腹园蛛所织圆网的捕丝长度和捕食面面积均与个体体重呈显著正相关,体重小于103.8 mg的大腹园蛛所织圆网的半径丝根数与个体体重呈显著正相关,而大于这一体重分界值的个体中,这种关系同样不显著.大腹园蛛圆网的平均捕丝间距与体重未呈现出相关关系,体重大于85.4 mg的个体中,其网的上、下部捕丝长度比与体重呈显著负相关.两种蜘蛛圆网结构特征的变化及圆网结构特征与个体体重关系变化的不同体重分界值,可能反映了它们在不同生境下不同生长阶段的捕食投入与捕食策略.     

Structural Features of Human Monoamine Oxidase A Elucidated from cDNA and Peptide Sequences

Monoamine oxidase (MAO), an important enzyme for the degradation of amine neurotransmitters, has been implicated in neuropsychiatric illness. The amino acid sequence for one form of the enzyme, MAO-A, has been deduced from human cDNA clones and verified against proteolytic peptides. The covalent binding site for the flavin adenine dinucleotide (FAD) cofactor is near the C-terminal region. The presence of features characteristic of the ADP-binding fold suggests that the N-terminal region is also involved in the binding of FAD. These cDNAs should facilitate the study of the structure, function, and intracellular targeting of MAO, as well as the analysis of its expression in normal and pathological states.     

The plasticity of a structural motif in RNA: Structural polymorphism of a kink turn as a function of its environment

The k-turn is a widespread structural motif that introduces a tight kink into the helical axis of double-stranded RNA. The adenine bases of consecutive G•A pairs are directed toward the minor groove of the opposing helix, hydrogen bonding in a typical A-minor interaction. We show here that the available structures of k-turns divide into two classes, depending on whether N3 or N1 of the adenine at the 2b position accepts a hydrogen bond from the O2′ at the −1n position. There is a coordinated structural change involving a number of hydrogen bonds between the two classes. We show here that Kt-7 can adopt either the N3 or N1 structures depending on environment. While it has the N1 structure in the ribosome, on engineering it into the SAM-I riboswitch, it changes to the N3 structure, resulting in a significant alteration in the trajectory of the helical arms.     

A dimeric RNA quadruplex architecture comprised of two G:G(:A):G:G(:A) hexads,G:G:G:G tetrads and UUUU loops

Using CD and NMR, we determined the structure of an RNA oligomer, r(GGAGGUUUUGGAGG) (R14), comprising two GGAGG segments joined by a UUUU segment. A modified quadruplex structure was observed for r(GGAGGUUUUGGAGG) in solution even in the absence of K(+). An unusually stable dimeric RNA quadruplex architecture formed from two strands of r(GGAGGUUUUGGAGG) at low K(+) concentration is reported here. In each strand of r(GGAGGUUUUGGAGG), two sets of successive turns in the GGAGG segments and turns at both ends of the UUUU loops drive four G-G steps to align in a parallel manner, a core with two stacked G-tetrads being formed. Two adenine bases bind to two edges of one G:G:G:G tetrad through the sheared G:A mismatch augmenting the tetrad into a G:G(:A):G:G(:A) hexad. Thus, one molecule of r(GGAGGUUUUGGAGG) folds into a modified quadruplex comprising a G:G:G:G tetrad, a UUUU double-chain reversal loop and a G:G(:A):G:G(:A) hexad. Two such molecules further associate by stacking through the dimeric hexad-hexad interface with a rotational symmetry. The ribose rings of most nucleotides take S (close to C2'-endo) puckering, which is unusual for an RNA. K(+) can increase the stability of this quadruplex structure; the number of bound K(+) was estimated from the results of the titration experiment. Besides G:G and G:A mismatches, a network of hydrogen bonds including O4'-NH(2) and C-H..O hydrogen bonds, and the extensive base stacking contribute to the high thermodynamic stability of R14. Our results could provide the stereochemical and thermodynamic basis for elucidating the biological role of the GGAGG-containing RNA segments abundantly existing in various RNAs. Relevance to quadruplex-mediated mRNA-FMRP binding and HIV-1 genome RNA dimerization is discussed.     

Fuzzy Kernel Clustering of RNA Secondary Structure Ensemble Using a Novel Similarity Metric

Abstract

Measuring the (dis)similarity between RNA secondary structures is critical for the study of RNA secondary structures and has implications to RNA functional characterization. Although a number of methods have been developed for comparing RNA structural similarities, their applications have been limited by the complexity of the required computation. In this paper, we present a novel method for comparing the similarity of RNA secondary structures generated from the same RNA sequence, i.e., a secondary structure ensemble, using a matrix representation of the RNA structures. Relevant features of the RNA secondary structures can be easily extracted through singular value decomposition (SVD) of the representing matrices. We have mapped the feature vectors of the singular values to a kernel space, where (dis)similarities among the mapped feature vectors become more evident, making clustering of RNA secondary structures easier to handle. The pair-wise comparison of RNA structures is achieved through computing the distance between the singular value vectors in the kernel space. We have applied a fuzzy kernel clustering method, using this similarity metric, to cluster the RNA secondary structure ensembles. Our application results suggest that our fuzzy kernel clustering method is highly promising for classifications of RNA structure ensembles, because of its low computational complexity and high clustering accuracy.     

银耳碱提孢子多糖A-BTF的分离与结构研究

银耳孢子发酵粉(Tremella fuciformisBerk),用热水煮提后,除去水溶性多糖。其沉淀用1M NaOH提取,Sevage法除去蛋白,用乙醇沉淀得到粗多糖。粗多糖经DEAE-32-cellulose和sephadex G-200反复分离后,纯化得到分布均一的多糖A-BTF。HPGPC测定A-BTF的分子量为67000,糖组成分析显示主要由葡萄糖组成。多糖A-BTF的甲基化产物,经水解、还原、乙酰化,通过GC-MS分析表明,主要含有1,6连接的葡萄糖和1,3,6连接的甘露糖,另外还有1,4连接的葡萄糖少量的半乳糖和1-NH2-来苏糖,末端为端基连接的葡萄糖。