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 sp2 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. This journal is © The Royal Society of Chemistry.

Introduction

The Bergman cycloaromatization[1] (Scheme 1) involves the conversion of enediynes 1 into highly reactive para-benzynes, mainly present in the benzenoid 1,4-diradical form 2. The p-benzynes are capable of abstracting two hydrogen atoms from H-donors like the sugar backbones of complementary strands of DNA[2] leading to double strand DNA cleavage followed by self-programmed cell death (apoptosis). Bergman cyclization has also been extensively applied to the development of polymeric materials with valuable thermal/electrical properties and to the synthesis of polynuclear aromatic compounds.[3] The enormous utility of this cycloaromatization in medicine and materials science has driven a series of continuous research efforts towards exploring enediyne chemistry including synthesis, reactivity, quantum chemical theory, thermodynamics and kinetics studies.[4] A literature survey reveals that Bergman cylization still remains a dynamic field for further exploration as manifested through many of the on-gonging endeavours. Luxon et. al. theoretically studied the singlet and triplet energy surfaces for the Bergman cyclization of (Z)-hex-3-ene-1,5-diyne using coupled-cluster methods.[5] The cyclization of the triplet enediyne was exothermic with a C2-symmetric transition state while the singlet enediyne showed endothermic cyclization with a C2v-symmetric transition state. The analysis of the frontier orbitals of each stationary point led to the conclusion that the triplet pathway is energetically unfavoured as compared to the singlet pathway. This was in-line with the conclusion drawn by Turro and workers regarding the non-feasibility of the triplet pathway.[6] Our group recently reported[7] the highly regioselective halide addition to the para-benzynes obtained from the Bergman cyclization of ortho substituted naphtho-fused cyclic azaenediyne to give halo naphtho tetrahydroisoquinolines in high yields. Steric hindrance to the trajectory of the nucleophile imposed by the peripheral hydrogen atom (close to the site of nucleophilic attack) was found to be the major driving force behind the observed regioselectivity. Also, the differential electrostatic potential as well as distortions at reactive centres played a minor role in controlling the regioselectivity. Binder and co-workers reported[8] the synthesis of main-chain diamino enediyne embedded polymers and showed that these polymers in solution exert a chain-length dependent DNA cleavage activity under physiological conditions, which is additionally tuneable by modulating the stereoelectronic environment via their substitution patterns. Photochemical activation of such systems led to the formation of long-lived free radicals.

Scheme 1

Bergman cyclization.

The core structure of the enediyne 1 is constituted of a double bond flanked by two acetylenic moieties for which the out-of-plane orbitals are in conjugation with the double bond. The intrinsic reactivity of enediynes, controlled in the general sense, through either strain or electronic effects, has intrigued the chemists for decades and that has driven a lot of experimental and theoretical studies towards delineating the structural parameters of enediynes that govern the stability and reactivity of enediynes towards Bergman cyclization.[9] One such seminal work (cd-distance theory)[9] was reported by Nicolaou and co-workers where it was empirically shown that the distance between carbon atoms c and d (c and d are the acetylenic carbons remote from double bond of the enediyne, also designated as the alkyne termini as shown in Fig. 1) plays a crucial role in determining the rate of p-benzyne formation.

Fig. 1

Labelling of the crucial carbon atoms and angles in the enediyne scaffold.

Upon correlating the cd-distances between the termini of the triple bonds of the enediyne systems and the rate of cycloaromatization, it was revealed that a decreased cd-distance and a consequent greater proximity between the acetylenic termini made the systems more prone towards Bergman reaction in order to relieve the torsional strain arising out of the alkyne electron cloud repulsion. A critical upper limit for the cd-distance of around 3.2–3.3 Å in the ground state, was proposed for the reaction to occur at a measurable rate at ambient temperatures.[9] These studies provide an explanation of why naturally occurring enediyne antibiotics contain either 9 or 10 membered cyclic enediyne cores as their cd-distances fall below or within the critical limit and thus they can be easily cycloaromatized under ambient conditions after removal of a so-called locking device. Subsequently, Snyder,[10] Tipsword,[11] Magnus and co-workers,[12] based on molecular mechanics calculations and experimental data pointed out the limitation of the cd-distance rule and suggested that the reactivity of an enediyne is dependent not only on the strain energy of the ground state but also on that of the transition state – the ease of cycloaromatization of an enediyne is critically determined by the difference in the strain energy between the ground state and the transition state. An enediyne with a smaller strain energy difference between the ground and transition states undergoes easier cycloaromatization.

