Journal of the Chilean Chemical Society
versión On-line ISSN 0717-9707
J. Chil. Chem. Soc. v.54 n.1 Concepción 2009
J. Chil. Chem. Soc, 54, N° 1 (2009); págs: 99-104
THE EFFECT OF C2 AND C5 SUBSTITUTIONS ON THE NITROGEN LONE PAIR INVERSION IN 1,3-OXAZINES"
J. GUILLERMO CONTRERAS*, S. MARCELA HURTADO
Facultad de Ciencias Químicas, Universidad de Concepción, casilla 160-C, Concepción, Chile. * e-mail: email@example.com
A series of N-methyl-2-R1,R2-5-R3,R4- 1,3-oxazines (where R are either H or p-nitro phenyl group and R3,4 = H or CH3- group) have been theoretically studied by ab initio methods in the frame of MO theory. Optimized geometries were obtained at HF and B3LYP/6-31G" levels, whereas energies were derived at MP2 and B3LYP levels employing the same basis set. Geometries obtained at DFT level are needed to obtain consistent values and good correlations between the experimental and theoretical 1H-NMR data. The nitrogen lone pair would be axial orientated in I and II, since the N-methyl group lies in the equatorial position in the gas phase, whereas in compounds III and IV the methyl axial conformers increase their concentration, implying that inversion of the nitrogen lone pair takes place in an appreciable extent. PCM method predicts that solvents (CCl4, CHCl3, CH2Cl2) produce a significant effect on conformers preference, but just a modest influence on the equilibrium geometries and 1H-NMR properties.
Keywords: 1,3-Oxazines, Nitrogen Inversion, Lone-Pair Inversion Barrier, Substituents Effects
The study of the spatial requirements of nitrogen nonbonding pairs of electrons imply the determination of the preferred orientation of alkyl group attached to nitrogen in the heterocycle being studied. Using dipole moment measurements, Urbanski1 assigned a chair form to N-methyl or ethyl-5-nitro-5-methyl-1,3-oxazines. The analysis ofthe 1H-NMR multiplet pattern and the position ofthe methylene signals in nitrocyclohexano suggest that the only species present at room temperature possesses a chair conformation with the nitro group being in an equatorial position2. Studying the kinetics ofthe cycle inversion process in some tetrahydro-1,3-oxazines by variable temperature NMR spectra, Lehn et al3 propose that nitrogen inversion process is faster than the ring inversion so that the 1H-NMR spectra will show an AB system for C2 protons, whereas the N-CH3 group will show up as single averaged signal. The experimental data support a conformational equilibrium between two equivalent chairs3,4. Allingham et al5 extended the Urbinski's series1,6 by including the n-propyl, i-propyl and cyclohexyl derivatives. The 1H-NMR data allow to determine the preferential conformation by studying the Jgem for C2, C4 and C6 protons. These coupling constants would be sensitive to orientation of the nitrogen lone pair. Katritzky et al7 found that ΔG°138 is ca. -0.16 kcal/mol for the equatorial ↵ axial equilibrium in 3-methyl-l-oxa-3-azacyclohexane. The mayor component is that having an axial methyl group. The calculated nitrogen inversion barrier is ca. 7.6 kcal/mol, whereas in pyperidine, a value of 8.7 kcal/mol has been reported8.
In the present work a series of N-methyl-2,5-substituted-1,3-oxazines have been studied with the aim to determine the effect of substitution on C5 and/or C2 on electronic delocalization, steric factors that might be important in the nitrogen inversion process. On the other hand, it looks promising to study the influence of these substitutions on the NMR properties. In fact, we have paid attention to relations between chemical shifts, 1J and Jgem spin-spin coupling constants with energy barriers and hyperconjugative interactions. The inversion barriers for some pairs of species have also been determined and compared with the experimental values when available. Important part of this work is devoted to the search the TSs connecting both conformers in the equatorial ↵ axial reaction. For this purpose we have assumed that, a nitrogen atom with three substituents and a nonbonding lone pair of electrons has a pyramidal geometry capable of inverting its configuration. The ground state has a nominally a sp3-hybridized nitrogen with the lone pair occupying an orbital that is approximately a sp3 hybrid, and the transition states has an sp2-hybridized nitrogen with the lone pair in a p orbital9. In the present work the preferred orientation of the methyl group on nitrogen has been studied and compared with the scarce experimental data. The C5 and C2 substitution by electron releasing or withdrawing groups allow inferring the effect on the nitrogen inversion barriers. The calculated protón chemical shifts and coupling constants have also been used to derive the relative % concentration ofthe axial and equatorial species.
Geometry optimizations were carried out at HF and B3LYP levels using the 6-31G** basis set. In all the cases, GAUSSIAN 03 series of programs were used10. The initial geometries were those previously reported for other oxazines11. Harmonic vibrational frequencies and IR intensities were calculated at the equilibrium geometries. In all cases, the equilibrium geometries yielded all real frequencies and thereby the 8 derived structures at either HF or B3LYP/6-31G** levels are local minima. Energy calculations were carried out using various basis sets in order to obtain the optimum level to be used in the calculation of the NMR properties. GIAO (Gauge-Independent Atomic Orbital) method12,13 was applied to derive the magnetic tensors. Chemical shits (in ppm) were referred to tetramethylsilane (TMS) standard (isotropic magnetic tensor = 32.1 ppm for 1H)14. Natural Bond Orbital15,16 analysis was performed to study the most prominent interactions in each molecular structure. NBO analysis allows the separation ofthe molecular energy into two fundamental contributions: the total energy (where delocalization is included), and the energy of the hypothetical Lewis structure. Accordingly, the electron are strictly located in bonds and lone pairs. The interactions between filled and antibonding (or Rydberg) orbitals represent the deviation ofthe molecule from the ideal Lewis structure and can be used as a measure of delocalizations (Edeloc = Etotal - ELewis). The energy of the hypothetical Lewis structure with strictly localized bonds is obtained by removing all off-diagonal elements from the Fock matrix and computing one SCF cycle. The hyperconjugation interaction energies were obtained from the second-order perturbation approach: , where the integral corresponds to the Fij. matrix element between interacting orbitals i and j whose energies are εσ and . 2 stands for the population ofthe donor orbital. There is good linear correlation between the second-order perturbation hyperconjugative energies and the deletion procedure17. All the NBO calculations were carried out at B3LYP/6-31G**. To obtain the theoretical 1H and 13C chemical shifts (δ), the magnetic shielding tensors were calculated at B3LYP/ 6-311+G**7/B3LYP/6-31G** level of theory14 using the Gauge- Independent Atomic Orbital (GIAO) method15. Chemical shifts were derived regarding to the NMR isotropic magnetic shielding tensors (in ppm) from the corresponding standard tensors (TMS, 31.904 ppm for :H and 183.436 ppm for 13C).
Since TS are local máxima in the reaction pathways, they are no accessible to experimental studies and thereby no direct information on their structures are available. Theoretical methods are the only way to infer some knowledge about transition status at present. Hammond has suggested20 that the energetic barrier is closely related with the reactants and products structures. Thus, Hammond's postúlate is a useful tool to connect the energy barrier position in the potential energy surface (PES) with the exo or endothermic character of the reaction. Accordingly, if the TS is close in energy to a stable species, its structure will be alike. Brónsted coefficient (eq. 1) and Marcus equation (eq.2) allow to locate a TS.
where: ΔE≠ is the activation energy, and ΔE0 corresponds to the reaction energy. According to these equations an exothermic reaction will possess a TS similar to the reactant and ß will be < ½ .. When ß is > ½ the TS resembles the products , and the reaction will be an endothermic one..
To model the solvent-solute interaction the polarizable continuum model (PCM)21, was employed. PCM considers the solvent as a continued unstructured medium whose polarization can be represented by a macroscopic constant such as the dielectric constant,e. The cavity hosting the solute is constructed with a number of overlapping spheres centered at the solute nuclei and possessing radii proportional to the Van Der Waals ones. The solute induces a polarization of the medium that brings on a charge density in the cavity surface in turn. Our large experience in our laboratory using this method to derived solvation free energies, allows us to be confident in the solvent modelation.
RESULTS AND DISCUSSION
1,3-oxazine is a saturated six-member heterocycle containing both nitrogen and oxygen with a distorted chair-like structure. Oxazines derivatives studied here are distinguished by their substituents as well as the orientation of the -CH3 group. The designation axial or equatorial refers to the orientation of the methyl group attached to N atom, whereas I, II ,III and IV refers to the type of substituents on C2 and C5 atoms (see figure 1).
From table 1, it can be inferred that difference between the dihedral angles that define each side of the ring, C201C6C5 and C2N3C4C5 (At), show that I ax and II ax conformers are more distorted at the O than at the N-side. In fact, the presence of two methyl groups on C5 leads Δτ values of ca. 11-14°, whereas in the corresponding equatorial conformers Δτ are ca 3-4° only. The presence of a large substituent on C2 seems to exert no much influence on the geometry of the oxazine ring since Δτ is just ca. 4° for both axial and equatorial IV derivatives. III ax is an exception as no substituents on C2 or C5 are present and behaves like I and II axial conformers, though to a lesser extend.
Important C-H bond distances are C2-H2a and C4-H4a for being the ones varying more within this series of conformers. These distances are longer in the equatorial conformers, implying a delocalization of the type (LP)N3 → σ* C2-H2a and (LP)N3 → σ* C4-H4a that contribute to the molecular stability. Larger C-H bond lengths dependence with C2 and C5 substitutions is observed in the equatorial conformers. Shorter C6Ha distance occurs in all cases but Illeq. The plot of ring relative (with respect conformer III) bond distances vs. the corresponding distances for all other conformers studied here, allow to remo ve any other effect, such as lone-pairs delocalization or other related to the substituent position but the one produced by the methyl groups on C5 (or C2 for the p-nitrophenyl group). Fig.2 shows that the increase of the C4C5 AND C5C6 bond lengths are due to the presence of two CH3- groups in conformers I and II. The p-nitrophenyl group in conformers II and IV lead to larger C2N3 distances. Certain degree of linearity in the C2N3 distances throughout this series of conformers: I < III < II < IV.
Since p-nitrophenyl is a strong electron-withdrawing group, weakens and enlarges the C2N3 distance, whereas methyl being an electron-releasing group counterbalances the above effect in conformer II. The 01C2 bond distance found for the series of conformers of these substituted oxazines are longer for II and IV most probably due to the p-nitrophenyl group presence. The substituent on N3 in an equatorial position facilitates the chain/ring tautomerization in 1,3-oxazines18. The chain/ring tautomerism have been extensively studied by means IR a UV spectroscopies19 and shown that the tautomeric equilibrium takes place due to the substituent located in the aromatic ring. Thus, electron-withdrawing groups increase the cyclic form contribution, whereas electron-releasing substituents would favor the chain forms. In the axial structures, the 01C2N3 AND C5C4N3 angles are larger than the corresponding angles in the equatorial conformers. These angles are related to the anomeric effect so essential for the axial conformations. Considering these angles in structure Illax and taking into account that 1,3-diaxial repulsions between the N-substituent and Ha on C5 that are very small, the methyl (ax) group effect on C5 can be determined, as well as the effect on angle 01C2N3 exerted by the p-nitrophenyl group. The variation of the above mentioned angles for al species studied here are shown in table 2.
The value shown by angle C5C4N3 reflects the contribution of the anomeric effect and the 1,3-diaxial repulsive interactions producing an increase of ca. 2 and 3 degrees in the equatorial an ↵ axial structures, respectively. Angle 01C2N3 shows a rather modest influence of the concomitant effects already mentioned. The increase in this angle is due to the anomeric effect whereas the p-nitrophenyl decreases such angles to ca 1.5 degrees.
The 1H-NMR parameters for all conformers were calculated applying the GIAO method and compared with the experimental values previously reported by Lenh et al.3 Good agreement of calculated and experimental values support well or model for the local minima found for all conformers. The results are given in table 3
The theoretical NMR chemical shifts (8) were calculated at the B3LYP/6-311+G*//B3LYP/6-31G** level, whereas the coupling constants (J) were derived at B3LYP/6-31G*//B3LYP/6-31G** in CH2Cl2 or CHCl3 solutions. Experimental values of δ and J were found to correlate well with the theoretical results obtained at the B3LYP/6-311+G**//B3LYP/6-31G** level of theory. Table 3 shows this agreement for the 8 and J values. 1H-NMR spectra4,7 in CF2Cl2 (e = 3.50) suggest that the conformational equilibrium Illa Illeq is partially displaced to the equatorial form. In fact, Lehn et al4 reported that at the equilibrium, a 45% of the axial conformer is present, whereas the equatorial is in a 55% concentration. For the IV species, a 57% of equatorial and a 43% of the axial conformer was found. Analysis of the equilibrium between the two forms of III by Katrizky7 indicates that ΔG°138 = -0.16 kcal/mol implying a 56.72% of the axial conformer, i.e. reversing the results of Lehn et al. The strong 1,3 interactions in axial conformers I and II, allow to expect, confidently, that equatorial form will be the predominant one at the equilibrium. Comparing the average value between CCl4 and CHCl3, for compound III a value of ca. 51.11% of the axial component, can be predicted. This is a value halfway of the reported by Lehn and Katrizky. Closer values to Lehn can be obtained by employing solvents with larger dielectric constants. axial conformer can achieve concentration of ca. 60% in compound IV. Table 4 shows relative concentration% of the axial component in various solvents.
From table 4, it can inferred that the axial and equatorial conformers are in equilibrium with a slightly larger concentration of the equatorial conformer in compound III, whereas in IV the axial conformer dominates. The equilibrium in compounds I and III are strongly equatorial in the gas phase and in solution of the solvents Usted in table 4. The equatorial components in I and III shows the larger dipole moments and thereby their stabihzation in solvents of higher dielectric constants can not be rule out and it is likely that the % of axial conformers in compounds II and IV could become larger as e grows since their dipole moments is also larger (see fig.3)
Reasoning out, the % of the axial conformer in equilibrium II and IV Would also increase as the medium polarity increases. However, the dipole moments of the axial component of II and IV systems are larger than the equatorial one (see fig.3) and thereby such effect cannot be observed, most likely due to problems of PCM to handle the non-electrostatic contributions , such as cavitation and dispersion. It is also possible to think that large non-stabilizing interactions in the axial structures cannot be overéame for other dipole/dipole interactions with the solvent. In system II, the 1,3-interactions between N-methyl group and the axial methyl on C5 would also play an important role.
The transition states show , after NBO treatment at B3LYP/6-3 IG** level, nitrogen atom that reverses its equatorial character to an axial one through planar state as demonstrated by the increase of its "p" character in the lone pair of electrons (see table 5).
The TSs found exploring the potential energy surface describes well the inversion-rotation processes showing imaginary frequencies at ca 180 cm-1. Two out of four TSs for the equatorial B axial interconversion were also characterized by following their reaction paths using IRC (intrinsic reaction coordínate) calculations. In both cases IRC (I and II) demonstrated that the corresponding TSs conduct to the reactants and products. The results of these calculations are shown in figure 4.
The instability of TS for the equilibrium equatorial B axial in conformers I and II .indicates shows a strong repulsive effect between C5-CH3(Δτ) y N3-CH3 subtituents. Similar to that observed at the axial minima I and II. In this sense a calculation of Brönstead coefficient yield a ß = 0.564 for the equilibrium in I and ß = 0.547 for II. In this way the structural similarities, imply that in both cases the corresponding TSs would show up the products alike. I and II cases follow well the postúlate of Hammond and their endothermic character. For equilibrium III and IV Brönstead coefficients do not allow to infer if the TSs are reactant or product alike. In fact, values of p of ca. 0.494 (III) and 0.490(IV) indícate that neither the reactants ñor products are good model for the TSs.
From the above results, some interesting conclusions can be outlined:
|1.-||The ring O and N sides dihedrals allow us to conclude that axial I and II conformers are more distorted at the O than at the N-side.|
The presence of a large substituent on C2 exerts no much influence on the geometry of the oxazine ring. In fact, Δτ is just ca. 4o for both axial and equatorial IV derivatives.
IIIax is an exception as no substituents on C2 or C5 are present and behaves like I and Ilaxial conformers, though to a lesser extend
|3.-||Bulky p-nitrophenyl is a strong electron-withdrawing group and weakens and enlarges C2N3 distance, whereas methyl being an electron-releasing group counterbalances the above effect in conformer II. The 01C2 bond distances found in the series are longer for II and IV due to the p-nitrophenyl group presence.|
Substituents on N3 in an equatorial position facilitates the chain/ring tautomerization in 1,3-oxazines. Electron-with drawing groups increase the cyclic form contribution, whereas electron-releasing substituents would favor the chain forms.
|5.-||Comparison of the calculated and experimental values of chemical shifts and coupling constants indicates that the theoretical level used in the calculations is good enough to apply it in other cases to be studied.|
|6.-||The instability of TS for the axial ℑ equatorial I and II species suggests a strong repulsive effect between C5-CH3(ax) y N3-CH3 subtituents. A Brönstead coefficient calculation yield a ß = 0.564 for the equilibrium in I and ß = 0.547 for II, implying that in both cases the corresponding TSs would look the products alike. I and II cases follow well Hammond postulate.|
This work was partially supported by a grant (N° 203,021.019-1.0) from the Universidad de Concepción, Chile. S.M.H.O thanks CONICYT for a scholarship.
1.- D. Gürne, T. J. Urbanski, J. Chem. Soc. 1912, (1959) [ Links ]
2.- W. Hofman, L. Stepaniak, T. J. Urbanski, M. J. Witanowski, J. Am. Chem. Soc. 86, 554, (1964) [ Links ]
3.- J. M. Lehn, P. Linscheid, F. G. Riddell, Bull. Soc. Chim., France, 3, (1964) 1172,(1964) [ Links ]
4.- F. G. Riddell, J. M. Lehn, J. Chem. Soc. (B) 1224, (1968). [ Links ]
5.- Y. Allingham, R. C. Cookson, T. A. Crabb, S. Vary, Tetrahedron 24, 4625 (1968). [ Links ]
6.- D. Gürne, T. J. Urbanski, L. Stepamak, M. J. Witanowski, Tetrahedron Suplemento 6, 211, (1963). [ Links ]
7.- A. R. Katrizky, V. J. Baker, F. M. S. Brito-Palma, J. Chem. Soc. (Perkin II) 1739, (1980). [ Links ]
8.- M. Tafazzoli, C. Suarez, N. S.Trae, C. B. LeMaster, C. L. LeMaster, J. Phys. Chem. 96, 10201, (1992). [ Links ]
9.- C. L. Perrin, J. D. Thoburn, S. Elsheimer, J, Org. Chem. 56, 7034, (1991). [ Links ]
10.- M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, J. A. Montgomery, Jr., T. Vreven, K. N. Kudin, J. C. Burant, J. M. Millam, S. S. Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. Cossi, G. Scalmani, N. Rega, G. A. Petersson, H. Nakatsuji, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, M. Klene, X. Li, J. E. Knox, H. P. Hratchian, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, P. Y. Ayala, K. Morokuma, G. A. Voth, P. Salvador, J. J. Dannenberg, V. G. Zakrzewski, S. Dapprich, A. D. Daniels, M. C. Strain, O. Farkas, D. K. Malick, A. D. Rabuck, K. Raghavachari, J. B. Foresman, J. V. Ortiz, Q. Cui, A. G. Baboul, S. Clifford, J. Cioslowski, B. B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, I. Komaromi, R. L. Martin, D. J. Fox, T. Keith, M. A. Al-Laham, C. Y. Peng, A. Nanayakkara, M. Challacombe, P. M. W. Gilí, B. Johnson, W. Chen, M. W. Wong, C. González and J. A. Pople. Gaussian 03, Revisión D.01,Gaussian, Inc., Wallingford CT, 2004. [ Links ]
11.- J. G. Contaras, S. M. Hurtado, M. M. C. Ferreira, L. Gerli and R. Castillo. J.Chil. Chem. Soc. 50, 731, (2005). [ Links ]
12.- K. Wolisky, J. F. Hilton, P. Pulay, J. Am. Chem. Soc, 112, 8251, (1990) [ Links ]
13.- J. R. Cheeseman, G. W. Trucks, T. A. Keith, M . J. Frisch, J. Chem. Phys., 104, 5497, (1996). [ Links ]
14.- J. G. Contreras, L. A. Gerli, J. Chil. Chem, Soc, 52, 1271, (2007) [ Links ]
15.-A. E. Reed, R. B. Weinstock, F. Weinhold, J. Chem. Phys., 83, 735 (1985). [ Links ]
16.- A. E. Reed, L. A. Curtis, F. Weinhold, Chem. Rev. 88, 899, (1988). [ Links ]
17.- I. V. Alabugin, T.A. Zeidan, J. Am. Chem. Soc, 124, 3175, (2002). [ Links ]
18.-(a)Cope, A. C; Hancock, E. M. J. Am. Chem. Soc. 64, 1503, (1966), [ Links ] (b) Bergmann E.D., Gil-Av E., Pinchas S., J. Am. Chem Soc. 75, 358, (1953) [ Links ] (c) Fülöp F., Pihlaja K., Mattinen J., Bernáth G., J. Org. Chem. 52, 3821, (1987) [ Links ] (d) Fülöp F., Pihlaja K., Mattinen I, Neuvonen K., Bernáth G, Argay G., Calman A., J. Org. Chem. 58, 1967, (1993) [ Links ] (e) Neuvonen K., Fülöp F., Neuvonen H, Koch A., Kleinpeter E., Pihlaja K., J. Org. Chem 66, 4132, (2001) [ Links ] (f) Szatmári I., Martinek T. A., Lázár L., Koch A., Kleinpeter E., Neuvonen K, Fülöp F., J. Org. Chem 69, 3645, (2004). [ Links ]
21.- (a) S. Miertus, E. Scrocco and J. Tomasi., J. Chem Phys. 55, 117, (1981) [ Links ] (b) J. Tomasi, R. Bonaccorsi, R. Cammi and F. J. Olivares del Valle., J. Molec. Struct. Theochem. 234, 401, (1991). [ Links ]
22.- M. Hurtado, J. G. Contreras, A. Matamala, O. Mo, M. Yañez. New J. Chem. 32, 2209 (2008). [ Links ]
(Received 1 October 2008 - Accepted 16 January 2008)