Studies on the structure and stabilizing factor of the CUUCGG hairpin RNA using chemically synthesized oligonucleotides

Accurate Modeling of RNA Hairpins Through the Explicit Treatment of Electronic Polarizability with the Classical Drude Oscillator Force Field

Molecular dynamics (MD) simulations play a crucial role in modeling biomolecular systems in which the electrostatic interactions are critical in dictating the structural and dynamical properties. Thus, the treatment of the electrostatic interactions defined in the underlying force field (FF) strongly affects the simulation accuracy. Most FFs use fixed partial atomic charges to include electrostatic interactions, and therefore lack the electronic polarization response, representing an intrinsic limitation. To address this limitation, polarizable FFs have been developed that treat atomic polarizabilities explicitly. Here we present the application of the all-atom polarizable (Drude) and non-polarizable (CHARMM) nucleic acid FFs in RNA hairpin systems to investigate the impact of polarization on structural properties, dipole moment distributions, and cation interactions. Results show that the presence of polarizability in the FF significantly improves the stabilization of RNA hairpin structure. As expected, the distributions of dipole moments show more fluctuations when simulated using the polarizable FF, with the variation in dipoles contributing to the stabilization of the structures of the loop regions of the RNAs. Contact map analyses between the bases and cations show that the variation of the ion distribution around the entire hairpin is larger for the polarizable FF and the cations occupy the outer hydration shell to a greater extent. The presented results indicate the importance of the explicit treatment of electronic polarizability in molecular simulations of RNA, including in non-canonical regions.

Hairpins under tension: RNA versus DNA

We use optical tweezers to control the folding and unfolding of individual DNA and RNA hairpins by force. Four hairpin molecules are studied in comparison: two DNA and two RNA ones. We observe that the conformational dynamics is slower for the RNA hairpins than for their DNA counterparts. Our results indicate that structures made of RNA are dynamically more stable. This difference might contribute to the fact that DNA and RNA play fundamentally different biological roles in spite of chemical similarity.

Propensities for loop structures of RNA & DNA backbones

RNA oligonucleotides exhibit a large tendency to bend and form a loop conformation which is a major motif contributing to their complex three-dimensional structure. This is in contrast to DNA molecules that predominantly form the double-helix structure. In this paper we investigate by molecular dynamics simulation, as well as, by its combination with the replica-exchange method, the propensity of RNA chains containing the GCUAA pentaloop to form spontaneously a hairpin conformation. The results were then compared with those of analogous hybrid oligonucleotides in which the ribose groups in the loop-region were substituted by deoxyriboses. We find that the RNA oligomers exhibit a marginal excess stability to form loop structures. The equilibrium constant for opening the loop to an extended conformation is twice as large in the hybrid than it is in the RNA chain. Analyses of the hydrogen bonds indicate that the excess stability for forming a hairpin is a result of hydrogen bonds the 2'-hydroxyls in the loop region form with other groups in the loop. Of these hydrogen bonds, the most important is the hydrogen bond donated from the 2'-OH at the first position of the loop to N7 of adenine at the forth position. RNA and DNA backbones are characterized by different backbone dihedral angles and sugar puckering that can potentially facilitate or hamper the hydrogen bonds involving the 2'-OH. Nevertheless, the sugar puckerings of all the pentaloop nucleotides were not significantly different between the two chains displaying the C3'-endo conformation characteristic to the A-form double helix. All of the other backbone dihedrals also did not show any considerable difference in the loop-region except of the δ-dihedral. In this case, the RNA loop exhibited bimodal distributions corresponding to, both, the RNA and DNA backbones, whereas the loop of the hybrid chain behaved mostly as that of a DNA backbone. Thus, it is possible that the behavior of the δ-dihedrals in the loop-region of the RNA adopts conformations that facilitate the intra-nucleotide hydrogen bondings of the 2'-hydroxyls, and consequently renders loop structures in RNA more stable.

Ribose 2'-Hydroxyl Groups Stabilize RNA Hairpin Structures Containing GCUAA Pentaloop

The chemical structure of RNA and DNA is very similar; however, the three-dimensional conformation of these two nucleic acids is very different. Whereas the DNA adopts a repetitive structure of a double-stranded helix, RNA is primarily single stranded with a complex three-dimensional structure in which the hairpin is the most common secondary structure. Apart from the difference between uracil and thymine, the difference in the chemical structure between RNA and DNA is the presence of a hydroxyl group at position 2' of the sugar (ribose) instead of a hydrogen (deoxyribose). In this paper, we present molecular dynamics simulations addressing the contribution of 2'-hydroxyls to the stability of a GCUAA pentaloop motif. The results indicate that the 2'-hydroxyls stabilize the hairpin conformation of the GCUAA pentaloop relative to an analogous oligonucleotide in which the ribose sugars in the loop region were substituted with deoxyriboses. The magnitude of the stabilization was found to be 23.8 ± 4.1 kJ/mol using an alchemical mutations free energy method and 4.2 ± 6.5 kJ/mol using potential of mean force calculations. The latter indicates that in addition to its larger thermodynamic stability the RNA hairpin is also kinetically more stable. We find that the excess stability is a result of intrahairpin hydrogen bonds in the loop region between the 2'-hydroxyls and sugars, bases, and phosphates. The hydrogen bonds with the sugars and phosphates involve predominantly interactions with adjacent nucleotides. However, the hydrogen bonds with the bases involve also interactions between groups on opposite sides of the loop or with the middle base of the loop and are therefore likely to contribute significantly to the stability of the loop. Of these hydrogen bonds, the most frequent is observed between the 2'-hydroxyl at the first position of the pentaloop with N6/N7 of adenine at the forth position, as well as between the 2'-hydroxyl at position -1 with N6 of adenine at the fifth position. Our results contribute to the notion that one of the important roles of the ribose sugars in RNA is to facilitate hairpin formation.

