Role of Heavy Atom Tunneling in Myers–Saito Cyclization of Cyclic

Jan 19, 2016 - (43-46) An understanding of the role of heavy atom tunneling in Myers–Saito cyclization for 1a and 3a would therefore shed light on N...
0 downloads 0 Views 419KB Size
Article pubs.acs.org/JPCB

Role of Heavy Atom Tunneling in Myers−Saito Cyclization of Cyclic Enyne-Cumulene Systems Sharmistha Karmakar and Ayan Datta* Department of Spectroscopy, Indian Association for the Cultivation of Science, 2A and 2B Raja S. C. Mullick Road, Jadavpur − 700032, Kolkata, West Bengal, India S Supporting Information *

ABSTRACT: Direct dynamics calculation using canonical variational transtition state theory (CVT) inclusive of small curvature tunneling (SCT) reveals heavy atom tunneling in Myers− Saito cyclization of 10- and 9-membered cyclic enyne−cumulene systems like 1,6didehydro[10]annulene and derivative of neocarzinostatin, respectively. The pure density functional theory functional, BLYP at a 6-31+G (d,p) basis set reproduce the observed reaction energies and barriers within 1.0 kcal/mol. The calculated rate constants of cyclization inclusive of heavy atom tunneling (kCVT+SCT = 3.26 × 10−4 s−1 at 222 K; t1/2 = 35 min) are in excellent agreement with experiments (t1/2 ∼ 21−31 min). Both primary and secondary kinetic isotope effect (KIE) become enhanced significantly upon inclusion of quantum mechanical tunneling. An Arrhenius plot of KIE shows measurable curvature at the experimental temperature of 222 K. The translation vector for the cyclization reactions in the transition-states (TS) show significant motion of primary and secondary carbon atoms explaining the origin of large KIE.



INTRODUCTION Quantum mechanical tunneling (QMT) is being realized as a key factor in many chemical1−7 and biological8−10 transformations. Many H atoms mediated reactions are expected to involve QMT effects;11−15 examples of reactions governing by heavy atom tunneling are comparatively rare. The automerisation of cyclobutadiene is known to involve carbon atom tunneling.16−18 Other important reactions where carbon atom tunneling has been documented are ring closure of cyclopentane-1,3-diyl and cyclobutane-1,3-diyl,19−21 rearrangement of cyclopropylcarbenes,22 ring opening of the cyclopropylcarbinyl radical,23,24 rearrangement of semibullvalene,25 Bergman cyclization of a 10-membered-ring enediyne,26 and ring expansion of noradamantyl carbenes and their analogs.27 A salient feature of the mentioned list of reactions is a rather low reaction barrier accompanied by a small displacement of carbon atoms from reactant to product. Recently, Doubleday and co-workers have predicted strong signatures of heavy atom tunneling in the Bergman cyclization of a cyclic 10-membered ene−diyne moiety, (3Z)-cyclodec-3en-1,5-diyne.26 Many anticancer antibiotics such as esperamicin,28 calicheamicin,29,30 dynemicin31 and other related compounds possess a common structural motif, namely, an ene−diyne fragment within a ring that converts to an aromatic diradical via Bergman cyclization.32 These benzenoid type diradicals are known to abstract hydrogen atoms from sugar phosphate backbone of DNA leading to DNA strand cleavage.33,34 Another similar class of antitumor antibiotics neocarzinostatin (NCS)35 contains a diene-diyne moiety as an active chromophore.36 NCS chromophore is known to rearrange by thiol activation to a 9-membered cyclic (Z)enyne−cumulene 5a (Scheme 1) that then follows the © XXXX American Chemical Society

Scheme 1. Myers−Saito Cyclization for 1a and 3a

cyclization pathway.37 The 9- and 10-membered cyclic enyne−cumulene38−41 systems, derivative of NCS(3a) and 1,6-didehydro[10]annulene (1a) respectively, follow the similar reactivity pattern like enyne−allene42 and undergo Myers− Saito (MS) cyclization to generate highly reactive biradical accounting for its cytotoxics behavior.43−46 An understanding of the role of heavy atom tunneling in Myers−Saito cyclization for 1a and 3a would therefore shed light on NCS-based antitumor agents. Received: December 21, 2015 Revised: January 15, 2016

A

DOI: 10.1021/acs.jpcb.5b12465 J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B





