Dynamical Effects along the Bifurcation Pathway Control

7 days ago - Post-transition-state dynamics during the deazetization of 3 resulting in two degenerate semibullvalenes (4 and 5) have been investigated...
1 downloads 18 Views 22MB Size
Subscriber access provided by Northwestern Univ. Library

Article

Dynamical Effects Along Bifurcation Pathway Control Semibullvalene Formation in Deazetization Reactions Nilangshu Mandal, and Ayan Datta J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.7b09533 • Publication Date (Web): 09 Jan 2018 Downloaded from http://pubs.acs.org on January 10, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 7 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Dynamical Effects Along Bifurcation Pathway Control Semibullvalene Formation in Deazetization Reactions Nilangshu Mandal, Ayan Datta*

Department of Spectroscopy, Indian Association for the Cultivation of Science, 2A and 2B Raja S.C. Mallick Road, Kolkata – 700032, India. Email: [email protected]

ABSTRACT:

Post-transition-state dynamics during the deazetization of 3 resulting in two degenerate semibullvalenes (4 and 5) have been investigated with density functional theory (DFT) and quasi-classical trajectory (QCT) calculations. Removal of N2 from 3 occurs through a synchronous and concerted pathway through an ambimodal transition-state (TS1). In addition to TS2, the exclusively anticipated product from minimum energy pathway (MEP) calculations, trajectories initiated from TS1 produce 4, TS2 and 5 in 1:1:1 ratio. Isotopic substitutions (12C(13C/14C)-H(D) at 1-2 positions) result in purely Newtonian kinetic isotope effects (4:5 ~1.4 for 13C1-13C2), an unequivocal evidence for dynamics controlled product formation.

Tetrazine (1) is known to act as an electron deficient diene for inverse electron demand Diels-Alder reactions with electron-rich dienophiles and has been extensively utilized in total synthesis for molecules like streptonigrin.12-15 Souer et al. have demonstrated an efficient method for the production of semibullvalenes (4 and 5) by the extrusion of N2 from tetrazine (1) and bis-cyclopropene (2).16 Such efficient and fast cycloaddition reactions are widely used in bio-orthogonal chemistry.17-19 Tetrazine cycloaddition reactions have been successfully utilized in live cell imaging with various dienophiles.21-22 Cyclopropene derivatives are one of the best dienephiles as they are highly strained as well less hydrophobic in nature. Recently, using molecular dynamic study, Houk and co-workers established the nature of Diels-Alder reactions of tetrazines with alkenes alongwith N2 extrusions from adducts.23 Singleton and co-workers reported cycloadditions of cyclopentadiene with diphenylketene and dichloroketene and found nonstatistical recrossing to be essential to describe the correct isotope effects.24 Using dynamic effect, Thomas et al. elucidated the mechanism through which product selectivity is controlled in the Diels-Alder cycloadditions of 3methoxycarbonylcyclopentadienone with 1,3-dienes.25 Recently, Houk and co-workers explored the origin of periselectivity in the ambimodal [6+4] cycloadditions of tropone to demethylfulvene using dynamical calculations.26 Tantillo and co-workers showed that carbocation-π interactions can manipulate dynamic behaviors of post-transition-state bifurcation reactions.27 Stopped flow kinetic studies of tetrazine cycloaddition toward sequential transition states have been reported by Sadasivam et al.28 Similarly, Yamataka and co-workers have established dynamic path bifurcation in the Beckmann reaction.29 Recently, Lόpez and co-workers have theoretically studied dynamical effects in [3,3] sigmatropic rearrangement along bifurcation pathways.30 Houk31 and Tantillo32 have reviewed several examples of post-TS bifurcations in organic reactions.

INTRODUCTION One of the most popular models to understand kinetic preferences of chemical reactions has been the transition-state theory (TST) or its improved variant considering statistical recrossing (the variational transition-state theory, VTST).1-2 However, several evidences have emerged recently wherein post-transition state non-statistical dynamics and recrossing become essential to qualitatively describe the outcome of organic reactions.3-8 For such cases, a static picture based on potential energy surface (PES) can give only the thermodynamic picture accurately while, product selectivity, regioselectivity are controlled by dynamics of nuclear motions and momenta of the atoms along non-internal reaction coordinate (non-IRC) pathways for bifurcation reactions.9-10 Reaction between tetrazine (1) and bis-cyclopropene (2) producing two degenerate semibullvalenes is one such example (Scheme 1).11