High-resolution NMR structure of an RNA model system: the 14-mer cUUCGg tetraloop hairpin RNA

We present a high-resolution nuclear magnetic resonance (NMR) solution structure of a 14-mer RNA hairpin capped by cUUCGg tetraloop. This short and very stable RNA presents an important model system for the study of RNA structure and dynamics using NMR spectroscopy, molecular dynamics (MD) simulations and RNA force-field development. The extraordinary high precision of the structure (root mean square deviation of 0.3 A) could be achieved by measuring and incorporating all currently accessible NMR parameters, including distances derived from nuclear Overhauser effect (NOE) intensities, torsion-angle dependent homonuclear and heteronuclear scalar coupling constants, projection-angle-dependent cross-correlated relaxation rates and residual dipolar couplings. The structure calculations were performed with the program CNS using the ARIA setup and protocols. The structure quality was further improved by a final refinement in explicit water using OPLS force field parameters for non-bonded interactions and charges. In addition, the 2'-hydroxyl groups have been assigned and their conformation has been analyzed based on NOE contacts. The structure currently defines a benchmark for the precision and accuracy amenable to RNA structure determination by NMR spectroscopy. Here, we discuss the impact of various NMR restraints on structure quality and discuss in detail the dynamics of this system as previously determined.

Retention of function in the DNA homolog of the RNA dopamine aptamer

While it is generally accepted that the functional tertiary structures formed by RNA cannot be replicated by a deoxy version of the same sequence, here we demonstrate conservation of function for a DNA homolog of an RNA aptamer. Using fluorescence anisotropy experiments, this work demonstrates that the all-DNA version of the RNA dopamine aptamer is able to bind dopamine with improved affinity and similar specificity relative to the RNA aptamer. Mutation studies suggest that the binding site is maintained in both structure types. These findings will help to elucidate what sequences and secondary structures allow for retention of function in both RNA and DNA.

Conformationally restricted nucleotides as a probe of structure-function relationships in RNA

We introduce the use of commercially available locked nucleic acids (LNAs) as a functional probe in RNA. LNA nucleotides contain a covalent linkage that restricts the pseudorotation phase of the ribose to C3'-endo (A-form). Introduction of an LNA at a single site thus allows the role of ribose structure and dynamics in RNA function to be assessed. We apply LNA probing at multiple sites to analyze self-cleavage in the lead-dependent ribozyme (leadzyme), thermodynamic stability in the UUCG tetraloop, and the kinetics of recognition of U1A protein by U1 snRNA hairpin II. In the leadzyme, locking a single guanosine residue into the C3'-endo pucker increases the catalytic rate by a factor of 20, despite the fact that X-ray crystallographic and NMR structures of the leadzyme ground state reported a C2'-endo conformation at this site. These results strongly suggest that a conformational change at this position is critical for catalytic function. Functional insights obtained in all three systems demonstrate the highly general applicability of LNA probing in analysis of the role of ribose orientation in RNA structure, dynamics, and function.

Preference for ribose over deoxyribose in loop-closing base pairs of extra stable nucleic acid hairpins

We have investigated the effect of switching ribose to deoxyribose at the closing base-pair of an extra-stable RNA hairpin. Specifically, we studied the sequence 5'-GGAC(UUCG)GUCC, a dodecanucleotide that folds into a well-characterized, "extra stable" RNA hairpin structure. Recently, we showed that hairpins containing a 2',5'-linked (UUCG) loop instead of the native 3',5'-linked loop also exhibit extra-stability (Hannoush and Damha, J. Am. Chem. Soc., 2001, 123, 12368-12374). In this article, we show that the ribose units located at the loop-closing positions (i.e., rC4 and rG9) contribute significantly to the stabilization of RNA hairpins, particularly those containing the 3',5'-UUCG loop. Interestingly, the requirement of rC4 and rG9 is more relaxed for DNA hairpins containing the 2',5'-UUCC loop and, in fact, they may be replaced altogether (ribose--> deoxyribose) without affecting stability. The results broaden our understanding of the behavior of highly stable (UUCG) hairpin loops and how they respond to structural perturbation of the loop-closing base pairs.

Imino proton NMR analysis of HDV ribozymes: nested double pseudoknot structure and Mg2+ ion-binding site close to the catalytic core in solution

Minimized trans-acting HDV ribozyme systems consisting of three (Rz-3) and two (Rz-2) RNA strands were prepared and their folding conformations were analyzed by NMR spectroscopy. The guanosine residues in one of the enzyme components of Rz-3 were labeled with 13C and 15N. Imino proton signals were assigned by analysis of NOESY and HSQC spectra. The results are consistent with the nested double pseudoknot model, which contains novel base pairs (P1.1), as observed in the crystal structure of a genomic HDV ribozyme. The NOE connectivities suggest an additional G:G pair at the bottom of P1.1 and at the top of P4. The effects of temperature and Mg2+ ions on base pairs for Rz-3 were examined. The temperature variation experiment on Rz-3 showed that P3 is the most stable and that P1.1 is as stable as P1 and P2. The imino proton signals of the G:U pair at the bottom of P1 and the top of P1.1, which are close to the cleavage site, showed the largest changes upon Mg2+ titration of Rz-3. The results suggest that the catalytic Mg2+ ion binds to the pocket formed by P1 and L3.

