Molecular modelling and footprinting studies of DNA minor groove binders: bisquaternary ammonium heterocyclic compounds

We report new quantitative footprinting data which reveal differences in binding constants of bisquaternary ammonium heterocyclic compounds (BQA) with AT-rich DNA sites depending on the ligand structure and on the size and sequence of the DNA binding site. In an attempt to understand the dependence of binding affinity on the ligand structure we have performed quantum-chemical AM1 calculations on the BQA compounds and on subunits to explore the conformational space and to calculate the electronic and structural features of individual ligand conformations. Due to the properties of the rotatable backbone bonds, there is a large number of possible conformations with almost equal energy. We present a new method for the calculation of the radius of curvature of molecular structures. Assuming that strong binders should have a shape complementary to the DNA minor groove, this measure is used to select the optimum conformations for DNA-drug binding. The approach yields the correct ligand conformation for SN6999, for which an X-ray DNA-drug structure is known. The curvature of the optimum conformations of all ligands is compared with the experimental binding constants. A correlation is found between curvature and binding constant provided other structural factors do not vary. Therefore, we conclude that within structurally similar BQA compounds the extent of curvature is the relevant quantity which modulates the binding affinity.