Energetics of the hairpin to mismatched duplex transition of d(GCCGCAGC) on NaCl solution

Influence of cationic molecules on the hairpin to duplex equilibria of self-complementary DNA and RNA oligonucleotides

A self-complementary nucleotide sequence can form both a unimolecular hairpin and a bimolecular duplex. In this study, the secondary structures of the self-complementary DNA and RNA oligonucleotides with different sequences and lengths were investigated under various solution conditions by gel electrophoresis, circular dichroism (CD) and electron paramagnetic resonance (EPR) spectroscopy and a ultraviolet (UV) melting analysis. The DNA sequences tended to adopt a hairpin conformation at low cation concentrations, but a bimolecular duplex was preferentially formed at an elevated cationic strength. On the other hand, fully matched RNA sequences adopted a bimolecular duplex regardless of the cation concentration. The thermal melting experiments indicated a greater change in the melting temperature of the bimolecular duplexes (by approximately 20 degrees C) than that of the hairpin (by approximately 10 degrees C) by increasing the NaCl concentration from 10 mM to 1 M. Hairpin formations were also observed for the palindrome DNA sequences derived from Escherichia coli, but association of the complementary palindrome sequences was observed when spermine, one of the major cationic molecules in a cell, existed at the physiological concentration. The results indicate the role of cations for shifting the structural equilibrium toward a nucleotide assembly and implicate nucleotide structures in cells.

Methylation of the nucleobases in RNA oligonucleotides mediates duplex-hairpin conversion

We have systematically investigated the duplex to hairpin conversion of oligoribonucleotides under the aspect of nucleobase methylation. The first part of our study refers to the self-complementary sequence rCGCGAAUUCGCGA, which forms a stable Watson-Crick base paired duplex under various buffer conditions. It is shown that this sequence is forced to adopt a hairpin conformation if one of the central 6 nt is replaced by the corresponding methylated nucleotide, such as 1-methylguanosine N(2),N(2)-dimethylguanosine, N(6),N(6)-dimethyladenosine (m(6)(2)A) or 3-methyluridine. On the other hand, the duplex structure is retained and even stabilized by replacement of a central nucleotide with N(2)-methylguanosine (m(2)G) or N(4)-methylcytidine. A borderline case is represented by N(6)-methyladenosine (m(6)A). Although generally a duplex-preserving modification, our data indicate that m(6)A in specific strand positions and at low strand concentrations is able to effectuate duplex-hairpin conversion. Our studies also include the ssu ribosomal helix 45 sequence motif, rGACCm(2)GGm(6)(2)Am(6)(2)AGGUC. In analogy, it is demonstrated that the tandem m(6)(2)A nucleobases of this oligoribonucleotide prevent duplex formation with complementary strands. Therefore, it can be concluded that nucleobase methylations at the Watson-Crick base pairing site provide the potential not only to modulate but to substantially affect RNA structure by formation of different secondary structure motifs.

A zipper-like duplex in DNA: the crystal structure of d(GCGAAAGCT) at 2.1 A resolution

BACKGROUND

The replication origin of the single-stranded (ss)DNA bacteriophage G4 has been proposed to fold into a hairpin loop containing the sequence GCGAAAGC. This sequence comprises a purine-rich motif (GAAA), which also occurs in conserved repetitive sequences of centromeric DNA. ssDNA analogues of these sequences often show exceptional stability which is associated with hairpin loops or unusual duplexes, and may be important in DNA replication and centromere function. Nuclear magnetic resonance (NMR) studies indicate that the GCGAAAGC sequence forms a hairpin loop in solution, while centromere-like repeats dimerise into unusual duplexes. The factors stabilising these unusual secondary structure elements in ssDNA, however, are poorly understood.

RESULTS

The nonamer d(GCGAAAGCT) was crystallised as a bromocytosine derivative in the presence of cobalt hexammine. The crystal structure, solved by the multiple wavelength anomalous dispersion (MAD) method at the bromine K-edge, reveals an unexpected zipper-like motif in the middle of a standard B-DNA duplex. Four central adenines, flanked by two sheared G.A mismatches, are intercalated and stacked on top of each other without any interstrand Watson-Crick base pairing. The cobalt hexammine cation appears to participate only in crystal cohesion.

CONCLUSIONS

