A Theoretical Investigation on the Sequence Selective Binding of Daunomycin to Double-Stranded Polynucleotides

Abstract
Theoretical computations are performed on the structural and energetical factors involved in the sequence selective binding of daunomycin (DNM) to six representative self- complementary double-stranded hexanucleotides: d(CGTACG)2, d(CGATCG)2, d(CITACI)2, d(TATATA)2, d(CGCGCG)2 and d(TACGTA)2. The conformational angles of the hexanucleotides are fixed in values found in the representative crystal structure of the d(CGTACG)2- DNM complex. The intermolecular DNM-hexanucleotide interaction energies and the conformational energy changes of DNM upon binding are computed and optimized in the framework of the SIBFA procedure, which uses empirical formulas based on ab initio SCF computations. Among the two regularly alternating hexanucleotides, d(TATATA)2 and d(CGCGCG)2, a stronger binding is predicted for the former, in agreement with experimental results obtained with poly(dA-dT)-poly(dA-dT) and poly(dG-dC)<2 and d(CGTACG)2, containing the intercalation site between two CG base pairs and an adjacent TA base pair. This situation may be related to the increased affinity of DNM for GC rich DNA's and to the situation in the crystal structure of the DNM-d(CGTACG)2 complex. Analysis of the intrinsic base sequence preferences expressed by the individual constituents of DNM, namely the daunosamine side chain, the chromophore ring and its two 9-hydroxyl and 9-acetoxy substituents, reveals that the overall sequence preference found is the result of a rather intricate interplay of intrinsic sequence preferences, in particular at the level of daunosamine and the 9-hydroxyl substituent. Altogether, it is seen that the selective base pair recognition by daunomycin cannot, in general, be defined in terms of the two base pairs implicated in the intercalation site alone (with the exception of homogeneous AT or GC base sequences) but must be expressed in terms of a triplet of base pairs.