Abstract: | The binding positions and relative minimum binding energies are calculated for complexes of 9-aminoacridine, proflavine, N-methylphenanthridinium, and ethidium in theoretically determined intercalation sites in B-DNA (sites I and II) and in unconstrained dimer-duplex sites. The selection of site I in B-DNA by these compounds agrees with the theoretical interpretation of studies of unwinding angles in closed circular DNA in all cases but ethidium, which is predicted to select site II. The most stable binding positions of the acridines and ethidium in unconstrained dimer-duplex units agree with experimental results of intercalation complexes of dinucleoside monophosphate units. Base-pair specificity for Watson-Crick pairing is examined. The energy of an intercalation complex is partitioned into ΔE23, the energy required to open base pairs BP2 and BP3 in B-DNA to a site, and ΔEIn, the energy change when a free molecular intercalates. ΔE23 depends strongly on the base-pair sequence, whereas ΔEIn for the four molecules studied does not. The three most stable sequences contain (pyrimidine)p(purine) units, and this provides a rationale for the exclusive formation of crystals of intercalation complexes with these units. In spite of this selectivity, the distribution of G?C and A?T base pairs is equal for these three units and persists as the more unstable sequences are included. Therefore, specificity arises from the interaction between the base pairs and the 2′-deoxyribose 5′-monophosphate backbone for the opening of B-DNA to an intercalation site and not from the interaction between the chromophore and the DNA. |