Higher order structure in the 3'-minor domain of small subunit ribosomal RNAs from a gram negative bacterium, a gram positive bacterium and a eukaryote

An experimental approach was used to determine and compare the highest order structure within the 150 to 200 nucleotides at the 3'-ends of the RNAs from the small ribosomal subunits of Escherichia coli, Bacillus stearothermophilus and Saccharomyces cerevisiae. Chemical reagents were employed to establish the degree of stacking and/or accessibility of each adenosine, guanosine and cytidine. The double helices were probed with a cobra venom ribonuclease from Naja naja oxiana, and the relatively unstructured and accessible sequences were localized with the single strand-specific ribonucleases A, T1, T2 and S1. The data enabled the various minimal secondary structural models, proposed for the 3'-regions of the E. coli and S. cerevisiae RNAs, to be critically examined, and to demonstrate that the main common features of these models are correct. The results also reveal the presence and position of additional higher order structure in the renatured free RNA. It can be concluded that a high level of conservation of higher order structure has occurred during the evolution of the gram negative and gram positive eubacteria and the eukaryote in both the double helical regions and the "unstructured" regions. Several unusual structural features were detected. Multiple G X A pairings in two of the putative helices, which are compatible with phylogenetic sequence comparisons, are strongly supported by the occurrence of cobra venom ribonuclease cuts adjacent to, and in one case between, these pairings. Evidence is also provided for the stacking of an A X A pair within a double helix of the yeast RNA. Other special structural features include adenosines bulged out from double helices; such nucleotides, which are hyper-reactive, have been implicated in protein recognition in 5 S ribosomal RNA. The 3'-terminal regions of the RNAs are particularly important for the functioning of the ribosome. They are involved in mRNA, tRNA and ribosomal factor binding. The results reveal that while the functionally important RNA sequences tend to be conserved, they are not always accessible in the free RNA; the pyrimidine-rich "Shine and Dalgarno" sequence, for example, which is involved in mRNA recognition, occurs in a double helix in both eubacterial RNAs.     

Different conformational families of pyrimidine.purine.pyrimidine triple helices depending on backbone composition.

Different helical conformations of DNA (D), RNA (R), and DNA.RNA (DR) hybrid double and triple helices have been detected using affinity cleavage analysis. Synthetic methods were developed to attach EDTA.Fe to a single nucleotide on RNA as well as DNA oligonucleotides. Cleavage patterns generated by a localized diffusible oxidant in the major groove on the pyrimidine strand of four purine.pyrimidine double helices consisting of all DNA, all RNA, and the corresponding hybrids reveal that the relative cleavage intensity shifts to the 5' end of the purine strand increasingly in the order: DD < DR < RD < RR. These results are consistent with models derived from structural studies. In six pyrimidine.purine.pyrimidine triple helices, the altered cleavage patterns of the Watson-Crick pyrimidine strands reveal at least two conformational families: (i) D + DD, R + DD, D + DR, and R + DR and (ii) R + RD and R + RR.     

The fastest global events in RNA folding: electrostatic relaxation and tertiary collapse of the Tetrahymena ribozyme

Large RNAs can collapse into compact conformations well before the stable formation of the tertiary contacts that define their final folds. This study identifies likely physical mechanisms driving these early compaction events in RNA folding. We have employed time-resolved small-angle X-ray scattering to monitor the fastest global shape changes of the Tetrahymena ribozyme under different ionic conditions and with RNA mutations that remove long-range tertiary contacts. A partial collapse in each of the folding time-courses occurs within tens of milliseconds with either monovalent or divalent cations. Combined with comparison to predictions from structural models, this observation suggests a relaxation of the RNA to a more compact but denatured conformational ensemble in response to enhanced electrostatic screening at higher ionic concentrations. Further, the results provide evidence against counterion-correlation-mediated attraction between RNA double helices, a recently proposed model for early collapse. A previous study revealed a second 100 ms phase of collapse to a globular state. Surprisingly, we find that progression to this second early folding intermediate requires RNA sequence motifs that eventually mediate native long-range tertiary interactions, even though these regions of the RNA were observed to be solvent-accessible in previous footprinting studies under similar conditions. These results help delineate an analogy between the early conformational changes in RNA folding and the "burst phase" changes and molten globule formation in protein folding.     

Common structural features of different viroids: serial arrangement of double helical sections and internal loops.