Remarkable stability of hairpins containing 2',5'-linked RNA loops

We report here the results of a comparative study of hairpin loops that differ in the connectivity of phosphodiester linkages (3',5'- versus 2',5'-linkages). In addition, we have studied the effect of changing the stem composition on the thermodynamic stability of hairpin loops. Specifically, we constructed hairpins containing one of six stem duplex combinations, i.e., DNA:DNA ("DD"), RNA:RNA ("RR"), DNA:RNA ("DR"), 2',5'-RNA:RNA ("RR"), 2',5'-RNA:DNA ("RD"), and 2',5'-RNA:2',5'-RNA ("RR"), and one of three tetraloop compositions, i.e., 2',5'-RNA ("R"), RNA ("R"), and DNA ("D"). All hairpins contained the conserved and well-studied loop sequence 5'-...C(UUCG)G...-3' [Cheong et al. Nature 1990, 346, 680-682]. We show that the 2',5'-linked loop C(UUCG)G, i.e.,...C(3'p5')U(2'p5')U(2'p5')C(2'p5')G(2'p5')G(3'p5')..., like its "normal" RNA counterpart, forms an unusually stable tetraloop structure. We also show that the stability imparted by 2',5'-RNA loops is dependent on base sequence, a property that is shared with the regioisomeric 3',5'-RNA loops. Remarkably, we find that the stability of the UUCG tetraloop is virtually independent of the hairpin stem composition (DD, RR, RR, etc.), whereas the native RNA tetraloop exerts extra stability only when the stem is duplex RNA (R:R). As a result, the relative stabilities of hairpins with a 2',5'-linked tetraloop, e.g. ggac(UUCG)gtcc (T(m) = 61.4 degrees C), are often superior to those with RNA tetraloops, e.g. ggac(UUCG)gtcc (T(m) = 54.6 degrees C). In fact, it has been possible to observe the formation of a 2',5'-RNA:DNA hybrid duplex by linking the hybrid's strands to a (UUCG) loop. These duplexes (RD), which are not stable enough to form in an intermolecular complex [Wasner et al. Biochemistry 1998, 37, 7478-7486], were stable at room temperature (T(m) approximately 50 degrees C). Thus, 2',5'-loops have potentially important implications in the study of nucleic acid complexes where structural data are not yet available. Furthermore, they may be particularly useful as structural motifs for synthetic ribozymes and nucleic acid "aptamers".

Molecular dynamics simulations of the denaturation and refolding of an RNA tetraloop

Tetraloops are very abundant structural elements of RNA that are formed by four nucleotides in a hairpin loop which is closed by a double stranded helical stem with some Watson-Crick base pairs. A tetraloop r(GCGAAGGC) was identified from the crystal structure of the central domain of 16S rRNA (727-730) in the Thermus thermophilus 30S ribosomal complex. The crystal structure of the 30S complex includes a total of 104 nucleotides from the central domain of the 16S rRNA and three ribosomal proteins S6, S15 and S18. Independent biochemical experiments have demonstrated that protein S15 plays the role in initiating the formation of the central domain of this complex. In the crystal, the tetraloop interacts with the protein S15 at two sites: one of them is associated with hydrogen bond interactions between residue His50 and nucleotide G730, and the other is associated with the occurrence of residue Arg53 beside A728. This paper uses molecular dynamics (MD) simulations to investigate the protein-dependent structural stability of the tetraloop and demonstrates the folding pathway of this tetraloop via melting denaturation and its subsequent refolding. Three important results are derived from these simulations: (i) The stability of nucleotide A728 appears to be protein dependent. Without the interaction with S15, A728 flips away from stacking with A729. (ii) The melting temperature demonstrated by the simulations is analogous to the results of thermodynamic experiments. In addition, the simulated folding of the tetraloop is stepwise: the native shape of the backbone is formed first; this is then followed by the formation of the Watson- Crick base pairs in the stem; and finally the hydrogen bonds and base stacking in the loop are formed. (iii) The tetraloop structure is similar to the crystal structure at salt concentrations of 0.1 M and 1.0 M used for the simulations, but the refolded structure at 0.1 M salt is more comparable to the crystal structure than at 1.0 M. The results from the simulations using both the Generalized Born continuum model and explicit solvent model (Particle Mesh Ewald) generate a similar pathway for unfolding/refolding of the tetraloop.

Extra stable 2',5'-linked RNA loops

We prepared hairpins that differ in the connectivity of phosphodiester linkages in the loop (RNA vs 2', 5'-RNA). We find that the stability of the extra stable RNA hairpin 5'-rGGAC(UUCG)GUCC-3' is the same as that observed for the hairpin containing a 2',5'RNA loop, i.e. 5'-rGGAC(UUCG)GUCC-3' (where UUCG = U2'p5'U2'p5' C2'p5'G2'p5'). Also significant is the finding that when the stem is duplex DNA, duplex 2',5'-RNA, or DNA:2',5'-RNA, hairpins with the UUCG loop are more stable than those with UUCG loop.

Thermodynamics of 2'-ribose substitutions in UUCG tetraloops

