Computation of writhe in modeling of supercoiled DNA

Equilibrium fluctuations of DNA plectonemes

Plectonemes are intertwined helically looped domains which form when a DNA molecule is supercoiled, i.e., over- or underwound. They are ubiquitous in cellular DNA, and their physical properties have attracted significant interest both from the experimental side and from the modeling side. In this paper, we investigate fluctuations of the end-point distance z of supercoiled linear DNA molecules subject to external stretching forces. Our analysis is based on a two-phase model, which describes the supercoiled DNA as composed of a stretched phase and a plectonemic phase. A variety of mechanisms are found to contribute to extension fluctuations, characterized by the variance 〈Δz^{2}〉. We find the dominant contribution to 〈Δz^{2}〉 to originate from phase-exchange fluctuations, the transient shrinking and expansion of plectonemes, which is accompanied by an exchange of molecular length between the two phases. We perform Monte Carlo simulations of the twistable wormlike chain and analyze the fluctuation of various quantities, the results of which are found to agree with the two-phase model predictions. Furthermore, we show that the extension and its variance at high forces are very well captured by the two-phase model, provided that one goes beyond quadratic approximations.

Computer simulations of melts of ring polymers with nonconserved topology: A dynamic Monte Carlo lattice model

We present computer simulations of a dynamic Monte Carlo algorithm for polymer chains on a fcc lattice which explicitly takes into account the possibility to overcome topological constraints by controlling the rate at which nearby polymer strands may cross through each other. By applying the method to systems of interacting ring polymers at melt conditions, we characterize their structure and dynamics by measuring, in particular, the amounts of knots and links which are formed during the relaxation process. In comparison with standard melts of unknotted and unconcatenated rings, our simulations demonstrate that the mechanism of strand crossing makes polymer dynamics faster provided the characteristic timescale of the process is smaller than the typical timescale for chain relaxation in the unperturbed state, in agreement with recent experiments employing solutions of DNA rings in the presence of the type II topoisomerase enzyme. In the opposite case of slow rates the melt is shown to become slower, and this prediction may be easily validated experimentally.

Characterization of entanglements in glassy polymeric ensembles using the Gaussian linking number

We propose a method for enumerating entanglements between long chained, linear polymers that is based on the Gaussian linking number. The linking number is calculated between closely approaching segments of the macromolecular chains. Topological features of an entanglement, i.e., the extent to which one open segment winds around another, are reflected by the linking number. We show that using this measure, we can track disentanglement events through a deformation history and gain insights into how large scale disentanglements lead to failure. Incorporating an additional step where the topological entanglements identified along each chain are optimally clustered using standard clustering algorithms, we can also obtain a measure of the average number of rheological constraints that exist along each chain in an ensemble. Comparisons with other methods of enumerating entanglements, especially the primitive path analysis, are also made. Our results indicate that the linking number between two entangled segments in the undeformed state is a good indicator of the strength of the entanglement. Also, disentanglements occurring overwhelmingly around chain ends are an important cause of failure when a triaxial stress state exists in the polymer.

Energy required to pinch a DNA plectoneme

DNA supercoiling plays an important role from a biological point of view. One of its consequences at the supramolecular level is the formation of DNA superhelices named plectonemes. Normally separated by a distance on the order of 10 nm, the two opposite double strands of a DNA plectoneme must be brought closer if a protein or protein complex implicated in genetic regulation is to be bound simultaneously to both strands, as if the plectoneme was locally pinched. We propose an analytic calculation of the energetic barrier, of elastic nature, required to bring closer the two loci situated on the opposed double strands. We examine how this energy barrier scales with the DNA supercoiling. For physically relevant values of elastic parameters and of supercoiling density, we show that the energy barrier is in the k_{B}T range under physiological conditions, thus demonstrating that the limiting step to loci encounter is more likely the preceding plectoneme slithering bringing the two loci side by side.

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.