The loop E-loop D region of Escherichia coli 5S rRNA: the solution structure reveals an unusual loop that may be important for binding ribosomal proteins

Structure of an internal loop motif with three consecutive U•U mismatches from stem-loop 1 in the 3'-UTR of the SARS-CoV-2 genomic RNA

The single-stranded RNA genome of SARS-CoV-2 is highly structured. Numerous helical stem-loop structures interrupted by mismatch motifs are present in the functionally important 5'- and 3'-UTRs. These mismatches modulate local helical geometries and feature unusual arrays of hydrogen bonding donor and acceptor groups. However, their conformational and dynamical properties cannot be directly inferred from chemical probing and are difficult to predict theoretically. A mismatch motif (SL1-motif) consisting of three consecutive U•U base pairs is located in stem-loop 1 of the 3'-UTR. We combined NMR-spectroscopy and MD-simulations to investigate its structure and dynamics. All three U•U base pairs feature two direct hydrogen bonds and are as stable as Watson-Crick A:U base pairs. Plasmodium falciparum 25S rRNA contains a triple U•U mismatch motif (Pf-motif) differing from SL1-motif only with respect to the orientation of the two closing base pairs. Interestingly, while the geometry of the outer two U•U mismatches was identical in both motifs the preferred orientation of the central U•U mismatch was different. MD simulations and potassium ion titrations revealed that the potassium ion-binding mode to the major groove is connected to the different preferred geometries of the central base pair in the two motifs.

Multiple Non-Canonical Base-Stacking Interactions as One of the Major Determinants of RNA Tertiary Structure Organization

Stacking interactions of heterocyclic bases of ribonucleotides are one of the most important factors in the organization of RNA secondary and tertiary structure. Most of these (canonical) interactions are formed between adjacent residues in RNA polynucleotide chains. However, with the accumulation of data on the atomic tertiary structures of various RNAs and their complexes with proteins, it has become clear that nucleotide residues that are not adjacent in the polynucleotide chains and are sometimes separated in the RNA primary structure by tens or hundreds of nucleotides can interact via (non-canonical) base stacking. This paper presents an exhaustive database of such nonadjacent base-stacking elements (NA-BSEs) and their environment in the macromolecules of natural and synthetic RNAs. Analysis of these data showed that NA-BSE-forming nucleotides, on average, account for about a quarter of all nucleotides in a particular RNA and, therefore, should be considered as bona fide motifs of the RNA tertiary structure. We also classified NA-BSEs by their location in RNA macromolecules. It was shown that the structure-forming role of NA-BSEs involves compact folding of single-stranded RNA loops, transformation of double-stranded bulges into imperfect helices, and binding of RNA regions distant in the primary and secondary RNA structure.

FARFAR2: Improved De Novo Rosetta Prediction of Complex Global RNA Folds

Predicting RNA three-dimensional structures from sequence could accelerate understanding of the growing number of RNA molecules being discovered across biology. Rosetta's Fragment Assembly of RNA with Full-Atom Refinement (FARFAR) has shown promise in community-wide blind RNA-Puzzle trials, but lack of a systematic and automated benchmark has left unclear what limits FARFAR performance. Here, we benchmark FARFAR2, an algorithm integrating RNA-Puzzle-inspired innovations with updated fragment libraries and helix modeling. In 16 of 21 RNA-Puzzles revisited without experimental data or expert intervention, FARFAR2 recovers native-like structures more accurate than models submitted during the RNA-Puzzles trials. Remaining bottlenecks include conformational sampling for >80-nucleotide problems and scoring function limitations more generally. Supporting these conclusions, preregistered blind models for adenovirus VA-I RNA and five riboswitch complexes predicted native-like folds with 3- to 14 Å root-mean-square deviation accuracies. We present a FARFAR2 webserver and three large model archives (FARFAR2-Classics, FARFAR2-Motifs, and FARFAR2-Puzzles) to guide future applications and advances.

Influence of Na+ and Mg2+ ions on RNA structures studied with molecular dynamics simulations

The structure of ribonucleic acid (RNA) polymers is strongly dependent on the presence of, in particular Mg2+ cations to stabilize structural features. Only in high-resolution X-ray crystallography structures can ions be identified reliably. Here, we perform molecular dynamics simulations of 24 RNA structures with varying ion concentrations. Twelve of the structures were helical and the others complex folded. The aim of the study is to predict ion positions but also to evaluate the impact of different types of ions (Na+ or Mg2+) and the ionic strength on structural stability and variations of RNA. As a general conclusion Mg2+ is found to conserve the experimental structure better than Na+ and, where experimental ion positions are available, they can be reproduced with reasonable accuracy. If a large surplus of ions is present the added electrostatic screening makes prediction of binding-sites less reproducible. Distinct differences in ion-binding between helical and complex folded structures are found. The strength of binding (ΔG‡ for breaking RNA atom-ion interactions) is found to differ between roughly 10 and 26 kJ/mol for the different RNA atoms. Differences in stability between helical and complex folded structures and of the influence of metal ions on either are discussed.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Identification of 15 candidate structured noncoding RNA motifs in fungi by comparative genomics

BACKGROUND

With the development of rapid and inexpensive DNA sequencing, the genome sequences of more than 100 fungal species have been made available. This dataset provides an excellent resource for comparative genomics analyses, which can be used to discover genetic elements, including noncoding RNAs (ncRNAs). Bioinformatics tools similar to those used to uncover novel ncRNAs in bacteria, likewise, should be useful for searching fungal genomic sequences, and the relative ease of genetic experiments with some model fungal species could facilitate experimental validation studies.

RESULTS

We have adapted a bioinformatics pipeline for discovering bacterial ncRNAs to systematically analyze many fungal genomes. This comparative genomics pipeline integrates information on conserved RNA sequence and structural features with alternative splicing information to reveal fungal RNA motifs that are candidate regulatory domains, or that might have other possible functions. A total of 15 prominent classes of structured ncRNA candidates were identified, including variant HDV self-cleaving ribozyme representatives, atypical snoRNA candidates, and possible structured antisense RNA motifs. Candidate regulatory motifs were also found associated with genes for ribosomal proteins, S-adenosylmethionine decarboxylase (SDC), amidase, and HexA protein involved in Woronin body formation. We experimentally confirm that the variant HDV ribozymes undergo rapid self-cleavage, and we demonstrate that the SDC RNA motif reduces the expression of SAM decarboxylase by translational repression. Furthermore, we provide evidence that several other motifs discovered in this study are likely to be functional ncRNA elements.

