New information content in RNA base pairing deduced from quantitative analysis of high-resolution structures

RNABPDB: Molecular Modeling of RNA Structure-From Base Pair Analysis in Crystals to Structure Prediction

The stable three-dimensional structure of RNA is known to play several important biochemical roles, from post-transcriptional gene regulation to enzymatic action. These structures contain double-helical regions, which often have different types of non-canonical base pairs in addition to Watson-Crick base pairs. Hence, it is important to study their structures from experimentally obtained or even predicted ones, to understand their role, or to develop a drug against the potential targets. Molecular Modeling of RNA double helices containing non-canonical base pairs is a difficult process, particularly due to the unavailability of structural features of non-Watson-Crick base pairs. Here we show a composite web-server with an associated database that allows one to generate the structure of RNA double helix containing non-canonical base pairs using consensus parameters obtained from the database. The database classification is followed by an evaluation of the central tendency of the structural parameters as well as a quantitative estimation of interaction strengths. These parameters are used to construct three-dimensional structures of double helices composed of Watson-Crick and/or non-canonical base pairs. Our benchmark study to regenerate double-helical fragments of many experimentally derived RNA structures indicate very high accuracy. This composite server is expected to be highly useful in understanding functions of various pre-miRNA by modeling structures of the molecules and estimating binding efficiency. The database can be accessed from http://hdrnas.saha.ac.in/rnabpdb .

Revisiting DNA Sequence-Dependent Deformability in High-Resolution Structures: Effects of Flanking Base Pairs on Dinucleotide Morphology and Global Chain Configuration

DNA carries more than the list of biochemical ingredients that drive the basic functions of living systems. The sequence of base pairs includes a multitude of structural and energetic signals, which determine the degree to which the long, threadlike molecule moves and how it responds to proteins and other molecules that control its processing and govern its packaging. The chemical composition of base pairs directs the spatial disposition and fluctuations of successive residues. The observed arrangements of these moieties in high-resolution protein–DNA crystal structures provide one of the best available estimates of the natural, sequence-dependent structure and deformability of the double-helical molecule. Here, we update the set of knowledge-based elastic potentials designed to describe the observed equilibrium structures and configurational fluctuations of the ten unique base-pair steps. The large number of currently available structures makes it possible to characterize the configurational preferences of the DNA base-pair steps within the context of their immediate neighbors, i.e., tetrameric context. Use of these knowledge-based potentials shows promise in accounting for known effects of sequence in long chain molecules, e.g., the degree of curvature reported in classic gel mobility studies and the recently reported sequence-dependent responses of supercoiled minicircles to nuclease cleavage.

Effects of Noncanonical Base Pairing on RNA Folding: Structural Context and Spatial Arrangements of G·A Pairs

Noncanonical base pairs play important roles in assembling the three-dimensional structures critical to the diverse functions of RNA. These associations contribute to the looped segments that intersperse the canonical double-helical elements within folded, globular RNA molecules. They stitch together various structural elements, serve as recognition elements for other molecules, and act as sites of intrinsic stiffness or deformability. This work takes advantage of new software (DSSR) designed to streamline the analysis and annotation of RNA three-dimensional structures. The multiscale structural information gathered for individual molecules, combined with the growing number of unique, well-resolved RNA structures, makes it possible to examine the collective features deeply and to uncover previously unrecognized patterns of chain organization. Here we focus on a subset of noncanonical base pairs involving guanine and adenine and the links between their modes of association, secondary structural context, and contributions to tertiary folding. The rigorous descriptions of base-pair geometry that we employ facilitate characterization of recurrent geometric motifs and the structural settings in which these arrangements occur. Moreover, the numerical parameters hint at the natural motions of the interacting bases and the pathways likely to connect different spatial forms. We draw attention to higher-order multiplexes involving two or more G·A pairs and the roles these associations appear to play in bridging different secondary structural units. The collective data reveal pairing propensities in base organization, secondary structural context, and deformability and serve as a starting point for further multiscale investigations and/or simulations of RNA folding.

How Does Mg2+ Modulate the RNA Folding Mechanism: A Case Study of the G:C W:W Trans Basepair