Motivated by such growing interest in reactions of tetrazenes, we have computationally studied the reaction between tertrazine (1) and bis-cyclopropene (2) resulting in 3 which subsequently forms degenerate semibullvalenes (4 and 5) through a bispericyclic transition state. Herein, we represent static density functional theory (DFT) calculations for the 3 → semibullvalenes (4 and 5). We have incorporated various functional groups on biscyclopropene (R = H, F, Cl, Br, OH, CN, NH2, Me and Et) to investigate the effects on their activation barriers, relative energies as well as in their valley ride inflection points (VRIs). We also perform detailed quasi-classical trajectory simulations from the

Scheme 1: Synthetic pathway for semibullvalene formation through deazetization of 3. 1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ambimodal transition state structure to examine the dynamical behavior. Newtonian kinetic isotope effects (KIEs) were computed to describe the origin of kinetic selectivity through isotope labeling.

COMPUTATIONAL METHODS All the structures are optimized at B3LYP/6-31G(d) and M062X/6-31G(d) levels of theory. For estimation of the barrier heights a higher basis set, TZVP was also utilized. All the static electronic structure calculations were performed with Gaussian 09.33 For a comprehensive understanding concerning “on-the-fly” behavior of the bifurcation pathways, the trajectories were propagated from the ambimodal transition state structure (TS1) (shown in Figure 1(a) and Figure 2(b)) at B3LYP/6-31G(d) level of theory. Normal mode sampling was performed by addition of the zero-point energy for each real normal mode of the TS followed by Boltzmann sampling to get a set of coordinates for trajectory propagation at 300K. The usage of TS distribution is necessary to get an unbiased product distribution. A velocity-Verlet algorithm was used to integrate the classical equations of motion using the Singleton’s Progdyn24 program with time step of 1 fs. The rate constants were computed using TheRate34 program. Rice– Ramsperger–Kassel–Marcus (RRKM) calculations were carried out using the Beyer–Swinehart direct count algorithm35-36. Product 4 and 5 were used as the starting points for lifetime calculations.

Page 2 of 7

Substitutions

TS1

VRI

TS2

4

5

I (R = H) II (R = F) III (R = Cl) IV (R = Br) V (R = CN) VI (R = OH) VII (R = NH2) VIII (R = Me) IX (R = Et)

12.5 11.9 11.8 11.7 11.9 11.7 11.5 12.2 12.3

6.8 5.9 6.1 6.0 2.0 6.0 6.0 6.8 6.8

-41.6 -48.3 -47.4 -47.7 -40.6 -46.4 -44.2 -42.1 -42.0

-46.4 -48.8 -48.9 -48.8 -47.3 -47.4 -47.9 -46.8 -47.4

-46.4 -52.1 -50.4 -50.7 -44.7 -49.5 -47.4 -45.6 -45.4

Table 1: Zero-point energy (ZPE) corrected relative free energies at B3LYP/6-31G(d) level of theory at 298.15K for I-IX with respect to 3. Energies are in kcal/mol. The intrinsic reaction coordinates (IRC) connect 3 with TS2 via TS1 in I (R = H). For all other reactants (II-IX), IRC connects 3 with the most stable of the semibullvalene (4 or 5) through TS1 without a detour into TS2 (shown in Figure S23–S28 and Table S3 in Supp. Info. file). Clearly, electronic perturbation by functional groups that break the high symmetry in TS1 can tune the reaction pathway. For I, a C2v symmetry is maintained all along the IRC from 3 to TS2, whereas for the substituted analogues, only a Cs symmetry is preserved along 3 to 4/5. The structures at TS1 and VRI-point and TS1 and 5 for I and III (R = Cl) are respectively shown in Figure 1 (a) – 1(d). For I, the VRIpoint is responsible for the bifurcation into 4 and 5 alongwith the formation of TS2. Specifically, for either C-N bond length of 2.13 Å, νimaginary(1)= 367.5i cm-1 results in the formation of TS2 while νimaginary(2)=149.6i cm-1 results in the formation of the two degenerate semibullvalene products namely, 4 and 5 at B3LYP/631G(d) level in I (Figure S22, Supp. Info. file).

