The complex formation of S- and R-ofloxacin with the self-complementary oligonucleotides, namely d[ATAGCGCTAT]₂, d[GCGATATCGC]₂ and d[ATAICICTAT]₂, were investigated by the molecular dynamics (MD) simulation. Four starting positions, including two intercalation positions with different insertion directions and two minor groove binding positions, were considered. The total energy of both S- and R-ofloxacin-d[ATAGCGCTAT]₂ complex, in which ofloxacin binds in the minor groove of the oligonucleotide, were lower than any intercalation binding mode. For both enantiomers, formation of the complex with GC oligonucleotide is more favorable than AT and IC oligonucleotides. When S- and R-ofloxacin are compared, the S-enantiomer exhibits more favorable total energy and torsion angles in the complex formation. This result is in agreement with the experimental observation [Hwangbo et al., Eur J Pharm Sci 18, 197 (2003)]. In the complex, both enantiomers form two hydrogen bonds: one between the carbonyl group of ofloxacin and the amine group of G16 and the other between the fluorine group and the G6 amine for S-ofloxacin. However, only one hydrogen bond is formed between endocyclic hydrogen atom at the C2 position of adenine and inosine base and carbonyl group of ofloxacin, which may be the reason for the GC preferentiality of ofloxacin.

Key words: Molecular dynamics; Ofloxacin; Quinolone; DNA; Binding Mode; and Enantioselective binding.