On loop folding in nucleic acid hairpin-type structures

In a series of studies, combining NMR, optical melting and T-jump experiments, it was found that DNA hairpins display a maximum stability when the loop part of the molecule comprises four or five nucleotide residues. This is in contrast with the current notion based on RNA hairpin studies, from which it had been established that a maximum hairpin stability is obtained for six or seven residues in the loop. Here we present a structural model to rationalize these observations. This model is based on the notion that to a major extent base stacking interactions determine the stability of nucleic acid conformations. The model predicts that loop folding in RNA is characterized by an extension of the base stacking at the 5'-side of the double helix by five or six bases; the remaining gap can then easily be closed by two nucleotides. Conversely, loop folding in DNA is characterized by extending base stacking at the 3'-side of the double helical stem by two or three residues; again bridging of the remaining gap can then be achieved by one or two nucleotides. As an example of loop folding in RNA the anticodon loop of yeast tRNAPhe is discussed. For the DNA hairpin formed by d(ATCCTAT4TAGGAT) it is shown that the loop structure obtained from molecular mechanics calculations obeys the above worded loop folding principles.     

Slow conformational exchange in DNA minihairpin loops: A conformationalstudy of the circular dumbbell d〈pCGC-TT-GCG-TT〉

In recent years various examples of highly stable two-residue hairpin loops (miniloops) in DNA have been encountered. As the detailed structure and stability of miniloops appear to be determined not only by the nature and sequence of the two bases in the loop, but also by the closing base pair, it is desirable to carry out in-depth studies of especially designed small model DNA compounds. Therefore, a circular DNA dumbbell-like molecule is tailored to consist of a stem of three Watson–Crick base pairs, flanked on each side by a minihairpin loop. The resulting circular DNA decanter 5′-d〈pCGC- TT-GCG- TT〉 -3′ ( I ) is studied in solution by means of nmr spectroscope. At a temperature of 269 K the molecule occurs in a 50/50 mixture of two dumbbell structures (denoted L2L2 and L2L4). L2L2 contains three Watson–Crick C-G base pairs and two two-residue loops (H2-family type) in opposite parts of the molecule. On raising the temperature from 269 to 314 K. The L2L4 conformer becomes increasingly dominant (95% at 314 K). This conformer has a partially disrupted closing G-C base pair in the 5′-GTTC-3′ loop with only one remaining solvent-accessible hydrogen bond between NH_α of the cytosine C(1) and O6 of the guanine G(8), whereas the opposite 5′-CTTG-3′ loop remains stable. The disruption of the C(1)-G(8) base pair in the L2L4 form is correlated with the presence of a syn orientation for the C(1) base at the 5′-3′ loop-stem junction in the 5′-GTTC-3′ loop. The two conformers. L2L2 and L2L4, occur in slow equilibrium (2–20 s^?1). Moderate line broadening of specific ¹H, ¹³C, and ³¹P resonances of residues C(1), G(8), T(9), and T(10) at low temperatures, due to chemical exchange between L2L2 and L2L4, show that the interconversion from an anti to syn conformer in residue C(1) has a small local effect on the structure of the dumbbell. T₁ relaxation measurements, chemical-shift considerations, and complete hand-shape calculations of the exchange process of the G(8) imino proton reveal a possibility for the existence of multiconformational slates in the anti–syn equilibrium. © 1995 John Wiley & Sons, Inc.     

Thermodynamics of melting of the circular dumbbell d〈pCGC-TT-GCG-TT〉