The GAAA consensus sequence can dimerise into a stable zipper-like duplex as well as forming a hairpin loop. The arrangement closes the minor groove and exposes the intercalated, unpaired, adenines to the solvent and DNA-binding proteins. Such a motif, which can transform into a hairpin, should be considered as a structural option in modelling DNA and as a potential binding site, where it could have a role in DNA replication, nuclease resistance, ssDNA genome packaging and centromere function.

Three-dimensional model of a selective theophylline-binding RNA molecule

A three-dimensional (3D) model for an RNA molecule that selectively binds theophylline but not caffeine is proposed. This RNA, which was found using SELEX (Jenison et al., 1994), is 10,000 times more specific for theophylline (Kn = 320 nM) than for caffeine (KD = 3.5 mM), although the two ligands are identical except for a methyl group substituted at N7 (present only in caffeine). The binding affinity for ten xanthine-based ligands was used to derive a comparative molecular field analysis model (R2 = 0.93 for three components, with cross-validated R2 of 0.73), using the SYBYL and GOLPE programs. A pharmacophoric map was generated to locate steric and electrostatic interactions between theophylline and the RNA binding site. This information was used to identify putative functional groups of the binding pocket and to generate distance constraints. On the basis of a model for the secondary structure (Jenison et al., 1994), the 3D structure of this RNA was then generated using the following method: each helical region of the RNA molecule was treated as a rigid body; single-stranded loops with specific end-to-end distances were generated. The structures of RNA-xanthine complexes were studied using a modified Monte Carlo algorithm. The detailed structure of an RNA-ligand complex model, as well as possible explanations for the theophylline selectivity are discussed.

Thermodynamic investigation of the association of ethidium, propidium and bis-ethidium to DNA hairpins

We have used a combination of calorimetric and spectroscopic techniques to investigate the association of the bis-intercalator ethidium homodimer (bis-ethidium) to short DNA hairpins with sequences: d(GCGCT5GCGC) and d(CGCGT5CGCG). The helix-coil transition of each hairpin, investigated by UV and calorimetric melting protocol, takes place in monomolecular two-state transitions with characteristic enthalpies of approximately 37 kcal mol-1 for disrupting the four dG-dC base pairs of the hairpin stems. Deconvolution of the bis-ethidium-hairpin calorimetric titration curves indicate that each hairpin contains two distinct binding sites for the ligand: a high affinity site in the stem (Kb approximately 10(7)) that accommodates one bis-ethidium molecule and a lower affinity site (Kb approximately 10(6)) located probably at the loop that accommodates two bis-ethidium molecules. The overall stoichiometries of three ligands per hairpin are in agreement with those obtained in continuous variation experiments using visible spectroscopy. The interaction of bis-ethidium for each type of sites results in enthalpy driven reactions, with average binding enthalpies, delta Hb, of -13.1 and -12.1 kcal mol-1 for the stem and loop sites, respectively. Comparison to the thermodynamic profiles of ethidium and propidium binding reveals that the bis-ethidium binding to the stem site of each hairpin has a more favorable free energy term of -1.4 kcal mol-1 and more favorable enthalpy of -4.2 kcal mol-1. These suggest that only one phenanthridine ring of bis-ethidium intercalates in the stem, while the second planar ring is exposed to solvent or weakly associated to the surface of DNA.

Kinetic and thermodynamic characterization of DNA duplex-hairpin interconversion for two DNA decamers: d(CAACGGGTTG) and d(CAACCCGTTG)

The duplex-hairpin interconversion of two DNA decamers, d(CAACGGGTTG) and d(CAACCCGTTG), has been characterized thermodynamically and kinetically by using uv-melting and nmr relaxation methods. Separately, each decamer shows slow exchange between hairpin and duplex conformations. The hairpin conformations have melting points of 47 and 50 degrees C, respectively, and exhibit similar thermodynamic stabilities. The enthalpies of duplex formation, measured by nmr, were found to be very similar (delta HDH = 26 +/- 3 kcal/mole) for both decamers at low salt concentrations (< 50 mM NaCl). However, as the salt concentration was increased the behavior of delta HDH and kinetics is significantly different for each decamer. The d(CAACGGGTTG) decamer forms a duplex containing two central G.G mismatches at high salt and DNA concentration. Based upon the measurement of high interconversion activation energies and a decrease in hairpin formation rate with increasing salt, the interconversion between hairpin and duplex was concluded to proceed by complete strand dissociation. In contrast, the d(CAACCCGTTG) decamer was determined to form a duplex with two centrally located C.C mismatches at pH values less than 6.2, consistent with the formation of a hemiprotonated C+.C mismatch. At pH values greater than 6.4, the hairpin-duplex equilibrium is almost completely shifted toward the hairpin conformation at DNA concentrations of 0.5-7.0 mM and salt concentrations of 10-100 mM. The interconversion of duplex and hairpin conformations was ascertained by means of both kinetic and thermodynamic measurements to proceed by a slightly different mechanism than its complementary decamer. Although the interconversion proceeds by complete strand separation as suggested by high duplex-hairpin interconversion activation enthalpies, the increasing hairpin formation rate with increasing ionic strength as well as the delta HDH dependence on salt indicate that an intermediate internally bulged duplex (no C+.C formation) is stabilized by increasing ionic strength. These data support an interconversion mechanism where an intermediate internally bulged duplex may be the rate limiting step before strand separation.