CONCLUSIONS

Systematic screening of fungal genomes using a computational discovery pipeline has revealed the existence of a variety of novel structured ncRNAs. Genome contexts and similarities to known ncRNA motifs provide strong evidence for the biological and biochemical functions of some newly found ncRNA motifs. Although initial examinations of several motifs provide evidence for their likely functions, other motifs will require more in-depth analysis to reveal their functions.

Numerical integration methods and layout improvements in the context of dynamic RNA visualization

BACKGROUND

RNA visualization software tools have traditionally presented a static visualization of RNA molecules with limited ability for users to interact with the resulting image once it is complete. Only a few tools allowed for dynamic structures. One such tool is jViz.RNA. Currently, jViz.RNA employs a unique method for the creation of the RNA molecule layout by mapping the RNA nucleotides into vertexes in a graph, which we call the detailed graph, and then utilizes a Newtonian mechanics inspired system of forces to calculate a layout for the RNA molecule. The work presented here focuses on improvements to jViz.RNA that allow the drawing of RNA secondary structures according to common drawing conventions, as well as dramatic run-time performance improvements. This is done first by presenting an alternative method for mapping the RNA molecule into a graph, which we call the compressed graph, and then employing advanced numerical integration methods for the compressed graph representation.

RESULTS

Comparing the compressed graph and detailed graph implementations, we find that the compressed graph produces results more consistent with RNA drawing conventions. However, we also find that employing the compressed graph method requires a more sophisticated initial layout to produce visualizations that would require minimal user interference. Comparing the two numerical integration methods demonstrates the higher stability of the Backward Euler method, and its resulting ability to handle much larger time steps, a high priority feature for any software which entails user interaction.

CONCLUSION

The work in this manuscript presents the preferred use of compressed graphs to detailed ones, as well as the advantages of employing the Backward Euler method over the Forward Euler method. These improvements produce more stable as well as visually aesthetic representations of the RNA secondary structures. The results presented demonstrate that both the compressed graph representation, as well as the Backward Euler integrator, greatly enhance the run-time performance and usability. The newest iteration of jViz.RNA is available at https://jviz.cs.sfu.ca/download/download.html .

How Mg2+ ion and water network affect the stability and structure of non-Watson-Crick base pairs in E. coli loop E of 5S rRNA: a molecular dynamics and reference interaction site model (RISM) study

The non-Watson-Crick (non-WC) base pairs of Escherichia coli loop E of 5S rRNA are stabilized by Mg²⁺ ions through water-mediated interaction. It is important to know the synergic role of Mg²⁺ and the water network surrounding Mg²⁺ in stabilizing the non-WC base pairs of RNA. For this purpose, free energy change of the system is calculated using molecular dynamics (MD) simulation as Mg²⁺ is pulled from RNA, which causes disturbance of the water network. It was found that Mg²⁺ remains hexahydrated unless it is close to or far from RNA. In the pentahydrated form, Mg²⁺ interacts directly with RNA. Water network has been identified by two complimentary methods; MD followed by a density-based clustering algorithm and three-dimensional-reference interaction site model. These two methods gave similar results. Identification of water network around Mg²⁺ and non-WC base pairs gives a clue to the strong effect of water network on the stability of this RNA. Based on sequence analysis of all Eubacteria 5s rRNA, we propose that hexahydrated Mg²⁺ is an integral part of this RNA and geometry of base pairs surrounding it adjust to accommodate the [Formula: see text]. Overall the findings from this work can help in understanding the basis of the complex structure and stability of RNA with non-WC base pairs.

Structural features of a 3' splice site in influenza a

Influenza A is an RNA virus with a genome of eight negative sense segments. Segment 7 mRNA contains a 3′ splice site for alternative splicing to encode the essential M2 protein. On the basis of sequence alignment and chemical mapping experiments, the secondary structure surrounding the 3′ splice site has an internal loop, adenine bulge, and hairpin loop when it is in the hairpin conformation that exposes the 3′ splice site. We report structural features of a three-dimensional model of the hairpin derived from nuclear magnetic resonance spectra and simulated annealing with restrained molecular dynamics. Additional insight was provided by modeling based on ¹H chemical shifts. The internal loop containing the 3′ splice site has a dynamic guanosine and a stable imino (cis Watson–Crick/Watson–Crick) GA pair. The adenine bulge also appears to be dynamic with the A either stacked in the stem or forming a base triple with a Watson–Crick GC pair. The hairpin loop is a GAAA tetraloop closed by an AC pair.

Regulation of the rplY gene encoding 5S rRNA binding protein L25 in Escherichia coli and related bacteria

Ribosomal protein (r-protein) L25 is one of the three r-proteins (L25, L5, L18) that interact with 5S rRNA in eubacteria. Specific binding of L25 with a certain domain of 5S r-RNA, a so-called loop E, has been studied in detail, but information about regulation of L25 synthesis has remained totally lacking. In contrast to the rplE (L5) and rplR (L18) genes that belong to the polycistronic spc-operon and are regulated at the translation level by r-protein S8, the rplY (L25) gene forms an independent transcription unit. The main goal of this work was to study the regulation of the rplY expression in vivo. We show that the rplY promoter is down-regulated by ppGpp and its cofactor DksA in response to amino acid starvation. At the level of translation, the rplY expression is subjected to the negative feedback control. The 5'-untranslated region of the rplY mRNA comprises specific sequence/structure features, including an atypical SD-like sequence, which are highly conserved in a subset of gamma-proteobacterial families. Despite the lack of a canonical SD element, the rplY'-'lacZ single-copy reporter showed unusually high translation efficiency. Expression of the rplY gene in trans decreased the translation yield, indicating the mechanism of autogenous repression. Site-directed mutagenesis of the rplY 5' UTR revealed an important role of the conserved elements in the translation control. Thus, the rplY expression regulation represents one more example of regulatory pathways that control ribosome biogenesis in Escherichia coli and related bacteria.

De novo tertiary structure prediction using RNA123--benchmarking and application to Macugen

