Abstract Molecular docking studies of two duplex DNA sequences as target fragments and allopsoralen as ligand were performed. The calculated interaction energies showed that the ligand can be docked into the minor groove as well as become intercalated. However, unlike psoralen, allopsoralen preferred binding mode for non-poly-TA sequences is minor groove binding. Calculated energies for intercalation between different base pairs suggest that the predicted sequence selectivity for allopsoralen is analogous to that observed for psoralen. Intercalation is favored in 5′-TpA sites in poly-TA sequences.