Exact algorithms for the calculation of melting curves of heterogeneous DNA with N base pairs apparently require computer time proportional to N2. However, it is shown that a decomposition of the loop entropy factor into a sum of I exponential functions (1) gives an extremely accurate approximation to the loop entropy factor for small values of I and (2) makes the computer time for the exact algorithms proportional to I·N. In effect, exact results for melting curves and lengths of helix or coil stretches are obtained with computer time comparable to that required for the Frank‐Kamenetskii approximation. The remarkable accuracy of the latter for the fraction of helical content (errors of 0.01–0.05) is confirmed, but appreciably larger errors are found for the lengths of helix or coil stretches (typical errors of 30–100%).
[1]
D. Crothers.
Calculation of melting curves for DNA
,
1968,
Biopolymers.
[2]
M. Frank-Kamenetskii,et al.
Comparison of different theoretical descriptions of helix–coil transitions in DNA
,
1976
.
[3]
D. Poland,et al.
Recursion relation generation of probability profiles for specific‐sequence macromolecules with long‐range correlations
,
1974,
Biopolymers.
[4]
M. Fixman.
Helix–coil transition in heterogeneous DNA
,
1975,
Biopolymers.
[5]
Y. Lyubchenko,et al.
Fine structure of DNA melting curves
,
1976,
Biopolymers.
[6]
Douglas Poland,et al.
Theory of helix-coil transitions in biopolymers
,
1970
.