The four-state model of a linear lattice: application to ligand-controlled hybridization of short duplex DNAs

An algorithm for calculating the order-disorder transition of a four-state linear lattice is presented. Recursion schemes for the probabilities of each site along the lattice occupying one of the four possible states were derived following the work of D. Poland (Biopolymers, 1974, Vol. 13, pp. 1859-1871) and H. Tachibana and A. Wada (Biopolymers, 1982, Vol. 21, pp. 1873-1885). The algorithm was parameterized to consider melting of short duplex DNA in the presence of duplex and single-strand binding ligands. Model calculations were performed for two 31 base-pair duplex DNAs with very different percent guanine-cytosine base pairs, and thus very different thermodynamic stabilities. In the absence of ligands, calculated melting curves of the two DNAs under identical solvent conditions and identical concentrations were separated by over 15 degrees C. In the presence of either duplex or single-strand binding ligands, if a sequence dependence to ligand binding is assumed, the melting curves of the two DNAs can be made to coalesce, i.e., stability of the two DNAs can be normalized! This example demonstrates the feasibility of controlling hybridization of short DNAs by sequence specific ligand binding.     

DNA strand exchange and selective DNA annealing promoted by the human immunodeficiency virus type 1 nucleocapsid protein.

Nucleocapsid protein (NC) of human immunodeficiency virus type 1 (HIV-1) was expressed in Escherichia coli and purified. The protein displayed a variety of activities on DNA structure, all reflecting an ability to promote transition between double-helical and single-stranded conformations. We found that, in addition to its previously described ability to accelerate renaturation of complementary DNA strands, the HIV-1 NC protein could substantially lower the melting temperature of duplex DNA and could promote strand exchange between double-stranded and single-stranded DNA molecules. Moreover, in the presence of HIV-1 NC, annealing of a single-stranded DNA molecule to a complementary DNA strand that would yield a more stable double-stranded product was favored over annealing to alternative complementary DNA strands that would form less stable duplex products (selective annealing). NC thus appears to lower the kinetic barrier so that double-strand <==> single-strand equilibrium is rapidly reached to favor the lowest free-energy nucleic acid conformation. This activity of NC may be important for correct folding of viral genomic RNA and may have practical applications.     

Inverted repeat sequences can influence the melting transitions of linear DNAs

The influence of inverted repeat sequences on the melting transitions of linear DNAs has been examined. Derivative melting curves (DMC) of a 514 base pair (bp) DNA, seven subfragments of this DNA, and four other DNAs have been compared to predictions of DNA melting theory. The 514-bp DNA contains three inverted repeat sequences that can form cruciform structures in supercoiled DNA. We refer to these sequences as c-inverted repeats. Previous work showed that the DMC of this DNA, unlike a number of other DNAs, is not accurately predicted by DNA melting theory. Since the theoretical model does not include hairpin-like structures, it was suggested that hairpin or cruciform formation in these inverted repeats may be responsible for this discrepancy. Our results support this hypothesis. Predicted DMCs are in good agreement with DNAs with no inverted repeats, or inverted repeats not evident in supercoiled DNA. Differences between the theoretical and experimental Tm's are less than or equal to 0.3 degrees C. DNA molecules that contain one or more of the three c-inverted repeats are not as accurately predicted. Experimental Tm values are lower than predicted values by 0.7-3.8 degrees C. It is concluded that some inverted repeat sequences can form hairpin-like structures during the melting of linear DNAs. These structures appear to lower overall DNA stability.     

Proton exchange and base pair opening in a DNA triple helix

Nuclear magnetic resonance spectroscopy has been used to characterize opening reactions and stabilities of individual base pairs in two related DNA structures. The first is the triplex structure formed by the DNA 31-mer 5'-AGAGAGAACCCCTTCTCTCTTTTTCTCTCTT-3'. The structure belongs to the YRY (or parallel) family of triple helices. The second structure is the hairpin double helix formed by the DNA 20-mer 5'-AGAGAGAACCCCTTCTCTCT-3' and corresponds to the duplex part of the YRY triplex. The rates of exchange of imino protons with solvent in the two structures have been measured by magnetization transfer from water and by real-time exchange at 10 degrees C in 100 mM NaCl and 5 mM MgCl2 at pH 5.5 and in the presence of two exchange catalysts. The results indicate that the exchange of imino protons in protonated cytosines is most likely limited by the opening of Hoogsteen C+G base pairs. The base pair opening parameters estimated from imino proton exchange rates suggest that the stability of individual Hoogsteen base pairs in the DNA triplex is comparable to that of Watson-Crick base pairs in double-helical DNA. In the triplex structure, the exchange rates of imino protons in Watson-Crick base pairs are up to 5000-fold lower than those in double-helical DNA. This result suggests that formation of the triplex structure enhances the stability of Watson-Crick base pairs by up to 5 kcal/mol. This stabilization depends on the specific location of each triad in the triplex structure.     

