Designing nucleosomal force sensors.

About three quarters of our DNA is wrapped into nucleosomes: DNA spools with a protein core. It is well known that the affinity of a given DNA stretch to be incorporated into a nucleosome depends on the geometry and elasticity of the basepair sequence involved, causing the positioning of nucleosomes. Here we show that DNA elasticity can have a much deeper effect on nucleosomes than just their positioning: it affects their "identities". Employing a recently developed computational algorithm, the mutation Monte Carlo method, we design nucleosomes with surprising physical characteristics. Unlike any other nucleosomes studied so far, these nucleosomes are short-lived when put under mechanical tension whereas other physical properties are largely unaffected. This suggests that the nucleosome, the most abundant DNA-protein complex in our cells, might more properly be considered a class of complexes with a wide array of physical properties, and raises the possibility that evolution has shaped various nucleosome species according to their genomic context.

Stiffer double-stranded DNA in two-dimensional confinement due to bending anisotropy.

Using analytical approach and Monte Carlo (MC) simulations, we study the elastic behavior of the intrinsically twisted elastic ribbons with bending anisotropy, such as double-stranded DNA (dsDNA), in two-dimensional (2D) confinement. We show that, due to the bending anisotropy, the persistence length of dsDNA in 2D conformations is always greater than three-dimensional (3D) conformations. This result is in consistence with the measured values for DNA persistence length in 2D and 3D in equal biological conditions. We also show that in two dimensions, an anisotropic, intrinsically twisted polymer exhibits an implicit twist-bend coupling, which leads to the transient curvature increasing with a half helical turn periodicity along the bent polymer.

Extreme bendability of DNA double helix due to bending asymmetry.

Experimental data of the DNA cyclization (J-factor) at short length scales exceed the theoretical expectation based on the wormlike chain (WLC) model by several orders of magnitude. Here, we propose that asymmetric bending rigidity of the double helix in the groove direction can be responsible for extreme bendability of DNA at short length scales and it also facilitates DNA loop formation at these lengths. To account for the bending asymmetry, we consider the asymmetric elastic rod (AER) model which has been introduced and parametrized in an earlier study [B. Eslami-Mossallam and M. R. Ejtehadi, Phys. Rev. E 80, 011919 (2009)]. Exploiting a coarse grained representation of the DNA molecule at base pair (bp) level and using the Monte Carlo simulation method in combination with the umbrella sampling technique, we calculate the loop formation probability of DNA in the AER model. We show that the DNA molecule has a larger J-factor compared to the WLC model which is in excellent agreement with recent experimental data.

Confinement dynamics of a semiflexible chain inside nano-spheres.

We study the conformations of a semiflexible chain, confined in nano-scaled spherical cavities, under two distinct processes of confinement. Radial contraction and packaging are employed as two confining procedures. The former method is performed by gradually decreasing the diameter of a spherical shell which envelopes a confined chain. The latter procedure is carried out by injecting the chain inside a spherical shell through a hole on the shell surface. The chain is modeled with a rigid body molecular dynamics simulation and its parameters are adjusted to DNA base-pair elasticity. Directional order parameter is employed to analyze and compare the confined chain and the conformations of the chain for two different sizes of the spheres are studied in both procedures. It is shown that for the confined chains in the sphere sizes of our study, they appear in spiral or tennis-ball structures, and the tennis-ball structure is more likely to be observed in more compact confinements. Our results also show that the dynamical procedure of confinement and the rate of the confinement are influential parameters of the structure of the chain inside spherical cavities.

Definition of the persistence length in the coarse-grained models of DNA elasticity.

By considering the detailed structure of DNA in the base pair level, two possible definitions of the persistence length are compared. One definition is related to the orientation of the terminal base pairs, and the other is based on the vectors which connect two adjacent base pairs at each end of the molecule. It is shown that although these definitions approach each other for long DNA molecules, they are dramatically different on short length scales. We show analytically that the difference mostly comes from the shear flexibility of the molecule and can be used to measure the shear modulus of DNA.

Contribution of nonlocal interactions to DNA elasticity.

A nonlocal harmonic elastic rod model is proposed to describe the elastic behavior of short DNA molecules. We show that the nonlocal interactions contribute to effective bending energy of the molecule and affect its apparent persistence length. It is also shown that the anomalous behavior which has been observed in all-atom molecular dynamic simulations [A. K. Mazur, Biophys. J. 134, 4507 (2006)] can be a consequence of both nonlocal interactions between DNA base pairs and the intrinsic curvature of DNA.

Asymmetric elastic rod model for DNA.

In this paper we consider the anharmonic corrections to the anisotropic elastic rod model for DNA. Our model accounts for the difference between the bending energies of positive and negative rolls, which comes from the asymmetric structure of the DNA molecule. We will show that the model can explain the high flexibility of DNA at small length scales, as well as kink formation at high deformation limit.

Stretching an anisotropic DNA.

We present a perturbation theory to find the response of an anisotropic DNA to the external tension. It is shown that the anisotropy has a nonzero but small contribution to the force-extension curve of the DNA. Thus an anisotropic DNA behaves like an isotropic one with an effective bending constant equal to the harmonic average of its soft and hard bending constants.