Molecular mechanics model of supercoiled DNA

Abstract

We describe a pseudo-atomic model of supercoiled DNA. Each base-pair of the DNA is represented in the model by three particles placed in a plane. The particle triplets are stacked to model stacked base-pairs in double-helical DNA, and closed circular conformations are generated to investigate supercoiling. This model is less detailed than allatom models, which are too computationally demanding to be used to study supercoiling. On the other hand, this model contains details at the base-pair level and is therefore more elaborate than elastomechanical models. A potential energy function is written in terms of a set of internal co-ordinates defined to resemble a limited number of helical parameters. The modeled helical parameters, helical twist, base-roll, tilt and rise, are the most important parameters of the global shape of DNA. Experimentally measured mechanical properties of DNA are used to define the forces holding the particles together. We then use a procedure incorporating energy minimization and molecular dynamics to locate low energy conformations of the model DNA. The model was found to behave very much like rubber-tubing and elastomechanical models. The conformations and the effects of supercoiling pressure (a number proportional to the degree to which the total twist of the DNA has been altered from its natural value) on these conformations are all very similar to those observed in the latter two models. We also used this model to examine the effects of supercoiling pressure, base-sequence and mechanical properties on the conformations and energies of five sequences. The sequences studied include models of naturally straight DNA and DNA with static or natural bends.