Correcting for heat capacity and 5'-TA type terminal nearest neighbors improves prediction of DNA melting temperatures using nearest-neighbor thermodynamic models

Nearest-neighbor thermodynamic (NNT) models currently provide some of the most accurate predictions of melting thermodynamics, including melting temperature (T(m)) values, for short DNA duplexes. Inherent to all existing NNT models is the assumption that ΔH° and ΔS° for the helix-to-coil transition are temperature invariant. Here we investigate the impact that this zero-ΔC(p) assumption has on the accuracy of T(m) predictions for 128 DNA duplexes. Previous and new melting thermodynamic data are analyzed to establish an estimate of ΔC(p)(bp), the heat capacity change per base pair, of 42 ± 16 cal mol(-1) K(-1) bp(-1), as well as an optimal thermodynamic reference temperature (T(ref)) of 53 ± 5 °C. These results were used to modify the unified NNT model to properly account for the temperature dependence of ΔH° and ΔS° and thereby extend the range over which T(m) is accurately predicted. This new approach is shown to be especially useful for duplexes that melt at a T(m) greater than 70 °C. Thermodynamic data collected by differential scanning calorimetry (DSC) for 16 duplexes designed to melt over a broad temperature range were used to verify the values of ΔC(p)(bp) and T(ref) and to show that ΔC(p)(bp) is essentially constant above 37 °C. Additional DSC analysis of 12 duplex sequences containing all 10 nearest neighbors allowed for errors associated with different terminal nearest neighbors to be examined and showed that duplexes containing one or more terminal 5'-TA groups are significantly more stable than predicted by the unified NNT model. A correction to improve prediction of the hybridization thermodynamics of duplexes with terminal 5'-TA groups is provided.     

Effect of locked nucleic acid modifications on the thermal stability of noncanonical DNA structure

We studied the kinetic and thermodynamic effects of locked nucleic acid (LNA) modifications on parallel and antiparallel DNA duplexes. The LNA modifications were introduced at cytosine bases of the pyrimidine strand. Kinetic parameters evaluated from melting and annealing curves showed that the association and dissociation rate constants for the formation of the LNA-modified parallel duplex at 25.0 °C were 3 orders of magnitude larger and 6 orders of magnitude smaller, respectively, than that of the unmodified parallel duplex. The activation energy evaluated from the temperature-dependent rate constants was largely altered by the LNA modifications, suggesting that the LNA modifications affected a prenucleation event in the folding process. Moreover, thermodynamic parameters showed that the extent of stabilization by the LNA modification for parallel duplexes (3.6 kcal mol(-1) per one modification) was much more significant than that of antiparallel duplexes (1.6 kcal mol(-1)). This large stabilization was due to the decrease in ΔH° that was more favorable than the decrease in TΔS°. These quantitative parameters demonstrated that LNA modification specifically stabilized the noncanonical parallel duplex. On the basis of these observations, we succeeded to stabilize the parallel duplex by LNA modification at the physiological pH. These results can be useful in the rational design of functional molecules such as more effective antisense and antigene strands, more sensitive strands for detection of target DNA and RNA strands, and molecular switches responding to solution pH.     

Synthesis and studies on the effect of 2-thiouridine and 4-thiouridine on sugar conformation and RNA duplex stability.

In order to understand the effect of 2-thiouridine (s2U) substitution on RNA structure and the potential for stabilization of tRNA codon-anticodon interactions through s2U-34 modification, a pentamer RNA sequence, Gs2UUUC, was synthesized and characterized by NMR spectroscopy. The single strand contains the UUU anticodon sequence of tRNALys with flanking GCs to increase duplex stability. Regiochemical effects of uridine thiolation were determined by comparing the structure and stability of the 2-thiouridine containing oligonucleotide with an identical sequence containing 4-thiouridine (s4U) and also the normal uridine nucleoside. Circular dichroism spectrum indicated an A-form helical conformation for Gs2UUUC which was further confirmed by 2D ROESY NMR experiments. The duplex stability of the three pentamers complexed with a 2'-O-methyl-ribonucleotide complementary strand, GmAmAmAmCm, was determined by UV thermal melting studies and by 1H NMR spectroscopy. The duplex containing s2U has a T m of 30.7 degrees C compared to 19. 0 degrees C for the unmodified control and 14.5 degrees C for the s4U containing duplex. The results from UV experiments were corroborated by imino proton NMR studies that show proton exchange rates, chemical shift differences, and NH proton linewidths indicative of the stability order s2U >U >s4U. The magnitude of the effect of s2U in our model system is comparable to the 20 degrees C stabilization observed by Grosjean and co-workers for 2-thiolation in a codon-anticodon model system composed of two tRNAs with complementary anticodon sequences [Houssier, C., Degee, P., Nicoghosian, K. and Grosjean, H. (1988) J. Biomol. Struct. Dyn., 5, 1259-1266].     

