Effect of undulations on spontaneous braid formation

Nucleation of Multiple Buckled Structures in Intertwined DNA Double Helices

We study the statistical-mechanical properties of intertwined double-helical DNAs (DNA braids). In magnetic tweezers experiments, we find that torsionally stressed stretched braids supercoil via an abrupt buckling transition, which is associated with the nucleation of a braid end loop, and that the buckled braid is characterized by a proliferation of multiple domains. Differences between the mechanics of DNA braids and supercoiled single DNAs can be understood as an effect of the increased bulkiness in the structure of the former. The experimental results are in accord with the predictions of a statistical-mechanical model.

Torque and buckling in stretched intertwined double-helix DNAs

We present a statistical-mechanical model for the behavior of intertwined DNAs, with a focus on their torque and extension as a function of their catenation (linking) number and applied force, as studied in magnetic tweezers experiments. Our model produces results in good agreement with available experimental data and predicts a catenation-dependent effective twist modulus distinct from what is observed for twisted individual double-helix DNAs. We find that buckling occurs near the point where experiments have observed a kink in the extension versus linking number, and that the subsequent "supercoiled braid" state corresponds to a proliferation of multiple small plectoneme structures. We predict a discontinuity in extension at the buckling transition corresponding to nucleation of the first plectoneme domain. We also find that buckling occurs for lower linking number at lower salt; the opposite trend is observed for supercoiled single DNAs.

Statistical mechanical model for a closed loop plectoneme with weak helix specific forces

We develop a statistical mechanical framework, based on a variational approximation, to describe closed loop plectonemes. This framework incorporates weak helix structure dependent forces into the determination of the free energy and average structure of a plectoneme. Notably, due to their chiral nature, helix structure dependent forces break the symmetry between left and right handed supercoiling. The theoretical approach, presented here, also provides a systematic way of enforcing the topological constraint of closed loop supercoiling in the variational approximation. At large plectoneme lengths, by considering correlation functions in an expansion in terms of the spatial mean twist density about its thermally averaged value, it can be argued that topological constraint may be approximated by replacing twist and writhe by their thermal averages. A Lagrange multiplier, containing the sum of average twist and writhe, can be added to the free energy to conveniently inforce this result. The average writhe can be calculated through the thermal average of the Gauss' integral in the variational approximation. Furthermore, this approach allows for a possible way to calculate finite size corrections due to the topological constraint. Using interaction energy terms from the mean-field Kornyshev-Leikin theory, for parameter values that correspond to weak helix dependent forces, we calculate the free energy, fluctuation magnitudes and mean geometric parameters for the plectoneme. We see a slight asymmetry, where interestingly, left handed supercoils have a looser structure than right handed ones, although with a lower free energy, unlike what the previous ground state calculations would suggest.

Collapse and coexistence for a molecular braid with an attractive interaction component subject to mechanical forces

Dual mechanical braiding experiments provide a useful tool with which to investigate the nature of interactions between rod-like molecules, for instance actin and DNA. In conditions close to molecular condensation, one would expect an appearance of a local minimum in the interaction potential between the two molecules. We investigate this situation, introducing an attractive component into the interaction potential, using a model developed for describing such experiments. We consider both attractive interactions that do not depend on molecular structure and those which depend on a DNA-like helix structure. In braiding experiments, an attractive term may lead to certain effects. A local minimum may cause molecules to collapse from a loosely braided configuration into a tight one, occurring at a critical value of the moment applied about the axis of the braid. For a fixed number of braid pitches, this may lead to coexistence between the two braiding states, tight and loose. Coexistence implies certain proportions of the braid are in each state, their relative size depending on the number of braid pitches. This manifests itself as a linear dependence in numerically calculated quantities as functions of the number of braid pitches. Also, in the collapsed state, the braid radius stays roughly constant. Furthermore, if the attractive interaction is helix dependent, the left-right handed braid symmetry is broken. For a DNA like charge distribution, using the Kornyshev-Leikin interaction model, our results suggest that significant braid collapse and coexistence only occurs for left handed braids. Regardless of the interaction model, the study highlights the possible qualitative physics of braid collapse and coexistence; and the role helix specific forces might play, if important. The model could be used to connect other microscopic theories of interaction with braiding experiments.

Undulations in a weakly interacting mechanically generated molecular braid under tension

We consider mechanically generated molecular braids composed of two molecules where long range interactions between them can be considered to be very weak. We describe a model that takes account of the thermal fluctuations of the braid, steric interactions between the molecules, and external mechanical forces. In this model, both sets of ends, of the two molecules, are considered to be separated by a fixed distance much larger than the radius of the braid. One set of ends is rotated to generate a braid of a certain number of pitches (or turns), while the other set remains fixed. This model may describe the situation in which the ends of each molecule are attached to a substrate and a magnetic bead; to the latter a pulling force and rotational torque can be applied. We discuss various aspects of our model. Most importantly, an expression for the free energy is given, from which equations, determining the various geometric parameters of the braid, can be obtained. By numerically solving these equations, we give predictions from the model for the external torque needed to produce a braid with a certain number of turns per bending persistence length, as well as the end to end extension of the two molecules for a given pulling force. Other geometric parameters, as well as the lateral force required to keep the ends of the two molecules apart, are also calculated.

Self-consistent treatment of electrostatics in molecular DNA braiding through external forces

In this paper we consider a physical system in which two DNA molecules braid about each other. The distance between the two molecular ends, on either side of the braid, is held at a distance much larger than supercoiling radius of the braid. The system is subjected to an external pulling force, and a moment that induces the braiding. In a model, developed for understanding such a system, we assume that each molecule can be divided into a braided and unbraided section. We also suppose that the DNA is nicked so that there is no constraint of the individual linking numbers of the molecules. Included in the model are steric and electrostatic interactions, thermal fluctuations of the braided and unbraided sections of the molecule, as well as the constraint on the braid linking (catenation) number. We compare two approximations used in estimating the free energy of the braided section. One is where the amplitude of undulations of one molecule with respect to the other is determined only by steric interactions. The other is a self-consistent determination of the mean-squared amplitude of these undulations. In this second approximation electrostatics should play an important role in determining this quantity, as suggested by physical arguments. We see that if the electrostatic interaction is sufficiently large there are indeed notable differences between the two approximations. We go on to test the self-consistent approximation-included in the full model-against experimental data for such a system, and we find good agreement. However, there seems to be a slight left-right-handed braid asymmetry in some of the experimental results. We discuss what might be the origin of this small asymmetry.