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.     

Conformational feasibility of a hairpin with two purines in the loop. 5'-d-GGTACIAGTACC-3'

Structural feasibility and conformational requirements for the sequence 5'-d-GGTACIAGTACC-3' to adopt a hairpin loop with I6 and A7 in the loop are studied. It is shown that a hairpin loop containing only two nucleotides can readily be formed without any unusual torsional angles. Stacking is continued on the 5'-side of the loop, with the I6 stacked upon C5. The base A7, on the 3'-side of the loop, can either be partially stacked with I6 or stick outside without stacking. Loop closure can be achieved for both syn and anti conformations of the glycosidic torsions for G8 while maintaining the normal Watson-Crick base pairing with the opposite C5. All torsional angles in the stem fall within the standard B-family of DNA helical structures. The phosphodiesters of the loop have trans,trans conformations. Loop formation might require the torsion about the C4'-C5' bond of G8 to be trans as opposed to the gauche+ observed in B-DNA. These results are discussed in relation to melting temperature studies [Howard et al. (1991) Biochemistry (preceding paper in this issue)] that suggest the formation of very stable hairpin structures for this sequence.     

Role of the closing base pair for d(GCA) hairpin stability: free energy analysis and folding simulations

Hairpin loops belong to the most important structural motifs in folded nucleic acids. The d(GNA) sequence in DNA can form very stable trinucleotide hairpin loops depending, however, strongly on the closing base pair. Replica-exchange molecular dynamics (REMD) were employed to study hairpin folding of two DNA sequences, d(gcGCAgc) and d(cgGCAcg), with the same central loop motif but different closing base pairs starting from single-stranded structures. In both cases, conformations of the most populated conformational cluster at the lowest temperature showed close agreement with available experimental structures. For the loop sequence with the less stable G:C closing base pair, an alternative loop topology accumulated as second most populated conformational state indicating a possible loop structural heterogeneity. Comparative-free energy simulations on induced loop unfolding indicated higher stability of the loop with a C:G closing base pair by ~3 kcal mol(-1) (compared to a G:C closing base pair) in very good agreement with experiment. The comparative energetic analysis of sampled unfolded, intermediate and folded conformational states identified electrostatic and packing interactions as the main contributions to the closing base pair dependence of the d(GCA) loop stability.     

Structure and dynamics of the DNA hairpins formed by tandemly repeated CTG triplets associated with myotonic dystrophy.

Anomalous expansion of the DNA triplet (CTG)n causes myotonic dystrophy. Structural studies have been carried out on (CTG)n repeats in an attempt to better understand the molecular mechanism of repeat expansion. NMR and gel electrophoretic studies demonstrate the presence of hairpin structures for (CTG)5 and (CTG)6 in solution. The monomeric hairpin structure remains invariant over a wide range of salt concentrations (10-200 mM NaCl), DNA concentrations (micromolar to millimolar in DNA strand) and pH (6.0-7.5). The (CTG)n hairpin contains three bases in the loop when n is odd and four bases when n is even. For both odd and even n the stacking and pairing in the stem remain the same, i.e, two hydrogen bond T.T pairs stack with the neighboring G.C pairs. All the nucleotides in (CTG)5 and (CTG)6 adopt C2'-endo, anti conformations. Full-relaxation matrix analysis has been performed to derive the NOE distance constraints from NOESY experiments at seven different mixing times (25, 50, 75, 100, 125, 200 and 500 ms). NOESY-derived distance constraints were subsequently used in restrained molecular dynamics simulations to obtain a family of structures consistent with the NMR data. The theoretical order parameters are computed for H5-H6(cytosines) and H2'-H2" dipolar correlations for both (CTG)5 and (CTG)6 by employing the Lipari-Szabo formalism. Experimental data show that the cytosine in the loop of the (CTG)5 hairpin is slightly more flexible than those in the stem. The cytosine in the loop of the (CTG)6 hairpin is extremely flexible, implying that the dynamics of the four base loop is intrinsically different from that of the three base loop.     

DNA hairpin structures in solution: 500-MHz two-dimensional 1H NMR studies on d(CGCCGCAGC) and d(CGCCGTAGC)

