Restrained molecular dynamics studies on complex of adriamycin with DNA hexamer sequence d-CGATCG

Abstract

Adriamycin is an anthracycline anticancer drug used widely for solid tumors in spite of its adverse side effects. The solution structure of 2:1 adriamycin-d-(CGATCG)₂ complex has been studied by restrained molecular dynamics simulations. The restraint data set consists of several intramolecular and intermolecular nuclear Overhauser enhancement cross-peaks obtained from two-dimensional nuclear magnetic resonance spectroscopy data. The drug is found to intercalate between CG and GC base pairs at two d-CpG sites. The drug-DNA complex is stabilized via specific hydrogen bonding and van der Waal's interactions involving 4OCH₃, O5, 6OH, and NH₃⁺ moiety of daunosamine sugar, and rings A protons. The O-glycosidic bond C7-O7-C1'-C2' lies in the range 138°-160° during the course of simulations. The O6-H6...O5 hydrogen bond is stable while O11-H11...O12 hydrogen bond is not favored. The intercalating base pairs are buckled and minor groove is wider in the complex. The phosphate on one strand at intercalation site C1pG2 is in B_I conformation and the phosphates directly lying on opposite strand is in B_II conformation. The phosphorus on adjacent site G2pA3 is in B_II conformation and hence a distinct pattern of B_I and B_II conformations is induced and stabilized. The role of various functional groups by which the molecular action is mediated has been discussed and correlated to the available biochemical evidence.