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.     

Preference for Ribose Over Deoxyribose in Loop-Closing Base Pairs of Extra Stable Nucleic Acid Hairpins

We have investigated the effect of switching ribose to deoxyribose at the closing base-pair of an extra-stable RNA hairpin. Specifically, we studied the sequence 5′-GGAC(UUCG)GUCC, a dodecanucleotide that folds into a well-characterized, “extra stable” RNA hairpin structure. Recently, we showed that hairpins containing a 2′,5′-linked (UUCG) loop instead of the native 3′,5′-linked loop also exhibit extra-stability (Hannoush and Damha, J. Am. Chem. Soc., 2001, 123, 12368–12374). In this article, we show that the ribose units located at the loop-closing positions (i.e., rC ₄ and rG ₉ ) contribute significantly to the stabilization of RNA hairpins, particularly those containing the 3′,5′-UUCG loop. Interestingly, the requirement of rC4 and rG9 is more relaxed for DNA hairpins containing the 2′,5′-UUCG loop and, in fact, they may be replaced altogether (ribose → deoxyribose) without affecting stability. The results broaden our understanding of the behavior of highly stable (UUCG) hairpin loops and how they respond to structural perturbation of the loop-closing base pairs.     

ISFOLD: Structure Prediction of Base Pairs in Non-Helical RNA Motifs from Isostericity Signatures in Their Sequence Alignments

Abstract

The existence and identity of non-Watson-Crick base pairs (bps) within RNA bulges, internal loops, and hairpin loops cannot reliably be predicted by existing algorithms. We have developed the Isfold (Isosteric Folding) program as a tool to examine patterns of nucleotide substitutions from sequence alignments or mutation experiments and identify plausible bp interactions. We infer these interactions based on the observation that each non-Watson-Crick bp has a signature pattern of isosteric substitutions where mutations can be made that preserve the 3D structure. Isfold produces a dynamic representation of predicted bps within defined motifs in order of their probabilities. The software was developed under Windows XP, and is capable of running on PC and MAC with Matlab 7.1 (SP3) or higher. A PC standalone version that does not require Matlab also is available. This software and a user manual are freely available at www.ucsf.edu/frankel/isfold.     

Post Hartree-Fock Studies of the Canonical Watson-Crick DNA Base Pairs: Molecular Structure and the Nature of Stability

Abstract

Gas-phase gradient optimization was carried out on the canonical Watson-Crick DNA base pairs using the second-order Møller-Plesset perturbation method at the 6–31G(d) and 6- 31G(d,p) basis sets. It is detected that full geometry optimization at the MP2 level leads to an intrinsically nonplanar propeller-twisted and buckled geometry of G-C and A-T base pairs; while HF and DFT methods predict perfect planar or almost planar geometry of the base pairs. Supposedly the nonplanarity of the pairs is caused by pyramidalization of the amino nitrogen atoms, which is underestimated by the HF and DFT methods. This justifies the importance of geometry optimization at the MP2 level for obtaining reliable prediction of the charge distribution, molecular dipole moments and geometrical structure of the base pairs. The Morokuma-Kitaura and the Reduced Variational Space methods of the decomposition for molecular HF interaction energies were used for investigation of the hydrogen bonding in the Watson-Crick base pairs. It is shown that the HF stability of the hydrogen-bonded DNA base pairs originates mainly from electrostatic interactions. At the same time, the calculated magnitude of the second order intramolecular correlation correction to the Coulomb energy showed that electron correlation reduces the contribution of the electrostatic term to the attractive interaction for the A-T and G-C base pairs. Polarization, charge transfer and dispersion interactions also make considerable contribution to the attraction energy of bases.     

A New Bias Site for Epigenetic Modifications: How Non-Canonical GC Base Pairs Favor Mechanochemical Cleavage of DNA

Properties of non-canonical GC base pairs and their relations with mechanochemical cleavage of DNA are analyzed. A hypothesis of the involvement of the transient GC wobble base pairs both in the mechanisms of the mechanochemical cleavage of DNA and epigenetic mechanisms involving of 5-methylcytosine, is proposed. The hypothesis explains the increase in the frequency of the breaks of the sugar-phosphate backbone of DNA after cytosines, the asymmetric character of these breaks, and an increase in break frequency in CpG after cytosine methylation. As an alternative hypothesis, probable implication of GC⁺ Hoogsteen base pairs is considered, which now exemplify the best-studied non-canonical GC base pairs in the DNA double helix. Also see the video abstract here https://youtu.be/EUunVWL0ptw.     