A hairpin structure contains two conformationally distinct domains: a double-helical stem with Watson-Crick base pairs and a single-stranded loop that connects the two arms of the stem. By extensive 1D and 2D 500-MHz 1H NMR studies in H2O and D2O, it has been demonstrated that the DNA oligomers d(CGCCGCAGC) and d(CGCCGTAGC) form hairpin structures under conditions of low concentration, 0.5 mM in DNA strand, and low salt (20 mM NaCl, pH 7). From examination of the nuclear Overhauser effect (NOE) between base protons H8/H6 and sugar protons H1' and H2'/H2", it was concluded that in d(CGCCGCAGC) and d(CGCCGTAGC) all the nine nucleotides display average (C2'-endo,anti) geometry. The NMR data in conjunction with molecular model building and solvent accessibility studies were used to derive a working model for the hairpins.     

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.     

The Hairpin Structure of a Topoisomerase II Site DNA Strand Analyzed by Combined NMR and Energy Minimization Methods

We report on the structural study of the single-stranded 19mer oligonucleotide d(AGCTTATC-ATC-GATAAGCT) 22(+). This corre sponds to the 15-to-33(+) strand of pBR322 DNA belonging to a strong cleavage site (site 22) for topoisomerase II coupled to antitumor drugs VP-16 or ellipticine. The partially self-complementary nature of this oligonucleotide makes likely its folding into a hairpin structure. To assess this property we carried out a quantitative analysis based on joint calculations and NMR experiments. The latter required two-dimensional (NOESY, P-COSY, TOCSY and proton-detected¹H-³¹P), and three- dimensional (NOESY-TOCSY) spectra to achieve the assignment of the overcrowded sugar H4′ ad H5′/H5′′ proton region. For molecular modeling, the JUMNA program was used together with NMR constraints; namely, the distances and the backbone torsion angles provided by NOEs and homo- and heteronuclear coupling constants. Experimental results proved that the 19mer oligonucleotide adopted a stable hairpin structure characterized by an eight base-pair stem and a three-membered loop (central-ATC-segment). Homonuclear¹H-¹H and heteronuclear¹H-³¹P coupling constant measurements provided information on the confor mational heterogeneity of the sugar and phosphate groups within both the stem and the loop. Restrained energy minimizations starting with different structures resulted in a family of closely related structures. All low-energy molecules presented the same, rather compact, folded structure with the base-stacking continuing into the loop, a sharp turn occurring between residues T10 and C11, and strong backbone distortions at the loop-stem junction.     

Hairpin opening by single-strand-specific nucleases.

DNA molecules with covalently sealed (hairpin) ends are probable intermediates in V(D)J recombination. According to current models hairpin ends are opened to produce short single-stranded extensions that are thought to be precursors of a particular type of extra nucleotides, termed P nucleotides, which are frequently present at recombination junctions. Nothing is known about the activities responsible for hairpin opening. We have used two single-strand-specific nucleases to explore the effects of loop sequence on the hairpin opening reaction. Here we show that a variety of hairpin ends are opened by P1 nuclease and mung bean nuclease (MBN) to leave short, 1-2 nt single-stranded extensions. Analysis of 22 different hairpin sequences demonstrates that the terminal 4 nt of the hairpin loop strongly influence the sites of cleavage. Correlation of the nuclease digestion patterns with structural (NMR) data for some of the hairpin loops studied here provides new insights into the structural features recognized by these enzymes.     

Strong binding of single-stranded DNA by stem-loop oligonucleotides.

We report that oligodeoxynucleotides which form stem-loop hairpin structures and which have pyrimidine-rich loops can form strong complexes with complementary single-stranded DNA sequences. Stem-loop oligonucleotides were constructed with a 25-nt T-rich loop and with variable Watson-Crick stems. The complexes of these oligomers with the sequence dA8 were studied by thermal denaturation. Evidence is presented that the complexes are one-to-one, bimolecular complexes in which the pyrimidine loop bases comprise the outer strands in a pyr.pur.pyr triplex, in effect chelating the purine strand in the center of the loop. Melting temperatures for the loop complexes are shown to be up to 29 degrees C higher than Watson-Crick duplex of the same length. It is shown that the presence of a stem increases stability of the triplex relative to an analogous oligomer without a stem. The effect of stem length on the stability of such a complex is examined. Such hairpin oligomers represent a new approach to the sequence-specific binding of single-stranded RNA and DNA. In addition, the finding raises the possibility that such a complex may exist in natural RNA folded sequences.     

Solution structure of a DNA double helix incorporating four consecutive non-Watson-Crick base-pairs