RESULTS AND DISSCUSSION Table 1 represents the relative energy change in terms of ∆G (at 298 K) for the conversion of 3 to semibullvalenes (4 and 5) with loss of N2 for I-IX. It is clear that the activation energy for the first-step (3 → TS2 or 3 → 4 or 3 → 5) is significantly higher than that for the degenerate cope-rearrangement (4 ↔ 5 via TS2). ‡ For example, in I (R = H), ∆G1 (3 → TS2 /3 → 4/ 3 → 5) =12.5 kcal/mol (10.2 kcal/mol at B3LYP/aug-cc-pVTZ level) while ‡ ∆G2 (4 → 5) = 4.8 kcal/mol (3.8 kcal/mol at B3LYP/aug-ccpVTZ level). Calculations at a higher basis set clearly do not alter the basic picture of the potential energy surface (PES), see Supp. Info. for further details. The computed reaction energies and barriers are in excellent agreement by previous reports on I by Birney and co-workers.11 Substitution by functional groups reduces the barrier for the rate-determining step with the ‡ maximum lowering, ∆∆G1 =1.0 kcal/mol occurring for R = NH2. Due to substitution, the degeneracy of 4 and 5 is lost and the relative stability of 4 and 5 can be understood based of the electron withdrawing/donating nature of the substituents. For IIIV and VI, preference for 5 over 4 arises as a consequence of the electron-withdrawing nature of the groups that increases the scharacter of the adjacent bonds (See Supp. Info. File for NBO analyses). This results in an increment in internal angle between the bonds (within the cyclopropyl ring) that makes 4 unstable. For VII-IX presence of electron-donating groups results in reduced scharacter of the adjacent bonds and lower internal angle within the ring. This stabilizes 4. In case of V, inspite of the electron withdrawing nature of –CN, its conjugation with the cyclopropane ring makes product 4 relatively more stable.

Figure 1: (a) and (b) are optimized TS1 structures for I (R = H) and III (R = Cl). (c) Valley ridge inflection (VRI) point for I (R = H). (d) Optimized chloro substituted semibullvalene (5). These structures are optimized at B3LYP/6-31G(d) level of theory. Interestingly, for II-IX, the IRC did not encounter the VRI-point close to the TS1 region and is located further right from TS1 (along the mass-weighted reaction coordinate, s). Hence, the steepest descent path from the saddle point structure (TS1) along IRC produces the thermodynamically more stable product 2

ACS Paragon Plus Environment

Page 3 of 7 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry quantum chemical calculations (using B3LYP/6-31G(d) level of theory) (c) Plot for the scatter in the C1-C2 and C4-C5 bondlengths across the 60 sampled TS1 structures.

amongst 4 or 5. A generalized potential energy surface schematically represents the bifurcation pathway for deazetization of I (Scheme 2).

