On the origin of the temperature dependence of the supercoiling free energy

Monte Carlo simulations using temperature-invariant torsional and bending rigidities fail to predict the rather steep decline of the experimental supercoiling free energy with increasing temperature, and consequently fail to predict the correct sign and magnitude of the supercoiling entropy. To illustrate this problem, values of the twist energy parameter (E(T)), which governs the supercoiling free energy, were simulated using temperature-invariant torsion and bending potentials and compared to experimental data on pBR322 over a range of temperatures. The slope, -dE(T)/dT, of the simulated values is also compared to the slope derived from previous calorimetric data. The possibility that the discrepancies arise from some hitherto undetected temperature dependence of the torsional rigidity was investigated. The torsion elastic constant of an 1876-bp restriction fragment of pBR322 was measured by time-resolved fluorescence polarization anisotropy of intercalated ethidium over the range 278-323 K, and found to decline substantially over that interval. Simulations of a 4349-bp model DNA were performed using these measured temperature-dependent torsional rigidities. The slope, -dE(T)/dT, of the simulated data agrees satisfactorily with the slope derived from previous calorimetric measurements, but still lies substantially below that of Duguet's data. Models that involve an equilibrium between different secondary structure states with different intrinsic twists and torsion constants provide the most likely explanation for the variation of the torsion constant with T and other pertinent observations.