Molecular models of neocarzinostatin damage of DNA: analysis of sequence dependence in 5'GAGCG:5'CGCTC.

Model building and molecular mechanics and dynamics calculations have been performed on a number of complexes of the post-activated form of the neocarzinostatin chromophore (NCS) with the B-DNA oligomer 5'GAGCG:5'CGCTC. Stable structures with the naphthoic acid moiety intercalated at all base pairs can be constructed. The observed bistranded lesions consisting of an abasic site at the Cyt residue in AGC and a direct break at the Thy residue on the complementary strand can be explained by assuming that NCS in the (R,R) form intercalates between the Ade2-Thy9/Gua3-Cyt8 base step with its 'diradical' core oriented towards the 3'-end of the (+) strand. Sites at C5', C4' and C1' in the minor groove are within a short enough distance from the two radical centers on NCS to permit hydrogen atom abstraction and the formation of the bistranded lesions. Strand cleavage at Thy9 may occur as a single lesion if NCS is intercalated into the Gua3-Cyt8/Cyt4-Gua7 base step with its active core towards the 3'-end of the (-) strand. The results are analyzed, and the utility and limitations of this type of model building are discussed.