Conformation and dynamics of a left-handed Z-DNA hairpin: studies of d(CGCGCGTTTTCGCGCG) in solution

DNA-based assay for calorimetric determination of protein concentrations in pure or mixed solutions

It was recently reported that values of the transition heat capacities, as measured by differential scanning calorimetry, for two globular proteins and a short DNA hairpin in NaCl buffer are essentially equivalent, at equal concentrations (mg/mL). To validate the broad applicability of this phenomenon, additional evidence for this equivalence is presented that reveals it does not depend on DNA sequence, buffer salt, or transition temperature (Tm). Based on the equivalence of transition heat capacities, a calorimetric method was devised to determine protein concentrations in pure and complex solutions. The scheme uses direct comparisons between the thermodynamic stability of a short DNA hairpin standard of known concentration, and thermodynamic stability of protein solutions of unknown concentrations. Sequences of two DNA hairpins were designed to confer a near 20°C difference in their Tm values. In all cases, evaluated protein concentrations determined from the DNA standard curves agreed with the UV-Vis concentration for monomeric proteins. For multimeric proteins evaluated concentrations were greater than determined by UV-Vis suggesting the calorimetric approach can also be an indicator of molecular stoichiometry.

Calorimetric analysis using DNA thermal stability to determine protein concentration

It was recently reported for two globular proteins and a short DNA hairpin in NaCl buffer that values of the transition heat capacities, Cp,DNA and Cp,PRO, for equal concentrations (mg/mL) of DNA and proteins, are essentially equivalent (differ by less than 1%). Additional evidence for this equivalence is presented that reveals this phenomenon does not depend on DNA sequence, buffer salt, or Tm. Sequences of two DNA hairpins were designed to confer a near 20°C difference in their Tm's. For the molecules, in NaCl and CsCl buffer the evaluated Cp,PRO and Cp,DNA were equivalent. Based on the equivalence of transition heat capacities, a calorimetric method was devised to determine protein concentrations in pure and complex solutions. The scheme uses direct comparisons between the thermodynamic stability of a short DNA hairpin standard of known concentration, and thermodynamic stability of protein solutions of unknown concentrations. In all cases, evaluated protein concentrations determined from the DNA standard curve agreed with the UV-Vis concentration for monomeric proteins. For samples of multimeric proteins, streptavidin (tetramer), Herpes Simplex Virus glycoprotein D (trimer/dimer), and a 16 base pair DNA duplex (dimer), evaluated concentrations were greater than determined by UV-Vis by factors of 3.94, 2.65, and 2.15, respectively.

Oxidative Modification of Guanine in a Potential Z-DNA-Forming Sequence of a Gene Promoter Impacts Gene Expression

One response to oxidation of guanine (G) to 8-oxo-7,8-dihydroguanine (OG) in a gene promoter is regulation of mRNA expression suggesting an epigenetic-like role for OG. A proposed mechanism involves G oxidation within a potential G-quadruplex-forming sequence (PQS) in the promoter, enabling a structural shift from B-DNA to a G-quadruplex fold (G4). When OG was located in the coding vs template strand, base excision repair led to an on/off transcriptional switch. Herein, a G-rich, potential Z-DNA-forming sequence (PZS) comprised of a d(GC) _n repeat was explored to determine whether oxidation in this motif was also a transcriptional switch. Bioinformatic analysis found 1650 PZSs of length >10 nts in the human genome that were overrepresented in promoters and 5'-UTRs. Studies in human cells transfected with a luciferase reporter plasmid in which OG was synthesized in a PZS context in the promoter found that a coding strand OG increased expression and a template strand OG decreased expression. The initial base excision repair product of OG, an abasic site (AP), was also found to yield similar expression changes as OG. Biophysical studies on model Z-DNA strands found OG favored a shift in the equilibrium to Z-DNA from B-DNA, while an AP disrupted Z-DNA to favor a hairpin, placing AP in the loop where it is a poor substrate for the endonuclease APE1. Overall, the impact of OG and AP in a PZS on gene expression was similar to that in a PQS but reduced in magnitude.