The ribose 2'-hydroxyl group confers upon RNA many unique molecular properties. To better appreciate its contribution to structure and stability and to monitor how substitutions of the 2' hydroxyl can alter an RNA molecule, each loop pyrimidine ribonucleotide in the UUCG tetraloop was substituted with a nucleotide containing either a fluorine (2'-F), hydrogen (2'-H), amino (2'-NH2), or methoxy (2'-OCH3) group, in the context of both the C:G and G:C loop-closing base pair. The thermodynamic parameters of these tetraloop variants have been determined and NMR experiments used to monitor the structural changes resulting from the substitutions. The modified riboses are better tolerated in the G[UUCG]C tetraloop, which may be due to its increased loop flexibility relative to the C[UUCG]G loop. Even for these simple substitutions, the free-energy change reflects a complex interplay of hydrogen bonding, solvation effects, and intrinsic pucker preferences of the nucleotides.

The crystal structure of UUCG tetraloop

All large structured RNAs contain hairpin motifs made of a stem closed by several looped nucleotides. The most frequent loop motif is the UUCG one. This motif belongs to the tetraloop family and has the peculiarity of being highly thermodynamically stable. Here, we report the first crystal structure of two UUCG tetraloops embedded in a larger RNA-protein complex solved at 2.8 A resolution. The two loops present in the asymmetric unit are in a different crystal packing environment but, nevertheless, have an identical conformation. The observed structure is globally close to that obtained in solution by nuclear magnetic resonance. However, subtle differences point to a more detailed picture of the role played by 2'-hydroxyl groups in stabilising this tetraloop.

Experimental and theoretical studies of the effects of deoxyribose substitutions on the stability of the UUCG tetraloop