The thermodynamic parameters of five different highly purified viroid "species" were determined by applying UV-absorption melting analysis and temperature jump methods. Their thermal denaturation proved to be a highly cooperative process with midpoint-temperatures (Tm) between 48.5 and 51 degrees C in 0.01 M sodium cacodylate, 1 mM EDTA, pH 6.8. The values of the apparent reaction enthalpies of the different viroid species range between 3,140 and 3,770 kJ/mol. Although the cooperativity is as high as found in homogeneous RNA double helices the Tm-value of viroid melting is more than 30 degrees C lower than in the homogeneous RNA. In order to explain this deviation, melting curves were simulated for different models of the secondary structure of viroids using literature values of the thermodynamic parameters of nucleic acids. Our calculations show that the following refinement of our earlier model is in complete accordance with the experimental data: In their native conformation viroids exist as an extended rodlike structure characterized by a series of double helical sections and internal loops. In the different viroid species 250-300 nucleotides out of total 350 nucleotides are needed to interprete the thermodynamic behaviour.     

Effect of single-residue bulges on RNA double-helical structures: crystallographic database analysis and molecular dynamics simulation studies

Asymmetric bulge loop motifs are widely dispersed in all types of functional RNAs. They are frequently occurring structural motifs in folded RNA structures and appear commonly in pre-microRNA and ribosomes, where they are involved in specific RNA–RNA and RNA–protein interactions. It is therefore necessary to understand such motifs from a structural point of view. We analyzed all available RNA structures and identified quite a few fragments of double helices that contain bulges. We found that these discontinuities often introduce kinks into the double helices, which also affects the stacking overlap between the base pairs across the irregularity. In order to understand the influence of these bulges on stability and flexibility, we carried out molecular dynamics simulations of three different single-residue bulge-containing RNA helices using the CHARMM36 force field. The structural variability at the junctions of RNA bulges is expected to differ from that in continuous double-helical stretches. The structural features of the junction region were observed to vary noticeably depending on the orientation of the bulge residue. When the base of the bulge residue is looped out, the RNA stretch behaves like a standard long A-form RNA double helix, whereas the entire RNA behaves differently when the base of the bulge residue is intercalated between base pairs inside the RNA stem. Such single-base intercalation was found to introduce a permanent kink into the composite double helix, which could be a recognition element for Dicer during the maturation of miRNA.     

The template RNAs of RNA polymerases can have compact secondary structure, formed by long double helices with partial violations of the complementarity

A new method of contextual analysis was used to search the long non-random inverted repeats and the complementary palindromes in the genes of E. coli and T7 RNA polymerases. These genes were found to contain from 25% to 50% of all the nucleotides involved in such helices. The 5' -and 3' -ends of mRNA can be protected by neighbouring double helices from the nuclease attack. Some double helices are competing and very similar to the attenuator of E. coli trp-operon.     

The relationship between internal chain length of amylopectin and crystallinity in starch

Molecular models of amylopectin were created and investigated by computer simulation. First, single and double helices of various lengths were constructed. The 1 → 6 branching in double and single helices of amylopectin was studied. Subunits of single helices, double helices, and branch points were used as building blocks of larger systems. The possible makeup of amylopectin unit clusters was investigated via a series of models, including single–single, double–single, and double–double helix systems. The lengths of the single helix section that linked two branch points (internal chains) was systematically varied between values of 0–10 glucose residues. It was found that certain internal chain lengths lead to parallel double helices. Thus, it was postulated that the length of internal chains may determine the degree of local crystallinity. Furthermore, it was noted that some of the low‐energy arrangement of double helices could be superimposed on either the two adjacent and nonadjacent double helices of crystalline A and B starch polymorphs. In other cases, the distance between the double helices is so large that it may in fact be a model for branching between two amylopectin crystals or unit clusters. Results obtained through this work were corroborated, where possible, with information available from crystallographic, branching, and enzymatic studies. © 1999 John Wiley & Sons, Inc. Biopoly 50: 381–390, 1999     

DbpA is a region‐specific RNA helicase

DbpA is a DEAD‐box RNA helicase implicated in RNA structural rearrangements in the peptidyl transferase center. DbpA contains an RNA binding domain, responsible for tight binding of DbpA to hairpin 92 of 23S ribosomal RNA, and a RecA‐like catalytic core responsible for double‐helix unwinding. It is not known if DbpA unwinds only the RNA helices that are part of a specific RNA structure, or if DbpA unwinds any RNA helices within the catalytic core's grasp. In other words, it is not known if DbpA is a site‐specific enzyme or region‐specific enzyme. In this study, we used protein and RNA engineering to investigate if DbpA is a region‐specific or a site‐specific enzyme. Our data suggest that DbpA is a region‐specific enzyme. This conclusion has an important implication for the physiological role of DbpA. It suggests that during ribosome assembly, DbpA could bind with its C‐terminal RNA binding domain to hairpin 92, while its catalytic core may unwind any double‐helices in its vicinity. The only requirement for a double‐helix to serve as a DbpA substrate is for the double‐helix to be positioned within the catalytic core's grasp.     

