  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.