Free-energy landscape of a hyperstable RNA tetraloop

We report the characterization of the energy landscape and the folding/unfolding thermodynamics of a hyperstable RNA tetraloop obtained through high-performance molecular dynamics simulations at microsecond timescales. Sampling of the configurational landscape is conducted using temperature replica exchange molecular dynamics over three isochores at high, ambient, and negative pressures to determine the thermodynamic stability and the free-energy landscape of the tetraloop. The simulations reveal reversible folding/unfolding transitions of the tetraloop into the canonical A-RNA conformation and the presence of two alternative configurations, including a left-handed Z-RNA conformation and a compact purine Triplet. Increasing hydrostatic pressure shows a stabilizing effect on the A-RNA conformation and a destabilization of the left-handed Z-RNA. Our results provide a comprehensive description of the folded free-energy landscape of a hyperstable RNA tetraloop and highlight the significant advances of all-atom molecular dynamics in describing the unbiased folding of a simple RNA secondary structure motif.

Single C8-Arylguanine modifications render oligonucleotides in the Z-DNA conformation under physiological conditions

Z-DNA is the only DNA conformation that has a left-handed helical twist. Although Z-DNA has been implicated in both carcinogenesis and mutagenesis, its specific biological role remains uncertain. We have demonstrated that the formation of C8-arylguanine DNA adducts, derived from arylhydrazines, shifts the B/Z-DNA equilibrium toward the Z-DNA conformation in d(CG)5 sequences. However, our previous work examined the effect of two adducts in the duplex, and it was unclear whether the two base modifications were working together to cause the equilibrium shift toward the Z-DNA conformation. Here we report the synthesis and characterization of a hairpin oligonucleotide sequence (d(CG)5T4(CG)5) containing only one C8-arylguanine modified base. The unmodified hairpin and the previously studied unmodified double-stranded oligonucleotide were conformationally similar, and each required ∼3 M NaCl to yield a B-/Z-DNA ratio of 1:1. The introduction of a single C8-arylguanine modification significantly reduced the NaCl concentration needed to produce a 1:1 B-/Z-DNA ratio in the hairpin. Further, the addition of MgCl2 and spermine to the C8-arylguanine-modified hairpin shifts the B/Z-DNA equilibrium such that the Z form predominated under physiological conditions. NMR and molecular modeling indicated the conformational effects produced by the C8-arylguanine modification occurred locally at the site of modification while CD data demonstrated that the C8-arylguanine-modified base destabilized the B form. Additionally, our data show that adopting the Z-DNA conformation is preferred over denaturation to the single-stranded form. Finally, the conformational effects of the C8-arylguanine modifications were not additive and the introduction of any such modifications drive Z-DNA formation under physiological conditions, which may provide a novel carcinogenesis mechanism where DNA adducts confer their carcinogenicity through a Z-DNA-mediated mechanism.

The Zα domain of fish PKZ converts DNA hairpin with d(GC)(n) inserts to Z-conformation

PKZ, protein kinase containing Z-DNA domains, is a novel member of the vertebrate eIF2α kinase family. Containing a catalytic domain in C-terminus and two Z-DNA binding domains (Zα1 and Zα2) in N-terminus, PKZ can be activated through the binding of Zα to Z-DNA. However, the regulatory function of PKZ Zα remains to be established. Here, to understand the impact of PKZ Zα on DNA conformational transition, wild-type Zα1Zα2 and 11 mutant proteins were expressed and purified. At the same time, several different lengths of DNA hairpins-d(GC)nT4(GC)n (n = 2-6) and an RNA hairpin-r(GC)6T4(GC)6 were synthesized. The effects of Zα1Zα2 and mutant proteins on the conformation of these synthetic DNA or RNA hairpins were investigated by using circular dichroism spectrum and gel mobility shift assays. The results showed that DNA hairpins retained a conventional B-DNA conformation in the absence of Zα1Zα2, while some of the DNA hairpins (n≥3) were converted to Z-conformation under Zα1Zα2 induction. The tendency was proportionally associated with the increasing amount of GC repeat. In comparison with Zα1Zα2, Zα1Zα1 rather than Zα2Zα2 displayed a higher ability in converting d(GC)6T4(GC)6 from B- to Z-DNA. These results demonstrated that Zα1 sub-domain played a more essential role in the process of B-Z conformational transition than Zα2 sub-domain did. Mutant proteins (K34A, N38A, R39A, Y42A, P57A, P58A, and W60A) could not convert d(GC)6T4(GC)6 into Z-DNA, whereas S35A or K56A retained some partial activities. Interestingly, Zα1Zα2 was also able to induce r(GC)6T4(GC)6 RNA from A-conformation to Z-conformation under appropriate conditions.

