Effects of single-base bulges on intercalator binding to small RNA and DNA hairpins and a ribosomal RNA fragment

The way in which a single-base bulge might affect the structure of an RNA helix has been examined by preparing a series of six RNA hairpins, all with seven base pairs and a four-nucleotide loop. Five of the hairpins have single-base bulges at different positions. The intercalating cleavage reagent (methidiumpropyl)-EDTA-Fe(II) [MPE-Fe(II)] binds preferentially at a CpG sequence in the helix lacking a bulge and in four of the five hairpins with bulges. Hairpins with a bulge one or two bases to the 3' side of the CpG sequence bind ethidium 4-5-fold more strongly than the others. V1 RNase, which is sensitive to RNA backbone conformation in helices, detects a conformational change in all of the helices when ethidium binds; the most dramatic changes, involving the entire hairpin stem, are in one of the two hairpins with enhanced ethidium affinity. Only a slight conformational change is detected in the hairpin lacking a bulge. A bulge adjacent to a CpG sequence in a 100-nucleotide ribosomal RNA fragment enhances MPE-Fe(II) binding by an order of magnitude. These results extend our previous observations of bulges at a single position in an RNA hairpin [White, S. A., & Draper, D.E. (1987) Nucleic Acids Res. 15, 4049] and show that (1) a structural change in an RNA helix may be propagated for several base pairs, (2) bulges tend to increase the number of conformations available to a helix, and (3) the effects observed in small RNA hairpins are relevant to larger RNAs with more extensive structure. A bulge in a DNA hairpin identical in sequence with the RNA hairpins does not enhance MPE-Fe(II) binding affinity, relative to a control DNA hairpin. The effects of bulges on ethidium intercalation are evidently modulated by helix structure.     

Bulge loops used to measure the helical twist of RNA in solution.

Bulge loops are commonly found in helical segments of cellular RNAs. When incorporated into long double-stranded RNAs, they may introduce points of flexibility or permanent bend that can be detected by the altered electrophoretic gel mobility of the RNA. We find that a single An or Un bulge loop near the middle of a long RNA helix significantly retards the RNA during polyacrylamide gel electrophoresis if n greater than or equal to 2. The mobility of an RNA containing two A2 bulges various periodically with the number of base pairs between the bulges. We interpret this to mean that A2 bulges varies periodically with the number of base pairs between the bulges. We interpret this to mean that Z2 bulges form torsionally stiff bends in the helix; the gel mobility reaches a minimum when the total helical twist between the bulges rotates the arms of the molecule into a cis conformation. The gel mobilities are proportional to the predicted end-to-end distance of the RNA if the average RNA helical repeat is 11.8 +/- 0.2 bp/turn and there is no helical twist (3 +/- 9 degrees) associated with the bulge (data obtained in 0.15 M Na+). Other sizes and sequences of bulges have very different effects on RNA helix conformation and flexibility. U2 bulges bend the helix to a much smaller degree than A2 bulges, while longer A or U bulge sequences probably allow bends of 90 degrees or more; all of these may be fairly flexible joints.(ABSTRACT TRUNCATED AT 250 WORDS)     

Bulged adenosine influence on the RNA duplex conformation in solution

The RNA single bulge motif is an unpaired residue within a strand of several complementary base pairs. To gain insight into structural changes induced by the presence of the adenosine bulge on RNA duplex, the solution structures of RNA duplex containing a single adenine bulge (5'-GCAGAAGAGCG-3'/5'-CGCUCUCUGC-3') and a reference duplex with all Watson-Crick base pairs (5'-GCAGAGAGCG-3'/5'-CGCUCUCUGC-3') have been determined by NMR spectroscopy. The reference duplex structure is a regular right-handed helix with all of the attributes of an A-type helix. In the bulged duplex, single adenine bulge stacks into the helix, and the bulge region forms a well-defined structure. Both structures were analyzed by the use of calculated helical parameters. Distortions induced by the accommodation of unpaired residue into the helical structure propagate over the entire structure and are manifested as the reduced base pairs inclination and x-displacement. Intrahelical position of bulged adenine A5 is stabilized by efficient stacking with 5'-neighboring residues G4.     

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.     

Substrate structure requirements of the Pac1 ribonuclease from Schizosaccharmyces pombe.