In 2010, Houk and co-workers proposed an interesting theoretical model known as the aryne distortion model[13] based on experimental and theoretical studies which explained the regioselectivity of nucleophilic addition reactions to unsymmetrical o-benzynes. According to this model, the nucleophilic addition to unsymmetrical benzynes was found to be favoured for attack at that aryne carbon which requires the minimum geometrical and energetic change in going from the aryne to the TS structure. Thus, closer the aryne bond angle to linearity, greater was the reactivity observed for that centre. This observation may be exemplified through the results of the computational and synthetic studies on the regioselectivity of nucleophilic attack on 5,6-indolyne (Scheme 2). The aryne carbon (as shown in 5,6-indolyne) undergoing nucleophilic addition gets flattened in the corresponding transition state (5a – due to attack at C-5; 6a – due to attack at C-6). The slightly greater linearity of C-5 (interior angle = 130°) than at C-6 (interior angle = 127°) in 4, increased the susceptibility of C-5 towards the nucleophilic attack as that would require lesser geometric and energetic change upon going to the transition state 5a than C-6 nucleophilic attack leading to the transition state 6a. Flattening of the aryne carbon (C-5) also induces significant p-character of the in-plane hybrid orbital at the site of attack and a slight positive charge (attributed to electron withdrawal by indolyl nitrogen at para-position, as shown for structure 5b), while the adjacent angle is compressed (C-6), leading to increased s-character of the in-plane hybrid orbital to stabilize the developing carbanion. These promising results obtained from correlating the angular distortion with the observed regioselectivity in the aryne systems led to the inception of our current work.

Scheme 2

Aryne distortion model to explain the regioselectivity in nucleophilic addition to indolyne systems.

It is well known that in a multistep reaction, the activation energy of the rate-determining step controls the over-all rate of the process.[14] Generally, the first step of Bergman cyclization involving formation of the 1,4-benzenoid diradical is the slow rate determining step and is an endothermic process.[1] Subsequent hydrogen atom abstraction step follows a much faster kinetics. So, the transition state for the rate determining step would resemble the diradical species (late transition state) in accordance to Hammond's postulate.[15] Theoretical studies[16] on the transition state of Bergman cyclization support this postulation as well. Thus, it is a reasonable assumption that the interior C–C–C bond angles in the transition state will get reduced in order to facilitate the attainment of an of 120° as the hybridization state of the acetylenic carbons change from sp to sp[2] upon formation of the p-benzynes. Extrapolation of the distortion model to enediyne system leads to the hypothesis that greater the deviation of the alkyne interior angle from linearity and closer it is towards trigonal geometry, greater will be the extent of mimicking the transition state, consequently lesser will be the difference in energy between the ground state and the transition state (lower activation energy) and hence faster should be the rate of cycloaromatization (Scheme 3). One may however argue that decreasing the internal acetylene angles should accompany a concomitant decrease in cd-distance, however this is not always true as shown by example 14 and 15 described later. It is interesting to note the for regioselective nucleophilic addition to o-arynes, distortion model predicts that greater the angular shift from trigonal to linear geometry, higher the reactivity; whereas in case of Bergman cyclization of enediynes, greater the angular shift from linear to trigonal geometry, higher should be the reactivity.

Scheme 3

The changes in the interior angles Pa and Pb during the formation of the p-benzyne diradical via BC.

Basis of the present work

Based on the above preamble, we carried out geometry optimization for a series of enediynes using DFT-method and an attempt was made to establish a correlation between their experimentally obtained half-lives with their respective angular distortions. The hypothesis was indeed supported by the kinetic data and could be reasonably extrapolated to systems for which the kinetic studies have not been done yet. It was found that among structurally similar enediynes, increase in the interior angular distortions (less than 180° and closer to 120°) of the acetylenic carbons at positions a and b of the enediyne motif, leads to decrease in half-life indicating an increase in reactivity. Distortions in enediyne core were primarily quantified through the proximal interior angles at alkyne carbons a and b mainly due to the following reasons: (i) in Bergman cyclization, the transition state is diradical-like in nature and atoms a and b become the diradical centres in the diradical-product; (ii) these atoms are not bonded to any substituent which is not a part of the enediyne core; (iii) unlike atoms c and d, a and b do not participate in any new bond formation. To further consider the combined effect of angular distortions over both the proximal alkyne carbons (a and b), we calculated the average interior proximal angle [(Pa + Pb)/2].[1,17] Subsequently, a threshold-value was proposed for this average parameter which has been successfully used to distinguish enediynes which cycloaromatize under ambient conditions from those which do not. It was gratifying to note that for all the studied systems, this threshold value was found to be consistent in categorizing the enediynes into either of the two categories (discussed later) as per their reactivity under ambient conditions.