COMPUTATIONAL DETAILS It is computationally challenging to study MS reaction because a closed shell reactant transforms to a trans-biradical product along the reaction coordinate. Biradical systems with small S-T gap can lead to significant configurational mixing. Hence, multireference techniques with dynamic correlation, such as CASPT2 (CASSCF method with inclusion of second order perturbation) and multireference configuration interaction (MRCI) methods should be ideal to handle them.47 However, these methods are known to suffer from systematic errors with respect to the active space dimension.48 Other ab initio method involving coupled cluster methods, such as CCSD(T)49 are computationally intractable for direct dynamics calculations. Though the hybrid density functional theory (DFT) functionals are known to overestimate the barrier height,50,51 the pure DFT functional BLYP gives excellent agreement with experiment.52 Schreiner and co-workers have studied Myers−Saito and Schmittel cyclization for a series of parent enyen−allene sytems and described BLYP level of theory as a reasonable choice for these purpose.42 We also calibrated various functionals including empirically corrected dispersion methods and BLYP indeed appears to be most suitable to describe the experimental energy barriers (see Table 1) for the MS reactions.

RESULT AND DISCUSSION 1. Suitability of BLYP Functional. To get accurate tunneling correction and rate constants it is necessary to describe the molecular potential energy surface properly. Hence, it is important to verify the performance of pure BLYP DFT functional with other dispersion corrected functional. To check the performance of pure DFT functional BLYP in describing the molecular potential energy surface, we have performed additional single point calculation on BLYP geometries at B97D3 and BLYPD3 level using 6-31+G(d,p) basis set. Table 1 lists the barrier height at different levels for the reaction 1a → 2a and 3a → 4a. The BLYP functional results show close agreement with experimental results; hence, all further calculations, namely, the tunneling calculations, were performed at BLYP/6-31+G(d,p) level. 2. Myers−Saito versus Schmittel Cyclization. Scheme 1 represents the MS cyclization for 10- and 9-membered enyne− cumulene ring systems; 1a and 3a.They undergo C3−C8 cyclization generating trans (σ−σ) biradical that eventually abstracts hydrogen atom from sugar phosphate backbone of DNA. Another alternative mode of cyclization, namely, Schmittel cyclization,65 was computationally proposed by Engels et al. for some model enyne−cumulene systems where five-membered ring cyclization occurs between C3−C7 centers leading to syn biradical product.66 Although Schmittel cyclization is very well-known for enyne−allene system, there is no experimental evidence for this C3−C7 cyclization in case of enyne−cumulene systems. Any attempt to locate the corresponding TS at BLYP level of theory reverts back to the initial reactant. Hence, we proceed with MS cyclization as the only mode of action for these NCS-based chromophores. Table 2 represents the reaction energy and activation energy at 298 K for 1a and 3a. For both the cases, the reactions are

Table 1. Enthalpy of Activation and Free Energy of Activation (at 298.15 K) in kcal/mol Calculated with Different Functional Using 6-31+G(d,p) Basis Set for 1a and 3aa density functional BLYP B97D3//BLYP BLYPD3//BLYP

1a 3a 1a 3a 1a 3a

→ → → → → →

2a 4a 2a 4a 2a 4a

ΔH‡

ΔG‡ (expt)

16.1 16.8 14.2 15.9 17.0 17.5

17.3(16.3 ± 0.1) 17.9(18.0 ± 0.1) 15.3 16.9 18.1 18.6

Article

Table 2. Enthalpy of Activation, Free Energy of Activation, and the Associated Reaction Energies (at 298.15 K) in kcal/ mol Calculated at BLYP/6-31+G (d,p) Level for 1a and 3aa

a

Corresponding experimental values are also mentioned in parentheses.

1a → 2a 3a → 4a

Therefore, pure DFT functional BLYP53−55 along with 631+G(d,p)56 basis set was employed for the estimation of reaction energies, barrier height, and direct dynamics calculations for MS reaction of the studied systems. The classical rate constants were computed using canonical variational transition state theory (CVT).57 All the calculations were performed for the singlet energy surface. Quantum effects on the reaction dynamics were incorporated, using the small curvature tunneling (SCT) approximation.58,59 All electronic structure calculations were performed with Gaussian09.60 Direct dynamics calculations were carried out with GAUSSRATE61 as the interface between Gaussian 09 and POLYRATE.62 Intrinsic reaction coordinate (IRC) calculations were performed to verify that the transition states are connected to the corresponding reactant and product. Harmonic approximation was used to determine vibrational levels of the reactant state for quantized reactant state tunneling (QRST) calculations.63 The reorientation of the dividing surface (RODS) algorithm has been applied to get accurate VaG (vibrational ground-state adiabatic potential energy) surfaces.64

ΔH

ΔG

ΔH‡

ΔG‡ (expt)

9.1 8.0

10.9 9.3

16.1 16.8

17.3(16.3 ± 0.1) 17.9(18.0 ± 0.1)

a

Corresponding experimental values are also mentioned in parentheses. The experimental temperatures for 1a and 3a are 222 and 235 K, respectively.