Lagrangian Coherent Structures in Tandem Flapping Wing Hovering

Lagrangian Coherent Structures (LCS) of tandem wings hovering in an inclined stroke plane is studied using ImmersedBoundary Method (IBM) by solving two dimensional (2D) incompressible Navier-Stokes equations.Coherent structures responsible for the force variation are visualized by calculating Finite Time Lyapunov Exponents (FTLE),and vorticity contours.LCS is effective in determining the vortex boundaries,flow separation,and the wing-vortex interactions accurately.The effects of inter-wing distance and phase difference on the force generation are studied.Results show that in-phase stroking generates maximum vertical force and counter-stroking generates the least vertical force.In-phase stroking generates a wake with swirl,and counter stroking generates a wake with predominant vertical velocity.Counter stroking aids the stability of the body in hovering.As the hindwing operates in the wake of the forewing,due to the reduction in the effective Angle of Attack (AoA),the hindwing generates lesser force than that of a single flapping wing.     

Crystal Structures and RNA-binding Properties of the RNA Recognition Motifs of Heterogeneous Nuclear Ribonucleoprotein L: INSIGHTS INTO ITS ROLES IN ALTERNATIVE SPLICING REGULATION

Heterogeneous nuclear ribonucleoprotein L (hnRNP L) is an abundant RNA-binding protein implicated in many bioprocesses, including pre-mRNA processing, mRNA export of intronless genes, internal ribosomal entry site-mediated translation, and chromatin modification. It contains four RNA recognition motifs (RRMs) that bind with CA repeats or CA-rich elements. In this study, surface plasmon resonance spectroscopy assays revealed that all four RRM domains contribute to RNA binding. Furthermore, we elucidated the crystal structures of hnRNP L RRM1 and RRM34 at 2.0 and 1.8 Å, respectively. These RRMs all adopt the typical β1α1β2β3α2β4 topology, except for an unusual fifth β-strand in RRM3. RRM3 and RRM4 interact intimately with each other mainly through helical surfaces, leading the two β-sheets to face opposite directions. Structure-based mutations and surface plasmon resonance assay results suggested that the β-sheets of RRM1 and RRM34 are accessible for RNA binding. FRET-based gel shift assays (FRET-EMSA) and steady-state FRET assays, together with cross-linking and dynamic light scattering assays, demonstrated that hnRNP L RRM34 facilitates RNA looping when binding to two appropriately separated binding sites within the same target pre-mRNA. EMSA and isothermal titration calorimetry binding studies with in vivo target RNA suggested that hnRNP L-mediated RNA looping may occur in vivo. Our study provides a mechanistic explanation for the dual functions of hnRNP L in alternative splicing regulation as an activator or repressor.     

The complete nucleotide sequence of a 16S ribosomal RNA gene from tobacco chloroplasts: Recombinant plasmid; sequence homology; stem and loop structure

The complete nucleotide sequence of a 16S ribosomal RNA gene from tobacco chloroplasts has been determined. This nucleotide sequence has 96% homology with that of maize chloroplast 16S rRNA gene and 74% homology with that of Escherichia coli16S gene.The 3′ terminal region of this gene contains the sequence ACCTCC which is complementary to sequences found at the 5′ termini of prokaryotic mRNAs.The large stem and loop structure can be constructed from the sequences surrounding the 5′ and 3′ ends of the 16S gene. These observations demonstrate the prokaryotic nature of chloroplast 16S rRNA.     

基于古菌酪氨酰tRNA合成酶非天然氨基酸插入的研究进展

构筑蛋白质的编码信息存在于高度保守的密码子表中,而生物体仅利用20种天然氨基酸,就能排列组合出不同的蛋白质来行使多种生物学功能。通过合成生物学的飞速发展,使得在蛋白质合成中可控地引入非天然氨基酸成为可能。这极大地拓展了蛋白质的结构和功能,并为生物学工具的开发和生物生理过程的研究提供了便利。具有活性基团的非天然氨基酸可以广泛地应用于蛋白质结构研究、蛋白质功能调控以及新型生物材料构建和医药研发等诸多领域。基因密码子拓展技术利用正交翻译系统,通过重新分配密码子改造中心法则,可以在蛋白质的指定位点引入非天然氨基酸。系统地介绍了目前提升密码子拓展技术插入非天然氨基酸效率的方法,包括tRNA以及氨酰tRNA合成酶的各种突变方法和翻译辅助因子的改造。汇总了利用古细菌酪氨酰tRNA合成酶插入的非天然氨基酸和突变位点并总结了密码子拓展技术在生物医药领域的前沿进展。最后讨论了该项技术目前所面临的挑战,如可利用的密码子数量不多、正交翻译系统的种类有限和非天然氨基酸多插效率低下。希望能够帮助研究者建立适合的非天然氨基酸插入方法并推动密码子拓展技术进一步发展。     