UDP-glucose dehydrogenase: substrate binding stoichiometry and affinity

Precise structural parameters of polyribonucleotides single stranded helices are determined as well as those of double stranded helices of poly 2′-O-methyl A and of poly A at neutral and acid pH. Infrared linear dichroism investigations indicate the similarity of the conformation of the sugar-phosphate backbone of these single and double stranded helices. The angles of the phosphate group for single stranded helix at neutral pH is found to be oriented at 48° for the 02P02 bisector and at about 65° for the 02–03 line to the helix axis. Similar values were found for double stranded poly A helix at acid pH. These structural parameters obtained for the first time on single stranded polynucleotide helices are proposed to be valid for other similar helical chains such as poly A segments of nuclear or messenger RNA and single stranded CCA acceptor end of transfer RNA.     

Parallel nucleic acid helices with hoogsteen base pairing: Symmetry and structure

Molecular structures for parallel DNA and RNA double helices with Hoogsteen pairing are proposed for the first time. The DNA helices have sugars in the C2′-endo region and the phosphodiester conformations are (trans, gauche^?), and the RNA helices have sugars in the C3′-endo region and the phosphodiester conformations are (gauche^?, gauche^?). A pseudorotational symmetry relates the two parallel strands of DNA helices and a screw symmetry relates the two strands of RNA helices, which have an associated tilt of the The conformational space of parallel helices with Hoogsteen base pairing, unlike the Watson-Crick duplex, is highly restricted due to the unique positioning of the symmetry axis in the former case. The features of the parallel double helix with Hoogsteen pairing are compared with the Watson-Crick duplex and the corresponding triple helix. © 1994 John Wiley & Sons, Inc.     

Studies on the secondary structure of ribosomal ribonucleic acid components of rabbit reticulocytes

The conformation in solution of fractionated 30 S and 19 S ribosomal RNA from rabbit reticulocytes has been studied by optical rotatory dispersion, analysis of thermal melting profiles and their derivatives, and spectrophotometric acid-base titration. From a consideration of the limitations of these methods, it has been possible to set limiting values on the degree of base-pairing and the lengths of the double helices: between 60 and 80% of the bases in 19 S and 30 S RNA are estimated to be paired. The paired segments are not shorter than 4 base pairs, and evidence from other sources is available which indicates that they are not longer than 8–16 base pairs. The spread of helix lengths is greater in the 30 S than in 19 S RNA; and other differences are noted. Several distinct populations of double helices, differing in their thermal stability, are present. Estimates are presented from spectrophotometric and titration data for the base compositions of the paired and unpaired regions.     

Prediction of three-dimensional structure of Escherichia coli ribosomal RNA

A model for the tertiary structure of 23S, 16S and 5S ribosomal RNA molecules interacting with three tRNA molecules is presented using the secondary structure models common to E. coli, Z. mays chloroplast, and mammalian mitochondria. This ribosomal RNA model is represented by phosphorus atoms which are separated by 5.9 A in the standard A-form double helix conformation. The accumulated proximity data summarized in Table 1 were used to deduce the most reasonable assembly of helices separated from each other by at least 6.2 A. Straight-line approximation for single strands was adopted to describe the maximum allowed distance between helices. The model of a ribosome binding three tRNA molecules by Nierhaus (1984), the stereochemical model of codon-anticodon interaction by Sundaralingam et al. (1975) and the ribosomal transpeptidation model, forming an alpha-helical nascent polypeptide, by Lim & Spirin (1986), were incorporated in this model. The distribution of chemically modified nucleotides, cross-linked sites, invariant and missing regions in mammalian mitochondrial rRNAs are indicated on the model.     

Protein binding sites on Escherichia coli 16S ribosomal RNA; RNA regions that are protected by proteins S7, S9 and S19, and by proteins S8, S15 and S17.

Selected groups of isolated 14C-labelled proteins from E. coli 30S ribosomal subunits were reconstituted with 32P-labelled 16S RNA, and the reconstituted complexes were partially digested with ribonuclease A. RNA fragments protected by the proteins were separated by gel electrophoresis and subjected to sequence analysis. Complexes containing proteins S7 and S19 protected an RNA region comprising helices 29 to 32, part of helix 41, and helices 42 and 43 of the 16S RNA secondary structure. Addition of protein S9 had no effect. When compared with previous data for proteins S7, S9, S14 and S19, these results suggest that S14 interacts with helix 33, and that S9 and S14 together interact with the loop-end of helix 41. Complexes containing proteins S8, S15 and S17 protected helices 7 to 10 as well as the "S8-S15 binding site" (helices 20, 22 and parts of helices 21 and 23). When protein S15 was omitted, S8 and S18 showed protection of part of helix 44 in addition to the latter regions. The results are discussed in terms of our model for the detailed arrangement of proteins and RNA in the 30S subunit.     