However, the extent of angular deviations, either for individual proximal alkyne carbon atoms (Pa and Pb) or the average parameter [(Pa + Pb)/2], should not be correlated to their individual reactivity order for the dissimilar enediynes. For such structurally different systems, the correlation may or may not hold good.

Thus, in addition to other major theories like the cd-distance model and transition-state model, this one offers a very simple yet alternative theoretical approach towards predicting the reactivity of enediynes. The results of our study are described in this paper.

Results and discussion

To study the effect of angular distortions, present in the enediynes on their reactivity, several different enediynes (10-33) were selected (Fig. 2) including two truncated structures of natural enediyne retaining the core enediyne structure along the trigging device in antibiotics calicheamicin γ1 (24) and dynemicin A (26).

Fig. 2

Enediynes geometry optimized using DFT-based calculations.

The ground state equilibrium structures of all the molecules were computed with DFT using B3LYP/6-31G+(d,p) and B3LYP/6-311G++(d,p) levels of theory, in gas phase, considering the singlet state (since we are dealing with the thermal mode of Bergman cyclization),[5,6,18] using the Gaussian 09 (Revision-A.02)[19] program package. It is to be noted that the output obtained from using B3LYP/6-31G+(d,p) optimization was used as the input geometry while using B3LYP/6-311G++(d,p).[20]

For each molecule, the basis set 6-311G++(d,p) furnished lowest energy values (included in ESI†). The following parameters were computed (Table 1 – acyclic and monocyclic enediyne; Table 2 – polycyclic enediyne) for each of the geometry optimized structures: (a) non-bonded distances between carbons a–b, c–d; (b) angles Pa, Pb, Dc and Dd; (c) energy values.

Distance between the carbon atoms c–d and angles Pa, Pb computed for the acyclic (10–12) and monocyclic (13–23) enediynes using DFT-calculation

Entry	Enediyne molecule (along with their reactivities)	c–d (distance in Å)	Bond angle (in degrees)		Entry	Enediyne molecule (along with their reactivities)	c–d (distance in Å)	Bond angle (in degrees)
			P a	P b				P a	P b
A	10	4.488a	182.26a	182.27a	H	17	3.763a	171.88a	170.56a
	Stable at 25 °C[1]	4.463b	182.06b	182.04b		Stable at 25°C[9a,21]	3.748b	171.76b	170.29b
B	11	4.171a	183.62a	183.62a	I	18	4.168a	177.45a	177.87a
	Stable at 25 °C[1]	4.152b	183.39b	183.39b		Stable at 25°C[9a,21,23]	4.149b	177.07b	177.62b
C	12	4.518a	182.41a	182.41a	J	19	3.339a	164.24a	166.38a
	Not reported, expected to cyclize only at high temperature	4.490b	182.06b	182.06b		t 1/2 = 72 h at 23 °C[24]	3.338b	164.72b	166.53b
D	13	2.918a	156.80a	156.51a	K	20	3.348a	166.64a	167.68a
	Expected to cyclize at ambient temperature[9b]	2.930b	157.55b	157.52b		t 1/2 = 52 h at 65 °C[24]	3.345b	167.02b	167.77b
E	14	3.376a	165.40a	165.40a	L	21	3.298a	163.81a	164.23a
	t 1/2 = 18 h at 25°C21,9a-b	3.383b	165.89b	165.89b		t 1/2 = 11.8 h at 37°C[9a,b,21]	3.310b	164.54b	164.97b
F	15	3.353a	166.90a	166.90a	M	22	3.441a	166.33a	166.25a
	ca t1/2 = 10 h at 84 °C[22]	3.357b	167.34b	167.34b		t 1/2 = 4 h at 50 °C[21,25]	3.436b	166.60b	166.43b