Experimental and theoretical thermodynamic studies of the consequences of 2'-hydroxyl substitution in the RNA UUCG tetraloop show distinct position dependence consistent with the diverse structural contexts of the four-loop ribose hydroxyls in this motif. The results suggest that even for simple substitutions, such as the replacement of the ribose hydroxyl (2'-OH) with hydrogen (2'-H), the free energy change reflects a complex interplay of hydrogen bonding and solvation effects and is influenced by the intrinsic pucker preferences of the nucleotides. Furthermore, theoretical studies suggest that the effect of these mutations in the single-strand state is sequence dependent, in contrast to what is commonly assumed. Free energy perturbation simulations of ribose-deoxyribose mutations in a single-strand dodecamer and in trinucleotide models suggest that in the denatured state, the magnitude of the free energy change for deoxyribose substitutions is determined to a larger extent by the identity of the nucleotide (A, C, G or U) rather than its structural context. Single-strand mutational effects must be considered when interpreting mutational studies in molecular terms.

Unrestrained stochastic dynamics simulations of the UUCG tetraloop using an implicit solvation model

Three unrestrained stochastic dynamics simulations have been carried out on the RNA hairpin GGAC[UUCG] GUCC, using the AMBER94 force field (Cornell et al., 1995. J. Am. Chem. Soc. 117:5179-5197) in MacroModel 5.5 (Mohamadi et al., 1990. J. Comp. Chem. 11:440-467) and either the GB/SA continuum solvation model (Still et al., 1990. J. Am. Chem. Soc. 112:6127-6129) or a linear distance-dependent dielectric (1/R) treatment. The linear distance-dependent treatment results in severe distortion of the nucleic acid structure, restriction of all hydroxyl dihedrals, and collapse of the counterion atmosphere over the course of a 5-ns simulation. An additional vacuum simulation without counterions shows somewhat improved behavior. In contrast, the two GB/SA simulations (1.149 and 3.060 ns in length) give average structures within 1.2 A of the initial NMR structure and in excellent agreement with results of an earlier explicit solvent simulation (Miller and Kollman, 1997. J. Mol. Biol. 270:436-450). In a 3-ns GB/SA simulation starting with the incorrect UUCG tetraloop structure (Cheong et al., 1990. Nature. 346:680-682), this loop conformation converts to the correct loop geometry (Allain and Varani, 1995. J. Mol. Biol. 250:333-353), suggesting enhanced sampling relative to the previous explicit solvent simulation. Thermodynamic effects of 2'-deoxyribose substitutions of loop nucleotides were experimentally determined and are found to correlate with the fraction of time the ribose 2'-OH is hydrogen bonded and the distribution of the hydroxyl dihedral is observed in the GB/SA simulations. The GB/SA simulations thus appear to faithfully represent structural features of the RNA without the computational expense of explicit solvent.

Observation of an A-DNA to B-DNA transition in a nonhelical nucleic acid hairpin molecule using molecular dynamics

One of the truly challenging problems for molecular dynamics (MD) simulations is demonstrating that the trajectories can sample not only in the vicinity of an experimentally determined structure, but also that the trajectories can find the correct experimental structure starting from some other structure. Frequently these transitions to the correct structure require that the simulations overcome energetic barriers to conformational change. Here we present unrestrained molecular dynamics simulations of the DNA analogs of the RNA 5'-GGACUUCGGUCC-3' hairpin tetraloop. In one simulation we have used deoxyuracil residues, and in the other we have used the native DNA deoxythymine residues. We demonstrate that, on a nanosecond time scale, MD is able to simulate the transitions of both of the A-DNA stems to B-DNA stems within the constraints imposed by the four-base loop that caps the helix. These results suggest that we are now in a position to use MD to address the nature of sequence-dependent structural effects in nonduplex DNA structures.

Theoretical studies of an exceptionally stable RNA tetraloop: observation of convergence from an incorrect NMR structure to the correct one using unrestrained molecular dynamics

We report on the results of five independent and unrestrained molecular dynamics simulations of an RNA tetraloop, r(GGACUUCGGUCC), and its related structures with the loop UUCG sugars changed to deoxyribose. Two separate NMR structures have been reported for the loop portion of this molecule, with the second refinement resulting in a slightly different and more accurate conformation for the loop. The root-mean-square deviation (RMSd) between the two NMR structures, for the loop portions only, is 2.5 A. Our simulations, starting from the two NMR structures, demonstrate that this tetraloop is a very stable and rigid structure with both nanosecond length simulations staying very close to the initial structures. Additionally, both simulations preserved most, if not all, of the NMR-derived interactions and violated very few of the nuclear Overhauser effect (NOE)-derived distances used in the structure refinements. However, when the two NMR structures were simulated with deoxyriboses in the loops instead of the native riboses, the flexibility of the systems increased and we observed a conversion from the incorrect to the correct loop conformation in the simulation which started in the incorrect loop conformation. When the riboses were subsequently re-introduced back into the structure which underwent the conversion, the agreement between this simulation and the one starting from the correct NMR structure was a remarkably low 0.5 A, demonstrating an almost complete convergence from the incorrect to the correct structure using unrestrained molecular dynamics.

Magnesium-mediated conversion of an inactive form of a hammerhead ribozyme to an active complex with its substrate. An investigation by NMR spectroscopy

The effects of magnesium ions on a 32-mer ribozyme (R32) were examined by high resolution NMR spectroscopy. In solution, R32 (without its substrate) consisted of a GAAA loop, stem II, a non-Watson-Crick 3-base pair duplex and a 4-base pair duplex that included a wobble G:U base pair. When an uncleavable substrate RNA (RdC11) was added to R32 without Mg2+ ions, a complex did not form between R32 and RdC11 because the substrate recognition regions of R32 formed intramolecular base pairs (the recognition arms were closed). By contrast, in the presence of Mg2+ ions, the R32-RdC11 complex was formed. Moreover, titration of mixtures of R32 and RdC11 with Mg2+ ions also induced the ribozyme-substrate interaction. Elevated concentrations (1.0 M) of monovalent Na+ ions could not induce the formation of the R32-RdC11 complex. These data suggest that Mg2+ ions are not only important as the true catalysts in the function of ribozyme-type metalloenzymes, but they also induce the structural change in the R32 hammerhead ribozyme that is necessary for establishment of the active form of the ribozyme-substrate complex.

An RNA fragment consisting of the P7 and P9.0 stems and the 3'-terminal guanosine of the Tetrahymena group I intron

On the basis of the nucleotide sequence of Tetrahymena group I intron, we constructed a 31 residue RNA that has the P7 stem and the 3'-terminal guanosine residue (3'-G) with a putative stem-loop structure (P9.0) intervening between them. For this model RNA (P7/P9.0/G), four residues around the guanosine binding site (GBS) in the P7 stem were found to exhibit much lower sensitivities to ribonuclease V1 than those of a variant RNA having adenosine in place of the 3'-G, suggesting that the 3'-G contacts around the GBS. NMR analyses of the imino proton resonances of the P7/P9.0/G RNA indicated that the base pairing in the GBS is retained on the interaction with the 3'-G, and that the two base pairs of the putative P9.0 stem-loop are definitely formed. Comparison of the RNA with its variants with either A (3'-A) or a deletion in place of the 3'-G suggested that the stability of the P9.0 stem-loop is affected by the GBS-3'-G interaction. The melting temperatures of the P9.0 stem-loop were determined from the UV absorbances of these RNAs, which quantitatively indicated that the P9.0 stem-loop is significantly stabilized by the interaction of the GBS with the 3'-G, rather than the 3'-A, and also by direct interaction with divalent cations (Mg2+, Ca2+ or Mn2+). Upon replacement of the G-C base pair by C-G in the GBS of the P7/P9.0/G RNA, the specificity was switched from 3'-G to 3'-A, as in the case of the intact intron.

Conformation of the circular dumbbell d<pCGC-TT-GCG-TT>: structure determination and molecular dynamics

The circular DNA decamer 5'-d<pCGC-TT-GCG-TT>-3' was studied in solution by means of NMR spectroscopy and molecular dynamics in H2O. At a temperature of 269 K, a 50/50 mixture of two dumbbell structures (denoted L2L2 and L2L4) is present. The L2L2 form contains three Watson-Crick C-G base pairs and two two-residue loops is opposite parts of the molecule. On raising the temperature from 269 K to 314 K, the L2L4 conformer becomes increasingly dominant (95% at 314 K). This conformer has a partially disrupted G(anti)-C(syn) closing base pair in the 5'-GTTC-3' loop with only one remaining (solvent-accessible) hydrogen bond between NH alpha of the cytosine dC(1) and O6 of the guanine dG(8). The opposite 5'-CTTG-3' loop remains stable. The two conformers occur in slow equilibrium (rate constant 2-20 s-1). Structure determination of the L2L2 and L2L4 forms was performed with the aid of a full relaxation matrix approach (IRMA) in combination with restrained MD. Torsional information was obtained from coupling constants. Coupling constant analysis (3JHH, 3JHP, 3JCP) gave detailed information about the local geometry around backbone torsion angles beta, gamma, delta, and epsilon, revealing a relatively high flexibility of the 5'-GTTC-3' loop. The values of the coupling constants are virtually temperature-independent. 'Weakly constrained' molecular dynamics in solvent was used to sample the conformational space of the dumbbell. The relaxation matrices from the MD simulation were averaged over <r-3> to predict dynamic NOE volumes. In order to account for the 1:1 conformational mixture of L2L2 and L2L4 present at 271 K, we also included S2 factors and <r-6> averaging of the <r-3>-averaged relaxation matrices. On matrix averaging, the agreement of NOE volumes with experiment improved significantly for protons located in the thermodynamically less stable 5'-GTTC-3' loop. The difference in stability of the 5'-CTTG-3' and 5'-GTTC-3' loops is mainly caused by differences in the number of potential hydrogen bonds in the minor groove and differences in stacking overlap of the base pairs closing the minihairpin loops. The syn conformation for dC(1), favored at high temperature, is stabilized by solvation in the major groove. However, the conformational properties of the dC(1) base, as deduced from R-factor analysis and MD simulations, include a large flexibility about torsion angle chi.