Activity, folding and Z-DNA formation of the 8-17 DNAzyme in the presence of monovalent ions

The effect of monovalent ions on both the reactivity and global folding of the 8-17 DNAzyme is investigated, and the results are compared with those of the hammerhead ribozyme, which has similar size and secondary structure. In contrast to the hammerhead ribozyme, the 8-17 DNAzyme activity is not detectable in the presence of 4 M K(+), Rb(+), or Cs(+) or in the presence of 80 mM, [Co(NH(3))(6)](3+). Only 4 M Li(+), NH(4)(+) and, to a lesser extent, Na(+) conferred detectable activity. The observed rate constants (k(obs) approximately 10(-3) min(-1) for Li(+) and NH(4)(+)) are approximately 1000-fold lower than that in the presence of 10 mM Mg(2+), and approximately 200,000-fold slower than that in the presence of 100 microM Pb(2+). Since the hammerhead ribozyme displays monovalent ion-dependent activity that is often within approximately 10-fold of divalent metal ion-dependent activity, these results suggest that the 8-17 DNAzyme, obtained by in vitro selections, has evolved to have a more stringent divalent metal ion requirement for high activity as compared to the naturally occurring ribozymes, making the 8-17 DNAzyme an excellent choice as a Pb(2+) sensor with high selectivity. In contrast to the activity data, folding was observed in the presence of all the monovalent ions investigated, although those monovalent ions that do not support DNAzyme activity have weaker binding affinity (K(d) approximately 0.35 M for Rb(+) and Cs(+)), while those that confer DNAzyme activity possess stronger affinity (K(d) approximately 0.22 M for Li(+), Na(+) and NH(4)(+)). In addition, a correlation between metal ion charge density, binding affinity and enzyme activity was found among mono- and divalent metal ions except Pb(2+); higher charge density resulted in stronger affinity and higher activity, suggesting that the observed folding and activity is at least partially due to electrostatic interactions between ions and the DNAzyme. Finally, circular dichroism (CD) study has revealed Z-DNA formation with the monovalent metal ions, Zn(2+) and Mg(2+); the K(d) values obtained using CD were in the same range as those obtained from folding studies using FRET. However, Z-DNA formation was not observed with Pb(2+). These results indicate that Pb(2+)-dependent function follows a different mechanism from the monovalent metal ions and other divalent metal ions; in the presence of latter metal ions, metal-ion dependent folding and structural changes, including formation of Z-DNA, play an important role in the catalytic function of the 8-17 DNAzyme.

Thermal stability, structural features, and B-to-Z transition in DNA tetraloop hairpins as determined by optical spectroscopy in d(CG)(3)T(4)(CG)(3) and d(CG)(3)A(4)(CG)(3) oligodeoxynucleotides

NMR and CD data have previously shown the formation of the T(4) tetraloop hairpin in aqueous solutions, as well as the possibility of the B-to-Z transition in its stem in high salt concentration conditions. It has been shown that the stem B-to-Z transition in T(4) hairpins leads to S (south)- to N (north)-type conformational changes in the loop sugars, as well as anti to syn orientations in the loop bases. In this article, we have compared by means of UV absorption, CD, Raman, and Fourier transform infrared (FTIR), the thermodynamic and structural properties of the T(4) and A(4) tetraloop hairpins formed in 5'-d(CGCGCG-TTTT-CGCGCG)-3' and 5'-d(CGCGCG-AAAA-CGCGCG)-3', respectively. In presence of 5M NaClO(4), a complete B-to-Z transition of the stems is first proved by CD spectra. UV melting profiles are consistent with a higher thermal stability of the T(4) hairpin compared to the A(4) hairpin. Order-to-disorder transition of both hairpins has also been analyzed by means of Raman spectra recorded as a function of temperature. A clear Z-to-B transition of the stem has been confirmed in the T(4) hairpin, and not in the A(4) hairpin. With a right-handed stem, Raman and FTIR spectra have confirmed the C2'-endo/anti conformation for all the T(4) loop nucleosides. With a left-handed stem, a part of the T(4) loop sugars adopt a N-type (C3'-endo) conformation, and the C3'-endo/syn conformation seems to be the preferred one for the dA residues involved in the A(4) tetraloop.