G	16	3.350a	167.10a	167.11a	N	23	3.482a	167.01a	167.01a
	Not reported, but expected to be stable because of benzannulation	3.352b	167.53b	167.54b		t 1/2 = 22 h at 50 °C[21,25]	3.470b	167.10b	167.05b

B3LYP/6-31G+(d,p) theory level.

B3LYP/6-311G++(d,p) theory level. Refer to the ESI for the corresponding values for a–b distance, angles Dc, Dd and the energy (at both the theory levels).

Distance between the carbon atoms c–d and angles Pa, Pb computed for the polycyclic enediynes (24–33) using DFT-calculation

Entry	Enediyne molecule (along with their reactivities)	c–d (distance in Å)	Bond angle (in degrees)		Entry	Enediyne molecule (along with their reactivities)	c–d (distance in Å)	Bond angle (in degrees)
			P a	P b				P a	P b
O	24	3.504a	168.15a	168.11a	T	29	3.418a	168.62a	164.93a
	Not reported, but expected to be stable at physiological (ambient) temperature	3.490b	168.14b	167.83b		t 1/2 = 2.1 h at 71 °C[12]	3.414b	168.71b	165.29b
P	25	3.171a	161.92a	162.59a	U	30	3.468a	168.99a	166.68a
	Not reported, but expected to cycloaromatize at physiological (ambient) temperature	3.182b	162.68b	163.41b		Reported to cycloaromatize 650 times slower than 29 at 124 °C[12]	3.456b	168.73b	166.87b
Q	26	3.654a	170.17a	168.98a	V	31	3.434a	167.53a	166.67a
	Not reported, but expected to be stable at physiological (ambient) temperature	3.628b	169.46b	168.69b		t 1/2 = 10.9 h at 81 °C[26]	3.429b	167.58b	166.93b
R	27	3.300a	163.21a	164.90a	W	32	3.420a	167.17a	166.49a
	Not reported, but expected to cycloaromatize at physiological (ambient) temperature	3.309b	163.87b	165.60b		Reported to cycloaromatize 1.6 times faster than 31, and at 65 °C[26]	3.416b	167.22b	166.73b
S	28	3.252a	163.36a	162.98a	X	33	3.315a	165.33a	164.73a
	Not reported, but expected to cycloaromatize at physiological (ambient) temperature	3.260b	164.23b	163.47b		Reported to cycloaromatize 500 times faster than 31 at 37 °C, t1/2 = 15.3 h at 37 °C[26]	3.317b	165.62b	165.29b

B3LYP/6-31G+(d,p) theory level.

B3LYP/6-311G++(d,p) theory level. Refer to the ESI for the corresponding values for a–b distance, angles Dc, Dd and the energy (at both the theory levels).

Observing the trend in the theoretically calculated structural parameters from the geometry optimized structures revealed that along with the cd-distance, the interior angles subtended at carbon atoms a, b, c and d (Pa, Pb, Dc and Dd respectively) also undergo a gradual change corresponding to the change in the rate of cycloaromatization. Correlating the literature reported experimental half-lives with the calculated Pa and Pb values depicted that for geometrically similar enediynes, there occurs a decrease in these angles (Pa and Pb) with increase in their reactivity (Tables 1 and 2). Consequently, the average internal angles were also found to decrease with increase in reactivity (Table 3).

Classification of enediynes based on their average interior proximal angles

Molecule	Theory level: B3LYP/6-31G+(d,p) average interior proximal angle [(Pa + Pb)/2] (in degrees)	Theory level: B3LYP/6-311G++(d,p) average interior proximal angle [(Pa + Pb)/2] (in degrees)
CATEGORY-I: average interior proximal angle is less than 166° (expected to cyclize at ambient temperatures)
13	156.66	157.54
14	165.40	165.89
19	165.31	165.63
21	164.02	164.76
25	162.26	163.05
27	164.06	164.74
28	163.17	163.85
33	165.03	165.46
CATEGORY-II: average interior proximal angle is less than 166° (expected to be stable at ambient temperatures)
10	182.27	182.05
11	183.62	183.39
12	182.41	182.06
15	166.90	167.34
16	167.11	167.54
17	171.22	171.03
18	177.66	177.35
20	167.16	167.40
22	166.29	166.52
23	167.01	167.08
24	168.13	167.99
26	169.58	169.08
29	166.78	167.00
30	167.84	167.80
31	167.10	167.26
32	166.83	166.98