The present benchmarking study utilizes the RNA123 program for de novo prediction of tertiary structures of a set of 50 RNA molecules for which X-ray/NMR structures are available, based on the nucleic acid sequence only. All molecules contain a hairpin loop motif and a helical structure of canonical and non-canonical base pairs, interrupted by bulges and internal loops to various degrees. RNA molecules with double helices made up purely by canonical base pairing, and molecules containing symmetric internal loops of non-canonical base pairing are, overall, very well predicted. Structures containing bulges and asymmetric internal loops, and more complex structures containing multiple bulges and internal loops in the same molecule, result in larger deviations from their X-ray/NMR predicted structures due to higher degree of flexibility of the nucleotide bases in these regions. In a majority of the molecules included herein, the RNA123 program was, however, able to predict the tertiary structure with a heavy atom RMSD of less than 5 Å to the X-ray/NMR structure, and the models were in most cases structurally closer to the X-ray/NMR structures than models predicted by MC-Fold and MC-Sym. A set of RNA molecules containing pseudoknot tertiary structure motifs were included, but neither of the programs was able to predict the folding of the single-stranded stem onto the helix without additional structural input. The RNA123 program was then applied to predict the tertiary structure of the RNA segment of Macugen®, the first RNA aptamer approved for clinical use, and for which no tertiary structure has yet been solved. Four possible tertiary structures were predicted for this 27-nucleic-acid-long RNA molecule, which will be used in constructing a full model of the PEGylated aptamer and its interaction with the vascular endothelial growth factor target.

The influenza A segment 7 mRNA 3' splice site pseudoknot/hairpin family

The 3′ splice site of the influenza A segment 7 transcript is utilized to produce mRNA for the critical M2 ion-channel protein. In solution a 63 nt fragment that includes this region can adopt two conformations: a pseudoknot and a hairpin. In each conformation, the splice site, a binding site for the SF2/ASF exonic splicing enhancer and a polypyrimidine tract, each exists in a different structural context. The most dramatic difference occurs for the splice site. In the hairpin the splice site is between two residues that are involved in a 2 by 2 nucleotide internal loop. In the pseudoknot, however, these bases are canonically paired within one of the pseudoknotted helices. The conformational switching observed in this region has implications for the regulation of splicing of the segment 7 mRNA. A measure of stability of the structures also shows interesting trends with respect to host specificity: avian strains tend to be the most stable, followed by swine and then human.

The 3' splice site of influenza A segment 7 mRNA can exist in two conformations: a pseudoknot and a hairpin

The 3′ splice site of influenza A segment 7 is used to produce mRNA for the M2 ion-channel protein, which is critical to the formation of viable influenza virions. Native gel analysis, enzymatic/chemical structure probing, and oligonucleotide binding studies of a 63 nt fragment, containing the 3′ splice site, key residues of an SF2/ASF splicing factor binding site, and a polypyrimidine tract, provide evidence for an equilibrium between pseudoknot and hairpin structures. This equilibrium is sensitive to multivalent cations, and can be forced towards the pseudoknot by addition of 5 mM cobalt hexammine. In the two conformations, the splice site and other functional elements exist in very different structural environments. In particular, the splice site is sequestered in the middle of a double helix in the pseudoknot conformation, while in the hairpin it resides in a two-by-two nucleotide internal loop. The results suggest that segment 7 mRNA splicing can be controlled by a conformational switch that exposes or hides the splice site.

Clustering RNA structural motifs in ribosomal RNAs using secondary structural alignment

RNA structural motifs are the building blocks of the complex RNA architecture. Identification of non-coding RNA structural motifs is a critical step towards understanding of their structures and functionalities. In this article, we present a clustering approach for de novo RNA structural motif identification. We applied our approach on a data set containing 5S, 16S and 23S rRNAs and rediscovered many known motifs including GNRA tetraloop, kink-turn, C-loop, sarcin–ricin, reverse kink-turn, hook-turn, E-loop and tandem-sheared motifs, with higher accuracy than the state-of-the-art clustering method. We also identified a number of potential novel instances of GNRA tetraloop, kink-turn, sarcin–ricin and tandem-sheared motifs. More importantly, several novel structural motif families have been revealed by our clustering analysis. We identified a highly asymmetric bulge loop motif that resembles the rope sling. We also found an internal loop motif that can significantly increase the twist of the helix. Finally, we discovered a subfamily of hexaloop motif, which has significantly different geometry comparing to the currently known hexaloop motif. Our discoveries presented in this article have largely increased current knowledge of RNA structural motifs.

A quantitative measure of chirality inside nucleic acid databank

We show the capability of a chirality index (Pietropaolo et al., Proteins 2008;70:667-677) to investigate nucleic acid structures because of its high sensitivity to helical conformations. By analyzing selected structures of DNA and RNA, we have found that sequences rich in cytosine and guanine have a tendency to left-handed chirality, in contrast to regions rich in adenine or thymine which show strong negative, right-handed, chirality values. We also analyze RNA structures, where specific loops and hairpin motifs are characterized by a well-defined chirality value. We find that in nucleosome the chirality is exalted, whereas in ribosome it is reduced. Our results illustrate the sensitivity of this descriptor for nucleic acid conformations.

Solution structure of the K-turn and Specifier Loop domains from the Bacillus subtilis tyrS T-box leader RNA

In Gram-positive bacteria, the RNA transcripts of many amino acid biosynthetic and aminoacyl tRNA synthetase genes contain 5' untranslated regions, or leader RNAs, that function as riboswitches. These T-box riboswitches bind cognate tRNA molecules and regulate gene expression by a transcription attenuation mechanism. The Specifier Loop domain of the leader RNA contains nucleotides that pair with nucleotides in the tRNA anticodon loop and is flanked on one side by a kink-turn (K-turn), or GA, sequence motif. We have determined the solution NMR structure of the K-turn sequence element within the context of the Specifier Loop domain. The K-turn sequence motif has several noncanonical base pairs typical of K-turn structures but adopts an extended conformation. The Specifier Loop domain contains a loop E structural motif, and the single-strand Specifier nucleotides stack with their Watson-Crick edges displaced toward the minor groove. Mg(2+) leads to a significant bending of the helix axis at the base of the Specifier Loop domain, but does not alter the K-turn. Isothermal titration calorimetry indicates that the K-turn sequence causes a small enhancement of the interaction between the tRNA anticodon arm and the Specifier Loop domain. One possibility is that the K-turn structure is formed and stabilized when tRNA binds the T-box riboswitch and interacts with Stem I and the antiterminator helix. This motif in turn anchors the orientation of Stem I relative to the 3' half of the leader RNA, further stabilizing the tRNA-T box complex.

