Chloroethylnitrosoureas (CENUs) are an important family of alkylating agents used in the clinical treatment of cancer. Their anticancer mechanism primarily involves the formation of DNA interstrand crosslinks (ICLs) induced by the chloroethyldiazonium ion derived from the decomposition of CENUs. In this work, the mechanism for the formation of ICLs was investigated by density functional theory (DFT) with B3LYP, wB97XD, and M062X functinoals using conductor-like polarizable continuum model solvent model. Three pathways leading to the formation of three types of G-C crosslinks were compared. G(N1)-C(N3) crosslink is predicted to be the dominant crosslinking product other than G(O 6 )-C(N 4 ) and G(N 2 )-C(O 2 ) crosslinks, which is consistent with the previous results obtained from QM/MM computations. The results indicate that the formation of the G(N1)-C(N3) crosslink via pathway A is the most favorable mechanism from both kinetic and thermodynamic standpoints. In this pathway, the chloroethyldiazonium ion alkylates guanine on the O 6 site followed by intramolecular cyclization to form O 6 ,N1ethanoguanine (4). The cytosine then reacts with intermediate 4 on the C a atom to yield the G(N1)-C(N3) crosslink. This work provides reasonable explanations for the supposed mechanism of CENUs-induced ICLs formation obtained from experimental investigations.the reaction was plotted to confirm the rate-determining step and to provide explanations for the experimental observations.
Models and ComputationsAs demonstrated in previous research efforts, [20] DNA ICLs induced by N-nitroso compounds (NNCs) were formed between complementary guanine and cytosine, because the ethidenes generated from NNCs exactly matched with the distance between the paired negative atoms. Therefore, three pathways (see Fig. 2) leading to three kinds of DNA ICLs, that is, G(O 6 )-C(N 4 ), G(N1)-C(N3), and G(N 2 )-C(O 2 ), were compared by DFT calculations. The geometric structures of all reactants, transition states (TSs), intermediates, and products were optimized with the DFT methods, using B3LYP, [21,22] wB97XD, [23] and M062X [24] functionals, respectively, with a 6-31þG(d,p) basis set. As all DNA damage events involved in anticancer or carcinogenic processes were supposed to take place in aqueous solution, the solvent effect of water on all reactions was taken into account. On the basis of the optimized geometries obtained in the gas phase, all structures were further optimized using self-consistent reaction field computations using the conductor-like polarizable continuum model (CPCM) [25,26] using the B3LYP, wB97XD, and M062X functionals, respectively, at the 6-31þG(d,p) theoretical level. Universal force field radii was used for all structures to evaluate the solvent effects of water, which is the default for the Gaussian 09 program package. The TSs were located by the Synchronous Transit-Guided Quasi-Newton method [27,28] with the QST2 or QST3 options to the Opt keyword. The analytical computations of the vibrational frequencies were performe...