Reverse Watson-Crick G:C basepairs (G:C W:W Trans) occur frequently in different functional RNAs. This is one of the few basepairs whose gas-phase-optimized isolated geometry is inconsistent with the corresponding experimental geometry. Several earlier studies indicate that through post-transcriptional modification, direct protonation, or coordination with Mg²⁺, accumulation of positive charge near N7 of guanine can stabilize the experimental geometry. Interestingly, recent studies reveal significant variation in the position of putatively bound Mg²⁺. This, in conjunction with recently raised doubts regarding some of the Mg²⁺ assignments near the imino nitrogen of guanine, is suggestive of the existence of multiple Mg²⁺ binding modes for this basepair. Our detailed investigation of Mg²⁺-bound G:C W:W Trans pairs occurring in high-resolution RNA crystal structures shows that they are found in 14 different contexts, eight of which display Mg²⁺ binding at the Hoogsteen edge of guanine. Further examination of occurrences in these eight contexts led to the characterization of three different Mg²⁺ binding modes: 1) direct binding via N7 coordination, 2) direct binding via O6 coordination, and 3) binding via hydrogen-bonding interaction with the first-shell water molecules. In the crystal structures, the latter two modes are associated with a buckled and propeller-twisted geometry of the basepair. Interestingly, respective optimized geometries of these different Mg²⁺ binding modes (optimized using six different DFT functionals) are consistent with their corresponding experimental geometries. Subsequent interaction energy calculations at the MP2 level, and decomposition of its components, suggest that for G:C W:W Trans , Mg²⁺ binding can fine tune the basepair geometries without compromising with their stability. Our results, therefore, underline the importance of the mode of binding of Mg²⁺ ions in shaping RNA structure, folding and function.

Approach to the unfolding and folding dynamics of add A-riboswitch upon adenine dissociation using a coarse-grained elastic network model

Riboswitches are noncoding mRNA segments that can regulate the gene expression via altering their structures in response to specific metabolite binding. We proposed a coarse-grained Gaussian network model (GNM) to examine the unfolding and folding dynamics of adenosine deaminase (add) A-riboswitch upon the adenine dissociation, in which the RNA is modeled by a nucleotide chain with interaction networks formed by connecting adjoining atomic contacts. It was shown that the adenine binding is critical to the folding of the add A-riboswitch while the removal of the ligand can result in drastic increase of the thermodynamic fluctuations especially in the junction regions between helix domains. Under the assumption that the native contacts with the highest thermodynamic fluctuations break first, the iterative GNM simulations showed that the unfolding process of the adenine-free add A-riboswitch starts with the denature of the terminal helix stem, followed by the loops and junctions involving ligand binding pocket, and then the central helix domains. Despite the simplified coarse-grained modeling, the unfolding dynamics and pathways are shown in close agreement with the results from atomic-level MD simulations and the NMR and single-molecule force spectroscopy experiments. Overall, the study demonstrates a new avenue to investigate the binding and folding dynamics of add A-riboswitch molecule which can be readily extended for other RNA molecules.

DSSR: an integrated software tool for dissecting the spatial structure of RNA

Insight into the three-dimensional architecture of RNA is essential for understanding its cellular functions. However, even the classic transfer RNA structure contains features that are overlooked by existing bioinformatics tools. Here we present DSSR (Dissecting the Spatial Structure of RNA), an integrated and automated tool for analyzing and annotating RNA tertiary structures. The software identifies canonical and noncanonical base pairs, including those with modified nucleotides, in any tautomeric or protonation state. DSSR detects higher-order coplanar base associations, termed multiplets. It finds arrays of stacked pairs, classifies them by base-pair identity and backbone connectivity, and distinguishes a stem of covalently connected canonical pairs from a helix of stacked pairs of arbitrary type/linkage. DSSR identifies coaxial stacking of multiple stems within a single helix and lists isolated canonical pairs that lie outside of a stem. The program characterizes 'closed' loops of various types (hairpin, bulge, internal, and junction loops) and pseudoknots of arbitrary complexity. Notably, DSSR employs isolated pairs and the ends of stems, whether pseudoknotted or not, to define junction loops. This new, inclusive definition provides a novel perspective on the spatial organization of RNA. Tests on all nucleic acid structures in the Protein Data Bank confirm the efficiency and robustness of the software, and applications to representative RNA molecules illustrate its unique features. DSSR and related materials are freely available at http://x3dna.org/.

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.

Nucleic acid helix structure determination from NMR proton chemical shifts

We present a method for de novo derivation of the three-dimensional helix structure of nucleic acids using non-exchangeable proton chemical shifts as sole source of experimental restraints. The method is called chemical shift de novo structure derivation protocol employing singular value decomposition (CHEOPS) and uses iterative singular value decomposition to optimize the structure in helix parameter space. The correct performance of CHEOPS and its range of application are established via an extensive set of structure derivations using either simulated or experimental chemical shifts as input. The simulated input data are used to assess in a defined manner the effect of errors or limitations in the input data on the derived structures. We find that the RNA helix parameters can be determined with high accuracy. We finally demonstrate via three deposited RNA structures that experimental proton chemical shifts suffice to derive RNA helix structures with high precision and accuracy. CHEOPS provides, subject to further development, new directions for high-resolution NMR structure determination of nucleic acids.

