An effective mesoscopic model of double-stranded DNA

Base pair fluctuations in helical models for nucleic acids

A statistical method is developed to estimate the maximum amplitude of the base pair fluctuations in a three dimensional mesoscopic model for nucleic acids. The base pair thermal vibrations around the helix diameter are viewed as a Brownian motion for a particle embedded in a stable helical structure. The probability to return to the initial position is computed, as a function of time, by integrating over the particle paths consistent with the physical properties of the model potential. The zero time condition for the first-passage probability defines the constraint to select the integral cutoff for various macroscopic helical conformations, obtained by tuning the twist, bending, and slide motion between adjacent base pairs along the molecule stack. Applying the method to a short homogeneous chain at room temperature, we obtain meaningful estimates for the maximum fluctuations in the twist conformation with ∼10.5 base pairs per helix turn, typical of double stranded DNA helices. Untwisting the double helix, the base pair fluctuations broaden and the integral cutoff increases. The cutoff is found to increase also in the presence of a sliding motion, which shortens the helix contour length, a situation peculiar of dsRNA molecules.

Ionic effects on the temperature-force phase diagram of DNA

Double-stranded DNA (dsDNA) undergoes a structural transition to single-stranded DNA (ssDNA) in many biologically important processes such as replication and transcription. This strand separation arises in response either to thermal fluctuations or to external forces. The roles of ions are twofold, shortening the range of the interstrand potential and renormalizing the DNA elastic modulus. The dsDNA-to-ssDNA transition is studied on the basis that dsDNA is regarded as a bound state while ssDNA is regarded as an unbound state. The ground state energy of DNA is obtained by mapping the statistical mechanics problem to the imaginary time quantum mechanics problem. In the temperature-force phase diagram the critical force F _c (T) increases logarithmically with the Na⁺ concentration in the range from 32 to 110 mM. Discussing this logarithmic dependence of F _c (T) within the framework of polyelectrolyte theory, it inevitably suggests a constraint on the difference between the interstrand separation and the length per unit charge during the dsDNA-to-ssDNA transition.

The "sugar" coarse-grained DNA model

More than 20 coarse-grained (CG) DNA models have been developed for simulating the behavior of this molecule under various conditions, including those required for nanotechnology. However, none of these models reproduces the DNA polymorphism associated with conformational changes in the ribose rings of the DNA backbone. These changes make an essential contribution to the DNA local deformability and provide the possibility of the transition of the DNA double helix from the B-form to the A-form during interactions with biological molecules. We propose a CG representation of the ribose conformational flexibility. We substantiate the choice of the CG sites (six per nucleotide) needed for the "sugar" GC DNA model, and obtain the potentials of the CG interactions between the sites by the "bottom-up" approach using the all-atom AMBER force field. We show that the representation of the ribose flexibility requires one non-harmonic and one three-particle potential, the forms of both the potentials being different from the ones generally used. The model also includes (i) explicit representation of ions (in an implicit solvent) and (ii) sequence dependence. With these features, the sugar CG DNA model reproduces (with the same parameters) both the B- and A- stable forms under corresponding conditions and demonstrates both the A to B and the B to A phase transitions. Graphical Abstract The proposed coarse-grained DNA model allows to reproduce both the B- and A- DNA forms and the transitions between them under corresponding conditions.

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.

Physical origin of DNA unzipping

In DNA transcription, the base pairs are unzipped in response to the enzymatic forces, separating apart two intertwined nucleotide strands. Consequently, the double-stranded DNA (dsDNA), in which two nucleotide strands wind about each other, transits structurally to the single-stranded DNA (ssDNA) in which two nucleotide strands are completely unwound and separated. The large interstrand separation is intimately related to the softening nucleotide strands. This conceptual framework is reinforced with the flow of the bending modulus toward zero under recursion relations derived from the momentum shell renormalization group. Interestingly, the stretch modulus remains the same under recursion relations. The renormalization of the bending modulus to zero has a profound implication that ssDNA has the shorter bending persistence length than does dsDNA in accordance with experiments.