found to be endergonic by 9−10 kcal/mol. The free energy of activation for the compounds 1a and 3a matches well with experiment (within ∼1 kcal/mol).41 For endergonic reaction, a particle can tunnel only if it has vibrational adiabatic ground-state energy equal or greater than product valley. Hence, tunneling correction at low temperature for endergonic reaction can be erroneous. Hence, we have calculated the rate with and without tunneling at the exergonic direction, that is, 2a → 1a, then scaled them by a factor of exp(−ΔVad/kBT)67 where ΔVad represents the difference in ground-state vibrational adiabatic energy between product and reactant side, kB is the Boltzmann constant, and T symbolizes temperature. 3. Role of QMT in MS Cyclization for 9- and 10Membered Enyne−Cumulene Systems. Table 3 lists the kCVT, kCVT+SCT rate constants, and tunneling transmission coefficient (κSCT) for 1a and 3a at different temperatures (T = 150, 200, and 298 K). The rate of cyclization for 3a with and B

DOI: 10.1021/acs.jpcb.5b12465 J. Phys. Chem. B XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry B Table 3. CVT, CVT+SCT Rate Constants, and Tunneling Transmission Coefficient (κSCT) Calculated at BLYP/631+G (d,p) Level for 1a and 3a at 150, 200, and 298 Ka T(K) 1a→2a

3a→4a

a

150 200 298 150 200 298

CVT 2.34 3.59 1.97 2.39 2.57 1.69

× 10−12 × 10−6 × 10−13 × 10−7

CVT+SCT 4.17 4.93 2.52 4.43 3.72 1.95

× 10−12 × 10−6 × 10−13 × 10−7

κSCT 1.8 1.4 1.2 1.9 1.4 1.2

The rate constants, k, are reported in s−1.

without inclusion of tunneling is lower than 1a because of large activation energy. At low temperature (∼100 K), the contribution from heavy atom tunneling is significant, resulting a 5−10 times rate enhancement in both the cases, which gradually decreases with an increase in temperature. Despite the higher barrier, κSCT for 3a is greater than 1a at low temperature and they eventually become similar for T ≥ 200 K. A comparison of barrier height and width for 1a → 2a and 3a → 4a is worthwhile. The free energy of activation for 3a lies 0.6 kcal/mol above than that of 1a. The average change in C3−C8 distance on traversing from reactant to TS is 1.26 Å for 1a and 1.14 Å for 3a. Hence the barrier width for 3a → 4a is narrower by ∼0.12 Å. The variation of tunneling probability (T) with the barrier height (V0), width (w) and the mass of moving particle (m) can be approximately written as T ∼ e−w V0m .1 Hence, an alteration of barrier width affects the tunneling probability more in comparison to a similar change in barrier height. As a result, the contribution from QMT is more in 3a for which the barrier is larger yet narrower. With the increase in temperature, a larger contribution of thermal activation makes the κSCT factors similar for both of the systems. Inclusion of tunneling enhances the classical rate by a factor 1.32 for 1a, giving kCVT+SCT = 3.26 × 10−4 s−1 at 222 K. At reaction temperature, the activation energy (Ea) and the preexponential factor (A) are Ea= 17.1 kcal/mol and log A = 13.3 s−1 at CVT level, and Ea = 16.8 kcal/mol and log A= 13.0 s−1 at CVT+SCT level. QMT increases the rate of cyclization by lowering the effective barrier height. The representative tunneling energy (RTE) lies ∼0.14 kcal/mol below the maxima of the VaG surface at 222 K. This is indicative of QMT from thermally activated states rather than the lowest vibrational level (v = 0).68 The half live (t1/2) of transformation was calculated to be 45 and 35 min at CVT and CVT+SCT level, respectively, at 222 K. The experimental t1/2 lies in the range 21−31 min.41 Hence, the t1/2 calculated at CVT+SCT level shows a better agreement with the experimental result, which is also evident from the Arrhenius plot of 1a. Figure 1 shows the CVT+SCT Arrhenius plots for the transformation 1a→ 2a in the temperature range 150−300 K along with the experimentally measured rate. A similar rate enhancement was also predicted for 3a resulting in kCVT = 2.56 × 10−5 s−1 and kCVT+SCT = 3.36 × 10−5 s−1 at 235 K. The calculated t1/2 for this reaction was 7.5 and 5.7 h at CVT and CVT+SCT level, respectively, while the experimental t1/2 value was found to be ∼2 h at 235 K.38−40 It is important to mention that even though the CVT+SCT t1/2 is in relatively better agreement with the experimental rate constant, the difference in their absolute magnitude probably arises due to structural modification in our computational model compared to the experimental structure (Scheme 1).

Figure 1. Arrhenius plot of the CVT+SCT rate constants for 1a from 150 to 300 K.