The conformational behavior of DNA minihairpin loops is sensitive to the directionality of the base pair that closes the loop. Especially tailored circular dumbbells, consisting of a stem of three Watson–Crick base pairs capped on each side with a minihairpin loop, serve as excellent model compounds by means of which deeper insight is gained into the relative stability and melting properties of hairpin loops that differ only in directionality of the closing pair: C-G vs G-C. For this reason the thermodynamic properties of the circular DNA decamers 5′-d〈pCGC-TT-GCG-TT〉-3′( I ) and reference compounds 5′-d〈pGGC-TT-GCC-TT≤-3′( II ) and 5′-d(GCG-TC-CGC)-3′( III ) are studied by means of nmr spectroscopy. Molecules I and II adopt dumbbell structures closed on both sides by a two-membered hairpin hop. At low temperature I consists of a mixture of two slowly exchanging forms, denoted L2L2 and L2L4 . The low-temperature L2L2 form is the fully intact minihairpin structure with three Watson–Crick C-G base pairs. The high-temperature form, L2L4 ,contains a partially disrupted closing G-C base pair in the 5′-GTTC-3′ loop, with the cytosine base placed in a syn orientation. The opposite 5′-CTTG-3′ loop remains stable. A study of the noncircular hairpin structure III shows similar conformational behavior for the 5′-GTTC-3′ loop as found in I a syn orientation for C(6) and two slowly exchanging imino proton signals for G(3). The melting point T_m of II was estimated to lie above 365 K. The T_m value of the duplex stem and the 5′-CTTG-3′ loop of the L2L4 form ofIis 352 ± 2 K. The ΔH° is calculated as ?89 ± 10 kJ/mol. The T_m value determined for the individual residues of the 5′-GTTC-3′ loop lies 4°–11° lower. The enthalpy ΔH° of melting the thymine residues in the 5′-GTTC-3′ loop is calculated to be -61± 7 kJ/mol. Thermodynamic data of the equilibrium between the slowly exchanging two- and four-membered loop conformers of I reveal an upper limit for ΔH° of +30 kJ/mol in going from a two-memberedto a four-membered loop, in agreement with the enthalpy difference of +28 k.j/mol between the two loops at the T_m midpoint. For hairpin III the upper limit for ΔH° going from a two-membered to a four-membered loop amounts to ±21 kJ/mol. The mutual exchange rate between the L2 and L4 form in III is estimated as 13.6 s^?1. Our results clearly suggest that small four-way DNA junctions(model for immobilized Holliday junctions) can be designed that consist of a single DNA strandthat features -CTTG-caps on three of the four arms of the junction. © 1995 John Wiley & Sons, Inc.     

Multinuclear NMR studies of DNA hairpins. 2. Sequence-dependent structural variations

The solution conformation of three related DNA hairpins, each with five bases in the loop, is investigated by proton and phosphorus 2D NMR methods. The sequences of the three oligomers are d(CGCGTTGTTCGCG), d(CGCGTTTGTCGCG), and d(CTGCTCTTGTTGAGCAG). One pair of hairpins shares the same stem sequence but differs in the loop, and the appearance of an unusual phosphate torsion in the stem is found to depend on the sequence in the loop of the hairpin. The second pair of hairpins shares the same loop region but differs in the stem sequence in that the base pair which closes the loop is a C-G or G-C pair. The pattern of NOEs reveals that the stacking arrangement in the loop region depends on the base pair that closes the stem. These results suggest that hairpin loop conformation and dynamics are sensitive to small changes in the loop and adjacent stem sequences. These findings are discussed in relation to sequence-dependent thermodynamic changes that have been observed in RNA hairpins.     

Folding of a DNA hairpin loop structure in explicit solvent using replica-exchange molecular dynamics simulations

Hairpin loop structures are common motifs in folded nucleic acids. The 5'-GCGCAGC sequence in DNA forms a characteristic and stable trinucleotide hairpin loop flanked by a two basepair stem helix. To better understand the structure formation of this hairpin loop motif in atomic detail, we employed replica-exchange molecular dynamics (RexMD) simulations starting from a single-stranded DNA conformation. In two independent 36 ns RexMD simulations, conformations in very close agreement with the experimental hairpin structure were sampled as dominant conformations (lowest free energy state) during the final phase of the RexMDs ( approximately 35% at the lowest temperature replica). Simultaneous compaction and accumulation of folded structures were observed. Comparison of the GCA trinucleotides from early stages of the simulations with the folded topology indicated a variety of central loop conformations, but arrangements close to experiment that are sampled before the fully folded structure also appeared. Most of these intermediates included a stacking of the C(2) and G(3) bases, which was further stabilized by hydrogen bonding to the A(5) base and a strongly bound water molecule bridging the C(2) and A(5) in the DNA minor groove. The simulations suggest a folding mechanism where these intermediates can rapidly proceed toward the fully folded hairpin and emphasize the importance of loop and stem nucleotide interactions for hairpin folding. In one simulation, a loop motif with G(3) in syn conformation (dihedral flip at N-glycosidic bond) accumulated, resulting in a misfolded hairpin. Such conformations may correspond to long-lived trapped states that have been postulated to account for the folding kinetics of nucleic acid hairpins that are slower than expected for a semiflexible polymer of the same size.     

