Catalytic Enantioselective [2,3]-Rearrangements of Allylic Ammonium Ylides: A Mechanistic and Computational Study.

A mechanistic study of the isothiourea-catalyzed enantioselective [2,3]-rearrangement of allylic ammonium ylides is described. Reaction kinetic analyses using 19F NMR and density functional theory computations have elucidated a reaction profile and allowed identification of the catalyst resting state and turnover-rate limiting step. A catalytically relevant catalyst-substrate adduct has been observed, and its constitution elucidated unambiguously by 13C and 15N isotopic labeling. Isotopic entrainment has shown the observed catalyst-substrate adduct to be a genuine intermediate on the productive cycle toward catalysis. The influence of HOBt as an additive upon the reaction, catalyst resting state, and turnover-rate limiting step has been examined. Crossover experiments have probed the reversibility of each of the proposed steps of the catalytic cycle. Computations were also used to elucidate the origins of stereocontrol, with a 1,5-S···O interaction and the catalyst stereodirecting group providing transition structure rigidification and enantioselectivity, while preference for cation-π interactions over C-H···π is responsible for diastereoselectivity.

Introduction

The [2,3]-rearrangement of allylic ammonium ylides is a direct and elegant method toward the synthesis of α-amino acid derivatives containing multiple stereocenters.[1] The mechanism of this process, and that of the competitive [1,2]-Stevens rearrangement, has been much discussed and disputed within the literature. A concerted thermally allowed sigmatropic process is thought to be operative in the [2,3]-rearrangement, while a radical mechanism involving bond cleavage and recombination is usually favored for [1,2]-rearrangement (Scheme A).[2]

Scheme 1

Rearrangements of Allylic Ammonium Ylides

To date, few mechanistic analyses of [2,3]-rearrangements of allylic ammonium ylides have been conducted, although Jacobsen and co-workers have recently reported a detailed mechanistic investigation into the related thiourea-catalyzed [2,3]-Wittig rearrangement.[3] The elegant experimental and computational work of Singleton and co-workers concerning the competitive [2,3]- and [1,2]-rearrangements of allylic ammonium ylides promoted by DBU represents the current state-of-the-art (Scheme B). Through 13C kinetic isotope effects, crossover experiments, and computation, these studies demonstrate that the origin of competitive [1,2]- and [2,3]-rearrangement is the common loose transition state leading to dynamic bond cleavage.[4] The development of both catalytic and stereoselective variants of the [2,3]-rearrangement of allylic ammonium ylides has been a significant synthetic challenge. Tambar and co-workers have reported a tandem ammonium salt formation and diastereoselective [2,3]-rearrangement process, exploiting Pd-catalyzed allylic substitution to form the reactive ammonium salt in situ, giving (±)-anti-α-amino acid derivatives with excellent diastereocontrol via proposed transition state 13 (Scheme A).[5] The observed diastereoselectivity of this [2,3]-rearrangement process, and indeed most [2,3]-rearrangements, can be rationalized through the exo- or endo- transition states 15 and 17 initially described by Houk and Marshall for the related [2,3]-Wittig rearrangement (Scheme B).[6]

Scheme 2

Stereochemical Models for [2,3]-Rearrangements

Prior to our studies within this area, only limited methods capable of imparting enantiocontrol in the [2,3]-rearrangement of allylic ammonium ylides had been developed. Sweeney first demonstrated a chiral auxiliary approach to allow access to enantiomerically enriched α-amino acid derivatives,[7] while the use of a superstoichiometric chiral Lewis acid promoter was subsequently reported by Somfai.[8] Catalytic enantioselective variants were unknown until 2014, when our laboratory reported an isothiourea-catalyzed[9] [2,3]-rearrangement of allylic quaternary ammonium salts to give syn-α-amino acid derivatives with excellent levels of diastereo- and enantiocontrol (Scheme ).[10] Treatment of quaternary ammonium salts 19 bearing an activated p-nitrophenol ester, either isolated or generated in situ, with catalytic (+)-benzotetramisole ((+)-BTM) 20, cocatalytic hydroxybenzotriazole (HOBt), and iPr2NH, gave stereoselective [2,3]-rearrangement into syn-α-amino acid derivatives with excellent levels of stereocontrol. This process can be performed in the absence of HOBt; however, its addition provides a subtle enhancement in both diastereo- and enantioselectivity. We tentatively proposed a Lewis base catalytic cycle, initiated by nucleophilic addition of (+)-BTM 20 into the activated ester substrate to form an acyl ammonium intermediate prior to the formation of ammonium ylide 22. However, alternative mechanistic pathways using either Lewis or Brønsted base catalysis proceeding via different intermediates can be envisaged (Scheme ). For example, assuming Lewis base catalysis is operative, the reaction could proceed through initial formation of a ketene intermediate 23 en route to acyl ammonium ylide 22. Furthermore, the origin of the observed diastereo- and enantiocontrol in the rearrangement process is currently unknown.

