J. Org. Chem. 1999, 64, 4725-4732
4725
Theoretical Investigations of the Influence of Pressure on the Selectivity of the Michael Addition of Diphenylmethaneamine to Stereogenic Crotonates F. Dumas,†,‡ C. Fressigne´,§ J. Langlet,† and C. Giessner-Prettre*,† Laboratoire de Chimie The´ orique, UMR 7616 CNRS, Universite´ P. & M. Curie, Case Courrier 137, 4, place Jussieu, 75252 Paris Cedex 05, France, and Centre Re´ gional Universitaire de Spectroscopie UPRESA 6014 CNRS, IRCOF, Universite´ de Rouen, 76821 Mont St Aignan Cedex, France Received December 30, 1998
SIBFA (sum of interactions between fragments ab initio computed) molecular mechanics systematics has been modified to take into account the effect of pressure on intra- and intermolecular energies. The van der Waals radii are related to the pressure, using Bridgman results on the variation of crystal volume, on one hand, and the relation between the volume of an atom and its van der Waals radius on the other. This procedure produces a decrease of the volume of the systems considered. The modified systematics is used for the study of the conformation at 1 atm and 15 kbar of two stereogenic crotonates and of the complexes formed by these two molecules with the diphenylmethaneamine and the three solvent molecules present in the experiment. The results obtained show that in the case of NMECC 1a the diastereoselectivity induced by high pressure and by the presence of methanol proceeds from an important stabilization of the pro-R reactive complex in which the crotonate has a stacked-transoid conformation. This stabilization is mainly due to intermolecular interactions. In the case of the second crotonate considered, NMCC 1b, our results indicate that pressure induces a stabilization of the pro-R and pro-S complexes having the axial conformation for which the reaction exhibits little diastereoselectivity in qualitative agreement with experimental data. This study tends to show that it is possible to account theoretically for the influence of pressure on molecular conformation and/or complex structure, using a molecular mechanics method that is able to take into account the variation of volumes of the different entities present in the system studied. Introduction A previous molecular mechanics study1 has found a close relationship between the theoretical stability of the stacked conformation of 8-arylmenthol crotonates 1 and the efficiency as a chiral auxiliary of the corresponding alcohols 3c for the diastereoselective addition of diphenylmethaneamine, leading to β-aminoesters (Chart 1), under high pressure (15 kbar).2 These calculations have even been able to propose a new chiral auxiliary, 8-omethoxyphenylmenthol (3d), the efficiency of which has been experimentally confirmed.3 However, recent experimental results suggest that the agreement between theory and experiment could seem partially fortuitous and require a need to take into account some additional features to be fulfilled in all cases.4 Thus, the computations have considered only the transoid conformation (Chart 2) of the crotonate moiety since the stereochemical outcome of this reaction is consistent with the addition of the amine to the transoid conformation of the crotonate reactant.2,3 Quite unexpectedly, the X-ray study of the
Chart 1
Chart 2
†
Universite´ P. &. M. Curie. Permanent address: Laboratoire de Chimie Organique, UPRESA 8076 CNRS Universite´ de Paris Sud, Centre d’Etudes Pharmaceutiques, 5, rue J.-B. Cle´ment, 92296 Chatenay-Malabry, France. § Universite ´ de Rouen. (1) Maddaluno, J. F.; Gresh, N.; Giessner-Prettre, C. J. Org. Chem. 1994, 59, 793. (2) d’Angelo, J.; Maddaluno, J. J. Am. Chem. Soc. 1986, 108, 8112. (3) Dumas, F.; Mezrhab, B.; d’Angelo, J.; Riche, C.; Chiaroni, A. J. Org. Chem. 1996, 61, 2293. (4) Mezrhab, B.; Dumas, F.; d’Angelo, J.; Riche, C. J. Org. Chem. 1994, 59, 500. ‡
related trans-2-[1-(2-naphthyl)-1-methylethyl]cyclohexyl(E)-crotonate (NMECC, 1a) has shown that, in the crystal, the crotonate moiety adopts a cisoid conformation.4 More recent computations carried out on four chiral cyclohexyl-based crotonates bearing a β-naphthyl group have shown that the energy of their cisoid and transoid
10.1021/jo9825308 CCC: $18.00 © 1999 American Chemical Society Published on Web 05/29/1999
4726 J. Org. Chem., Vol. 64, No. 13, 1999
Dumas et al.
Table 1. Diastereoselectivities in the Michael Addition of Ph2CHNH2 to NMEC- and NMC-based Chiral Crotonates as a Function of Experimental Conditions (values from Ref 3) de in the Michael addition of Ph2CHNH2 to chiral crotonates derived from alcohols:
conditions MeOH, 15 kbar MeOH, 1 atm CH3CN-Et2O, 15 kbar CH3CN-Et2O, 1 atm a
98 10 8 a
18 7 5 a
No reaction in these conditions.
conformations differs by 0.1-0.6 kcal/mol.5 These values are not large enough to justify some of the selectivities obtained experimentally (Table 1).3 Furthermore, the diastereoisomeric excess (de) appears to be very sensitive to the solvent; it decreases considerably when going from CH3OH to CH3CN/Et2O.2,3 Finally, under atmospheric pressure, the reaction is not only considerably slower even in CH3OH (5 days at 40 °C) but also much less diastereoselective as can be seen from the values reported in Table 1.3 This set of data tends to show that the stability of the stacked conformers with respect to the other possible conformations (Chart 2) of such chiral crotonates is only one (even if it is the leading as suggested by previous results)1,3 of the factors at work in such high-pressure diastereoselective Michael additions. Because of the importance of activation volumes in high-pressure chemistry6 and the pressure dependence of the stereoselectivity that is assumed to originate from a difference in the activation volume of the possible diastereomeric transition states, the volumes of the different conformers of chiral crotonates under experimental study have been calculated. The results obtained show, provided that under high-pressure conditions the conformers with the smallest volumes and relative energies up to 1-2 kcal/ mol are favored, a good correlation between the conformational molecular volume and their π-facial discrimation in the high-pressure-induced Michael addition of diphenylmethaneamine.5 However, if these results can explain the influence of the pressure on the de, they do not provide any explanation for the role of the solvent on the selectivity. At this point, we underline that the role of the solvent on the kinetics of such a reaction is not considered in the present work since previous studies have provided satisfactory accounts of the contribution of protic solvents to related Michael additions.7 (5) Dumas, F.; Fressigne´, C.; Langlet, J.; Giessner-Prettre, C. In Chemical processes and reactions under extreme or nonclassical conditions, COST action D6, Luche, J. L., Balny, C., Benefice, S., Denis, J.-M., Pe´trier, C., Eds.; Office for Official Publications of the European Communities: Luxembourg, 1998; p 71. (6) (a) Isaacs, N. S. In Liquid-Phase High-Pressure Chemistry; John Wiley and Sons: Chichester, New York, Brisbane, Toronto, 1981; pp 182-201. (b) Van Eldik, R.; Asano, T.; Le Noble, W. J. Chem. Rev. 1989, 89, 549. (c) Tietze, L. F.; Ott, C.; Gerke, K.; Buback, M. Angew. Chem., Int. Ed. Engl. 1993, 32, 1485. (d) Jenner, G. Tetrahedron 1997, 53, 2669. (e) Kla¨rner, F. G.; Ruster, V.; Zimny, B.; Hochstrate, D. HighPressure Res. 1991, 7, 133. (7) (a) Abbot, D. J.; Colonna, S.; Stirling, C. J. M. J. Chem. Soc., Perkin Trans. 1 1976, 492. (b) Pardo, L.; Osman, R.; Weinstein, H.; Rabinowitz, J. R. J. Am. Chem. Soc. 1993, 115, 8263.
To get some insight into the different factors that could contribute to the diastereoselectivity of this particular reaction, we undertook a molecular mechanics study using a version of the SIBFA (sum of interactions between fragments ab initio computed) method8 modified to take into account pressure. Previous studies concerned by modeling molecular conformation/interactions under high pressure have been carried out using Monte Carlo9 or molecular dynamics10 methods. Because of the size of the systems to be studied (98 atoms in the largest ones), we could not consider the use of an ab initio method.11 Furthermore, the number of conformational domains (12), depicted in Chart 2, to be explored for the different situations to be considered would require computing facilities beyond thoses currently in use in organic chemistry modeling if molecular dynamics computations were undertaken. On the other hand, SIBFA, which has given good results not only for the correlation stacked conformer energy/diastereoselectivity of the high-pressure-induced Michael addition in these systems1 but also for precursors of the transition states of an intramolecular Diels-Alder reaction,12 can be expected to give a reasonable first estimate of the influence of pressure on the arrangement of the reactants as well as their relative energies. The calculations will be carried out at 15 kbar, as in the experiments, and at atmospheric pressure. The largest systems considered will include a β-naphthylmethylethylcyclohexyl crotonate (NMECC, 1a), the diphenylmethaneamine, and the number (three) of methanol molecules corresponding to the experimental conditions. The search of the most stable geometrical arrangement of such complexes, carried out under both conditions, will provide indications on the conformational changes induced by high pressure. These computations will unfortunately be dealing not with the transition state itself but only with what can be considered as the “pretransition state” (PTS) or the “near attack conformation ” (NAC) as proposed by Lightstone and Bruice.13 But under high pressure, the possibilities of conformational equilibrium and of diffusion are considerably reduced; the stability of the PTS should, therefore, be able to reflect to some extent that of the transition state. These computations are carried out for two auxiliaries, trans2-[1-(2-naphthyl)-1-methylethyl]cyclohexanol (NMEC, 3a) and trans-2-(2-naphthylmethyl)cyclohexanol (NMC, 3b), which give significantly different π-facial differentiation of the corresponding crotonates (Table 1). A delineation of the role of the different forces at work in the reaction will be obtained from the variation, with the pressure, of the intra-/intermolecular energies (and of their different contributions). Furthermore, similar calculations carried out on the isolated chiral crotonates on one hand and on a model of the “reacting” system on the other will (8) Gresh, N. J. Chim. Phys. 1997, 94, 1365. (9) Jorgensen, W. L. J. Am. Chem. Soc. 1981, 103, 421. (10) (a) Root, L. J.; Berne, B. J. J. Chem. Phys. 1997, 107, 4350. (b) Paci, E.; Marchi, M. J. Phys. Chem. 1996, 100, 4314. (c) Kowall, T.; Foglia, F.; Hlem, L.; Merbach, A. E. Chem. Eur. J. 1996, 2, 285. (d) Ferrario, M.; Rychaert, J. P. Mol. Phys. 1985, 54, 587. (e) Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, A.; Haak, J. R. J. Chem. Phys. 1984, 81, 3684. (11) (a) Remler, D. K.; Madden, P. A. Mol. Phys. 1990, 70, 921. (b) Galli, G.; Parinello, M. In Computer Simulations in Material Science; Meyer, M., Pontikis, V., Eds.; Klower: Dodrecht, 1991; p 1. (12) Giessner-Prettre, C.; Hu¨ckel, S.; Maddaluno, J.; Jung, M. E. J. Org. Chem. 1997, 62, 1439. (13) Lightstone, F. C.; Bruice, T. J. Am. Chem. Soc. 1994, 116, 10789.
Pressure on the Selectivity of a Michael Addition
provide information of the role of the methanols in the reaction through the hydrogen bonding chart found in the complex. This analysis is possible since molecular mechanics computations give separately the inter- and intramolecular contributions to the energy as well as their different components (electrostatic, dispersion, etc.). The first part of this study is concerned by the cisoidtransoid energy difference in methyl crotonate, a small enough molecule to be studied by sophisticated methods. Such calculations are necessary to validate the corresponding molecular mechanics results. This quantity is of prime importance for the problem studied since no de is expected for isoenergetic cisoid and transoid conformers.
Methodology With SIBFA systematics, the total energy is calculated as the sum of inter- and intramolecular contributions
E ) EInter +
∑E
J. Org. Chem., Vol. 64, No. 13, 1999 4727 intermolecular energy within the SIBFA framework. Isaacs15 has suggested that the compressibility (molecular volume/ interatomic distance decrease with pressure) proceeds from a reduction of the “effective” radius of atoms in solids and of the “void” between molecules in liquids. Following this suggestion, we have introduced the influence of pressure in the inter-/ intramolecular energies by varying the van der Waals radii, using the above relations and Bridgman results. Bridgman16 has determined, from measurements of volume change, due to pressure P, of isotropic crystals a pressure/volume relationship
-∆V ) aP - bP2 V0 it is possible to take into account the pressure effects on the van der Waals radii, which are the “effective” radii of the atoms in the crystal. Following the definition of the compressibility β at temperature T
1 ∂V β)V0 ∂P
( )
T
Intra
and taking into account that
the summation running the different molecules of the system considered. Einter is calculated as the sum of 5 terms8
∂∆V )( (∂V ∂P) ∂P ) T
EInter ) EMTP + ERep + EDisp + EPol + ECT In the case of EIntra, there is no charge-transfer term but a torsion contribution, ETor, is taken into account, and we have
EIntra ) EMTP + ERep + EDisp + EPol + ETor
it appears that the coefficients a and b are related to the compressibility by
β ) a - 2bP and
(∂P∂β)
Because EMTP, ERep, and EDisp are pair-wise additive terms, E can also be written
E)
∑E
ij
+ EPol + ECT + ETor
i,j
where Eij is the interaction energy between nonbonded atoms i and j belonging to fragments F and F′ (belonging to the same or to two different molecules of the system studied), respectively. The van der Waals contribution to Eij is given by
Eijvw ) EDisp + ERep ) Aij/Rij6 + B exp(-RijRij) where Rij is the interatomic distance and B a constant, while Aij and Rij have values that depend on the atomic species of i and j. Following Kitaygorodski,14 Aij is taken as:
Aij ) A′ro6ij where A′ is a universal constant and roij, corresponds to the best interatomic distance measured in reference crystals, such as graphite for carbon. Atoms i and j van der Waals radii, rvwi and rvwj, are deduced from roij using relation
roij ) rvwi + rvwj Similarly for Rij, we can write
Rij ) R′/roij The calculation of the van der Waals contribution (EDisp + ERep) is carried out following Kitaygorodski’s chart, which relates this quantity to equilibrium interatomic distances in crystals. This expression for EDisp and ERep is of particular adequacy for taking into account the effect of pressure on intra- and (14) Kitaygorodsky, A. I. Tetrahedron 1961, 14, 230.
T
) -2b
T
On the other hand, each atom can be considered as a sphere with a radius equal to its van der Waals radius. Thus, at atmospheric pressure its volume is equal to
Vo )
4 πr 3 3 o vw
and the volume variation with pressure is given by
∆V ) Vp - Vo )
4π 3 - ro3vw) (r 3 p vw
Therefore, the relation 3
3 vw
∆V rp vw - ro ) Vo ro3vw
)
rp3vw ro3vw
-1
gives the van der Waals radius rpvw at pressure P provided the volume variation ∆V be known for that pressure. Bridgman has tabulated for a number of atoms either the coefficients a and b which relate ∆V to the pressure or ∆V at different pressures. In this last case, the a and b coefficients are obtained from a least-squares treatment.16 As precedingly,5,17 the computation of the molecular volume Ω is based on van der Waals volume and not on electron density envelopes.18 With that approximation, it can be obtained using Ostrogradski’s formula from a surface integration where Σ is the volume Ω integration boundary, that is (15) Isaacs, N. S. Tetrahedron 1991, 40, 8463. See also: Matsumoto, K.; Acheson, R. M. In Organic Synthesis at High Pressure; Matsumoto, K., Acheson, R. M., Eds.; John Wiley: New York, 1991; pp 3-73. (16) Bridgman, P. W. In The Physics of High Pressure; G. Bell: London, U.K., 1949. (17) Huron, M.-J.; Claverie, P. J. Phys. Chem. 1972, 76, 2123. (18) Wong, M. W.; Wiberg, K. B.; Frisch, M. J. J. Comput. Chem. 1995, 16, 385.
4728 J. Org. Chem., Vol. 64, No. 13, 1999 Chart 3
Dumas et al. Table 2. Conformational Energy Minima (kcal/mol) and Conformational Molecular Volume (mL) of NMEC-based Chiral Crotonatea 1 atm
Chart 4
15 kbar
conf.
E
∆E
V
∆V
E
∆E
V
∆V
Sc St Tc Tt Ac At
606.4 606.6 607.4 607.6 610.0 610.3
0.0 0.2 1.0 1.2 3.6 3.9
190.7 186.6 194.1 194.3 192.4 191.8
4.1 0.0 7.5 7.7 5.8 5.2
412.7 413.0 414.9 415.2 413.6 413.8
0.0 0.3 2.2 2.5 0.9 1.1
161.9 159.5 168.4 168.0 166.3 166.1
2.4 0.0 8.9 8.5 6.8 6.6
a Under its threshold “Stacked” (S), “Trans” (T) and “Axial” (A) (position of the aryl group), as well as transoid (t) and cisoid (c) [only the most stable of the two orientations of the naphthyl group (1 or 2, Chart 2) are given] conformation.
role of methanol might be interpreted in terms of hydrogenbond formation between this solvent and the crotonate moiety, a phenomenon favored under high-pressure, because of the decrease in atomic distance and suppression of thermal motions.3 The ab initio calculations have been carried out using the Gaussian9419 series of programs. The geometry optimizations were performed with the 6-31G** basis and the B3P8620 functional.
Chart 5
Results and Discussion
the molecular surface and b ne is the normal vector to Σ directed toward the outside of Ω.
Ω)
∫ ∫∫ dv ) ∫ 31 br bn (Ω) dσ Ω
Σ
e
The integration surface Σ is obtained from the summation of (parts of) spheres centered on each atom i of the molecule (Chart 3). For atoms i and j, if the interatomic distance Rij is larger than the sum of their van der Waals radii, the surface of both spheres is included in Σ. If the interatomic distance Rij is smaller than the sum of the van der Waals radii, part of the spheres have to be skipped. The surface element dσi is kept if and only if the distance between dσi and atom j is larger than atom j’s van der Waals radius by running over sphere i using a Korobov grid.17 Points dσj of sphere j are kept according to the same criteria. The values of the van der Waals radii used for the volumes computations take into account the pressure. Therefore, it is possible to have an estimation of the pressure on this quantity for the different possible conformers. Other procedures have been proposed to calculate molecular volumes.18 The one retained herein presents the advantage of being fully consistent with the method used to take into account the effect of the pressure. The fragments (Chart 4) used as input for the chiral crotonates are the same as those discussed in ref 1 except for the crotonate entity. In the present work, the crotonaldehyde is made of three fragments: formaldehyde, ethylene, and methane in order to allow the cisoid-transoid conformational variation. Diphenylmethaneamine is made of ammonia, methane, and two benzene fragments, while water and methane are those of the methanols. The multipoles (monopoles, dipoles, and quadrupoles) are calculated from the MP2 density matrixes in order to take into account the contribution of electron correlation to the molecular electron distribution. Optimizations for the “reacting” systems were carried out through the hydrogen bonding chart depicted in Chart 5 since the critical
Methyl Crotonate. The ab initio complete optimization of the cisoid and transoid conformations, using the B3P86 hybrid functional, gives a 0.7 kcal/mol energy difference in favor of the cisoid form. This value is close to those reported by Houk et al.21 on one hand and RuizLopez et al.22 on the other for methyl acrylate at the SCF level (0-1.8 kcal/mol depending on the basis set and computational level used). SIBFA gives a cisoid conformer more stable than the transoid one by 0.4 kcal/ mol. This result shows that SIBFA is not only able to give the correct order of stabilities but also reasonable values of energy differences between quasi-isoenergetic conformations. This systematics should, therefore, account correctly, in the case of the compounds under study, for the energy variation with this particular conformational parameter. trans-2-[1-(2-Naphthyl)-1-methylethyl]cyclohexyl(E)-crotonate (NMECC) 1a. This compound corresponds to the case of an excellent de (98%) when the addition takes place at 15 kbar in methanol and a poor one in other conditions as can be seen from Table 1. The values of Table 2, which concern the isolated NMECderived chiral crotonate at 1 atm, and Figure 1a show that the SIBFA computations give results in excellent agreement with experiment since the most stable conformation is the stacked cisoid (Sc) form with orientation 1 (Chart 2) of the naphthyl group as observed in the crystal.4 We see also from Table 2 that the energy (19) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Gill, P. M. W.; Johnson, B. G.; Robb, M. A.; Cheeseman, J. R.; Keith, T. A.; Petterson, G. A.; Montgomery, J. A.; Raghavachari, K.; Al-Laham, M. A.; Zarzewski, V. G.; Ortiz, J. V.; Foresman, J. B.; Cioslowski, J.; Stefanov, B. R.; Nanayakkara, A.; Challacombe, M.; Peng, C. Y.; Ayala, P. Y.; Chen, W.; Wong, M. W.; Andres, J. L.; Replogle, E. S.; Gomperts, R.; Martin, R. L.; Fox, D. J.; Binkley, J. S.; Defrees, D. J.; Baker, J.; Stewart, J. P.; Head-Gordon, M.; Gonzales, C.; Pople, J. A. GAUSSIAN 94 (Revision B.1); Gaussian Inc.: Pittsburgh, PA, 1995. (20) (a) Becke, A. D. J. Chem. Phys. 1993, 98, 5648. (b) Perdew, J. P. Phys. Rev. B 1986, 33, 8822. (21) Loncharich, R. J.; Schwartz, T. R.; Houk, K. N. J. Am. Chem. Soc. 1987, 109, 14. (22) Ruiz-Lopez, M. F.; Assfeld, X.; Garcia, J. I.; Mayoral, J. A.; Salvatella, L. J. Am. Chem. Soc. 1993, 115, 8780.
Pressure on the Selectivity of a Michael Addition
J. Org. Chem., Vol. 64, No. 13, 1999 4729 Table 3. Conformational Energy Minima (kcal/mol) and Conformational Molecular Volume (mL) of the Reactive Complexa 1 atm (de ) 18%)
15 kbar (de ) 98%)
conf.
E
∆E
V
∆V
E
∆E
V
∆V
Sc St Tc Tt Ac At
1254.0 1253.6 1256.1 1254.8 1258.7 1256.6
0.4 0.0 2.5 1.2 5.1 3.0
382.2 376.3 386.4 383.6 379.8 381.0
5.9 0.0 10.1 7.3 3.5 4.7
864.8 862.3 868.5 865.0 865.6 864.1
2.5 0.0 6.2 2.3 3.3 1.7
314.7 314.6 325.9 326.9 322.1 321.6
0.1 0.0 11.3 12.3 7.5 7.0
a Including the NMEC based chiral crotonate under relevant conformations (see Table 2), diphenylmethaneamine, and three methanol molecules.
Figure 1. Most stable geometrical arrangement of the cisoid form of NMECC at 1 atm (A) and under 15 kbar (B).
difference between the cisoid and transoid conformations is very small (0.2-0.3 kcal/mol) and of the same order of magnitude for the stacked, trans, and axial conformers. This feature, which is not modified by a pressure of 15 kbar as can be seen from the corresponding values of Table 2, is in agreement with the low de obtained in absence of protic solvent since the stacked cisoid (Sc) and transoid (St) conformers, the most stable ones, give opposite diastereoisomers. The values of Table 2 show also that pressure decreases the difference of volume between the cisoid and transoid forms of the stacked conformers, the St conformer remaining, however, by far the more compact. In addition, pressure inverts the stability of the trans and axial conformations. At 15 kbar, the latters have an energy that differs by less than 1 kcal/mol from that of the stacked conformations. This variation, to be related to the molecular volumes at 1 atm of the different conformers (the axial being more compact than the trans), confirms our hypothesis on the importance of this factor for their stabilization under high pressure since the volume of the axial arrangement is more compact than the trans one. From Figure 1, which shows the optimized geometries of the Sc conformation at both pressures, we see that the only significant variation concerns the orientation of the crotonate moiety that rotates about the C-O bond, getting closer to the naphthyl rings upon pressure increase. The same variation takes place for the St conformer (not shown). When the “complex”, including the diphenylmethaneamine and three methanols in addition to NMECC 1a (Chart 5) which can reasonably be considered as the entity that undergoes the chemical reaction, is optimized, the results obtained at 1 atm and 15 kbar differ significantly. From Table 3, we see that under both conditions the St conformation of NMECC 1a is the most stable. However, the tabulated values show clearly that at 1 atm the energy difference between the cisoid and the transoid forms (0.4 kcal/mol) is not large enough to give rise to a large de. Under a pressure of 15 kbar, we find that the St conformation is more stable by 2.5 and 1.7 kcal/mol than the Sc and axial-transoid (At) conformers, a result able to account correctly for the large de (98%), which is obtained in these conditions.3 This last value, which shows that the energy difference between the stacked and At conformers is considerably reduced by the pressure (1.7 kcal/mol instead of 3.0 kcal/mol), suggests that under high pressure the second diastereoisomer formed would not be due to the cisoid/transoid equilibrium but to the
Figure 2. Most stable geometrical arrangements of the cisoid (A & B) and transoid (C & D) of the complex formed by NMECC, three methanols and the diphenylmethaneamine at 1 atm (A & C) and under 15 kbar (B & D).
stacked/axial one. The values of Table 2 concerning the volume of the complex (Figure 2) show that the geometrical arrangement of the “complex” is significantly less modified when going from 1 atm to 15 kbar in the case of the Sc conformer than in the case of the St form. The examination of the different contributions to the energy of the complex for the Sc and St conformations of NMECC 1a shows that the intramolecular energy is in favor of the cisoid form (1279.4/1279.8 kcal/mol at 1 atm and 901.3/905.1 (Table 4) under 15 kbar), while the intermolecular contribution is larger for the transoid conformation of the crotonate (-25.3/-26.3 kcal/mol at 1 atm and -35.7/-42.8 kcal/mol under 15 kbar). From the values of Table 4, we see that the larger value of EInter in the system with the St crotonate is mainly due to the electrostatic forces (EMTP), while the van der Waals contribution (EDisp + ERep) is in favor of the Sc form. From the values of the interaction energies (and of their contributions) of the different “dimers” of the complex, we see that the interactions of NMECC 1a with the methanol 3 (hydrogen bonded to the carbonyl) on one hand and methanol 4 with the amine on the other are responsible for the greater stability of the complex with St conformation (see Chart 5 for the molecules numbering). The values of Table 4 show also that the interaction
4730 J. Org. Chem., Vol. 64, No. 13, 1999
Dumas et al.
Table 4. Contributions (kcal/mol) to the Intermolecular Interaction Energy in the “Pretransition State” Complex and Between the Different Possible “Dimers” Presentb intermolecular EDisp
ERep
-19.3 41.4
Sc 1-2 1-3 1-4 1-5 2-3 2-4 2-5 3-4 3-5 4-5 St 1-2 1-3 1-4 1-5 2-3 2-4 2-5 3-4 3-5 4-5
EPol
intramolecular
ECT
EMTP
EInter
EaIntra
Etot
-9.9 -5.6 -42.5 -35.7
901.3
864.8
-1.4 -2.2 -0.2 -0.1 0.0 0.0 -0.4 -3.2 -0.1 -0.4
-0.8 -8.1 -10.7 -2.1 -9.1 -8.4 -0.0 -0.5 -1.6 -0.1 -0.9 -2.3 0.0 0.1 0.1 0.0 0.1 0.1 -0.2 -4.4 -4.4 -2.3 -16.5 -15.4 0.0 -1.0 -1.7 -0.1 -2.3 -2.6
539.6 539.6 539.7 521.7 230.8 230.9 212.9 230.9 212.9 212.9
528.9 531.2 538.1 519.4 230.9 231.0 208.5 215.5 211.2 210.3
-24.1 56.7 -14.2 -7.9 -53.5 -42.8
905.1
862.3
540.0 540.0 540.0 522.6 231.0 230.9 215.4 230.9 215.5 215.4
530.1 528.6 539.0 519.4 231.2 231.0 213.2 215.0 212.9 207.5
-3.5 3.1 -4.5 9.4 -1.0 0.2 -3.2 1.9 0.0 0.0 0.0 0.0 -1.1 1.9 -4.1 10.5 -1.0 0.4 -0.9 1.0
-4.3 3.6 -7.0 18.8 -0.5 0.0 -4.2 2.5 0.0 0.0 0.0 0.0 -0.7 0.5 -4.4 11.8 -0.7 0.2 -2.3 5.3
-1.3 -4.2 -0.1 -0.4 0.0 0.0 -0.3 -3.5 -0.1 -1.4
-0.7 -7.3 -9.9 -4.0 -15.2 -11.4 -0.0 -0.4 -1.0 -0.1 -1.0 -3.2 0.0 0.3 0.2 0.0 0.1 0.1 0.0 -1.7 -2.2 -2.5 -17.5 -15.9 0.0 -2.0 -2.6 -0.6 -8.9 -7.9
aE Intra is the sum of the five intramolecular energies or of those of the two monomers when dealing with dimers. b For the Sc and St conformations of NMECC under a pressure of 15 kbar (see Chart 5 for the molecule numbering).
Table 5. Conformational Energy Minima (kcal/mol) and Conformational Molecular Volume (mL) of NMC Based Chiral Crotonate under Relevant Conformationsa 1 atm (de ) 7%)
15 kbar (de ) 18%)
conf.
E
∆E
V
∆V
E
∆E
V
∆V
Sc St Tc Tt Ac At
503.5 502.1 504.8 504.4 504.0 503.7
1.4 0.0 2.7 2.3 1.9 1.6
180.8 178.8 184.4 184.2 180.2 180.1
2.0 0.0 5.6 5.4 1.4 1.3
364.3 363.8 365.8 365.5 364.5 364.2
0.5 0.0 2.0 1.7 0.7 0.4
151.0 153.9 159.5 159.4 157.1 157.8
0.0 2.9 8.5 8.4 6.1 6.8
a
(See Table 2.)
energy between methanols 3 and 4 is larger than that between NMECC 1a and methanol 3; this feature explains the preference for the methanol-methanol hydrogen bond over the methanol NMECC 1a one in the case of methanol 4. trans-2-(2-Naphthylmethyl)cyclohexyl-(E)-crotonate (NMCC) 1b. For this compound, the values of Table 5 show that at 1 atm the St conformation is the most stable by 1.4 kcal/mol and that the axial arrangements have the same energy as the Sc form. At this point, we have to notice that the energy difference between the St conformer on one hand and the Sc and axial arrangements on the other does not appear in complete agreement with NMR3 results concerning this compounds. These data, obtained in solution, tend to show that the NMCC stacked conformations are significantly less populated than in the case of NMECC 1a, while the present computations give for both compounds similar populations of the stacked forms. However, our results give the correct overall trend when going from 1a to 1b since the average energy difference between the stacked and nonstacked conformers is 2.7 kcal/mol for the former and 1.7 kcal/mol for the latter. Moreover, our calculations do
Figure 3. Most stable geometrical arrangements of stacked and axial conformations of the transoid form of NMCC under 15 kbar.
Chart 6
not take into account the solvent effect on the stability of the different conformers in solution, a contribution that can destabilize the stacked conformations with respect to axial and trans arrangements. In addition, the ThorpeIngold23 effect, which produces variations of the valence angles at the tethering carbon between the cyclohexyl ring and the naphthyl group, is not taken into account in the present computations since with the version of SIBFA used in this study, only the torsion angles are varied. Since it has been shown24 that these variations can be of some importance, we carried out some computations, at 1 atm, to get an estimation of the (de)stabilizing contribution of such variations on the energy difference between the St and At conformers. In a first step, we carried out a complete optimization, using AM1, of 2,3-dimethyl-2-phenylbutane (4) to model compound 1a and of 1-phenyl-2-methylpropane (5), to model compound 1b (Chart 6). These computations show that replacing the methyl groups of 4 by hydrogens produces an increase of 3° of the angle between the plane of the phenyl and isopropyl C-C bond. The calculations carried out for the transoid stacked and axial conformations of NMCC with an angle of 112.5° instead of 109.5° for the C-C-C valence angle concerned by the Thorpe-Ingold effect, the energy difference between the St and At conformations is of 1.1 kcal/mol instead of 1.6 kcal/mol as reported in Table 5. This result shows that, for NMCC, the difference between the theoretical and experimental results concerning the most stable conformation in solution, as deduced from NMR, is not due to inacuracy of SIBFA method/parametrization. Under high pressure, the two axial and the two stacked conformations have energies which differ by 0.4-0.7 kcal/ mol with again the St as the most stable. The stabilization by pressure of the axial conformers can be related to the smaller difference of volume between the axial and stacked conformers than in the case of the preceding compound. Figure 3 shows that with the At conformation the two pro-R and pro-S approaches of the crotonate (23) (a) Beesley, R. M.; Ingold, C. K.; Thorpe, J. F. J. Chem Soc. 1915, 107, 1080. (b) Ingold, C. K. J. Chem. Soc. 1921, 119, 305. (24) von R. Schleyer, P. J. Am. Chem. Soc. 1961, 83, 1368.
Pressure on the Selectivity of a Michael Addition
J. Org. Chem., Vol. 64, No. 13, 1999 4731 Table 7. Contributions (kcal/mol) to the “Pretransition State” Complex Energy for the Sc, St, and At Conformations of NMCCa,b intermolecular EDisp ERep
Figure 4. Most stable geometrical arrangements of pro-R (A) and pro-S (B) complexes formed by axial transoid NMCC, three methanols and diphenylmethaneamine under 15 kbar. Table 6. Conformational Energy Minima (kcal/mol) and Conformational Molecular volume (mL) of the Reactive Complexa 1 atm (de ) 7%)
15 kbar (de ) 18%)
conf.
E
∆E
V
∆V
E
∆E
V
∆V
Sc pro-S St pro-R Tc pro-S Tt pro-R Ac pro-S At pro-R At pro-S
1133.8 1133.1 1137.4 1136.1 1137.0 1133.8 1135.0
0.7 0.0 4.3 3.0 3.9 0.7 1.9
372.5 368.2 374.9 377.5 369.7 370.9 372.5
4.3 0.0 2.4 9.3 1.5 2.7 4.3
800.6 796.7 804.4 800.2 802.1 797.6 798.1
3.9 0.0 7.7 3.5 5.4 0.9 1.4
309.3 309.9 316.3 316.3 312.3 314.3 314.3
0.0 0.6 7.0 7.0 3.0 5.0 5.0
a Including the NMC based chiral crotonate under relevant conformations (see Table 2), one diphenylmethaneamine and three methanol molecules.
moiety by the amine are possible. The values of Table 6, which concern the “pretransition state” complex, show that for this crotonate at 1 atm the At conformation corresponds to arrangements that have a stability close to that in which it adopts the St form and equal to that of the Sc one. However, one should notice that the most stable axial conformation is pro-R as is the St arrangement. The values corresponding to the complex under 15 kbar show that in this case the pressure tends to stabilize the axial-pro-S conformation that now has an energy differing only 0.5 kcal/mol from that of the axial-pro-R arrangement. The stabilization of the axial-pro-S form is in qualitative agreement with the low de increase induced by pressure in this case.3 The values reported in Table 7 show that for the complex with the stacked conformation of the chiral crotonate the difference of stability between the cisoid and transoid conformation has the same origin as in the case of NMECC 1a since the variation of the value of the different contributions follows the same trends. The values concerning the complex with an axial crotonate conformer show that the intramolecular energy is the factor responsible for the stabilization of the axial-pro-S complex. Conclusion Using a relation between compressibility and atomic van der Waals radii, based on previous work by Bridgman16 and Isaacs,15 SIBFA systematics has been modified to take into account the influence of pressure on interand intramolecular energies. With such a procedure, the parametrization of the method is not modified. Therefore, it can reasonably be expected that the quality of the results8 is the same at all pressures. The results obtained
Sc St At pro-R 1-2 1-3 1-4 1-5 2-3 2-4 2-5 3-4 3-5 4-5 At pro-S 1-2 1-3 1-4 1-5 2-3 2-4 2-5 3-4 3-5 4-5
-18.8 -21.0 -20.8 -3.5 -5.2 -1.7 -2.5 0.0 0.0 -1.4 -4.4 -1.3 -0.8 -20.0 -3.7 -6.9 -1.4 -2.2 0.0 0.0 0.0 -4.7 -0.1 -1.6
41.3 49.1 49.9 3.1 13.8 0.6 1.5 0.0 0.0 2.5 11.8 0.7 0.5 50.2 3.2 16.0 0.5 1.0 0.0 0.0 0.0 13.4 0.0 3.1
intramolecular
EPol
ECT
EMTP
EInter EaIntra
-10.0 -12.0 -12.2 -1.3 -3.2 -0.5 -0.1 0.0 0.0 -0.7 -3.5 -0.2 -0.3 -12.5 -1.3 -3.8 -0.5 -0.2 0.0 0.0 0.0 -3.7 0.0 -0.8
-5.6 -6.7 -6.7 -0.8 -3.1 -0.1 -0.1 0.0 0.0 -0.3 -2.5 0.0 0.0 -7.2 -0.8 -3.5 0.0 0.0 0.0 0.0 0.0 -2.5 0.0 -0.4
-42.9 -49.7 -50.2 -7.6 -13.0 -3.9 -0.9 0.0 0.0 -5.0 -17.4 -1.1 -1.4 -49.7 -8.2 -15.1 -3.7 -0.7 0.2 0.0 0.1 -17.6 -0.2 -5.1
-35.9 -40.1 -40.1 -10.2 -10.7 -5.5 -2.0 0.0 0.0 -5.0 -15.8 -1.9 -1.9 -38.9 -10.8 -13.2 -5.1 -2.2 0.2 0.0 0.1 -14.9 -0.3 -4.7
836.5 836.8 837.7 475.1 475.1 475.0 457.2 231.1 231.2 213.2 231.1 213.2 213.5 837.0 474.8 474.7 474.7 456.4 230.4 230.4 212.2 230.5 212.1 212.1
Etot 800.6 796.7 797.6 464.9 464.4 469.5 455.1 231.1 231.2 208.2 215.3 211.3 211.3 798.1 463.9 461.5 469.6 454.2 230.6 230.4 212.3 215.6 211.8 207.4
a E Intra is the sum of the five intramolecular energies or of the sum of those of the two monomers when dealing with dimers. b Under a pressure of 15 kbar and for the different possible dimers in the pro-R and pro-S at forms (see Chart 5 for the molecule numbering).
for two stereogenic crotonates show that the energy variation with pressure depends significantly upon molecular conformation. Thus, for one of the compounds studied, NMCC 1b, the most stable conformation is different at 1 atm and 15 kbar. Our results show also that pressure modifies intermolecular interactions since for NMECC 1a the complex with diphenylmethaneamine plus three methanols is more stabilized when the crotonate has a St conformation than when it has any of the other possible arrangements. In this case, it is a feature that broadens the range of energy of the different possible complexes. As a consequence, pressure makes the addition more diastereoselective in agreement with experiment.3 In the second case, NMCC 1b, this variation does not occur because several conformations undergo a stabilization of comparable magnitude. This result prevents the increased stereoselectivity obtained for NMECC 1a, a feature in agreement with experiment. At a pressure of 15 kbar, the molecular and complexes volumes are calculated to be much smaller than at 1 atm but the variation of the molecular volumes with molecular conformation is not noticeably modified by pressure. On the other hand, our results show clearly that the Tc and Tt conformations with the largest volumes at 1 atm are those which have their energy less modified by pressure. This feature is also found but to a smaller extend in the case of the complexes. This finding is in line with those of Diedrich and Kla¨rner25 concerning the role of molecular volume and packing coefficients in the effect of pressure on pericyclic reactions. Our results tend to show that a computational method taking into account the compressibility of the system studied through the (25) Diedrich, M. K.; Kla¨rner, F.-K. J. Am. Chem. Soc. 1998, 120, 6212.
4732 J. Org. Chem., Vol. 64, No. 13, 1999
volume variation appears to be able to give, at least at the qualitative level, structural modifications due to pressure. Since SIBFA has provided pertinent results for peptides and proteins,8,26 it should, with the modifications that we have introduced, be able to account for the conformational changes taking place in pressure-denatured proteins.27 (26) (a) Gresh, N.; Roques B.-P. Biopolymers 1997, 41, 145. (b) Garmer, D. R.; Gresh, N.; Roques, B.-P. Proteins, Struct., Funct, Genet. 1998, 31, 41. (27) Hummer, G.; Garde, S.; Garcia, A. E.; Paulaitis, M. E.; Pratt, L. R. Proc. Natl. Acad. Sci. 1998, 95, 1552.
Dumas et al.
Acknowledgment. The computations have been carried out at the CRIHAN (Centre de Ressources Informatiques de Haute Normandie), Mont Saint Aignan, within the framework of the “Re´seau Normand pour la Mode´lisation Mole´culaire” and the C.C.R. (Centre de Calcul pour la Recherche, Universite´ P. & M. Curie), Paris. We thank Dr. Jean d’Angelo and Dr. Nohad Gresh for fruitfull discussions.
JO9825308