NMR structure and dynamics of the Specifier Loop domain from the Bacillus subtilis tyrS T box leader RNA

Gram-positive bacteria utilize a tRNA-responsive transcription antitermination mechanism, designated the T box system, to regulate expression of many amino acid biosynthetic and aminoacyl-tRNA synthetase genes. The RNA transcripts of genes controlled by this mechanism contain 5′ untranslated regions, or leader RNAs, that specifically bind cognate tRNA molecules through pairing of nucleotides in the tRNA anticodon loop with nucleotides in the Specifier Loop domain of the leader RNA. We have determined the solution structure of the Specifier Loop domain of the tyrS leader RNA from Bacillus subtilis. Fifty percent of the nucleotides in the Specifier Loop domain adopt a loop E motif. The Specifier Sequence nucleotides, which pair with the tRNA anticodon, stack with their Watson–Crick edges rotated toward the minor groove and exhibit only modest flexibility. We also show that a Specifier Loop domain mutation that impairs the function of the B. subtilis glyQS T box RNA disrupts the tyrS loop E motif. Our results suggest a mechanism for tRNA–Specifier Loop binding in which the phosphate backbone kink created by the loop E motif causes the Specifier Sequence bases to rotate toward the minor groove, which increases accessibility for pairing with bases in the anticodon loop of tRNA.

Bacterial 5S rRNA-binding proteins of the CTC family

The presence of CTC family proteins is a unique feature of bacterial cells. In the CTC family, there are true ribosomal proteins (found in ribosomes of exponentially growing cells), and at the same time there are also proteins temporarily associated with the ribosome (they are produced by the cells under stress only and incorporate into the ribosome). One feature is common for these proteins - they specifically bind to 5S rRNA. In this review, the history of investigations of the best known representatives of this family is described briefly. Structural organization of the CTC family proteins and their occurrence among known taxonomic bacterial groups are discussed. Structural features of 5S rRNA and CTC protein are described that predetermine their specific interaction. Taking into account the position of a CTC protein and its intermolecular contacts in the ribosome, a possible role of its complex with 5S rRNA in ribosome functioning is discussed.

Intergenic regions of Borrelia plasmids contain phylogenetically conserved RNA secondary structure motifs

BACKGROUND

Borrelia species are unusual in that they contain a large number of linear and circular plasmids. Many of these plasmids have long intergenic regions. These regions have many fragmented genes, repeated sequences and appear to be in a state of flux, but they may serve as reservoirs for evolutionary change and/or maintain stable motifs such as small RNA genes.

RESULTS

In an in silico study, intergenic regions of Borrelia plasmids were scanned for phylogenetically conserved stem loop structures that may represent functional units at the RNA level. Five repeat sequences were found that could fold into stable RNA-type stem loop structures, three of which are closely linked to protein genes, one of which is a member of the Borrelia lipoprotein_1 super family genes and another is the complement regulator-acquiring surface protein_1 (CRASP-1) family. Modeled secondary structures of repeat sequences display numerous base-pair compensatory changes in stem regions, including C-G-->A-U transversions when orthologous sequences are compared. Base-pair compensatory changes constitute strong evidence for phylogenetic conservation of secondary structure.

CONCLUSION

Intergenic regions of Borrelia species carry evolutionarily stable RNA secondary structure motifs. Of major interest is that some motifs are associated with protein genes that show large sequence variability. The cell may conserve these RNA motifs whereas allow a large flux in amino acid sequence, possibly to create new virulence factors but with associated RNA motifs intact.

Analyzing the flexibility of RNA structures by constraint counting

RNA requires conformational dynamics to undergo its diverse functional roles. Here, a new topological network representation of RNA structures is presented that allows analyzing RNA flexibility/rigidity based on constraint counting. The method extends the FIRST approach, which identifies flexible and rigid regions in atomic detail in a single, static, three-dimensional molecular framework. Initially, the network rigidity of a canonical A-form RNA is analyzed by counting on constraints of network elements of increasing size. These considerations demonstrate that it is the inclusion of hydrophobic contacts into the RNA topological network that is crucial for an accurate flexibility prediction. The counting also explains why a protein-based parameterization results in overly rigid RNA structures. The new network representation is then validated on a tRNA(ASP) structure and all NMR-derived ensembles of RNA structures currently available in the Protein Data Bank (with chain length >/=40). The flexibility predictions demonstrate good agreement with experimental mobility data, and the results are superior compared to predictions based on two previously used network representations. Encouragingly, this holds for flexibility predictions as well as mobility predictions obtained by constrained geometric simulations on these networks. Potential applications of the approach to analyzing the flexibility of DNA and RNA/protein complexes are discussed.

Solution structure of an rRNA substrate bound to the pseudouridylation pocket of a box H/ACA snoRNA

Base pairing between the RNA components of box H/ACA small nucleolar ribonucleoproteins (snoRNPs) and sequences in other eukaryotic RNAs target specific uridines for pseudouridylation. An RNA called HJ1 has been developed that interacts with the rRNA sequence targeted by the 5' pseudouridylation pocket of human U65 snoRNA the same way as intact U65 snoRNA. Sequences on both strands of the analog of the U65 snoRNP pseudouridylation pocket in HJ1 pair with its substrate sequence, and the resulting complex, called HJ3, is strongly stabilized by Mg(2+). The solution structure of HJ3 reveals an Omega-shaped RNA interaction motif that has not previously been described, which is likely to be common to all box H/ACA snoRNP-substrate complexes. The topology of the complex explains why the access of substrate sequences to snoRNPs is facile and how uridine selection may occur when these complexes form.

Solution structure of DNA/RNA hybrid duplex with C8-propynyl 2'-deoxyadenosine modifications: Implication of RNase H and DNA/RNA duplex interaction