The three-dimensional structure of a DNA hairpin in solution two-dimensional NMR studies and structural analysis of d(ATCCTATTTATAGGAT)

The hairpin formed by d(ATCCTATTTATAGGAT) was studied by means of two-dimensional NMR spectroscopy and conformational analysis. Almost all 1H resonances of the stem region could be assigned, while the 1H and 31P spectra of the loop region were interpreted completely; this includes the stereospecific assignment of the H5' and H5" resonances. The derivation of the detailed loop structure was carried out in a stepwise fashion including some improved and new methods for structure determination from NMR data. In the first step, the mononucleotide structures were examined. The conformational space available to the mononucleotide was scanned systematically by varying the glycosidic torsion angle and pseudorotational parameters. Each generated conformer was tested against the experimental J coupling constants and NOE parameters. In the following stage, the structures of dinucleotides and longer fragments were derived. Inter-residue distances between protons were calculated by means of a procedure in which the simulated NOEs, obtained via a relaxation-matrix approach, were fitted to the experimental NOEs without the introduction of a molecular model. In addition, the backbone torsion angles beta, gamma and epsilon were deduced from homocoupling and heterocoupling constants. These data served as constraints in the next step, in which the loop sequence was subjected to a multi-conformer generation procedure. The resulting structures were tested against the mentioned constraints and disregarded if these constraints were violated. This yielded a family of structures for the loop region, confined to a relatively narrow conformational space. A representative conformation was subsequently docked on a B-type stem which fulfilled the structural constraints (derived from the NMR experiments for the stem region) to yield the hairpin structure. Results obtained from subsequent restrained-molecular-mechanics as well as free-molecular-mechanics calculations are in accordance with those obtained by means of the analysis described above. The structure of the hairpin loop is a compactly folded conformation and the first base of the central TTTA region forms a Hoogsteen T-A pair with the fourth base. This Hoogsteen base pair is stacked upon the sixth base pair of the B-type double-helical stem. The second base of the loop is folded into the minor groove, whereas the third base of the loop is partly stacked on the first and fourth bases. The phosphate backbone exhibits a sharp turn between the third and fourth nucleotides of the loop. The peculiar structure of this hairpin loop is discussed in relation to loop folding in DNA and RNA hairpins and in relation to a general model for loop folding.     

Construction of a DNA four-way junction: Design and NMR spectroscopy

In 1964 Holliday postulated the formation of cruciform structures (four-way junctions) in duplex DNA as intermediate in genetic recombination. Since then, many biochemical and biophysical investigations were directed at solving questions concerning structural details of stable four-way junctions. Thus far, NMR spectroscopy played a minor part in these investigations on account of the minimum size of the molecule (expressed as the number of nucleotide residues) that was thought necessary to produce a stable cruciform structure. Indeed, the smallest four-way junction studied thus far by NMR methods was built from four separate DNA strands, each containing 16 nucleotides, a total of 64. Obviously, with such a large structure one runs into assignment problems. We considered the possibility of constructing a stable four-way junction from a single strand of DNA. The underlying idea was to make use of our detailed knowledge of the building principles of stable minihairpin loops. These loops, containing only two nucleotides to bridge the gap between antiparallel strands, are maximally stable in DNA sequences like 5-d(-C-TT-G-), where C and G form a normal Watson-Crick base pair and the two T residues cross the minor groove to form the minihairpin loop. Three of such miniloops could in principle cap three arms of the cruciform. The fourth arm would have an open end. The problem to be solved is to find the minimum length that is required to insure stability of the three closed arms and of the fourth open arm. We were successful with a structure that has three short stems (four base pairs each) and an open-end stem consisting of eight base pairs, a total of 46 residues. NMR experiments, carried out on this molecule in the presence of Mg²⁺, showed details of folding which have never been observed before.     