Scheme 3

Catalytic Enantioselective [2,3]-Rearrangement

Herein we report experimental and computational investigations into the mechanism and origins of stereocontrol in the isothiourea-catalyzed [2,3]-rearrangement of allylic ammonium ylides (Scheme C). In situ NMR analysis has allowed a reaction profile to be elucidated, while isotopic-labeling studies have unambiguously identified a genuine productive catalytic intermediate. Kinetic analysis has given insight into the overall process, and crossover studies have provided information about the reversibility of each step. Kinetic isotope effects have also been used to probe the stereodetermining [2,3]-rearrangement step of the process. Computational reaction coordinate modeling provides deeper insight into the catalytic cycle, and transition state modeling reveals the origins of stereochemical control.

Results and Discussion

Mechanistic Studies

Temporal Concentration Profiles

Initial studies aimed to establish the kinetics of the [2,3]-rearrangement and identify any reaction intermediate(s) or catalyst resting states. Ammonium salt 25a (34 mM) rearranges to 26a in d3-MeCN/d6-DMSO (9:1) at −20 °C, catalyzed by (+)-BTM 20 (20 mol %) (Scheme ). The presence of various salts in the reaction medium causes extensive line broadening in the 1H NMR spectrum, making it unsuitable for in situ analysis of the [2,3]-rearrangement. However, the 19F{1H} and 13C{1H} NMR spectra were tractable, and the 4-fluoro substituent in 25a allows quantitative monitoring of the process by in situ19F{1H} NMR (δF = −113.1 ppm), with PhCF3 as an internal standard.[11] After an initial burst phase (<1000 s), ammonium salt 25a is converted into 26a (>80% 26a, δF = −117.6 ppm) with pseudo-first-order kinetics, vide infra, over a period of 4 h. During the reaction evolution, a transient species (δF = −117.0 ppm) was detected, accumulating to a maximum concentration of ∼5.2 mM in the early stages of catalysis and then depleting as substrate 25a was consumed. Rearrangement in the absence of HOBt resulted in a similar reaction profile, but afforded higher concentrations of the same species (δF = −117.0 ppm), over a longer period, assisting in its analysis.

Scheme 4

System Chosen for in Situ19F NMR Study

Identification of an On-Cycle Catalytic Intermediate

Catalyst Speciation

To determine if the transient species (δF = −117.0 ppm) involves the catalyst, a fluorinated variant,[12] (+)-F-BTM 27, was prepared.[11]In situ monitoring of the [2,3]-rearrangement of 25a, catalyzed by 27 (Figure ), confirmed conversion of free (+)-F-BTM 27 (δF = −123.3 ppm) into a catalyst-derived species (δF = −113.4 ppm), which was also transient, reaching a maximum concentration of 4.6 mM at ∼2500 s, and decaying as the reaction proceeds to completion (see inset graph to Figure ). The comparable temporal intensities of the two signals (δF = −117.0 ppm and −113.4 ppm) strongly suggest they arise from a single transient intermediate containing both 25a and 27 in a formal 1:1 combination. Based on reference to isothiouronium salt 28 (δF = −113.3 ppm), the 19F chemical shift of the catalyst-derived component in transient species 29 suggested it to be an N-acylated isothiourea.

Figure 1

Temporal concentration data for [2,3]-rearrangement of 25a using (+)-F-BTM 27. Conditions: 25a (29.5 mM), (+)-F-BTM 27 (6.8 mM), iPr2NH (47 mM), d3-MeCN/d6-DMSO (9:1), −20 °C. Inset: monitoring of catalyst-derived species.

