Jens H.
Aasheim
a,
Heike
Fliegl
b,
Einar
Uggerud
b,
Tore
Bonge-Hansen
*a and
Odile
Eisenstein
*bc
aDepartment of Chemistry, University of Oslo, PO Box 1033, Blindern, 0315 Oslo, Norway. E-mail: tore.hansen@kjemi.uio.no; Fax: +47 228 55441; Tel: +47 228 55446
bCentre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, P. O. Box 1033, Blindern, 0315 Oslo, Norway
cInstitut Charles Gerhardt, UMR 5253 CNRS Université Montpellier 2, cc1501, Place E., Bataillon, 34095 Montpellier, France. E-mail: odile.eisenstein@univ-montp2.fr; Fax: +33 467144839; Tel: +33 467143306
First published on 26th September 2014
A prospective novel organocatalyst scaffold has been investigated by computational studies for its potential to induce stereoselective reactions. Density functional theory calculations (DFT) of an organocatalysed Diels–Alder reaction between cinnamaldehyde and cyclopentadiene, show that it is possible to induce selectivity by placing appropriate chemical groups on the scaffold. In the present case small substituents with strong electronic effects are found to be more effective in controlling the stereochemical induction than sterically bulky substituents. We show, using noncovalent interaction analysis (NCI plots), that the CF3 group is especially efficient in improving the stereochemical control as well as increasing the reaction rate.
One of the first studies of organocatalysis was carried out by MacMillan who used MM3 Monte Carlo calculations.2 Following this study, computational chemistry, has been considered as a helpful tool in both describing, explaining and also validating organocatalytic structures, selectivity and mechanisms. Additional relevant and notable examples are by the groups of Houk8 and Tomkinson.9 It is thus attractive to consider the possibility of using computational methods to go beyond rationalization and envisage the prediction and rational catalyst design as shown for example by Tomkinson,9 Houk,10 Moitessier11 and also one of us.12 One essential aspect of organocatalysis is the representation of the weak interactions that exist between the various groups of the reactants. This requires the use of a methodology that properly represents dispersion forces. The group of Truhlar has been especially active in this field, and has developed several levels of the Minnesota functional that is appropriate for these kinds of studies.13 Among the various functionals of this family, M06-2X, parameterized for chemical species without metal atoms, is well adapted for applications involving main group thermochemistry, kinetics and non covalent interactions.13b It has been used successfully for studying organocatalytic and cycloaddition reactions.14,15 This functional has been selected for the present study.
The main purpose of this work is to study the possibility of reactivity manipulation, by creating a synthetically viable organocatalytic system that includes adjustable electronic and steric effects. The core of the system is a rigid bicyclic structure, to limit the degrees of freedom and thus provide a better control of the conformational space in addition to significantly simplifying the calculations. Conformational control will also provide a means to shield one side of the reactive site, hence yielding stereochemical control.
A possibly interesting system is the bicyclic Geissman–Waiss lactone (BL) (Fig. 1) as a candidate for this exploratory study. This lactone provides a bicyclic skeleton with a sp2-hybridised carbonyl,16 providing the needed structural rigidity. The reaction that is considered is shown in Scheme 1. The amine group is transformed into an iminium through a condensation reaction. Then, the reaction between the iminium carrying the reactive group R and the substrate (reaction between R and L giving R′) occurs in a chiral environment that could induce stereocontrol, which can be further influenced by the nature of the group at the 4 position of the lactone. Hydrolysis of the iminium liberates the end product and regenerates the amine form of BL, which re-enters the catalytic cycle.
With these basic structural elements in place we would have a tunable system that could be modified relatively easy by changing the aromatic shielding group. The reaction that we consider is a Diels–Alder reaction in which the iminium is conjugated to a CC double bond that serves as a dienophile.
In accordance with other theoretical studies of the Diels–Alder reaction,8e,23 all transition states correspond to an asynchronous concerted reaction pathway.
The two new carbon–carbon bond lengths differ by about 0.2 Å for both the endo and the exo reactions, as shown in Fig. 2. The calculations show a slight preference for the exo reaction, in contrast to the normal endo preference usually attributed to secondary orbital interactions between the incoming CP and the carbonyl group.24 A possible reason for this difference is the presence of additional secondary interactions between the incoming CP and the π-system of the phenyl ring in the exo reaction.
The catalysed version of the reaction, following the proposal presented in Scheme 1, is shown in Scheme 3. The condensation of the cinnamaldehyde and the bicyclic lactone (BL) catalyst yields the conjugated iminium ion (Scheme 2), which can serve as an activated dienophile in the Diels–Alder reaction. Hydrolysis recovers the carbonyl group, in the product, and the amine functionality of BL. This study focuses on the cycloaddition step, as it is the rate determining step of the catalytic cycle.9c
Scheme 3 The catalytic cycle of the organocatalysed Diels–Alder reaction between CP and the E-isomer of the iminium ion formed from cinnamaldehyde and BL. |
To induce stereoselectivity in Diels–Alder reactions (Scheme 3) one needs to control three aspects: (i) the formation of the E- or Z-iminium ion isomer; (ii) endo- or exo-attack; (iii) attack on either side of the iminium ion, which will be referred to as top and bottom.
The condensation reaction between cinnamaldehyde and the bicyclic lactone (BL) catalyst (1) forms an iminium ion with either an E (H on imine carbon cisoid to the lactone) or a Z conformation. The calculations show a preference of 2 kcal mol−1 for the 1-E over the 1-Z isomer, indicating that the condensation is conformationally controlled at the CN double bond period (Fig. 3). The reason for this slight energy preference is likely to come from the close proximity of the phenyl ring at C4 of the lactone (see Fig. 1 for labelling) and the π system of the imimium moiety in the Z-conformation. It remains to be determined if this preference for the E over Z isomer translates into the transition states, and contributes to differentiate the top from bottom faces of the iminium and the endo from the exo pathways.
Fig. 3 E/Z iminium ions formed by the condensation of 1 and cinnamaldehyde with relative Gibbs energies in kcal mol−1. |
Fig. 4 displays the preferred transition states TS1 for the Diels–Alder reaction between the 1-E iminium ion and CP. They all correspond to asynchronous concerted reactions and have activation barriers that are of similar magnitude to that of the uncatalyzed reaction (Fig. 2). The Gibbs energies of activation are the differences in Gibbs energies of the transition state and of the most stable adduct formed between 1-E and CP. The lowest transition state is only 2 kcal mol−1 lower (19.2 versus 21.6 kcal mol−1 for the catalyzed and uncatalyzed systems, respectively). Thus the iminium and the aldehyde forms of the activated dienophile have similar reactivities.
Fig. 4 Transition states for the reaction between CP and iminium ion 1-E with Gibbs energies of activation in kcal mol−1 and relevant distances in Å. |
The exo attack is slightly preferred over the endo approach by 2 kcal mol−1 on average. The preference for the exo attack is thus somewhat larger in the presence of BL (compare Fig. 2 and 4). Unfortunately, the exo attack is not associated with a significant differentiation for approach on one or another of the two faces. This is illustrated by a negligible difference in energy of only 0.1 kcal mol−1 between the top and bottom transition states. In contrast, the endo attack differentiates between the two faces, and the bottom approach is preferred by 3.3 kcal mol−1 over the top approach (Fig. 4). This indicates that the phenyl group is quite effective in directing the reaction towards one side of the iminium moiety for the endo attack, but unfortunately not for the preferred exo approach. This is due to the inherent flexibility of the iminium chain and some limited flexibility at the C4 phenyl group.
Before considering the effects of further modification of the BL catalysts, the following points have been established: (i) the BL-structure contributes to differentiate addition to the top and bottom faces for the endo transition states, but not sufficiently for the preferred exo transition states; (ii) the energy preference for the E- over Z-iminium ion isomer is conserved in the transition states; (iii) the exo attack of CP on the E isomer is preferred over the endo one. Finally, the catalyst has only a small impact on the rate of reaction since the activation barriers are lowered by only 2 kcal mol−1 compared to the uncatalyzed systems.
The next topic concerns the effect of substituents R1, R2, R3 and R4 (Fig. 5) on the activation energy and selectivity of the reaction (endo vs. exo and top vs. bottom). The Gibbs energies of the eight transition states relative to the most stable adduct formed between CP and each of the relevant iminium ions are given in Table 1. The associated endo:exo ratio and ee for both the endo and exo addition, obtained from a Boltzman analysis carried out at 298 K, is given in Table 2.
Entry | TS-endo-top | TS-endo-bottom | TS-exo-top | TS-exo-bottom | |
---|---|---|---|---|---|
1 | TS1-E | 23.3 | 20.0 | 19.3 | 19.2 |
TS1-Z | 22.5 | 21.9 | 23.2 | 21.2 | |
2 | TS2-E | 22.7 | 20.2 | 19.1 | 18.9 |
TS2-Z | 22.5 | 23.3 | 23.3 | 23.1 | |
3 | TS3-E | 18.7 | 19.0 | 17.2 | 17.5 |
TS3-Z | 20.5 | 19.8 | 19.4 | 18.8 | |
4 | TS4-E | 23.4 | 20.6 | 19.7 | 19.9 |
TS4-Z | 24.2 | 22.8 | 24.5 | 22.6 | |
5 | TS5-E | 25.0 | 20.7 | 21.4 | 20.2 |
TS5-Z | 19.8 | 18.6 | 20.6 | 17.6 | |
6 | TS6-E | 23.8 | 19.8 | 19.5 | 18.5 |
TS6-Z | 23.1 | 21.3 | 23.2 | 20.7 | |
7 | TS7-E | 22.1 | 19.5 | 19.1 | 18.5 |
TS7-Z | 21.8 | 20.9 | 22.0 | 20.5 | |
8 | TS8-E | 18.4 | 18.6 | 15.5 | 17.5 |
TS8-Z | 18.4 | 19.7 | 20.6 | 19.6 |
endo:exo | ee – endo | ee – exo | |
---|---|---|---|
a S and R refer to the absolute configuration at the two carbons of the dienophile. | |||
1 | 1:6.9 | 91.5 S,Sa | 2.1 S,S |
2 | 1:17 | 96.1 S,S | 16.7 S,S |
3 | 1:13 | 26.7 R,R | 23.0 R,R |
4 | 1:7.9 | 94.0 S,S | 22.1 R,R |
5 | 1:4.9 | 71.2 R,R | 96.4 R,R |
6 | 1:9.5 | 86.7 S,S | 62.6 S,S |
7 | 1:6.3 | 82.6 S,S | 39.4 S,S |
8 | 1:54 | 22.6 S,S | 93.6 R,R |
The results of Table 1 show that the Gibbs energies of activation for the Diels–Alder reaction vary between 15 and 25 kcal mol−1, which indicates that the nature of the substituents influence the global rate of the reaction. An electron-withdrawing group close to the iminium moiety (R2 = CF3, entry 8) significantly lowers the energy of all transition states. This is consistent with the well-known accelerating effect of electron acceptor groups on the dienophile.2,9,25 Interestingly, changing R3 from Ph to tBu (entry 3) also gives relatively low activation barriers, especially for the addition to the preferred E-isomer. In total, the influence of R3 on the activation energy is relatively modest, indicating that the BL catalyst does not significantly affect the reaction rate.
Table 2 shows that the exo addition is always preferred, but the endo:exo ratio varies between 1:4.9 (entry 5) to 1:54 (entry 8). These two entries giving the lowest and highest endo:exo ratio, correspond to a CH3 substituent at R2 for entry 5 and CF3 for entry 8. Since these two groups have similar effective vdW radii,26 this indicates that other factors influence the selectivity. The significant lowering of the activation barriers with CF3, and notably for the exo addition, could be responsible for this result. Changes in positions R3 and R4 have less influence on the endo:exo ratio. Even replacing Me with F at R4 gives similar endo:exo ratio (entries 1 and 7). However, increasing the bulk of R4, as for entry 2 with R4 = m-xylene and entry 3 with R4 = tBu, gives larger exo:endo ratio than R4 = phenyl.
Entry 1 is rather unique in giving such a small ee for the exo approach. For all other entries, the ee ratio is larger and can be even very large like for entries 5 and 8. The ee from the endo addition is always very large with the exception of entries 3 and 8. However, we do not need to discuss further the endo addition, as it is not preferred.
The most attractive case appears to be entry 8. The structure of the lowest transition state is given in Fig. 6. As mentioned earlier, this entry is associated with an increase in the reaction rate and a high exo selectivity, in addition to high enantiomeric selectivity. Lowering of the energy barrier is an expected result of the presence of an electron withdrawing group on the dienophile, as known in the Alder rule.2,9,25 This does not however, account for the high selectivity in the exo addition and the significant preference for adding to the top face of the iminium ion. In order to obtain a deeper insight into this issue, an analysis of the noncovalent interactions was performed for the four E iminium ion transition states of entry 8 with NCIplot.21
Fig. 6 The preferred transition states (TS2 to TS8) for entries 2 to 8 (Table 1) with Gibbs energies of activation in kcal mol−1, and selected distances in Å. |
Fig. 7 shows the result of the NCI analysis for the lowest transition state (TS8-E-exo-top) associated with entry 8. The weak interactions revealed by the NCI analysis are color coded from blue, for the strongest attractive weak interactions (hydrogen bonds, for example), to red for repulsive ones (steric clashes). In between, the very weak van der Waals interactions appear in green. Focus is mostly placed on the spaces that define the selectivity of the reaction. Fig. 7 shows that large green surfaces occur in the space between the incoming CP and both the R4 phenyl group and the R2 CF3 group of BL. They are clearly associated with weak interaction between the CH bonds of CP and both the electron density of the phenyl ring and the electron density at the fluorine atoms. These interactions are representative of non-classical CH⋯X interactions, whose importance in determining the structures of large systems have been established.27 Unfortunately, there is not yet any quantitative relation between the size of the surface combined with its colour, and the energy of the interaction. However, a decrease of the size of the surface and a change of colour from darker to lighter green is indicative of a decrease in the magnitude of the weak interactions.
Fig. 7 Color-coded noncovalent interactions (NCI) surfaces (attractive decreasing from blue to green, repulsive increasing from yellow to red) for TS8-E-exo-top. |
Thus, going to the less preferred transition states for entry 8 is associated with a decrease, or disappearance, of the green surfaces. In particular going from TS8-E-exo-top to TS8-E-endo-top, which corresponds to a raise of 2 kcal mol−1, (Table 1), shows a clear decrease of the green surfaces between CP and the phenyl and CF3 groups (Fig. 8). It seems that in the exo orientation, the CH bonds of CP creates the best network of dispersion interactions with the two electron rich groups. In the case of the bottom attacks, these weak interactions between CP and the phenyl and CF3 groups are lost and the interaction between the R4 phenyl group and the iminium skeleton does not fully compensate (Fig. 9).
Fig. 8 Color-coded noncovalent interactions (NCI) surfaces (attractive decreasing from blue to green, repulsive increasing from yellow to red) for TS8-E-endo-top. |
Fig. 9 Color-coded noncovalent interactions (NCI) surfaces (attractive decreasing from blue to green, repulsive increasing from yellow to red) for TS8-E-exo-bottom. |
The NCI analysis suggests the paradoxical result that it could be preferable to attack the iminium from the most hindered side because approach on this side is in fact dominated by weak interactions. One can thus understand the influence of the CF3 group that modifies the electron distribution and increases the electron density in the appropriate part of the space and enhances the attractive interactions.
Proper substitutions on the BL catalyst (for example by introducing an electron withdrawing group such as CF3) highlight the possibility to direct selectivity by creating an electronic environment that could fit the approaching reagent, an effect somewhat related to enzymatic control. The ideal system may have to be adapted to the reaction and even to the two substrates. This is not surprising in regard of the very high specificity of enzymatic catalysis.
Footnote |
† Electronic supplementary information (ESI) available: Cartesian coordinates of all reported structures as well as the total electronic and zero-point energies. See DOI: 10.1039/c4nj01460j |
This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2014 |