4. Kinetic Isotope Effect (KIE) Study on 1a. Substitution of 12C by its heavier analog 13C is expected to reduce tunneling. Therefore, a ratio of rates involving 12C and 13C, that is, k12C/ k13C (KIE) should magnify the tunneling effects of 12C with respect to 13C. For light atom mediated reactions like hydrogen atom migration, the tunneling probability is majorly governed by the movement of the light atom along the reaction coordinate. However, for reactions where the movement of the heavy atom takes part, the tunneling probability depends strongly on the motion of different parts of the reacting system participating in the transition state imaginary mode. Hence a detailed KIE study has been carried out where both the carbon atom C-2 and C-3 has been substituted by 13C. Replacement of 12 C-3, 12C-2, and 12C-2 + 12C-3 by 13C leads to primary (1°), secondary (2°), and both (1° + 2°) KIE, respectively. Table 4 lists the KIE values for 1a → 2a transformation at different temperatures. Table 4. 1°, 2°, and 1° + 2° KIE for 1a at Three Different Temperatures; T = 200, 220, and 298 K 1a→2a primary (1°)

secondary (2°)

both (1° + 2°)

T (K)

CVT

CVT+SCT

CVT

CVT+SCT

CVT

CVT+SCT

200 220 298

1.027 1.024 1.017

1.109 1.099 1.085

1.006 1.005 1.003

1.035 1.033 1.025

1.034 1.030 1.020

1.148 1.135 1.109

As can be seen from Table 4, the computed KIE at CVT +SCT level are larger than CVT KIE. Though both CVT and CVT+ SCT KIE decreases with increase in temperature, the variation in CVT KIE is small while the CVT+SCT KIE value shows appreciable decrease with temperature. For 1a → 2a, the PKIE at CVT level is 1.024, which enhances significantly upon inclusion of tunneling to 1.099 at 222 K. The difference in activation energy, Ea (1a′ → 2a′) − Ea (1a → 2a), is 24.3 cal/ mol at CVT level and 41.3 cal/mol at CVT+SCT level, respectively, at 222 K. The κSCT factor is 1.22 for 13C (1a′ → 2a′) and 1.32 for 12C (1a → 2a) at the same temperature. Hence, thermally activated tunneling by the lighter isotope is C

DOI: 10.1021/acs.jpcb.5b12465 J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B



facilitated by not only mass effect but also by a decrease in activation energy, resulting higher KIE value at CVT+SCT level. The 2° KIE is nominal at CVT level that increases modestly with inclusion of tunneling. At 222 K, the 2° KIE for 1a is 1.005 at CVT level and 1.033 at CVT+SCT level. Although the primary carbon atom (C-3) moves by the maximum extent during the bond formation, the secondary carbon atom (C-2) also involves significant displacement along the reaction coordinate. It leads to non-negligible SKIE at CVT+SCT level.69 This is verified by animating the translational vectors of TS structure for 1a → 2a, which shows contribution from C2 as well (see Supporting Information). Substitution of both primary and secondary carbon atom by 13C gives the CVT +SCT (1° + 2°) KIE = 1.134, which is approximately 9% greater than its classical value, T = 222 K. An Arrhenius plot of ln KIE versus 1/T for 1a (Figure 2) shows appreciable curvature in the CVT+SCT curve that can clearly be seen for T = 160−200 K.

