Molecular dynamics simulations provide a structural basis for the experimentally observed nucleotide preferences for DNA interstrand cross-links induced by aziridinylbenzoquinones

Spermine inhibition of the 2,5-diaziridinyl-1,4-benzoquinone (DZQ) crosslinking reaction with DNA duplexes containing poly(purine). poly(pyrimidine) tracts

Upon reduction, 2,5-diaziridinyl-1,4-benzoquinone (DZQ) can form an interstrand guanine to guanine crosslink with DNA duplexes containing a d(GC).d(GC) dinucleotide step. The reaction is enhanced by a thymine positioned 5[prime] to each guanine [i.e. in a d(TGCA). d(TGCA) duplex fragment]. Here we show that spermine can inhibit DZQ crosslink formation in duplexes of sequence d[C(N6)TGCA(M6)C]. d[G(M[prime]6)TG-CA(N[prime]6)G]. For N6= M6= GGGGGG, N6= M6= a 'random' sequence and N6= GGGGGG and M6= a 'random' sequence, spermine concentrations of 20, 1 and 3 microM, respectively, were required for 50% inhibition of the DZQ crosslink. This suggests that spermine is more strongly bound to the polyguanosine tract than the random sequence, making it less available for crosslink inhibition. When the polyguanosine tract is interrupted by N 7-deazaguanine (D) located three bases, d(CGGGDGGTGCAGGDGGGC), and four bases, d(CG-GDGGGTGCAGGGDGGC), from the d(TGCA).d(TGCA) site, 30 and 3 microM spermine, respectively, were required for 50% crosslink inhibition. We suggest that this difference is due to the relative proximity of the three-guanosine tract to the d(TGCA).d(TGCA) site. We were able to confirm these conclusions with further experiments using duplexes containing three-guanosine and two-guanosine tracts and from computer simulations of the spermine-DNA complexes.

Specificity and kinetics of interstrand and intrastrand bifunctional alkylation by nitrogen mustards at a G-G-C sequence

Previous work showed that melphalan-induced mutations in the aprt gene of CHO cells are primarily transversions and occur preferentially at G-G-C sequences, which are potential sites for various bifunctional alkylations involving guanine N-7. To identify the DNA lesion(s) which may be responsible for these mutations, an end-labeled DNA duplex containing a frequent site of melphalan-induced mutation in the aprt gene was treated with melphalan, mechlorethamine or phosphoramide mustard. The sequence specificity and kinetics of formation of both interstrand and intrastrand crosslinks were determined. All mustards selectively formed two base-staggered interstrand crosslinks between the 5'G and the G opposite C in the 5'G-G-C sequence. Secondary alkylation was much slower for melphalan than for the other mustards and the resulting crosslink was more stable. Mechlorethamine and phosphoramide mustard induced intrastrand crosslinks between the two contiguous Gs in the G-G-C sequence in double-stranded DNA, but melphalan did not. Molecular dynamic simulations provided a structural explanation for this difference, in that the monofunctionally bound intermediates of mechlorethamine and phosphoramide mustard assumed thermodynamically stable conformations with the second arm in a position appropriate for intrastrand crosslink formation, while the corresponding melphalan monoadduct did not.

Molecular dynamics simulations of chlorambucil/DNA adducts. A structural basis for the 5'-GNC interstrand DNA crosslink formed by nitrogen mustards

The alkylation of DNA by chlorambucil has been studied using a computational approach. Molecular dynamics simulations were performed on the fully solvated non-covalent complex, two monoadducts and a crosslinked diadduct of chlorambucil with the d(CGG3G2CGC).-d(GCG1CCCG) duplex, in which the N7 atoms of G1, G2 and G3 are potential alkylation sites. The results provide a structural basis for the preference of nitrogen mustards to crosslink DNA duplexes at a 5'-GNC site (a 1,3 crosslink, G1-G3) rather than at a 5'-GC sites (a 1,2 crosslink, G1-G2). In the non-covalent complex simulation the drug reoriented from a non-interstrand crosslinking location to a position favorable for G1-G3 diadduct formation. It proved possible to construct a G1-G3 diadduct from a structure from the non-covalent simulation, and continue the molecular dynamics calculation without further disruption of the DNA structure. A crosslinked diadduct developed with four BII conformations on the 3' side of each alkylated guanine and of their respective complementary cytosine. In the first monoadduct simulation the starting point was the same DNA conformation used in the crosslinked diadduct simulation with alkylation at G1. In this simulation the DNA deformation was reduced, with the helix returning to a more canonical form. A second monoadduct simulation was started from a canonical DNA conformation alkylated at G3. Here, no significant motion towards a potential crosslinking conformation occurred. Collectively, the results suggest that crosslink formation is dependent upon the drug orientation prior to alkylation and the required deformation of the DNA to permit 1,3 crosslinking can largely be achieved in the non-covalent complex.