Solution structures of DNA/RNA hybrid duplexes, d(GCGCA*AA*ACGCG): r(cgcguuuugcg)d(C) (designated PP57), containing two C8-propynyl 2'-deoxyadenosines (A*) and unmodified hybrid (designated U4A4) are solved. The C8-propynyl groups on 2'-deoxyadenosine perturb the local structure of the hybrid duplex, but overall the structure is similar to that of canonical DNA/RNA hybrid duplex except that Hoogsteen hydrogen bondings between A* and U result in lower thermal stability. RNase H is known to cleave RNA only in DNA/RNA hybrid duplexes. Minor groove widths of hybrid duplexes, sugar puckerings of DNA are reported to be responsible for RNase H mediated cleavage, but structural requirements for RNase H mediated cleavage still remain elusive. Despite the presence of bulky propynyl groups of PP57 in the minor groove and greater flexibility, the PP57 is an RNase H substrate. To provide an insight on the interactions between RNase H and substrates we have modeled Bacillus halodurans RNase H-PP57 complex, our NMR structure and modeling study suggest that the residue Gly(15) and Asn(16) of the loop residues between first beta sheet and second beta sheet of RNase HI of Escherichia coli might participate in substrate binding.

Facilitating RNA structure prediction with microarrays

Determining RNA secondary structure is important for understanding structure-function relationships and identifying potential drug targets. This paper reports the use of microarrays with heptamer 2'-O-methyl oligoribonucleotides to probe the secondary structure of an RNA and thereby improve the prediction of that secondary structure. When experimental constraints from hybridization results are added to a free-energy minimization algorithm, the prediction of the secondary structure of Escherichia coli 5S rRNA improves from 27 to 92% of the known canonical base pairs. Optimization of buffer conditions for hybridization and application of 2'-O-methyl-2-thiouridine to enhance binding and improve discrimination between AU and GU pairs are also described. The results suggest that probing RNA with oligonucleotide microarrays can facilitate determination of secondary structure.

Secondary structure, high variability and conserved motifs for domain III of 12S rRNA in the Arthropleona (Hexapoda; Collembola)

The robustness of different hypothetical structural models for ribosomal RNAs can be established by comparing sequences from widely different taxonomic groups. Secondary structure models of the two mitochondrial ribosomal subunits have been proposed for a large number of pterygote insects, but little information is available for the more ancient entognath hexapods. We have investigated the most frequently sequenced of the four domains of the mitochondrial small subunit rRNA in twenty-two collembolan species in the taxon Arthropleona. We provide secondary structure models of this third domain for representative species, and a consensus based on all species studied. This consensus is in partial agreement with previous models. High levels of variation, in terms of length and of primary sequence, are present in the peripheral portions of the domain. Some of the structural elements differ more among genera of springtails than among insect orders. Conversely, a limited number of structural differences compared with other taxa have been detected in the core regions.

NMR spectroscopy of RNA

NMR spectroscopy is a powerful tool for studying proteins and nucleic acids in solution. This is illustrated by the fact that nearly half of all current RNA structures were determined by using NMR techniques. Information about the structure, dynamics, and interactions with other RNA molecules, proteins, ions, and small ligands can be obtained for RNA molecules up to 100 nucleotides. This review provides insight into the resonance assignment methods that are the first and crucial step of all NMR studies, into the determination of base-pair geometry, into the examination of local and global RNA conformation, and into the detection of interaction sites of RNA. Examples of NMR investigations of RNA are given by using several different RNA molecules to illustrate the information content obtainable by NMR spectroscopy and the applicability of NMR techniques to a wide range of biologically interesting RNA molecules.

The structure of a ribosomal protein S8/spc operon mRNA complex

In bacteria, translation of all the ribosomal protein cistrons in the spc operon mRNA is repressed by the binding of the product of one of them, S8, to an internal sequence at the 5' end of the L5 cistron. The way in which the first two genes of the spc operon are regulated, retroregulation, is mechanistically distinct from translational repression by S8 of the genes from L5 onward. A 2.8 A resolution crystal structure has been obtained of Escherichia coli S8 bound to this site. Despite sequence differences, the structure of this complex is almost identical to that of the S8/helix 21 complex seen in the small ribosomal subunit, consistent with the hypothesis that autogenous regulation of ribosomal protein synthesis results from conformational similarities between mRNAs and rRNAs. S8 binding must repress the translation of its own mRNA by inhibiting the formation of a ribosomal initiation complex at the start of the L5 cistron.

Three-dimensional motifs from the SCOR, structural classification of RNA database: extruded strands, base triples, tetraloops and U-turns

Release 2.0.1 of the Structural Classification of RNA (SCOR) database, http://scor.lbl.gov, contains a classification of the internal and hairpin loops in a comprehensive collection of 497 NMR and X-ray RNA structures. This report discusses findings of the classification that have not been reported previously. The SCOR database contains multiple examples of a newly described RNA motif, the extruded helical single strand. Internal loop base triples are classified in SCOR according to their three-dimensional context. These internal loop triples contain several examples of a frequently found motif, the minor groove AGC triple. SCOR also presents the predominant and alternate conformations of hairpin loops, as shown in the most well represented tetraloops, with consensus sequences GNRA, UNCG and ANYA. The ubiquity of the GNRA hairpin turn motif is illustrated by its presence in complex internal loops.

Symmetric K+ and Mg2+ Ion-binding Sites in the 5S rRNA Loop E Inferred from Molecular Dynamics Simulations

Potassium binding to the 5S rRNA loop E motif has been studied by molecular dynamics at high (1.0 M) and low (0.2 M) concentration of added KCl in the presence and absence of Mg2+. A clear pattern of seven deep groove K+ binding sites or regions, in all cases connected with guanine N7/O6 atoms belonging to GpG, GpA, and GpU steps, was identified, indicating that the LE deep groove is significantly more ionophilic than the equivalent groove of regular RNA duplexes. Among all, two symmetry-related sites (with respect to the central G.A pair) were found to accommodate K+ ions with particularly long residence times. In a preceding molecular dynamics study by Auffinger et al. in the year 2003, these two sites were described as constituting important Mg2+ binding locations. Altogether, the data suggest that these symmetric sites correspond to the loop E main ion binding regions. Indeed, they are located in the deep groove of an important ribosomal protein binding motif associated with a fragile pattern of non-Watson-Crick pairs that has certainly to be stabilized by specific Mg2+ ions in order to be efficiently recognized by the protein. Besides, the other sites accommodate monovalent ions in a more diffuse way pointing out their lesser significance for the structure and function of this motif. Ion binding to the shallow groove and backbone atoms was generally found to be of minor importance since, at the low concentration, no well defined binding site could be characterized while high K+ concentration promoted mostly unspecific potassium binding to the RNA backbone. In addition, several K+ binding sites were located in positions equivalent to water molecules from the first hydration shell of divalent ions in simulations performed with magnesium, indicating that ion binding regions are able to accommodate both mono- and divalent ionic species. Overall, the simulations provide a more precise but, at the same time, a more intricate view of the relations of this motif with its ionic surrounding.

