Melting profile and temperature dependent binding constant of an anticancer drug daunomycin-DNA complex

We calculate thermal fluctuational base pair opening probability and the drug binding constant of a daunomycin-bound Poly d(CGTA).Poly d(TACG) at temperatures from room temperature to its melting temperature. For comparison we also carry out a calculation on a drug-free DNA with the same sequence. Our calculations are carried out by means of a statistical approach using microscopic structures and established force fields and with cooperative effects incorporated into the algorithm. Both hydrogen bond disruption probabilities and drug unstacking probability are determined self-consistently. These probabilities are then used to determine temperature dependent base pair opening probabilities and the drug binding constant. The calculated base pair opening probabilities and drug binding constant are found to be in fair agreement with experiments carried out at room temperature. Our calculation shows cooperative base pair disruption and drug dissociation at certain critical temperatures close to the observed melting temperatures for similar helices. We find that the temperature dependence of the drug binding constant fits well to the van't Hoff relation, in agreement with observations. Our calculation indicates the occurrence of a premelting transition in the drug-bound DNA helix. Some comments are made about this premelting transition.