Calculating sequence-dependent melting stability of duplex DNA oligomers and multiplex sequence analysis by graphs

The analytical methods for characterizing DNA sequence-dependent thermodynamic stability have been reviewed. A set of n-n sequence stability parameters is presented. Examples in which these values are used to calculate the thermodynamic stability of short duplex DNA oligomers are presented. The problem of determining sets of isothermal sequences is addressed by representing DNA sequences as graphs. Representing DNA sequences by a graph descriptor with special mathematical properties minimizes the computational difficulty of determining the number of DNA sequences with identical predicted thermodynamic stability. This is achieved by replacement of a whole set of sequences by a single representative. Applications of this concept were demonstrated for sequences assembled from individual bases and sequences assembled from oligomeric blocks.

Cisplatin (cis-Pt(NH(3))(2)Cl(2)) and cis-[Pt(NH(3))(2)(H(2)O)(2)](2+) Intrastrand Cross-Linking Reactions at the Telomere GGGT DNA Sequence Embedded in a Duplex, a Hairpin, and a Bulged Duplex: Use of Mg(2+) and Zn(2+) to Convert a Hairpin to a Bulged Duplex

In the past, we showed that metal species have a high affinity for the central G in the GGG sequence of the duplex d(A(1)T(2)G(3)G(4)G(5)T(6)A(7)C(8)C(9)C(10)A(11)T(12))(2) (G3-D) and that cisplatin (cis-Pt(NH(3))(2)Cl(2)) and G3-D formed an N7-Pt-N7 G(4),G(5) intrastrand cross-link preferentially over the G(3),G(4) adduct ( approximately 25:1). Thus, a putative G(4) monoadduct was postulated to cross-link in the 3'- rather than the normally more favorable 5'-direction. To evaluate this hypothesis and also to explore why the G3-D G(4),G(5) adduct had an unusual hairpin structure, we have now introduced the use of N,N'-dimethylthiourea (DMTU) as a monoadduct trap and have extended the study to a G3-D analogue with a hairpin form, d(A(1)T(2)G(3)G(4)G(5)T(6)T(7)C(8)C(9)C(10)A(11)T(12)) (G3-H). Chemical shift and 2D (1)H and (13)C NMR data indicated that the G3-H hairpin has a stem region with B-form structure and a nonhelical loop region. Zn(2+) or Mg(2+) ions transformed G3-H into a bulged duplex. Downfield shifts of G(4)H8 and G(4)C8 NMR signals indicated that Zn(2+) binds preferentially to G(4)N7. Reaction of cisplatin or cis-[Pt(NH(3))(2)(H(2)O)(2)](2+) with the bulged duplex and hairpin forms of G3-H gave a similar intrastrand cross-link ratio, G(4),G(5):G(3),G(4) = 7:3. This ratio is insensitive to DNA form or Pt leaving group. For G3-D this ratio is lower in the cis-[Pt(NH(3))(2)(H(2)O)(2)](2+) reaction ( approximately 1:1) than in the cisplatin reaction (25:1), indicating that the leaving group influences the cross-linking step for G3-D. The G(4) monoadducts of the cis-Pt(NH(3))(2)Cl(2)-G3-H and -G3-D and the cis-[Pt(NH(3))(2)(H(2)O)(2)](2+)-G3-D reactions were trapped with DMTU, but no monoadduct was trapped in the cis-[Pt(NH(3))(2)(H(2)O)(2)](2+)-G3-H reaction. The results suggest that the respective monoadducts are more long-lived for G3-D. We postulate that the G(5) in the G3-D Cl-G(4) monoadduct is placed in a favorable position to form the cross-link because of a prior conformational change induced by G(4)-A(7) stacking. This accounts for the very high selectivity for 3'-cross-linking. Nevertheless, in all other cases, regardless of the form or conformation, 3'-direction cross-linking is unusually favored at GGGT sequences, suggesting that the sequence itself contributes greatly to the 3'-cross-linking preference; since telomeres have multiple repeats of this GGGT sequence, this finding may have biological relevance.