Solution structure of an LNA hybridized to DNA: NMR study of the d(CT(L)GCT(L)T(L)CT(L)GC):d(GCAGAAGCAG) duplex containing four locked nucleotides

We have used two-dimensional (1)H NMR spectroscopy at 750 MHz to determine a high-resolution solution structure of an oligonucleotide containing restricted nucleotides with a 2'-O, 4'-C-methylene bridge (LNA) hybridized to the complementary DNA strand. The LNA:DNA duplex examined contained four thymidine LNA modifications (T(L), d(C1T(L)2G3C4T(L)5T(L)6C7T(L)8G9C10):d( G11C12A13G14A15A16G17C 18A19G20). A total relaxation matrix approach was used to obtain interproton distance bounds from NOESY cross-peak intensities. These distance bounds were used as restraints in molecular dynamics (rMD) calculations. Forty final structures were generated for the duplex from A-form and B-form DNA starting structures. The root-mean-square deviation (RMSD) of the coordinates for the 40 structures of the complex was 0.6 A. The sugar puckerings are averaged values of a dynamic interchange between N- and S-type conformation except in case of the locked nucleotides that were found to be fixed in the C3'-endo conformation. Among the other nucleotides in the modified strand, the furanose ring of C7 and G9 is predominantly in the N-type conformation whereas that of G3 is in a mixed conformation. The furanose rings of the nucleotides in the unmodified complementary strand are almost exclusively in the S-type conformation. Due to these different conformations of the sugars in the two strands, there is a structural strain between the A-type modified strand and the B-type unmodified complementary strand. This strain is relaxed by decreasing the value of rise and compensating with tip, buckle, and propeller twist. The values of twist vary along the strand but for a majority of the base pairs a value even lower than that of A-DNA is observed. The average twist over the sequence is 32+/-1 degrees. On the basis of the structure, we conclude that the high stability of LNA:DNA duplexes is caused by a local change of the phosphate backbone geometry that favors a higher degree of stacking.     

Sequence-dependent thermodynamic parameters for locked nucleic acid (LNA)-DNA duplex formation

The design of modified nucleic acid probes, primers, and therapeutics is improved by considering their thermodynamics. Locked nucleic acid (LNA) is one of the most useful modified backbones, with incorporation of a single LNA providing a substantial increase in duplex stability. In this work, the hybridization DeltaH(o), DeltaS(o), and melting temperature (T(M)) were measured from absorbance melting curves for 100 duplex oligonucleotides with single internal LNA nucleotides on one strand, and the results provided DeltaDeltaH(o), DeltaDeltaS(o), DeltaDelta, and DeltaT(M) relative to reference DNA oligonucleotides. LNA pyrimidines contribute more stability than purines, especially A(L), but there is substantial context dependence for each LNA base. Both the 5' and 3' neighbors must be considered in predicting the effect of an LNA incorporation, with purine neighbors providing more stability. Enthalpy-entropy compensation in DeltaDeltaH(o) and DeltaDeltaS(o) is observed across the set of sequences, suggesting that LNA can stabilize the duplex by either preorganization or improved stacking, but not both simultaneously. Singular value decomposition analysis provides predictive sequence-dependent rules for hybridization of singly LNA-substituted DNA oligonucleotides to their all-DNA complements. The results are provided as sets of DeltaDeltaH(o), DeltaDeltaS(o), and DeltaDelta parameters for all 32 of the possible nearest neighbors for LNA+DNA:DNA hybridization (5' MX(L) and 5' X(L)N, where M, N, and X = A, C, G, or T and X(L) represents LNA). The parameters are applicable within the standard thermodynamic prediction algorithms. They provide T(M) estimates accurate to within 2 degrees C for LNA-containing oligonucleotides, which is significantly better accuracy than previously available.     

Characterization of oligonucleotides with LNA-monomers for PCR detection of point mutations in <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis> genome

Point mutations associated with isoniazid resistance in Mycobacterium tuberculosis (MTB) have been analyzed in codon 315 of the katG gene by conventional polymerase chain reaction (PCR) using primers containing locked nucleic acid (LNA) modified nucleotides. Purity and structure of primers containing 5 LNA monomers of 17 nucleotides in length were characterized by matrix assisted laser desorption ionization time of flight mass spectrometry (MALDI-TOF MS) and a 17-mer duplex formed by two complementary oligonucleotides was characterized by the method of thermal denaturation. The duplex containing five LNA monomers per each strand was characterized by a higher melting temperature than it was expected using extrapolation of theoretical calculation for nucleotide modification of one strand of the duplex. Detection of any of six possible mutations in katG codon 315 (i.e. discrimination between sensitive and resistant MTB) requires just one PCR employing a set of two primers with one LNA-modified primer; this is an important advantage of oligonucleotides containing LNA over unmodified nucleotides: employment of multiplex PCR would require up to 12 primers. Problems of control of oligonucleotide modification by LNA monomers are discussed.     

Evidence for a DNA triplex in a recombination-like motif: I. Recognition of Watson-Crick base pairs by natural bases in a high-stability triplex

Data are presented on a triplex type with two parallel homologous strands for which triplex formation is almost as strong as duplex formation at least for some sequences and even at pH 7 and 0.2 M NaCl. The evidence mainly rests upon comparing thermodynamic properties of similar systems. A paperclip oligonucleotide d(A12C4T12C4A12) with two linkers C4 obviously can form a triplex with parallel back-folded adenine strand regions, because the single melting transition of this complex splits in two transitions by introducing mismatches only in the third strand region. Respectively, a hairpin duplex d(A12C4T12) and a single strand d(A12) form a triplex as a 1:1 complex in which the second adenine strand is parallel oriented to the homologous one in the Watson-Crick paired duplex. In this system the melting temperature T(m) of the triplex is practically the same as that of the duplex d(A12)-d(T12), at least within a complex concentration range of 0.2-4.0 microM. The melting behaviour of complexes between triplex stabilizing ligand BePI and the system hairpin duplex plus single strand supports the triplex model. Non-denaturing gel electrophoresis suggests the existence of a triplex for a system in which five of the twelve A-T*A base triads are substituted by C-G*C base triads. The recognition between any substituted Watson-Crick base pair (X-Y) in the hairpin duplex d(A4XA7C4T7YT4) and the correspondingly replaced base (Z) in the third strand d(A4ZA7) is mutually selective. All triplexes with matching base substitutions (Z = X) have nearly the same stability (T(m) values from 29 to 33.5 degrees C), whereas triplexes with non-matching substitutions (Z not equal X) show a clearly reduced stability (T(m) values from 15 to 22 degrees C) at 2microM equimolar oligonucleotide concentration. Most nucleic acid triple helices hitherto known are limited to homopurine-homopyrimidine sequences in the target duplex. A stable triplex formation is demonstrated for inhomogeneous sequences tolerating at least 50% pyrimidine content in the homologous strands. On the basis of the surprisingly similar thermodynamic parameters for duplex and triplex, and of the fact that this triplex type seems to be more stable than many other natural DNA triplexes known, and on the basis of semiempirical and molecule mechanical calculations, we postulate bridging interactions of the third strand with the two other strands in the triplex according to the recombination motif. This triplex, denoted by us 'recombination-like form', tolerates heterogeneous base sequences.     

Thermodynamic, counterion, and hydration effects for the incorporation of locked nucleic acid nucleotides into DNA duplexes

A locked nucleic acid (LNA) monomer is a conformationally restricted nucleotide analogue with an extra 2'-O, 4'-C-methylene bridge added to the ribose ring. LNA-modified oligonucleotides are known to exhibit enhanced hybridization affinity toward complementary DNA and RNA. In this work, we have evaluated the hybridization thermodynamics of a series of LNA-substituted DNA octamers, modified to various extents by one to three LNA substitutions, introduced at either adenine (5'-AGCACCAG) or thymine (5'-TGCTCCTG) nucleotides. To understand the energetics, counterion effects, and the hydration contribution of the incorporation of LNA modification, a combination of spectroscopic and calorimetric techniques was used. The CD spectra of the corresponding duplexes showed that the modified duplexes adopt an A-type conformation. UV and DSC melting studies revealed that each type of duplex unfolds in a two-state transition. A complete thermodynamic profile at 5 degrees C indicated that the net effect of modification on thermodynamic parameters might be positional and that the neighboring bases flanking the modification might influence the favorable formation of the modified duplexes. Furthermore, relative to the formation of the unmodified reference duplexes, the formation of modified duplexes is accompanied by a higher uptake of counterions and a lower uptake of water molecules.     

Site-specifically located 8-amino-2'-deoxyguanosine: thermodynamic stability and mutagenic properties in Escherichia coli

2-Nitropropane (2-NP), an important industrial solvent and a component of cigarette smoke, is mutagenic in bacteria and carcinogenic in rats. 8-Amino-2'-deoxyguanosine (8-amino-dG) is one of the types of DNA damage found in liver, the target organ in 2-NP-treated rats. To investigate the thermodynamic properties of 8-amino-dG opposite each of the four DNA bases, we have synthesized an 11mer, d(CCATCG*CTACC), in which G* represents the modified base. By annealing a complementary DNA strand to this modified 11mer, four sets of duplexes were generated each containing one of the four DNA bases opposite the lesion. Circular dichroism studies indicated that 8-amino-dG did not alter the global helical properties of natural right-handed B-DNA. The thermal stability of each duplex was examined by UV melting measurements and compared with its unmodified counterpart. For the unmodified 11mer, the relative stability of the complementary DNA bases opposite G was in the order C > T > G > A, as determined from their -DeltaG degrees values. The free energy change of each modified duplex was lower than its unmodified counterpart, except for the G*:G pair that exhibited a higher melting transition and a larger -DeltaG degrees than the G:G duplex. Nevertheless, the stability of the modified 11mer duplex also followed the order C > T > G > A when placed opposite 8-amino-dG. To explore if 8-amino-dG opposite another 8-amino-dG has any advantage in base pairing, a G*:G* duplex was evaluated, which showed that the stability of this duplex was similar to the G*:G duplex. Mutagenesis of 8-amino-dG in this sequence context was studied in Escherichia coli, which showed that the lesion is weakly mutagenic (mutation frequency approximately 10(-3)) but still can induce a variety of targeted and semi-targeted mutations.     

Thermodynamics of DNA-RNA heteroduplex formation: effects of locked nucleic acid nucleotides incorporated into the DNA strand

A locked nucleic acid (LNA) monomer is a conformationally restricted nucleotide analogue exhibiting enhanced hybridization efficiency toward complementary strand. The potential of LNA-based oligonucleotides has been sought to improve the selectivity and specificity of probe sets employed in detection and specific targeting of nucleic acids. We have evaluated the influence of "locked nucleic acid" residues on hybridization thermodynamics, counterions and hydration of DNA.RNA heteroduplex using spectroscopic and calorimetric techniques. One to three LNA substitutions have been introduced either at the adenine (5'-AGCACCAG) or thymine (5'-TGCTCCTG) residues of the DNA strand. A complete thermodynamic profile for heteroduplex formation suggested that LNA-induced stabilization results from a favorable increase in the enthalpy of hybridization that compensates for the unfavorable entropy change. Analysis of differential scanning calorimetry data indicated a nonzero heat capacity change, DeltaCp, accompanying the heteroduplex formation. Isothermal titration calorimetry measurements indicated an increase in binding affinity of the two strands as the LNA content of the heteroduplex is increased. Overall our result demonstrated that the effect of LNA-substitution at the thymine residue is more pronounced compared to the adenine residue. Furthermore, optical melting studies showed that, compared to an unmodified duplex, the formation of LNA-modified duplex is accompanied by a higher uptake of counterions and a lower uptake of water molecules. Our result, thus, presents a preliminary attempt toward the characterization of hybridization thermodynamics of the LNA-based probe-target sets, which will in turn aid in the selection of optimal conditions for hybridization experiments, and evaluation of the minimum probe-length required for hybridization and cloning experiments.     

Stability and mismatch discrimination of locked nucleic acid-DNA duplexes

Locked nucleic acids (LNA; symbols of bases, +A, +C, +G, and +T) are introduced into chemically synthesized oligonucleotides to increase duplex stability and specificity. To understand these effects, we have determined thermodynamic parameters of consecutive LNA nucleotides. We present guidelines for the design of LNA oligonucleotides and introduce free online software that predicts the stability of any LNA duplex oligomer. Thermodynamic analysis shows that the single strand-duplex transition is characterized by a favorable enthalpic change and by an unfavorable loss of entropy. A single LNA modification confines the local conformation of nucleotides, causing a smaller, less unfavorable entropic loss when the single strand is restricted to the rigid duplex structure. Additional LNAs adjacent to the initial modification appear to enhance stacking and H-bonding interactions because they increase the enthalpic contributions to duplex stabilization. New nearest-neighbor parameters correctly forecast the positive and negative effects of LNAs on mismatch discrimination. Specificity is enhanced in a majority of sequences and is dependent on mismatch type and adjacent base pairs; the largest discriminatory boost occurs for the central +C·C mismatch within the +T+C+C sequence and the +A·G mismatch within the +T+A+G sequence. LNAs do not affect specificity in some sequences and even impair it for many +G·T and +C·A mismatches. The level of mismatch discrimination decreases the most for the central +G·T mismatch within the +G+G+C sequence and the +C·A mismatch within the +G+C+G sequence. We hypothesize that these discrimination changes are not unique features of LNAs but originate from the shift of the duplex conformation from B-form to A-form.     

Stabilizing effect of propionic acid derivative of anthraquinone--polyamine conjugate incorporated into α-β chimeric oligonucleotides on the alternate-stranded triple helix

Two types of anthraquinone conjugates were synthesized as non-nucleosidic oligonucleotide components. These include an anthraquinone derivative conjugated with 2,2-bis(hydroxymethyl)propionic acid and an anthraquinone--polyamine derivative conjugated with 2,2-bis(hydroxymethyl)propionic acid. The conjugates were successfully incorporated into the "linking-region" of the α-β chimeric oligonucleotides via phosphoramidite method as non-nucleosidic backbone units. The resultant novel α-β chimeric oligonucleotides possessed two diastereomers that were generated by the introduction of the anthraquinone conjugate with a stereogenic carbon atom. The isomers were successfully separated by a reversed-phase HPLC. UV-melting experiments revealed that both stereoisomers formed a substantially stable alternate-strand triple helix, irrespective of the stereochemistry of the incorporated non-nucleosidic backbone unit. However, the enhancing effect on thermal stability depended on the length of the alkyl linker connecting anthraquinone moiety and the propionic acid moiety. The sequence discrimination ability of the chimeric oligonucleotides toward mismatch target duplex was also examined. The T(m) values of the triplexes containing the mismatch target were substantially lower than the T(m) values of those containing the full-match target. The thermodynamic parameters (ΔH°, ΔS°, and ΔG°) required for the dissociation of the triplexes into the third strand and target duplex were also measured.     

Solution structure of a dsDNA:LNA triplex

We have determined the NMR structure of an intramolecular dsDNA:LNA triplex, where the LNA strand is composed of alternating LNA and DNA nucleotides. The LNA oligonucleotide binds to the dsDNA duplex in the major groove by formation of Hoogsteen hydrogen bonds to the purine strand of the duplex. The structure of the dsDNA duplex is changed to accommodate the LNA strand, and it adopts a geometry intermediate between A- and B-type. There is a substantial propeller twist between base-paired nucleobases. This propeller twist and a concomitant large propeller twist between the purine and LNA strands allows the pyrimidines of the LNA strand to interact with the 5′-flanking duplex pyrimidines. Altogether, the triplex has a regular global geometry as shown by a straight helix axis. This shows that even though the third strand is composed of alternating DNA and LNA monomers with different sugar puckers, it forms a seamless triplex. The thermostability of the triplex is increased by 19°C relative to the unmodified DNA triplex at acidic pH. Using NMR spectroscopy, we show that the dsDNA:LNA triplex is stable at pH 8, and that the triplex structure is identical to the structure determined at pH 5.1.     

The solution structure of a locked nucleic acid (LNA) hybridized to DNA

LNA (Locked Nucleic Acids) is a novel oligonucleotide analogue containing a conformationally restricted nucleotide with a 2'-O, 4'-C-methylene bridge that induces unprecedented thermal affinities when mixed with complementary single stranded DNA and RNA. We have used two-dimensional 1H NMR spectroscopy obtained at 750 and 500 MHz to determine a high resolution solution structure of an LNA oligonucleotide hybridized to the complementary DNA strand. The determination of the structure was based on a complete relaxation matrix analysis of the NOESY cross peaks followed by restrained molecular dynamics calculations. Forty final structures were generated for the duplex from A-type and B-type dsDNA starting structures. The root-mean-square deviation (RMSD) of the coordinates for the forty structures of the complex was 0.32A. The structures were analysed by use of calculated helix parameters. This showed that the values for rise and buckle in the LNA duplex is markedly different from canonical B-DNA at the modification site. A value of twist similar to A-DNA is also observed at the modification site. The overall length of the helix which is 27.3 A. The average twist over the sequence are 35.9 degrees +/- 0.3 degrees. Consequently, the modification does not cause the helix to unwind. The bis-intercalation of the thiazole orange dye TOTO to the LNA duplex was also investigated by 1H NMR spectroscopy to sense the structural change from the unmodified oligonucleotide. We observed that the bis-intercalation of TOTO is much less favourable in the 5'-CT(L)AG-3' site than in the unmodified 5'-CTAG-3' site. This was related to the change in the base stacking of the LNA duplex compared to the unmodified duplex.     

Properties of 2'-fluorothymidine-containing oligonucleotides: interaction with restriction endonuclease EcoRV

2'-Fluorothymidine (Tf) was synthesized via an improved procedure with (diethylamino)sulfur trifluoride. The compatibility of the analogue with DNA synthesis via the phosphoramidite method was demonstrated after complete enzymatic digestion of the oligonucleotides d(Tf11T) and d(Tf3T), the sole products of which were 2'-fluorothymidine and thymidine in the expected ratio. The 2'-fluorothymidine was also incorporated into the EcoRV recognition sequence (underlined), within the complementary oligonucleotides d(CAAACCGATATCGTTGTG) and d(CACAACGATATCGGTTTG). Thermal melting characteristics of these duplexes showed a significant decrease in stability only when both of the thymidine residues in one of the strands were replaced. In contrast, when all of one strand of a duplex contained 2'-fluorothymidine, as in d(Tf11T).d(A12), a substantially higher Tm and cooperativity of melting was observed relative to the unmodified structure. EcoRV cleaved a duplex that contained a 2'-fluorothymidine at the scissile linkage in each strand at two-thirds of the rate obtained for the unmodified structure. A duplex containing two 2'-fluorothymidine residues in one strand and none in the other was cleaved at one-third of the rate in its unsubstituted strand, whereas the cleavage rate was reduced to 22% in its modified strand.     

Cytarabine-induced destabilization of a model Okazaki fragment.

Cytarabine is a potent anticancer drug that interferes with elongation of the lagging strand at the replication fork during DNA synthesis. The effects of cytarabine substitution on the structural and thermodynamic properties of a model Okazaki fragment were investigated using UV hyperchromicity and 1H NMR spectroscopy to determine how cytarabine alters the physicochemical properties of Okazaki fragments that are intermediates during DNA replication. Two model Okazaki fragments were prepared corresponding to a primary initiation site for DNA replication in the SV40 viral genome. One model Okazaki fragment consisted of five ribo- and seven deoxyribonucleotides on the hybrid strand, together with its complementary (DNA) strand. The second model Okazaki fragment was identical to the first with the exception of cytarabine substitution for deoxycytidine at the third DNA nucleotide of the hybrid strand. Thermodynamic parameters for the duplex to single strand transition for each model Okazaki fragment were calculated from the concentration dependence of the T m at 260 nm. Cytarabine significantly decreased the stability of this model Okazaki fragment, decreasing the melting temperature from 46.8 to 42.4 degrees C at a concentration of 1.33 x 10(-5) M. The free energy for the duplex to single strand transition was 1.2 kcal/mol less favorable for the cytarabine-substituted Okazaki fragment relative to the control at 37 degrees C. Analysis of the temperature dependence of the imino1H resonances for the two duplexes demonstrated that cytarabine specifically destabilized the DNA:DNA duplex portion of the model Okazaki fragment. These results are consistent with inhibition of lagging strand DNA synthesis by cytarabine substitution resulting from destabilization of the DNA:DNA duplex portion of Okazaki fragments in vivo .     

Melting studies of short DNA hairpins: Influence of loop sequence and adjoining base pair identity on hairpin thermodynamic stability

Spectroscopic and calorimetric melting studies of 28 DNA hairpins were performed. These hairpins form by intramolecular folding of 16 base self‐complementary DNA oligomer sequences. Sequence design dictated that the hairpin structures have a six base pair duplex linked by a four base loop and that the first five base pairs in the stem are the same in every molecule. Only loop sequence and identity of the duplex base pair closing the loop vary for the set of hairpins. For these DNA samples, melting studies were carried out to investigate effects of the variables on hairpin stability. Stability of the 28 oligomers was ascertained from their temperature‐induced melting transitions in buffered 115 mM Na⁺ solvent, monitored by ultraviolet absorbance and differential scanning calorimetry (DSC). Experiments revealed the melting temperatures of these molecules range from 32.4 to 60.5°C and are concentration independent over strand concentrations of 0.5 to 260 μM; thus, as expected for hairpins, the melting transitions are apparently unimolecular. Model independent thermodynamic transition parameters, ΔH_cal, ΔS_cal, and ΔG_cal, were determined from DSC measurements. Model dependent transition parameters, ΔH_vH, ΔS_vH, and ΔG_vH were estimated from a van't Hoff (two‐state) analysis of optical melting transitions. Results of these studies reveal a significant sequence dependence to DNA hairpin stability. Thermodynamic parameters evaluated by either procedure reveal the transition enthalpy, ΔH_cal (ΔH_vH) can differ by as much as 20 kcal/mol depending on sequence. Similarly, values of the transition entropy ΔS_cal (ΔS_vH) can differ by as much as 60 cal/Kmol (eu) for different molecules. Differences in free energies ΔG_cal (ΔG_vH) are as large as 4 kcal/mol for hairpins with different sequences. Comparisons between the model independent calorimetric values and the thermodynamic parameters evaluated assuming a two‐state model reveal that 10 of the 28 hairpins display non‐two‐state melting behavior. The database of sequence‐dependent melting free energies obtained for the hairpins was employed to extract a set of n‐n (nearest‐neighbor) sequence dependent loop parameters that were able to reproduce the input data within error (with only two exceptions). Surprisingly, this suggests that the thermodynamic stability of the DNA hairpins can in large part be reasonably represented in terms of sums of appropriate nearest‐neighbor loop sequence parameters. © 1999 John Wiley & Sons, Inc. Biopoly 50: 425–442, 1999     

Studies of DNA dumbbells. V. A DNA triplex formed between a 28 base-pair DNA dumbbell substrate and a 16 base linear single strand

CD spectra and melting curves were collected for a 28 base-pair DNA fragment in the form of a DNA dumbbell (linked on both ends by T₄ single-strand loops) and the same DNA sequence in the linear form (without end loops). The central 16 base pairs (bp) of the 28-bp duplex region is the poly(pu) sequence: 5′-AGGAAGGAGGAAAGAG-3′. Mixtures of the dumbbell and linear DNAs with the 16-base single-strand sequence 5′-TCCTTCCTCCTTTCTC-3′ were also prepared and studied. At 22°C, CD measurements of the mixtures in 950 mM NaCl, 10 mM sodium acetate, 1 mM EDTA, pH 5.5, at a duplex concentration of 1.8 μM, provided evidence for triplex formation. Spectroscopic features of the triplexes formed with either a dumbbell or linear substrate were quite similar. Melting curves of the duplex molecules alone and in mixtures with the third strand were collected as a function of duplex concentration from 0.16 to 2.15 μM. Melting curves of the dumbbell alone and mixtures with the third strand were entirely independent of DNA concentration. In contrast, melting curves of the linear duplex alone or mixed with the third strand were concentration dependent. At identical duplex concentrations, the dumbbell alone melts ～20°C higher than the linear duplex. The curve of the linear duplex displayed a significant pretransition probably due to end fraying. On melting curves of mixtures of the dumbbell or linear duplex with the third strand, a low temperature transition with much lower relative hyperchromicity change (～ 5%) was observed. This transition was attributed to the melting of a new molecular species, e.g., the triplex formed between the duplex and single-strand DNA molecules. In the case of the dumbbell/single-strand mixture, these melting transitions of the triplex and the dumbbell were entirely resolvable. In contrast, the melting transitions of the linear duplex and the triplex overlapped, thereby preventing their clear distinction. To analyze the data, a three-state equilibrium model is presented. The analysis utilizes differences in relative absorbance vs temperature curves of dumbbells (or linear molecules) alone and in mixtures with the third strand. From the model analysis a straightforward derivation of f_T(T), the fraction of triplex as a function of temperature, was obtained. Analysis of f_T vs temperature curves, in effect melting curves of the triplexes, provided evaluation of thermodynamic parameters of the melting transition. For the triplex formed with the dumbbell substrate, the total transition enthalpy is ΔH_T = 118.4 ± 12.8 kcal/mol (7.4 ± 0.8 kcal/mol per triplet unit) and the total transition entropy is ΔS_T = 344 ± 36.8 cal/K · mol (eu) (21.5 ± 2.3 eu per triple unit). The transition curves of the triplex formed with the linear duplex substrate displayed two distinct regions. A broad pretransition region from f_T = 0 to 0.55 and a higher, sharper transition above f_T = 0.55. The transition parameters derived from the lower temperature region of the curve are ΔH′_T = 44.8 ± 9.6 kcal/mol and ΔS′_T = 112 ± 33.6 eu (or ΔH′ = 2.8 ± 0.6 kcal/mol and ΔS′ = 7.0 ± 2.1 eu per triplet). These values are probably too small to correspond to actual melting of the triplex but instead likely reveal effects of end fraying of the duplex substrate on triplex stability. Transition parameters of the upper transition are ΔH′_T = 128.0 ± 2.3 kcal/mol and ΔS′_T = 379.2 ± 6.4 eu (ΔH′ = 8.0 ± 0.2 kcal/mol and ΔS′ = 23.7 ± 0.4 eu per triplet) in good agreement (within experimental error) with the transition parameters of the triplex formed with the dumbbell substrate. Supposing this upper transition reflects actual dissociation of the third strand from the linear duplex substrate this triplex is comparable in thermodynamic stability to the triplex formed with a dumbbell substrate. Even so, the biphasic melting character of the linear triplex obscures the whole analysis, casting doubt on its absolute reliability. Apparently triplexes formed with a dumbbell substrate offer technical advantages over triplexes formed from linear or hairpin duplex substrates for studies of DNA triplex stability. © 1993 John Wiley & Sons, Inc.     

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°C, respectively, and exhibit similar thermodynamic stabilities. The enthalpies of duplex formation, measured by nmr, were found to be very similar (ΔH_DH = 26 ± 3 kcal/mole) for both decanters at low salt concentrations (< 50 mM NaCl). However, as the salt concentration was increased the behavior of ΔH_DH, 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(CAAC-CCGTTG) 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 ΔH_DH, dependence on sail 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. © 1995 John Wiley & Sons, Inc.     

Impact of the oxidized guanine lesion spiroiminodihydantoin on the conformation and thermodynamic stability of a 15-mer DNA duplex

Spiroiminodihydantoin (Sp) is a hyperoxidized guanine base produced from oxidation of the mutagenic DNA lesion 7,8-dihydro-8-oxo-2'-deoxguanosine (8-oxoG) by a variety of species including peroxynitrite, singlet oxygen, and the high-valent metals Ir(IV) and Cr(V). In this study, the conformation and thermodynamic stability of a 15-mer DNA duplex containing an Sp lesion are examined using spectroscopic techniques and differential scanning calorimetry (DSC). The Sp lesion does not alter the global B-form conformation of the DNA duplex as determined by circular dichroism spectroscopy. Thermal denaturation experiments find that Sp significantly lowers the thermal stability of the duplex by approximately 20 degrees C. The enthalpies, entropies, and free energies of duplex formation for 15-mers containing guanine, 8-oxoG, and Sp were determined by performing DSC experiments as well as van't Hoff analysis of UV melting spectroscopic data. The thermodynamic stability of the Sp duplex is significantly reduced compared to that of both the 8-oxoG and parent G duplexes, with the thermodynamic destabilization being enthalpic in origin. The thermodynamic impact of the Sp lesion is compared to what is found for other types of DNA base damage and discussed in relation to how the presence of this lesion could affect cellular processes, in particular the recognition and repair of these adducts by the base excision repair enzymes.