Effets non linéaires en catalyse asymétrique

Un effet non linéaire positif, (+)-NLE, est présent dans une réaction asymétrique qui démontre un ee du produit (ee_Produit) plus élevé que prévu par rapport à une situation linéaire idéale (cas A au lieu du cas B, figure 1). Il s'agira d'une amplification asymétrique. Un exemple d'effet non linéaire positif est observé dans le cas de l'époxydation de Sharpless avec le substrat géraniol[5].

Un effet non linéaire négatif (ou épuisement asymétrique) est présent lorsque le ee_produit est inférieur à celui prédit par une situation linéaire idéale (cas C au lieu du cas B, figure 1). Contrairement à un (+)-NLE, un (-)-NLE entraîne une vitesse de réaction globale plus rapide et une diminution de l'énantiosélectivité. Un exemple intéressant d'un effet (-)-NLE a été rapporté dans les oxydations asymétriques de sulfure.

Si le ligand chiral utilisé possède un excès énantiomérique ee_L (0 < ee_L < 100) pour générer un catalyseur, il est possible qu’une partie du catalyseur soit piégé sous forme d’espèces inactives (le réservoir). Soit α, la fraction de ces complexes non impliqués dans le cycle catalytique. Cette fraction α possédera un excès énantiomérique noté ee_res (_res pour réservoir). Si α ≠ 0 et ee_res ≠ ee_L, cela modifiera l’excès énantiomérique du ligand engagé dans le cycle catalytique. Le ee_L sera modifié, il donnera un excès énantiomérique efficace noté ee_eff. La figure 2 résume la situation.

Figure 2. Mécanisme de l'effet réservoir. Une partie du ligand engagé dans la réaction ayant un excès énantiomérique ee_L est piégé dans des espèces inactives ayant pour excès énantiomérique ee_res. Cela modifie l'excès énantiomérique de l'espèce active (ee_eff). α représente la fraction de ligand piégée dans le réservoir.

La relation entre ee_eff et ee_L sera donnée par l’équation suivante :

${\displaystyle ee_{\text{eff}}={\frac {ee_{\text{L}}-\alpha ee_{\text{res}}}{1-\alpha }}}$ et ${\displaystyle ee_{\text{produit}}=ee_{\text{0}}\ ee_{\text{eff}}}$.

Un exemple de représentation graphique ee_eff = f(ee_L) est donné sur la figure 3 pour α = 0, α = 0,5 et α = 0,75 en considérant que le réservoir piège une espèce racémique (c’est-à-dire ee_res = 0 %). Une partie du ligand sous sa forme racémique étant piégé dans le réservoir, l’excès énantiomérique de l’espèce active (ee_eff) sera augmentée par rapport à sa valeur de départ (ee_L). Il y aura une amplification de la chiralité qui sera augmentée avec la fraction du réservoir α.