The Pac1 ribonuclease of Schizosaccharomyces pombe is a member of the RNase III family of double-strand-specific ribonucleases. To examine RNA structural features required for efficient cleavage by the Pac1 RNase, we tested a variety of double-stranded and hairpin RNAs as substrates for the enzyme. The Pac1 RNase required substrates that have a minimal helix length of about 20 base pairs. The enzyme cut both strands of the helix at sites separated by two base pairs. However, Pac1 was also able to make a single-stranded cleavage within an internal bulge of an authentic Escherichia coli substrate at the same site chosen by RNase III. Pac1 efficiently degraded the structurally complex adenovirus VA RNA(I), but was inactive against the short HIV-1 TAR RNA hairpin. These results indicate that the Pac1 RNase prefers straight, perfect helices, but it can tolerate internal bulges that do not distort the helix severely. Like its homologue from Saccharomyces cerevisiae, the Pac1 RNase cleaved at two in vivo RNA processing sites in a hairpin structure in the 3' external transcribed spacer of the S. pombe pre-rRNA, suggesting a role for the enzyme in rRNA maturation.     

Conformational transitions in RNA single uridine and adenosine bulge structures: a molecular dynamics free energy simulation study

Extra unmatched nucleotides (single base bulges) are common structural motifs in folded RNA molecules and can participate in RNA-ligand binding and RNA tertiary structure formation. Often these processes are associated with conformational transitions in the bulge region such as flipping out of the bulge base from an intrahelical stacked toward a looped out state. Knowledge of the flexibility of bulge structures and energetics of conformational transitions is an important prerequisite to better understand the function of this RNA motif. Molecular dynamics simulations were performed on single uridine and adenosine bulge nucleotides at the center of eight basepair RNA molecules and indicated larger flexibility of the bulge bases compared to basepaired regions. The umbrella sampling method was applied to study the bulge base looping out process and accompanying conformational and free energy changes. Looping out toward the major groove resulted in partial disruption of adjacent basepairs and was found to be less favorable compared to looping out toward the minor groove. For both uridine and adenosine bulges, a positive free energy change for full looping out was obtained which was approximately 1.5 kcal mol-1 higher in the case of the adenosine compared to the uridine bulge system. The simulations also indicated stable partially looped out states with the bulge bases located in the RNA minor groove and forming base triples with 5'-neighboring basepairs. In the case of the uridine bulge this state was more stable than the intrahelical stacked bulge structure. Induced looping out toward the minor groove involved crossing of an energy barrier of approximately 3.5 kcal mol-1 before reaching the base triple state. A continuum solvent analysis of intermediate bulge states indicated that electrostatic interactions stabilize looped out and base triple states, whereas van der Waals interactions and nonpolar contributions favor the stacked bulge conformation.     

Twist-joints and double twist-joints in RNA structure

Analysis of available RNA crystal structures has allowed us to identify a new family of RNA arrangements that we call double twist-joints, or DTJs. Each DTJ is composed of a double helix that contains two bulges incorporated into different strands and separated from each other by 2 or 3 bp. At each bulge, the double helix is over-twisted, while the unpaired nucleotides of both bulges form a complex network of stacking and hydrogen-bonding with nucleotides of helical regions. In total, we identified 14 DTJ cases, which can be combined in three groups based on common structural characteristics. One DTJ is found in a functional center of the ribosome, another DTJ mediates binding of the pre-tRNA to the RNase P, and two more DTJs form the sensing domains in the glycine riboswitch.     

Binding of the 60-kDa Ro autoantigen to Y RNAs: evidence for recognition in the major groove of a conserved helix.

The 60-kDa Ro autoantigen is normally complexed with small cytoplasmic RNAs known as Y RNAs. In Xenopus oocytes, the Ro protein is also complexed with a large class of variant 5S rRNA precursors that are folded incorrectly. Using purified baculovirus-expressed protein, we show that the 60-kDa Ro protein binds directly to both Y RNAs and misfolded 5S rRNA precursors. To understand how the protein recognizes these two distinct classes of RNAs, we investigated the features of Y RNA sequence and structure that are necessary for protein recognition. We identified a truncated Y RNA that is stably bound by the 60-kDa Ro protein. Within this 39-nt RNA is a conserved helix that is proposed to be the binding site for the Ro protein. Mutagenesis of this minimal Y RNA revealed that binding by the 60-kDa Ro protein requires specific base pairs within the conserved helix, a singly bulged nucleotide that disrupts the helix, and a three-nucleotide bulge on the opposing strand. Chemical probing experiments using diethyl pyrocarbonate demonstrated that, in the presence of the two bulges, the major groove of the conserved helix is accessible to protein side chains. These data are consistent with a model in which the Ro protein recognizes specific base pairs in the conserved helix by binding in the major groove of the RNA. Furthermore, experiments in which dimethyl sulfate was used to probe a naked and protein-bound Y RNA revealed that a structural alteration occurs in the RNA upon Ro protein binding.     

Raman and Raman optical activity (ROA) analysis of RNA structural motifs in Domain I of the EMCV IRES

