Bifunctional binding of cisplatin to DNA: Why does cisplatin form 1,2-intrastrand cross-links with AG but not with GA?

Abstract

The bifunctional binding of the anticancer drug cisplatin to two adjacent nucleobases in DNA is modeled using density functional theory. Previous experimental studies revealed that cisplatin binding to adjacent guanine and adenine is sensitive to nucleobase sequence. Whereas AG 1,2-intrastrand cross-links are commonly observed, the analogous GA adducts are not known. This study focuses on understanding this directional preference by constructing a full reaction profile using quantum chemical simulation methods. Monofunctional and bifunctional cisplatin adducts were generated, and the transition states that connect them were located for the dinucleotides d(pApG) and d(pGpA), assuming that initial platination takes place at the guanine site. Our computer simulations reveal a significant kinetic preference for formation of the AG over the GA adduct. The activation free energies of similar to 23 kcal/mol for AG and similar to 32 kcal/mol for GA suggest that bifunctional closure is similar to 6 orders of magnitude faster for AG than for GA. A strong hydrogen bond between one of the ammine ligands of cisplatin and the 5' phosphate group of the DNA backbone is responsible for the stabilization of the transition state that affords the AG adduct. This interaction is absent in the transition state that leads to the GA adduct because the right-handed helix of the DNA backbone places the phosphate out of reach for the ammine ligand. We found only an insignificant thermodynamic difference between AG and GA adducts and conclude that the preference of AG over GA binding is largely under kinetic control. The puckering of the deoxyribose ring plays an important role in determining the energetics of the bifunctional platination products. Whereas the 3'-nucleoside remains in the native C2'-endo/C3'-exo form of B-DNA, the deoxyribose of the 5'-nucleoside always adopts the C2'-exo/C3'-endo puckering in our simulations. A detailed analysis of the energies and structures of the bifunctional adducts revealed that the observed sugar puckering patterns are necessary for platinum to bind in a relaxed coordination geometry.