High resolution solution structure of a DNA duplex alkylated by the antitumor agent duocarmycin SA

Single Mutation on Trastuzumab Modulates the Stability of Antibody-Drug Conjugates Built Using Acetal-Based Linkers and Thiol-Maleimide Chemistry

Antibody-drug conjugates (ADCs) are a class of targeted therapeutics used to selectively kill cancer cells. It is important that they remain intact in the bloodstream and release their payload in the target cancer cell for maximum efficacy and minimum toxicity. The development of effective ADCs requires the study of factors that can alter the stability of these therapeutics at the atomic level. Here, we present a general strategy that combines synthesis, bioconjugation, linker technology, site-directed mutagenesis, and modeling to investigate the influence of the site and microenvironment of the trastuzumab antibody on the stability of the conjugation and linkers. Trastuzumab is widely used to produce targeted ADCs because it can target with high specificity a receptor that is overexpressed in certain breast cancer cells (HER2). We show that the chemical environment of the conjugation site of trastuzumab plays a key role in the stability of linkers featuring acid-sensitive groups such as acetals. More specifically, Lys-207, located near the reactive Cys-205 of a thiomab variant of the antibody, may act as an acid catalyst and promote the hydrolysis of acetals. Mutation of Lys-207 into an alanine or using a longer linker that separates this residue from the acetal group stabilizes the conjugates. Analogously, Lys-207 promotes the beneficial hydrolysis of the succinimide ring when maleimide reagents are used for conjugation, thus stabilizing the subsequent ADCs by impairing the undesired retro-Michael reactions. This work provides new insights for the design of novel ADCs with improved stability properties.

Structural evolution of a DNA repair self-resistance mechanism targeting genotoxic secondary metabolites

Microbes produce a broad spectrum of antibiotic natural products, including many DNA-damaging genotoxins. Among the most potent of these are DNA alkylating agents in the spirocyclopropylcyclohexadienone (SCPCHD) family, which includes the duocarmycins, CC-1065, gilvusmycin, and yatakemycin. The yatakemycin biosynthesis cluster in Streptomyces sp. TP-A0356 contains an AlkD-related DNA glycosylase, YtkR2, that serves as a self-resistance mechanism against yatakemycin toxicity. We previously reported that AlkD, which is not present in an SCPCHD producer, provides only limited resistance against yatakemycin. We now show that YtkR2 and C10R5, a previously uncharacterized homolog found in the CC-1065 biosynthetic gene cluster of Streptomyces zelensis, confer far greater resistance against their respective SCPCHD natural products. We identify a structural basis for substrate specificity across gene clusters and show a correlation between in vivo resistance and in vitro enzymatic activity indicating that reduced product affinity-not enhanced substrate recognition-is the evolutionary outcome of selective pressure to provide self-resistance against yatakemycin and CC-1065.

Inverse Electron Demand Diels-Alder Reactions of Heterocyclic Azadienes, 1-Aza-1,3-Butadienes, Cyclopropenone Ketals, and Related Systems. A Retrospective

A summary of the investigation and applications of the inverse electron demand Diels-Alder reaction is provided that have been conducted in our laboratory over a period that now spans more than 35 years. The work, which continues to provide solutions to complex synthetic challenges, is presented in the context of more than 70 natural product total syntheses in which the reactions served as a key strategic step in the approach. The studies include the development and use of the cycloaddition reactions of heterocyclic azadienes (1,2,4,5-tetrazines; 1,2,4-, 1,3,5-, and 1,2,3-triazines; 1,2-diazines; and 1,3,4-oxadiazoles), 1-aza-1,3-butadienes, α-pyrones, and cyclopropenone ketals. Their applications illustrate the power of the methodology, often provided concise and nonobvious total syntheses of the targeted natural products, typically were extended to the synthesis of analogues that contain deep-seated structural changes in more comprehensive studies to explore or optimize their biological properties, and highlight a wealth of opportunities not yet tapped.

The Difference a Single Atom Can Make: Synthesis and Design at the Chemistry-Biology Interface