Introduction to genome biology: features, processes, and structures

Genomic analyses increasingly make use of sophisticated statistical and computational approaches in investigations of genomic function and evolution. Scientists implementing and developing these approaches are often computational scientists, physicists, or mathematicians. This article aims to provide a compact overview of genome biology for these scientists. Thus, the article focuses on providing biological context to the genomic features, processes, and structures analysed by these approaches. Topics covered include (1) differences between eukaryotic and prokaryotic cells; (2) the physical structure of genomes and chromatin; (3) different categories of genomic regions, including those serving as templates for RNA and protein synthesis, regulatory regions, repetitive regions, and "architectural" or "organisational" regions, such as centromeres and telomeres; (4) the cell cycle; (5) an overview of transcription, translation, and protein structure; and (6) a glossary of relevant terms.

Complete modification maps for the cytosolic small and large subunit rRNAs of Euglena gracilis: functional and evolutionary implications of contrasting patterns between the two rRNA components

In the protist Euglena gracilis, the cytosolic small subunit (SSU) rRNA is a single, covalently continuous species typical of most eukaryotes; in contrast, the large subunit (LSU) rRNA is naturally fragmented, comprising 14 separate RNA molecules instead of the bipartite (28S+5.8S) eukaryotic LSU rRNA typically seen. We present extensively revised secondary structure models of the E. gracilis SSU and LSU rRNAs and have mapped the positions of all of the modified nucleosides in these rRNAs (88 in SSU rRNA and 262 in LSU rRNA, with only 3 LSU rRNA modifications incompletely characterized). The relative proportions of ribose-methylated nucleosides and pseudouridine (∼60% and ∼35%, respectively) are closely similar in the two rRNAs; however, whereas the Euglena SSU rRNA has about the same absolute number of modifications as its human counterpart, the Euglena LSU rRNA has twice as many modifications as the corresponding human LSU rRNA. The increased levels of rRNA fragmentation and modification in E. gracilis LSU rRNA are correlated with a 3-fold increase in the level of mispairing in helical regions compared to the human LSU rRNA. In contrast, no comparable increase in mispairing is seen in helical regions of the SSU rRNA compared to its homologs in other eukaryotes. In view of the reported effects of both ribose-methylated nucleoside and pseudouridine residues on RNA structure, these correlations lead us to suggest that increased modification in the LSU rRNA may play a role in stabilizing a 'looser' structure promoted by elevated helical mispairing and a high degree of fragmentation.

All-atom knowledge-based potential for RNA structure prediction and assessment

MOTIVATION

Over the recent years, the vision that RNA simply serves as information transfer molecule has dramatically changed. The study of the sequence/structure/function relationships in RNA is becoming more important. As a direct consequence, the total number of experimentally solved RNA structures has dramatically increased and new computer tools for predicting RNA structure from sequence are rapidly emerging. Therefore, new and accurate methods for assessing the accuracy of RNA structure models are clearly needed.

RESULTS

Here, we introduce an all-atom knowledge-based potential for the assessment of RNA three-dimensional (3D) structures. We have benchmarked our new potential, called Ribonucleic Acids Statistical Potential (RASP), with two different decoy datasets composed of near-native RNA structures. In one of the benchmark sets, RASP was able to rank the closest model to the X-ray structure as the best and within the top 10 models for ∼93 and ∼95% of decoys, respectively. The average correlation coefficient between model accuracy, calculated as the root mean square deviation and global distance test-total score (GDT-TS) measures of C3' atoms, and the RASP score was 0.85 and 0.89, respectively. Based on a recently released benchmark dataset that contains hundreds of 3D models for 32 RNA motifs with non-canonical base pairs, RASP scoring function compared favorably to ROSETTA FARFAR force field in the selection of accurate models. Finally, using the self-splicing group I intron and the stem-loop IIIc from hepatitis C virus internal ribosome entry site as test cases, we show that RASP is able to discriminate between known structure-destabilizing mutations and compensatory mutations.

AVAILABILITY

RASP can be readily applied to assess all-atom or coarse-grained RNA structures and thus should be of interest to both developers and end-users of RNA structure prediction methods. The computer software and knowledge-based potentials are freely available at http://melolab.org/supmat.html.

CONTACT

fmelo@bio.puc.cl; mmarti@cipf.es

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Site-directed mutants of 16S rRNA reveal important RNA domains for KsgA function and 30S subunit assembly

