Before revision:  Computing the melting profiles, I discuss the relation between nonlinear character of the base pair interactions and twisting geometry .
After revision:  Computing the equilibrium thermodynamics, we focus on the interplay between twisting of the complementary strands around the molecule axis and nonlinear stacking potential: it is shown that the latter affects the melting profiles only if the rotational degrees of freedom are included in the Hamiltonian. The use of ladder Hamiltonian models for the DNA complementary strands in the pre-melting regime is questioned .