Atom Connectivity: 13C/15N Labeling

Isotopically labeled substrates and catalyst (1-[13C1]-25a, 1,2-[13C2]-25a, 1,2,3′-[13C3]-25b (Figure A) and (±)-[15N1]-27) were prepared[11] to deduce connectivity between various atoms in the intermediate and probe for reversibility in its generation. Rearrangement of 1,2-[13C2]-25a (34 mM) catalyzed by (±)-[15N1]-27 (6.8 mM, 20 mol %) was monitored by 13C{1H} NMR (Figure B). Reducing the iPr2NH concentration to 23.8 mM prolonged the lifetime of the intermediate, [15N1,13C2]-29a, allowing detailed analysis of the 13C=O region, Figure C. The characteristic doublets (1JCC = 52 Hz) arising from the adjacent 13C labels, C(1)–C(2), in the substrate (1,2-[13C2]-25a) and the product (1,2-[13C2]-26a) are replaced by a double–doublet in [15N1, 13C2]-29a (δC = 173.1 ppm).[13] The magnitude of the C–C coupling (1JCC = 52 Hz) confirms that C(2) remains sp3-hybridized. The magnitude of the additional coupling constant (1JCN = 5.2 Hz)[13] indicates that C(1) is directly bound to the 15N-labeled atom in the catalyst, confirming 29a to be an N-acylated isothiourea (Figure ).

Figure 2

(A) Isotopically labeled substrates. (B) Generation of [15N1, 13C2]-labeled intermediate 29a. (C) 13C{1H} NMR subspectrum (101 MHz, d3-MeCN/d6-DMSO (9:1), 273 K) of the 13C=O region. Conditions: 1,2-[13C2]-25a (34 mM), (±)-15N-F-BTM [15N1-27] (6.8 mM), iPr2NH (23.8 mM).

Attempts to identitify the atom adjacent to C(2) in [15N1, 13C2]-29a from its 1H-coupled 13C NMR signal were thwarted by line broadening.[13] Instead, the 13C{1H} NMR spectrum of intermediate 1,2,3-[13C3]-29b, generated from 1,2,3′-[13C3]-25b, was analyzed (Figure ). The C(2) signal of the resulting [2,3]-rearrangement product 1,2,3-[13C3]-26b (δC = 70.5 ppm) displayed the expected double–doublet coupling arising from C(2)–C(3) bond formation. However, this coupling pattern was also evident in intermediate 1,2,3-[13C3]-29b (δC = 68.0 ppm, 1JCC 51.4 Hz, 1JCC 35.8 Hz) with the magnitude of the C(2)–C(3) coupling indicative of sp3-hybdrization at both centers.[13] Overall the data confirms that the intermediate (29) is a catalyst-bound, post-[2,3]-rearrangement, acyl ammonium salt (Figure ).

Figure 3

Catalyst-bound [2,3]-rearrangement product, 13C{1H} NMR subspectrum (101 MHz, d3-MeCN/d6-DMSO (9:1), 293 K) of the 13C(2)-H region. Conditions: 1,2,3′-[13C3]-25b (34 mM), (+)-F-BTM 27 (6.8 mM), iPr2NH (23.8 mM).

Productivity of the Intermediate