Solution structures of stem-loop RNAs that bind to the two N-terminal RNA-binding domains of nucleolin

Nucleolin, a multi-domain protein involved in ribosome biogenesis, has been shown to bind the consensus sequence (U/G)CCCG(A/G) in the context of a hairpin loop structure (nucleolin recognition element; NRE). Previous studies have shown that the first two RNA-binding domains in nucleolin (RBD12) are responsible for the interaction with the in vitro selected NRE (sNRE). We have previously reported the structures of nucleolin RBD12, sNRE and nucleolin RBD12-sNRE complex. A comparison of free and bound sNRE shows that the NRE loop becomes structured upon binding. From this observation, we hypothesized that the disordered hairpin loop of sNRE facilitates conformational rearrangements when the protein binds. Here, we show that nucleolin RBD12 is also sufficient for sequence- specific binding of two NRE sequences found in pre-rRNA, b1NRE and b2NRE. Structural investigations of the free NREs using NMR spectroscopy show that the b1NRE loop is conformationally heterogeneous, while the b2NRE loop is structured. The b2NRE forms a hairpin capped by a YNMG-like tetraloop. Comparison of the chemical shifts of sNRE and b2NRE in complex with nucleolin RBD12 suggests that the NRE consensus nucleotides adopt a similar conformation. These results show that a disordered NRE consensus sequence is not a prerequisite for nucleolin RBD12 binding.

The Mg2+ binding sites of the 5S rRNA loop E motif as investigated by molecular dynamics simulations

Molecular dynamics simulations have been used to investigate the binding of Mg(2+) ions to the deep groove of the eubacterial 5S rRNA loop E. The simulations suggest that long-lived and specific water-mediated interactions established between the hydrated ions and the RNA atoms lining up the binding sites contribute to the stabilization of this motif. The Mg(2+) binding specificity is modulated by two factors: (i) a required electrostatic complementarity and (ii) a structural correspondence between the hydrated ion and its binding pocket that can be estimated by its degree of dehydration and the resulting number and lifetime of the intervening water-mediated contacts. Two distinct binding modes for pentahydrated Mg(2+) ions that result in a significant freezing of the tumbling motions of the ions are described, and mechanistic details related to the stabilization of nucleic acids by divalent ions are provided.

Non-Watson-Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E

Explicit solvent and counterion molecular dynamics simulations have been carried out for a total of >80 ns on the bacterial and spinach chloroplast 5S rRNA Loop E motifs. The Loop E sequences form unique duplex architectures composed of seven consecutive non-Watson-Crick basepairs. The starting structure of spinach chloroplast Loop E was modeled using isostericity principles, and the simulations refined the geometries of the three non-Watson-Crick basepairs that differ from the consensus bacterial sequence. The deep groove of Loop E motifs provides unique sites for cation binding. Binding of Mg(2+) rigidifies Loop E and stabilizes its major groove at an intermediate width. In the absence of Mg(2+), the Loop E motifs show an unprecedented degree of inner-shell binding of monovalent cations that, in contrast to Mg(2+), penetrate into the most negative regions inside the deep groove. The spinach chloroplast Loop E shows a marked tendency to compress its deep groove compared with the bacterial consensus. Structures with a narrow deep groove essentially collapse around a string of Na(+) cations with long coordination times. The Loop E non-Watson-Crick basepairing is complemented by highly specific hydration sites ranging from water bridges to hydration pockets hosting 2 to 3 long-residing waters. The ordered hydration is intimately connected with RNA local conformational variations.

5 S rRNA: structure and interactions

5 S rRNA is an integral component of the large ribosomal subunit in all known organisms. Despite many years of intensive study, the function of 5 S rRNA in the ribosome remains unknown. Advances in the analysis of ribosome structure that have revealed the crystal structures of large ribosomal subunits and of the complete ribosome from various organisms put the results of studies on 5 S rRNA in a new perspective. This paper summarizes recently published data on the structure and function of 5 S rRNA and its interactions in complexes with proteins, within and outside the ribosome.

The solution structure of the loop E region of the 5S rRNA from spinach chloroplasts

A structure has been obtained for the loop E region of the 5S rRNA from Spinacia oleracia chloroplast ribosomes using residual dipolar coupling data as well as NOE, J coupling and chemical shift information. Even though the loop E sequence of this chloroplast 5S rRNA differs from that of Escherichia coli loop E at approximately 40% of its positions, its conformation is remarkably similar to that of E.coli loop E. Consistent with this conclusion, ribosomal protein L25 from E.coli, which binds to the loop E region of both intact E.coli 5S rRNA and to oligonucleotides containing that sequence, also binds to the chloroplast-derived oligonucleotide discussed here.

The lonepair triloop: a new motif in RNA structure

The lonepair triloop (LPTL) is an RNA structural motif that contains a single ("lone") base-pair capped by a hairpin loop containing three nucleotides. The two nucleotides immediately outside of this motif (5' and 3' to the lonepair) are not base-paired to one another, restricting the length of this helix to a single base-pair. Four examples of this motif, along with three tentative examples, were initially identified in the 16S and 23S rRNAs with covariation analysis. An evaluation of the recently determined crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits revealed the authenticity for all of these proposed interactions and identified 16 more LPTLs in the 5S, 16S and 23S rRNAs. This motif is found in the T loop in the tRNA crystal structures. The lonepairs are positioned, in nearly all examples, immediately 3' to a regular secondary structure helix and are stabilized by coaxial stacking onto this flanking helix. In all but two cases, the nucleotides in the triloop are involved in a tertiary interaction with another section of the rRNA, establishing an overall three-dimensional function for this motif. Of these 24 examples, 14 occur in multi-stem loops, seven in hairpin loops and three in internal loops. While the most common lonepair, U:A, occurs in ten of the 24 LPTLs, the remaining 14 LPTLs contain seven different base-pair types. Only a few of these lonepairs adopt the standard Watson-Crick base-pair conformations, while the majority of the base-pairs have non-standard conformations. While the general three-dimensional conformation is similar for all examples of this motif, characteristic differences lead to several subtypes present in different structural environments. At least one triloop nucleotide in 22 of the 24 LPTLs in the rRNAs and tRNAs forms a tertiary interaction with another part of the RNA. When a LPTL containing the GNR or UYR triloop sequence forms a tertiary interaction with the first (and second) triloop nucleotide, it recruits a fourth nucleotide to mediate stacking and mimic the tetraloop conformation. Approximately half of the LPTL motifs are in close association with proteins. The majority of these LPTLs are positioned at sites in rRNAs that are conserved in the three phylogenetic domains; a few of these occur in regions of the rRNA associated with ribosomal function, including the presumed site of peptidyl transferase activity in the 23S rRNA.

