Sequence directed flexibility of DNA and the role of cross-strand hydrogen bonds

Persistence length and torsional rigidity for different B-DNA sequences have been calculated by analysing crystal structure database. The values of these parameters for mixed sequence DNA are in good agreement with those estimated by others. Persistence lengths for the homopolymeric sequences, namely poly(da).poly(dT) and poly(dG).poly(dC), are significantly large compared to those of others as expected from the inability of these sequences to form nucleosome under normal conditions. The heteropolymeric sequences poly(dA-dC).poly(dc-dT) and poly(dG-dC).poly(dC-dC), on the other hand, have smaller persistence lengths. This implies larger flexibility of the d(AC).d(GT), d(CA).d(TG), d(GC).d(GC) and d(CG).d(CG) doublers, some of which constitute: the genetic disease forming triplet repeats d(CTG).d(CAG) and d(CGG).d(CCG). Thus it is expected that these triplet repeat sequences are also flexible and wrap around the histone octamer efficiently. Persistence length calculations also indicate larger flexibility for these tripler repeat sequences. Furthermore, our computations reveal that the rigidity of a given DNA sequence is controlled by its ability to form cross-strand bifurcated hydrogen bonds between the successive base pairs. Molecular orbital calculations suggest that these hydrogen bonds are generally extended with bond lengths around 3 Angstrom.