NMR evidence for dynamic secondary structure in helices II and III of the RNA of Escherichia coli

The MC-Fold and MC-Sym pipeline infers RNA structure from sequence data

The classical RNA secondary structure model considers A.U and G.C Watson-Crick as well as G.U wobble base pairs. Here we substitute it for a new one, in which sets of nucleotide cyclic motifs define RNA structures. This model allows us to unify all base pairing energetic contributions in an effective scoring function to tackle the problem of RNA folding. We show how pipelining two computer algorithms based on nucleotide cyclic motifs, MC-Fold and MC-Sym, reproduces a series of experimentally determined RNA three-dimensional structures from the sequence. This demonstrates how crucial the consideration of all base-pairing interactions is in filling the gap between sequence and structure. We use the pipeline to define rules of precursor microRNA folding in double helices, despite the presence of a number of presumed mismatches and bulges, and to propose a new model of the human immunodeficiency virus-1 -1 frame-shifting element.

Bistable secondary structures of small RNAs and their structural probing by comparative imino proton NMR spectroscopy

We investigate 25-34 nucleotide RNA sequences, that have been rationally designed to adopt two different secondary structures that are in thermodynamic equilibrium. Experimental evidence for the co-existence of the two conformers results from the NH...N 1H NMR spectra. When compared to the NH...N 1H NMR spectra of appropriate reference sequences the equilibrium position is easily quantifiable even without the assignment of the individual NH resonances. The reference sequences represent several Watson-Crick base-paired double helical segments, each encountered in either of the two conformers of the bistable target sequence. In addition, we rationalize the influence of nucleotide mutations on the equilibrium position of one of the bistable RNA sequences. The approach further allows a detailed thermodynamic analysis and the evaluation of secondary structure predictions for multistable RNAs obtained by computational methods.

Mg2+ dependence of the structure and thermodynamics of wheat germ and lupin seeds 5S rRNA

The formation and stability of structural elements in two 5S rRNA molecules from wheat germ (WG) and lupin seeds (LS) as a function of Mg2+ concentration in solution was determined using the adiabatic differential scanning microcalorimetry (DSC). The experimentally determined thermodynamic parameters are compared with calculations using thermodynamic databases used for prediction of RNA structure. The 5S rRNA molecules which show minor differences in the nucleotide sequence display very different thermal unfolding profiles (DSC profiles). Numerical deconvolution of DSC profiles provided information about structural transformations that take place in both 5S rRNA molecules. A comparative analysis of DSC data and the theoretical thermodynamic models of the structure was used to establish a relationship between the constituting transitions found in the melting profiles and the unfolding of structural domains of the 5S rRNA and stability of its particular helical elements. Increased concentrations of Mg2+ ions induces additional internal interactions stabilising 5S rRNA structures found at low Na+ concentrations. Observed conformational transitions suggest a structural model in which the extension of helical region E dominates over the postulated tertiary interaction between hairpin loops. We propose that helix E is stabilised by a sequence of non-standard pairings extending this helix by the formation of tetra loop e and an almost total reduction of loop d between helices E and D. Two hairpin structures in both 5S rRNA molecules: the extended C-C' and the extended E-E'-E" hairpins appear as the most stable elements of the structure. The cooperativity of the unfolding of helixes in these 5S rRNA molecules changes already at 2 mM Mg2+.

Structural investigation of helices II, III, and IV of B. megaterium 5S ribosomal RNA by molecular dynamics calculations

The structures of the helices II-III region and the helix IV region of B. megaterium 5S rRNA have been examined by means of energy minimization and molecular dynamics calculations. Calculated distances between neighboring hydrogen-bonded imino protons in helices II, III, and IV were between 3.5 and 4.5 A. The overall axis for the helices II-III region is warped rather than straight. Formation of additional Watson-Crick base pairs in loop B and loop C was not evident from the atomic positions calculated by molecular dynamics. Bases in loop C are well stacked, showing no significant change during dynamics. Bulge migration in helix III does not seem to be possible; the helices II-III region prefers one conformation. Helix II is more stable than helix III. Five base pairs in helix IV were sufficiently stable to establish that helix IV is terminated by a hairpin loop of three nucleotides. U87 protrudes from loop D. Structures of the helices II-III segment and the helix IV segment of B. megaterium 5S rRNA obtained by molecular dynamics were generally consistent with the solution structure inferred from high-field proton nmr spectroscopy.

Structural analysis of three prokaryotic 5S rRNA species and selected 5S rRNA--ribosomal-protein complexes by means of Pb(II)-induced hydrolysis