Free energy of imperfect nucleic acid helices. 3. Small internal loops resulting from mismatches

Physical studies of enzymioally synthesized oligoribonucleotides of defined sequence are used to evaluate quantitatively the destabilizing influence of mismatched bases in a double helix. The series (A-)₄G(-C)_n(-U)₄, N = 1 to 6, exist as imperfect dimer helices when N is equal to or less than 4, and as monomolecular hairpin helices when N is 5 and 6. Internal loops become progressively more destabilizing as their size increases from 2 to 4 to 6 nucleotides resulting from 1, 2 and 3 consecutive mismatched base pairs. However, the stability of a helix will generally be greater if a given number of mismatched pairs occur consecutively rather than in isolation from one another.These data may be used for improved calculations of stability of RNA secondary structure, to estimate the frequency of structural fluctuations in a double helix and to assess the stability of modified polynucleotide helices. An unmodified double helix of one million randomly arranged base pairs should contain on the time average approximately 10 G.C and 500 A.U pairs in non-hydrogen bonded, unstacked conformations at 25 °C. Our estimate of the effect of mismatching on T_m values of high polymers is less precise because of the long temperature extrapolation required. However, we estimate that DNA or RNA treated with mutagens which interrupt up to 20% of the nucleotide pairs will show a drop of about 1.2 deg. C in melting temperature with each unit per cent of modification.     

Defining topological equivalences in macromolecules.

We describe a completely automated and objective method for defining topological equivalents in macromolecules. The method is based on well established techniques for identifying topologically and topographically equivalent atoms in small molecules and has been used in structural alignment of proteins and RNA molecules, and to extract fragments of molecules (protein secondary structures and RNA and DNA double helices) from structural databases consistent with some specified template structure.     

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.     

Phylogenetic relationships among Myrceugenia, Blepharocalyx, and Luma (Myrtaceae) based on paired-sites models and the secondary structures of ITS and ETS sequences

To establish relationships among the genera Blepharocalyx, Luma, and Myrceugenia, the phylogeny of tribe Myrteae was reconstructed based on secondary structures of the sequences of ITS (ITS1-5.8S-ITS2) and ETS regions from 93 taxa belonging to 29 genera. The DNA sequences were aligned according to their secondary structure and analysed with Bayesian inference (BI) using base substitution models for RNA, which can differentiate between the regions consisting of base pairs (helices) and unpaired bases (loops). The best-fit models were RNA7C for ITS and RNA7B for ETS, which differ in how they predict double substitutions and symmetry of the frequency of base pairs. The phylogenetic hypothesis was compared with results obtained using classical 4-state models, the results of maximum parsimony, and using sequence alignment without considering the secondary structure. Analyses show low support along the backbone of the consensus trees, those using substitution models of paired sites showing more highly supported derived clades than substitution models of independent sites. All tests were consistent and differed in the positioning of certain groups, but they also showed that Blepharocalyx, Luma, and Myrceugenia arise from three independent lineages that are not closely related. The substitution model for RNA shows that Luma is a monophyletic group and sister to the other genera of Myrteae, except for Myrtus communis. Myrceugenia appears as a monophyletic group, except for M. fernandeziana, which is related to Blepharocalyx, a genus that appears to be biphyletic.     

Molecular Modelling of Helix Stability in Carrageenans

Molecular models of disaccharides, and single and double helices up to eight monomers in length have been constructed of the two types of glycosidic linkage in the carrageenan chain. These links are a galactose to anhydrogalactose link (GA link), and an anhydrogalactose to galactose link (AG link). These models are also based on 3-carrageenan, which contains a 4-sulphate galactose ring. The effects of the sulphate groups on the conformation of the helices may be seen by the angles of $ϕ$ϕ and N explored during the simulations by the AG and GA linkages. It has been observed that the molecule can explore a greater area of conformational space about the GA link than the AG link. This could be due to steric hindrance caused by the bulky sulphate group near the AG link. The sulphate group is further away from the GA link than from the AG link, and this may provide a possible explanation for the relatively unhindered movement about the GA link compared to the AG link. The results have also shown that the conformational space for the AG linkages, as well as the GA linkages vary between different lengths of the polysaccharide chain. Single helix models show little stability in molecular dynamics simulation, whereas the eight monomer double helix model is more stable than a six monomer double helix model.