An algorithm for studying cooperative transitions in DNA

Cooperative transitions in DNA (B to Z, B to A, helix to coil, etc.) are known to depend strongly on nucleotide sequence. In general the change in free energy involved in the transition can be expressed as: delta G(seq) = 2RT log (sigma) where sigma is a factor arising from the free energy associated with boundaries of different conformations along the molecule. This formula allows to infer a general algorithm with which DNA sequences can be partitioned into well defined domains in which, under suitable conditions, base pairs change state cooperatively. The different partitions of the sequence that can be generated by varying the values of the physical parameters involved in the above formula, are shown to be embedded into a binary tree hierarchy. Application to a reliable prediction of Z-DNA antibody binding sites will be illustrated for the 0X174 genome. Possible biological implications are briefly discussed.