Soumi
Tribedi
a,
Christopher M.
Hadad
b and
Raghavan B.
Sunoj
*a
aDepartment of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India. E-mail: sunoj@chem.iitb.ac.in
bDepartment of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
First published on 25th June 2018
Asymmetric catalysis using two chiral catalysts in combination using one-pot reaction conditions is in its initial stages of development and understanding. We employ density functional theory (SMD(toluene)/M06/6-31G**,SDD(Ir)) computations to shed light on the action of chiral phosphoric acid and a chiral Cp*Ir(diamine) in stereoinduction in an asymmetric amination reaction of an alcohol. First, the protonation of the Ir–diamine complex by the phosphoric acid forms an ion-pair of the active catalytic dyad. Both chiral catalysts are involved throughout the catalytic cycle, thus constituting an important example of true cooperative catalysis. A borrowing hydrogen mechanism operates, wherein the phosphate abstracts the hydroxyl proton of the alcohol while the electrophilic Ir(III) simultaneously extracts the α-hydrogen to form a [Ir]–H species. The ketone thus derived from the alcohol through dehydrogenation condenses with aniline to form an imine. In the diastereocontrolling transition state, the hydride adds to the activated iminium, held in position in the chiral pocket of the catalytic dyad through a network of noncovalent interactions (C–H⋯π, N–H⋯O and C–H⋯O). The enantioselectivity in this DYKAT process is identified as taking place at an earlier stage of the catalytic cycle prior to the diastereo-determining transition state.
The mechanism of such Brønsted acid catalysis is expected to become increasingly more intricate when it is used in conjunction with other chiral transition metal catalysts.1b Phosphoric acid and transition metal catalysts when used together under one-pot reaction conditions can present different mechanistic scenarios wherein the former can act as (a) a ligand coordinated to the metal, (b) a Brønsted acid hydrogen bonded to a basic center of the transition metal complex, or (c) a counteranion to the in situ formed positively charged intermediate.4 Deeper mechanistic understanding is required to ascertain the role of chiral Brønsted acids in such forms of catalytic reactions.
In a recent demonstration of dual catalysis, Zhao and coworkers have employed the borrowing hydrogen method for the asymmetric amination of alcohols. Importantly, both the organocatalyst (phosphoric acid) and the transition metal catalyst (Cp*Ir(diamine)) used as the catalytic dyad are chiral.5 The borrowing hydrogen method is an encouraging strategy to directly engage the readily available and low-cost alcohols in C–C bond formation reactions.6 The method can help oxidize an alcohol to a carbonyl compound, which in turn, can open up a gamut of transformations such as condensation with amines. The resulting imines can then be reduced to respective amines by the transfer of the “borrowed” hydrogen from the catalyst. Chiral amines can also be synthesized through the transfer hydrogenation strategy using small molecules such as alcohols and formic acids or Hantzsch esters as a sacrificial species.7 However, the borrowing hydrogen method offers a major advantage over transfer hydrogenation in that it does not require a sacrificial hydrogen donor; rather, the source of hydrogen is the substrate alcohol itself. While the idea showed early promise,8 the enantioselective variants of such reactions became available only recently by employing catalyst combinations involving compatible organocatalysts and transition metal catalysts.5 Zhao and coworkers reported an elegant asymmetric dual chiral catalytic (ADCC) method for the stereoconvergent amination between four isomers of a chiral alcohol and aryl amine (Scheme 1).9 The action of two chiral catalysts is not a well-understood concept owing to the scarcity of reports focusing on the mechanism and origin of stereoinduction in dual chiral catalytic systems.10
Scheme 1 A representative example of the asymmetric amination of an alcohol using two chiral catalysts.9 |
The fascinating feature in this reaction is that although all of the four stereoisomers of the parent alcohol (a pair of enantiomers and a pair of diastereomers) are involved in the reaction as a racemic mixture, only one stereoisomer of the product amine is obtained in excellent enantio- and diastereo-selectivity. The yield of the major diastereomer is more than 80%, suggesting that all of the stereoisomers of the reactant convert to the product. In this article, we wish to present the first molecular insights of stereoinduction in the asymmetric amination of a racemic alcohol using the borrowing hydrogen method under dual chiral catalytic conditions. In particular, we wish to probe whether both catalysts act in a cooperative/synergistic manner or in independent steps in a sequential manner in the catalytic cycle.
Several attempts to optimize any other conformer of D that maintains at least one of the stronger N–H⋯O interactions (a and b) did not succeed. Such geometries converged to the minima, wherein the free oxygen of the covalently bound phosphate participates in a number of C–H⋯O interactions. A comparison of the atoms-in-molecule (AIM) features of C and D revealed a well-woven network of noncovalent interactions such as N–H⋯O (a, b, and n), C–H⋯O (c, d, o, p, and q), and C–H⋯π (e, f, g, h, i, j, k, l, and m) and a unique N–H⋯π (r) interaction between the two catalysts. The ion-pair complex C has a larger number of C–H⋯π contacts, which are lost in D. Due to the loss of the stronger electrostatic N–H⋯O interactions, the free oxygen of the phosphate in D engages in a larger number of C–H⋯O contacts (c, o, p, and q). The overall stabilization of C by 16.6 kcal mol−1 over D can be rationalized by (i) the absence of the stronger interactions in D (a and b) due to conformational constraints and (ii) an enthalpic advantage of the ionic over the covalent interaction between the two catalysts.
The first important catalytic step is the dehydrogenation of the alcohol to form the corresponding ketone. The dual catalytic combination can act on the substrate alcohol, in two mechanistically distinct modes: inner-sphere or outer-sphere.24 The most preferred dehydrogenation is found to involve a concerted outer-sphere pathway,25 as shown as TS1 in Scheme 2. Interestingly, dehydrogenation in this mode enjoys a cooperative action of both phosphate and Ir(III) catalysts. The phosphate abstracts the proton of the hydroxyl group at one end while the electrophilic Ir(III) center simultaneously extracts the α-hydrogen in the form of a hydride from the other end of the substrate. The whole process can be considered as ‘borrowing’ of a molecule of hydrogen from the substrate in a heterolytic manner, thereby forming a ketone and a hydrogenated catalyst complex, i1. The dehydrogenation of the alcohol directly by the iridium catalyst B, in the absence of the phosphate counterion, is found to be much higher in energy than the cooperative mode involving the ion-pair C.26
Of the four stereoisomers of the substrate alcohol bearing two chirality centres, we focus on a pair of diastereomers (1(2S,3S) and 1′(2R,3S)) to probe the mechanism of the reaction as shown in Scheme 2. Here, the first and second configurational descriptors are respectively for the hydroxyl and the phenyl ring bearing chirality centres in the alcohol. The important issue at this juncture is to examine whether or not the catalyst dyad will dehydrogenate one or both diastereomers. The question is relevant as both catalysts and the substrate are chiral, giving rise to a chiral recognition or a match-mismatch situation. Interestingly, the relative Gibbs free energies (with respect to the separated reactants) of the transition states for concerted dehydrogenation of the two diastereomers of alcohol 1 are found to vary from 19.5 kcal mol−1 for 1(2S,3S) to 23.6 kcal mol−1 in the case of 1′(2R,3S).27
The ketone can now condense with aniline in the presence of the chiral phosphoric acid catalyst to form the respective imine. More details on the condensation steps are provided in Scheme S1 and Fig. S11 in the ESI.† After the loss of chirality at the α-C center due to the formation of the ketone and then the imine, only the configuration of the β-C remains unaltered as (S). Such imine intermediates can be protonated by the chiral phosphoric acid, leading to a second ion-pair species i2 with the chiral phosphate as the counterion (Scheme 2).28 The phosphate-bound iminium can then be intercepted by the [Ir]–H at the imino carbon. In other words, the hydrogen “borrowed” from the substrate by the catalyst dyad in the initial step of the reaction is effectively returned to the developing product in the final step. The hydride transfer to the ion-pair intermediate i2, as shown in Fig. 3, can lead to two diastereomeric products 3(2S,3S) and 3′(2R,3S) depending on the hydride addition to the re or the si prochiral face respectively. While both (E) and (Z)-imines can give both the products, the most favoured hydrogenation TS is found to be when the (E)-imine is involved. The relative energies and geometries of the (Z)-imine are provided in Fig. S15 in the ESI.†
The Gibbs free energy profile for the formation of the preferred stereoisomer 3(2S,3S) is provided in Fig. 2. It can be noticed that the hydride transfer to the iminium viaTS2SS is the rate-determining step with a barrier of 35 kcal mol−1 with respect to the exoergic PRC. Such a barrier can be regarded as a modest over-estimation for a reaction temperature of 110 °C. We have refined the energies of various species involved in the profile using larger basis sets, the details of which are given in Table S9 in the ESI.† The best estimate of the barrier for the rate determining step is found to be about 29.5 kcal mol−1 at the SMD(toluene)/B3LYP-D3/6-31G**,def2-TZVP (Ir)//SMD(toluene)/M06/6-31G**,SDD(Ir) level of theory. This energetic estimate is in line with the reaction conditions and further endorses the action of the two chiral catalysts in the critical hydride transfer step of the reaction. The transition state for the attack of [Ir]–H on the si-face of S-iminium is 3.6 kcal mol−1 higher than that on the re-face. Such a large difference in the energies between the transition states for the re face and si face hydrogenation corresponds to a de of 99.5% in favor of 3(2S,3S), which is in excellent agreement with the experimentally observed de of 94%. The origin of the energy difference can be analyzed at this stage, by probing the molecular features more closely.
Since the chirality at the second carbon chirality centre in the product gets fixed through this hydride transfer transition state, we have carefully analyzed the geometric features to understand more about the counterion-induced chirality transfer. The key interactions are electrostatic and are denoted by using letters a, b, and t in Fig. 3(i) and (ii). In both of these transition states, one of the phosphate oxygen atoms is involved in two relatively stronger N–H⋯O interactions (b and t) with both the electrophilic and nucleophilic partners. The other oxygen atom of the phosphate interacts strongly with one of the hydrogen atoms of the –NH2 group of the iridium catalyst (a). A network of weak non-covalent interactions, such as C–H⋯O (c, d, and v′), C–H⋯π (e, f, g, h, j, l, m, and m′), lpO⋯π (s) and π⋯π (z), between the 3,3′-substituents of the chiral phosphoric acid and the metal catalyst is noticed. The chiral counterion with its bulky aryl arms provides a chiral space for the electrophilic and nucleophilic partners to interact. The topological analysis reveals that the noncovalent interactions offered by the counterion help bring the two reacting partners together and hold them in a suitable geometric position. The mapping of the catalyst–catalyst interactions in both the diastereomeric transition states, as shown in Fig. 3, reveals that the C–H⋯O (c and d) and C–H⋯π (e, f, l, m, and m′) interactions between the chiral phosphate and the metal hydride are better in the case of the hydride transfer to the re face. An effective ‘communication’ between the catalysts, enabled by a series of noncovalent interactions, is noticed in the stereocontrolling transition state, in addition to the expected catalyst–substrate interactions (vide infra). We believe that such camaraderie between the catalysts is vital to their cooperative action in dual chiral catalysis.
Fig. 3 The general representation of the catalyst–catalyst interactions between the chiral phosphate counterion and Ir(III)–H complex (i) in TS2SS and (ii) in TS2RS. Two different orientations of the transition states are shown to improve the clarity of various weak interactions operating between the two chiral catalysts. The noncovalent contacts are in Å and the corresponding electron densities (ρ × 10−2 au) at the bcps are given in parentheses. Hydrogen atoms have been omitted wherever possible, to ensure better clarity. Only the interactions of electrophiles and nucleophiles with the phosphate are shown. Refer to Fig. 1 for atom colour codes. |
A comparison of the space-filling models of the transition states for stereoselective hydrogenation is provided in Fig. 4. The chiral phosphate counterion is represented using a grey shade and the Ir-catalyst in red. It can be noticed that both the catalysts together create a chiral pocket into which the substrate iminium fits in. An improved fit of iminium occurs when its re face is exposed for the [Ir]–H addition rather than the si face. In the most favoured transition state TS2SS for the re face addition of [Ir]–H, all the substituents of the iminium such as the N-phenyl and that on the β-chirality centre participate in noncovalent interactions with the catalyst dyad. However, in the higher energy TS2RS, the si face of the iminium should point toward the [Ir]–H. In such an arrangement the β-phenyl substituent on the chirality centre is found to remain away from the catalyst in a perpendicular arrangement. An interesting difference between these diastereomeric transition states in the region where the iminium fits in is also noticed. In the chiral pocket of the catalyst dyad, the upper region is relatively more compact than the lower region, as shown in Fig. 4. In TS2SS, the tetrahedral β-chirality centre occupies the spacious lower region of the catalyst leading to a relatively better fit of the substituents. In such a fit, better N–H⋯O and C–H⋯O interactions (t and w) between the iminium and the chiral phosphate are noticed. In the higher energy TS2RS, the larger tetrahedral carbon would have to fit into the compact upper region of the catalyst dyad. In order to accommodate the tetrahedral center, distortion in the dihedral angle ϕ1(C1–C2–C3–C4), up to −102° is noticed in the higher energy TS2RS as compared to −94° in the most preferred TS2SS.
Fig. 4 Space-filling models (left) of the diastereomeric transition states (i) TS2RS and (ii) TS2SS for the asymmetric hydride transfer to (S)-iminium. Given on the right side are all the important interatomic distances (Å) of noncovalent interactions and the corresponding electron densities (ρ × 10−2 au) at the bond critical points (bcps). Hydrogen atoms have been omitted wherever required, to ensure better clarity. Relative free energies are given in parentheses in kcal mol−1. Refer to Fig. 1 for atom color codes. |
The noncovalent interactions only between the catalysts and substrate are shown in Fig. 4.29 The most important interaction is an N–H⋯O electrostatic interaction (denoted as t) between the iminium and phosphate. In addition, the catalyst dyad exhibits a network of weak contacts such as C–H⋯O (u, v, w, w′, w′′) and C–H⋯π (x, x′, x′′, y, y′, y′′) with the substrate through (i) 3,3′-aryl substituents of the phosphate and (ii) the substituents on the Ir-catalyst. In the hydride transfer transition state, there are only two or three interactions between the substrate and the [Ir]–H whereas six interactions are noted with the chiral phosphate. This led us to infer that the chiral phosphate is mostly, if not entirely, responsible for the stereochemical outcome in this dual catalytic reaction. Since chiral induction by the phosphate is recognized as more impactful than that by the chiral iridium complex, we propose that further modifications to this dual catalytic strategy could focus on the phosphoric acid 3,3′-aryl substituents. It would also be of interest to see the stereochemical outcome with a substrate alcohol bearing one more chirality centre.
The difference in the catalyst–catalyst and catalyst–substrate non-covalent interactions in the two diastereomeric hydride addition TSs from the same (S)-iminium phosphate ion-pair intermediate, i2(S), is found to be an important factor impacting the diastereoselectivity. However, enantiomers of the product amine are formed from different intermediates, i.e., 3(2S,3S) from i2(S) and 3′′(2R,3R) from i2(R). In other words, intermediate i2(S) cannot give the minor enantiomer 3′′, neither can i2(R) form the major enantiomer. It should also be noted that i2(S) and i2(R) are diastereomeric intermediates and have different energies. Therefore, the standard approach for calculating ee, on the basis of the difference in the relative free energies of the hydride addition TSs, does not directly apply in this situation. In most asymmetric reactions, enantioselectivity is achieved through a crucial step where one of the prochiral faces of the substrate reacts faster than the other. In such cases, all elementary steps leading up to the stereodetermining TS are expected to show no particular preference towards the formation of any one enantiomer. In the present example, the diastereoselectivity in the formation of 3(2S,3S) is due to the difference in the free energies of the corresponding hydride transfer TS only. However, the question of enantioselectivity presents a non-trivial scenario. The initial dehydrogenation and ensuing condensation as well as the hydride addition might impact the enantioselectivity, as the R or S configuration at the β-C of the substrate is preserved in these steps. Since the yield of the major product under the optimised reaction conditions is more than 80%, one can infer the involvement of a DYKAT mechanism,30 in which case, at some point during the course of the reaction, the R substrate is expected to convert to the S substrate prior to the final hydride transfer step.
We have considered an epimerization of the β-R stereocentre to β-S at the ketone stage, i.e., in intermediate i1′ (Scheme 3).31 The final product is formed through a condensation of the (S)-ketone with aniline and a hydride addition to the resulting (S)-iminium intermediate. The two key steps in condensation are (a) nucleophilic attack of aniline on the carbonyl carbon viaTSi-k(R) giving intermediate k(R), and (b) H2O elimination to form (R)-iminium l(R) viaTSk-l(R). The epimerisation to the (S)-ketone also involves two steps such as (a) a deprotonation of the β-C–H mediated by the chiral phosphoric acid, and (b) re-protonation through the other prochiral face of the enol j to give the (S)-ketone. The epimerization of the (R) to the (S) ketone via a common enol intermediate is found to be energetically more favoured than its condensation to form (R)-iminium (Fig. 5).
Scheme 3 Two plausible pathways for the involvement of the (R)-ketone in the reaction. The solid arrows represent a faster pathway as compared to the dotted arrow. |
Fig. 5 Comparison of the free energy profiles (kcal mol−1) for epimerisation of the (R) to the (S) ketone and the condensation of the (R)-ketone to form the (R)-iminium intermediate. |
As the calculation of enantiomeric excess in a DYKAT situation is not quite straightforward, we approached the problem as follows. Two critical mechanistic routes that can impact the configuration of the products are (a) the ease of conversion of i1′(R) to i1′(S) followed by condensation to form (S)-iminium and subsequent hydride addition to give the (2S,3S) product, or (b) a direct condensation of i1′(R) to form the respective iminium and then the hydride transfer to (R)-iminium to form the minor stereoisomer 3′′(2R,3R). In the present example, the configuration of the product is determined in the phosphoric acid catalysed epimerization of the (R)-ketone to the (S)-ketone through TSR-j and TSj-S (Scheme 3). The Gibbs free energy profile, as given in Fig. 5, indicates that the most energy demanding elementary step in the direct condensation path-(b), which leads to the 3R configuration in the product, is the expulsion of water viaTSk-l(R). On the other hand, ketone epimerization viaTSj-S in path-(a) would result in the 3S configuration. It should be noted that the configuration at the C2 centre would depend on the preferred prochiral face involved in the hydride addition. Hence, the energetic origin of enantioselectivity can be considered as the difference in the relative Gibbs free energies of TSs with the highest energy in the epimerization and condensation pathways for the (R)-ketone, as shown in Scheme 3. This difference amounts to (26.2–22.7) 3.5 kcal mol−1, which is equivalent to an ee of 99.4%, in very good agreement with the experimental value of 98%.32
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c8sc01433g |
This journal is © The Royal Society of Chemistry 2018 |