The data presented so far does not discriminate between the N-acylated isothiourea species (29) being peripheral to the productive catalytic cycle (Figure , case A) or an integral part of it (cases B–D). The following isotopic entrainment test distinguishes these four possibilities. A catalytic reaction employing 1-[13C1]-25a (17 mM) was allowed to evolve until 1-[13C1]-29a had reached its maximum concentration (∼5 mM). A further 1.0 equiv (17 mM) of a differently labeled substrate, 1,2-[13C2]-25a, was then rapidly added, resulting in an isotopic perturbation of the system. At the point that 1,2-[13C2]-25a is added, there has been 32% net conversion of 25a ([13C1]- and [13C2]-). However, neither of the [2,3]-rearrangement products (29a and 26a) yet contain any of the [13C2]-label: all of this resides in unreacted [13C1,2]-25a, which comprises 0.28 [13C1]/0.72 [13C2]. The key features are the changes in 13C2-populations in the substrate 25a, intermediate 29a and product 26a, as the reaction evolves. Irrespective of the pathway (A–D) the population in the final product (26a) must ultimately rise from 0% to 50%, as dictated by the equal proportions of 1-[13C1]-25a and 1,2-[13C2]-25a added overall. For case A, where 29a is not productive, the isotope population in 29a will depend only on that of the final product (26a; max 50% 13C2) and at all stages will be lower or equal to it. For case B, where the intermediate is productive, but is in equilibrium with 25a, the 13C2-population in 25a will be reduced, in the limit from 72% to 55%. For cases C and D, where the [2,3]-rearrangement to 29a is irreversible, the isotope population in 25a is constant (72% 13C2-) and the 13C2-content in 29a rises from 0% to a maximum of 72% as it is repopulated from 25a.[14] However, for case C, equilbration of 29a with product 26a will attenuate the rise in 13C2-population in 29a, in the limit to 50%. Only for case D will the 13C2 isotope population in 29a rise, in advance of 26a, to reach a maximum 72% 13C2.[14] Comparison of the predicted and experimentally determined 13C2-populations as a function of net conversion (Figure ) confirms that 26 arises from two irreversible sequential first-order interconversions (25 → 29 → 26) where 29 is the productive catalytic intermediate (case D).[14] Kinetic modeling confirms that the impact of heavy-atom (12C/13C) KIEs on the isotope-entrainment are negligible.[11]

Figure 4

13C2 Entrainment into the catalytic cycle. Conditions: 1-[13C1]-25a (17 mM), (+)-BTM (6.8 mM), iPr2NH (47.6 mM), d3-MeCN/d6-DMSO (9:1), with 1,2-[13C2]-25a (17 mM) added at 32% net conversion of 25a. Open circles: experimental (13C{1H} NMR) data for 13C2-incorporation (%) versus net conversion (%). Dashed lines: kinetic simulation where 29 is a productive intermediate in two irreversible sequential pseudo-first-order interconversions (25 → 29 → 26, Case D; with rate ratio 0.407).[14]

HOBt Cleavage of Acyl Ammonium Intermediate

HOBt provides optimal diastereo- and enantiocontrol (Scheme ); however, its role within the catalytic cycle is unclear. To probe if the HOBt enhances stereocontrol via suppression of the base-mediated background reaction, the iPr2NH-mediated rearrangement of 25a was examined. Reaction of 25a with iPr2NH (47 mM) in the absence of the BTM catalyst resulted in slow formation of racemic 26a with low diastereocontrol (79:21 dr) and a kobs of 2.74 × 10–5 s–1. The addition of a catalytic amount of HOBt (6.8 mM) resulted in no change in rate (kobs 2.70 × 10–5 s–1). This rules out the role of HOBt as improving stereocontrol through direct suppression of the rate of the background reaction.

The (+)-BTM-catalyzed rearrangement of 25a in the presence of stoichiometric HOBt (34 mM) was studied by 19F{1H} NMR. The presence of HOBt strongly suppressed accumulation of acyl ammonium intermediate 29c ((+)-BTM replaces (+)-F-BTM) and resulted in the formation of the corresponding HOBt ester 30 (δF = −117.2 ppm, confirmed by comparison with an authentic sample),[11] in addition to the PNPO ester product 26a. Addition of HOBt once the acyl ammonium intermediate 29c reached the pseudo steady-state (5 mM, t = 3920 s, Figure ) resulted in immediate formation of HOBt ester 30 and consumption of acyl ammonium intermediate 29c. HOBt thus shifts the catalyst speciation to be strongly dominated by free (+)-BTM (Figure ). The background iPr2NH-mediated reaction that converts 25 into racemic (±)-26 presumably involves the generation of ylide intermediate 24. Interception of this species by free (+)-BTM would also generate intermediate 29, thus leading to nonracemic 26. The higher the concentration of free BTM, the more effective this interception.

Figure 5

Addition of stoichiometric HOBt at t = 3920 s. Conditions: 25a (34 mM), (+)-BTM (6.8 mM), iPr2NH (47.6 mM), d3-MeCN/d6-DMSO (9:1), HOBt (34 mM). Total Product refers to both HOBt ester 30 and 26a.

Overall, the beneficial effect of HOBt on the selectivity may arise from both a change in catalyst speciation to favor free BTM and by diversion of the background reaction onto the enantioselective pathway (Scheme ).