A Perspective of work in our laboratory on the examination of biologically active compounds, especially natural products, is presented. In the context of individual programs and along with a summary of our work, selected cases are presented that illustrate the impact single atom changes can have on the biological properties of the compounds. The examples were chosen to highlight single heavy atom changes that improve activity, rather than those that involve informative alterations that reduce or abolish activity. The examples were also chosen to illustrate that the impact of such single-atom changes can originate from steric, electronic, conformational, or H-bonding effects, from changes in functional reactivity, from fundamental intermolecular interactions with a biological target, from introduction of a new or altered functionalization site, or from features as simple as improvements in stability or physical properties. Nearly all the examples highlighted represent not only unusual instances of productive deep-seated natural product modifications and were introduced through total synthesis but are also remarkable in that they are derived from only a single heavy atom change in the structure.

Toxicity and repair of DNA adducts produced by the natural product yatakemycin

Yatakemycin (YTM) is an extraordinarily toxic DNA alkylating agent with potent antimicrobial and antitumor properties and the most recent addition to the CC-1065 and duocarmycin family of natural products. While bulky DNA lesions the size of those produced by YTM are normally removed from the genome by the nucleotide excision repair (NER) pathway, YTM adducts are also a substrate for the bacterial DNA glycosylases AlkD and YtkR2, unexpectedly implicating base excision repair (BER) in their elimination. The reason for the extreme toxicity of these lesions and the molecular basis for how they are eliminated by BER have been unclear. Here, we describe the structural and biochemical properties of YTM adducts responsible for their toxicity, and define the mechanism by which they are excised by AlkD. These findings delineate an alternative strategy for repair of bulky DNA damage and establish the cellular utility of this pathway relative to that of NER.

A fundamental relationship between hydrophobic properties and biological activity for the duocarmycin class of DNA-alkylating antitumor drugs: hydrophobic-binding-driven bonding

Two systematic series of increasingly hydrophilic derivatives of duocarmycin SA that feature the incorporation of ethylene glycol units (n = 1-5) into the methoxy substituents of the trimethoxyindole subunit are described. These derivatives exhibit progressively increasing water solubility along with progressive decreases in cell growth inhibitory activity and DNA alkylation efficiency with the incremental ethylene glycol unit incorporations. Linear relationships of cLogP with -log IC50 for cell growth inhibition and -log AE (AE = cell-free DNA alkylation efficiency) were observed, with the cLogP values spanning the productive range of 2.5-0.49 and the -log IC50 values spanning the range of 11.2-6.4, representing IC50 values that vary by a factor of 10(5) (0.008 to 370 nM). The results quantify the fundamental role played by the hydrophobic character of the compound in the expression of the biological activity of members in this class (driving the intrinsically reversible DNA alkylation reaction) and define the stunning magnitude of its effect.

Synthesis and evaluation of a series of C5'-substituted duocarmycin SA analogs

The synthesis and evaluation of a key series of analogs of duocarmycin SA, bearing a single substituent at the C5' position of the DNA binding subunit, are described.

Investigation of the transformations of a novel anti-cancer agent combining HPLC, HPLC-MS and direct ESI-HRMS analyses

One of the main problems of anti-cancer therapy is an insufficient differentiation between normal and malignant cells by the known anti-proliferant agents. The antibody-directed enzyme prodrug therapy is a promising approach for a selective treatment of cancer, in which a non-toxic prodrug is enzymatically converted into a highly cytotoxic drug at the surface of malignant cells by a targeted antibody-enzyme conjugate. The transformations and the stability of a very promising novel prodrug and its corresponding cytotoxic derivative were now investigated in detail by high-performance liquid chromatography (HPLC)-mass spectrometry (MS). In order to determine the time-dependent DNA alkylation efficiency and the sequence selectivity of the novel compounds, DNA binding studies using direct electrospray-Fourier transform ion cyclotron resonance-MS (ESI-FTICR-MS) have been performed. These measurements were accompanied by HPLC analyses followed by MS of the separated species to confirm the results of the direct ESI-FTICR-MS measurements. The sites of DNA alkylation could be identified unambiguously by the mass spectrometric fragmentation pattern of the alkylated oligodeoxynucleotides as well as by the results of HPLC followed by MS. A combination of all techniques applied led to a better understanding of the mode of action of the new therapeutics and might be used for an estimation of the cytotoxicity of different prodrugs and drugs since the alkylation efficiency correlates with the bioactivity of the compounds in cell culture investigations.

Sliding of alkylating anticancer drugs along the minor groove of DNA: new insights on sequence selectivity

