The solvation contribution to the binding energy of DNA with non-intercalating antibiotics

The influence of the solvent on the binding energies to DNA of six non-intercalating antibiotics - netropsin, distamycin-3, distamycin-2, SN 18071, berenil and stilbamidine - is evaluated by combining the effect of the first hydration shell with that of bulk water. The first effect is computed by a methodology based on a spherical/point dipole model of water and limited to electrostatic interaction energies. Hydration shells are obtained which are energy optimized with respect to both water-solute and water-water interactions for the complexes and for the isolated DNA oligomers and ligands. The method allows even very large complexes to be studied in reasonable computation times. The second effect is introduced via a cavity treatment. It is shown that if the vacuum interaction energies already predict correctly the preference of the ligands for the minor groove of AT sequences of B-DNA, the introduction of the solvation effect is indispensable for reproducing the order of affinity of the ligands and for bringing the values of the complexation energies into close agreement with experimental data.