Mechanical properties of nucleic acids and the non-local twistable wormlike chain model

Abstract

Mechanical properties of nucleic acids play an important role in many biological processes that often involve physical deformations of these molecules. At sufficiently long length scales (say, above [Formula: see text] base pairs), the mechanics of DNA and RNA double helices is described by a homogeneous Twistable Wormlike Chain (TWLC), a semiflexible polymer model characterized by twist and bending stiffnesses. At shorter scales, this model breaks down for two reasons: the elastic properties become sequence-dependent and the mechanical deformations at distal sites get coupled. We discuss in this paper the origin of the latter effect using the framework of a non-local Twistable Wormlike Chain (nlTWLC). We show, by comparing all-atom simulations data for DNA and RNA double helices, that the non-local couplings are of very similar nature in these two molecules: couplings between distal sites are strong for tilt and twist degrees of freedom and weak for roll. We introduce and analyze a simple double-stranded polymer model that clarifies the origin of this universal distal couplings behavior. In this model, referred to as the ladder model, a nlTWLC description emerges from the coarsening of local (atomic) degrees of freedom into angular variables that describe the twist and bending of the molecule. Different from its local counterpart, the nlTWLC is characterized by a length-scale-dependent elasticity. Our analysis predicts that nucleic acids are mechanically softer at the scale of a few base pairs and are asymptotically stiffer at longer length scales, a behavior that matches experimental data.