Non-Watson-Crick base associations in DNA and RNA revealed by single crystal x-ray diffraction methods: Mismatches,modified bases,and nonduplex DNA

Single crystal x-ray diffraction methods have been used to characterize numerous oligonucleotide structures, providing valuable information on the fine structure of DNA, oligonucleotide hydration, interactions with small molecule ligands and proteins. There has been a particular focus on nonstandard base associations and a number of groups have sought to characterize different non-Watson-Crick base pairs to further the understanding of their influence on the structure of duplex DNA and RNA, and to investigate which structural features might be utilized by enzymes in recognition and repair of these errors in DNA. Bases that have been chemically damaged by mutagenic or carcinogenic agents have distinctive modified hydrogen-bonding patterns and these have been investigated. The structure determination of a series of nonduplex DNA structures including examples of a triplex, quadruplexes, and a novel DNA loop have recently been published. In this article we survey the structures of a series of non-Watson-Crick base associations in duplex DNA and RNA. We show how nonstandard base pairs, base triads, and tetrads play an important role in stabilizing nonduplex structures. © 1997 John Wiley & Sons, Inc. Biopoly 44: 91–103, 1997     

Detailed physical mapping of the ribosomal RNA genes of Bacillus subtilis

Characterization of patterns of ribosomal RNA (rRNA) homology with restriction digests of Bacillus subtilis 168 chromosomal DNA and with cloned DNA sequences has resulted in the construction of a physical map of the rRNA gene sets. There are two types of gene sets which differ in the size of "spacer" DNA sequences separating the 16S and 23S rRNA determinants. It was estimated that there are ten rRNA gene sets on the B. subtilis chromosome.     

文法推断RNA二级结构的研究进展

生物序列可看成是一种语言,通过计算语言学的方法理解生物序列的内涵是近年来研究的热点,本文综述了文法推断RNA二级结构的基本原理,研究历史和现状,阐述了文法推断RNA二级结构的理论模型和算法,列举了一些有代表性的预测方法,总结了存在的问题并展望了研究的趋势。     

On Spectral Identification of DNA-Base Pairs Polymorphism

The results of MNDO-PM3 theoretical study of H-bonds liability of watson-crick base pairs are discussed. Some microwave and IR spectral criteria are suggested for identification of hidden polymorphism of DNA base pairing.     

MBNL binds similar RNA structures in the CUG repeats of myotonic dystrophy and its pre-mRNA substrate cardiac troponin T

Myotonic dystrophy (DM) is a genetic disorder with multisystemic symptoms that is caused by expression (as RNA) of expanded repeats of CTG or CCTG in the genome. It is hypothesized that the RNA splicing factor muscleblind-like (MBNL) is sequestered to the expanded CUG or CCUG RNAs. Mislocalization of MBNL results in missplicing of a subset of pre-mRNAs that are linked to the symptoms found in DM patients. We demonstrate that MBNL can bind short structured CUG and CCUG repeats with high affinity and specificity. Only 6 base pairs are necessary for MBNL binding: two pyrimidine mismatches and four guanosine-cytosine base pairs in a stem. MBNL also has a preference for pyrimidine mismatches, but many other mismatches are tolerated with decreased affinity. We also demonstrate that MBNL binds the helical region of a stem-loop in the endogenous pre-mRNA target, the cardiac troponin T (cTNT) pre-mRNA. The stem-loop contains two mismatches and resembles both CUG and CCUG repeats. In vivo splicing results indicate that MBNL-regulated splicing is dependent upon the formation of stem-loops recognized by MBNL. These results suggest that MBNL may bind all of its RNA substrates, both normal and pathogenic, as structured stem-loops containing pyrimidine mismatches.     

Organization of the ribosomal RNA gene cluster in Aspergillus nidulans

DNA coding for ribosomal RNA in Aspergillus nidulans was found to consist of a unit 7.8 kb in size which is tandemly repeated in the genome and codes for 5.8S, 18S and 26S rRNA. The repeat unit has been cloned, and its restriction map and the location of the individual rRNA coding sequences within the unit have been established.