Article

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcb.5b12465. Optimized geometries, energies, thermal corrections, and harmonic frequencies for 1a−4a and the transition structure connecting them, CVT and CVT+SCT rate constants from 100 to 300 K, image of imaginary mode translational vector of TS connecting 1a and 2a, plot of VaG and VMEP surfaces for 1a and 3a and the complete lists of authors for refs 60 and 62. (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS S.K. thanks CSIR India for SRF. A.D. thanks DST, INSA, and BRNS for partial funding. We thank Professor A. Fernandez Ramos for discussion on calculation of CVT+SCT rates for endothermic reactions.



REFERENCES

(1) Bell, R. P. Tunnel Effect in Chemistry; Chapman and Hall: London, 1980. (2) Kästner, J. Theory and Simulation of Atom Tunneling in Chemical Reactions. WIREs Comput. Mol. Sci. 2014, 4, 158−168. (3) Ley, D.; Gerbig, D.; Wagner, J. P.; Reisenauer, H. P.; Schreiner, P. R. Cyclopropylhydroxycarbene. J. Am. Chem. Soc. 2011, 133, 13614− 13621. (4) Ley, D.; Gerbig, D.; Schreiner, P. R. Tunnelling Control of Chemical Reactions − The Organic Chemist’s Perspective. Org. Biomol. Chem. 2012, 10, 3781−3790. (5) Patureau, F. W. Atom Tunneling in Organic Transformations. Angew. Chem., Int. Ed. 2012, 51, 4784−4786. (6) Pan, Z.; Horner, J. H.; Newcomb, M. Tunneling in C-H Oxidation Reactions by an Oxoiron(IV) Porphyrin Radical Cation: Direct Measurements of Very Large H/D Kinetic Isotope Effects. J. Am. Chem. Soc. 2008, 130, 7776−7777. (7) Karmakar, S.; Datta, A. Tunneling Assists the 1, 2-Hydrogen Shift in N-Heterocyclic Carbenes. Angew. Chem., Int. Ed. 2014, 53, 9587− 9591. (8) Quantum Tunneling in EnzymeCatalysed Reactions; Allemann, R. K., Scrutton, N. S., Eds.; RSC Publishers: Cambridge, U.K., 2009. (9) Dybala-Defratyka, A.; Paneth, P.; Banerjee, R.; Truhlar, D. G. Coupling of Hydrogenic Tunneling to Active-Site Motion in the Hydrogen Radical Transfer Catalyzed by a Coenzyme B12-dependent Mutase. Proc. Natl. Acad. Sci. U. S. A. 2007, 104, 10774−10779. (10) Pu, J.; Gao, J.; Truhlar, D. G. Multidimensional Tunneling, Recrossing, and the Transmission Coefficient for Enzymatic Reactions. Chem. Rev. 2006, 106, 3140−3169. (11) Hynes, J. T.; Klinman, J. P.; Limbach, H.-H.; Schowen, R. L. Hydrogen Transfer Reactions; Wiley-VCH: Weinheim, Germany, 2007; Vols. 1−4. (12) Garrett, B. C.; Truhlar, D. G.; Bowman, J. M.; Wagner, A. F.; Robie, D.; Arepalli, S.; Presser, N.; Gordon, R. J. Ab Initio Predictions and Experimental Confirmation of Large Tunneling Contributions to Rate Constants and Kinetic Isotope Effects for Hydrogen Atom Transfer Reactions. J. Am. Chem. Soc. 1986, 108, 3515−3516. (13) Liu, Y. − P.; Lynch, G. C.; Truong, T. N.; Lu, D. − h.; Truhlar, D. G.; Garrett, B. C. Molecular Modeling of the Kinetic Isotope Effect for the [ 1,5] Sigmatropic Rearrangement of cis- 1,3-Pentadiene. J. Am. Chem. Soc. 1993, 115, 2408−2415.

Figure 2. Arrhenius plots of CVT and CVT+SCT KIE for 1a from 100 to 300 K.



SUMMARY In summary, calculations inclusive of small curvature tunneling reveal the role of QMT in the Myers−Saito cyclization for cyclic enyne−cumulene systems. Compared to 1a the effect of QMT is pronounced for 3a at low temperature for which the reaction passes through a narrower yet larger barrier. For both 1a and 3a, the calculated t1/2 at CVT+SCT level shows better agreement with experiments. Incorporation of tunneling enhances the KIE significantly and is shown to magnify the importance of QMT. QMT also leads to a significant curvature in Arrhenius plot of 12C/13C KIE at temperatures much higher than the cryogenic limit. Clearly, examples of heavy atom tunneling need not always be to confined to cryogenic temperatures and for small molecules. The present case shows that tunneling indeed can be relevant as well as observable at ambient conditions for important chemical transformations beyond just academic interests. D

DOI: 10.1021/acs.jpcb.5b12465 J. Phys. Chem. B XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry B

Dependent or Phototriggering Devices. Chem. Rev. 2007, 107, 2861− 2890. (34) Smith, A. L.; Nicolaou, K. C. The Enediyne Antibiotics. J. Med. Chem. 1996, 39, 2103−2117. (35) Ishida, N.; Miyazaki, K.; Kumagai, K. M.; Ribikimaru, M. Neocarzinostatin, an Antitumor Antibiotic of High Molecular Weight. Isolation, Physiochemical Properties and Biological Activities. J. Antibiot. 1965, 18, 68−76. (36) Grissom, J. W.; Gunawardena, G. U.; Klingberg, D.; Huang, D. The Chemistry of Enediynes, Enyne Allenes and Related Compounds. Tetrahedron 1996, 52, 6453−6518. (37) Suffert, J.; Abraham, E.; Raeppel, S.; Brückner, R. Synthesis of 5Ring/10-Ring Analogs of the Dienediyne Core of Neocarzinostatin Chromophore Through Palladium(0)-mediated Ring-closure Reactions. Liebigs Ann. 1996, 1996, 447−456. (38) Myers, A. G. Proposed Structure of the Neocarzinostatin Chromophore-Methyl Thioglycolate Adduct; a Mechanism for the Nucleophilic Activation of Neocarzinostatin. Tetrahedron Lett. 1987, 28, 4493−4496. (39) Myers, A. G.; Proteau, P. J.; Handel, T. M. Stereochemical Assignment of Neocarzinostatin Chromophore. Structures of Neocarzinostatin Chromophore-Methyl Thioglycolate Adducts. J. Am. Chem. Soc. 1988, 110, 7212−7214. (40) Myers, A. G.; Proteau, P. J. Evidence for Spontaneous, LowTemperature Biradical Formation from a Highly Reactive Neocarzinostatin Chromophore-Thiol Conjugate. J. Am. Chem. Soc. 1989, 111, 1146−1147. (41) Myers, A. G.; Finney, N. S. Synthesis of 1, 6-Didehydro[10]annulene. Observation of its Exceptionally Facile Rearrangement to form the Biradical 1, 5-Dehydronaphthalene. J. Am. Chem. Soc. 1992, 114, 10986−10987. (42) Schreiner, P. R.; Prall, M. Myers−Saito versus C2-C6 (“Schmittel”) Cyclizations of Parent and Monocyclic Enyne-Allenes: Challenges to Chemistry and Computation. J. Am. Chem. Soc. 1999, 121, 8615−8627. (43) Myers, A. G.; Kuo, E. Y.; Finney, N. S. Thermal Generation of α,3-Dehydrotoluene from (Z)-1,2,4-Heptatrien-6-yne. J. Am. Chem. Soc. 1989, 111, 8057−8059. (44) Myers, A. G.; Dragovich, P. S.; Kuo, E. Y. Studies on the Thermal Generation and Reactivity of a Class of (σ, π)-1,4-Biradicals. J. Am. Chem. Soc. 1992, 114, 9369−9386. (45) Saito, K.; Watanabe, T.; Takahashi, K. [4 + 2]-Type Cycloadditions of Tropone and Heptafulvene Derivatives with a Tautomeric Mixture of Cycloheptatrienylidene and Cycloheptatetraene. Chem. Lett. 1989, 18, 2099−2102. (46) Nagata, R.; Yamanaka, H.; Murahashi, E.; Saito, I. DNA Cleavage by Acyclic Eneyne-Allene Systems Related to Neocarzinostatin and Esperamicin-Calicheamicin. Tetrahedron Lett. 1990, 31, 2907−2910. (47) Andersson, K.; Malmqvist, P.-Ȧ .; Roos, B. O.; Sadlej, A. J.; Wolinski, K. Second-Order Perturbation Theory with a CASSCF Reference Function. J. Phys. Chem. 1990, 94, 5483−5488. (48) Andersson, K.; Roos, B. O. Multiconfigurational Second-Order Perturbation Theory: A Test of Geometries and Binding Energies. Int. J. Quantum Chem. 1993, 45, 591−607. (49) Scuseria, G. E. The Open-Shell Restricted Hartree-Fock Singles and Doubles Coupled-Cluster Method Including Triple Excitations CCSD(T): Application to C3+. Chem. Phys. Lett. 1991, 176, 27−35. (50) Chen, W.-C.; Zou, J.-W.; Yu, C.-H. Density Functional Study of the Ring Effect on the Myers−Saito Cyclization and a Comparison with the Bergman Cyclization. J. Org. Chem. 2003, 68, 3663−3672. (51) Wenthold, P.G.; Lipton, M. A. A Density Functional Molecular Orbital Study of the C2-C7 and C2-C6 Cyclization Pathways of 1,2,4Heptatrien-6-ynes. The Role of Benzannulation. J. Am. Chem. Soc. 2000, 122, 9265−9270. (52) Schreiner, P. R. Monocyclic Enediynes: Relationships between Ring Sizes, Alkyne Carbon Distances, Cyclization Barriers, and Hydrogen Abstraction Reactions. Singlet-Triplet Separations of

(14) Shelton, G. R.; Hrovat, D. A.; Borden, W. T. Tunneling in the 1,5-Hydrogen Shift Reactions of 1,3-Cyclopentadiene and 5-Methyl1,3-Cyclopentadiene. J. Am. Chem. Soc. 2007, 129, 164−168. (15) Karmakar, S.; Datta, A. Role of Quantum Mechanical Tunneling on the γ-Effect of Silicon on Carbenes in 3-Trimethylsilylcyclobutylidene. J. Phys. Chem. B 2014, 118, 2553−2558. (16) Whitman, D. W.; Carpenter, B. K. Experimental Evidence for Nonsquare Cyclobutadiene as a Chemically Significant Intermediate in Solution. J. Am. Chem. Soc. 1980, 102, 4272−4274. (17) Carpenter, B. K. Heavy-Atom Tunneling as the Dominant Pathway in a Solution-Phase Reaction? Bond Shift in Antiaromatic Annulenes. J. Am. Chem. Soc. 1983, 105, 1700−1701. (18) Whitman, D. W.; Carpenter, B. K. Limits on the Activation Parameters for Automerization of Cyclobutadiene-1,2-d2. J. Am. Chem. Soc. 1982, 104, 6473−6474. (19) Buchwalter, S. L.; Closs, G. L. An Electron Spin Resonance Study of Matrix Isolated 1,3-Cyclopentadiyl, a Localized 1,3-Carbon Biradical. J. Am. Chem. Soc. 1975, 97, 3857−3858. (20) Buchwalter, S. L.; Closs, G. L. Electron Spin Resonance and CIDNP Studies on 1,3-Cyclopentadiyls. A Localized 1,3 Carbon Biradical System with a Triplet Ground State Tunneling in CarbonCarbon Bond Formation. J. Am. Chem. Soc. 1979, 101, 4688−4694. (21) Sponsler, M. B.; Jain, R.; Coms, F. D.; Dougherty, D. A. MatrixIsolation Decay Kinetics of Triplet Cyclobutanediyls. Observation of Both Arrhenius Behavior and Heavy-Atom Tunneling in C-C BondForming Reactions. J. Am. Chem. Soc. 1989, 111, 2240−2252. (22) Gerbig, D.; Ley, D.; Schreiner, P. R. Light- and Heavy-Atom Tunneling in Rearrangement Reactions of Cyclopropylcarbenes. Org. Lett. 2011, 13, 3526−3529. (23) Datta, A.; Hrovat, D. A.; Borden, W. T. Calculations Predict Rapid Tunneling by Carbon from the Vibrational Ground State in the Ring Opening of Cyclopropylcarbinyl Radical at Cryogenic Temperatures. J. Am. Chem. Soc. 2008, 130, 6684−6685. (24) Gonzalez-James, O. M.; Zhang, X.; Datta, A.; Hrovat, D. A.; Borden, W. T.; Singleton, D. A. Experimental Evidence for HeavyAtom Tunneling in the Ring-Opening of Cyclopropylcarbinyl Radical from Intramolecular 12C/13C Kinetic Isotope Effects. J. Am. Chem. Soc. 2010, 132, 12548−12549. (25) Zhang, X.; Hrovat, D. A.; Borden, W. T. Calculations Predict That Carbon Tunneling Allows the Degenerate Cope Rearrangement of Semibullvalene to Occur Rapidly at Cryogenic Temperatures. Org. Lett. 2010, 12, 2798−2801. (26) Greer, E. M.; Cosgriff, C. V.; Doubleday, C. Computational Evidence for Heavy-Atom Tunneling in the Bergman Cyclization of a 10-Membered-Ring Enediyne. J. Am. Chem. Soc. 2013, 135, 10194− 10197. (27) Kozuch, S. The Reactivity Game: Theoretical Predictions for Heavy Atom Tunneling in Adamantyl and Related Carbenes. Phys. Chem. Chem. Phys. 2014, 16, 7718−7727. (28) Konishi, M.; Ohkuma, H.; Saitoh, K.-I.; Kawaguchi, H.; Golik, J.; Dubay, G.; Groenewold, G.; Krishnan, B.; Doyle, T. W.; Esperamicins, A. Novel Class of Potent Antitumor Antibiotics I. Physico-Chemical Data and Partial Structure. J. Antibiot. 1985, 38, 1605−1609. (29) Lee, M. D.; Dunne, T. S.; Chang, C. C.; Ellestad, G. A.; Siegel, M. M.; Morton, G. O.; McGahren, W. J. Borders, D. B. Calichemicins, a Novel Family of Antitumor Antibiotics. 2. Chemistry and Structure of Calichemicin γlI. J. Am. Chem. Soc. 1987, 109, 3466−3468. (30) Lee, M. D.; Dunne, T. S.; Siegel, M. M.; Chang, C. C.; Morton, G. O.; Borders, D. B. Calichemicins, a Novel Family of Antitumor Antibiotics. 1. Chemistry and Partial Structure of Calichemicin γlI. J. Am. Chem. Soc. 1987, 109, 3464−3466. (31) Konishi, M.; Ohkuma, H.; Matsumoto, K.; Tsuno, T.; Kamei, H.; Miyaki, T.; Oki, T.; Kawaguchi, H.; VanDuyne, G. D.; Clardy, J.; Dynemicin, A. a Novel Antibiotic with the Anthraquinone and 1, 5Diyn-3-ene Subunit. J. Antibiot. 1989, 42, 1449−1452. (32) Bergman, R. G. Reactive 1,4-Dehydroaromatics. Acc. Chem. Res. 1973, 6, 25−31. (33) Kar, M.; Basak, A. Design, Synthesis, and Biological Activity of Unnatural Enediynes and Related Analogues Equipped with pHE

DOI: 10.1021/acs.jpcb.5b12465 J. Phys. Chem. B XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry B Methyl-Substituted p-Benzynes. J. Am. Chem. Soc. 1998, 120, 4184− 4190. (53) Lee, C.; Yang, W.; Parr, R. G. Development of the Colic-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B: Condens. Matter Mater. Phys. 1988, 37, 785−789. (54) Becke, A. D. A New Mixing of Hartree−Fock and Local Density-Functional Theories. J. Chem. Phys. 1993, 98, 1372−1377. (55) Becke, A. D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A: At., Mol., Opt. Phys. 1988, 38, 3098−3100. (56) Clark, T.; Chandrasekhar, J.; Spitznagel, G. W.; Schleyer, P. v. R. Efficient Diffuse Function-Augmented Basis Sets for Anion Calculations. III.* The 3-21+G Basis Set for First-Row Elements, Li-F. J. Comput. Chem. 1983, 4, 294−301. (57) Truhlar, D. G.; Garrett, B. C. Variational Transition State Theory. Annu. Rev. Phys. Chem. 1984, 35, 159−189. (58) Hu, W.-P.; Liu, Y.-P.; Truhlar, D. G. Variational Transition-State Theory and Semiclassical Tunnelling Calculations with Interpolated Corrections: A New Approach to Interfacing Electronic Structure Theory and Dynamics for Organic Reactions. J. Chem. Soc., Faraday Trans. 1994, 90, 1715−1725. (59) Fernandez-Ramos, A.; Ellingson, B. A.; Garrett, B. C.; Truhlar, D. G. Variational Transition State Theory with Multidimensional Tunneling; Reviews in Computational Chemistry; Lipkowitz, K. B., Cundari, T. R., Eds.; Wiley-VCH: Hoboken, NJ, 2007; Vol. 23, pp 125−232. (60) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, J. A.; Vreven, Jr. T.; Kudin, K. N.; Burant, J. C.; et al. Gaussian 09, revision A.01; Gaussian, Inc.: Wallingford, CT, 2012. (61) Zheng, J.; Zhang, S.; Corchado, J. C.; Chuang, Y. Y.; Coitino, E. L.; Ellingson, B. A.; Truhlar, D. G. GAUSSRATE version 2010; University of Minnesota: Minneapolis, MN, 2010. (62) Zheng, J.; Zhang, S.; Lynch, B. J.; Corchado, J. C.; Chuang, Y.Y.; Fast, P. L.; Hu, W.-P.; Liu, Y.-P.; Lynch, G. C.; Nguyen, K. A. et al. POLYRATE, version 2010; University of Minnesota: Minneapolis, MN, 2010. (63) Chen, Y.-L.; Hu, W.-P. Rate Constant Calculation for HArF→ Ar + HF and HKrF →Kr + HF Reactions by Dual-Level Variational Transition State Theory with Quantized Reactant State Tunneling. J. Phys. Chem. A 2004, 108, 4449−4454. (64) González-Lafont, A.; Villà, J.; Lluch, J. M.; Bertrán, J.; Steckler, R.; Truhlar, D. G. Variational Transition State Theory and Tunneling Calculations with Reorientation of the Generalized Transition States for Methyl Cation Transfer. J. Phys. Chem. A 1998, 102, 3420−3428. (65) Schmittel, M.; Vavilala, C.; Cinar, M. E. The Thermal C2−C6 (Schmittel)/Ene Cyclization of Enyne−Allenes − Crossing the Boundary between Classical and Nonstatistical Kinetics. J. Phys. Org. Chem. 2012, 25, 182−197. (66) Musch, P. W.; Engels, B. Which Structural Elements Are Relevant for the Efficacy of Neocarzinostatin? Angew. Chem., Int. Ed. 2001, 40, 3833−3836. (67) Arnaut, L. G.; Formosinho, S. J.; Barroso, M. Tunnelling in Low-Temperature Hydrogen-Atom and Proton Transfers. J. Mol. Struct. 2006, 786, 207−214. (68) Zuev, P. S.; Sheridan, R. S.; Albu, T. V.; Truhlar, D. G.; Hrovat, D. A.; Borden, W. T. Carbon Tunneling from a Single Quantum State. Science 2003, 299, 867−870. (69) Datta, A.; Hrovat, D. A.; Borden, W. T. Calculations Find That Tunneling Plays a Major Role in the Reductive Elimination of Methane from Hydridomethylbis(trimethylphosphine)platinum: How to Confirm This Computational Prediction Experimentally. J. Am. Chem. Soc. 2008, 130, 2726−2727.



NOTE ADDED AFTER ASAP PUBLICATION This paper was published ASAP on January 28, 2016. The abstract was updated. The revised paper was reposted on January 29, 2016. F

DOI: 10.1021/acs.jpcb.5b12465 J. Phys. Chem. B XXXX, XXX, XXX−XXX