A series of DNA 21-mers containing a variety of the 4 x 4 internal loop sequence 5'-CAAG-3'/3'-ACGT-5' were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers exhibit excellent resolution in the NMR spectra and reveal many unusual NOEs (nuclear Overhauser effect) that allow for the detailed characterization of a DNA hairpin incorporating a track of four different non-Watson-Crick base-pairs in the stem. These include a wobble C.A base-pair, a sheared A.C base-pair, a sheared A.G base-pair, and a wobble G.T base-pair. Significantly different twisting angles were observed between the base-pairs in internal loop that results with excellent intra-strand and inter-strand base stacking within the four consecutive mismatches and the surrounding canonical base-pairs. This explains why it melts at 52 degrees C even though five out of ten base-pairs in the stem adopt non-Watson-Crick pairs. However, the 4 x 4 internal loop still fits into a B-DNA double helix very well without significant change in the backbone torsion angles; only zeta torsion angles between the tandem sheared base-pairs are changed to a great extent from the gauche(-) domain to the trans domain to accommodate the cross-strand base stacking in the internal loop. The observation that several consecutive non-canonical base-pairs can stably co-exist with Watson-Crick base-pairs greatly increases the limited repertoire of irregular DNA folds and reveals the possibility for unusual structural formation in the functionally important genomic regions that have potential to become single-stranded.     

A DNA hairpin with a single residue loop closed by a strongly distorted Watson-Crick G x C base-pair

Our previous NMR and modeling studies have shown that the single-stranded 19mer oligonucleotides d(AGCTTATC-ATC-GATAA GCT) -ATC- and d(AGCTTATC-GAT-GATAAGCT) -GAT- encompassing the strongest topoisomerase II cleavage site in pBR322 DNA could form stable hairpin structures. A new sheared base-pair, the pyrimidine-purine C x A, was found to close the single base -ATC- loop, while -GAT- displayed a flexible loop of three/five residues with no stabilizing interactions. Now we report a structural study on -GAC-, an analog of -GAT-, derived through the substitution of the loop residue T by C. The results obtained from NMR, non-denaturing PAGE, UV-melting, circular dichroism experiments and restrained molecular dynamics indicate that -GAC- adopts a hairpin structure folded through a single residue loop. In the -GAC- hairpin the direction of the G9 sugar is reversed relative to the C8 sugar, thus pushing the backbone of the loop into the major groove. The G9 x C11 base-pair closing the loop is thus neither a sheared base-pair nor a regular Watson-Crick one. Although G9 and C11 are paired through hydrogen bonds of Watson-Crick type, the base-pair is not planar but rather adopts a wedge-shaped geometry with the two bases stacked on top of each other in the minor groove. The distortion decreases the sugar C1'-C1' distance between the paired G9 and C11, to 8 A versus 11 A in the standard B-DNA. The A10 residue at the center of the loop interacts with the G9 x C11 base-pair, and seems to contribute to the extra thermal stability displayed by -GAC- compared to -GAT-. Test calculations allowed us to identify the experimental NOEs critical for inducing the distorted G.C Watson-Crick base-pair. The preference of -GAC- for a hairpin structure rather than a duplex is confirmed by the diffusion constant values obtained from pulse-field gradient NMR experiments. All together, the results illustrate the high degree of plasticity of single-stranded DNAs which can accommodate a variety of turn-loops to fold up on themselves.     

Hairpin-loop formation by inverted repeats in supercoiled DNA is a local and transmissible property

Short inverted repeat sequences adopt hairpin stem-loop type structures in supercoiled closed circular DNA molecules, demonstrated by S1 nuclease cleavage. Fine mapping of cleavage frequencies is in good agreement with expected cleavage patterns based upon the interaction between an unpaired loop and a sterically bulky enzyme molecule. Whilst the topological properties of underwound DNA circles depend ultimately upon reduced linkage, necessarily a global molecular property, hairpin loop formation is an essentially local property. Thus molecular size is unimportant for the S1 hypersensitivity of the Co1E1 inverted repeat. Furthermore, a 440 bp Sau3AI, EcoRI fragment of Co1E1 which contains the inverted repeat has been cloned into pBR322 whereupon it exhibits S1 cleavage similar to Co1E1 in the supercoiled recombinant molecule. The effect is therefore both local and transmissible. Direct competition, between inverted repeats in the recombinant, coupled with close examination of flanking sequences, enables some simple 'rules' for base pairing in hairpin loops to be formulated. Whilst limited G-T and A-C base pairing appears not to be destabilising, A-G, T-C or loop outs are highly destabilising.     