Solution structure of mRNA hairpins promoting selenocysteine incorporation in Escherichia coli and their base-specific interaction with special elongation factor SELB.

On the basis of chemical probing data, the solution structures of RNA hairpins within fdhF and fdnG mRNAs in Escherichia coli, which both promote selenocysteine incorporation at UGA codons, were derived with the help of computer modeling. We find that these mRNA hairpins contain two separate structural domains that possibly also exert two different functions. The first domain is comprised of the UGA codon, which is included within a complex and distorted double-stranded region. Thereby, release factor 2 might be prevented from binding to the UGA codon to terminate protein synthesis. The second domain is located within the apical loop of the mRNA hairpin structures. This loop region exhibits a defined tertiary structure in which no base is involved in Watson-Crick interactions. The structure of the loop is such that, following a sharp turn after G22 (A22 in fdnG mRNA), bases G23 and U24 are exposed to the solvent on the deep groove side of the supporting helix. Residues C25 and U26 close the loop with a possible single H-bonding interaction between the first and last residues of the loop, 04(U26) and N6(A21). The bulge residues U17 and U18 (in fdhF mRNA), or Ul7 only in fdnG mRNA, point their Watson-Crick positions in the same direction as loop residues G23 and U24 do, and at the same time open up the deep groove at the top of the hairpin helix. Chemical probing data demonstrate that bases G23 and U24 in both mRNA hairpins, as well as residues U17 and Ul7/U18 (for fdhF mRNA) located in a bulge 5' to the loop, are involved directly in binding to special elongation factor SELB in both mRNAs. Therefore, SELB recognizes identical bases within both mRNA hairpins despite differences in their primary sequence, consistent with the derived 3D models for these mRNAs, which exhibit similar tertiary structures. Binding of SELB to the fdhF mRNA hairpin was estimated to proceed with an apparent Kd of 30 nM.     

Structural characterization of three RNA hexanucleotide loops from the internal ribosome entry site of polioviruses.

Structural characteristics of three RNA hairpins from the internal ribosome entry site of poliovirus mRNAs have been determined in solution by NMR. Complete proton, phosphorus and carbon resonance assignments were made for the three 16 nt hairpins. The loop sequences, 5'-AAUCCA , AAACCA and GAACCA, have been shown to be essential for viral mRNA translation. NOESY spectra for the three oligomers were very similar indicating a common three dimensional structure. Stems were A-type duplexes with C3'-endo sugar pucker. In the loops, sequential base stacking interactions were detected for all bases except between U8/A8 and C9, indicating a turn in the phosphodiester backbone at this point. Only one nucleotide, U8/A8, had a sugar pucker which deviated appreciably from C3'-endo. The final base in the loop, A11, exhibited an unusual gauche (-) gamma angle. An ensemble of 10 structures calculated for one hairpin using restrained molecular dynamics shows that the first three bases of the loop are turned so as to be exposed to the exterior of the molecule, while the remaining three bases are in an orientation approximating a continuation of the stem helix. Structure calculations and NMR relaxation measurements indicate that the loop apex is subject to considerable local dynamics.     

NOVEL DNA NANOPARTICLES AND NETWORKS

Joining the thrombin-binding aptamer 5′-d(GGTTGGTGTGGTTGG) and the minihairpin 5′-d(GCGAAGC) leads to new DNA nanoparticles, which are different from rod-like helical double-stranded DNA. Covalent interstrand cross-links in DNA duplexes generated by bifunctional alkadiyne chains were used to build-up the DNA networks.     

Eukaryote specific RNA and protein features facilitate assembly and catalysis of H/ACA snoRNPs

