Thermodynamics of DNA containing very stable chemically modified base pairs

DNA chemical modifications caused by the binding of some antitumor drugs give rise to a very strong local stabilization of the double helix. These sites melt at a temperature that is well above the melting temperatures of ordinary AT and GC base pairs. In this work we have examined the melting behavior of DNA containing very stable sites. Analytical expressions were derived and used to evaluate the thermodynamic properties of homopolymer DNA with several different distributions of stable sites. The results were extended to DNA with a heterogeneous sequence of AT and GC base pairs. The results were compared to the melting properties of DNA with ordinary covalent interstrand cross-links. It was found that, as with an ordinary interstrand cross-link, a single strongly stabilized site makes a DNA's melting temperature (T(m)) independent of strand concentration. However in contrast to a DNA with an interstrand cross-link, a strongly stabilized site makes the DNA's T(m) independent of DNA length and equal to T(infinity), the melting temperature of an infinite length DNA with the same GC-content and without a stabilized site. Moreover, at a temperature where more than 80% of base pairs are melted, the number of ordinary (non-modified) helical base pairs (n) is independent of both the DNA length and the location of the stabilized sites. For this condition, n(T) = (2 omega-a)S/(1-S) and S = expDeltaS(T(infinity)-T)/(RT)] where omega is the number of strongly stabilized sites in the DNA chain, a is the number of DNA ends that contain a stabilized site, and DeltaS, T, and R are the base pair entropy change, the temperature, and the universal gas constant per mole. The above expression is valid for a temperature interval that corresponds to n<0.2N for omega=1, and n<0.1N for omega>1, where N is the number of ordinary base pairs in the DNA chain.