Wen-Jie
Wei
,
Kalman J.
Szabo
* and
Fahmi
Himo
*
Department of Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91 Stockholm, Sweden. E-mail: kalman.j.szabo@su.se; fahmi.himo@su.se
First published on 27th February 2025
The detailed reaction mechanism of diphenyl selenide-catalyzed sulfenofunctionalization of chiral α-CF3 allylboronic acids is investigated by means of density functional theory calculations. It is demonstrated that the reaction starts with the transfer of the SCF3 group from the (PhSO2)2NSCF3 reagent to the Ph2Se catalyst, a process that is shown to be assisted by the presence of Tf2NH acid. After a proton transfer step, the SCF3 group is transferred to the CC double bond of the substrate to generate a thiiranium ion. Concerted deborylative opening of the thiiranium ion yields then the final product. Several representative substrates are considered by the calculations, and the origins of the stereoselectivity of the reactions are analyzed by comparing the energies and geometries of the transition states leading to the different products.
Very recently, Szabó and co-workers presented an efficient method to introduce the SCF3 group by using (PhSO2)2NSCF3 reagent362 with α-CF3 allylboronic acid derivatives 1 to form a chiral alkenyl SCF3 compound 5 (Scheme 1).37 The reaction relies on the application of selenium-based Lewis base Ph2Se 3 as a catalyst in the presence of Tf2NH 4 as the activator. These reactions have a high degree of functional group tolerance and proceed with excellent stereo-, diastereo- and site-selectivity.
![]() | ||
Scheme 1 Se-catalyzed sulfenofunctionalization of allylboronic acid to form chiral allyl SCF3 compound investigated in the present study. |
The application of electrophilic sulfenofunctionalization in the presence of selenium catalysis has been documented by the pioneering studies of the groups of Denmark,38 Zhao,26 and others.39 However, the application of allylboronate substrates for trifluoromethylthiolation revealed a couple of new, interesting mechanistic aspects.37 An important feature is the excellent stereochemistry of the reaction, including chirality transfer and high E-selectivity for the formation of the new double bond. The studies also pointed out the benefits of using allylboronic acid type of substrates. Replacement of the B(OH)2 group of 1 with Bpin leads to a significant decrease in the yield of the product, indicating that allyl-Bpin species have a substantially lower reactivity than allylboronic acids.
In the present work, we have performed density functional theory (DFT) calculations to elucidate the mechanism for the formation of allyl SCF3 compounds via the selenium-catalyzed sulfenofunctionalization of allylboronic acids. We consider the reactions of several representative substrates (1a–1d in Scheme 1) and discuss the origins of selectivity for each of them.
To evaluate the effect of performing the geometry optimizations with the smaller basis set and in the gas phase, the geometries of the first step of the reaction (2 + 3 + 4 → Int1, see below) were re-optimized twice, in implicit DCM solvent and with the larger basis set. The calculations showed that the geometries and energies were not affected significantly (see ESI†).
Intrinsic reaction coordinate (IRC) calculations were performed on all transition states to confirm the nature of the connecting intermediates. A thorough manual conformation search was carried out on all stationary points to ensure that the structures with the lowest energy were located.
![]() | ||
Fig. 1 Optimized structures of the transition states and intermediates of TS1, Int1, and TS2. Selected bond distances are indicated in Å. |
Next step is a proton transfer from 4 to (PhSO2)2N− anion 6viaTS2, generating (PhSO2)2NH 8 and Tf2N− anion 9 (Scheme 2b). The energy barrier of this step is only 0.6 kcal mol−1 relative to Int1.
We also calculated the alternative pathway with the reversed step order, i.e. in which a proton is first transferred from acid 4 to (PhSO2)2NSCF3, followed by the transfer of the SCF3 group from the generated (PhSO2)2NHSCF3 cation to the Ph2Se catalyst 3. However, the energy barrier of this scenario was found to be very high, 52.9 kcal mol−1 (see ESI†).
In the following step of the mechanism, the SCF3 group of the SePh2SCF3 cation 7 is transferred to the CC double bond of the substrate to generate a thiiranium ion (Scheme 2c). Two competing stereoselective pathways are possible and were investigated, in which the SCF3 group is transferred either to the Re-face of the substrate viaTS3a(3R) to form the (R)-configuration Int2a(3R), or to the Si-face viaTS3a(3S) to generate the (S)-configuration Int2a(3S). The optimized structures of the transition states and intermediates are shown in Fig. 3. The energies of TS3a(3R) and TS3a(3S) are calculated to be 18.4 and 18.6 kcal mol−1, respectively, relative to the previous intermediate, and the resulting Int2a(3R) and Int2a(3S) are 10.6 and 11.2 kcal mol−1 higher, respectively.
![]() | ||
Fig. 3 Optimized structures of the transition states and intermediates for the generation of the thiiranium ion Int2a. Selected bond distances are indicated in Å. |
For comparison, we also investigated the uncatalyzed reaction, in which the SCF3 group is transferred directly from reagent 2 to substrate 1a, generating Int2a without the participation of catalyst 3 (see ESI†). The energy barrier was calculated to be 29.9 kcal mol−1 for both the 3R- and 3S-pathways, which is significantly higher than the case with catalyst 3. Accordingly, using PhSePh catalyst (3) increases the reactivity of 2. In the positively charged 7 the electrophilicity of SCF3 is substantially increased compared to 2. In addition, the cleavage of the weak Se–S bond in 7 is also easier than the cleavage of the N–S bond of 2. The stability of 7 is poor under ambient conditions, and therefore 2 is converted to 7in situ (in the presence of 4) under the applied reaction conditions.
The final step of the mechanism (Scheme 2d) involves the Tf2N− anion 9 performing a nucleophilic attack on the boronate group of Int2a(3R) or Int2a(3S), triggering the concerted deborylative opening of the thiiranium ion viaTS4 to yield the four possible forms of the final product. The (3R)-configured products are formed from Int2a(3R) and can result in either the E-configuration through the anti-elimination pathway viaTS4anti(3R), or alternatively, rotation of the Cα–Cβ bond leads to the syn-elimination pathway viaTS4syn(3R), resulting in the Z-configuration. Similarly, the E-(3S)- or Z-(3S)-configured products can be achieved viaTS4syn(3S) or TS4anti(3S), respectively. The optimized structures of these transition states are shown in Fig. 4, while the optimized structures of the products are given in the ESI.†
![]() | ||
Fig. 4 Optimized structures of the transition states for the step of deborylative elimination in model substrate 1a. Selected bond distances are indicated in Å. |
Although the chiral center is formed at TS3, the calculations show that the energy barriers for the eliminations of the boronate group viaTS4 are higher and irreversible, which means that the latter step is the stereoselectivity-determining step of the reaction.
By comparing the energy profiles of the above four pathways (Fig. 2), the calculations of the model substrate 1a show that the energy barriers of anti-eliminations are considerably lower than the syn-eliminations, 21.7 kcal mol−1 in TS4anti(3R)vs. 27.8 kcal mol−1 in TS4syn(3R), and 22.0 kcal mol−1 in TS4anti(3S)vs. 29.1 kcal mol−1 in TS4syn(3S). Inspection of the optimized structures in Fig. 4 shows that the reason for this energy difference is mainly the steric repulsion between the SCF3 group and the leaving boronate group, as these two moieties point toward each other in the syn-elimination, TS4syn(3R) and TS4syn(3S), while in the anti-elimination, TS4anti(3R) and TS4anti(3S), they point away from each other.
The difference in energy between TS4anti(3R) and TS4anti(3S), which lead to the E-(3R)-5a and Z-(3S)-5a products, respectively, is calculated to be only 0.3 kcal mol−1 in favor of the former. The calculations show thus that already the model substrate 1a, with the small methyl substituent, qualitatively reproduces the experimental selectivity trend, in that the formation of the E-(3R)-configured product is associated with the lowest-energy pathway, albeit with a small energy. However, as will be shown below, the calculations using the experimentally employed substrates, with bulkier substituents, yield a more quantitative agreement with the experiments.
To summarize, the mechanism proposed on the basis of the current calculations is given in Scheme 3. The obtained overall energy profile (Fig. 2) indicates that the first step, i.e. the acid-assisted transfer of the SCF3 group viaTS1, is the rate-determining step (RDS) for the model substrate, with a barrier of 24.2 kcal mol−1. However, the final step, i.e. deborylative opening of the thiiranium ring viaTS4, has an overall barrier of 22.2 kcal mol−1, which is quite close in energy, and it is therefore not possible to determine confidently the nature of the RDS based only on the calculations. In particular, various substituents on the substrate may lead to significant changes in the energy of the final step (see below).
Here, it is interesting to mention two previous mechanistic studies on sulfenofunctionalizations of alkenes catalyzed by selenides, where DFT calculations were employed to investigate various aspects of the reactions. However, none of these studies involved a deborylation step, which is a novel aspect of the present study. Denmark and co-workers analyzed the geometries and energies of the transition states for the thiiranium ion formation step, which was assumed to be the enantio-determining step of the reaction,46 while Zhao and co-workers investigated the mechanism of selenide-catalyzed trifluoromethylthiolation of gem-diaryl tethered alkenes to synthesize trifluoromethylthiolated tetrahydronaphthalenes.27
In addition to the results discussed above, we have also considered some other mechanistic alternatives that turned out to have higher energy barriers. As seen from Fig. 2 and 4, catalyst 3 does not participate in the deborylative elimination step in TS4. We have considered whether it can assist this step, but the energy barriers for this scenario were found to be higher (see ESI†). We also considered whether the (PhSO2)2NH species 8 could act as the nucleophile to attack the boronate group of Int2a, but the calculated energy barriers for this pathway were calculated to be very high as compared to when the anionic Tf2N−9 is the nucleophile (see ESI†). Finally, the experiments reveal that replacing the B(OH)2 group of the substrate with Bpin significantly decreases the yield of the product.37 Consistently with this result, the calculations show that the barrier for the case of Bpin is 3.7 kcal mol−1 higher than the case of the B(OH)2 group (see ESI†).
For substrate 1b, with the benzyl substituent, the mechanism was calculated to be very similar to that of the model substrate 1a shown in Scheme 3 (see calculated energy profile in Fig. 5). One small difference is that the formation of the Z-(3R)-configured product through the syn-elimination was found to occur in a stepwise manner (see ESI† for detailed results). The calculations show that the overall barrier for substrate 1b is ca. 2 kcal mol−1 lower than for 1a, and very importantly, the extent of the stereo-differentiation is well-reproduced.
Similarly to substrate 1a, the barriers of anti-eliminations for 1b are considerably lower than the syn-eliminations due to steric repulsion between the SCF3 group and the leaving boronate group. In addition, the pathway leading to the E-(3R)-5b product is now 3.7 kcal mol−1 lower than that leading to the Z-(3S)-5b product (Fig. 5), due to a steric repulsion between the benzyl substituent of the substrate and the SCF3 group (see ESI†). The calculated energy difference is in good agreement with the experimentally observed ee of 93% in favor of the E-(3R) product.
For substrate 1c with the tBu substituent (see ESI†), the situation is very similar to substrate 1b, with both syn-eliminations found to take place in a stepwise manner. The overall barrier was calculated to be ca. 3 kcal mol−1 higher than for 1b, and the selectivity is determined by the same factors, with an energy difference of 4.2 kcal mol−1, in good agreement with the experimental outcome of 98% ee.
Substrate 1d, with the phthalimide substituent, represents an interesting case, because the carbonyl group of 1d may perform a nucleophilic attack at the Cβ atom of the thiiranium ion through an intramolecular mechanism,25,47,48 leading to the opening of the thiiranium ion and yielding a six-membered ring intermediate Int3d′ (see Fig. 6). The calculations show that the energy barrier for this competing intramolecular nucleophilic attack viaTS4d′ is much lower than for the intermolecular reaction of the thiiranium ion with the Tf2N− anion 9viaTS4d, which was found for the other substrates. However, from Int3d′, the barriers for the following steps, which would be the nucleophilic attack of Tf2N− and the deborylative opening of the six-membered ring, were found to be higher in energy compared to the intermolecular reaction (see ESI†), indicating that the intramolecular nucleophilic attack is a reversible process. Thus, formation of Int3d′ can be regarded as an unproductive dead-end for the deborylative trifluoromethylthiolation process.
The calculations show thus that the reaction of substrate 1d also follows the mechanism of the model substrate 1a. However, as seen from Fig. 6, the energy of Int3d′ is the lowest point on the energy profile before TS4d, which means that the overall barrier should be calculated relative to Int3d′, resulting in a slightly higher barrier as compared to the other substrates (24.9 kcal mol−1 compared to 22.2, 20.3, and 23.2 kcal mol−1, for 1a, 1b and 1c, respectively). Importantly, the stereoselectivity is reproduced also for substrate 1d, with a selectivity-determining energy difference of 2.5 kcal mol−1, in good agreement with the 97% ee observed experimentally. The origins of the selectivity are found to be the same as for the other substrates.
The reaction starts with the generation of the catalytically active SePh2SCF3 cation, a process that takes place is in a stepwise manner, with a transfer of the SCF3 group from the (PhSO2)2NSCF3 reagent 2 to the Ph2Se catalyst 3, followed by a proton transfer from acid Tf2NH 4 to the formed (PhSO2)2N− anion 6. Interestingly, the Tf2NH acid stabilizes the negative charge that develops on the nitrogen anion of the (PhSO2)2N− species 6 through a hydrogen bond interaction, lowering thus the barrier for the SCF3 transfer.
Next, the SCF3 group of the formed SePh2SCF3 cation is transferred to the CC bond of the α-CF3 allylboronic acid, generating the thiiranium ion species Int2. Two different pathways are possible, depending on whether the SCF3 group is transferred to the Si or Re face of the C
C bond, which eventually lead to the S- or R-configurations of the product, respectively.
Finally, the Tf2N− anion 9 performs a nucleophilic attack at the boronate group of Int2, triggering the opening of the thiiranium ion and the leaving of the boronate group. Rotation around the Cα–Cβ bond of Int2 leads to either the syn- or anti-elimination, generating the E- or Z- configurations of the product. This step constitutes the selectivity-determining step of the reaction, and the calculations show a clear preference for formation of the E-form, in excellent agreement with the high E-selectivity reported in the experimental studies.37 The enantioselectivity is also very well reproduced by the calculations, and analysis of the transition state structures shows that the selectivity is mainly controlled by steric repulsions between the SCF3 group and both the leaving boronate group and the substituent of the α-CF3 allylboronic acid.
The mechanistic insights provided by the current study will be valuable for the development of new regio-, stereo- and diastero-selective selenide-catalyzed, deborylative electrophilic sulfenofunctionalization reactions.
Footnote |
† Electronic supplementary information (ESI) available: Additional computational results discussed in the text, absolute energies and energy corrections, and Cartesian coordinates. See DOI: https://doi.org/10.1039/d4qo02170c |
This journal is © the Partner Organisations 2025 |