Melting of the solvent structure around a RNA duplex: a molecular dynamics simulation study

From three 2.4-ns molecular dynamics simulations of the r(CpG)(12) duplex conducted at 5, 25 and 37 degrees C, a strong temperature dependence of the dynamics of the water molecules and ions located in the first nucleic acid coordination shell is observed. At 5 degrees C, the highest residence times of bound water molecules exceed 1 ns while, at 37 degrees C, they decrease to 0.5 ns in agreement with available NMR data. Similar temperature dependencies are observed for the potassium ions bound to the duplex. In this temperature range, the structure of the RNA helix remains essentially unchanged. Thus, the observed alterations correspond to a 'premelting' of the solvent structure around the duplex. It is proposed that, before the nucleic acid structure melts, the entropy of the solvent increases to a point where it is no longer compensated by the enthalpic contribution of solute-solute and solute-solvent interactions. At this stage, the weakest structural elements start to melt. In other terms, the experimentally observed melting processes are preceded by a melting of the more labile solvent structure.     

Structural characteristics of DNA of Bacillus thuringiensis phages

The spectral characteristics of DNA from two phages of the polylysogenic Bacillus thuringiensis subsp. galleriae culture 1-97 were studied. The typical parameters of melting and the negative reaction with formaldehyde are indicative of the double-helical structure of these DNAs. The phage DNAs differ in the molar content of nitrogen bases (32 and 38 mole% of GC) and in their distribution along the molecule. This distribution is uniform in the DNA of one phage whereas the other phage DNA is composed of heterological segments with a different nucleotide composition.     

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in four 40 base pair deoxyoligonucleotides

Effects of different end sequences on melting, circular dichroism spectra (CD), and enzyme binding properties were investigated for four 40 base pair, non-self-complementary duplex DNA oligomers. The center sequences of these oligoduplexes have either of two 22 base pair modules flanked on both sides by sequences differing in AT content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from a van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming that the melting transition is two-state. Melting free energies (20 degrees C) evaluated from DSC melting experiments on the four duplex DNAs ranged from -52.2 to -77.5 kcal/mol. Free energies based on the van't Hoff analysis were -37.9 to -58.8 kcal/mol. Although the values are different, trends in the melting free energies of the four duplex DNAs as a function of sequence were identical in both DSC and optical analyses. Subject to several assumptions, values for the initiation free energy were estimated for each duplex, defined as DeltaG(int) = DeltaG(cal) - DeltaG(pred), where DeltaG(cal) is the experimental free energy at 20 degrees C determined from the experimentially measured values of the transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal). The predicted free energy of the sequence, DeltaG(pred)(20 degrees C), is obtained using published nearest-neighbor sequence stability values. For three of the four duplexes, values of DeltaG(int) are essentially nil. In contrast, the duplex with 81.8% GC has a considerably higher estimate of DeltaG(int) = 7.1 kcal/mol. The CD spectra for the six duplexes collected over the wavelength range from 200 to 320 nm are also sequence-dependent. Factor analysis of the CD spectra by singular value decomposition revealed that the experimental CD spectra could be reconstructed from linear combinations of two minor and one major subspectra. Changes in the coefficients of the major subspectrum for different sequences reflect incremental sequence-dependent variations of the CD spectra. Equilibrium binding by BamHI restriction endonuclease to the 40 base pair DNAs whose central eight base pairs contain the recognition sequence for BamHI restriction enzyme bounded by A.T base pairs, 5'-A-GGATCC-A-3' was investigated. Binding assays were performed by titering BamHI against a constant concentration of each of the duplex DNA substrates, in the absence of Mg(2+), followed by analysis by gel retardation. Under the conditions employed, the enzyme binds but does not cleave the DNAs. Results of the assays revealed two binding modes with retarded gel mobilities. Binding isotherms for the fraction of bound DNA species versus enzyme concentration for each binding mode were constructed and analyzed with a simple two-step equilibrium binding model. This analysis provided semiquantitative estimates on the equilibrium binding constants for BamHI to the four DNAs. Values obtained for the binding constants varied only 7-fold and ranged from 6 x 10(-)(8) to 42 x 10(-)(8) M, with binding free energies from -8.6 to -9.7 (+/- 0.2) kcal/mol depending on the sequence that flanks the enzyme binding site. Unlike what was found earlier in binding studies of the 22 base pair duplexes that constitute the core modules of the present 40-mers [Riccelli, P. V., Vallone, P. M., Kashin, I., Faldasz, B. D., Lane, M. J., and Benight, A. S. (1999) Biochemistry 38, 11197-11208], no obvious relationship between binding and stability was found for these longer DNAs. Apparently, effects of flanking sequence stability on restriction enzyme binding may only be measurable in very short duplex deoxyoligonucl     