H/ACA Box ribonucleoprotein complexes (RNPs) play a major role in modification of rRNA and snRNA, catalyzing the sequence specific pseudouridylation in eukaryotes and archaea. This enzymatic reaction takes place on a substrate RNA recruited via base pairing to an internal loop of the snoRNA. Eukaryotic snoRNPs contain the four proteins Nop10, Cbf5, Gar1 and Nhp2, with Cbf5 as the catalytic subunit. In contrast to archaeal H/ACA RNPs, eukaryotic snoRNPs contain several conserved features in both the snoRNA as well as the protein components. Here, we reconstituted the eukaryotic H/ACA RNP containing snR81 as a guide RNA in vitro and report on the effects of these eukaryote specific features on complex assembly and enzymatic activity. We compare their contribution to pseudouridylation activity for stand-alone hairpins versus the bipartite RNP. Using single molecule FRET spectroscopy, we investigated the role of the different eukaryote-specific proteins and domains on RNA folding and complex assembly, and assessed binding of substrate RNA to the RNP. Interestingly, we found diverging effects for the two hairpins of snR81, suggesting hairpin-specific requirements for folding and RNP formation. Our results for the first time allow assessing interactions between the individual hairpin RNPs in the context of the full, bipartite snoRNP.     

Unusual DNA duplex and hairpin motifs

Single-stranded DNA or double-stranded DNA has the potential to adopt a wide variety of unusual duplex and hairpin motifs in the presence (trans) or absence (cis) of ligands. Several principles for the formation of those unusual structures have been established through the observation of a number of recurring structural motifs associated with different sequences. These include: (i) internal loops of consecutive mismatches can occur in a B-DNA duplex when sheared base pairs are adjacent to each other to confer extensive cross- and intra-strand base stacking; (ii) interdigitated (zipper-like) duplex structures form instead when sheared G·A base pairs are separated by one or two pairs of purine·purine mismatches; (iii) stacking is not restricted to base, deoxyribose also exhibits the potential to do so; (iv) canonical G·C or A·T base pairs are flexible enough to exhibit considerable changes from the regular H-bonded conformation. The paired bases become stacked when bracketed by sheared G·A base pairs, or become extruded out and perpendicular to their neighboring bases in the presence of interacting drugs; (v) the purine-rich and pyrimidine-rich loop structures are notably different in nature. The purine-rich loops form compact triloop structures closed by a sheared G·A, A·A, A·C or sheared-like G_anti·C_syn base pair that is stacked by a single residue. On the other hand, the pyrimidine-rich loops with a thymidine in the first position exhibit no base pairing but are characterized by the folding of the thymidine residue into the minor groove to form a compact loop structure. Identification of such diverse duplex or hairpin motifs greatly enlarges the repertoire for unusual DNA structural formation.     

Melting studies of short DNA hairpins containing the universal base 5-nitroindole.

Effects of the universal base 5-nitroindole on the thermodynamic stability of DNA hairpins having a 6 bp stem and four base loops were investigated by optical absorbance and differential scanning calorimetry techniques. Melting studies were conducted in buffer containing 115 mM Na(+). Five different modified versions of DNA hairpins containing a 5-nitroindole base or bases substituted at different positions in the stem and loop regions were examined. Thermo-dynamic parameters of the melting transitions estimated from a two-state analysis of optical melting curves and measured directly by calorimetry revealed that the presence of 5-nitroindole bases in the duplex stem or loop regions of short DNA hairpins significantly affects both their enthalpic and entropic melting components in a compensating manner, while the transition free energy varies linearly with the transition temperature. The calorimetrically determined enthalpy and entropy values of the modified hairpins were considerably smaller (43-53%) than the two-state optical parameters, suggesting that solvent effects may be significant in the melting processes of these hairpins. Results of circular dichroism measurements also revealed slight differences between the modified hairpins and the control in both the duplex and melted states, suggesting subtle structural differences between the control and DNA hairpins containing a 5-nitroindole base or bases.     

Analysis of melting transitions of the DNA hairpins formed from the oligomer sequences d[GGATAC(X)4GTATCC] (X = A, T, G, C)