A 6 bp Z-DNA hairpin binds two Z alpha domains from the human RNA editing enzyme ADAR1

The Z alpha domain of the human RNA editing enzyme double-stranded RNA deaminase I (ADAR1) binds to left-handed Z-DNA with high affinity. We found by analytical ultracentrifugation and CD spectroscopy that two Z alpha domains bind to one d(CG)3T4(CG)3 hairpin which contains a stem of six base pairs in the Z-DNA conformation. Both wild-type Z alpha and a C125S mutant show a mean dissociation constant of 30 nM as measured by surface plasmon resonance and analytical ultracentrifugation. Our data suggest that short (> or = 6 bp) segments of Z-DNA within a gene are able to recruit two ADAR1 enzymes to that particular site.

Studies of DNA dumbbells. IV. Preparation and melting of a DNA dumbbell with the 16 base-pair sequence 5'G-T-A-T-C-C-C-T-C-T-G-G-A-T-A-C3' linked on the ends by dodecyl chains

The preparation and melting of a 16 base-pair duplex DNA linked on both ends by C12H24 (dodecyl) chains is described. Absorbance vs temperature curves (optical melting curves) were measured for the dodecyl-linked molecule and the same duplex molecule linked on the ends instead by T4 loops. Optical melting curves of both molecules were measured in 25, 55, and 85 mM Na+ and revealed, regardless of [Na+], the duplex linked by dodecyl loops is more stable by at least 6 degrees C than the same duplex linked by T4 loops. Experimental curves in each salt environment were analyzed in terms of the two-state and multistate theoretical models. In the two-state, or van't Hoff analysis, the melting transition is assumed to occur in an all-or-none manner. Thus, the only possible states accessible to the molecule throughout the melting transition are the completely intact duplex and the completely melted duplex or minicircle. In the multistate analysis no assumptions regarding the melting transition are required and the statistical occurrence of every possible partially melted state of the duplex is explicitly considered. Results of the analysis revealed the melting transitions of both the dodecyl-linked molecule and the dumbbell with T4 end loops are essentially two state in 25 and 55 mM Na+. In contrast, significant deviations from two-state behavior were observed in 85 mM Na+. From our previously published melting data of DNA dumbbells with Tn end loops where n = 2, 3, 4, 6, 8, 10, 14 [T. M. Paner, M. Amaratunga, and A. S. Benight, (1992) Biopolymers, Vol. 32, pp. 881-892] and the dumbbell with T4 end loops of this study, a plot of d(Tm)/d ln [Na+] was constructed. Extrapolation of this data to n = 1 intersects with the value of d (Tm)/d ln [Na+] obtained for the alkyl-linked dumbbell, suggesting the salt-dependent stability of the alkyl-linked molecule behaves as though the duplex of this molecule were linked by end loops comprised of a single T residue.

Studies of DNA dumbbells. II. Construction and characterization of DNA dumbbells with a 16 base-pair duplex stem and Tn end loops (n = 2, 3, 4, 6, 8, 10, 14)