Structural characterization of d(CAACCCGTTG) and d(CAACGGGTTG) mini-hairpin loops by heteronuclear NMR: the effects of purines versus pyrimidines in DNA hairpins

The DNA decamers, d(CAACCCGTTG) and d(CAACGGGTTG) were studied in solution by proton and heteronuclear NMR. Under appropriate conditions of pH, temperature, salt concentration and DNA concentration, both decamers form hairpin conformations with similar stabilities [Avizonis and Kearns (1995) Biopolymers, 35, 187-200]. Both decamers adopt mini-hairpin loops, where the first and last four nucleotides are involved in Watson-Crick hydrogen bonding and the central two nucleotides, CC or GG respectively, form the loop. Through the use of proton-proton, proton-phosphorus and natural abundance proton-carbon NMR experiments, backbone torsion angles (beta, gamma and epsilon), sugar puckers and interproton distances were measured. The nucleotides forming the loops of these decamers were found to stack upon one another in an L1 type of loop conformation. Both show gamma tr and unusual beta torsion angles in the loop-closing nucleotide G7, as expected for mini-hairpin loop formation. Our results indicate that the beta and epsilon torsion angles of the fifth and sixth nucleotides that form the loop and the loop-closing nucleotide G7 are not in the standard trans conformation as found in B-DNA. Although the loop structures calculated from NMR-derived constraints are not well defined, the stacking of the bases in the two different hairpins is different. This difference in the base stacking of the loop may provide an explanation as to why the cytosine-containing hairpin is thermodynamically more stable than the guanine-containing hairpin.

Most compact hairpin-turn structure exerted by a short DNA fragment, d(GCGAAGC) in solution: an extraordinarily stable structure resistant to nucleases and heat

The three-dimensional structure of a short DNA fragment, d(GCGAAGC) exhibiting an extraordinarily stable hairpin structure was determined by nuclear magnetic resonance spectroscopy. Two possible models were obtained by molecular modelling using distance and torsion constraints. Only one of these two models is the correct structure, which can clearly explain all the 1H chemical shifts. d(GCGAAGC) is folded back on itself between A4 and A5, and all the sugars in the fragment adopt the C2'-endo conformation. This compact molecule is stabilized by regular extensive base-stacking interaction within each B-form helical strand of G1C2G3A4 and A5G6C7, and by two G-C and one G3-A5 base pairs. The molecule is hard to differentiate into stem and loop regions, so that we classify it as a turn (hairpin-turn) structure experted by a single-stranded DNA. This highly stacked structure shows high thermostability and strong resistance against nucleases contained in E. coli extracts and in human serum.

Synthetic oligoribonucleotides carrying site-specific modifications for RNA structure-function analysis

Synthetic oligoribonucleotides have become increasingly valuable in studies of RNA structure and function. A range of nucleotide analogues is available which carry modifications in the base, sugar or phosphate moieties. Such analogues have been incorporated into synthetic RNA structures to eliminate or alter individual functional groups in the RNA which potentially can take part in hydrogen-bonding or other non-covalent interactions. Comparisons of the properties of the modified RNAs with unmodified RNA models allow conclusions to be drawn concerning the importance or otherwise of specific functional groups within the RNA. These methods have been applied to studies of RNA interactions with proteins, RNA catalysis and RNA structure.

High-resolution NMR study of a synthetic oligoribonucleotide with a tetranucleotide GAGA loop that is a substrate for the cytotoxic protein, ricin

Ricin is a cytotoxic protein that inactivates ribosomes by hydrolyzing the N-glycosidic bond at position A4324 in eukaryotic 28S rRNA. Its substrate domain forms a double helical stem and a 17-base loop that includes the sequence GAGA, the second adenosine of which corresponds to A4324. Recently, studies of mutant RNAs have shown that the four-nucleotide loop, GAGA, can function as a substrate for ricin. To investigate the structure that is recognized by ricin, we studied the properties of a short synthetic substrate, the dodecaribonucleotide r-CUCAGAGAUGAG, which forms a RNA hairpin structure with a GABA loop and a stem of four base pairs. The results of NMR spectroscopy allowed us to construct the solution structure of this oligonucleotide by restrained molecular-dynamic calculations. We found that the stem region exists as an A-form duplex. 5G and 8A in the loop region form an unusual G:A base pair, and the phosphodiester backbone has a turn between 5G and 6A. This turn seems to help ricin to gain access to 6A which is the only site of depurination in the entire structure. The overall structure of the GAGA loop is similar to those of the GAAA and GCAA loops that have been described but that are not recognized by ricin. Therefore, in addition to the adenosine at the depurination site, the neighboring guanosine on the 3' side (7G) may also play a role in the recognition mechanism together with 5G and 8A.

