We propose a statistical mechanics model for the melting transition of DNA. Base pairing and stacking are treated as separate degrees of freedom, and the interplay between pairing and stacking is described by a set of local rules which mimic the geometrical constraints in the real molecule. This microscopic mechanism intrinsically accounts for the cooperativity related to the free energy penalty of bubble nucleation. The model describes both the unpairing and unstacking parts of the spectroscopically determined experimental melting curves. Furthermore, the model explains the observed temperature dependence of the effective thermodynamic parameters used in models of the nearest neighbor type. We compute the partition function for the model through the transfer matrix formalism, which we also generalize to include nonlocal chain entropy terms.

• Received 13 September 2004

DOI:https://doi.org/10.1103/PhysRevE.71.041909