Figure 3. Simulation ee_eff = f(ee_L) avec ee_res = 0 % et α = 0, 0,5 ou 0,75 (modèle de l'effet réservoir).

Ce modèle simple ne prenant pas en compte les paramètres cinétiques de la réaction permet de comprendre comment une amplification asymétrique est possible par piégeage d'une partie des espèces catalytiques actives. Le modèle peut être appliqué à des résultats expérimentaux, la figure 4 superpose le modèle « effet réservoir » avec les résultats obtenus pour la réaction d'alkylation de ZnEt₂ sur le benzaldéhyde en présence du ligand chiral DAIB[6]. Ce modèle reste cependant très simple et ne peut donner que des droites comme simulation ce qui est éloigné de la réalité expérimentale. Il ne prend pas en compte l'aspect cinétique des réactions et la fraction d'espèce inactive α est constante quelle que soit la valeur de ee_L, ce qui ne correspond pas à la réalité des systèmes étudiés.

Figure 4. Simulation des données expérimentales de la réaction d'alkylation du benzaldéhyde en présence du ligand DAIB avec le modèle « effet réservoir ». La simulation a été obtenue avec ee₀ = 98 %, ee_res = 0 % et α = 0,88.

Les modèles communément appelés « ML_n » impliquent la présence d’espèces catalytiques qui associent un métal à plusieurs ligands : deux ligands = modèle ML₂, trois ligands = modèle ML₃, etc.[1]. Dès que l’on associe deux ligands (ou plus) sur un même métal, plusieurs diastéréomères sont possibles. Toutes ces espèces peuvent avoir une concentration et une réactivité différente et donc induire des perturbations. Ces modèles prennent en compte la cinétique de la réaction. Ces modèles s'appliquent également lorsque l'on a des oligomères : dimère (ML)₂, trimère (ML)₃, etc.

Le cas le plus simple serait l'assemblage de deux ligands chiraux avec un métal pour former un complexe, ce qui donne le modèle ML₂. Dans un tel système, trois stéréoisomères sont possibles comme le montre la figure 5. Deux complexes homochiraux (ML_RL_R et ML_SL_S, qui donneront les produits avec les excès énantiomériques ee₂ et -ee₂, respectivement) et un complexe hétérochiral (ML_RL_S) (souvent appelé « composé méso ») qui donnera le produit racémique (ee = 0 %). Le produit de la catalyse est de ce fait obtenu à partir des trois catalyseurs en équilibre et qui fonctionnent simultanément. Par conséquent, l'excès énantiomérique du produit (ee_P) dépendra de (i) l’ee du ligand (ee_L), (ii) de la réactivité relative g, avec k₂ et k₂' étant les constantes cinétiques correspondantes et de (iii) la concentration relative en catalyseur β (avec [RR], [SS] et [RS] définis comme les concentrations individuelles des espèces ML_RL_R, ML_SL_S et ML_RL_S. K est la constante d'équilibre entre les espèces homochirales et hétérochirales.

Figure 5. Schéma général pour le modèle ML₂ qui permet de simuler les courbes ee_P = f(ee_L). L'ee_P global dépendra de l'ee du ligand (ee_L), de l'énantiosélectivité du dimère homochiral (ee₂), des quantités relatives des catalyseurs présents ([RR], [SS] et [RS], exprimées par la constante d'équilibre K et de leur réactivité relative g.

En prenant en compte tous ces paramètres, il est possible de relier ee_P à ee_L :

${\displaystyle ee_{\text{P}}=ee_{\text{2}}\ ee_{\text{L}}\ {\frac {1+\beta }{1+g\beta }}}$

avec ${\displaystyle \beta ={\frac {-K{ee_{\text{L}}}^{2}+{\sqrt {-4{ee_{\text{L}}}^{2}+K(4+K{ee_{\text{L}}}^{2})}}\ }{4+K{ee_{\text{L}}}^{2}}}}$.

• Si β = 0, aucun complexe hétérochiral (ML_RL_S) n'est présent ; l'équation ML₂ se simplifie. Idem si g = 1, les complexes auront la même réactivité. Par conséquent, dans ces deux cas, les propriétés additives simples devraient s'appliquer pour établir une relation linéaire entre l'énantiosélectivité du produit et l'énantiopureté du ligand chiral.

• Si l'on a une distribution statistique des trois espèces chirales, alors la constante d'équilibre sera égale à 4. La figure 6 décrit des simulations de courbes à différentes valeurs de g. Si g = 0 (c'est-à-dire que l'espèce hétérochirale n'est pas active), un effet (+)-NLE est observé. Si g > 1 alors l'espèce hétérochirale est plus active que l'espèce homochirale, un effet (-)-NLE est observé.

Figure 6. Modélisation de courbes pour le modèle ML₂ avec K = 4 (c'est-à-dire une distribution statistique des trois espèces) et pour différentes valeurs de g (0, 0,01, 0,1, 0,33, 1, 3, 10 et 100).

• Une situation devient intéressante lorsque la constante K est élevée, autrement dit lorsqu'il y a une prédominance du composé hétérochiral. La figure 7 représente des exemples à différentes valeurs de g pour K = 2 500. Si g = 0 (espèce hétérochirale inactive), l'effet non linéaire positif devient extrêmement important. Cette situation rejoint celle de l'effet réservoir. L'espèce hétérochirale (racémique) joue le rôle de réservoir et induit un fort effet (+)-NLE. Inversement, si g est grand (par exemple g = 100), l'espèce hétérochirale prédomine et comme elle génère le produit sous sa forme racémique, une forte déplétion est observée (effet non linéaire négatif).

Figure 7. Modélisation de courbes pour le modèle ML₂ avec K = 2 500 et pour différentes valeurs de g (0, 0,01, 0,1, 0,33, 1, 3, 10 et 100).

Sur le même schéma, le modèle peut être étendu à des systèmes plus complexes, par exemple lorsqu'un complexe métallique avec 3 ligands donnant une espèce de type ML₃ est formé, ou une espèce trimérique (ML)₃. Dans le cas ML₃, quatre espèces sont possibles : (ML_RL_RL_R), (ML_SL_SL_S), (ML_RL_RL_S) et (ML_SL_SL_R) (figure 8). Comme dans le modèle ML₂, il y a coexistence d'une espèce homochirale et d'une espèce hétérochirale, mais avec une différence majeure : l'espèce hétérochirale n'est pas méso et va générer un produit de réaction qui ne sera pas racémique. Il s'agit de ce fait d'une situation où les deux espèces vont générer le produit avec un certain ee, les courbes des effets non linéaires peuvent être très caractéristiques. La forme devient indicative d'un mécanisme.

Un exemple intéressant avec le modèle ML₃ est décrit sur la figure 9. Si l'espèce hétérochirale est plus réactive que l'espèce homochirale (grande valeur de g avec g = k₃'/k₃) et qu'elle induit un produit avec une meilleure énantiosélectivité (c'est-à-dire que ee'₃ >> ee₃₎ alors, il serait possible d'avoir un catalyseur plus efficace lorsqu'il n'est pas énantiopure. Si l'on prend les courbes de la figure 9, les paramètres ont été fixés à : ee₃ = 50 % ee et ee'₃ = 100 % ee, donc l'espèce hétérochirale induit une meilleure énantiosélectivité que l'espèce homochirale. Dès que g est significatif, la courbe montre une zone où le catalyseur est plus efficace sous une forme non énantiopure.

Afin de nommer ce type d'effet non linéaire, Henri Kagan a proposé le terme d'« effet non linéaire hyperpositif »[1]. Le premier exemple expérimental d'effet non linéaire hyperpositif a été observé 26 ans après les hypothèses de Kagan, cependant le mécanisme s'est avéré différent[7]

Figure 9. Exemple de simulation pour le modèle ML₃.

Le modèle ML₄ a également été traité par Henri Kagan, il suit le même schéma que les deux modèles ci-dessus. Des modélisations de ce cas ont donné des courbes très caractéristiques.

Les effets non linéaires de la réaction de Sharpless ont été étudiés avec le géraniol comme substrat par Kagan dès 1986 (figure 10a)[5]. L'effet non linéaire observé est un effet non linéaire positif avec une légère amplification asymétrique comme le montre la figure 10b. En appliquant le modèle de type ML₂, il est possible d'avoir une bonne prédiction de la courbe avec les paramètres indiqués sur la figure. Les résultats sont cohérents avec les études menées par Barry Sharpless qui suggéra une espèce dimérique comme espèce active (donc de type (ML)₂)[8]. La valeur de g dans la simulation est de 0,35 c'est-à-dire que k'₂/k₂ ≈ 1/3 ce qui suggère que l'espèce dimère (ML)₂ homochirale (figure 10c) est trois fois plus réactive que l'espèce homochirale.

Figure 10. a) la réaction d'époxydation asymétrique de Sharpless avec géraniol comme substrat et le (R,R) diéthyltartrate comme ligand chiral. b) Effet non linéaire observé, superposition des résultats expérimentaux et du modèle ML₂. c) Structure moléculaire suggérée de l'espèce homochirale active.

La réaction asymétrique du dialkylzinc et du benzaldéhyde en présence du ligand chiral (-)-3-exo-(diméthylamino)isobornéol (nommé DAIB) présente des phénomènes non linéaires inhabituels ayant fait l'objet de nombreuses études mécanistiques. En présence d'une quantité catalytique de DAIB, la réaction du dialkylzinc et d'un aldéhyde aromatique est nettement accélérée pour donner, après hydrolyse, l'alcool correspondant avec une très grande pureté énantiomérique (figure 11a). L'étude de l'excès énantiomérique du produit en fonction de l'excès énantiomérique du ligand [ee_P = f(ee_L)] donne un effet non linéaire positif très marqué comme le montre la figure 11b. Par exemple pour R = Et, un excès énantiomérique de 95 % pour le produit est obtenu avec un ligand ayant uniquement 15 % d'excès énantiomérique.

Figure 11. a) La réaction d'alkylation du benzaldéhyde, catalysée en présence du ligand DAIB. b) Un effet non linéaire positif (+)-NLE est observé avec ZnEt₂ ou ZnMe₂ comme réactif.

Le mécanisme de la réaction a été étudié par Ryoji Noyori et l'origine de l'effet non linéaire a pu être déterminé[6] - [9]. Celui-ci dérive du modèle ML₂ de Kagan. Ici, les (+)-NLE proviennent de la formation de dimères homochiraux et hétérochiraux, qui jouent le rôle de réservoir. Contrairement au modèle ML₂ de Kagan, la réaction est catalysée par des complexes monomériques tandis que les dimères sont inactifs et en équilibre dynamique avec les monomères. Ici, le NLE est régie purement par les stabilités relatives des espèces dimériques, qui à leur tour modifient le rapport R/S des monomères catalytiquement actifs. Un (+)-NLE se produit lorsque le dimère hétérochiral est plus stable que le dimère homochiral (et inversement). La figure 12 donne un aperçu global du mécanisme.

Figure 12. Mécanisme expliquant l'effet non linéaire de la figure 11. L'ee_P global est une conséquence du rapport R/S des catalyseurs monomériques actifs (les dimères sont inactifs) qui dépend de l'ee du ligand (ee_L), de l'énantiosélectivité du catalyseur et des constantes de dimérisation pour former les dimères. Si l'espèce dimérique hétérochirale est plus stable, alors un effet non linéaire positif sera observé.