Design and specificity of hammerhead ribozymes against calretinin mRNA

We obtained a partial sequence of mouse calretinin mRNA from cDNA clones, and designed hammerhead ribozymes to cleave positions within it. With a view to optimising hammerhead ribozymes for eliminating the mRNA in vivo, we varied the length and sequence of the three duplex 'arms' and measured the cleavage of long RNA substrates in vitro at 37 degrees C (as well as 50 degrees C). Precise cleavage occurred, but it could only go to completion with a large excess of ribozyme. The evidence suggests that the rate-limiting step with a large target is not the cleavage, but the formation of the active ribozyme: substrate complex. The efficiency varied unpredictably according to the target site, the length of the substrate RNA, and the length of the ribozyme; secondary structure in vitro may be responsible. We particularly investigated the degree of sequence-specificity. Some mismatches could be tolerated, but shortening of the total basepairing with the substrate to less than 14 bp drastically reduced activity, implying that interaction with weakly-matched RNAs is unlikely to be a serious problem in vivo. These results suggest that specific and complete cleavage of a mRNA in vivo should be possible, given high-level expression of a ribozyme against a favourable target site.

Hammerhead ribozymes: importance of stem-loop II for activity

The activity of several hammerhead ribozyme constructs with constant lengths of stems I and III of 5 nt each but with variously shortened stems II is reported. Stems with 2 bp rather than the conventional 4 bp show essentially unaltered catalytic activity, independent of the composition of the tetraloop. Further reduction in size to 1 bp or 0 bp decreases activity drastically. Inversion of the G10.1.C11.1 bp next to the invariant core leads to a loss in activity, even when the stem consists of 4 bp. Thus, the minimal structural requirement for stem-loop II is a 2-bp stem with a conserved G.C bp. The reduction in catalytic activity is predominantly a result of a decrease of catalytic constant kcat, whereas Km is only slightly affected. Thus, the structural requirement for optimal activity in these constructs where the chemical-cleavage step is rate limiting is determined by the stabilization of the transition state.

The solution structure of a d[C(TTCG)G] DNA hairpin and comparison to the unusually stable RNA analogue

The solution structure of the DNA analogue of the unusually stable r[C(UUCG)G] RNA hairpin, 5'-d[GGA-C(TTCG)GTCC]-3', has been determined by NMR spectroscopy, and its structure has been compared to that of the RNA molecule. The RNA molecule is compact and rigid with a highly structured loop. However, the DNA molecule is much less structured. The DNA hairpin contains a B-form stem of four base pairs. The terminal base pair frays, and the 3'-terminal nucleotides, C11 and C12, are in equilibrium between 2'-endo and 3'-endo conformations. Unlike the RNA loop, the DNA loop contains no syn nucleotides, and there is no evidence for base-base or base-phosphate hydrogen bonding in the loop. The loop is flexible, and reveals no specific internucleotide interactions.

Large scale chemical synthesis, purification and crystallization of RNA-DNA chimeras

RNA-DNA chimeras, in which both DNA and RNA monomers are site-specifically substituted in the same strand, may be prepared only by chemical synthesis. Biochemical studies have revealed a number of surprising and subtle effects resulting from the insertion of either a ribonucleotide into a DNA strand or a deoxyribonucleotide into an RNA strand. The availability of large quantities of these chimeras allows for their crystallization and subsequent x-ray structure determination. We describe a flexible and efficient method for the large-scale preparation of these compounds, their purification, and their crystallization. The methodology is based on a combination of existing DNA phosphoramidite synthons and those recently introduced for the preparation of biochemically active RNA1. We demonstrate that these two different synthons are compatible, produce large quantities of nucleic acid needed for physical studies, and that high resolution diffraction quality crystals may be grown from these chimeras. Of the duplex chimeras synthesized and crystallized, [r(G)d(CGTATACGC)]2, [d(GCGT)r(A)d(TACGC)]2 and [r(GCG)d(TATACCC) + d(GGGTATACGC)] form A-helices and d(CG)r(CG)d(CG)]2 forms a left-handed Z-helix.

Extraordinarily stable mini-hairpins: electrophoretical and thermal properties of the various sequence variants of d(GCGAAAGC) and their effect on DNA sequencing

A small DNA fragment having a characteristic sequence d(GCGAAAGC) has been shown to form an extraordinarily stable mini-hairpin structure and to have an unusually rapid mobility in polyacrylamide gel electrophoresis, even when containing 7M urea. Here, we have studied the stability of the various sequence variants of d(GCGAAAGC) and the corresponding RNA fragments. Many such sequence variants form stable mini-hairpins in a similar manner to the d(GCGAAAGC) sequence. The RNA fragment, r(GCGAAAGC) also forms a mini-hairpin structure with less stability. The DNA mini-hairpins with GAAA or GAA loop are much more stable than DNA and RNA mini-hairpins with other loop sequence so far as has been examined. The stability difference between DNA and RNA mini-hairpins may be deduced to the stem structures formed by DNA (B form) and RNA (A form). The stable hairpins consisting of the GCGAAAGC sequence cause strong band compression on the sequencing gel. This phenomenon should be carefully considered in DNA sequencing.