Structure and stability of wild-type and mutant RNA internal loops from the SL-1 domain of the HIV-1 packaging signal

The packaging signal (Psi) of the human immunodeficiency virus type 1 (HIV-1) enables encapsidation of the full-length genomic RNA against a background of a vast excess of cellular mRNAs. The core HIV-1 Psi is approximately 109 nucleotides and contains sequences critical for viral genomic dimerisation and splicing, in addition to the packaging signal. It consists of a series of stem-loops (termed SL-1 to SL-4), which can be arranged in a cloverleaf secondary structure. Using a combination of NMR spectroscopy, UV melting experiments, molecular modeling and phylogenetic analyses, we have explored the structure of two conserved internal loops proximal to the palindromic sequence of SL-1. Internal loop A, composed of six purines, forms a flexible structure that is strikingly similar to the Rev responsive element motif when bound to Rev protein. This result suggests that it may function as a protein-binding site. The absolutely conserved four-purine internal loop B is instead conformationally and thermodynamically unstable, and exhibits multiple conformations in solution. By introducing a double AGG to GGA mutation within this loop, its conformation is stabilised to form a new intra-molecular G:A:G base-triplet. The structure of the GGA mutant explains the relative instability of the wild-type loop. In a manner analogous to SL-3, we propose that conformational flexibility at this site may facilitate melting of the structure during Gag protein capture or genomic RNA dimerisation.

Structural changes of tRNA and 5S rRNA induced with magnesium and visualized with synchrotron mediated hydroxyl radical cleavage

The structure of native yeast tRNA(Phe) and wheat germ ribosomal 5S RNA induced by different magnesium ion concentrations was studied in solution with a synchrotron mediated hydroxyl radical RNA cleavage reaction. We showed that very small amounts of Mg+2 can induce significant changes in the hydroxyl radical cleavage pattern of tRNA(Phe). It also turned out that a reactivity of tRNAz(Phe) towards *OH coincides with the strong metal binding sites. Because of the Mg ions are heavily hydrated one can suggest the strong correlation of the observed nucleosides reactivity in vicinity of Mg2+ binding sites with availability of water molecules as a source of hydroxyl radical. On the other hand the structure of wheat germ 5S rRNA is less sensitive to the hydroxyl radical reaction than tRNA(Phe) although some changes are visible at 4 mM Mg ions. It is probably due to the lack of strong Mg+2 binding sites in that molecule. The reactivity of nucleotides in loops C and D of 5S rRNA is not effected, what suggests their flexibility or involvement in higher order structure formation. There is different effect of magnesium on tRNA and 5S rRNA folding. We found that nucleotides forming strong binding sites for magnesium are very sensitive to X-ray generated hydroxyl radical and can be mapped with *OH. The results show, that guanine nucleotides are preferentially hydrated. X-ray footprinting mediated hydroxyl radical RNA cleavage is a very powerful method and has been applied to studies of stable RNAs for the first time.

5S ribosomal RNA

Ribosomes have been visualized in electron micrographs in 1943 but 5S rRNA was discovered 20 years later. The next four decades witnessed big advances in our understanding of the ribosome using biochemical, genetic and low resolution structural approaches. During those times many experimental data accumulates also on 5S rRNA, but its precise function remains unknown. To understand the role of this RNA in ribosome a high-resolution structure is urgently needed. Because the ribosome is a dynamic machine, details on the interaction of 5S rRNA with proteins within entire ribosome are required. Big progress in the structural analysis of ribosome will stimulate further understanding of 5S rRNA.

The structure of helix III in Xenopus oocyte 5 S rRNA: an RNA stem containing a two-nucleotide bulge

The solution structure of an oligonucleotide containing the helix III sequence from Xenopus oocyte 5 S rRNA has been determined by NMR spectroscopy. Helix III includes two unpaired adenosine residues, flanked on either side by G:C base-pairs, that are required for binding of ribosomal protein L5. The consensus conformation of helix III in the context provided by this oligonucleotide has the two adenosine residues located in the minor groove and stacked upon the 3' flanking guanosine residue, consistent with biochemical studies of free 5 S rRNA in solution. A distinct break in stacking that occurs between the first adenosine residue of the bulge and the flanking 5' guanosine residue exposes the base of the adenosine residue in the minor groove and the base of the guanosine residue in the major groove. The major groove of the helix is widened at the site of the unpaired nucleotides and the helix is substantially bent; nonetheless, the G:C base-pairs flanking the bulge are intact. The data indicate that there may be conformational heterogeneity centered in the bulge region. The corresponding adenosine residues in the Haloarcula marismortui 50 S ribosomal subunit form a dinucleotide platform, which is quite different from the motif seen in solution. Thus, the conformation of helix III probably changes when 5 S rRNA is incorporated into the ribosome.

Zipper-like Watson-Crick base-pairs

A series of DNA heptadecamers containing the DNA analogues of RNA E-like 5'-d(GXA)/(AYG)-5' motifs (X/Y is complementary T/A, A/T, C/G, or G/C pair) were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers reveal excellent resolution in NMR spectra and exhibit many unusual nuclear Overhauser effects (NOEs) that allow for good characterization of an unusual zipper-like conformation with zipper-like Watson-Crick base-pairs; the potential canonical X.Y H-bonding is not present, and the central X/Y pairs are transformed instead into inter-strand stacks that are bracketed by sheared G.A base-pairs. Such phenomenal structural change is brought about mainly through two backbone torsional angle adjustments, i.e. delta from C2'-endo to C3'-endo for the sugar puckers of unpaired residues and gamma from gauche(+) to trans for the following 3'-adenosine residues. Such motifs are analogous to the previously studied (GGA)(2) motif presumably present in the human centromeric (TGGAA)(n) tandem repeat sequence. The novel zipper-like motifs are only 4-7 deg. C less stable than the (GGA)(2) motif, suggesting that inter-strand base stacking plays an important role in stabilizing unusual nucleic acid structures. The discovery that canonical Watson-Crick G.C or A.T hydrogen-bonded pairs can be transformed into stacking pairs greatly increases the repertoire for unusual nucleic acid structural motifs.