Currently, little is known about the molecular recognition pathways between DNA-alkylating anticancer drugs and their targets despite their pharmacological relevance. In the framework of classical molecular dynamics simulations, here we use umbrella sampling to map the potential of mean force (PMF) associated with sliding along the DNA minor groove of two of these compounds. These are an indole derivative of duocarmycin (DSI) and the putative reactive form of anthramycin (anhydro-anthramycin, IMI). Twenty-three configurations were considered for each drug/DNA complex, corresponding to a movement along approximately 3 basepairs. The alkylation site turns out to be the most favorable for DSI, while a barrier of approximately 6 kcal/mol separates the reactive configuration of IMI.DNA from the absolute minimum. An analysis of various contributions to the PMF reveals that solvent effects play an important role for the largest and more flexible drug DSI. Instead, the PMF of IMI.DNA overall correlates with changes in the binding enthalpy. Implications of these results on the sequence selectivity of the two drugs are discussed.

Binding free energy calculation for duocarmycin/DNA complex based on the QPLD-derived partial charge model

The 3ns unrestrained MD simulations were carried out on the DNA/duocarmycin complex based on (1) the classic RESP charge model, and (2) the QM-polarized ligand docking (QPLD)-based charge model. The RMSDs of the trajectories and the DeltaG(bind) of the QPLD model perform much better than the RESP model, with the DeltaG(bind) estimation for QPLD model (-16.11 kcal/mol) versus DeltaG(bind) estimation for RESP model (-10.05 kcal/mol).

Modeling anticancer drug–DNA interactions via mixed QM/MM molecular dynamics simulations

The development of anticancer drugs started over four decades ago, with the serendipitous discovery of the antitumor activity of cisplatin and its successful use in the treatment of various cancer types. Despite the efforts made in unraveling the mechanism of the action of cisplatin, as well as in the rational design of new anticancer compounds, in many cases detailed structural and mechanistic information is still lacking. Many of these drugs exert their anticancer activity by covalently binding to DNA inducing a distortion or simply impeding replication, thus triggering a cellular response, which eventually leads to cell death. A detailed understanding of the structural and electronic properties of drug-DNA complexes and their mechanism of binding is the key step in elucidating the principles of their anticancer activity. At the theoretical level, the description of covalent drug-DNA complexes requires the use of state-of-the-art computer simulation techniques such as hybrid quantum/classical molecular dynamics simulations. In this review we provide a general overview on: drugs which covalently bind to DNA duplexes, the basic concepts of quantum mechanics/molecular mechanics (QM/MM), molecular dynamics methods and a list of selected applications of these simulations to the study of drug-DNA adducts. Finally, the potential and the limitations of this approach to the study of such systems are critically evaluated.

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 Acid and Cationic Catalysis on the Reactivity of Duocarmycins

Several catalyzed alkylation reactions of 9-methyladenine by a model [CPI, cyclopropa[c]pyrrolo[3,2-e]indol-4(5H)-one (1)] of duocarmycin anticancer drugs have been compared to the uncatalyzed reaction in gas phase and in water solvent bulk, using density functional theory at the B3LYP level with the 6-31+G(d,p) basis set and C-PCM solvation model. The effect on the CPI reactivity induced by water, formic and phosphoric acids (general acid catalysis), H3O+ (specific acid catalysis), sodium, and ammonium cation complexation (cationic catalysis) has been investigated. The calculations indicate that the specific acid catalysis and the catalysis induced by sodium cation complexation are strong in the gas phase, but solvation reduces them dramatically by electrostatic effects. The specific acid catalysis is still operative, but strongly reduced in water solution, where the reaction barrier is reduced by 8.6 kcal mol(-1) in comparison to the uncatalyzed reaction. The general acid catalysis induced by phosphoric acid (-7.3 kcal mol(-1)) and the catalysis induced by Na+ and NH4+ complexation become competitive, with a catalytic effect of -3.6 and -4.1 kcal mol(-1) in water, respectively. With the specific acid catalysis, the high acidity (low pK(a) value) of the conjugated acid of CPI (CPIH+), computed in water solution using both C-PCM (pK(a) = +2.6) and PCM-B3LYP/6-31+G(d,p) (pK(a) = +2.4) solvation models, suggests that the catalytic effects induced by NH4+ complexation could become more important than the specific acid catalysis and the general catalysis by H3PO4 under physiological conditions, due to concentration effects of the catalysts.

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.