Normal mode sampling was conducted on the dynamical bottleneck structure at 300 K to obtain an ensemble of TS1 structures (See computational section for details). An overlay of 25 sampled structures for TS1 is shown in Figure 2(a). To verify the suitability of the sampling method, we have plotted the C1-C2 and C4-C5 bonds (See Figure 2(b) for atom numbering in the TS1 structure) of 60 sampled TS1 structures, which is represented in Figure 2(c). The transition zone is depicted schematically within an ellipsoid in which both bond lengths are 1.63 Å ± 0.01 Å. This region conserves > 98% of the transition state geometry. Total 303 trajectories were propagated from TS1 in both forward and backward directions of which 269 were productive. Reactant (3) and products (4, 5 and TS2) were defined by the choice of specific cut-offs: C1–C24.0 Å for 5 and C1–C2>2.1 Å, C4–C5>2.1 Å, C3–N16>4.0 Å and C6–N17>4.0 Å for TS2. Amongst the 269 productive trajectories, 91 trajectories reached 4, 88 trajectories reached 5 and rest (90) of them arrived at TS2. Out of the 34 unproductive trajectories, some of them recrossed back into the reactant or product side. Therefore, there is an approximately 1:1:1 ratio for the formations of 4, 5 and TS2 from the dynamical bottleneck structure (TS1), which might be due to equal dynamic racing of different modes while trajectory propagation.37 The plots of the trajectories (Figure 3) reveal that they move along TS1 → VRI and subsequently bifurcate to produce either 4 or 5 or TS2. From the IRC calculations one observes that though the minimum energy pathway connects the two saddle points (TS1 and TS2) yet it becomes unstable at the VRI point. Therefore, beyond VRI, the trajectories have a freedom to diverge towards the alternative product valleys (4 and 5). Snapshots for the formation of 4, 5 and TS2 from typical trajectories are shown in Figure S3. Interestingly, formations of these structures are concerted and proceed only through the critical structures (TS1, VRI point) without forming any intermediate. Figure 3(a) shows the trajectories that lead to TS2. Formation of “probable synthetic transition state” has been indicated in previous reports.11, 38-39 Figure 3(b) clearly depicts the production of the degenerate semibullvalenes (4 and 5) via a bifurcating pathway through VRI. Concomitant cleavage of both the C–N bonds (C3–N16 and C6– N17) occurs within ~15 fs during backward and forward trajectory propagations in a synchronous manner. Similarly, production of TS2 also follows a synchronous pathway wherein both C1–C2 and C4–C5 bonds break simultaneously. When both the C–N bonds are ~4.0 Å, formation of TS2, 4 and 5 were confirmed during trajectory simulations using the cut-off criteria listed above. Previous experimental studies by Carpenter et al. for the thermal deazetization of 2,3-Diazabicyclo[2.2.1]hept-2-ene also revealed a synchronous cleavage of two C-N bonds.40 Interestingly, the computed average time for the trajectories to arrive at TS2, 4 and 5 (Table 2) from 3 are similar. To explore the chemical importance of the VRI point, few trajectories were also simulated for V (R = CN), where two products (4 and 5) are inequivalent (activation energy 6.6 kcal/mol and 4.1 kcal/mol with respect to 4 and 5 respectively). Dynamic pictures clearly show the preference of 4 over 5. This dynamical preference can easily justifiable from their thermodynamic stability order (ΔΔG (4-5) = -2.6 kcal/mol). Nonetheless, by stereo-electronic switching one can manipulate the chameleonic behavior of the

Scheme 2: Bifurcation in the generalized potential energy surface for deazetization of I-IX producing degenerate semibullvalenes (4 and 5) and TS2. From the above IRC calculations in I-IX, it is evident that the balance between the three pathways (3 → TS1 → TS2 and 3 → TS1 → 4 or 3 → TS1 → 5) is actually very subtle and therefore, requires more detailed dynamical study. Additionally, since, neither of the species involved in this reaction are overtly polar (µ = 0–5 Debye), so polarity of the solvent(s) is unlikely to effect the overall mechanistic pathway (details are shown in Table S5–S6 of Supporting Information file). Hence, we believe that gas phase quasi-classical trajectory (QCT) analyses should indeed represent reasonable product distribution under ambient experimental conditions.

Figure 2: (a) Overlay of 25 sampled TS1 structures in I (R = H) (b) Optimized structure of TS1 for I (R = H) available from 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

and 14C isotopic substitution at various positions in 3 to tune the product selectivity in bifurcation (Figure 4 (A) – 4(D)). The IRCs are examined starting from the C2v saddle-point, TS1. For the dideuterated case (A), the IRC terminates at TS2 but for both B and C, the IRCs exclusively produce 4. The IRC expectedly prefers to propagate along the lighter atoms (C4 and C5) and hence, move more in this pathway to form the major product, 4. Even though this qualitative approach predicts the major product correctly yet unfortunately, it fails to indicate if the preference is modest or exclusive. To get a quantitative estimate for the kinetic selectivity, we have carried out quasi-classical trajectory calculations for all the four isotopomers. We have propagated 315, 359, 373 and 233 trajectories from TS1 of A, B, C and D respectively at 300 K (See Table 3). Products (4 and 5) were determined using the above mentioned cut-offs.

bifurcating point (See Supp. Info. for detailed dynamical outcomes).

Figure 3: Variations of C4-C5 vs C1-C2 bond lengths (in Å) of I (R = H), where (a) 30 representative trajectories for the formation of TS2. (b) 58 representative trajectories that result in formation of 4 and 5. The yellow and white asterisks represent TS1 and VRI respectively.

Species TS2 4 5

Page 4 of 7