Lead ions have been applied to the structural analysis of 5S rRNA from Thermus thermophilus, Bacillus stearothermophilus and Escherichia coli. Based on the distribution of Pb(II)-induced cleavages, some minor modifications of the consensus secondary structure model of 5S rRNA are proposed. They include the possible base pairing between nucleotides at positions 11 and 109, as well as changes in secondary interactions within the helix B region. The 'prokaryotic arm' region is completely resistant to hydrolysis in the three RNA species, suggesting that it is a relatively stable, highly ordered structure. Hydrolysis of E. coli 5S rRNA complexed with ribosomal protein L18 shows, besides the shielding effect of the bound protein, a highly enhanced cleavage between A108 and A109. It supports the concept that the major L18-induced conformational change involves the junction of helices A, B and D.

Three-dimensional model of Escherichia coli ribosomal 5 S RNA as deduced from structure probing in solution and computer modeling

The conformation of Escherichia coli 5 S rRNA was investigated using chemical and enzymatic probes. The four bases were monitored at one of their Watson-Crick positions with dimethylsulfate (at C(N-3) and A(N-1], with a carbodiimide derivative (at G(N-1) and U(N-3] and with kethoxal (at G(N-1, N-2]. Position N-7 of purine was probed with diethylpyrocarbonate (at A(N-7] and dimethylsulfate (at G(N-7]. Double-stranded or stacked regions were tested with RNase V1 and unpaired guanine residues with RNase T1. We also used lead(II) that has a preferential affinity for interhelical and loop regions and a high sensitivity for flexible regions. Particular care was taken to use uniform conditions of salt, magnesium, pH and temperature for the different enzymatic chemical probes. Derived from these experimental data, a three dimensional model of the 5 S rRNA was built using computer modeling which integrates stereochemical constraints and phylogenetic data. The three domains of 5 S rRNA secondary structure fold into a Y-shaped structure that does not accommodate long-range tertiary interactions between domains. The three domains have distinct structural and dynamic features as revealed by the chemical reactivity and the lead(II)-induced hydrolysis: domain 2 (loop B/helix III/loop C) displays a rather weak structure and possesses dynamic properties while domain 3 (helix V/region E/helix IV/loop D) adopts a highly structured and overall helical conformation. Conserved nucleotides are not crucial for the tertiary folding but maintain an intrinsic structure in the loop regions, especially via non-canonical pairing (A.G, G.U, G.G, A.C, C.C), which can close the loops in a highly specific fashion. In particular, nucleotides in the large external loop C fold into an organized conformation leading to the formation of a five-membered loop motif. Finally, nucleotides at the hinge region of the Y-shape are involved in a precise array of hydrogen bonds based on a triple interaction between U14, G69 and G107 stabilizing the quasi-colinearity of helices II and V. The proposed tertiary model is consistent with the localization of the ribosomal protein binding sites and possesses strong analogy with the model proposed for Xenopus laevis 5 S rRNA, indicating that the Y-shape model can be generalized to all 5 S rRNAs.

Conformation of an RNA pseudoknot

The structure of the 5' GCGAUUUCUGACCGCUUUUUUGUCAG 3' RNA oligonucleotide was investigated using biochemical and chemical probes and nuclear magnetic resonance spectroscopy. Formation of a pseudoknot is indicated by the imino proton spectrum. Imino protons are observed consistent with formation of two helical stem regions; nuclear Overhauser enhancements between imino protons show that the two stem regions stack to form a continuous helix. In the stem regions, nucleotide conformations (3'-endo, anti) and internucleotide distances, derived from two-dimensional correlated, spectroscopy and two-dimensional nuclear Overhauser effect spectra, are characteristic of A-form geometry. The data suggest minor distortion in helical stacking at the junctions of stems and loops. The model of the pseudoknot is consistent with the structure originally proposed by Pleij et al.

Dynamic structure of bacterial ribosomal 5S RNA helices II and III of B. megaterium 5S RNA

A possible switch between two conformations, previously observed in an enzymatically cleaved fragment of E. coli 5S ribosomal RNA (a Gram-negative bacterium) containing helices II and III, has been examined by means of proton nuclear magnetic resonance spectroscopy (10-15 ppm) as a function of [Mg2+] and temperature for an RNase-T1 digested fragment of Bacillus megaterium 5S rRNA (a Gram-positive bacterium) containing the same helices II and III. The conformational changes induced in the fragment are not accompanied by breakage of some base-pairs and formation of others, but rather consist simply of tightening or loosening of helices with retention of existing base-pairs. Helix III is found to be more flexible than helix II. Finally, the loop conformation is conserved over a wide range of Mg2+ concentration, suggesting that the loop may serve an important role in the biological function of 5S rRNA in ribosomes.

Ribosomal 5S RNA from Xenopus laevis oocytes: conformation and interaction with transcription factor IIIA

This review describes extensive studies on 5S rRNA from X laevis oocytes combining conformational analyses in solution (using a variety of chemical and enzymatic probes), computer modeling, site-directed mutagenesis, crosslinking and TFIIIA binding. The proposed 3-dimensional model adopts a Y-shaped structure with no tertiary interactions between the different domains of the RNA. The conserved nucleotides are not crucial for the tertiary folding but they maintain an intrinsic structure in the loop regions. The model was tested by the analysis of several 5S rRNA mutants. A series of 5S RNA mutants with defined block sequence changes in regions corresponding to each of the loop regions was constructed by in vitro transcription of the mutated genes. Our results show that none of the mutations perturbs the Y-shaped structure of the RNA, although they induce conformational changes restricted to the mutated regions. The interaction of the resulting 5S rRNA mutants with TFIIIA was determined by a direct binding assay. Only the mutations in the hinge region between the 3 helical domains have a significant effect on the binding for the protein. Finally, TFIIIA was crosslinked by the use of trans-diamminedichloroplatinum (II) to a region covering the fork region. Our results show that (i) the tertiary structure does not involve long-range interactions; (ii) the intrinsic structures in loops are strictly sequence-dependent; (iii) the hinge nucleotides govern the relative orientation of the 3 helical domains; (iv) TFIIIA recognizes essentially specific features of the tertiary structure of 5S rRNA.

500-MHz proton NMR evidence for two solution structures of the common arm base-paired segment of wheat germ 5S ribosomal RNA

The base-pair protons of the common arm duplex fragment of wheat germ (Triticum aestivum) ribosomal 5S RNA have been identified and assigned by means of 500-MHz proton NMR spectroscopy. The two previously reported extra base pairs within the fragment [Li et al. (1987) Biochemistry 26, 1578-1585] are now explained by the presence of two distinct solution structures of the common arm fragment (and its corresponding base-paired segment in intact 5S rRNA). The present conclusions are supported by one- and two-dimensional proton homonuclear Overhauser enhancements in H2O and by temperature variation and Mg2+ titration of the downfield 1H NMR spectrum. The difference between the two conformers is most likely due to difference in helical tightness. Some additional amino proton resonances have also been assigned.

Physical studies of 5S RNA variants at position 66

Two variants of the 5S RNA of E. coli have been examined by imino proton NMR spectroscopy, one of them a deletion of A66 (Christiansen, J., Douthwaite, S.R., Christensen, A. and Garrett, R.A. (1985) EMBO J. 4, 1019-1024) and the other a replacement of A66 with a C (Goringer, H.U. and Wagner, R. (1986) Biol. Chem. Hoppe-Seyler 367, 769-780). Both are of interest because the role the bulged A in helix II of 5S RNA is supposed to play in interactions with ribosomal protein L18. The data show that the structural perturbations that result from these mutations are minimal, and assign the resonances of some of the imino protons around position 66. Some mutations at or near position 66 greatly reduce the L18-dependent increase in the circular dichroism of 5S RNA at 267 nm first observed by Bear and coworkers (Bear, D.G., Schleich, T., Noller, H.F. and Garrett, R.A. (1977) Nucl. Acids Res. 4, 2511-2526).

Conformation of a bulge-containing oligomer from a hot-spot sequence by NMR and energy minimization

Two-dimensional nmr data on a bulge-containing oligodeoxyribonucleotide, 5'dGATGGGCAG.dCTGACCCATC, and a regular oligomer of similar sequence, 5'dGATGGCAG.dCTGCCATC, are presented. The nonexchangeable protons are assigned from sequential nuclear Overhauser effect spectroscopy (NOESY) connectivities. The two-dimensional NOE (NOESY) and correlated (COSY) spectra of the bulge-containing oligomer are compared to those of the perfect 8-mer. Experimental proton-proton distances are determined from NOESY spectra acquired with mixing times of 100, 150, and 200 ms, using comparable distances in the B-DNA region of the molecule as a calibration. With this approach, measured distances do not depend systematically on mixing time. Energy minimization techniques are used to calculate a three-dimensional structure for the bulge-containing oligomer in agreement with the nmr data. The helix is of the B family, with the extra adenine stacked into the helix, and the helix axis is bent by 20 degrees.

Protein L18 binds primarily at the junctions of helix II and internal loops A and B in Escherichia coli 5 S RNA. Implications for 5 S RNA structure

Ribonuclease and chemical probes were used to investigate the binding sites of ribosomal protein L18 on Escherichia coli 5 S RNA using both end-labelling and reverse transcriptase procedures. The results, together with earlier data, were superimposed on a cylindrical projection of RNA double helices and most of the protection effects were found to cluster in the major groove at two sites located on one side of the RNA at the junctions of helix II with the adjoining internal loops A and B. The loop A/helix II junction was investigated using 5 S RNA mutants, produced by site-directed mutagenesis, that exhibited altered binding properties to L18. These results, together with those from a circular dichroism study of L18 complexed with the wild-type and different mutant RNAs, enabled us to assign an L18-induced conformational change to loop A. We infer that this change contributes to the co-operative binding of L5 to helix I, which may be reinforced by the binding of the very basic N-terminal peptide of L18 within the minor groove of helix I. A psoralen derivative formed a mono-addition product with U25 within loop B in the free RNA but not in the L18 complex. Moreover, the modified molecules were selected against in L18 binding experiments. Protection effects that occurred within the adjoining helix III and loop C were compatible with a tertiary interaction between loop C and loop B/helix III that could be stabilized by the L18 binding to the junction of helix II and loop B. Further support for a bipartite binding site derived from the finding that ethidium bromide molecules that are displaced from E. coli 5 S RNA by L18 intercalate both at the loop A/helix II junction and in loop B at the binding site of the psoralen derivative.

Exploration of the L18 binding site on 5S RNA by deletion mutagenesis

Several deletion variants of E. coli 5S RNA have been constructed and produced either in vivo or in vitro using T7 RNA Polymerase. Their structures and ribosomal protein L18 binding properties have been examined. All of them are similar to wild-type 5S RNA in their helix II-III regions, where L18 binds [Huber, P.W. and Wool, I.G. (1984) Proc. Natl. Acad. Sci. (USA) 81, 322-326; Douthwaite, S., Christensen, A., and Garrett, R.A. (1982) Biochemistry 21, 2313-2320.], by NMR criteria. However, none of the molecules examined that lack the helix IV-helix V stem bind L18 efficiently, even though that portion of 5S RNA is outside the L18 footprint. The L18 binding site is clearly more than a simple hairpin loop.

A comparison of the solution structures and conformational properties of the somatic and oocyte 5S rRNAs of Xenopus laevis

The secondary and tertiary structures of Xenopus oocyte and somatic 5S rRNAs were investigated using chemical and enzymatic probes. The accessibility of both RNAs towards single-strand specific nucleases (T1, T2, A and S1) and a helix-specific ribonuclease from cobra venom (RNase V1) was determined. The reactivity of nucleobase N7, N3 and N1 positions towards chemical probes was investigated under native (5 mM MgCl2, 100 mM KCl, 20 degrees C) and semi-denaturing (1 mM EDTA, 20 degrees C) conditions. Ethylnitrosourea was used to identify phosphates not reactive towards alkylation under native conditions. The results obtained confirm the presence of the five helical stems predicted by the consensus secondary structure model of 5S rRNA. The chemical reactivity data indicate that loops C and D are involved in a number of tertiary interactions, and loop E folds into an unusual secondary structure. A comparison of the data obtained for the two types of Xenopus 5S rRNA indicates that the conformations of the oocyte and somatic 5S rRNAs are very similar. However, the data obtained with nucleases under native conditions, and chemical probes under semi-denaturing conditions, reveal that helices III and IV in the somatic 5S rRNA are less stable than the same structures in oocyte 5S rRNA. Using chimeric 5S rRNAs, it was possible to demonstrate that the relative resistance of oocyte 5S rRNA to partial denaturation in 4 M urea is conferred by the five oocyte-specific nucleotide substitutions in loop B/helix III. In contrast, the superior stability of oocyte 5S rRNA in the presence of EDTA is related to a single C substitution at position 79.

Secondary structure of 5S RNA: NMR experiments on RNA molecules partially labeled with nitrogen-15

A method has been found for reassembling fragment 1 of Escherichia coli 5S RNA from mixtures containing strand III (bases 69-87) and the complex consisting of strand II (bases 89-120) and strand IV (bases 1-11). The reassembled molecule is identical with unreconstituted fragment 1. With this technique, fragment 1 molecules have been constructed 15N-labeled either in strand III or in the strand II-strand IV complex. Spectroscopic data obtained with these partially labeled molecules show that the terminal helix of 5S RNA includes the GU and GC base pairs at positions 9 and 10 which the standard model for 5S secondary structure predicts [see Delihas, N., Anderson, J., & Singhal, R. P. (1984) Prog. Nucleic Acid Res. Mol. Biol. 31, 161-190] but that these base pairs are unstable both in the fragment and in native 5S RNA. The data also assign three resonances to the helix V region of the molecule (bases 70-77 and 99-106). None of these resonances has a "normal" chemical shift even though two of them correspond to AU or GU base pairs in the standard model. The implications of these findings for our understanding of the structure of 5S RNA and its complex with ribosomal protein L25 are discussed.