Raman and Raman optical activity (ROA) spectra were collected for four RNA oligonucleotides based on the EMCV IRES Domain I to assess the contributions of helix, GNRA tetraloop, U·C mismatch base pair and pyrimidine-rich bulge structures to each. Both Raman and ROA spectra show overall similarities for all oligonucleotides, reflecting the presence of the same base paired helical regions and GNRA tetraloop in each. Specific bands are sensitive to the effect of the mismatch and asymmetric bulge on the structure of the RNA. Raman band changes are observed that reflect the structural contexts of adenine residues, disruption of A-form helical structure, and incorporation of pyrimidine bases in non-helical regions. The ROA spectra are also sensitive to conformational mobility of ribose sugars, and verify a decrease in A-type helix content upon introduction of the pyrimidine-rich bulge. Several Raman and ROA bands also clearly show cooperative effects between the mismatch and pyrimidine-rich bulge motifs on the structure of the RNA. The complementary nature of Raman and ROA spectra provides detailed and highly sensitive information about the local environments of bases, and secondary and tertiary structures, and has the potential to yield spectral signatures for a wide range of RNA structural motifs.     

Molecular dynamics study of intrinsic stability in six RNA terminal loop motifs

Single stranded RNA molecules can assume a wide range of tertiary structures beyond the canonical A-form double helix. Certain sequences, termed motifs, are more common than a random distribution would suggest. The existence of such motifs can be rationalized in structural terms. In this study, we have investigated the intrinsic structural stability of RNA terminal loop motifs using multiple MD simulations in explicit water. Representative loops were chosen from the major tetraloop motifs, including also the U-turn motif. Not all loops retain their folded starting structure, but lowering the temperature to 277 K, or adding adjacent base pairs from the stem to which the motif is attached, helps stabilizing the folded loop structure.     

GU receptors of double helices mediate tRNA movement in the ribosome

A new RNA structural motif consisting of two double helices closely packed via minor grooves is found in many places in the ribosome structure. The packing requires that a GU base pair in one helix be packed against a Watson-Crick pair in the other helix. Two such motifs mediate the interaction of the P- and E-tRNA with the large ribosomal subunit. Analysis of the particular positions of these two motifs in view of the available data on occupancy of tRNA-binding sites and structural changes in the ribosome during the elongation cycle suggests a distinct role for each motif in tRNA translocation.     

Distortions induced in double-stranded oligonucleotides by the binding of cis- or trans-diammine-dichloroplatinum(II) to the d(GTG) sequence.

Conformational changes induced in double-stranded oligonucleotides by the binding of trans- or cis-diamminedichloro platinum(II) to the d(GTG) sequence have been characterized by means of melting temperatures, electrophoretic migrations in non-denaturing polyacrylamide gels, reactivities with the artificial nuclease Phenanthroline-copper and with chemical probes. The cis-platinum adduct behaves more as a centre of directed bend than as a hinge joint, the induced bend angle being of the order of 25-30 degrees. The double helix is locally denatured over 2 base pairs (corresponding to the platinated 5'G residue and the central T residue) and is distorted over 4-5 base pairs. The trans-platinum adduct behaves also more as a centre of directed bend than as a hinge joint, the induced bend angle being of the order of 60 degrees. The double helix is locally denatured over 4 base pairs (corresponding to the immediately 5'T residue adjacent to the adduct and to the three base residues of the adduct). Both the cis- and trans-platinum adducts decrease the thermal stability of the double helix.     

Molecular Dynamics Study of Intrinsic Stability in Six RNA Terminal Loop Motifs

Abstract

Single stranded RNA molecules can assume a wide range of tertiary structures beyond the canonical A-form double helix. Certain sequences, termed motifs, are more common than a random distribution would suggest. The existence of such motifs can be rationalized in structural terms. In this study, we have investigated the intrinsic structural stability of RNA terminal loop motifs using multiple MD simulations in explicit water. Representative loops were chosen from the major tetraloop motifs, including also the U-turn motif. Not all loops retain their folded starting structure, but lowering the temperature to 277 K, or adding adjacent base pairs from the stem to which the motif is attached, helps stabilizing the folded loop structure.     

Solution conformation of a bulged adenosine base in an RNA duplex by relaxation matrix refinement

Bulges are common structural motifs in RNA secondary structure and are thought to play important roles in RNA-protein and RNA-drug interactions. Adenosine bases are the most commonly occurring unpaired base in double helical RNA secondary structures. The solution conformation and dynamics of a 25-nucleotide RNA duplex containing an unpaired adenosine, r(GGCAGAGUGCCGC): r(GCGGCACCUGCC) have been studied by NMR spectroscopy and MORASS iterative relaxation matrix structural refinement. The results show that the bulged adenosine residue stacks into the RNA duplex with little perturbation around the bulged region. Most of the bases in the RNA duplex adopt C(3)'-endo conformation, exhibiting the N-type sugar pucker as found in the A form helices. The sugars of the bulged residue and the 5' flanking residue to it are found to exhibit C(2)'-endo conformation. None of the residues are in syn conformation.     

