Modeling DNA loops using the theory of elasticity

A versatile approach to modeling the conformations and energetics of DNA loops is presented. The model is based on the classical theory of elasticity, modified to describe the intrinsic twist and curvature of DNA, the DNA bending anisotropy, and electrostatic properties. All the model parameters are considered to be functions of the loop arclength, so that the DNA sequence-specific properties can be modeled. The model is applied to the test case study of a DNA loop clamped by the lac repressor protein. Several topologically different conformations are predicted for various lengths of the loop. The dependence of the predicted conformations on the parameters of the problem is systematically investigated. Extensions of the presented model and the scope of the model's applicability, including multi-scale simulations of protein-DNA complexes and building all-atom structures on the basis of the model, are discussed