Hybrid hexanucleotide duplex containing cyclonucleotides and deoxynucleotides: the d(TA) segment can adopt a high anti left-handed double-helical structure

It is known that oligonucleotides containing cyclonucleosides with a high anti (intermediate between anti and syn) glycosidic conformation adopt left-handed, single- and double-helical structures [Uesugi, S., Yano, J., Yano, E., & Ikehara, M. (1977) J. Am. Chem. Soc. 99, 2313-2323]. In order to see whether DNA can adopt the high anti left-handed double-helical structure or not, a self-complementary hexanucleotide containing 6,2'-O-cyclocytidine (C(o)), 8,2'-O-cycloguanosine (G(o)), thymidine, and deoxyadenosine, C(o)G(o)dTdAC(o)G(o), was synthesized. Imino proton NMR spectra and the results of nuclear Overhauser effect experiments strongly suggest that C(o)G(o)dTdAC(o)G(o) adopts a left-handed double-helical structure where the deoxynucleoside residues are involved in hydrogen bonding and take a high anti glycosidic conformation. A conformational model of the left-handed duplex was obtained by calculation with energy minimization. Thus it appears that DNA can form a high anti, left-handed double helix under some constrained conditions, which is quite different from that of Z-DNA.     

N(4)C-ethyl-N(4)C cross-linked DNA: synthesis and characterization of duplexes with interstrand cross-links of different orientations.

The preparation and physical properties of short DNA duplexes that contain a N(4)C-ethyl-N(4)C interstrand cross-link are described. Duplexes that contain an interstrand cross-link between mismatched C-C residues and duplexes in which the C residues of a -CG- or -GC- step are linked to give "staggered" interstrand cross-links were prepared using a novel N(4)C-ethyl-N(4)C phosphoramidite reagent. Duplexes with the C-C mismatch cross-link have UV thermal transition temperatures that are 25 degrees C higher than the melting temperatures of control duplexes in which the cross-link is replaced with a G-C base pair. It appears that this cross-link stabilizes adjacent base pairs and does not perturb the structure of the helix, a conclusion that is supported by the CD spectrum of this duplex and by molecular models. An even higher level of stabilization, 49 degrees C, is seen with the duplex that contains a -CG- staggered cross-link. Molecular models suggest that this cross-link may induce propeller twisting in the cross-linked base pairs, and the CD spectrum of this duplex exhibits an unusual negative band at 298 nm, although the remainder of the spectrum is similar to that of B-form DNA. Mismatched C-C or -CG- staggered cross-linked duplexes that have complementary overhanging ends can undergo self-ligation catalyzed by T4 DNA ligase. Analysis of the ligated oligomers by nondenaturing polyacrylamide gel electrophoresis shows that the resulting oligomers migrate in a manner similar to that of a mixture of non-cross-linked control oligomers and suggests that these cross-links do not result in significant bending of the helix. However, the orientation of the staggered cross-link can have a significant effect on the structure and stability of the cross-linked duplex. Thus, the thermal stability of the duplex that contains a -GC- staggered cross-link is 10 degrees C lower than the melting temperature of the control, non-cross-linked duplex. Unlike the -CG- staggered cross-link, in which the cross-linked base pairs can still maintain hydrogen bond contacts, molecular models suggest that formation of the -GC- staggered cross-link disrupts hydrogen bonding and may also perturb adjacent base pairs leading to an overall reduction in helix stability. Duplexes with specifically positioned and oriented cross-links can be used as substrates to study DNA repair mechanisms.     