The acyclic enediynes are highly unreactive as compared to their cyclic analogues. This is substantiated by the corresponding very high values of Pa and Pb (entries A–C in Table 1) for molecules 10–12. Among the simple unsubstituted cyclic enediynes, nine membered (13) is so reactive[9] that it has not been isolated (as it readily cyclizes under ambient condition) and it showed the lowest values (Fig. 3) for Pa and Pb (entry D in Table 1) among the molecules studied. A subsequent increase in the ring size in molecules 14 and 17[9] (ring sizes 10 and 11 respectively) decreases its reactivity manifested through a consequent increase in the values of Pa and Pb (entries E and H respectively, in Table 1). Experimental studies[22,24] suggested that fusion of aryl rings with the cyclic enediyne scaffold leads to a decrease in their reactivity. This trend was supported by the increased Pa and Pb values in presence of fused phenyl or naphthyl rings in the ten-membered carbocyclic enediyne 14 to form molecules 15 and 16 respectively (entries F and G in Table 1).[22] Similarly, increase in half-life was found to be accompanied by the increase in the values of Pa and Pb upon going from ten-membered aza-enediyne 19 to it's for the benzo-fused derivative 20 (entries J and K respectively, in Table 1).[24] It is interesting to note that the computed cd-distances decreased with subsequent fusion of aromatic ring(s) in molecules 15 and 16 (as compared to simple ten-membered cyclic enediyne 14) which contradicted the order of their reactivities as per Nicolaou's distance theory.

Fig. 3

Geometry optimized structures of a few representative enediynes (10,11, 13, 14, 15, 17, 19 and 20) at B3LYP/6-311G++(d,p) level.

Among the studied diols[9] (entries L–N, in Table 1), molecule 21 was found to be most reactive while 23 was found to be least reactive (as reflected from their experimentally reported half-lives). The proposed distortion theory as manifested through the proximal interior angles (Pa and Pb) supported their reactivity order as the angular distortion was found to be highest for 21 and lowest for 23.

The idea of correlating reactivity with extent of bond angle distortion was further validated on the simplified structural analogues of two well established enediyne antibiotics – calicheamicin γ1 and dynemicin A.[2] The enediyne functionality of both these molecules are locked within rigid bridged ring systems awaiting activation to undergo the Bergman cyclization. Calicheamicin γ1 has a bicyclic aglycon enediyne unit that is triggered for the cycloaromatization through a conjugate addition of thiolate to the proximate α, β-unsaturated carbonyl system. On the other hand, dynemicin A is activated through NADPH or glutathione induced epoxide ring opening reaction. The newly induced strain makes the triggered form much more reactive (activated) than the untriggered form. We selected the surrogates (truncated and simplified structural analogues due to computational constraints) of calicheamicin γ1 (represented by molecule 24, entry O in Table 2) and dynemicin A (represented by 26, entry Q in Table 2) as our model systems. The activated form of calicheamicin γ1 (conjugate addition product) was represented in its simplified form by molecule 25. For the dynemicin A analogue, two simplified analogues of its activated forms were considered: (a) epoxide ring opened by sulphur based nucleophilic attack (represented by molecule 27); (b) epoxide ring opened by hydride attack (represented by molecule 28). For both calicheamicin γ1 (24) and dynemicin A analogues (26), as expected, the values of angles Pa and Pb were found to be closer to trigonal geometry (Table 2) upon their activation (Fig. 4).

Fig. 4

Geometry optimized structures of natural product analogues 24, 25, 26 and 28 at B3LYP/6-311G++(d,p) level.

During the course of work, our attention was drawn to the interesting relative reactivity profiles of a few bicyclic enediynes 29 to 33 (entries T–X, in Table 2). Magnus reported[12] that although the cd-distance is lower for 12-ketobicyclo[7.2.1]enediyne (30) than 13-ketobicyclo[7.3.1] enediyne (29), yet 30 cycloaromatizes 650 times slower than 29 which has a half-life of 2.1 h at 71 °C. This observation by Magnus was in disagreement with the cd-distance theory proposed by Nicolaou. Our theoretical studies[28] revealed that the interior proximal angles Pa and Pb are slightly higher (Table 2, entries T and U) for 30 than that for 29 (Fig. 5). Also, the average interior angle (Table 3) is higher for 30 than 29. Thus, the angle distortion parameters are successfully able to rationalize the reactivity order. Interestingly, unlike Magnus' report, our optimization studies revealed a marginal increase in the cd-distance upon going from 29 to 30 (Table 2, entries T and U), which also extended support to the distance theory.[29]

