Andrea
Škvorcová
and
Radovan
Šebesta
*
Department of Organic Chemistry, Faculty of Natural Sciences, Comenius University, Mlynska dolina CH-2, SK-84215, Bratislava, Slovak Republic. E-mail: radovan.sebesta@fns.uniba.sk; Tel: +421 2 60296208
First published on 4th October 2013
Diastereoselective ortho-lithiations serve for the preparation of many important chiral ferrocenes, however diastereoselections of these lithiations are explained only by simple steric models. We elucidated ortho-lithiations of three important ferrocenes based on DFT calculations. The calculations showed that simple models of transition states involving ferrocenyl substrates and lithium bases can only in some cases account for the experimental results. Transition state models, which include solvent or coordinating additives as distinct entities bound to lithium, can satisfactorily explain diastereoselection of the ortho-lithiations of chiral ferrocenes.
Scheme 1 Formation of two possible lithium intermediates evolved from a chiral ferrocene derivative. |
Ugi and co-workers described the first diastereoselective ortho-lithiations of ferrocene derivative – N,N-dimethyl-1-ferrocenylethylamine (Ugi amine).10 Based on this reaction Hayashi developed chiral amino phosphanes, which later culminated in the development of ferrocenyl diphosphines such as the Josiphos family of chiral ligands.11 Apart from the Ugi amine, other amines can serve as ortho-directing groups.12–14 Nitrogen, as a stronger nucleophile, serves as a directing atom in oxazolines, which form the basis of another important group of chiral ligands. Three research groups described ortho-lithiations of ferrocenyl oxazolines almost simultaneously.15–17 Other nitrogen based ortho-directing groups are hydrazones,18,19 azepines,20 sulfoximines,21 and imidazolines.22 Oxygen performs the ortho-directing function in sulfoxides23 acetals,24 methoxy-derivatives,25 and alcohols.26 The structures of the major products of diastereoselective ortho-lithiations have been well established experimentally, either by X-ray or NMR; however, a detailed understanding of the reasons for the observed diastereoselectivity is lacking. It is usually explained by empirical consideration based on simple models of transition states or even intermediates. DFT calculations, only recently started to help explain the underlying reasons for diastereoselection in ortho-lithiations of ferrocene derivatives.27,28 We have studied diastereoselective ortho-lithiations of ferrocenophanes29 and have suggested the mechanism of dilithiation reaction leading to the Taniaphos family of ligands.30 A lack of understanding of the lithiation of ferrocenes prompted us to investigate diastereoselective ortho-lithiations of three relevant ferrocene compounds by DFT calculations. Although many factors, such as steric demands or aggregation of alkyllithium reagent, solvent, additive and size of the substituent on a chiral stereogenic centre, seem to influence the diastereoselectivity of ortho-lithiations we argue that even a rather simplistic model including only ferrocene substrate and a lithium base can account, at least on a qualitative level, for the observed diastereoselection. On the basis of these facts, we decided to provide a theoretical explanation for the experimentally observed diastereoselective lithiations on ferrocene derivatives 1–3, which are essential for the synthesis of useful chiral ferrocenyl ligands (Fig. 1).
For the compounds 1–3, Table 1 summarizes experimental results obtained by ortho-lithiation.3
Compound | Base | Additive | E+ | Solvent | dr | Major diastereomer |
---|---|---|---|---|---|---|
(R)-1 | n-BuLi | — | PPh2 | Et2O | 96:4 | (R,Rp) |
(S)-2 | s-BuLi | — | Me3Si | THF | 8:1 | (S,Sp) |
(S)-2 | s-BuLi | TMEDA | Me3Si | Hexane | 500:1 | (S,Sp) |
(S)-3 | LDA | — | Me3Si | THF | 99:1 | (S,Rp) |
Compound | ΔErela (kJ mol−1) | Compound | ΔErelb (kJ mol−1) | Compound | ΔErel (kJ mol−1) |
---|---|---|---|---|---|
a RI-DFT, BP86, def2-TZVP. b RI-DFT-COSMO (ε = 4.33), BP86, def2-TZVP single point. | |||||
(R)-1 + n-BuLi | 83.5 | (R)-1 + n-BuLi | 45.2 | (R)-1 + n-BuLi + 2Et2O | 95.5 |
— | — | — | — | (R)-1 + n-BuLi·2Et2O | 0.2 |
(R)-1·n-BuLi (C1a) | 0.0 | (R)-1·n-BuLi (C1a) | 0.0 | (R)-1·n-BuLi·2Et2O (C2a) | 0.0 |
(R)-1·n-BuLi (C1b) | 2.5 | (R)-1·n-BuLi (C1b) | 3.6 | (R)-1·n-BuLi·2Et2O (C2b) | 1.8 |
(R,Rp)-TS1 | 59.3 | (R,Rp)-TS1 | 72.7 | (R,Rp)-TS2 | 62.1 |
(R,Sp)-TS1 | 68.7 | (R,Sp)-TS1 | 77.1 | (R,Sp)-TS2 | 73.1 |
(R,Rp)-Int1 + BuH | −25.0 | (R,Rp)-Int1 + BuH | −38.0 | (R,Rp)-Int2 + BuH | −66.0 |
(R,Sp)-Int1 + BuH | −21.9 | (R,Sp)-Int1 + BuH | −33.9 | (R,Sp)-Int2 + BuH | −52.7 |
The complex C1a leads to transition state (R,Rp)-TS1 and C1b to transition state (R,Sp)-TS1. From the two calculated diastereomeric transition states, (R,Rp)-TS1 is energetically more favoured (ΔE = 12 kJ mol−1) than (R,Sp)-TS1. This fact is in accordance with the experimentally observed configurations of products. The fact that transition state TS1 indeed connects starting material with product was confirmed by performing a potential energy scan (see ESI†).
Experimental data show that the course of lithiation often heavily depends on the solvent. The importance of solvation of lithium compounds is well known. Moreover, solvation can influence the formation of transition states or stabilization of the resulting Li-intermediates.34 Taking this fact into account, we proposed to simulate the effect of solvation of Li-derivatives firstly as a dielectric medium and secondly by two molecules of Et2O, which are included not just as a dielectric medium, but as defined molecules coordinated to lithium.
The calculation performed with COSMO (Et2O, ε = 4.33) resulted in an energy difference of transition states TS-1 of just 4 kJ mol−1. This result suggests that the solvent in the lithiation did not have influence just as a dielectric medium, but that the solvent molecules needed to be explicitly coordinated on lithium with a direct effect on diastereoselectivity of the reaction.
Interestingly, the inclusion of solvent resulted in approximately the same energy difference between the two diastereomeric transition states. A significant change was observed in the relative energies of lithiated products (R,Rp)-Int2, (R,Sp)-Int2. The energy difference increased to ΔErel = 13 kJ mol−1. The inclusion of solvent in the calculation did not affect the formation of transition states but a stabilization of Li-intermediates.
By examination of electron density maps of the calculated transition states TS1, it was possible to see steric repulsions, which were responsible for differences between the transition states. In accordance with Ugi's explanation, there indeed seems to be repulsion of the methyl group with a ferrocene moiety in the non-preferred transition states (R,Sp)-TS1 (Fe–H(CH3) 3.216 Å) and (R,Sp)-TS2 (Fe–H(CH3) 3.173 Å) (Fig. 3). The methyl group, in preferred transition states, is located in the plane of the cyclopentadienyl (Cp) ring.
The repulsion of methyl group by the ferrocene skeleton can be visualized by a NCI plot (Fig. 4).35
The alkyllithium bases are present as aggregates in solution. As this can also influence the formation of transition states, we performed a calculation of the lithiation of amine (R)-1 with n-BuLi-dimer. Although in n-BuLi solutions aggregates up to a hexamer can exist, the presence of an amine coordinating atom, such as nitrogen, will likely decrease the aggregation (Fig. 5).
Therefore we used the n-BuLi-dimer as a model for the aggregated species. The energy difference between the two diastereomeric transition states increased to 18 kJ mol−1 and the correct diastereomer was predicted as being more favourable. The reason for the bigger energy difference between the transition states is probably the increased steric interactions in the presence of n-BuLi-dimer in complexes and transition states. The relative energies of the reactants, the transitions states and the resulting Li-intermediates of the lithiation of amine (R)-1 are summarized in Table 3.
Compound | ΔE (kJ mol−1) |
---|---|
(R)-1 + (n-BuLi·2Et2O)2 | 12.4 |
(R)-1·(n-BuLi)2·3Et2O (C3) + Et2O | 0.0 |
(R,Rp)-TS3 + 3Et2O | 71.7 |
(R,Sp)-TS3 + 3Et2O | 89.7 |
(R,Rp)-Int2 + BuH + (n-BuLi·2Et2O) | −5.9 |
(R,Sp)-Int2 + BuH + (n-BuLi·2Et2O) | 7.4 |
The calculations suggest that this simple model may not be sufficient to account for the observed diastereoselection. Diastereomeric transition states TS4 composed of oxazoline 2 and s-BuLi showed only a negligible energy difference ∼1 kJ mol−1. The situation manifests itself already in the calculation of complexes. The energy difference between complexes C4a and C4b was just 0.27 kJ mol−1. As is shown in Table 1, a major influence of THF was experimentally observed. The calculation of transition states, where two molecules of THF are included, confirmed the experimental results. Also, in this case it is not sufficient to include solvent only as the dielectric medium. The energies of the transition states were the same and neither transition state was favoured. However, the experimentally observed diastereoselectivity of the lithiation of compound 2 in THF is 8:1.
The calculation of lithiation, where two THF molecules were coordinated on the lithium atom observed a preferred formation of the less hindered transition state (S,Rp)-TS5 by 7.4 kJ mol−1. While in the transition states TS4, we did not observe any repulsion of the i-Pr group with base or ferrocene skeleton, in TS5, THF molecules showed repulsion of the i-Pr group with approaching base. This fact made one transition state more favourable and the diastereoselectivity changed. This theoretical result confirmed the experimental observations of Sammakia and co-workers.39
The inclusion of strongly chelating agent N,N,N′,N′-tetramethylethylenediamine (TMEDA) had a significant impact on the results. As this sterically bulky compound coordinated to lithium, it filled the space around lithium and thus helped increase steric interactions between the base and the substrate.
This effect was manifested both in complexation and was followed by the formation of transition states. It is confirmed by the large energy difference between complexes C6a and C6b. The formation of C6b was less favoured and therefore the formation of (S,Sp)-TS6 was not preferred. The calculated transition state (S,Rp)-TS6 is by 19 kJ mol−1 more stable than (S,Sp)-TS6. These calculated results correlate with the experimentally obtained major diastereomer (Fig. 7). The relative energies of the reactants, the transitions states and the resulting Li-intermediates of the lithiation reactions of oxazoline (S)-2, both with and without TMEDA and THF, are summarized in Tables 4 and 5.
Compound | ΔEa (kJ mol−1) | Compound | ΔEb (kJ mol−1) |
---|---|---|---|
a RI-DFT, BP86, def2-TZVP. b RI-DFT-COSMO (ε = 7.58), BP86, def2-TZVP single point. | |||
(S)-2 + s-BuLi | 86.9 | (S)-2 + s-BuLi | 37.5 |
(S)-2·s-BuLi (C4a) | 0.0 | (S)-2·s-BuLi (C4a) | 0.0 |
(S)-2·s-BuLi (C4b) | 0.3 | (S)-2·s-BuLi (C4b) | 0.1 |
(S,Rp)-TS4 | 45.6 | (S,Rp)-TS4 | 57.8 |
(S,Sp)-TS4 | 46.4 | (S,Sp)-TS4 | 57.9 |
(S,Rp)-Int3 + BuH | −65.0 | (S,Rp)-Int3 + BuH | −79.1 |
(S,Sp)-Int3 + BuH | −65.7 | (S,Sp)-Int3 + BuH | −80.7 |
Compound | ΔE (kJ mol−1) |
---|---|
(S)-2 + s-BuLi + 2THF | 124.7 |
(S)-2 + s-BuLi·2THF | 6.3 |
(S)-2·s-BuLi·2THF (C5a) | 0.0 |
(S)-2·s-BuLi·2THF (C5b) | 14.6 |
(S,Rp)-TS5 | 38.6 |
(S,Sp)-TS5 | 46.0 |
(S,Rp)-Int4 + BuH | −106.6 |
(S,Sp)-Int4 + BuH | −105.4 |
(S)-2 + s-BuLi + TMEDA | 119.4 |
(S)-2 + s-BuLi·TMEDA | 7.7 |
(S)-2·s-BuLi·TMEDA (C6a) | 0.0 |
(S)-2·s-BuLi·TMEDA (C6b) | 18.1 |
(S,Rp)-TS6 | 39.9 |
(S,Sp)-TS6 | 59.2 |
(S,Rp)-Int5 + BuH | −106.3 |
(S,Sp)-Int5 + BuH | −105.0 |
Inspection of the electron density maps of diastereomeric transition states TS4 shows that the distance between i-Pr group oriented up and alkyllithium base is sufficiently large (2.92 Å) and no repulsion is observed. On the other hand, the presence of TMEDA, as a bulky additive coordinated to the alkyllithium base, leads to larger repulsion with the i-Pr group in (S,Sp)-TS6 and thus the formation of (S,Rp)-TS6 is favoured (Fig. 8).
At the beginning of this part of the study, we have realized that LDA conformations play an important role and significantly influence the results. Energy differences of up to ΔErel = 5 kJ mol−1 were found (see ESI†). The formation of complex C7a is more favored (3.5 kJ mol−1) than the formation of complex C7b (Fig. 9 and Table 6).
Compound | ΔEa (kJ mol−1) | Compound | ΔEb (kJ mol−1) | Compound | ΔEa (kJ mol−1) |
---|---|---|---|---|---|
a RI-DFT, BP86, def2-TZVP. b RI-DFT-COSMO (ε = 7.58), BP86, def2-TZVP, single point. | |||||
(S)-3 + LDA | 88.5 | (S)-3 + LDA | 36.4 | (S)-3 + LDA + 2THF | 114.3 |
— | — | — | — | (S)-3 + LDA·2THF | 5.9 |
(S)-3·LDA (C7b) | 0.0 | (S)-3·LDA (C7b) | 0.0 | (S)-3 + LDA·2THF (C8b) | 0.0 |
(S)-3·LDA (C7a) | 3.5 | (S)-3·LDA (C7a) | 1.6 | (S)-3 + LDA·2THF (C8a) | 6.5 |
(S,Rp)-TS7 | 68.8 | (S,Rp)-TS7 | 74.8 | (S,Rp)-TS8 | 49.5 |
(S,Sp)-TS7 | 62.0 | (S,Sp)-TS7 | 61.4 | (S,Sp)-TS8 | 30.7 |
(S,Rp)-Int6 + (i-Pr)2NH | 55.5 | (S,Rp)-Int6 + (i-Pr)2NH | 26.5 | (S,Rp)-Int7 +(i-Pr)2NH | 1.1 |
(S,Sp)-Int6 + (i-Pr)2NH | 40.8 | (S,Sp)-Int6 + (i-Pr)2NH | 12.5 | (S,Sp)-Int7 +(i-Pr)2NH | −9.5 |
Diastereoselection seems to be induced by the orientation of a tosyl group. While the tosyl group in C7a and (S,Sp)-TS7 is oriented over the ferrocene skeleton, the position of this group in C7b and (S,Rp)-TS7 is near the ferrocene moiety (Fig. 10).
The inclusion of the solvent in the calculation increased the energy difference between the complexes and the transition states. The coordination of two molecules of THF on lithium increased the steric demands of the lithium base. The repulsion of LDA with a tosyl group is observed in non-favoured transition state (S,Rp)-TS8. On the other hand, this repulsion is not possible in the transition state (S,Sp)-TS8, because the tosyl group is oriented on the opposite side. The energetically favoured transition states lead to the product (S,Rp)-Int7, which is in accordance with experiment (Fig. 9 and Table 6).
The electron density maps of transition states help observe interactions in the diastereomeric transition states. Table 7 summarizes the bonding distances, energy differences for complexes, transition states and intermediates of the studied lithiations.
(R)-1 | (S)-2 | (S)-3 | ||||||
---|---|---|---|---|---|---|---|---|
Int | Int1 | Int2 | Int2 e | Int3 | Int4 | Int5 | Int6 | Int7 |
a Energy of calculated complexes, transition states and intermediates are displayed in kJ mol−1. b Distances between hydrogen and base in complexes displayed in Å. c Distance between hydrogen and ferrocenyl carbon displayed in Å. d Distance between ferrocenyl carbon and lithium displayed in Å. e Intermediates produced when BuLi-dimer was considered. | ||||||||
C | C1 | C2 | C3 | C4 | C5 | C6 | C7 | C8 |
ΔE C-a | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 3.5 | 6.5 |
d(H–C)b | 3.01 | 3.40 | 2.74 | 2.56 | 2.61 | 4.34 | 2.48 | |
ΔE C-b | 2.5 | 1.8 | — | 0.3 | 14.6 | 18.1 | 0.0 | 0.0 |
d(H–C)b | 2.97 | 3.48 | 2.72 | 2.73 | 3.85 | 4.36 | 5.69 | |
TS | TS1 | TS2 | TS3 | TS4 | TS5 | TS6 | TS7 | TS8 |
ΔE (Rp)-TS | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 6.8 | 18.7 |
d(H–C)c | 1.40 | 1.44 | 1.34 | 1.37 | 1.41 | 1.40 | 1.43 | 1.45 |
d(C–Li)d | 2.11 | 2.40 | 2.63 | 2.30 | 2.75 | 2.68 | 2.56 | 3.11 |
ΔE (Sp)-TS | 9.4 | 10.9 | 18.0 | 0.8 | 7.4 | 19.3 | 0.0 | 0.0 |
d(H–C)c | 1.43 | 1.44 | 1.35 | 1.36 | 1.41 | 1.39 | 1.45 | 1.43 |
d(C–Li)d | 2.14 | 2.38 | 2.60 | 2.26 | 2.77 | 2.64 | 2.90 | 3.45 |
Int | Int1 | Int2 | Int2e | Int3 | Int4 | Int5 | Int6 | Int7 |
ΔE (Rp)-Int | 0.0 | 0.0 | 0.00 | 0.7 | 0.0 | 0.0 | 14.7 | 10.6 |
ΔE (Sp)-Int | 3.1 | 13.3 | 13.3 | 0.0 | 1.2 | 1.3 | 0.0 | 0.0 |
The NCI plot software was use for visualization of non-covalent repulsions in transition states.
Footnote |
† Electronic supplementary information (ESI) available: Energies and Cartesian coordinates of all calculated structures, imaginary frequencies and NCI plots for transition states, table of LDA conformers, and potential energy scan for lithiations of amine (R)-1. See DOI: 10.1039/c3ob41339j |
This journal is © The Royal Society of Chemistry 2014 |