Sequence and temperature effect on hydrogen bond disruption in DNA determined by a statistical analysis

We use the modified self-consistent phonon approximation theory to calculate temperature dependent interbase hydrogen bond disruption profiles for a number of six base pair repeating sequence infinite B-DNA polymers with various guanine-cytosine/adenine-thymine ratios, For comparison we also include results we have obtained in our earlier work on several B-DNA homopolymers, copolymers and a four-base-pair repeating sequence polymer. Our theory gives a statistical estimate of thermal fluctuational disruption probability of individual hydrogen bonds in individual base pairs in DNA as a function of temperature. The calculated probabilities show no sequence dependence at premelting temperatures, in agreement with proton exchange measurements. These probabilities however become very sensitive to base sequence at temperatures close to the observed melting temperatures. Multi phasic critical transitions are found in which a portion of base pairs are disrupted at temperatures below the final disruption temperature. These transitions include localized as well as non-localized base pair opening. The localized transitions involve disruption of a few base-pairs at every other location without large scale base unstacking, and they may not appear in the observed UV curves with current resolution. On the other hand the overall disruption behavior is consistent with observations. The midpoint transition temperatures are close to the observed melting temperatures and these temperatures show the observed linear dependence on guanine-cytosine content. Our calculations indicate that our theory can be used effectively to calculate H-bond disruption behavior of different DNA sequences.