Interaction of DNA hairpin loops and a complementary strand by a triplet of base pairs

Hypothesis of non-enzymatic recognition of primordial tRNA and mRNA precursors is experimentally approached. DNA hairpins containing a different number of deoxyguanosine residues in the loop are analyzed for their binding ability to a chemically fixed single-strand of oligo(dC). In presence of small Mg2+ concentration a hairpin with five dG residues in the loop is adsorbed to affinity matrix. Comparison of elution temperatures of hairpin oligonucleotides with those of single-stranded oligoguanylic acids with length of the loop indicates, that smallest loop able to bind forms a triplet of base pairs.     

The formation of adjacent triplex-duplex domainsin DNA.

The ability of single-stranded DNA oligomers to form adjacent triplex and duplex domains with two DNA structural motifs was examined. Helix-coil transition curves and a gel mobility shift assay were used to characterize the interaction of single-stranded oligomers 12-20 nt in length with a DNA hairpin and with a DNA duplex that has a dangling end. The 12 nt on the 5'-ends of the oligomers could form a triplex structure with the 12 bp stem of the hairpin or the duplex portion of the DNA with a dangling end. The 3'-ends of the 17-20 nt strands could form Watson-Crick pairs to the five base loop of the hairpin or the dangling end of the duplex. Complexes of the hairpin DNA with the single-stranded oligomers showed two step transitions consistent with unwinding of the triplex strand followed by hairpin denaturation. Melting curve and gel competition results indicated that the complex of the hairpin and the 12 nt oligomer was more stable than the complexes involving the extended single strands. In contrast, results indicated that the extended single-stranded oligomers formed Watson-Crick base pairs with the dangling end of the duplex DNA and enhanced the stability of the adjacent triplex region.     

Strong Binding of Single-stranded DNA by Stem-loop Oligonucleotides

Abstract

We report that oligodeoxynucleotides which form stem-loop hairpin structures and which have pyrimidine-rich loops can form strong complexes with complementary single-stranded DNA sequences. Stem-loop oligonucleotides were constructed with a 25-nt T-rich loop and with variable Watson-Crick stems. The complexes of these oligomers with the sequence dA_gwere studied by thermal denaturation. Evidence is presented that the complexes are one-to-one, bimolecular complexes in which the pyrimidine loop bases comprise the outer strands in a pyr · pur · pyr triplex, in effect chelating the purine strand in the center of the loop. Melting temperatures for the loop complexes are shown to be up to 29 °C higher than Watson- Crick duplex of the same length. It is shown that the presence of a stem increases stability of the triplex relative to an analogous oligomer without a stem. The effect of stem length on the stability of such a complex is examined. Such hairpin oligomers represent a new approach to the sequence-specific binding of single-stranded RNA and DNA. In addition, the finding raises the possibility that such a complex may exist in natural RNA folded sequences.     

Properties and Anti-HIV Activity of Nicked Dumbbell Oligonucleotides

Abstract

We have designed a new type of oligodeoxyribonucleotide. These oligodeoxyribonucleotides form two hairpin loop structures with base pairs (sense and antisense) in the double helical stem at the 3′ and 5′-ends (nicked dumbbell oligonucleotides). The nicked dumbbell oligonucleotides are molecules with free ends that are more resistant to exonuclease attack. Furthermore, the nicked dumbbell oligonucleotide containing phosphorothioate (P=S) bonds in the hairpin loops has increased nuclease resistance, as compared to the unmodified nicked oligonucleotide. The binding of the nicked dumbbell oligonucleotide to RNA is lower than that of a single-stranded DNA. We also describe the anti-HIV activity of nicked dumbbell oligonucleotides.

     

The nonenzymatic hydrolysis of oligoribonucleotides. VII. Structural elements affecting hydrolysis

Several elements of oligoribonucleotide structure are important for efficient hydrolysis. We have found that the following factors influence oligoribonucleotide hydrolysis: (i) single-stranded structure of RNA flanking the scissile phosphodiester bond, (ii) the substituent on atom C-5 of the uridine adjacent to the cleaved internucleotide bond, (iii) the position of the scissile UA phosphodiester bond within a hairpin loop, (iv) the concentration of formamide, urea, ethanol and sodium chloride.