The HU protein from Thermotoga maritima: recombinant expression, purification and physicochemical characterization of an extremely hyperthermophilic DNA-binding protein.

The histone-like protein TmHU from the hyperthermophilic eubacterium Thermotoga maritima was cloned, expressed to high levels in Escherichia coli, and purified to homogeneity by heat precipitation and cation exchange chromatography. CD spectroscopical studies with secondary structure analysis as well as comparative modeling demonstrate that the dimeric TmHU has a tertiary structure similar to other homologous HU proteins. The Tm of the protein was determined to be 96 degrees C, and thermal unfolding is nearly completely reversible. Surface plasmon resonance measurements for TmHU show that the protein binds to DNA in a highly cooperative manner, with a KD of 73 nM and a Hill coefficient of 7.6 for a 56 bp DNA fragment. It is demonstrated that TmHU is capable to increase the melting point of a synthetic, double-stranded DNA (poly[d(A-T)]) by 47 degrees C, thus suggesting that DNA stabilization may be a major function of this protein in hyperthermophiles. The significant in vitro protection of double-helical DNA may be useful for biotechnological applications.     

Transforming Region of Group A, B, and C Adenoviruses: DNA Homology Studies with Twenty-Nine Human Adenovirus Serotypes

The 31 human adenovirus (Ad) serotypes form five groups based upon DNA genome homologies: group A (Ad12, 18, 31), group B (Ad3, 7, 11, 14, 16, 21), group C (Ad1, 2, 5, 6), group D (Ad8, 9, 10, 13, 15, 17, 19, 20, 22-30), and group E (Ad4) (M. Green, J. Mackey, W. Wold, and P. Rigden, Virology, in press). Group A Ads are highly oncogenic in newborn hamsters, group B Ads are weakly oncogenic, and other Ads are nononcogenic. However, most or all Ads transform cultured cells. We have studied the homology of Ad5, Ad7, and Ad12 transforming restriction endonuclease DNA fragments with DNAs of 29 Ad types. Ad5 HindIII-G (map position 0-7.3), Ad7 XhoI-C (map position 0-10.8), and Ad12 (strain Huie) EcoRI-C (map position 0-16) and SalI-C (map position 0-10.6) fragments were purified, labeled in vitro (nick translation), and annealed with DNAs of Ad1 to Ad16, Ad18 to Ad24, and Ad26 to Ad31. Hybrids were assayed by using hydroxylapatite. Ad5 HindIII-G hybridized 98 to 100% with DNAs of group C Ads, but only 1 to 15% with DNAs of other types. Ad7 XhoI-C fragment hybridized 85 to 99% with DNAs of group B Ads, but only 6 to 21% with DNAs of other types. Ad12 (Huie) EcoRI-C hybridized 53 to 68% with DNAs of five other Ad12 strains, 53% with Ad18 DNA, 56% with Ad31 DNA, but only 3 to 13% with DNAs of other types. In vitro-labeled Ad12 (Huie) SalI-C hybridized 35 to 71% with DNAs of 6 other Ad12 strains, 44% with Ad18 DNA, 52% with Ad31 DNA, but only 2 to 7% with DNAs Ad7, Ad2, Ad26, or Ad4. When assayed using S-1 nuclease, SalI-C annealed 17 to 44% with DNAs of group A Ads. The melting temperatures of the hybrids of Ad5 HindIII-G with all group C Ad DNAs were 84 degrees C in 0.12 M sodium phosphate (pH 6.8). The melting temperature of the Ad12 (Huie) EcoRI-C hybrid with Ad12 (Huie) DNA was 83 degrees C, but was only 71 to 77 degrees C with DNAs of other group A Ads. Thus, group C and group B Ads both have very homologous transforming regions that are not represented in DNAs of non-group C Ads or non-group B Ads, respectively. Similarily, group A Ads have unique but less homologous transforming regions. These different transforming nucleotide sequences may be reflected in the different oncogenic properties of group A, B, and C Ads.     