Optical melting transitions of the short DNA hairpins formed from the self-complementary DNA oligomers d[GGATACX4GTATCC] where X = A, T, G, or C measured in 100 mM NaCl are presented. A significant dependence of the melting transitions on loop sequence is observed and transition temperatures, tm, of the hairpins vary from 58.3 degrees C for the T4 loop hairpin to 55.3 degrees C for the A4 loop. A nearest-neighbor sequence-dependent theoretical algorithm for calculating melting curves of DNA hairpins is presented and employed to analyze the experimental melting transitions. Experimental melting curves were fit by adjustment of a single theoretical parameter, Fend(n), the weighting function for a hairpin loop comprised of n single-strand bases. Empirically determined values of Fend(n) provide an evaluation of the free-energy of hairpin loop formation and stability. Effects of heterogeneous nearest-neighbor sequence interactions in the duplex stem on hairpin loop formation were investigated by evaluating Fend(n) in individual fitting procedures using two of the published sets of nearest-neighbor stacking interactions in DNA evaluated in 100 mM NaCl and given by Wartell and Benight, 1985. In all cases, evaluated values of Fend(n) were obtained that provided exact theoretical predictions of the experimental transitions. Results of the evaluations indicate: (1) Evaluated free-energies of hairpin loop formation are only slightly dependent on loop sequences examined. At the transition temperature, Tm, the free-energy of forming a loop of four bases is approximately equal for T4, G4, or C4 loops and varies from 3.9 to 4.8 kcal/mole depending on the set of nearest-neighbor interactions employed in the evaluations. This result suggests, in light of the observed differences in stability between the T4, G4, and C4 loop hairpins, that sequence-dependent interactions between base residues of the loop are most likely not the source of the enhanced stability of a T4 loop.(ABSTRACT TRUNCATED AT 400 WORDS)     

Selective Pressures on RNA Hairpins In Vivo and In Vitro

Comparison of the most stable potential hairpins in the sequences of natural ribozymes with those in the randomized sequences has revealed that the hairpin loop energies are lower than expected by chance. Although these hairpins are not necessarily parts of functional structures, there is a selective pressure to diminish the destabilizing free energies of the hairpin loops. In contrast, no significant bias is observed in the stacking values of the most stable stems. In the ribozymes isolated in vitro the loops of potential hairpins are closer to random values, which can result in less efficient folding rates. Furthermore, the effects of kinetic traps seem to be more significant in the folding pathways of the in vitro isolates due to a potential to form stable stacks incompatible with the functional folds. Similarly to natural ribozyme sequences, the untranslated regions of viral RNAs also form hairpins with relatively low loop free energies. These evolutionary trends suggest ways for efficient engineering of improved RNA constructs on the basis of analysis of in vitro isolates and approaches for the search of regions coding for functional RNA structures in large genome sequences. Received: 12 January 2001 / Accepted: 21 May 2001     

A distinctive RNA fold: the solution structure of an analogue of the yeast tRNAPhe T Psi C domain.

The structure of an analogue of the yeast tRNAPhe T Psi C stem-loop has been determined by NMR spectroscopy and restrained molecular dynamics. The molecule contained the highly conserved modification ribothymidine at its naturally occurring position. The ribothymidine-modified T Psi C stem-loop is the product of the m5U54-tRNA methyltransferase, but is not a substrate for the m1A58-tRNA methyltransferase. Site-specific substitutions and 15N labels were used to confirm the assignment of NOESY cross-peaks critical in defining the global fold of the molecule. The structure is unusual in that the loop folds far over into the major groove of the curved stem. This conformation is stabilized by both stacking interactions and hydrogen bond formation. Furthermore, this conformation appears to be unique among RNA hairpins of similar size. There is, however, a considerable resemblance to the analogous domain in the crystal structure of the full-length yeast tRNAPhe. We believe, therefore, that the structure we have determined may represent an intermediate in the folding pathway during the maturation of tRNA.     

Exploring the complex folding kinetics of RNA hairpins: II. Effect of sequence, length, and misfolded states