RNA Structure and Stability

RNA molecules fold into specific base-paired conformations that contain single-stranded regions, A-form double helices, hairpin loops, internal loops, bulges, junctions, pseudoknots, kissing hairpins, and so forth. These structural motifs are recognized by proteins, other RNAs, and other parts of the same RNA. The interactions of these structural elements are crucial to the biological functions of the RNA molecules. We describe the different motifs and discuss their thermodynamic stabilities relative to single strands of RNA. The stabilities determine under what conditions they occur and whether they change when interacting with proteins or other ligands.     

Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model.

The analysis and prediction of non-canonical structural motifs in RNA is of great importance for an understanding of the function and design of RNA structures. A hierarchical method has been employed to generate a large variety of sterically possible conformations for a single-base adenine bulge structure in A -form DNA and RNA. A systematic conformational search was performed on the isolated bulge motif and neighboring nucleotides under the constraint to fit into a continuous helical structure. These substructures were recombined with double-stranded DNA or RNA. Energy minimization resulted in more than 300 distinct bulge conformations. Energetic evaluation using a solvation model based on the finite-difference Poisson-Boltzmann method identified three basic classes of low-energy structures. The three classes correspond to conformations with the bulge base stacked between flanking nucleotides (I), location of the bulge base in the minor groove (II) and conformations with a continuous stacking of the flanking helices and a looped out bulge base (III). For the looped out class, two subtypes (IIIa and IIIb) with different backbone geometries at the bulge site could be distinguished. The conformation of lowest calculated energy was a class I structure with backbone torsion angles close to those in standard A -form RNA. Conformations very close to the extra-helical looped out bulge structure determined by X-ray crystallography were also among the low-energy structures. In addition, topologies observed in other experimentally determined bulge structures have been found among low-energy conformers. The implicit solvent model was further tested by comparing an uridine and adenine bulge flanked by guanine:cytosine base-pairs, respectively. In agreement with the experimental observation, a looped out form was found as the energetically most favorable form for the uridine bulge and a stacked conformation in case of the adenine bulge. The inclusion of solvation effects especially electrostatic reaction field contributions turned out to be critically important in order to select realistic low-energy bulge structures from a large number of sterically possible conformations. The results indicate that the approach might be useful to model the three-dimensional structure of non-canonical motifs embedded in double-stranded RNA, in particular, to restrict the number of possible conformations to a manageable number of conformers with energies below a certain threshold.     

New insights into the fundamental role of topological constraints as a determinant of two-way junction conformation

Recent studies have shown that topological constraints encoded at the RNA secondary structure level involving basic steric and stereochemical forces can significantly restrict the orientations sampled by helices across two-way RNA junctions. Here, we formulate these topological constraints in greater quantitative detail and use this topological framework to rationalize long-standing but poorly understood observations regarding the basic behavior of RNA two-way junctions. Notably, we show that the asymmetric nature of the A-form helix and the finite length of a bulge provide a physical basis for the experimentally observed directionality and bulge-length amplitude dependence of bulge induced inter-helical bends. We also find that the topologically allowed space can be modulated by variations in sequence, particularly with the addition of non-canonical GU base pairs at the junction, and, surprisingly, by the length of the 5' and 3' helices. A survey of two-way RNA junctions in the protein data bank confirms that junction residues have a strong preference to adopt looped-in, non-canonically base-paired conformations, providing a route for extending our bulge-directed framework to internal loop motifs and implying a simplified link between secondary and tertiary structure. Finally, our results uncover a new simple mechanism for coupling junction-induced topological constraints with tertiary interactions.     

NMR spectroscopy of RNA duplexes containing pseudouridine in supercooled water

We have performed NMR experiments in supercooled water in order to decrease the temperature-dependent exchange of protons in RNA duplexes. NMR spectra of aqueous samples of RNA in bundles of narrow capillaries that were acquired at temperatures as low as -18 degrees C reveal resonances of exchangeable protons not seen at higher temperatures. In particular, we detected the imino protons of terminal base pairs and the imino proton of a non-base-paired pseudouridine in a duplex representing the eukaryotic pre-mRNA branch site helix. Analysis of the temperature dependence of chemical shift changes (thermal coefficients) for imino protons corroborated hydrogen bonding patterns observed in the NMR-derived structural model of the branch site helix. The ability to observe non-base-paired imino protons of RNA is of significant value in structure determination of RNA motifs containing loop and bulge regions.