An NMR study of the conformation and thermodynamics of the circular dumbbell d [formula: see text] Slow exchange between two- and four-membered hairpin loops

The circular DNA decamer 5'-d [formula: see text] 3' is studied in solution by means of NMR spectroscopy. At low temperature the molecule adopts a dumbbell structure with three Watson-Crick C-G base pairs and two two-residue loops in opposite parts of the molecule. On raising the temperature another conformer appears, in which the closing C-G base pair in the 5'-GTTC-3' loop is disrupted, whereas the opposite 5'-CTTG-3' loop remains stable. The two conformers are in slow equilibrium over a limited temperature range.

A thermodynamic study of unusually stable RNA and DNA hairpins

About 70% of the RNA tetra-loop sequences identified in ribosomal RNAs from different organisms fall into either (UNCG) or (GNRA) families (where N = A, C, G, or U; and R = A or G). RNA hairpins with these loop sequences form unusually stable tetra-loop structures. We have studied the RNA hairpin GGAC(UUCG)GUCC and several sequence variants to determine the effect of changing the loop sequence and the loop-closing base pair on the thermodynamic stability of (UNCG) tetra-loops. The hairpin GGAG(CUUG)CUCC with the conserved loop G(CUUG)C was also unusually stable. We have determined melting temperatures (Tm), and obtained thermodynamic parameters for DNA hairpins with sequences analogous to stable RNA hairpins with (UNCG), C(GNRA)G, C(GAUA)G, and G(CUUG)C loops. DNA hairpins with (TTCG), (dUdUCG), and related sequences in the loop, unlike their RNA counterparts, did not form unusually stable hairpins. However, DNA hairpins with the consensus loop sequence C(GNRA)G were very stable compared to hairpins with C(TTTT)G or C(AAAA)G loops. The C(GATA)G and G(CTTG)C loops were also extra stable. The relative stabilities of the unusually stable DNA hairpins are similar to those observed for their RNA analogs.

RNA structure and NMR spectroscopy

RNA molecules perform a wide variety of biological functions, from enzymic activity to storage and propagation of genetic information.

Structure of an unusually stable RNA hairpin

The structure of a very common RNA hairpin, 5'GGAC(UUCG)GUCC, has been determined in solution by NMR spectroscopy. The loop sequence, UUCG, occurs exceptionally often in ribosomal and other RNAs, and may serve as a nucleation site for RNA folding and as a protein recognition site. Reverse transcriptase cannot read through this loop, although it normally transcribes RNA secondary structure motifs. A hairpin with that loop displays unusually high thermodynamic stability; its stability decreases when conserved nucleotides are mutated. The three-dimensional structure for the hairpin was derived from interproton distances and scalar coupling constants determined by NMR using distance geometry, followed by restrained energy minimization. The structure was well-defined despite the conservative use of interproton distances, by constraining the backbone conformation by means of scalar coupling measurements. A mismatch G.U base pair, with syn-guanosine, closes the stem. This hairpin has a loop of only two nucleotides; both adopt C2'-endo sugar pucker. A sharp turn in the phosphodiester backbone is stabilized by a specific cytosine-phosphate contact, probably a hydrogen bond, and by stacking of the cytosine nucleotide on the G.U base pair. The structural features of the loop can explain the unusual thermodynamic stability of this hairpin and its sensitivity to mutations of loop nucleotides.

Non-ribozyme sequences enhance self-cleavage of ribozymes derived from Hepatitis delta virus

Analysis of the self-cleavage of ribozymes derived from the genomic RNA of Hepatitis delta virus (HDV) has revealed that certain co-transcribed vector sequences significantly affect the activity of the ribozyme. Specifically, the t1/2 of self-cleavage for a 135 nucleotide HDV RNA varied, at 42 degrees C, from 5 min to 88 min, depending on the vector-derived sequences flanking the 5' end of the ribozyme. Further analysis suggested that this phenomenon was most likely due to the interaction of vector-derived sequences with a 16 nucleotide region found at the 3' end of the ribozyme. These findings have implications for studies of ribozymes transcribed from cDNA templates, and may provide information regarding the catalytic structure of the HDV ribozyme.

Structural elements in RNA

This chapter describes the RNA structural characteristics that have emerged so far. Folded RNA molecules are stabilized by a variety of interactions, the most prevalent of which are stacking and hydrogen bonding between bases. Many interactions among backbone atoms also occur in the structure of tRNA, although they are often ignored when considering RNA structure because they are not as well-characterized as interactions among bases. Backbone interactions include hydrogen bonding and the stacking of sugar or phosphate groups with bases or with other sugar and phosphate groups. The interactions found in a three-dimensional RNA structure can be divided into two categories: secondary interactions and tertiary interactions. This division is useful for several reasons. Secondary structures are routinely determined by a combination of techniques discussed in chapter, whereas tertiary interactions are more difficult to determine. Computer algorithms that generate RNA structures can search completely through possible secondary structures, but the inclusion of tertiary interactions makes a complete search of possible structures impractical for RNA molecules even as small as tRNA. The division of RNA structure into building blocks consisting of secondary or tertiary interactions makes it easier to describe RNA structures. In those cases in which RNA studies are incomplete, the studies of DNA are described with the rationalization that RNA structures may be analogous to DNA structures, or that the techniques used to study DNA could be applied to the analogous RNA structures. The chapter focuses on the aspects of RNA structure that affect the three-dimensional shape of RNA and that affect its ability to interact with other molecules.