The flexibility of short DNA chains is investigated via computation of the average correlation function between dimers which defines the persistence length. Path integration techniques have been applied to confine the phase space available to base pair fluctuations and derive the partition function. The apparent persistence lengths of a set of short chains have been computed as a function of the twist conformation both in the over-twisted and the untwisted regimes, whereby the equilibrium twist is selected by free energy minimization. The obtained values are significantly lower than those generally attributed to kilo-base long DNA. This points to an intrinsic helix flexibility at short length scales, arising from large fluctuational effects and local bending, in line with recent experimental indications. The interplay between helical untwisting and persistence length has been discussed for a heterogeneous fragment by weighing the effects of the sequence specificities through the non-linear stacking potential.
1The DNA double helical structure is stable enough to preserve genetic information encoded in the Watson-Crick paired bases and also loose enough to allow for those transient base pair (bp) openings [1] which make the code accessible to enzymes during the processes of replication, transcription and repair. At physiological temperatures DNA molecules fluctuate between a variety of random coil conformations in which even distant segments along the helical axis can be brought in close proximity [2]. This points to an inherent flexibility of the DNA chain which has been widely probed over the last twenty five years [3]. While these experiments demonstrate that stretch and twist elasticity are intertwined [4], they also call for models in which bp fluctuations and stacking interactions are considered as dependent on the specific twist conformation of the molecule. Modeling of the helix and its conformational states can be carried out at different levels of resolution ranging from all-atom simulations to continuous worm-like chain (WLC) models which simply treat DNA as a homogeneous and inextensible rod [5], not accounting for the interplay between twist and bp fluctuations. This may explain the shortcomings of the WLC model emerged in the analysis of the cyclization properties [6,7] at those short length scales in which the details of the bp interactions matter. In this regard, mechanical models such as the Dauxois-Peyrard-Bishop (DPB) model