An ONIOM study of the Bergman reaction: a computationally efficient and accurate method for modeling the enediyne anticancer antibiotics

Angle distortion model for predicting enediyne activation towards Bergman cyclization: an alternate to the distance theory

The kinetics of Bergman cyclization (BC) of enediynes into 1,4-benzene diradicals (also known as p-benzynes) have attracted interest ever since the discovery of natural enediynes which pointed out a surprising reactivity profile difference across enediynes with varying structural architectures. From the analysis of experimental kinetic data, several models were proposed to have a structure-kinetics correlation, out of which, the cd-distance model and the transition state model are the most accepted ones. Recently, Houk et al. introduced a distortion model to explain the regioselectivity of nucleophilic addition to unsymmetrical o-benzynes based on the geometry of the transition state. In the case of BC, since the reaction is endothermic, the transition state geometrically resembles the product structure which implies that in the reaction pathway, the sp-carbons of enediynes are transformed into the trigonal sp² carbons of the benzenoid product. Thus, greater bending of the interior angles at the proximal alkyne carbons in the enediynes will lead to a lower activation barrier for the BC and hence faster cyclization. This hypothesis has been tested on a series of enediynes including natural product surrogates and the extent of deviation correlates well with the kinetic results. A cut-off value for the average internal proximal angles has been proposed to categorize enediynes as per their reactivity under ambient conditions. We believe that this distortion theory offers an alternative model in designing new unnatural enediynes with desired kinetic stabilities.

Structurally similar enediynes showed decrease in interior proximal bond angles (P_a and P_b) with increase in reactivity. Enediynes with average interior proximal angles [(P_a + P_b)/2] less than 166° undergo cyclization under ambient conditions.

Efficient and accurate characterization of the Bergman cyclization for several enediynes including an expanded substructure of esperamicin A1

Incorporation of enediynes into anticancer drugs remains an intriguing yet elusive strategy for the design of therapeutically active agents. Density functional theory was used to locate reactants, products, and transition states along the Bergman cyclization pathways connecting enediynes to reactive para-biradicals. Sum method correction to low-level calculations confirmed B3LYP/6-31G(d,p) as the method of choice in investigating enediynes. Herein described as MI:Sum, calculated reaction enthalpies differed from experiment by an average of 2.1 kcal x mol(-1) (mean unsigned error). A combination of strain energy released across the reaction coordinate and the critical intramolecular distance between reacting diynes explains reactivity differences. Where experimental and calculated barrier heights are in disagreement, higher level multireference treatment of the enediynes confirms lower level estimates. Previous work concerning the chemically reactive fragment of esperamcin, MTC, is expanded to our model system MTC2.

The Reactivity of Calicheamicin γ1I in the Minor Groove of DNA: The Decisive Role of the Environment

Triggering and Bergman cyclization of calicheamicin gamma(1) (I) outside and inside the minor groove of the duplex 9mer-B-DNA sequence d(CACTCCTGG).d(CCAGGAGTG) were investigated by using density functional theory and molecular mechanics (DFT and MM) descriptions in which the ligand is completely described at the DFT and the receptor at the MM level. The calculated docking energy of calicheamicin gamma(1) (I) (-12.5 kcal mol(-1)) is close to the measured value of -9.7 kcal mol(-1) and the site specificity is in line with experimental observations. Calicheamicin is triggered in the minor groove in such a way that out of a cyclohexenone by Michael addition an E rather than a Z form of a cyclohexanone is formed, which in turn adopts a chair rather than a twistboat form. Decisive for the stereochemistry of the Michael addition is the orientation of the carbamate substituent at the headgroup of calicheamicin. Triggered calicheamicin can undergo the Bergman cyclization at body temperature only if present in its E chair form (activation enthalpy 16.4 kcal mol(-1)). An intermediate biradical is formed (docking energy -13.6 kcal mol(-1)), which has a sufficient lifetime to abstract two hydrogen atoms. Hydrogen abstraction is a two- rather than one-step process and involves the C5(H5') atom first and then the T22(H4') atom in line with experimental observations. The decisive role of using a DFT rather than an MM description for the ligand is documented.

Tuning rate of the bergman cyclization of benzannelated enediynes with ortho substituents

[reaction: see text] The Bergman cyclization of benzannelated enediynes is highly sensitive to ortho substitution. This finding opens possibilities for the rational design of conformer-specific and pH-dependent DNA-cleaving agents.