The complexity of RNA hairpin folding arises from the interplay between the loop formation, the disruption of the slow-breaking misfolded states, and the formation of the slow-forming native base stacks. We investigate the general physical mechanism for the dependence of the RNA hairpin folding kinetics on the sequence and the length of the hairpin loop and the helix stem. For example, 1), the folding would slow down when a stable GC basepair moves to the middle of the stem; 2), hairpin with GC basepair near the loop would fold/unfold faster than the one with GC near the tail of the stem; 3), within a certain range of the stem length, a longer stem can cause faster folding; and 4), certain misfolded states can assist folding through the formation of scaffold structures to lower the entropic barrier for the folding. All our findings are directly applicable and quantitatively testable in experiments. In addition, our results can be useful for molecular design to achieve desirable fast/slow-folding hairpins, hairpins with/without specific misfolded intermediates, and hairpins that fold along designed pathways.     

Structural and thermodynamic analysis of the conformational states of self-complementary hexanucleotides 5′-d(GCATGC) and 5′-d(GCTAGC) in Aqueous Solution

The conformational states of hexanucleotides 5′-d(GCATGC) and 5′-d(GCTAGC) capable of forming hairpins in aqueous solution were studied by 1D and 2D ¹H NMR and molecular dynamics. The equilibrium thermodynamic parameters were determined for the formation of duplexes and hairpins, and the spatial structures were computed for the GCATGC and GCTAGC conformers. The mobility of the hexamer constituents was evaluated by nanosecond molecular dynamics simulation. The possible causes of the observed difference in the thermodynamic stability of the duplex and the hairpin are discussed.     

RNA simulations: probing hairpin unfolding and the dynamics of a GNRA tetraloop

Simulations of an RNA hairpin containing a GNRA tetraloop were conducted to allow the characterization of its secondary structure formation and dynamics. Ten 10 ns trajectories of the folded hairpin 5'-GGGC[GCAA]GCCU-3' were generated using stochastic dynamics and the GB/SA implicit solvent model at 300 K. Overall, we find the stem to be a very stable subunit of this molecule, whereas multiple loop conformations and transitions between them were observed. These trajectories strongly suggest that extension of the C6 base away from the loop occurs cooperatively with an N-type-->S-type sugar pucker conversion in that residue and that similar pucker transitions are necessary to stabilize other looped-out bases. In addition, a short-lived conformer with an extended fourth loop residue (A8) lacking this stabilizing 2'-endo pucker mode was observed. Results of thermal perturbation at 400 K support this model of loop dynamics. Unfolding trajectories were produced using this same methodology at temperatures of 500 to 700 K. The observed unfolding events display three-state behavior kinetically (including native, globular, and unfolded populations) and, based on these observations, we propose a folding mechanism that consists of three distinct events: (i) collapse of the random unfolded structure and sampling of the globular state; (ii) passage into the folded region of configurational space as stem base-pairs form and gain helicity; and (iii) attainment of proper loop geometry and organization of loop pairing and stacking interactions. These results are considered in the context of current experimental knowledge of this and similar nucleic acid hairpins.     

High-resolution nmr studies of A- and G-containing oligonucleotides

We report the temperature dependence of the H2 and H8 purine ring proton resonances of oligoriboadenylates up to chain length 11, with or without a single guanosine residue at the 5′-end, second position, or 3′-end. The results suggest the following generalizations: (1) Stacking of the bases in a right-handed single-stranded helix is more extensive in the interior of the chain than at the chain ends. (2) The tendency of the terminal base to unstack is greater at the 5′-end than at the 3′-terminus. (3) G stacks more weakly than A, as evidenced by weak stacking of 3′-terminal G. Anomalies were also observed in the unstacking profile of G at the second position in the chain, indicating a conformational anomaly such as looping out of G, thereby allowing adjacent A's to stack together, or adoption by G of some other alternative structure. (4) The results imply that the environment at a given base is influenced by effects of longer range than nearest- or next-nearest-neighbor. Increasing ion condensation as chain length increases may be responsible for the slow approach of oligomer behavior to the properties of the high polymer.