Correlation of Tm and sequence of DNA duplexes with ΔH computed by an improved empirical potential method

Abstract

T_m values for 20 DNA duplexes with different repeating base sequences provided the data base for developing a rational and relatively simple methodology for computing apparent enthalpies for the helix → coil transitions of DNA helices, ΔH_calc. The computational variables and their range of acceptable values were selected on the basis of physically plausible arguments. Over 350,000 different combinations of the variables were tested for degree of fit. It was therby possible to find a combination giving a high degree of linear fit between T_m and ΔH_calc (correlation coefficient, 0.99), with T_m values deviating (on average) from the regression line by only ±2.17°C. Most of this uncertainty is attributed to experimental limitations, although computational approximations also contribute. With ΔH_calc for the melting of each of the unique complementary dinucleotide fragments computed by the method developed, it is possible to estimate T_m and (relative) ΔH_calc reliable for the melting of any particular DNA [with base pairs G(I)·C and A·T] given only its base sequence. The ΔH_calc values for the complementary dinucleotide fragments, together with statistical considerations, make it apparent that T_m of DNAs with repeating base sequence show only an approximate linear dependence on G·C content because A·T and G·G pairs do not contribute to helix stability independently of the base‐pair sequence in which they occur. In fact, the nearestneighbor stacking interactions are so significant that certain complementary dinucleotide fragment sequences with 0,50, and 100% G·C content have the same stability.