Scheme 5

Effect of HOBt on Overall Catalytic Cycle

Reaction Kinetics and Impact of Additives

Having identified, by 19F{1H} and 13C{1H} NMR, the major reactant-derived and catalyst-derived components present in the reaction mixture, the empirical rate equation was established by analysis of the decay in substrate 25a during the pseudosteady-state phase of the catalysis (Scheme ).[15] The standard conditions [25a (34 mM), (+)-BTM (6.8 mM), HOBt (6.8 mM), and iPr2NH (47 mM)], afforded a pseudo-first-order rate constant, kobs = 1.37 × 10–4 s–1.

Scheme 6

Empirical Rate Equation,

d[26a]/dt = kobs[25a]; kobs = k[BTM][NR3].

Conditions: 25a (17–68 mM), (+)-BTM (1.7–13.6 mM), HOBt (0–34 mM), iPr2NH (23.8–95.2 mM), d3-MeCN/d6-DMSO (9:1), 253 K, 4 h.

There was a linear relationship between the enantiopurity of (+)-BTM and product 26a, consistent with predominant or exclusive speciation of the catalyst in an active, monomeric, form.[11] Systematic variation of the concentration of the reaction components afforded empirical first-order dependencies on (+)-BTM (1.7–13.6 mM) and diisopropylamine (23.8–95.2 mM), with no rate impact from the HOBt (0–34 mM). Addition of n-Bu4N 4-nitrophenoxide (27.2 mM) in the absence of HOBt resulted in a large increase in rate (kobs = 2.42 × 10–4 s–1) and reduced the accumulation of the [2,3]-rearrangement acyl ammonium intermediate 29c (≤2 mM). Control studies showed that n-Bu4NBr (27.2 mM) resulted in a similar rate enhancement (kobs = 2.35 × 10–4 s–1) but did not suppress the accumulation of 29c (5 mM). The observed increase in kobs may be rationalized by an increase in ionic strength of the medium.[11] The n-Bu4N 4-nitrophenoxide thus promotes catalyst turnover of acyl ammonium 29c into product 26a and (+)-BTM. Addition of 4-nitrophenol (0–34 mM) resulted in a decrease in rate, with a negative first-order dependency. Overall, this suggests that 4-nitrophenoxide or benzotriazololate is required for efficient turnover of acyl ammonium 29c into products 26a/30 and (+)-BTM.

Secondary Kinetic Isotope Effect (SKIE)

Although it is clear that the [2,3]-rearrangement to generate 29 is irreversible, it is not evident whether this is the product-determining step, i.e. the first irreversible step in the cycle. During the [2,3]-rearrangement step, the carbon that becomes C(3)–H in 29 undergoes rehybridization from sp2 to sp3. The process to generate a C(3)–D isotopologue of 29 would be thus expected to exhibit a 2H-SKIE. Alternatively, if deprotonation at C(2)–H to generate an allylic ammonium ylide (Scheme ) was the product-determining step, there would not be any significant SKIE, as C(3)–D is remote.

The SKIE was measured by competition using aryl-D1-C(3)-D0-25a and aryl-D0-C(3)-D1-25a, in a double-labeling method[16] in which the C(3)–D/C(3)–H ratio as a function of fractional conversion is determined by 19F NMR (ΔδF = 0.28 ppm; aryl-D0/aryl-D1). After correction for the effect of aryl deuteration,[17,18] a value of kH/kD = 1.031 was obtained. The presence of a small positive SKIE is consistent with a product-determining [2,3]-rearrangement transition state (Scheme ). A linear free energy relationship analysis of a range of C(3)-aryl substrates, against standard Hammett sigma values, showed the C(3) position to be relatively insensitive to electronic substituent effects.[11]

Scheme 7

C(3) SKIE Competition Experiment

Crossover and Reversibility Studies