Synthesis and Evaluation of Duocarmycin and CC-1065 Analogues Incorporating the 1,2,9,9a-Tetrahydrocyclopropa[c]benz[e]-3-azaindol-4-one (CBA) Alkylation Subunit

An efficient eight-step synthesis (53% overall) and the evaluation of 1,2,9,9a-tetrahydrocyclopropa[c]benz[e]-3-azaindol-4-one (CBA) and its derivatives containing an aza variant of the CC-1065/duocarmycin alkylation subunit are detailed. This unique deep-seated aza modification provided an unprecedented 2-aza-4,4-spirocyclopropacyclohexadienone that was characterized chemically and structurally (X-ray). CBA proved structurally identical with CBI, the carbon analogue, including the stereoelectronic alignment of the key cyclopropane, its bond lengths, and the bond length of the diagnostic C3a-N2 bond, reflecting the extent of vinylogous amide (amidine) conjugation. Despite these structural similarities, CBA and its derivatives were found to be much more reactive toward solvolysis and hydrolysis, much less effective DNA alkylating agents (1000-fold), and biologically much less potent (100- to 1000-fold) than the corresponding CBI derivatives.

DNA alkylation properties of yatakemycin

Yatakemycin represents the newest and now most potent member of a class of naturally occurring antitumor compounds that includes CC-1065 and the duocarmycins, which derive their biological properties from a characteristic DNA alkylation reaction. Herein, the first description of the yatakemycin DNA alkylation properties is detailed, constituting the first such study of a naturally occurring "sandwiched" member of this class. Thus, the event, sequence selectivity, relative rate and efficiency, and reversibility of the DNA alkylation reaction of yatakemycin are described.

Establishment of substituent effects in the DNA binding subunit of CBI analogues of the duocarmycins and CC-1065

An extensive series of CBI analogues of the duocarmycins and CC-1065 exploring substituent effects within the first indole DNA binding subunit is detailed. In general, substitution at the indole C5 position led to cytotoxic potency enhancements that can be >/=1000-fold providing simplified analogues containing a single DNA binding subunit that are more potent (IC(50)=2-3 pM) than CBI-TMI, duocarmycin SA, or CC-1065.

Structure of a beta-alanine-linked polyamide bound to a full helical turn of purine tract DNA in the 1:1 motif

Polyamides composed of N-methylpyrrole (Py), N-methylimidazole (Im) and N-methylhydroxypyrrole (Hp) amino acids linked by beta-alanine (beta) bind the minor groove of DNA in 1:1 and 2:1 ligand to DNA stoichiometries. Although the energetics and structure of the 2:1 complex has been explored extensively, there is remarkably less understood about 1:1 recognition beyond the initial studies on netropsin and distamycin. We present here the 1:1 solution structure of ImPy-beta-Im-beta-ImPy-beta-Dp bound in a single orientation to its match site within the DNA duplex 5'-CCAAAGAGAAGCG-3'.5'-CGCTTCTCTTTGG-3' (match site in bold), as determined by 2D (1)H NMR methods. The representative ensemble of 12 conformers has no distance constraint violations greater than 0.13 A and a pairwise RMSD over the binding site of 0.80 A. Intermolecular NOEs place the polyamide deep inside the minor groove, and oriented N-C with the 3'-5' direction of the purine-rich strand. Analysis of the high-resolution structure reveals the ligand bound 1:1 completely within the minor groove for a full turn of the DNA helix. The DNA is B-form (average rise=3.3 A, twist=38 degrees ) with a narrow minor groove closing down to 3.0-4.5 A in the binding site. The ligand and DNA are aligned in register, with each polyamide NH group forming bifurcated hydrogen bonds of similar length to purine N3 and pyrimidine O2 atoms on the floor of the minor groove. Each imidazole group is hydrogen bonded via its N3 atom to its proximal guanine's exocyclic amino group. The important roles of beta-alanine and imidazole for 1:1 binding are discussed.

The DNA phosphate backbone is not involved in catalysis of the duocarmycin and CC-1065 DNA alkylation reaction

