Abstract

Psoralen interaction with two models of DNA was investigated using molecular mechanics and molecular dynamics methods. Calculated energies of minor groove binding and intercalation were compared in order to define a preferred binding mode for the ligand. We found that both binding modes are possible, explaining the low efficiency for monoadduct formation from intercalated ligands. A comparison between the interaction energy for intercalation between different base pairs suggests that the observed sequence selectivity is due to favorable intercalation in 5-TpA in (AT)n sequences.