KsgA is an rRNA methyltransferase important to the process of small subunit biogenesis in bacteria. It is ubiquitously found in all life including archaea and eukarya, where the enzyme is referred to as Dim1. Despite the emergence of considerable data addressing KsgA function over the last several years, details pertaining to RNA recognition are limited, in part because the most accessible substrate for in vitro studies of KsgA is the 900000 Da 30S ribosomal subunit. To overcome challenges imposed by size and complexity, we adapted recently reported techniques to construct in vivo assembled mutant 30S subunits suitable for use in in vitro methyltransferase assays. Using this approach, numerous 16S rRNA mutants were constructed and tested. Our observations indicate that the 790 loop of helix 24 plays an important role in overall catalysis by KsgA. Moreover, the length of helix 45 also is important to catalysis. In both cases loss of catalytic function occurred without an increase in the production of N(6)-methyladenosine, a likely indication that there was no critical reduction in binding strength. Both sets of observations support a "proximity" mechanism of KsgA function. We also report that several of the mutants constructed failed to assemble properly into 30S subunits, while some others did so with reduced efficiency. Therefore, the same technique of generating mutant 30S subunits can be used to study ribosome biogenesis on the whole.

Unusual target site disruption by the rare-cutting HNH restriction endonuclease PacI

The crystal structure of the rare-cutting HNH restriction endonuclease PacI in complex with its eight-base-pair target recognition sequence 5'-TTAATTAA-3' has been determined to 1.9 A resolution. The enzyme forms an extended homodimer, with each subunit containing two zinc-bound motifs surrounding a betabetaalpha-metal catalytic site. The latter is unusual in that a tyrosine residue likely initiates strand cleavage. PacI dramatically distorts its target sequence from Watson-Crick duplex DNA base pairing, with every base separated from its original partner. Two bases on each strand are unpaired, four are engaged in noncanonical A:A and T:T base pairs, and the remaining two bases are matched with new Watson-Crick partners. This represents a highly unusual DNA binding mechanism for a restriction endonuclease, and implies that initial recognition of the target site might involve significantly different contacts from those visualized in the DNA-bound cocrystal structures.

Role of unsatisfied hydrogen bond acceptors in RNA energetics and specificity

RNA plays essential roles in much of biology. These functions are dictated by structures mediated by hydrogen bonding, stacking, electrostatics, and steric interactions. Roles of unsatisfied hydrogen bond functionalities in these structures are less well understood. Herein, we evaluated the energetic contributions of unsatisfied hydrogen bonding groups by placing chemically modified substituents in select internal positions in RNA helices and conducting thermodynamic studies. We find that unsatisfied carbonyl groups make exceptional contributions to structure formation (approximately 3 kcal/mol in free energy), most likely due to a combination of strain and dehydration effects. Thus, unsatisfied hydrogen bonding groups are likely key determinants in the folding energetics and specificity of many RNA and DNA molecules and may be especially important in tertiary structure interactions.

On the role of Hoogsteen:Hoogsteen interactions in RNA: ab initio investigations of structures and energies

We use a combination of database analysis and quantum chemical studies to investigate the role of cis and trans Hoogsteen:Hoogsteen (H:H) base pairs and associated higher-order structures in RNA. We add three new examples to the list of previously identified base-pair combinations belonging to these families and, in addition to contextual classification and characterization of their structural and energetic features, we compare their interbase interaction energies and propensities toward participation in triplets and quartets. We find that some base pairs, which are nonplanar in their isolated minimum energy geometries, attain planarity and stability upon triplet formation. A:A H:H trans is the most frequent H:H combination in RNA structures. This base pair occurs at many distinct positions in known rRNA structures, where it helps in the interaction of ribosomal domains in the 50S subunit. It is also present as a part of tertiary interaction in tRNA structures. Although quantum chemical studies suggest an intrinsically nonplanar geometry for this base pair in isolated form, it has the tendency to attain planar geometry in RNA crystal structures by forming higher-order tertiary interactions or in the presence of additional base-phosphate interactions. The tendency of this base pair to form such additional interactions may be helpful in bringing together different segments of RNA, thus making it suitable for the role of facilitator for RNA folding. This also explains the high occurrence frequency of this base pair among all H:H interactions.

Diversity of 16S rRNA genes within individual prokaryotic genomes