Effects of sodium ions on DNA duplex oligomers: improved predictions of melting temperatures

Melting temperatures, T(m), were systematically studied for a set of 92 DNA duplex oligomers in a variety of sodium ion concentrations ranging from 69 mM to 1.02 M. The relationship between T(m) and ln [Na(+)] was nonlinear over this range of sodium ion concentrations, and the observed melting temperatures were poorly predicted by existing algorithms. A new empirical relationship was derived from UV melting data that employs a quadratic function, which better models the melting temperatures of DNA duplex oligomers as sodium ion concentration is varied. Statistical analysis shows that this improved salt correction is significantly more accurate than previously suggested algorithms and predicts salt-corrected melting temperatures with an average error of only 1.6 degrees C when tested against an independent validation set of T(m) measurements obtained from the literature. Differential scanning calorimetry studies demonstrate that this T(m) salt correction is insensitive to DNA concentration. The T(m) salt correction function was found to be sequence-dependent and varied with the fraction of G.C base pairs, in agreement with previous studies of genomic and polymeric DNAs. The salt correction function is independent of oligomer length, suggesting that end-fraying and other end effects have little influence on the amount of sodium counterions released during duplex melting. The results are discussed in the context of counterion condensation theory.     

Characterization of ciprofloxacin binding to the linear single- and double-stranded DNA

The binding of ciprofloxacin to natural and synthetic polymeric DNAs was investigated at different solvent conditions using a combination of spectroscopic and hydrodynamic techniques. In 10 mM cacodylate buffer (pH 7.0) containing 108.6 mM Na(+), no sequence preferences in the interaction of ciprofloxacin with DNA was detected, while in 2 mM cacodylate buffer (pH 7.0) containing only 1.7 mM Na(+), a significant binding of ciprofloxacin to natural and synthetic linear double-stranded DNA was observed. At low ionic strength of solution, ciprofloxacin binding to DNA duplex containing alternating AT base pairs is accompanied by the largest enhancement in thermal stability (e.g. DeltaT(m) approximately 10 degrees C for poly[d(AT)].poly[d(AT)]), and the most pronounced red shift in the position of the maximum of the fluorescence emission spectrum (lambda(max)). Similar red shift in the position of lambda(max) is also observed for ciprofloxacin binding to dodecameric duplex containing five successive alternating AT base pairs in the row. On the other hand, ciprofloxacin binding to poly[d(GC)].poly[d(GC)], calf thymus DNA and dodecameric duplex containing a mixed sequence is accompanied by the largest fluorescence intensity quenching. Addition of NaCl does not completely displace ciprofloxacin bound to DNA, indicating the binding is not entirely electrostatic in origin. The intrinsic viscosity data suggest some degree of ciprofloxacin intercalation into duplex.     

Solvent effects on the stability of A7U7p

The thermodynamics of double-helix formation were measured spectrophotometrically for A7U7 in water at 1 M NaCl and for A7U7p in a variety of solvent mixtures and salt. Comparison of the A7U7 results with calorimetric measurements indicates duplex formation involves intermediate states. For A7U7p between 0.06 and 0.55 M Na+, dTm/d(log [Na+]) = 17.4 degrees C, similar to the value of 19.6 degrees C for poly-(A).poly(U) [Krakauer, H., & Sturtevant, J. M. (1968) Biopolymers 6, 491-512]. At 1 M NaCl, the A7U7p duplex is most stable in 100% water. For 10 mol % solutions, the order for A7U7p duplex stability is ethylene glycol greater than glycerol greater than ethanol greater than 2-propanol greater than dimethyl sulfoxide greater than 1-propanol greater than formamide greater than N,N-dimethylformamide greater than urea greater than dioxane. Comparison of changes in stability and thermodynamic parameters with literature results for proteins suggests proteins and A7U7p interact differently with solvent. The results suggest hydrophobic bonding is not a major contributor to the stability of the A7U7p duplex. Comparisons with bulk solvent surface tension suggest the energy of cavity formation is also not a major contributor to duplex stability.     

