Investigation of interactions of doxorubicin with purine nucleobases by molecular modeling

Doxorubicin, an anthracycline antibiotic with anti-tumor activity, is produced by the bacterium Streptomyces peucetius. The interactions between doxorubicin and genetic material and the details of the intercalation with DNA have been controversial issues. Thus, the interactions of doxorubicin with purine nucleobases were studied by quantum mechanical methods. Initially, conformer analyses of doxorubicin were performed with Spartan 08 software and 319 different conformers from 422 initial structures for doxorubicin were obtained. Geometry optimizations and frequency analyses were performed for each structure using density functional theory (DFT) at B3LYP/6-31G** level using Gaussian 09 software. The most stable 20 conformers of doxorubicin and tautomers of purine nucleobases were optimized again with ɷB97XD/6-31G** level and their interactions were also analyzed at the same level. The Discovery Studio 3.5 Visualizer was used to draw the initial and optimized structures of investigated geometries. The noncovalent interactions (NCIs) were visualized by calculating reduced density gradient (RDG) with Multiwfn program. The color-filled isosurfaces and RDG scatter maps of most stable interaction geometries were plotted by Visual Molecular Dynamics (VMD) software and Gnuplot 5.3 software, respectively. This study showed that adenine, guanine, and hypoxanthine nucleobases interact with doxorubicin by forming strong hydrogen bonds and π-π interactions. Considering the normal cellular conditions, the effect of solvent (water) on the interaction geometries were also analyzed and when compared to gas phase it was determined that the movements of the molecules were restricted and there was a minimal change between initial and optimized structures in the aqueous phase.