RNA tertiary interactions in the large ribosomal subunit: the A-minor motif

Analysis of the 2.4-A resolution crystal structure of the large ribosomal subunit from Haloarcula marismortui reveals the existence of an abundant and ubiquitous structural motif that stabilizes RNA tertiary and quaternary structures. This motif is termed the A-minor motif, because it involves the insertion of the smooth, minor groove edges of adenines into the minor groove of neighboring helices, preferentially at C-G base pairs, where they form hydrogen bonds with one or both of the 2' OHs of those pairs. A-minor motifs stabilize contacts between RNA helices, interactions between loops and helices, and the conformations of junctions and tight turns. The interactions between the 3' terminal adenine of tRNAs bound in either the A site or the P site with 23S rRNA are examples of functionally significant A-minor interactions. The A-minor motif is by far the most abundant tertiary structure interaction in the large ribosomal subunit; 186 adenines in 23S and 5S rRNA participate, 68 of which are conserved. It may prove to be the universally most important long-range interaction in large RNA structures.

Structure and function of the conserved 690 hairpin in Escherichia coli 16 S ribosomal RNA. II. NMR solution structure

The solution structure of the conserved 690 hairpin from Escherichia coli 16 S rRNA was determined by NMR spectroscopy. The 690 loop is located at the surface of the 30 S subunit in the platform region and has been implicated in interactions with P-site bound tRNA, E-site mRNA, S11 binding, IF3 binding, and in RNA-RNA interactions with the 790 loop of 16 S rRNA and domain IV of 23 S rRNA. The structure reveals a novel sheared type G690.U697 base-pair with a single hydrogen bond from the G690 amino to U697-04. G691 and A696 also form a sheared pair and U692 forms a U-turn with an H-bond to the A695 non-bridging phosphate oxygen. The sheared pairs and U-turn result in the continuous single-stranded stacking of five residues from 6693 to U697 with their Watson-Crick functional groups exposed in the minor groove. The overall fold of the 690 hairpin is similar to the anticodon loop of tRNA. The structure provides an explanation for chemical protection patterns in the loop upon interaction with tRNA, the 50 S subunit, and S11. In vivo genetic studies demonstrate the functional importance of the motifs observed in the solution structure of the 690 hairpin.

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA-RNA and RNA-protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals 'well-determined' single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit beta-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA-RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

Stability and structure of RNA duplexes containing isoguanosine and isocytidine

Isoguanosine (iG) and isocytidine (iC) differ from guanosine (G) and cytidine (C), respectively, in that the amino and carbonyl groups are transposed. The thermodynamic properties of a set of iG, iC containing RNA duplexes have been measured by UV optical melting. It is found that iG-iC replacements usually stabilize duplexes, and the stabilization per iG-iC pair is sequence-dependent. The sequence dependence can be fit to a nearest-neighbor model in which the stabilities of iG--iC pairs depend on the adjacent iG--iC or G--C pairs. For 5'-CG-3'/3'-GC-5' and 5'-GG-3'/3'-CC-5' nearest neighbors, the free energy differences upon iG-iC replacement are smaller than 0.2 kcal/mol at 37 degrees C, regardless of the number of replacements. For 5'-GC-3'/3'-CG-5', however, each iG--iC replacement adds 0.6 kcal/mol stabilizing free energy at 37 degrees C. Stacking propensities of iG and iC as unpaired nucleotides at the end of a duplex are similar to those of G and C. An NMR structure is reported for r(CiGCGiCG)(2) and found to belong to the A-form family. The structure has substantial deviations from standard A-form but is similar to published NMR and/or crystal structures for r(CGCGCG)(2) and 2'-O-methyl (CGCGCG)(2). These results provide benchmarks for theoretical calculations aimed at understanding the fundamental physical basis for the thermodynamic stabilities of nucleic acid duplexes.

Molecular dynamics simulation reveals conformational switching of water-mediated uracil-cytosine base-pairs in an RNA duplex

A 4 ns molecular dynamics simulation of an RNA duplex (r-GGACUUCGGUCC)(2 )in solution with Na+ and Cl- as counterions was performed. The X-ray structure of this duplex includes two water-mediated uracil-cytosine pairs. In contrast to the other base-pairs in the duplex the water-mediated pairs switch between different conformations. One conformation corresponds to the geometry of the water-mediated UC pairs in the duplex X-ray structure with water acting both as hydrogen-bond donor and acceptor. Another conformation is close to that of a water-mediated UC base-pair found in the X-ray structure of the 23 S rRNA sarcin/ricin domain. In this case the oxygen of the water molecule is linked to two-base donor sites. For a very short time also a direct UC base-pair and a further conformation that is similar to the one found in the RNA duplex structure but exhibits an increased H3(U)...N3(C) distance is observed. Water molecules with unusually long residence times are involved in the water-mediated conformations. These results indicate that the dynamic behaviour of the water-mediated UC base-pairs differs from that of the duplex Watson-Crick and non-canonical guanine-uracil pairs with two or three direct hydrogen bonds. The conformational variability and increased flexibility has to be taken into account when considering these base-pairs as RNA building blocks and as recognition motifs.

NMR spectroscopic evidence for Mn(2+)(Mg(2+)) binding to a precursor-tRNA microhelix near the potential RNase P cleavage site

The binding of Mg(2+)/Mn(2+) to acceptor stem microhelices as minimal models for precursor-tRNA(Gly) is demonstrated by NMR spectroscopy. From the evaluation of COSY and NOESY spectra, binding sites for Mg(2+)/Mn(2+) can be inferred. In particular, one binding site exists near the ribose moiety of nucleotide -1 at the position of cleavage by RNase P. From comparison with a variant possessing a deoxynucleotide at this position, it is concluded that the 2'-OH group of this nucleotide is indispensable for coordinating the divalent metal ion. Hence, this catalytically important metal ion is "pre-bound" to the precursor-tRNA before complexation with RNase P.