The preparation and characterization of DNA dumbbells that contain the 16 base-pair duplex sequences 5'G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C3' (set 1) and 5'G-C-A-T-C-A-T-C-G-A-T-G-A-T-G-C3' (set 2) are reported. The dumbbells of set 1 have the duplex stem nucleated on both ends by Tn (n = 2, 3, 4, 6, 8, 10, and 14) loops. The dumbbells of set 2 have Tn (n = 2, 4, 8, 10) end loops. For the molecules of set 1, effects of end loop size on the electrophoretic mobility, CD and UV absorbance spectra, and cleavage by restriction enzymes, were investigated. Effects of loop size on the CD spectra and restriction enzyme cleavage of the molecules of set 2 were also examined. Optical melting curves of the molecules of set 1 were collected as a function of sodium ion concentration from 30 to 120 mM. These investigations revealed that as loop size decreases, the electrophoretic mobilities, rates of enzyme cleavage, and optical melting temperatures increase. For end loops with at least three T's the observed increases are inversely proportional to loop size. The behavior of the dumbbell with T2 end loops departs from this linear dependence and is anomalous in every experimental context. For molecules with end loops comprised of at least four T's CD spectra were virtually indistinguishable. However, these spectra differed considerably from the CD spectrum of the T2-looped molecule. The CD spectrum of the dumbbell with T3 end loops displayed features common to the dumbbells with larger loops and T2 end loops. Thermodynamic evidence that the terminal G.C base pairs (bps) nucleating the T2 end loops were intact was obtained from a comparison of the melting temperature of this molecule with that of a DNA dumbbell containing the 14 central bps of the set 1 duplex sequence linked instead by end loops comprised of the four base sequence, C-T-T-C. The tm of this latter molecule was determined to be 9 degrees C less than that of the former dumbbell assumed to contain a 16-bp stem and T2 end loops.(ABSTRACT TRUNCATED AT 400 WORDS)

Coupling of sequential transitions in a DNA double hairpin: energetics, ion binding, and hydration

In an effort to evaluate the relative contributions of sequence, ion binding, and hydration to the thermodynamic stability of nucleic acids, we have investigated the melting behavior of a double hairpin and that of its component single hairpins. Temperature-dependent UV absorption and differential scanning calorimetry techniques have been used to characterize the helix-coil transitions of three deoxyoligonucleotides: d(GTACT5GTAC), d(GCGCT5GCGC), and d(GCGCT5GCGCGTACT5GTAC). The first two oligomers melt with transition temperatures equal to 28 and 69 degrees C, respectively, in 10 mM dibasic sodium phosphate at pH 7.0. The Tm's are independent of strand concentration, strongly indicating the presence of single-stranded hairpin structures at low temperatures. The third oligomer, with a sequence corresponding to the joined sequences of the first two oligomers, melts with two apparently independent monomolecular transitions with Tm's of 41 and 69 degrees C. These transitions correspond to the melting of a double hairpin. In the salt range of 10-100 mM in NaCl, we obtain average enthalpies of 24 and 38 kcal/mol for the transitions in the single-hairpin molecules. Each transition in the double hairpin has an enthalpy of 32 kcal/mol. In addition, dtm/d log [Na+] for the transitions are 4.1 and 4.7 degrees C for the single hairpins and 12.6 and 11.2 degrees C for each transition in the double hairpin. The differential ion binding parameter between the double hairpin and that of the sum of single hairpins is roughly equal to 1.1 mol of Na+ ions/mol of double hairpin and is consistent with an increase in the electrostatic behavior of the stem phosphates of this molecule.

DNA hairpin loops in solution. Correlation between primary structure, thermostability and reactivity with single-strand-specific nuclease from mung bean

Hairpin structures formed by seven DNA inverted repeats have been studied by PAGE, UV(CD)-spectroscopy and nuclease cleavage. The hairpins consisted of (CG)3 stems and loops of 2, 3 and 4 residues. Thermal stabilities (Tm) have been determined in low and high ionic strength buffers, where the hairpins were structured in the B- and Z-DNA form respectively. The thermodynamic parameters of hairpin formation have been obtained by a two-state analysis of the hairpin-coil transitions. It is found that, on increasing the number of bases in the loop from 2 to 3 and 4, the Tms of the B-hairpins decrease, whereas the Tms of the same hairpins in the Z-form increase. This confirms previous evidence (1,2) that in a hairpin molecule the size and structure of the loop are modulated by the conformation of the helical stem. Moreover, B-hairpins with loops comprising 2, 3 and 4 bases have been digested with the single-strand-specific nuclease from mung bean. In our experimental conditions (0 degrees C) the nuclease preferentially cleaves the unbonded nucleotides of the loops. However, the rates of loop hydrolysis, which roughly follow a first-order kinetics, markedly depend on the size of the loop. At a ratio of 3 enzyme units/micrograms DNA, the half-lives of hairpins which are expected to form loops of 4, 3 and 2 residues are 90, 145 and 440 minutes respectively. Thermostability and enzymatic digestion data suggest that two-membered loops can be formed in B-hairpins but not in Z-hairpins.

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)

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