Figure 4: Four isotopomers for 3 (R = H) on which direct dynamics calculations are performed.

Average time (fs) 116±14 112±14 116±11

Table 2: Average time-scales for producing either a product or TS2 starting from 3 (R = H). One simple test to experimentally determine the importance of post-transition state dynamics is to measure the kinetic isotope effects. 3 → TS1 being the rate-determining step would result in a kinetic isotope effect arising from ZPE and quantum mechanical tunneling (See Supp. Info. for computed rate constants, Table S7). Nevertheless, since the product selectivity occurs subsequent to TS1, kinetic isotope effects would be benign towards bifurcation7. Newtonian KIEs on the other hand, are indeed sensitive towards the motion of atoms along bifurcating potential surfaces.37, 41-42 With isotopic substitutions that make the products non-equivalent, the Newton’s second law would favor a product which arises through direct displacement along the steepest-descent path in mass-weighted coordinates (IRC). Herein, we used deuterium, 13C

Cases

A

B

C

D

No of trajectories

315

359

373

233

Trajectories affording 4 Trajectories affording 5 Isotopomer ratio (4:5) 95% confidence

141

179

196

134

137

146

143

84

1:1

1.2:1

1.4:1

1.6:1

1.3–0.8

1.4–1.1

1.5–1.2

1.8–1.4

Table 3: Results for the quasi-classical trajectories starting from TS1 of A, B, C and D respectively at B3LYP/6-31G(d) level. In case of A, product ratio between 4 and 5 still remains ~ 1:1 as that for the parent, non-isotopically labeled, I. This is expected as H/D replacement occurs only at the distal position. 4

ACS Paragon Plus Environment

Page 5 of 7 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry AUTHOR INFORMATION

Interestingly, isotopic effects are anticipated to be large for substitutions only at the C1-C2 and C4-C5 positions. For a single 12 13 C/ C isotopic replacement (B), 4:5 is 1.2:1 which enhances to 1.4:1 in C for two 12C/13C isotopic replacements (along C1-C2). As a proof of principle, two 12C/14C isotopic replacements in the C1-C2 bond makes 4:5 as high as 1.6.43 This clearly indicates that pure Newtonian effects govern the KIEs for these deazetization reactions (for further details see the supporting information file). RRKM lifetimes of each case are also analyzed (Table 4) from which we observe that, an increment in mass can retard the extent of rearrangement between 4 and 5. Since, the products are stable till ~ 300 fs, one can experimentally observe the non-statistical population ratios (4:5) using a femtosecond time resolved transient 2D-IR (t-2D-IR) technique to probe the initial product ratios. We believe that, high sensitivity of t-2D-IR spectroscopy can provide a useful basis for characterizing these dynamical effects generated from their nuclear motions44-48 (for further details see the supporting information file).

Corresponding Author Corresponding 24734971.

Author:

[email protected].

Phone:

+91-33-

Notes The authors declare no competing financial interests. ACKNOWLEDGMENT NM thanks CSIR for JRF. AD thanks DST and BRNS for partial funding. We thank Professor Daniel A. Singleton for useful discussions regarding the use and applications of the Progdyn program. Dedicated to Professor Roald Hoffmann on his 80th birthday.

REFERENCES 1.

Molecules

Lifetime of 4 (fs)

Lifetime of 5 (fs)

2.

Parent (R =H)

314

314

3.

A (R = H)

327

327

B (R = H)

316

316

C (R = H)

321

321

D (R = H)

328

326

4.

5. 6.

Table 4: Computed RRKM lifetimes (in fs) of 4 and 5 for different isotopomers (A, B, C, D respectively) at 300K using B3LYP/631G(d) level of theory.

CONCLUSIONS

7.

In summary, we have shown that semibullvalene formation from 3 is concerted, synchronous, featuring a bispericyclic transition state (TS1) and is dictated by dynamical effects. Experimental evidence for post-transition-state dynamic control can be readily obtained by the observation of Newtonian kinetic isotope effects in the ratio of 4 and 5. We hope that our computational prediction would be subject to experimental test soon.

8.

9. 10.