Analysis of intragenomic variation of 16S rRNA genes is a unique approach to examining the concept of ribosomal constraints on rRNA genes; the degree of variation is an important parameter to consider for estimation of the diversity of a complex microbiome in the recently initiated Human Microbiome Project (http://nihroadmap.nih.gov/hmp). The current GenBank database has a collection of 883 prokaryotic genomes representing 568 unique species, of which 425 species contained 2 to 15 copies of 16S rRNA genes per genome (2.22 +/- 0.81). Sequence diversity among the 16S rRNA genes in a genome was found in 235 species (from 0.06% to 20.38%; 0.55% +/- 1.46%). Compared with the 16S rRNA-based threshold for operational definition of species (1 to 1.3% diversity), the diversity was borderline (between 1% and 1.3%) in 10 species and >1.3% in 14 species. The diversified 16S rRNA genes in Haloarcula marismortui (diversity, 5.63%) and Thermoanaerobacter tengcongensis (6.70%) were highly conserved at the 2 degrees structure level, while the diversified gene in B. afzelii (20.38%) appears to be a pseudogene. The diversified genes in the remaining 21 species were also conserved, except for a truncated 16S rRNA gene in "Candidatus Protochlamydia amoebophila." Thus, this survey of intragenomic diversity of 16S rRNA genes provides strong evidence supporting the theory of ribosomal constraint. Taxonomic classification using the 16S rRNA-based operational threshold could misclassify a number of species into more than one species, leading to an overestimation of the diversity of a complex microbiome. This phenomenon is especially seen in 7 bacterial species associated with the human microbiome or diseases.

The RNA backbone plays a crucial role in mediating the intrinsic stability of the GpU dinucleotide platform and the GpUpA/GpA miniduplex

The side-by-side interactions of nucleobases contribute to the organization of RNA, forming the planar building blocks of helices and mediating chain folding. Dinucleotide platforms, formed by side-by-side pairing of adjacent bases, frequently anchor helices against loops. Surprisingly, GpU steps account for over half of the dinucleotide platforms observed in RNA-containing structures. Why GpU should stand out from other dinucleotides in this respect is not clear from the single well-characterized H-bond found between the guanine N2 and the uracil O4 groups. Here, we describe how an RNA-specific H-bond between O2′(G) and O2P(U) adds to the stability of the GpU platform. Moreover, we show how this pair of oxygen atoms forms an out-of-plane backbone ‘edge’ that is specifically recognized by a non-adjacent guanine in over 90% of the cases, leading to the formation of an asymmetric miniduplex consisting of ‘complementary’ GpUpA and GpA subunits. Together, these five nucleotides constitute the conserved core of the well-known loop-E motif. The backbone-mediated intrinsic stabilities of the GpU dinucleotide platform and the GpUpA/GpA miniduplex plausibly underlie observed evolutionary constraints on base identity. We propose that they may also provide a reason for the extreme conservation of GpU observed at most 5′-splice sites.

RNA structural motifs that entail hydrogen bonds involving sugar-phosphate backbone atoms of RNA

The growing number of high-resolution crystal structures of large RNA molecules provides much information for understanding the principles of structural organization of these complex molecules. Several in-depth analyses of nucleobase-centered RNA structural motifs and backbone conformations have been published based on this information, including a systematic classification of base pairs by Leontis and Westhof. However, hydrogen bonds involving sugar-phosphate backbone atoms of RNA have not been analyzed systematically until recently, although such hydrogen bonds appear to be common both in local and tertiary interactions. Here we review some backbone structural motifs discussed in the literature and analyze a set of eight high-resolution multi-domain RNA structures. The analyzed RNAs are highly structured: among 5372 nucleotides in this set, 89% are involved in at least one "long-range" RNA-RNA hydrogen bond, i.e., hydrogen bonds between atoms in the same residue or sequential residues are ignored. These long-range hydrogen bonds frequently use backbone atoms as hydrogen bond acceptors, i.e., OP1, OP2, O2', O3', O4', or O5', or as a donor (2'OH). A surprisingly large number of such hydrogen bonds are found, considering that neither single-stranded nor double-stranded regions will contain such hydrogen bonds unless additional interactions with other residues exist. Among 8327 long-range hydrogen bonds found in this set of structures, 2811, or about one-third, are hydrogen bonds entailing RNA backbone atoms; they involve 39% of all nucleotides in the structures. The majority of them (2111) are hydrogen bonds entailing ribose hydroxyl groups, which can be used either as a donor or an acceptor; they constitute 25% of all hydrogen bonds and involve 31% of all nucleotides. The phosphate oxygens OP1 or OP2 are used as hydrogen bond acceptors in 12% of all nucleotides, and the ribose ring oxygen O4' and phosphodiester oxygens O3' and O5' are used in 4%, 4%, and 1% of all nucleotides, respectively. Distributions of geometric parameters and some examples of such hydrogen bonds are presented in this report. A novel motif involving backbone hydrogen bonds, the ribose-phosphate zipper, is also identified.