We report Potential of Mean Force studies to describe the relative thermodynamic stabilities of d(GCCGCAGC) in a mismatched duplex and a hairpin monomer conformation in NaCl solution. The PMF calculations are combined with previous molecular mechanics and normal mode analysis in order to estimate the role of different components of the free energy in determining the relative stability of the duplex and hairpin structures. The high entropy associated with the loop region and the lack of minor groove phosphate-phosphate interactions in the hairpin compete against the gain in enthalpic contribution to the free energy due to base pairing in the mismatched duplex. The combined free energy calculations show that the hairpin is the most stable conformation at low salt and that a hairpin to duplex transition takes place at approximately 0.47 M NaCl. In addition, we studied the hairpin to partially stacked single helical conformation equilibrium at low salt. We found a small variation in transition temperature in salt concentration, delta Tm/delta log10(cs) approximately 2-3 degrees K/decade, in contrast to the duplex to hairpin or duplex to partially stacked single helix transition where the transition temperature exhibited marked dependence on salt concentration. This is in qualitative agreement with experimental data. Based on the Potential of Mean Force free energy calculation, the order of relative stability of the three-conformations studied varies with salt concentration. We observed the following orders of stability: stacked single helix greater than hairpin greater than duplex for cs less than 0.77 M NaCl; single helix greater than duplex greater than hairpin for 0.77 less than Cs less than 2.1 M; and duplex greater than hairpin greater than single strand for cs greater than 2.1 M. From the calculated PMF free energy curves in the NaCl concentration range, 0.012 less than cs less than 5.0 M, we can assign upper and lower bounds for the non-ionic differences in free energy between the duplex, hairpin, and stacked single helical states (at standard conditions: cs = 1.0 M, T = 25 degrees C, and 1 M oligomer concentration). We found that for delta G duplex single helix = G duplex - 2 x G single helix less than -7.38 Kcal/mol, the single helix is the least stable state. For the duplex-to-hairpin free energy difference in the range, -1.87 less than delta G duplex-hairpin less than 0.03 Kcal/mol, there will always be a salt-induced hairpin-to-duplex transition for 0.01 less than cs less than 1.6 M NaCl. If delta G duplex-hairpin less than -1.87, the duplex is always more stable than the hairpin; and for delta G duplex-hairpin greater than Kcal/mol, the hairpin state is always more stable than the duplex, for all salt concentrations.

B to Z transitions of the short DNA hairpins formed from the oligomer sequences: d[(CG)3X4(CG)3] (X = A, T, G, C)

Circular Dichroism (CD) spectra were collected as a function of sodium perchlorate concentration [NaClO4] for the set of DNA hairpins formed from the oligomer sequences d[(CG)3X4(CG)3] where X = A, T, G or C. Over the range in salt concentration from 0 to 4.0 M NaClO4, the CD spectra invert in a manner characteristic of the B to Z transition. A factor analysis routine is described and employed to determine the least number of basis spectra required to fit the measured spectra of each hairpin over the entire salt range examined. In every case, linear combinations of only two sub-spectra fit the experimental spectra of the hairpins with greater than 98% accuracy, indicating the spectrally monitored structural transitions are two-state. From the relative weights of the individual sub-spectra, B-Z transition curves are constructed. The transitions are analyzed in terms of a simple two-state equilibrium model which yields an evaluation of the transition free-energy, delta GB-Z, as a function of NaClO4 concentration. At 1.0 M NaClO4 and 21 degrees C, delta GB-Z = 5.4, 4.9, 3.6 and 2.3 kcal/mole for the G4, T4, A4 and C4 loop hairpins, respectively.

