Md Abdul Shafeeuulla
Khan
a,
Ji
Zhang
b,
Koushik Das
Sarma
*c and
Bishwajit
Ganguly
*a
aAnalytical Discipline and Centralized Instrument Facility, Central Salt & Marine Chemicals Research Institute (Council of Scientific and Industrial Research), Bhavnagar, Gujarat 364 002, India. E-mail: ganguly@csmcri.org; Fax: (+91)-120-2586776
bKing Kai-Ray Pharma, Mail Box 3016, Princeton, NJ-08540, USA
cJubilant Chemsys Ltd., B-34, Sector-58, Noida, UP 201301, India. E-mail: kdsarma@gmail.com; Fax: (+91)-278-2567562
First published on 17th July 2012
Oxygenated five-membered ring compounds are more widespread in nature than the corresponding aza-analogues due to their more interesting biological activity. In our earlier studies, we have reported a simple methodology to synthesize novel synthons like 5-(1′-hydroxy)-γ-butenolide and 5-(1′-hydroxy)-γ-butyrolactam through the direct vinylogous aldol addition reaction. The observed diastereoselectivities for γ-butyrolactones and γ-butyrolactams were not understood in that report. We have reported herein a detailed computational study to understand the observed diastereoselectivity of α,β-dichloro-γ-butenolides and α,β-dichloro-γ-butyrolactams. The transition state calculations performed with B3LYP/cc-pVDZ level of theory to examine the diastereoselectivity of α,β-dichloro-γ-butenolides and α,β-dichloro-γ-butyrolactams have been found to be in good agreement with the observed experimental results. The origin of selectivity is examined by an activation strain model and molecular electrostatic potential (MESP) surfaces. The distortion energies calculated using the activation strain model revealed the origin of diastereoselectivities for butenolides and butyrolactams. The reversal in stereochemical outcome of α,β-dichloro-γ-butenolides and α,β-dichloro-γ-butyrolactams is due to the presence of a large substituent at the nitrogen center of the γ-butyrolactams. We have further revealed the enhanced diastereoselectivity of these butenolides while reacting with the aromatic aldehydes bearing ortho substituents. The generated molecular electrostatic potential (MESP) surfaces have shown an additional Vmin of value = −115.0 kcal mol−1 responsible for enhanced diastereoselectivity with ortho-substituted benzaldehydes as compared to unsubstituted benzaldehyde. The difference in transition state distortion energies of syn and anti aldol products are found to correlate with their corresponding differences in activation energies.
One of our recent efforts explored a simple and facile method to synthesize novel synthons like 5-(1′-hydroxy)-γ-butenolide and 5-(1′-hydroxy)-γ-butyrolactam from a direct vinylogous aldol addition reaction between α,β-dichloro-γ-butenolide/α,β-dichloro-γ-butyrolactam and aromatic aldehydes.19 Diastereoselectivity has been observed and the effects of substitution on the benzaldehyde on diastereomeric ratio have also been ascertained. From 1H NMR spectra analysis, it has been observed that the syn/anti ratio of aldol product with unsubstituted benzaldehyde is 2/1 (Table 1). It has also been observed that the syn/anti ratio of aldol product with para- and meta-substituted benzaldehydes is same as observed in the case of unsubstituted benzaldehyde, whereas, the syn/anti ratio was enhanced to 5/1 with ortho-substituted benzaldehyde (Table 1).19 We have speculated the role of steric effect between the aromatic aldehydes and chlorine atoms of α,β-unsaturated γ-butenolide or α,β-unsaturated γ-butyrolactam in their transition state geometries to rationalize the observed diastereoselectivity (Scheme 1).19 Interestingly, α,β-unsaturated γ-butyrolactams have shown opposite diastereoselectivity with respect to α,β-unsaturated γ-butenolides (Table 1) under the same reaction conditions. These observations were intriguing and warrant further investigations. In the present article, we have undertaken a computational study to understand the experimentally observed diastereoselectivity for the formation of 5-(1′-hydroxy)-γ-butenolide and 5-(1′-hydroxy)-γ-butyrolactam derivatives by employing DFT calculations. The energetic data of transition state geometries predicted with B3LYP/cc-pVDZ showed a good agreement with the syn/anti diastereomeric ratios determined experimentally. The origin of diastereoselectivity for 5-(1′-hydroxy)-γ-butenolide/5-(1′-hydroxy)-γ-butyrolactam derivatives has been rationalized on the basis of a distortion/interaction energy model20 and steric effects. Molecular Electrostatic Potential (MESP) analysis performed on the transition state geometries revealed the enhanced selectivity of α,β-unsaturated γ-butenolide with ortho-substituted benzaldehyde.
Scheme 1 The Newman projections of possible transition state conformers for syn and anti products.19 |
Entry | Ar-CHO | Product | syn/anti |
---|---|---|---|
1 | 2/1 | ||
2 | 2/1 | ||
3 | 2/1 | ||
4 | 5/1 | ||
5 | 1/5 | ||
6 | 1/5 |
Scheme 2 α,β-Dichloro-γ-butenolide, benzaldehyde and its derivatives 1–4. |
For simplicity, we have adopted the following notations; 1S is the transition state formed between 1 and X; 2S for 2 and X; 3S for 3 and X; 4S for 4 and X; here S stands for syn orientation. Similarly, 1A, 2A, 3A and 4A are considered for the corresponding anti transition state geometries. All optimized transition state geometries (1S–4S and 1A–4A) are shown in Fig. 1 along with the corresponding solvated Gibbs free energy differences. The carbanion of X moves toward the carbonyl carbon center of 1 at a distance of 1.911 Å in syn orientation to form a transition state (1S). The single imaginary frequency of 1S is 237.1i cm−1, with the normal vibrational mode characterizing the –C⋯C– bond formation. The corresponding bond distance and single imaginary frequency for the anti oriented transition state (1A) are 1.900 Å and 228.3i cm−1, respectively. The B3LPY/cc-pVDZ calculated free energies show that the syn transition state (1S) is energetically more stable by 2.9 kJ mol−1 compared to the corresponding anti transition state (1A) (Fig. 1). The cyano substitution at the para position of benzaldehyde (2) with α,β-dichloro-γ-butenolide (X) also predicted the syn transition state 2S is energetically preferred by 2.1 kJ mol−1 than the anti transition state 2A. The transition states calculated with the meta methoxybenzaldehyde 3 and X showed the preference for the formation of syn transition state (Fig. 1). The enhancement in the formation of syn transition state was observed with 4 and X, where the nitro group is substituted at the ortho position. The syn transition state 4S was found to be 13.0 kJ mol−1 energetically more stable than the corresponding anti transition state 4A. The influence of ortho substitution was further exemplified with ortho-methoxybenzaldehyde. The calculated results suggest that the formation of syn product is preferred by 10.0 kJ mol−1 compared to the formation of anti product in terms of Gibb's free energy as observed in the case of ortho-nitrobenzaldehyde (ESI, Fig. S1†).
Fig. 1 B3LYP/cc-pVDZ optimized transition state geometries of syn (left column) and anti (right column) aldol products of γ-butenolides in methanol along with the corresponding Gibb's free energies () and electronic energies with ZPVE correction [] (grey = carbon, red = oxygen, blue = nitrogen, white = hydrogen and green = chlorine). |
The role of the methoxy(oxo)methanide group at the para position in 4 was also investigated towards the diastereoselectivity of X. We replaced the methoxy(oxo)methanide group by a hydrogen atom and determined the syn and anti transition states (ESI, Fig. S2†). The calculated free energy difference showed the larger preference for the syn-aldol product over anti-aldol product (ESI, Fig. S2†). Such preference for the formation of syn-aldol over anti-aldol is in good agreement with the observed diastereoselectivity for the aldol reaction between substituted benzaldehydes and α,β-dichloro-γ-butenolide (Table 1). Further, the enhanced syn/anti product ratio for 4 with X compared to 1–3 was also predicted by the transition state calculations (Fig. 1).
Turning to the diastereoselectivity exhibited by α,β-unsaturated γ-butyrolactams, we have considered 5 and 6 to partake in a vinylogous aldol addition reaction with 4-NO2-benzaldehyde (Y) as suggested in our previous19 experimental study (Scheme 3). The optimized syn and anti transition state geometries of (5S, 5A, 6S and 6A) are shown in Fig. 2 along with corresponding free energy differences. Here, 5S is the transition state formed between 5 and Y; 6S is the transition state formed between 6 and Y; again S stands for syn orientation. Similarly, 5A and 6A are considered for the transition state geometries in anti orientation. The carbanion of 5 moves toward the carbonyl group carbon center of Y at a distance of 1.877 Å in syn orientation to form a transition state 5S. The single imaginary frequency of 5S is 187.3i cm−1, with the normal vibrational mode characterizing the –C⋯C– bond formation. The corresponding bond distance and single imaginary frequency for the anti oriented transition state (1A) are 1.887 Å and 192.7i cm−1, respectively. Similarly, the bond distances of 6S and 6A are 1.875 Å and 1.892 Å and the imaginary frequencies are 186i cm−1 and 196.3i cm−1, respectively. The free energy differences calculated with B3LYP/cc-pVDZ level for the syn and anti transition states for 5 and 6 with Y showed that the anti forms are more stable by 5.5 kJ mol−1 and 7.6 kJ mol−1, respectively, than their corresponding syn forms (Fig. 2). This opposite trend in diastereoselectivity of α,β-unsaturated γ-butyrolactams compared to diastereoselectivity of α,β-unsaturated γ-butenolides is in well agreement with the experimentally observed diastereoselectivity of these compounds (Table 1).
Fig. 2 B3LYP/cc-pVDZ optimized transition state geometries of syn (left column) and anti (right column) aldol products of γ-butyrolactams in methanol along with the corresponding Gibb's free energies () and electronic energies with ZPVE correction [] (grey = carbon, red = oxygen, blue = nitrogen, white = hydrogen and green = chlorine). |
Scheme 3 α,β-unsaturated γ-Butyrolactams (deprotonated forms) 5, 6, 7 and the para-nitrobenzaldehyde Y chosen for the study. |
The reversal of the diastereoselectivity in the case of α,β-unsaturated γ-butyrolactams (5 and 6) may be due to the presence of large group attached to the nitrogen atom of the five membered ring which is absent in the α,β-unsaturated γ-butenolides. We have replaced the large group on the nitrogen atom in both syn and anti transition state geometries with a hydrogen atom. These model transition state geometries 7S and 7A were optimized at the same level of theory (Fig. 2). Interestingly, the computed free energies show that the syn transition state (7S) is energetically stable by 0.5 kJ mol−1 compared to the anti transition state (7A) geometry of α,β-unsaturated γ-butyrolactams as we have seen in the case of α,β-unsaturated γ-butenolides. Therefore, it is now clear that the stabilization of the anti transition state over the syn transition state is due to the presence of large groups on the nitrogen atoms of the five membered rings of α,β-unsaturated γ-butyrolactams, which causes a larger steric repulsive interaction for the syn orientation over the corresponding anti orientation (Fig. 2). To further examine the effect of ortho-substituted benzaldehyde on the diastereoselectivity of α,β-unsaturated γ-butyrolactams, we have also computed the transition states formed during the reaction of ortho-nitrobenzaldehyde with N-substituted α,β-unsaturated γ-butyrolactam (5) and N-unsubstituted α,β-unsaturated γ-butyrolactam (7). The difference in activation Gibb's free energy data shows that the syn transition state formed between 7 and ortho-nitrobenzaldehyde is 13.0 kJ mol−1 more stable than the corresponding anti transition state as observed in the case of α,β-unsaturated γ-butenolides (ESI, Fig. S3†). However, N-substituted α,β-unsaturated γ-butyrolactam 5 also shows that the syn transition state is 9.6 kJ mol−1 more stable than the anti form in contrast to the reaction with Y (para-substituted benzaldehyde) (ESI, Fig. S3†). Due to the steric effect of ortho-substitution of benzaldehyde, the syn preference is enhanced with the N-unsubstituted α,β-unsaturated γ-butyrolactam 7, whereas, the selectivity is reversed with the N-substituted α,β-unsaturated γ-butyrolactam 5 as compared to the case of para-substituted benzaldehyde. That means, the steric effect of the ortho-substituent outweighs the steric effect of the larger substituent at the N-atom of α,β-unsaturated γ-butyrolactam. Therefore, it appears that the positional isomers of benzaldehydes can drive the diastereoselectivity of α,β-unsaturated γ-butyrolactams.
The origin of diastereoselectivity for the aldol reactions of α,β-unsaturated γ-butenolides/α,β-unsaturated γ-butyrolactams with benzaldehyde derivatives was examined with distortion energy model calculations. The distortion/interaction energy or activation strain model has been applied for many reactions in literature viz. aldol addition, Diels–Alder reactions of cyclo additions and reactions involving chiral dienes and achiral dienophiles, Evans chiral auxiliary-based asymmetric aldol reactions and 1, 3-dipolar cycloadditions by considering transition state distortion energy or interaction energy as a controlling factor towards reactivity.20–25 For syn and anti transition states obtained between the α,β-unsaturated γ-butenolides or α,β-unsaturated γ-butyrolactams and benzaldehyde derivatives, activation strain model was employed to segregate the distortion and interaction energies for the respective systems. The activation strain model calculations suggest that the distortion energies for the syn transition states of the aldol reaction between α,β-unsaturated γ-butenolides and benzaldehyde derivatives are lower than the corresponding anti transition states (ESI, Table S1†). Similarly, the distortion energies calculated for the α,β-unsaturated γ-butyrolactams and benzaldehyde derivatives showed that anti-transition states are energetically less distorted than the syn ones (ESI, Table S1†). A linear correlation is found between the difference of activation energies and the corresponding distortion energies for all the studied reactions with R2 = 0.941 (Fig. 3). Analysis of the transition states employing the activation strain model indicates that the diastereoselectivity of these aldol reactions is mainly governed by the geometric distortions of reactant molecules to achieve corresponding transition state geometries.
Fig. 3 Plot of linear fit computed at B3LYP/cc-pVDZ for the difference of distortion energies between syn and anti transition states (1–7) versus the difference of activation energies between syn and anti transition states (1–7) with energies given in kJ mol−1. |
The transition state calculations and activation strain model shed light on the origin of diastereoselectivity for α,β-unsaturated γ-butenolides and α,β-unsaturated γ-butyrolactams with benzaldehyde derivatives. However, the reason for the enhanced syn/anti selectivity in the case of 4 with X is yet to be answered. To address this point, we have generated the MESP surfaces of syn and anti transition states of 1 and 4 with X. MESP is an important and widely used topographical quantity for understanding molecular reactivity, intermolecular interactions, molecular recognition, electrophilic reactions, substituent effects etc.26–37 MESP studies have also been used as tool for measuring electrostatic interactions (repulsions or attractions) due to the π-electron cloud of aromatic rings.38–40 MESP analysis gives the most negative valued point (Vmin) in electron rich regions obtained through topography calculations in any system. The MESP calculations show that in the case of 4A, there is an additional Vmin near to ortho-nitro group and chloro-substituents of the α,β-dichloro-γ-butenolide moiety (Fig. 4). Such large Vmin values can cause severe repulsion to the neighboring partner and hence would destabilize the anti aldol product more compared to 1A. The calculated free energy differences corroborate the MESP results.
Fig. 4 The B3LYP/cc-pVDZ generated MESP surfaces in methanol for 1S and 1A at an isosurface value of −410.0 kJ mol−1, and 4S and 4A at an isosurface value of −340.0 kJ mol−1. Red spots represent location of Vmin at the π-electron cloud of benzene and the black dot represents location of Vmin at the oxygen of the nitro group at the ortho-position. |
(1) |
The activation strain or distortion energy model is used to analyze the energy barriers. This method is a fragment-based approach used to understand chemical reactions and their associated barriers.20 The distortion energy depends on the ability of the reactants to reorganize into their deformed structures in the transition state, which is devoid of any interaction between the reactants.21 The distortion energy (ΔEd‡) or interaction energy (ΔEi‡) model is a powerful tool to explain the reactivity in bimolecular reactions.22 By definition, activation energy is ΔE‡ = ΔEd‡ + ΔEi‡. The starting point is the two separate reactants, which approach from infinity and begin to interact and deform each other. The distortion energy computed by this method is associated with the strain caused by deforming the individual reactants. Throughout this article, energies and distances are given in kJ mol−1 and angstroms, respectively.
Footnote |
† Electronic supplementary information (ESI) available: B3LYP/cc-pVDZ optimized Cartesian coordinates of all stationary points, including all the thermodynamic data generated from frequency calculation and imaginary frequencies of all the transition state geometries. See DOI: 10.1039/c2ra21203j |
This journal is © The Royal Society of Chemistry 2012 |