We have previously demonstrated the [2,3]-rearrangement step to be intramolecular and irreversible.[10] To distinguish which steps prior to the [2,3]-rearrangement are reversible, a crossover reaction between ammonium salts 25a and 31 (1:1) bearing two distinct activated ester groups (4-NO2C6H4 and 3,5-(CF3)2C6H3) and two distinct C(3)-aryl units (4-FC6H4 and 3-FC6H4) was monitored in situ under catalytic conditions (Scheme A). Complete equilibration with ammonium salts 32 and 33 was observed,[11] consistent with reversible generation of nonrearranged (+)-BTM acyl ammonium intermediates. To examine reversibility at the deprotonation step, the [2,3]-rearrangement of α-dideuterio ammonium salt α-[D2]-25 (75% D2) was monitored in situ. The product was obtained with significantly lower deuterium incorporation (29% D), consistent with a reversible deprotonation step (Scheme B).[19] Reaction of ammonium salt 33 in the presence of rearrangement product 26b (Scheme C) bearing distinct C(3)-aryl units and activated esters demonstrated no crossover, consistent with catalyst turnover being irreversible, confirming the conclusions deduced from isotopic entrainment.

Scheme 8

Crossover Studies

(i) (+)-BTM (20 mol %), iPr2NH (1.4 equiv), 25a (0.5 equiv), 31 (0.5 equiv), d3-MeCN/d6-DMSO (9:1), 253 K, 4 h; (ii) (+)-BTM (20 mol %), HOBt (20 mol %), iPr2NH (1.4 equiv), 26b (0.5 equiv), 33 (0.5 equiv), MeCN, 253 K.

Computational Studies

Computed Catalytic Cycle

We computed all intermediates, transition structures (TSs), and possible salt complexes involved in the catalytic cycle shown in Scheme . Exhaustive searches were performed to locate all pertinent conformations. Geometries and thermodynamic corrections were computed at the M06-2X[20]/6-31G(d)[21] level of theory.[22] Vibrational frequencies and thermal corrections to the Gibbs free energy were calculated at −20 °C and 1 atm to match the experimental conditions. Further energy refinements were completed using M06-2X/6-311++G(2df,p).[23] Implicit solvation corrections were applied using the polarized continuum model (PCM)[24] with UFF radii for acetonitrile in both the geometry optimizations and the single-point energy refinements. The hybrid meta-GGA functional M06-2X is generally more robust than B3LYP at accounting for dispersion and nonbonding interactions routinely found in organocatalytic reactions.[25] Kinetic isotope effects were calculated using the theory of Bigeleisen and Mayer[26] along with the rigid-rotor harmonic oscillator approach (ΔHΔS).[27] Quantum mechanical tunneling effects were also calculated for both methods using the one-dimensional parabolic approximation.[28] The calculation of the KIE was automated by use of the Onyx isotope effect program.[29] The catalytic cycle and the computed reaction coordinate are summarized in Scheme .

Scheme 9

Proposed Catalytic Cycle and Computed Reaction Coordinate

Direct acylation begins with BTM attack on allylic ammonium activated substrate (TS-II, ΔG‡ = 14.8 kcal·mol–1) to form tetrahedral intermediate III. Release of PNPO– (TS-IV, ΔG‡ = 12.0 kcal·mol–1) gives dication V. Indirect acylation through formation of the ammonium ketene III′ was ruled an unlikely reactive intermediate based on its unfavorable thermodynamics (ΔG‡ = 22.0 kcal·mol–1). The endergonicity of dication V (ΔG‡ = 3.3 kcal·mol–1) confirms the observed reversibility of catalyst acylation. Dication V is in equilibrium with ylide VII (ΔG = 1.6 kcal·mol–1) through deprotonation of the α-proton of V by PNPO– (ΔG‡ = 11.4 kcal·mol–1), also in agreement with the experimentally observed reversibility of the deprotonation step.[30,31]

NBO analyses reveal significant enolate character of ylide intermediate VII. Intermediate VII subsequently undergoes stereoselective and turnover-rate limiting [2,3]-rearrangement (TS-VIII-(23)-Major, ΔG‡ = 17.3 kcal·mol–1) to yield enantio- and diastereoenriched acyl ammonium product–catalyst complex IX. Catalyst turnover is found to be stepwise, begins with PNPO– attack (TS-X) and ends with catalyst and product release (TS-XII). The barrier for PNPO– attack as calculated from intermediate IX (ΔG‡ = 16.9 kcal·mol–1) indicates that, in the absence of HOBt, this step is highly competitive with rearrangement as turnover-rate limiting.

