Computer simulation of DNA supercoiling

We treat supercoiled DNA within a wormlike model with excluded volume. A modified Monte Carlo approach has been used, which allowed computer statistical-mechanical simulations of moderately and highly supercoiled DNA molecules. Even highly supercoiled molecules do not have a regular shape, though with an increase in writhing the chains look more and more like branched interwound helixes. The averaged writhing (Wr) approximately 0.7 delta Lk. The superhelical free energy F is calculated as a function of the linking number. Lk. The calculations have shown that the generally accepted quadratic dependence of F on Lk is valid for a variety of conditions, though it is by no means universal. Significant deviations from the quadratic dependence are expected at high superhelical density under ionic conditions where the effective diameter of DNA is small. The results are compared with the available experimental data.