Dynamics and relative stabilities of parallel- and antiparallel-stranded DNA duplexes

The dynamics and stability of four DNA duplexes are studied by means of molecular dynamics simulations. The four molecules studied are combinations of 4, 15 bases long, single-stranded oligomers, F1, F2, F3, and F4. The sequence of these single strand oligomers are chosen such that F1-F2 and F3-F4 form parallel (ps) DNA double helices, whereas F1-F4 and F2-F3 form antiparallel-stranded (aps) DNA double helices. Simulations were done at low (100 K) and room (300 K) temperatures. At low temperatures the dynamics are quasi-harmonic and the analysis of the trajectories gives good estimates of the low frequency vibrational modes and density of states. These are used to estimate the linear (harmonic) contribution of local fluctuations to the configurational entropy of the systems. Estimates of the differences in enthalpy between ps and aps duplexes show that aps double helices are more stable than the corresponding ps duplexes, in agreement with experiments. At higher temperatures, the distribution of the fluctuations around the average structures are multimodal and estimates of the configurational entropy cannot be obtained. The multi-basin, nonlinear character of the dynamics at 300 K is established using a novel method which extracts large amplitude nonlinear motions from the molecular dynamics trajectories. Our analysis shows that both ps DNA exhibit much larger fluctuations than the two aps DNA. The large fluctuations of ps DNA are explained in terms of correlated transitions in the beta, epsilon, and zeta backbone dihedral angles.

Sampling of the conformations of the d(CGCTGCGGC) hairpin in solution by two-dimensional nuclear magnetic resonance and theoretical methods

Most NMR studies of DNA oligomers have focused on rigid structures that show a strong preference for one or a small set of ground-state conformations. There is an increasing interest in extending NMR methods to investigate DNA systems in which this preference does not exist. A DNA hairpin is one such system where a large number of low-energy structures coexist in solution. In this article we show how 1D/2D NMR data of the d(C1-G2-C3-T4-G5-C6-G7-G8-C9) hairpin are used to map the conformational space of this molecule. First, we characterize the gross morphology of the hairpin by monitoring the exchangeable imino signals in the molecule. Second, we extract a set of inter-proton distances (i.e., the average values and the associated dispersions) for various pairwise interactions by performing full-matrix NOESY simulation with respect to the observed NOESY data for mixing times of 250 and 100 ms. Third, we use these distances as structural constraints to perform a 300-ps molecular dynamics simulation at 500 K. Fourth, we extract 600 snapshots (one after every 0.5 ps) from the MD trajectory and perform constrained energy minimization to map local minima on the sampled energy surface (we call this the rapid temperature quenching step). Fifth, we assign 600 structures to 14 disjoint clusters such that conformationally similar hairpins belong to the same cluster while conformationally distinct hairpins belong to different clusters. Finally, we interpret the NOESY data in terms of conformationally distinct structures by recalculating NOESY contributions taken from representative structures of different clusters. Our analyses clearly demonstrate that the NMR data correspond to an ensemble of distinct structures, i.e., a set of energetically stable but conformationally distinct structures that satisfies the constraints of loop folding in the d(C1-G2-C3-T4-G5-C6-G7-G8-C9) hairpin. Two types of loop folding consistent with NMR data are obtained: (i) a hairpin with two G.C pairs in the stem and four residues in the loop and (ii) a hairpin with two G.C pairs and a reverse wobble G.T pair in the stem plus two residues in the loop.