Density functional and Ab initio studies on N-acetylduocarmycin SA: insight into its DNA interaction properties

Duocarmycins binding to DNA investigated by molecular simulation

Duocarmycins are a potent class of antitumor agents, whose activity arises through their covalent binding to adenine nucleobases of DNA.(1-3) Here, we perform molecular dynamics (MD) and hybrid Car-Parinello QM/MM simulations to investigate aspects of duocarmycin binding to the d(pGpApCpTpApApTpTpGpApC) oligonucleotide. We focus on the derivatives (+)-duocarmycin SA (DSA) and (+)-duocarmycin SI (DSI), for which structural information of the covalent complex with the oligonucleotide is available, as well as on the related, but less reactive, NBOC-duocarmycin SA (NBOC-DSA), interacting with the same oligonucleotide. Comparison is made with adenine alkylation reaction in water performed by the smallest of these compounds (NBOC-DSA). The MD calculations suggest that, in noncovalent complexes, (i) drug binding causes a partial dehydration of the minor groove, without inducing a significant conformational changes, and (ii) DSA and DSI occupy a more favorable position for nucleophilic attack than NBOC-DSA, consistently with the lower reactivity of the latter. The QM/MM calculations, which are used to investigate the first step of the alkylation reaction, turn out to provide strongly underestimated free energy barriers. Within these approximations, our calculations suggest that an important ingredient for the experimentally observed DNA catalytic power is the polarization of the drug by the biomolecular scaffold.

On the role of stereo-electronic effects in tuning the selectivity and rate of DNA alkylation by duocarmycins

The role of local geometric and stereo-electronic effects in tuning the alkylation of DNA by duocarmycins has been analyzed by an integrated computational tool rooted in the density functional theory and the polarizable continuum model. Our study points out that together with steric accessibility, different electronic delocalisations also contribute to determine the higher reactivity of adenine with respect to guanine. Also the effect of the methyl ester group on the alkylating agent has an electronic origin. Furthermore, deviations from the planarity in the drug structure (conformational catalysis) could be less important than currently accepted since, according to our computations, compounds with strongly different reactivity have nearly constant and very similar out of plane distortions before and after the reaction. Model computations suggest, instead, that specific non covalent interactions could discriminate between different drugs selectively reducing some activation energies with respect to the corresponding processes in solution.

Modeling Substituent and Conformational Effects on the Reactivity of Antitumor Agents Containing a Cyclopropylcyclohexadienone Subunit

The uncatalyzed alkylation reactions of ammonia by the parent spirocyclopropylcyclohexadienone (6), its 3-amino analogue (7), the cyclic derivative (8), its N-formyl derivative (9), and a closer model (10) of the CPI (1-4) drugs have been investigated in gas phase and in water solvent bulk, using density functional theory at the B3LYP level with several basis sets and the C-PCM solvation model. The effect of several structural key features such as the vinylogous amide conjugation, the acylation of the 2-amino substituent, the ring constraint of the heterocyclic nitrogen atom at C(2) carbon in a ring, and the presence of a condensed pyrrole ring on the reaction activation energy have been investigated. Substrate 7, which is a flexible conformational model of the cyclopropylpyrroloindole moiety (CPI) contained in the duocarmycins, has been used to model the shape-dependent reactivity of these drugs, in gas phase and water solutions. The calculations indicate that shape dependence of reactivity is strongly operative both in gas phase and in polar solvents, since conformational effects are capable of reducing the reaction activation energy by -8.4 and -4.3 kcal mol(-1) in gas phase and in water solution, respectively, that is required to promote "conformational catalysis".

NMR Structure of the (+)-CPI-indole/d(GACTAATTGAC)-d(GTCAATTAGTC) Covalent Complex

We report the NMR solution structure of (+)-CPI-indole (CPI, 1,2,8,8a-tetrahydrocyclopropa[c]pyrrolo[3,2-e]indol-4(5H)-one), an agent belonging to the CC-1065/duocarmycin family of antitumor compounds. This (+)-CPI-indole structure is covalently bound to d(G(1)ACTAATTGTC(11))-d(G(12)TCAATTAGTC(22)), a synthetic DNA duplex containing a high-affinity binding site. The three-dimensional structure has been determined by several cycles of restrained molecular dynamics calculations with a total of 563 NMR-derived constraints, both in vacuo and by using the generalized Born solvent continuum model. In-depth analysis of the structure of this ligand-DNA complex led to a detailed knowledge of the bound state conformation of the CPI-indole, the most simplified agent related to CC-1065 and duocarmycins, the parent members of a family of extremely potent antitumor compounds. Comparison of the CPI-indole bound conformation with those previously found for (+)-duocarmycin SA (DSA), its unnatural enantiomer (-)-DSA, and the demethoxylated analogue (+)-DSI in their DNA complexes provided additional evidence of the tight correlation between the catalytic effect exerted by DNA on the alkylation reaction and the extent of angular twist between the two planar heteroaromatic subunits of these agents. Additionally, comparison of the structural features of the DNA-bound state of a "naked" ligand, such as CPI-indole, with those of various other duocarmycin agents provided useful information for the interpretation of the observed effects on chemical reactivity of the different substitution patterns at the hemispheres of these types of complex.

Effect of protonation of the N-acetyl neuraminic acid residue of sialyl Lewis(X): a molecular orbital study with insights into its binding properties toward the carbohydrate recognition domain of E-selectin

Semiempirical molecular orbital (MO) calculations with an implicit treatment of the water environment were employed in order to assess whether the sialyl Lewis(X) (sLe(X)) tetrasaccharide binds to E-selectin in the anionic or neutral (i.e., protonated) state. The analysis of the frontier molecular orbitals, electrostatic potential surfaces, and conformational behavior of the sugar indicates that its neutral form possesses the necessary characteristics for binding. In particular, the LUMO level of the neutral sLe(X) molecule is localized both on the carboxylic group of the N-acetyl neuraminic acid (NeuNAc) residue and on the nearby glycosidic linkage. These two moieties interact with the Arg97 residue of E-selectin, as revealed by a recent crystal structure analysis of the E-selectin/sLe(X) complex. The energetics of this specific interaction was investigated with the aid of ab initio Hartree-Fock MO calculations, which resulted in a BSSE-corrected binding energy of 16.63 kcal/mol. Our observations could open up new perspectives in the design of sLe(X) mimics.