ASSOCIATED CONTENT Supporting Information Method for quasi-classical trajectory (QCT) calculations, input parameters for Progdyn, normal mode sampling, snapshots from trajectories, variation of bond lengths obtained from QCTs, lifetimes, detailed analyses of dynamic selectivity study, uncertainties in the trajectory ratios, Simulations of V (R = CN), DFT calculations, VRI-point analysis, IRC for R = Cl from TS1, IRCs for substituted case, relative energies and stability analyses of different cases, implicit solvation study, computed rate constants, RRKM lifetimes, optimized cartesian coordinates and harmonic frequencies.

11. 12.

13.

Modern Physical Organic Chemistry, Anslyn, E. V.; Dougherty, D. A. University Science Books, USA (2006). Truhlar, D. G.; Garrett, B. C. Acc. Chem. Res, 1980, 13, 440448. Singleton, D. A.; Hang, C.; Szymanski, M. J.; Meyer, M. P.; Leach, A. G.; Kuwata, K. T.; Chen, J. S.; Greer, A.; Foote, C. S.; Houk, K. N. Mechanism of Ene Reactions of Singlet Oxygen. A Two-Step No-Intermediate Mechanism J. Am. Chem. Soc. 2003, 125, 1319−1328. Bekele, T.; Christian, C. F.; Lipton, M. A.; Singleton, D. A. “Concerted” Transition State, Stepwise Mechanism. Dynamics Effects in C2-C6 Enyne Allene Cyclizations J. Am. Chem. Soc. 2005, 127, 9216. Litovitz, A. E.; Keresztes, I.; Carpenter, B. K. Evidence for Nonstatistical Dynamics in the Wolff Rearrangement of a Carbene J. Am. Chem. Soc. 2008, 130, 12085−12094. Glowacki, D. R.; Marsden, S. P.; Pilling, M. J. Significance of Nonstatistical Dynamics in Organic Reaction Mechanisms: Time-Dependent Stereoselectivity in Cyclopentyne−Alkene Cycloadditions J. Am. Chem. Soc. 2009, 131, 13896−13897. Wang, Z.; Hirschi, J. S.; Singleton, D. A. Recrossing and Dynamic Matching Effects on Selectivity in a Diels–Alder Reaction Angew. Chem., Int. Ed. 2009, 48, 9156. Patel, A.; Chen, Z.; Yang, Z.; Gutierrez, O.; Liu, H.; Houk, ́ K. N.; Singleton, D. A. Dynamically Complex [6+4] and [4+2] Cycloadditions in the Biosynthesis of Spinosyn A J. Am. Chem. Soc. 2016, 138, 3631. Hornsby, C. E.; Paton, R. S. Natural product biosynthesis: It’s all downhill from here. Nat. Chem. 2014, 6, 88−89. Martín-Sómer, A.; Yáñez, M.; Hase, W. L.; Gaigeot, M.-P.; Spezia, R. Post-Transition State Dynamics in Gas Phase Reactivity: Importance of Bifurcations and Rotational Activation J. Chem. Theory Comput. 2016, 12, 974−982. Zhou, C.; Birney, D. Sequential Transition States and the Valley−Ridge Inflection Point in the Formation of a Semibullvalene Org. Lett. 2002, 4, 3279−3282. Boger, D. L.; Panek, J. S.; Duff, S. R. Inverse electron demand Diels-Alder reactions of heterocyclic azadienes: Formal total synthesis of streptonigrin J. Am. Chem. Soc. 1985, 107, 5745−5754. Benson, S. C.; Palabrica, C. A.; Snyder, J. K. Indole as a Dienophile in Inverse Electron Demand Diels-Alder Reactions. 5H-Pyridazino[4,5-b]indoles as Cycloadducts with 3,6-Dicarbomethoxy- 1,2,4,5-tetrazine J. Org. Chem. 1987, 52, 4610−4614.

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