The rates of DNA alkylation were established for the reaction of (+)-duocarmycin SA (1) with the native duplex d(G(1)TCAATTAGTC(11))*d(G(12)ACTAATTGAC(22)), an 11 bp deoxyoligonucleotide that contains a single high-affinity alkylation site that has been structurally characterized at exquisite resolution, and modified duplexes in which the four backbone phosphates proximal to the C4 carbonyl of bound 1 were replaced with methylphosphonates. All were found to react at comparable rates establishing that these backbone phosphates do not participate in catalysis of the DNA alkylation reaction.

Substituent effects within the DNA binding subunit of CBI analogues of the duocarmycins and CC-1065

A series of CBI analogues of the duocarmycins and CC-1065 exploring substituent effects within the first indole DNA binding subunit are detailed. Substitution at the indole C5 position led to cytotoxic potency enhancements that are > or =1000-fold, providing simplified analogues containing a single DNA binding subunit that are more potent (IC(50)=2-3 pM) than CBI-TMI, duocarmycin SA, or CC-1065.

Synthesis, chemical properties, and biological evaluation of CC-1065 and duocarmycin analogues incorporating the 5-methoxycarbonyl-1,2,9,9a-tetrahydrocyclopropa

The synthesis of 5-methoxycarbonyl-1,2,9,9a-tetrahydrocyclopropa[c]benz[e]indol-4-one (C5-CO2Me-CBI), a substituted CBI derivative bearing a C5 methoxycarbonyl group, and its corresponding 5-hydroxymethyl derivative are described in efforts to establish substituent electronic effects on the agents' functional reactivity and the resulting effect this has on their rate of DNA alkylation. Resolution of an immediate C5-CO2Me-CBI precursor and its incorporation into both enantiomers of 16 and 17, analogues of the duocarmycins, are also detailed. A study of the solvolysis reactivity and regioselectivity of N-BOC-C5-CO2Me-CBI (12) revealed that the introduction of a C5 methyl ester modestly slowed the rate of solvolysis (1.8x, pH 3) without altering the inherent reaction regioselectivity (>20:1). The comparison of the X-ray structures of the N-CO2Me derivatives of C5-CO2Me-CBI and CBI revealed correlations with the reaction regioselectivity and the relative reactivity of the compounds. The latter correlated well with the less reactive C5-CO2Me-CBI exhibiting a shortened N2-C2a bond length (1.386 vs 1.390 A) and smaller chi1 dihedral angle (8.1 degrees vs 21.2 degrees ) indicative of greater vinylogous amide conjugation and was accompanied by a diminished (cross-conjugated) cyclopropane conjugation (shorter bond lengths). Establishment of the DNA alkyation properties revealed that C5-CO2Me-CBI-based agents retained the identical alkylation selectivity of the natural products. More importantly, the C5 methyl ester was found to decrease the rate (0.77x) of DNA alkylation relative to CBI, consistent with its inherent lower reactivity. These results indicate that the previously observed increase in the rate of DNA alkylation for C7-substituted CBI analogues including CCBI (7-cyano-CBI) is contrary to expectations based on their inherent reactivities. Unlike 17, in which the C5 methyl ester does not bind in the minor groove, the C7 substituent lies in the minor groove extending the rigid length of the agents, further enhancing the DNA binding-induced conformational change responsible for activation toward nucleophilic attack and catalysis of the DNA alkylation reaction.

The structural basis for in situ activation of DNA alkylation by duocarmycin SA

Duocarmycin SA is a member of a growing class of interesting lead compounds for chemotherapy, distinguished by the manner in which they bind to and react with DNA substrates. The first three-dimensional structure of a DNA adduct of an unnatural enantiomer from this family has been determined by (1)H NMR methods. Comparison to the previously determined structure of the natural enantiomer bound in the same DNA-binding site provides unique insights into the similarities and critical distinctions producing the respective alkylation products and site selectivities. The results also support the hypothesis that the duocarmycin SA alkylation reaction is catalyzed by the binding to DNA, and provide a deeper understanding of the structural basis for this unique mode of activation.

Synthesis and evaluation of 1,2,8, 8a-Tetrahydrocyclopropa[c]pyrrolo[3,2-e]indol-4(5H)-one, the parent alkylation subunit of CC-1065 and the duocarmycins: impact of the alkylation subunit substituents and its implications for DNA alkylation catalysis