Structure of a T4 hairpin loop on a Z-DNA stem and comparison with A-RNA and B-DNA loops

The synthetic DNA oligomer C-G-C-G-C-G-T-T-T-T-C-G-C-G-C-G crystallizes as a Z-DNA hexamer, capped at one end by a T4 loop. The crystals are monoclinic, space group C2, with a = 57.18 A, b = 21.63 A, c = 36.40 A, beta = 95.22 degrees, and one hairpin molecule per asymmetric unit. The structure of the z-hexamer stem was determined by molecular replacement, and the T4 loop was positioned by difference map methods. The final R factor at 2.1 A resolution for hairpin plus 70 water molecules is 20% for 2 sigma data, with a root-mean-square error of 0.26 A. The (C-G)3 stem resembles the free Z-DNA hexamer with minor crystal packing effects. The T4 loop differs from that observed on a B-DNA stem in solution, or in longer loops in tRNA, in that it shows intraloop and intermolecular interactions rather than base stacking on the final base-pair of the stem. Bases T7, T8 and T9 stack with one another and with the sugar of T7. Two T10 bases from different molecules stack between the C1-G12 terminal base-pairs of a third and fourth molecule, to simulate a T.T "base-pair". Distances between thymine N and O atoms suggest that the two thymine bases are hydrogen bonded, and a keto-enol tautomer pair is favored over disordered keto-keto wobble pairs. The hairpin molecules pack in the crystal in herringbone columns in a manner that accounts well for the observed relative crystal growth rates in a, b and c directions. Hydration seems to be most extensive around the phosphate groups, with lesser hydration within the grooves.

Thermodynamic stability of the 5' dangling-ended DNA hairpins formed from sequences 5'-(XY)2GGATAC(T)4GTATCC-3', where X, Y = A, T, G, C

Expressions for the partition function Q (T) of DNA hairpins are presented. Calculations of Q (T), in conjunction with our previously reported numerically exact algorithm [T. M. Paner, M. Amaratunga, M. J. Doktycz, and A. S. Benight (1990) Biopolymers, 29, 1715-1734], yield a numerical method to evaluate the temperature dependence of the transition enthalpy, entropy, and free energy of a DNA hairpin directly from its optical melting curve. No prior assumptions that the short hairpins melt in a two-state manner are required. This method is then applied in a systematic manner to investigate the stability of the six basepair duplex stem 5'-GGATAC-3' having four-base dangling single-strand ends with the sequences (XY)2, where X, Y = A, T, G, C, on the 5' end and a T4 loop on the 3' end. Results show that all dangling ends of the sample set stabilize the hairpin against melting. Increases in transition temperatures as great as 4.0 degrees C above the blunt-ended control hairpin were observed. The hierarchy of the hairpin transition temperatures is dictated by the identity of the first base of the dangling end adjoining the duplex in the order: purine greater than T greater than C. Calculated melting curves of every hairpin were fit to experimental curves by adjustment of a single parameter in the numerically exact theoretical algorithm. Exact fits were obtained in all cases. Experimental melting curves were also calculated assuming a two-state melting process. Equally accurate fits of all dangling-ended hairpin melting curves were obtained with the two-state model calculation. This was not the case for the melting curve of the blunt-ended hairpin, indicating the presence of a four-base dangling-end drives hairpin melting to a two-state process. Q (T) was calculated as a function of temperature for each hairpin using the theoretical parameters that provided calculated curves in exact agreement with the experimentally obtained optical melting curves. From Q (T), the temperature dependence of the transition enthalpy delta H, entropy delta S, and free energy delta G were calculated for every hairpin providing a quantitative assessment of the effects of dangling ends on hairpin thermodynamics. Comparisons of our results are made with those of the Breslauer group [M. Senior, R. A. Jones, and K. J. Breslauer (1988) Biochemistry 27, 3879-3885] on the T2 5' dangling-ended d(GC)3 duplexes.(ABSTRACT TRUNCATED AT 400 WORDS)