Effect of Counterions on the Theoretical KIE

The presence of counterions posed a challenge to the accuracy of DFT and significantly increased the complexity of the conformational search and the number of relevant structures to consider.[32] Almost all species present in the catalytic cycle prior to catalyst turnover bear a positive charge, with intermediate V being dicationic. Species indicated to include a counterion in Scheme were optimized with the explicit ion shown. Given the charged nature of the species present, the identification of the structures that compose the free energy span[33] resulted from considering all possible counterion coordination combinations for all conformations of each charged species in the catalytic cycle. This exhaustive process led to the identification of acyl substrate I and TS-VIII as the two states contributing to the free energy span. Computed barriers show the rearrangement step as the first irreversible step of the mechanism, thereby allowing kinetic isotopic fractionation to occur. The computed KIE then depends on the vibrational frequencies of I and TS-VIII, and could then be utilized to corroborate the computed thermodynamics and barriers of this free energy span.[34]

In a multistep reaction with highly charged and zwitterionic species, leveraging KIE prescribes a means to identify not only the structures that compose the free energy span[33] but also which ions coordinate, the specific binding site of the counterion,[35] and the conformation.[27b,36] We sought to identify the coordination state of TS-VIII by leveraging both the KIE and computed barriers for TS-VIII-F, i.e. bearing the 4-fluoro substituent used for KIE determination (Scheme ). Coordination to TS-VIII-F, formation of byproduct salt complexes, and conformations all affect the barrier in going from I-F to TS-VIII-F. Two possible counterions, PNPO– and Br–, were considered as TS counterions, while iPr2NH2+ was evaluated as a component of the possible remaining complexes (Figure ). No coordination to the TS (Figure , left) leaves H-bond complex PNPOH···PrNH as the lowest-energy remaining complex, giving an overall ΔG‡ = 18.3 kcal·mol–1. Bromide ion binding to the TS (TS-VIII-F-(2,3)···Br, Figure , middle) also leaves complex PNPOH···PrNH (ΔG‡ = 18.7 kcal·mol–1). PNPO– binding and the complexation of iPr2NH2+ and Br– gives an even higher barrier (TS-VIII-F-(2,3)···PNPO, ΔG‡ = 22.8 kcal·mol–1). With no counterion coordination to the TS (TS-VIII-F()-Major), the KIEcomp of 1.028 matches well with experiment (KIEexp = 1.031). Bromide complexation, which is 0.4 kcal·mol–1 higher, also matches fairly closely, giving a KIEcomp of 1.041.[38] PNPO– complexation leads to an erroneously large magnitude of rate difference between kH/kD, yielding a KIEcomp of 1.050.

Figure 6

Computed TSs, ions, complexes, and KIEs involving the 4-fluoro substituted substrate I-F. The computed KIE depends on the coordination state of substrate I-F and TS-VIII-F. The violet highlighted atom is the isotopic proton (H/D).[37] All energies in kcal·mol–1. Shaded gray lines represent forming/breaking bonds. Green lines represent C–H electrostatic interactions and hydrogen bonds.

Stereocontrol Model

The computed diastereomeric [2,3]-rearrangement TSs are shown in Figure . All TSs feature concerted C–C bond formation and ammonium N–C bond cleavage.[4,6] The four main elements that control the stereochemical outcome of the reaction are as follows (Figure ):

Figure 8

Computed model systems (all energies in kcal·mol–1 and distances in Å). (A) Preference for Z over E enolate. Enolate-like character indicated by bond orders (B. O.) estimated from the Wiberg bond indices (bottom left inset). (B) Effect of S···O interaction on acylated catalyst conformation. (C) With the acylated catalyst conformation held rigid (S···O), the BTM stereodirecting Ph sterically biases open enolate face. (D) Endo rearrangement is favored. In the major TS, this preference is reinforced by a π–cation interaction (bottom right inset).

E vs Z Configuration of the Enolate in Ylide VII