The synthesis of 1,2,8,8a-tetrahydrocyclopropa[c]pyrrolo[3, 2-e]indol-4(5H)-one (CPI), the parent CC-1065 and duocarmycin SA alkylation subunit, is detailed. The parent CPI alkylation subunit lacks the C7 methyl substituent of the CC-1065 alkylation subunit and the C6 methoxycarbonyl group of duocarmycin SA, and their examination permitted the establishment of the impact of these natural product substituents. The studies revealed a CPI stability comparable to the CC-1065 alkylation subunit but which was 6x more reactive than the (+)-duocarmycin SA alkylation subunit, and it displayed the inherent reaction regioselectivity (4:1) of the natural products. The single-crystal X-ray structure of (+)-N-BOC-CPI depicts a near identical stereoelectronic alignment of the cyclopropane accounting for the identical reaction regioselectivity and a slightly diminished vinylogous amide conjugation relative to (+)-N-BOC-DSA suggesting that the stability distinctions stem in part from this difference in the vinylogous amide as well as alterations in the electronic nature of the fused pyrrole. Establishment of the DNA binding properties revealed that the CPI-based agents retain the identical DNA alkylation selectivities of the natural products. More importantly, the C6 methoxycarbonyl group of duocarmycin SA was found to increase the rate (12-13x) and efficiency (10x) of DNA alkylation despite its intrinsic lower reactivity while the CC-1065 C7 methyl group was found to slow the DNA alkylation rate (4x) and lower the alkylation efficiency (ca. 4x). The greater DNA alkylation rate and efficiency for duocarmycin SA and related analogues containing the C6 methoxycarbonyl is proposed to be derived from the extended length that the rigid C6 methoxycarbonyl provides and the resulting increase in the DNA binding-induced conformational change which serves to deconjugate the vinylogous amide and activate the alkylation subunit for nucleophilic attack. The diminished properties resulting from the CC-1065 C7 methyl group may be attributed to the steric impediment this substituent introduces to DNA minor groove binding and alkylation. Consistent with this behavior, the duocarmycin SA C6 methoxycarbonyl group increases biological potency while the CC-1065 C7 methyl group diminishes it.

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

Density functional (DF) and Møller-Plesset second order perturbation (MP2) calculations were carried out on N-acetylduocarmycin SA (N-Ac-DSA), an analogue of a series of potent antitumor antibiotics that include the duocarmycins. These computational methods were used to investigate the degree of ground state destabilization of duocarmycins that would result upon binding to DNA. Ground state destabilization has been proposed as the origin of the ligand's enhanced rate of alkylation by more than a millionfold. The conformations of the 'Unbound' and 'DNA-Bound' N-Ac-DSA were generated using available geometric data for duocarmycin SA. Specifically, the dihedral angles chi1/chi2 were locked at 6.9 degrees/4.5 degrees for the Unbound and 22.0 degrees/11.0 degrees for the Bound form. The structures were optimized using DF theory, with subsequent MP2 calculations to improve the electronic energies. All of the calculations were performed on the unprotonated (1) as well as the C6-carbonyl protonated form (2). The results showed that the ground state destabilization energies of the Unbound and Bound forms, for the unprotonated and protonated series, were fairly small (< 0.8 kcal/mol). Similarly, the difference in the electronic nature of the Unbound and Bound forms, as indicated by changes in bond lengths and charge density, were also small. In summary, it appears that twisting of two key torsional angles, the concomitant ground state destabilization, and C6-carbonyl protonation may not fully account for the significant rate increase of adenine-N3 alkylation upon binding to DNA.

Influence of 5-bromodeoxycytosine substitution on triplex DNA stability and conformation

Three triple-helical hairpin DNAs with substitution of 5-bromocytosine for cytosine in different strands have been investigated by molecular mechanics and Raman spectroscopy. The stability of the three substituted triplexes were compared with the corresponding unsubstituted triplex DNA by the molecular mechanics method. Base stacking interactions and strand--strand interactions of each triplex were analyzed in detail. Sugar conformations in these triplexes have been determined by both vibrational spectroscopy and molecular dynamics simulation. The hairpin triplexes with substitution occurring in strand I or both in strands I and III have the main sugar conformation of C3'-endo, while the triplex with substitution occurring in strand III is the combination of C3'-endo and C2'-endo sugar conformation. Theoretical results are basically in agreement with experiments.