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).poly(dG-dC). Altogether, however, among the six investigated sequences, the strongest complexes are computed for the mixed hexanucleotides d(CGATCG)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.(ABSTRACT TRUNCATED AT 250 WORDS)