A recent study described an allosteric effect in which the binding of a protein to DNA is influenced by another protein bound nearby. The effect shows a periodicity of ∼10 basepairs and decays with increasing protein-protein distance. As a mechanistic explanation, the authors reported a similar periodic, decaying pattern of the correlation coefficient between major groove widths inferred from a shorter molecular dynamics simulation. Here we show that in a state-of-the-art, microsecond-long simulation of the same DNA sequence, the periodicity of the correlation coefficient is not observed. To study the problem further, we extend an earlier mechanical model of DNA allostery based on constrained minimization of effective quadratic deformation energy of the DNA. We demonstrate that, if the constraints mimicking the bound proteins are properly applied, the periodicity in the binding energy is indeed recovered.

• MeSH
• Allosteric Regulation MeSH
• DNA chemistry   MeSH
• Nucleic Acid Conformation * MeSH
• Molecular Dynamics Simulation * MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA MeSH

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
• Biomechanical Phenomena MeSH
• DNA chemistry   MeSH
• RNA, Double-Stranded chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Base Pairing MeSH
• Molecular Dynamics Simulation MeSH
• Torque MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Comparative Study MeSH
• Names of Substances
• DNA MeSH
• RNA, Double-Stranded MeSH