14. Sauer, J.; Heldmann, D. K.; Hetzenegger, J.; Krauthan, J.; Sichert, H.; Schuster, 1,2,4,5-Tetrazine: Synthesis and Reactivity in [4+2] Cycloadditions. Eur. J. Org. Chem. 1998, 1998, 2885−2896. 15. Boger, D. L.; Boyce, C. W.; Labroli, M. A.; Sehon, C. A.; Jin, Q. Total Syntheses of Ningalin A, Lamellarin O, Lukianol A, and Permethyl Storniamide A Utilizing Heterocyclic Azadiene Diels−Alder Reactions J. Am. Chem. Soc. 1999, 121, 54−62. 16. Sauer, J.; Bäuerlein, P.; Ebenbeck, W.; Schuster, J.; Sellner, I.; Sichert, H.; Stimmelmayr, H. An One-Pot Synthesis of Semibullvalenes and Its Mechanism Eur. J. Org. Chem. 2002, 791-801. 17. Patterson, D. M.; Nazarova, L. A.; Xie, B.; Kamber, D. N.; Prescher, J. A. Functionalized cyclopropenes as bioorthogonal chemical reporters J. Am. Chem. Soc. 2012, 134, 18638−18643. 18. Yang, J.; Šeĉkutė, J.; Cole, C. M.; Devaraj, N. K. Live-Cell Imaging of Cyclopropene Tags with Fluorogenic Tetrazine Cycloadditions Angew. Chem., Int. Ed. 2012, 51, 7476−7479. 19. Thomas, J. D.; Cui, H.; North, P. J.; Hofer, T.; Rader, C.; Burke, T. R. Application of strain-promoted azide–alkyne cycloaddition and tetrazine ligation to targeted Fc-drug conjugates Bioconjugate Chem. 2012, 23, 2007−2013. 20. Devaraj, N. K. Advancing tetrazine bioorthogonal reactions through the development of new synthetic tools Synlett 2012, 23, 2147−2152. 21. Blackman, M. L.; Royzen, M.; Fox, J. M. Tetrazine ligation: fast bioconjugation based on inverse-electron-demand Diels− Alder reactivity J. Am. Chem. Soc. 2008, 130, 13518−13519. 22. Devaraj, N. K.; Weissleder, R.; Hilderbrand, S. A. Tetrazinebased cycloadditions: application to pretargeted live cell imaging Bioconjugate Chem. 2008, 19, 2297−2299. 23. Törk, L.; Jimenez-Osés, G.; Doubleday, C.; Liu, F.; Houk, K. N. Molecular Dynamics of the Diels–Alder Reactions of Tetrazines with Alkenes and N2 Extrusions from Adducts J. Am. Chem. Soc. 2015, 137, 4749. 24. Ussing, B. R.; Hang, C.; Singleton, D. A. Dynamic Effects on the Periselectivity, Rate, Isotope Effects, and Mechanism of Cycloadditions of Ketenes with Cyclopentadiene J. Am. Chem. Soc. 2006, 128, 7594−7607. 25. Thomas, J. B.; Waas, J. R.; Harmata, M.; Singleton, D. A. Control Elements in Dynamically Determined Selectivity on a Bifurcating Surface J. Am. Chem. Soc. 2008, 130, 14544−14555. 26. Yu, P.; Chen, T.Q.; Yang, Z.; He, C.Q.; Patel, A.; Lam, Y.; Liu, C.; Houk, K.N. Mechanisms and Origins of Periselectivity of the Ambimodal [6 + 4] Cycloadditions of Tropone to Dimethylfulvene J. Am. Chem. Soc. 2017, 139, 8251-8258. 27. Hare, S.R.; Pemberton, R.P.; Tantillo, D.J. Navigating Past a Fork in the Road: Carbocation–π Interactions Can Manipulate Dynamic Behavior of Reactions Facing PostTransition-State Bifurcations J. Am. Chem. Soc. 2017, 139, 7485–7493. 28. Sadasivam, D.V.; Prasad, E.; Flowers R.A.; Birney D.M. Stopped-Flow Kinetics of Tetrazine Cycloadditions; Experimental and Computional Studies toward Sequential Transition States J. Phys. Chem. A, 2006, 110, 1288–1294. 29. Yamamoto, Y.; Hasegawa, H.; Yamataka, H. Dynamic Path Bifurcation in the Beckmann Reaction: Support from Kinetic Analyses J. Org. Chem. 2011, 76, 4652. 30. López, R.V.; Faza, O.N.; López, C.S. Dynamic Effects Responsible for High Selectivity in a [3,3] Sigmatropic

31.

32.

33.

34.

35. 36.

37.

38.

39.

40.

41.

42.

43.

44.

45.

