Mechanics and dynamics of DNA
The bending stiffness of DNA at high curvatures is fundamental to its biological activity, yet this regime has been difficult to probe experimentally and literature results have not been consistent. We created a “molecular vise” in which base-pairing interactions generated a compressive force on sub-persistence length segments of dsDNA. Short dsDNA strands (< 41 base pairs) resisted this force and remained straight; longer strands became bent, a phenomenon called “Euler buckling”. We monitored the buckling transition via FRET between appended fluorophores. For low-to-moderate concentrations of monovalent salt (up to ~150 mM), our results are in quantitative agreement with the wormlike chain (WLC) model of DNA elasticity, without the need to invoke any “kinked” states. Greater concentrations of monovalent salts or 1 mM Mg2+ induced an apparent softening of the dsDNA which was best accounted for by a kink in the region of highest curvature. We tested the effects of all single-nucleotide mismatches on the DNA bending. Remarkably, the propensity to kink correlated with the thermodynamic destabilization of the mismatched DNA relative to the perfectly complementary strand, suggesting that the kinked state is locally melted. The molecular vise is exquisitely sensitive to the sequence-dependent linear and nonlinear elastic properties of dsDNA.
.The Figure shows molecular vises used to apply compressive forces to short strands of duplex DNA. (A) Possible conformations of a molecular vise. The base-pairing force in the hairpin stem (composed of all A-T base pairs) imparted a roughly constant compressive force on the ends of the target strand (top duplex). When the target strand was shorter than the buckling length (left and center images), it withstood the compressive force and remained rigid. The dye separation increased as the target strand grew longer, resulting in a decrease of FRET efficiency. Past the buckling transition (right image), the target strand bent under the compressive force and the FRET efficiency recovered. Molecular cartoons were generated using Nucleic Acid Builder (50) and PyMOL (51). (B,C) Combined measurements of electrophoretic mobility and FRET. Molecular vises of five loop sizes (30, 36, 40, 46, and 50 nt) were each bound with varying-length target strands, analyzed by native polyacrylamide gel electrophoresis, and imaged on a commercial scanner (B). FRET efficiency was quantified from gel images and plotted for each loop size as a function of the target strand length (C). The local minimum in FRET efficiency at a target strand length of 40 bp signified the buckling transition and was consistent with the predictions of a statistical mechanical model (shaded areas; see Supporting Materials and Methods for details). The buckling transition also manifested as a change in the dependence of electrophoretic mobility on target length (B; quantified in Fig. S2). Error bars in C are SEM from four independent replicates.
- A. Fields, E. Meyer, A.E. Cohen, Euler buckling and nonlinear kinking of double-stranded DNA, Nucl. Acids Res., 41, 9881-9890, 2013.
Single-molecule DNA dynamics
Many single-molecule experiments have focused on the shape and dynamics of DNA tethered to a surface or stretched between a pair of beads, but much less attention has been paid to the dynamics of single molecules in their free-solution, relaxed, conformation. We used the ABEL trap to immobilize single molecules l-DNA and to study the internal dynamics of trapped molecules. We labeled the DNA with the intercalating dye YOYO-1. This dye lights up the molecule like a string of Christmas lights. Then we took ~60,000 frames of high speed video of single molecules fluctuating in the trap.
We analyzed these images to learn about:
- The average shape of a molecule of DNA
- The deviations about this average shape
- The dynamics of how one shape becomes another shape
We learned that the average distribution of mass about the center of mass does not follow a Gaussian distribution. A simple analytical model quantitatively agrees with the observed distribution. The key concept is that each piece of the DNA follows a Gaussian distribution about the center of mass, but different pieces follow different distributions (the ends wander further from the center of mass than does the middle). The total density distribution is the sum of many Gaussians of different widths, and thus is not Gaussian--even in the complete absence of excluded volume or other non-ideal interactions.
The figure shows a cross-section of the time-averaged density distribution, a fit to a Gaussian, and a fit to the model that includes the contributions from each element of the DNA chain. The formula for the fit is given in our PRL paper on DNA.
|By performing Principal Components Analysis (PCA) we decomposed the shape fluctuations into a set of normal modes, akin to the vibrational modes of a drumhead. We developed an analytical model for these principal components that quantitatively agrees with our data. The eigenvalue associated with each principal component is a measure of the entropic stiffness of the DNA in that mode. The eigenvalues and the comparison to theory are discussed in our PNAS paper on DNA.|
- Adam Cohen, W. E. Moerner: Principal Components Analysis of shape fluctuations of single DNA molecules, PNAS, 104, 12622-12627, 31 July 2007. Journal link and Supporting Material and Movies.
- Adam Cohen, W. E. Moerner: Internal mechanical response of a polymer in solution, Phys. Rev. Lett., 98, 116001 (2007). Supplementary material.