Chunhui Liua,
Peilin Hana,
Zhizhong Xiec,
Zhihong Xu*a and
Donghui Wei*b
aSchool of Chemistry and Chemical Engineering, Xuchang University of China, Xuchang, Henan Province 461000, P. R. China
bThe College of Chemistry and Molecular Engineering, Center of Computational Chemistry, Zhengzhou University of China, Zhengzhou, Henan Province 450001, P. R. China
cDepartment of Chemistry, School of Chemistry, Chemical Engineering and Life Sciences, Wuhan University of Technology, Wuhan, Hubei Province 430070, P. R. China
First published on 4th December 2018
The competing mechanisms of Ag(I)-catalyzed chemoselective addition reactions of amino alcohols and electron-deficient olefins leading to the O-adduct or N-adduct products were systematically studied with density functional theory methods. Calculations indicate that the AgHMDS/dppe versus AgOAc/dppe catalytic systems can play different roles and thereby generate two different products. The AgHMDS/dppe system works as a Brønsted base to deprotonate the amino alcohol OH to form the Ag–O bond, which leads to formation of the O-adduct. In contrast, the AgOAc/dppe system mainly acts as a Lewis acid to coordinate with O and N atoms of the amino alcohol, but it cannot act as the Brønsted base to further activate the OH group because of its weaker basicity. Therefore, the AgOAc/dppe catalyzed reaction has a mechanism that is similar to the non-catalyzed reaction, and generates the same N-adduct. The obtained insights will be important for rational design of the various kinds of cooperatively catalyzed chemoselective addition reactions, including the use of the less nucleophilic hydroxyl groups of unprotected amino alcohols.
Recently, a number of studies were reported to have solved these questions.5–12 In 2006, Gunnoe and coauthors8,9 reported the determination method of reaction rates of O vs. N nucleophiles in the Cu(I) catalyzed Michael addition. Subsequently, the Ohshima group10–12 reported that a series of soft Lewis acids (Cu(I) and Ag(I) salt) catalyzed the chemoselective conjugate addition reactions of amino alcohols to electron-deficient olefins (an oxa-Michael reaction), which enabled the chemoselective deprotonation of the hydroxyl group (OH) of the amino alcohol. As shown in Scheme 1, the Ohshima group reported the first example of the Lewis acid and base cooperatively catalyzed conjugate addition reaction of amino alcohols to α,β-unsaturated nitriles in 2014.10 Generally, the nucleophilicity of the NH2 group is stronger than that of the OH group in amino alcohols, and N-addition should be the preferred process. Interestingly, the cooperative nature of the catalyst (Cu(MeSal)/dppe and LiHMDS) enabled chemoselective activation of alcohols over amines, leading to the O-addition reaction. Then the Ohshima group reported the addition reaction of amino alcohol to α,β-unsaturated sulfonyl derivatives,11 as shown in Scheme 2, gave varied reaction products according to different reaction conditions. The AgHMDS/dppe system and the cooperative catalyst system (Ag(OAc)/dppe and KHMDS) enabled chemoselective activation of alcohols over amines to induce the O-addition reaction, but the AgOAc/dppe catalyst cannot reverse the innate activity order of OH and NH2 groups of amino alcohols. In 2017, the Ohshima group reported the chemoselective conjugate addition of an amino alcohol to an α,β-unsaturated ester,12 as shown in Scheme 3, in which the cooperative catalyst achieved chemoselective addition of the hydroxyl group over the amino group, and favored chemoselective conjugate addition over transesterification. These reactions not only contribute to atom and step economy,13 but also reverse the innate activity order of OH and NH2 groups of amino alcohols.
Scheme 2 The different O- and N-additions of amino alcohols to α,β-unsaturated sulfonyl derivatives controlled by different catalysts. |
Scheme 3 The different O- and N-additions of amino alcohols to α,β-unsaturated esters controlled by different catalysts. |
The reactions discussed in Schemes 1–3 are very convenient and unique, therefore, the detailed reaction mechanisms should be worthy to disclose, and the origin of switchable chemoselectivities is urgent to confirm and explore deeply in theory. Prompted by the interesting experimental observations and our interest in catalysis mechanisms,14 we selected and investigated the addition reactions of amino alcohols and olefins based on the different catalytic systems in Scheme 2 as model reactions. In the present study, we sought to provide exact pictures of the detailed reaction mechanisms, explore the roles of the different catalysts, and find the origin of the chemoselectivity. All the calculations were carried out using density functional theory (DFT), which has been confirmed to be a powerful tool for investigating enzyme-catalyzed,15 organocatalytic,16 and transition metal-catalyzed17 reaction mechanisms and other theoretical studies.18
Fig. 1 Computed Gibbs free-energy profiles (298 K) for chemoselective addition reactions without catalyst. |
Fig. 2 Computed Gibbs free-energy profiles (253 K) for chemoselective addition reactions catalyzed by 3′ system (distances in Å). |
Fig. 3 The optimized 3-D structures of the transition states in Fig. 2 and 4. Some H atoms are omitted for clarity (distances in Å). |
As depicted in Fig. 2, phenyl vinyl sulfone 2 approaches ADD2-O to undergo the O- or N-addition process via the transition states TS2-O or TS2-N with the relative free energy of −2.0 or 18.1 kcal mol−1, respectively. For the O-addition process, the C2 atom of 2 undergoes nucleophilic attack by the O atom of the Ag–O group via the transition state TS2-O with a relative free-energy difference of −2.0 kcal mol−1. In TS2-O, the O–C bond distance is 1.80 Å, and HHMDS locates close to the carbanion of the CC bond. TS2-O evolves towards ADD4-O in which the C–O bond has been formed with a carbanion based on the CC bond. HHMDS then works as a Brønsted acid to return a proton to the carbanion of the CC bond via transition state TS3-O with a relative free-energy difference of −0.6 kcal mol−1. ADD5-O is 6.6 kcal mol−1 below the reactant energy level. Separation of the O-adduct and regeneration of catalyst 3′ is favorable on the free-energy scale; the overall exergonicity of the whole process is 19.3 kcal mol−1.
For the N-addition process, the CC bond of 2 undergoes nucleophilic attack by the NH2 group of ADD2-O via a four-membered ring transition state TS2-N with a relative free-energy difference of 18.1 kcal mol−1. In TS2-N, the HHMDS group is weakly bonded to the substrate with the hydrogen bond (1.91 Å). The N-addition process leads to ADD4-N, and then HHMDS works as a Brønsted acid to return a proton to the O atom of the substrate via four-membered ring TS3-N with an energy barrier of 9.1 kcal mol−1. In TS3-N, a proton is transferred from the N atom of HHMDS to the O atom with a N–H distance of 1.29 Å and an O–H distance of 1.23 Å, while breaking of the Ag–O bond (2.37 Å) is coupled with formation of Ag–N (2.40 Å). TS3-N evolves towards ADD5-N, in which the N-adduct is weakly bound to catalyst 3′ through a hydrogen bond (1.99 Å).
As shown in Fig. 2, catalyst 3′ works as a Brønsted base to capture a proton from the OH group and form the Ag–O bond to attain ADD2-O before phenyl vinyl sulfone 2 approaches ADD2-O to undergo the O- or N-addition process. We then considered whether the Ag–O group can attack the CC bond of 2 together to form the four-membered ring, which is a similar process to the N-addition process via TS2-N. This process is depicted in Fig. 4, and the optimized 3-D structures of the transition states are shown in Fig. 3. As shown in Fig. 4, the Ag–O bond adds to the CC bond of 2 to form the four-membered ring transition state TS2′-O. In TS2′-O, the HHMDS group is weakly bound to the substrate via the hydrogen bond (2.05 Å). This new process leads to ADD4′-O, and then HHMDS returns a proton to the C atom of the substrate through four-membered ring TS3′-O. In TS3′-O, a proton is transferred from the N atom of HHMDS to the C atom of the substrate (N–H, 1.29 Å; C–H, 1.48 Å), and the Ag–C bond is broken to form the Ag–N bond. This process is not favored because of the high energy barrier of 21.9 kcal mol−1 (energy difference between ADD4′-O and TS3′-O). Therefore, for the O-addition process, the four-membered ring TS2′-O formed by the Ag–O group direct attack of the CC bond of 2 is unsuitable.
Fig. 4 Computed Gibbs free-energy profiles (253 K) for the possible reaction pathways via TS2′-O and TS3′-O catalyzed by 3′ (distances in Å). |
If catalyst 3′ only acts as a Lewis acid and does not work as a Brønsted base, the mechanism should proceed in only one step as shown in Fig. S1 of ESI.† Initially, reactant 1 approaches catalyst 3′ to form ADD1′, and a hydrogen bond forms between the N atom of HMDS− and the OH group of 1 (1.90 Å). Afterwards, the approach of 2 to ADD1′ results in N-addition or O-addition through the respective transition states TS1′-N or TS1′-O. TS1′-N and TS1′-O are four-membered ring transition states with relative free energies of 17.4 and 33.3 kcal mol−1, respectively. In TS1′-N depicted in Fig. S1,† the NH2 group adds to the vinyl group to form the N-addition product (N-adduct), while for TS1′-O, the OH group adds to the vinyl group to form the O-addition product (O-adduct). The energy difference indicates that the Lewis acid-catalyzed N-addition pathway is more energetically favorable than the Lewis acid-catalyzed O-addition pathway, but this is not consistent with the experimental observations depicted in Scheme 2.11 Therefore, the mechanism via TS1′-N and TS1′-O for which catalyst 3′ only works as a Lewis acid cannot be the most favorable pathway for this system.
As shown in Fig. 2, the O-addition reaction pathway associated with transition states TS1-O, TS2-O, and TS3-O should be the most favorable pathway with an energy barrier of 9.5 kcal mol−1. In this case, the AgHMDS/dppe catalyst not only acts as a Lewis acid to coordinate with the O and N atom of the substrate, but also works as a Brønsted base to capture a proton from the OH group to form the Ag–O bond and as a Brønsted acid to return a proton to recycle itself. Comparing the structure of TS2-O with those of TS1′-N and TS2-N, the O atom of the Ag–O group is more suitable than the N atom of the NH2 group for nucleophilic addition to the CC bond of 2. This is probably because the O atom of Ag–O has a stronger nucleophilicity than the N atom of NH2 groups, which is discussed in the Parr function analysis. This factor is responsible for the chemoselectivity of the addition reaction and O-adduct formation.
In the overall Gibbs free-energy profiles (Fig. 2, 4, and S1†), the second step (conjugate addition process) determines the chemoselectivity of O- and N-addition through TS2-O, TS2-N, and TS1′-N. The energy of TS2-O is much lower than those of TS2-N and TS1′-N, which indicates that the chemoselectivity of O-addition is preferred in this system, and is consistent with the experimental observations shown in Scheme 2.11 To compare the effects of the real and model AgHMDS/dppe system (i.e., 3 and 3′) on the chemoselectivity of the reaction, we additionally constructed and optimized TS2-O-real, TS2-N-real, and TS1′-N-real with the actual dppe ligand by employing catalyst 3. As shown in Fig. 5, the energy of transition state TS2-O-real associated with O-addition is still lower than those of TS2-N-real and TS1′-N-real associated with the N-addition. We also reproduced Fig. 2 with the actual dppe ligand by using catalyst 3 (see the Fig. S2 of ESI†). The computed Gibbs free-energy profiles and energy barriers showed the same trend and conclusions as the calculated results of the model systems. This indicates that the selected models and the computed results should be suitable and reliable.
Fig. 5 The optimized 3-D structures and relative free energies (253 K, kcal mol−1) of TS2-O-real, TS2-N-real and TS1′-N-real. Some H atoms are omitted for clarity (distances in Å). |
In the favored O-addition reaction pathway, the deprotonation process is the rate-determining step. We computed the TS1-O–Cu-HMDS and TS1-O–Cu-Bu transition states with CuHMDS/dppe and CuBu/dppe as catalysts, in which the phenyl groups of dppe were replaced by methyl groups. As shown in Fig. 2, the free energies related to the separate reactants of TS1-O–Cu-HMDS and TS1-O–Cu-Bu are higher than that of TS1-O, indicating that the rate of Ag(I)-catalyzed addition would be faster than that of the Cu(I)-catalyzed addition. This is consistent with the experimental observations described in Ohshima's report.11 To test the impact of different computational methods on the calculated results, we additionally computed the rate-determining step by using different computational methods. As summarized in Table S1 (ESI),† the free energies calculated by the different methods are very close and have small differences, which indicate that the calculated results should be reliable by using the selected computational method.
Fig. 6 The energy scans of the deprotonation processes along with the distance changes and 3-D structures. Some H atoms are omitted for clarity (distances in Å). |
The detailed mechanisms of the addition reaction catalyzed by 4′ and the corresponding energy profiles are shown in Fig. 7. Initially, L-phenylalaninol 1 approaches catalyst 4′ to form ADD1-OAc, and two hydrogen bonds are formed between the oxygen atom of OAc− and the OH and NH2 group with distances of 1.79 Å and 1.92 Å. ADD1-OAc sits only 0.03 kcal mol−1 above the reactants in energy, and afterwards the approach of 2 to ADD1-OAc results in two separate conjugate addition processes (O-addition and N-addition) through transition states TS1-N–OAc and TS1-O–OAc. Both TS1-N–OAc and TS1-O–OAc are four-membered ring transition states. In TS1-N–OAc depicted in Fig. 7, the NH2 group adds to the vinyl group to form the N-addition product (N-adduct) via nucleophilic attack of the N atom on the C2 atom of 2 and electrophilic attack of the H atom on the C1 atom of 2. Similarly, in TS1-O–OAc (Fig. 7), the OH group adds to the vinyl group to form the O-addition product (O-adduct) via addition of the O atom to the C2 atom of 2 while the H atom adds to the C1 atom of 2. The free energies of transition states TS1-N–OAc and TS1-O–OAc are 17.0 and 33.8 kcal mol−1, respectively, relative the reactants. This energy difference indicates that the N-addition process is preferred, while the O-addition process is difficult because of the much higher free-energy barrier of 33.8 kcal mol−1. The N-addition process mainly leads to ADD2-N–OAc, in which the N-adduct is weakly bound to catalyst 4′ through hydrogen bonds (1.71 Å and 2.38 Å) and Ag–O through a coordination bond (2.49 Å). ADD2-N–OAc sits 22.5 kcal mol−1 below the reactants in energy terms. Separation of the main N-adduct product and regeneration of catalyst 4′ can proceed smoothly with an overall exergonicity of 29.8 kcal mol−1.
Fig. 7 Computed free-energy profiles (253 K) and optimized 3-D structures for the chemoselective addition reactions catalyzed by 4′. Some H atoms are omitted for clarity (distances in Å). |
Comparing the structure of TS1-N–OAc with that of TS1-N-non depicted in Fig. 1, the Ag atom is coordinated with O and N atoms of the substrate and two hydrogen bonds are formed (1.67 Å and 1.97 Å) between the OAc− group and the OH and NH2 groups in TS1-N–OAc. This interaction would cause the energy barrier to be lower than that via TS1-N-non. Therefore, the rate of reaction catalyzed by the AgOAc/dppe system is higher than that without the presence of the catalyst. According to the above discussion, the transition states TS1-N–OAc and TS1-O–OAc would be the key to determining the chemoselectivity of the preferred N-adduct.
In addition, we found that the chemoselectivity has the same trend by constructing and computing the key transition states TS1-N–OAc-real and TS1-O–OAc-real catalyzed by AgOAc/dppe 4. As shown in Fig. 8, transition state TS1-O–OAc-real associated with the O-addition process is higher than transition state TS1-N–OAc-real associated with the N-addition process by 15.4 kcal mol−1. This is close to the energy difference of 16.8 kcal mol−1 between TS1-N–OAc and TS1-O–OAc depicted in Fig. 7. Therefore, the N-addition process still has the energy advantage for the AgOAc/dppe 4 system, and the chemoselectivity of the N-addition product is consistent with the experimental observations shown in Scheme 2.11
Fig. 8 The optimized 3-D structures and relative free energies (253 K, kcal mol−1) of TS1-N–OAc-real and TS1-O–OAc-real. Some H atoms are omitted for clarity (distances in Å). |
As discussed above, the AgOAc/dppe system only acts as a Lewis acid to form the coordination bonds with the O and N atoms of the substrate during the addition, but does not work as a Brønsted base to capture a proton from the OH and NH2 group. When compared with the AgHMDS/dppe system, AgOAc/dppe does not enable the deprotonation of the OH group to form the Ag–O bond, and cannot change the innate reactivity of the OH and NH2 groups. This means that the origin of chemoselectivity in the addition reaction is the same as the non-catalyst reaction depicted in Fig. 1, which leads to the N-adduct.
In brief, we can conclude that the AgOAc/dppe only works as Lewis acid, and the strong base (such as KHMDS, NaHMDS, LiHMDS) works as Brønsted base. Therefore, the whole system should be a cooperatively catalytic system in the favorable O-addition pathway. Obviously, both of the Ag(I) salt and the strong base additives should be necessary, which is in agreement with the experimental observation that the O-adduct is the main product with the presence of AgOAc/dppe and the strong base (such as KHMDS, LiHMDS, NaHMDS).11
It is noteworthy that reactant 1, intermediate ADD2-O (Fig. 2), and ADD1-OAc (Fig. 7) all have two reactive nucleophilic sites: the O atom of OH (denoted as 1-O) and the N atom of NH2 (denoted as 1-N) in reactant 1, the O atom of Ag–O (denoted as ADD2-O–O) and the N atom of NH2 (denoted as ADD2-O–N) in ADD2-O, and the O atom of OH (denoted as ADD1-OAc-O) and the N atom of NH2 (denoted as ADD1-OAc-N) in ADD1-OAc. For these species, the question of relative site reactivity was addressed by a simple calculation of Parr functions (Pk+ and Pk−).
As summarized in Table 1, the order of nucleophilic Parr function (Pk−) in reactant 1 is 1-N (0.0359) > 1-O (0.0092). This assessment reflects how the innate reactivity of the NH2 group is stronger than the OH group of the amino alcohol, and why the N-addition process is preferred in the absence of catalyst. The order of nucleophilic Parr function (Pk−) in ADD2-O is ADD2-O–O (0.5570) > ADD2-O–N (0.1034). This comparison confirms that the nucleophilic reactivity of the O atom of Ag–O is stronger than the N atom of NH2 in ADD2-O, and helps to explain why the O-addition process is preferred with the 3′ catalyst system. For ADD1-OAc, the order of nucleophilic Parr function (Pk−) is ADD1-OAc-N (0.0014) > ADD1-OAc-O (0.0013), indicating that the nucleophilic reactivity order of NH2 and OH does not change with the 4′ catalyst system, and that the N-addition process is still preferred. These results suggest that we can successfully predict the most nucleophilic reactive site by the Parr function analysis using a simple single-point calculation.
Parr function | 1-O | 1-N | ADD2-O–O | ADD2-O–N | ADD1-OAc-O | ADD1-OAc-N |
---|---|---|---|---|---|---|
Pk+ | −0.0015 | −0.0061 | −0.0099 | 0.0237 | 0.0011 | 0.0235 |
Pk− | 0.0092 | 0.0359 | 0.5570 | 0.1034 | 0.0013 | 0.0014 |
However, the calculated results suggested that AgOAc/dppe catalyst can only act as the Lewis acid to coordinate with O and N atoms of the substrate, but cannot work as the Brønsted base to enable deprotonation of the OH group to form the Ag–O bond. This means that the reaction catalyzed by Lewis acid AgOAc/dppe has the same mechanism and chemoselectivity as the uncatalyzed reaction via the direct four-membered ring transition state. Through this original research of Ag(I)-catalyzed chemoselective addition reactions, we reasonably expect that the obtained insights will be useful for rational design of more efficient addition reactions of olefins with OH and NH2 groups with special selectivities.
Footnote |
† Electronic supplementary information (ESI) available: Table S1, Fig. S1–S3, Cartesian coordinates, absolute energies and imaginary frequency of the computed structures. See DOI: 10.1039/c8ra09065c |
This journal is © The Royal Society of Chemistry 2018 |