Rearrangement Featuring a Bispericyclic Transition State J. Org. Chem. 2017, 82 (9), 4758–4765. Ess, D. H.; Wheeler, S. E.; Iafe, R. G.; Xu, L.; CelebiOlcum, N.; Houk, K. N. Bifurcations on Potential Energy Surfaces of Organic Reactions Angew. Chem. Int. Ed. 2008, 47, 7592-7601. Hare, S. R.; Tantillo, D. J. Post-transition state bifurcations gain momentum – current state of the field Pure Appl. Chem. 2017, 89, 679-698. Gaussian 09, Revision E.01, Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G.E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A. et al. Gaussian, Inc., Wallingford CT, 2009. Duncan, W. T.; Bell, R. L.; Truong, T. N. TheRate: program for ab initio direct dynamics calculations of thermal and vibrational-state-selected rate constants. J. Comput. Chem. 1998, 19, 1039−1052. https://cdssim.chem.ttu.edu/htmlpages/licensemenu.jsp Zhu, L.; Hase, W. L., A General RRKM Program. Quantum Chemistry Program Exchange; Indiana University: Bloomington, IN, 1994; Program No. QCPE 644. Andujar-De Sanctis, I. L.; Singleton, D. A. Racing Carbon Atoms. Atomic Motion Reaction Coordinates and Structural Effects on Newtonian Kinetic Isotope Effects Org. Lett. 2012, 14, 5238−5241. Toma, L.; Romano, S.; Quadrelli, P.; Caramella, P. Merging of 4+2 and 2+4 cycloaddition paths in the regiospecific dimerization of methacrolein. A case of concerted cryptodiradical cycloaddition Tetrahedron Lett. 2001, 42, 50775080. Caramella, P.; Quadrelli, P.; Toma, L. An Unexpected Bispericyclic Transition Structure Leading to 4+2 and 2+4 Cycloadducts in the Endo Dimerization of Cyclopentadiene J. Am. Chem. Soc. 2002, 124, 1130-1131. Reyes, M. B.; Carpenter, B. K. Mechanism of the Thermal Deazetization of 2,3-Diazabicyclo[2.2.1]hept-2-ene and Its Reaction Dynamics in Supercritical Fluids J. Am. Chem. Soc. 2000, 122, 10163-10176. Singleton, D. A.; Hang, C.; Szymanski, M. J.; Greenwald, E. E. A New Form of Kinetic Isotope Effect. Dynamic Effects on Isotopic Selectivity and Regioselectivity J. Am. Chem. Soc. 2003, 125, 1176−1177. Kelly, K. K.; Hirschi, J. S.; Singleton, D. A Newtonian Kinetic Isotope Effects. Observation, Prediction, and Origin of Heavy-Atom Dynamic Isotope Effect. J. Am. Chem. Soc. 2009, 131, 8382–8383. 14 C is a radioactive isotope and hence, it is rather unlikely to be used for kinetic isotope effects studies. Nevertheless, the increasing ratio of 4:5 as a function of replacement of C1 or/and C2 by the heavier isotopes acts as a diagnostic proof for Newton’s 2nd Law governed kinetic isotope effects in I. Carpenter, B. K.; Harvey, J. N.; Orr-Ewing, A. J. The Study of Reactive Intermediate in Condensed Phases. J. Am. Chem. Soc. 2016, 138, 4695–4705. Schleif, T.; Mierez-Perez, J.; Henkel, S.; Ertelt, M.; Borden, W.T.; Sander, W. The Cope Rearrangement of 1,5Dimethylsemibullvalene-2(4)-d1: Experimental Evidence for Heavy-Atom Tunneling. Angew. Chem., Int. Ed. 2017, 56, 10746–10749.

6

ACS Paragon Plus Environment

Page 6 of 7

Page 7 of 7 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry 46. Asvany, O.; Kumar, P.; Redlich, B.; Hegemann, I.; Schlemmer, S.; Marx, D. Understanding the Infrared Spectrum of Bare CH5+. Science 2005, 309, 1219–1222. 47. Zewail, A. H. Femtochemistry: Atomic-Scale Dynamics of the Chemical Bonds. J. Phys. Chem. A. 2000, 104, 5660– 5694. 48. Bredenbeck, J.; Helbing, J.; Hamm, P. Transient twodimensional infrared spectroscopy: Exploring the polarization dependence. J. Chem. Phys. 2004, 121, 5943.

TOC Graphic

7

ACS Paragon Plus Environment