Transfer-matrix calculations of DNA polymer micromechanics under tension and torque constraints

Recent development of single-molecule manipulation technologies has made it possible to exert constant force and torque on individual DNA biopolymers to probe their elastic characteristics and structural stability. It has been previously shown that depending on the nature of applied mechanical constraints, DNA can exist in several forms including B-, L-, and P-DNA. However, there is still a lack of understanding of how structural heterogeneity of DNA, which may naturally arise due to sequence-dependent DNA properties, protein binding, or DNA damage, influences local stability of the above DNA states. To provide a more complete and detailed description of the DNA mechanics, we developed a theoretical framework based on transfer-matrix calculations and demonstrated how it can be used to predict the DNA behavior upon application of a wide range of force and torque constraints. The resulting phase diagram shows DNA structural transitions that are in good agreement with previous experimental and theoretical studies. We further discuss how the constructed formalism can be extended to include local inhomogeneities in the DNA physical properties, thus making it possible to investigate the effect of DNA sequence as well as protein binding on DNA structural stability.