Specific DNA binding by the homeodomain Nkx2.5(C56S): detection of impaired DNA or unfolded protein by isothermal titration calorimetry

Titrations of specific 18-bp duplex DNA with the cardiac-specific homeodomain Nkx2.5(C56S) have utilized an ultrasensitive isothermal titration calorimeter (ITC). As the free DNA nears depletion, we observe large apparent decreases in the binding enthalpy when the DNA is impaired or when the temperature is sufficiently high to produce some unfolding of the free protein. Either effect can be attributed to refolding of the biopolymer that occurs as a result of stabilization due to the large favorable change in free energy on the homeodomain binding to DNA (-49.4 kJ/mol at 298 K). In either case, thermodynamic parameters obtained in such ITC experiments are unreliable. By using a lower temperature (85 vs. 95 degrees C) during the annealing of complementary DNA strands, damage of the 18-bp duplex DNA (T(m) = 72 degrees C) is avoided, and titrations with the homeodomain are normal at temperatures from 10 to 40 degrees C when >95% of the protein is folded. Under the latter conditions, the heat capacity plot is linear with a DeltaC(p) value of -0.80 +/- 0.03 kJ K(-1) mol(-1), which is more negative than that calculated from the burial of solvent accessible surface areas (-0.64 +/- 0.05 kJ K(-1) mol(-1)), consistent with water structures being at the protein-DNA interfaces.     

Flexibility and curvature of duplex DNA containing mismatched sites as a function of temperature

The flexibility and curvature of duplex DNAs containing mismatched sites have been monitored as a function of temperature. The diffusion coefficients are dependent on the flexibility and the curvature of the DNA, and these have been determined by NMR-based methods. The diffusion coefficients, D, depend on a Boltzmann term and the viscosity of the solvent, eta, which is also temperature dependent. To analyze the temperature dependence of the diffusion results, the shape function, S(f) = etaD/T, is used. The shape functions do not have the viscosity and temperature dependence of the diffusion coefficients. The presence of mismatched sites significantly enhances the shape function of duplex DNA at all temperatures examined. The observed increases in the shape functions are attributed to the mismatched sites acting as localized flexible joints. The results on the temperature dependence of the shape functions, the optical absorbance, and the proton chemical shifts indicate that local melting at, and adjacent to, mismatched site occurs at a lower temperature than the overall melting of the duplexes. The localized melting gives rise to a considerable increase in the shape function. The contribution of the curvature of the mismatched sites to the enhanced diffusion has been examined. A DNA with mismatches that are in phase with respect to the helical repeat and a DNA which has the mismatches out of phase with respect to the helical repeat have been examined. The results indicate that mismatched sites have modest curvature.     

Sequences in the HSP90 promoter form G-quadruplex structures with selectivity for disubstituted phenyl bis-oxazole derivatives

The HSP90 protein is an important target in cancer. We report here that stable quadruplex DNAs can be formed from a promoter sequence in the HSP90 gene, on the basis of melting, circular and NMR studies, and show that these can be selectively targeted by non-macrocyclic quadruplex-stabilizing phenyl bis-oxazole derivatives. These do not bind significantly to duplex DNA and show low stabilization of the human telomeric quadruplex. These results suggest an approach to targeting HSP90 at the DNA level.     

Salt-concentration dependence of melting profiles of lambda phage DNAs: evidence for long-range interactions and pronounced end effects.

Melting profiles of DNAs from wild-type λ phage and a deletion mutant phage λb2 were examined in a wide range of salt concentration. The fine structure of the melting profiles changed sharply with salt concentration, especially in the range [Na⁺] ? 10 mM. A comparison of the melting profiles between the wild-type and the deletion mutant DNAs provided good evidence for extremely high melting cooperativity under low salt conditions, which is clearly manifested as the long-range interactions and the pronounced end effects; a large melting peak appeared as a result of the b2 deletion without any inserted sequence in the salt range [Na⁺] ? 2.8 mM. It was also suggested that in the further reduced salt range [Na⁺] ? 2.0 mM, melting of a λ DNA molecule starts from its right end rather than the most (A + T)-rich central region. The molecular basis of the high melting cooperativity at low salt concentrations can be explained in terms of the increased free energy associated with loop formation in the double-helical structure of DNA.