Fig. 5

Geometry optimized structures of a few representative enediynes (10,11, 13, 14, 15, 17, 19 and 20) at B3LYP/6-311G++(d,p) level.

For dynemicin core azabicyclo[7.3.1]enediynes (31, 32 and 33),[26] Magnus reported that the reactivity increase (as revealed from their kinetic studies) upon going from 31 to 33 in spite of having almost same cd-distance. Our optimization studies revealed that upon going from 31 to 33, increase in reactivity was accompanied by a concomitant decrease in the interior proximal angles (Table 2, entries V to X) as well as average internal proximal angles (Table 3). Thus, the angle distortion hypothesis was able to justify the order of experimental half-lives. The cd-distance obtained at B3LYP/6-31G+(d,p) and B3LYP/6-311G++(d,p) levels of theory, were found to decrease upon going from 31 to 33 (Table 2, entries V to X). Thus, like molecules 29 and 30, the theoretically obtained cd-distances for 31, 32 and 33 are in accordance to the cd-distance theory.

Based on our studies, we envisaged that the enediynes may be broadly classified into the following two categories depending on their reactivity: (i) category I – those which undergo Bergman cyclization under ambient conditions (at temperature ≤ 37 °C); (ii) category II – those which are relatively stable under ambient conditions and undergo Bergman cyclization at higher temperatures (at temperature > 37 °C). Keeping this in perspective, the scope of our hypothesis was further extended by proposing an empirical cut-off value of 166° for the average interior proximal angle [(Pa + Pb)/2]. It was found that the enediynes having this average value less than 166° belonged to category I, and the ones having this average value greater than 166° belonged to category II (Table 3). None of the studied enediynes was found to violate this empirical threshold value. Further, this threshold value was also successful in uniformly classifying the enediynes into the two categories, in accordance to their reactivities, irrespective of their structural similarity/dissimilarity.

It is interesting to note that the benzannulated enediynes 15, 16 and 20 (entries F, G and K respectively, in Table 1) defied the modified cd-distance range of 3.4–2.9 Å that was proposed by Schreiner[9] for spontaneous cyclization, as these systems showed a cd-distance within the proposed range, yet they cyclized only under higher temperature conditions. On the other hand, the threshold value of the average interior proximal angle (166°) as proposed in this work, was able to categorize (category II in Table 3) them as stable molecules under ambient conditions.[30]

Conclusion

A novel and simple theoretical approach has been proposed to correlate the reactivity of enediynes with their angular distortions at the proximal alkyne carbon atoms. Since the p-benzyne formation during Bergman cyclization adopts a late transition state, so substrates which mimic the corresponding p-benzyne follow a faster reaction kinetics. For the structurally similar kind of enediyne systems, it was observed that greater the angular deviation of the proximal alkyne carbons from linearity, lower is their experimental half-lives. A threshold value of 166° was empirically proposed for the composite parameter of average interior proximal angle to distinguish the potentially reactive enediynes from the unreactive ones. Although activation energy calculation for the respective cyclization pathway has the potential to probably give a more accurate correlation between the structure and reactivity than the results obtained via the proposed semi-quantitative distortion model, considering the advantage of the latter in terms of computational resources and time, the current hypothesis is a viable trade-off which may be used to design various enediyne scaffolds with controlled reactivities.

Computational details

All of the computations were performed in gas phase with the Gaussian 09 (Revision-A.02)[19] software package and calculation of the structural parameters were done using Gauss View 6.0.[31] Becke's three-parameter hybrid functional,[32] the correlation functional of Lee, Yang, and Parr[33] and split-valence basis sets 6-31G+(d,p)[34] and 6-311G++(d,p)[35] were used.[36] Electronic energies in atomic unit (Hartrees) and Cartesian coordinates of the optimized geometries of enediynes are reported in the ESI.† The nature of the stationary point was characterized by the absence of imaginary frequency in the vibrational frequency analysis.

Conflicts of interest

There are no conflicts to declare.