Double stranded helical DNA and RNA are flexible molecules that can undergo global conformational fluctuations. Their bending, twisting and stretching deformabilities are of similar magnitude. However, recent single-molecule experiments revealed a striking qualitative difference indicating an opposite sign for the twist-stretch couplings of dsDNA and dsRNA [Lipfert et al. 2014. Proc. Natl. Acad. Sci. U.S.A. 111, 15408] that is not explained by existing models. Employing unconstrained Molecular Dynamics (MD) simulations we are able to reproduce the qualitatively different twist-stretch coupling for dsDNA and dsRNA in semi-quantitative agreement with experiment. Similar results are also found in simulations that include an external torque to induce over- or unwinding of DNA and RNA. Detailed analysis of the helical deformations coupled to twist indicate that the interplay of helical rise, base pair inclination and displacement from the helix axis upon twist changes are responsible for the different twist-stretch correlations. Overwinding of RNA results in more compact conformations with a narrower major groove and consequently reduced helical extension. Overwinding of DNA decreases the size of the minor groove and the resulting positive base pair inclination leads to a slender and more extended helical structure.

• MeSH
• biomechanika MeSH
• DNA chemie   MeSH
• dvouvláknová RNA chemie   MeSH
• konformace nukleové kyseliny MeSH
• párování bází MeSH
• simulace molekulární dynamiky MeSH
• točivý moment MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• srovnávací studie MeSH
• Názvy látek
• DNA MeSH
• dvouvláknová RNA MeSH

A-tracts are functionally important DNA sequences which induce helix bending and have peculiar structural properties. While A-tract structure has been qualitatively well characterized, their mechanical properties remain controversial. A-tracts appear structurally rigid and resist nucleosome formation, but seem flexible in DNA looping. In this work, we investigate mechanical properties of symmetric AnTn and asymmetric A2n tracts for n = 3, 4, 5 using two types of coarse-grained models. The first model represents DNA as an ensemble of interacting rigid bases with non-local quadratic deformation energy, the second one treats DNA as an anisotropically bendable and twistable elastic rod. Parameters for both models are inferred from microsecond long, atomic-resolution molecular dynamics simulations. We find that asymmetric A-tracts are more rigid than the control G/C-rich sequence in localized distortions relevant for nucleosome formation, but are more flexible in global bending and twisting relevant for looping. The symmetric tracts, in contrast, are more rigid than asymmetric tracts and the control, both locally and globally. Our results can reconcile the contradictory stiffness data on A-tracts and suggest symmetric A-tracts to be more efficient in nucleosome exclusion than the asymmetric ones. This would open a new possibility of gene expression manipulation using A-tracts.

• MeSH
• adenin chemie   MeSH
• biomechanika MeSH
• DNA chemie   MeSH
• entropie MeSH
• konformace nukleové kyseliny MeSH
• nukleozomy chemie   MeSH
• poly A chemie   MeSH
• sekvence nukleotidů MeSH
• simulace molekulární dynamiky MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• adenin MeSH
• DNA MeSH
• nukleozomy MeSH
• poly A MeSH