Search for rigid sub domains in DNA from molecular dynamics simulations

A strategy is presented for searching which atoms can be regrouped within rigid sub-units during the time course of Molecular Dynamics simulations of biopolymers. The root mean square fluctuations of the interatomic distances are used as a criterion. The number of rigid sub-units which are found depends on the tolerance r(c) for the definition of a rigid body, i.e. until which value the fluctuations can be neglected. The method is applied to two self-complementary oligonucleotides belonging to the B-form family which give identical results. With r(c)=0.027 nm each nucleotide may be described as 3 rigid sub-units: the sugar ring, the base and the backbone (PO4 + C5' atoms). With r(c)=0.01 nm, the same sub-units are obtained except that C5' can no more be regrouped with the PO4 atoms. It is shown that the variation of the coulombic potential owing to the deformation of the sub-units during the time course of the simulation is on the same order of magnitude as the inaccuracy due to the choice of the force field parameters.