NBO analysis indicates that both ylide VII and TS-VIII(2,3)Major have significant enolate character.[39] Ylide VII displays a C–O bond order of 1.39 and a C–C bond order of 1.52 (Figure , bottom left inset), while TS-VIII(2,3)Major displays a C–O bond order of 1.54 and a C–C bond order of 1.21.[11] The computed bond order of 1.52 for ylide VII suggests partial C–C double bond character leading to distinct isomeric E and Z enolate configurations prior to rearrangement with the configurations set in place by the deprotonation step. The Z-configuration is heavily favored over the E, as shown in the model system /-35 where the Z is favored by >16 kcal·mol–1. All stable [2,3]-rearrangement transition structures feature the Z-enolate.

Anti vs syn Scatalyst to Osubstrate Orientation

In all the lowest energy conformations of the ylide-VII and rearrangement TS-VIII, the S–O relationship is syn. The syn distances (∼2.7–2.8 Å) are significantly below the sum of the van der Waals radii (3.4 Å), indicating close-contact S···O interactions (Figures and 8, orange lines).[40] Computed model systems show >4 kcal·mol–1 preference for the conformation which contains the 1,5-S···O interaction (/-36, Figure ). All [2,3]-rearrangement TSs that do not bear the S···O interaction are higher by >6 kcal·mol–1.[11] The conformational bias toward the S–O arrangement is proposed to result from nO to σ*C–S delocalization coupled with electrostatic attraction of the partially positive sulfur atom and partially negative oxygen atom.[41]

Figure 7

Stereodetermining [2,3]-rearrangement TSs.[34] All energies in kcal·mol–1, and distances in Å. Shaded gray lines represent forming/breaking bonds. Solid orange lines represent nonbonding S···O interactions. Dashed blue lines represent aromatic interactions.

Facial Selectivity of Rearrangement

The S···O interaction significantly rigidifies the ylide-VII structure, leaving conformational freedom only to the substrate cinnamyl group. With rearrangement possible from either face of the catalyst isothiourea plane, the facial selectivity is controlled by the catalyst Ph stereodirecting group. The most favorable [2,3]-TSs favor approach opposite to this group (TS-VIII-(2,3)-Major and TS-VIII-(2,3), Figure ). Approach on the same side as the stereodirecting Ph is disfavored by >6 kcal·mol–1 (TS-VIII-(2,3) and TS-VIII-(2,3)).

Endo vs exo [2,3]-TS

Rearrangement can occur either endo or exo with respect to the substrate C=O (Figure ). In the simple allyl model TS, the endo/exo preference is ∼1 kcal·mol–1. This preference is ∼2 kcal·mol–1 between TS-VIII-(2,3)-Major and TS-VIII-(2,3), and additional interactions contribute to the diastereoselectivity. In the major TS there is a π–cation interaction,[42] which is favored over the π-C–H interaction found in the minor.[43] Truncated fully optimized model systems probing the difference in energy between these interactions in the context of cationic BTM reveal a ∼1 kcal·mol–1 preference for 38 π-cation over 38 π-C–H (Figure , bottom right inset).[11] These two factors contribute to the computed 2 kcal·mol–1 preference for TS-VIII-(2,3)-Major over TS-VIII-(2,3), in good agreement with the experimental selectivity of 1.5 kcal·mol–1.

Conclusions

The experimental and computational investigation reported herein has provided mechanistic and stereochemical insight into the enantioselective isothiourea-catalyzed [2,3]-rearrangement of allylic ammonium ylides. Analysis by in situ19F NMR has allowed reaction profiles to be established and has identified an intermediate species.

Isotopic labeling of catalyst (15N) and substrate (13C) has confirmed the constitution of the catalytic intermediate as 29/IX by 13C NMR. Isotopic entrainment has shown 29/IX to be an irreversibly generated intermediate that is productive toward catalysis. A series of crossover experiments have provided detailed information regarding the reversibility of each individual step of the catalytic cycle. The turnover-rate limiting step of the process varies between product release and [2,3]-rearrangement, depending on substrate conversion. The effect of excess HOBt upon the reaction is to accelerate product release, thus generating a greater proportion of the free BTM catalyst (Figure ). This may then result in more effective interception of the background racemic reaction, and thus greater diversion onto the enantioselective pathway (Scheme ). Computational analysis has provided finer detail for the fundamental steps in the catalytic cycle as well as the key interactions that control the stereochemical outcome of the process. The insight gained into this process will have implications in a wider context, especially in the use of activated esters in Lewis base catalysis, which is currently under investigation in our laboratories.[44]