Timothy J.
Dudley
*,
Jennifer E.
Beck
,
Earl E. P.
Santos
,
Kathryn A.
Johnston
,
William S.
Kassel
,
William G.
Dougherty
,
Walter J.
Boyko
and
Deanna L.
Zubris
*
Department of Chemistry, Villanova University, Villanova, PA 19085, USA. E-mail: deanna.zubris@villanova.edu; timothy.dudley@villanova.edu
First published on 2nd May 2012
Sterically hindered (imino)pyridine 2-{(2,6-Me2-C6H3)NC(i-Pr)}C5H4N (1) was synthesized via addition of isolated imidoyl chloride to an in situ lithiated pyridine. Room temperature 1-D and 2-D NMR spectroscopy reveals two rapidly equilibrating isomers in solution. Interconversion of these two isomers was verified by 2D-EXSY NMR spectroscopy. Calculations at the B3LYP and MP2 levels of theory reveal four relevant isomers, with two atropisomers of E geometry (1-EA and 1-EB) and two atropisomers of Z geometry (1-ZA and 1-ZB). A simple carbon–carbon bond rotation to alter the orientation of the isopropyl group provides a fifth, related conformer, 1-ZB′, that is the most stable species at the MP2 level. The transition states for E/Z isomerization and the isomerization pathways between atropisomers have been characterized. Comparison of experimental and ab initio NMR chemical shifts in combination with NOE analysis suggests that isomers 1-EB and 1-ZB/1-ZB′ are the dominant species in our solution phase NMR studies. Our understanding of the isomerization behavior of 1 will help inform the future design of readily complexed, sterically hindered mono(imine) and bis(imine) ligands.
![]() | ||
| Fig. 1 Previously reported bis(imino)pyridine (A,4C,6 and D6), bis(imino)pyrimidine (B5), and mono(imino)pyridine ligands (E8 and F9) with evidence for the existence of the Z,Z- or Z-isomer, respectively. | ||
The two limiting mechanisms of E/Z isomerization for imines are inversion at nitrogen (lateral shift) or rotation about the C
N double bond.10–12 For N-aryl imines, inversion at nitrogen is the generally accepted mechanism.13 Identification of both the E and Z isomers for N-aryl imines has been made by X-ray crystallography14 and by 1H NMR spectroscopy, with NOE experiments used to assist with conformational assignments,15 and variable temperature NMR used to measure the free energies of activation for the E/Z isomerization.15,16 Calculations suggest that C
N rotation would have a free energy of activation that is twice the magnitude of that for inversion at nitrogen, with the latter typically falling in a range from 12–29 kcal mol−1.16,17and cited references
For select imines, atropisomerism18 is invoked in combination with E/Z isomerism with experimental evidence derived from NMR spectroscopy16 and X-ray crystallography.14 In other cases, atropisomerism is not invoked, presumably since no experimental evidence exists to suggest its relevance.15,17,19 Atropisomerism was implied by Small and Brookhart20 and Cámpora et al.21 for free asymmetric bis(imino)pyridine ligands and for select corresponding Fe(II) complexes; 13C NMR spectroscopy provided evidence of hindered rotation about the imino N-aryl C bonds. It is notable that these authors do not invoke E/Z isomerization as an explanation for the fluxional behavior observed by NMR spectroscopy; presumably, the primary alkyl substituents on the imino carbon positions (methyl20,21 or butyl21) render the E,E-isomer energetically favorable. This is in contrast to studies by Orrell et al. where a symmetric bis(imino)pyridine ligand serves a bidentate ligand for Pd(II) or Pt(II),22 Re(I)23,24 and Pt(IV),25 and E/Z isomerization of the non-coordinated imine arm is invoked in combination with atropisomerism (again, with methyl substituents on the imino carbon positions). Finally, an account from 1987 with the closest structural similarity to our work describes theoretical identification of two atropisomers of the E-isomer of mono(imino)pyridine 2-(HN
CH)C5H4N.26
Herein, we describe the preparation and conformational analysis of the sterically hindered mono(imino)pyridine, 2-{(2,6-Me2-C6H3)NC(i-Pr)}C5H4N (1). For 1, our calculations reveal two atropisomers of E geometry and two atropisomers of Z geometry that are energetically viable. We have obtained experimental evidence of facile isomerization via 1-D and 2-D NMR spectroscopy. Importantly, ligand 1 is readily complexed with Ni(II) (vide infra), providing further support for E/Z isomerization. We believe our studies will help direct future development of sterically hindered bidentate and tridentate ligands containing the imine functionality.
N infrared stretching frequency of 1651 cm−1 for 1 is consistent with our structural assignment.
![]() | ||
| Scheme 1 Reagents and conditions: (i) 2,6-dimethylaniline, NEt3, CH2Cl2, reflux, 1 h ; (ii) PCl5, toluene, 22 °C, 15 h ; (iii) n-BuLi, CH2Cl2, −78 °C, 15 min ; (iv) −78–22 °C, 21 h. | ||
![]() | ||
| Fig. 2 Compound 1 atom labeling scheme. | ||
![]() | ||
| Fig. 3 Upfield region of 1H NMR spectra of 1 at 20 °C (top), −40 °C (middle), and −60 °C (bottom), respectively, in methylene chloride-d2 solvent. | ||
Assignments of 1 at low temperature were simplified due to general weak coupling giving “first order” resonances. The 6-position of the pyridine ring (H6b) was the starting point for the assignment, because this position has a smaller ortho coupling to the 5-position (H5b) than the other ortho couplings in the pyridine ring.29 A gCOSY plot was especially useful for resolving the partial overlap of the minor isomer H5b and major isomer H4b protons (see ESI†). Final confirmation used an EXSY plot at 20 °C where each site pair for the two isomers gave an exchange crosspeak on the 2D map. A listing of 1H NMR chemical shift assignments and coupling constants at −60 °C is provided in Table 1, and the −40 °C and 20 °C 1H NMR spectral data are present in the ESI.†
| δ (ppm) | Major isomer | Minor isomer | J (Hz) | Major isomer | Minor isomer |
|---|---|---|---|---|---|
| H3,5-a | 6.86 | 7.05 | H3a, H4a | 7.59 | 7.38 |
| H4a | 6.75 | 6.90 | H3b, H4b | 8.01 | 7.89 |
| Mea | 1.91 | 2.02 | H3b, H5b | 1.16 | 1.00 |
| H3b | 6.67 | 8.04 | H3b, H6b | 0.98 | 0.85 |
| H4b | 7.39 | 7.81 | H4b, H5b | 7.59 | 7.76 |
| H5b | 7.17 | 7.36 | H4b, H6b | 1.75 | 1.76 |
| H6b | 8.57 | 8.60 | H5b, H6b | 4.80 | 4.94 |
| Hi-Pr | 3.51 | 2.64 | Hi-Pr, Mei-Pr | 6.81 | 7.00 |
| Mei-Pr | 1.22 | 1.26 | — | — | — |
For a given proton of 1, the chemical shift difference between the major and minor isomer at −60 °C is generally small (< 0.20 ppm), with three notable exceptions: Hi-Pr (0.87 ppm), H3b (1.37 ppm) and H4b (0.42 ppm). The magnitudes of these chemical shift differences are maintained across the range of temperatures examined.
In an attempt to rationalize the chemical shift differences for Hi-Pr, H3b and H4b for the major versus minor isomer of 1, NOE experiments were carried out at −60 °C using a degassed methylene chloride-d2 solution of 1. A 2D-NOESY pulse sequence was used, and NOE peaks were observed as negative cross peaks found on the two-dimensional map. The only apparent through-space NOE signal was observed for the Mea–H3b interaction for the major isomer, as shown in Fig. 4. One may consider the normalized volume integral data to support this through-space interaction. If the major isomer Mea–H3,5-a volume integral is set to 100, the minor isomer Mea–H3,5-a normalized volume integral has a value of 94 and the major isomer Mea–H3b normalized volume integral has a value of 49. Thus, qualitatively the two “a”-ring methyl groups (Mea) must make a close approach to the pyridine proton, H3b, in the major isomer.
![]() | ||
| Fig. 4 NOESY map of 1 at −60 °C in methylene chloride-d2. | ||
The involvement of H3b in this through-space interaction is notable, since this is one of the three protons with highly disparate chemical shift behavior for the major and minor isomers of 1. The structural implications of these results will be discussed in context with our computational analysis.
To determine the temperature dependence of the rate of exchange for the major and minor isomers of 1, we used dynamic NMR techniques. 1H NMR spectra were collected in methylene chloride-d2 solvent over a range of temperatures from 0 °C (273 K) to 30 °C (303 K). Line shape analysis was carried out by using the commercially available program WINDNMR. Line broadening of the Mea signals was modelled, and both k and ΔG‡ were calculated over the range of temperatures following the methods described in the recent microreview by Mazzanti and coworkers.30 At 10 °C, the value of k is 10.9 ± 0.1 s−1, corresponding to a ΔG‡ value of 15.3 ± 0.1 kcal mol−1. An isomer ratio of 1.42
:
1 (major
:
minor) was found at this temperature. Our experimental isomerization barrier is in line with a computed isomerization barrier, as described in the next section.
![]() | ||
| Fig. 5 Structures and relative energies associated with the isomerization pathway between 1-EA and 1-ZA. Geometries are optimized at the B3LYP/6-31G* level of theory and energies are determined using single point MP2/6-311+G** calculations. The relative energies are referenced to 1-ZB. | ||
![]() | ||
| Fig. 6 Structures and relative energies associated with the isomerization pathway between 1-EB and 1-ZB. Geometries are optimized at the B3LYP/6-31G* level of theory and energies are determined using single point MP2/6-311+G** calculations. The relative energies are referenced to 1-ZB. | ||
Though 1-EA resembles the structure of the ligand when it is complexed to a metal, it is the highest energy equilibrium structure presented (see Fig. 5). It is connected to its geometric isomer, 1-ZA, through a transition state (1-TSA) that corresponds to inversion at the imine nitrogen, as shown in Fig. 5. This inversion mechanism is consistent with other reports for the interconversion of geometric isomers of N-aryl substituted imines.13 While the 2,6-dimethylphenyl moiety (“a” ring, Fig. 2) changes orientation drastically between 1-EA and 1-ZA, the Nim–Cim–C2b–Npy dihedral angle changes little between structures (39.7° in 1-ZA and 38.0° in 1-EA). The most interesting features of transition state 1-TSA involve three atoms: Cim, Nim, and C1a. The angle between these atoms is nearly linear (175.3°), which is expected since the transition between 1-ZA and 1-EA is described as an inversion at Nim. Both the Cim–Nim and C1a–Nim bonds (see Table 2) are significantly shorter (between 0.03 and 0.06 Å) in 1-TSA than in either 1-EA or 1-ZA.
| 1-EA | 1-ZA | 1-TSA | |
|---|---|---|---|
| C1a–Nim | 1.405 | 1.414 | 1.350 |
| Nim–Cim | 1.275 | 1.276 | 1.248 |
| Cim–C2b | 1.511 | 1.506 | 1.521 |
| C2b–Npy | 1.345 | 1.344 | 1.345 |
Unlike structures 1-EA and 1-ZA, structures 1-EB and 1-ZB exhibit a significant change in the Nim–Cim–C2b–Npy dihedral angle (see Fig. 6). The atoms are nearly planar in 1-EB (∠Nim–Cim–C2b–Npy = 175.8°) while they deviate significantly from planarity in 1-ZB (∠Nim–Cim–C2b–Npy = −138.8°). The transition state (1-TSB) connecting 1-EB and 1-ZB resembles 1-EB in terms of the aforementioned dihedral (∠Nim–Cim–C2b–Npy = −174.7°), but differs from both 1-EB and 1-ZB in terms of the bonds involving the same atoms (see Table 3). The trends are similar to those observed in 1-TSA, where both bonds involving Nim in 1-TSB are significantly shorter (between 0.03 and 0.06 Å) than those in either intermediate structure (either 1-EB or 1-ZB). Intrinsic reaction coordinate (IRC) calculations suggest the transition from 1-EB to 1-ZB corresponds to inversion at Nim and a subsequent rotation around the Cim–C2b bond.
| 1-EB | 1-ZB | 1-TSB | |
|---|---|---|---|
| C1a–Nim | 1.411 | 1.411 | 1.350 |
| Nim–Cim | 1.283 | 1.281 | 1.253 |
| Cim–C2b | 1.506 | 1.509 | 1.516 |
| C2b–Npy | 1.345 | 1.346 | 1.344 |
The relative energies of the structures at the MP2/6-311+G** level of theory (all referenced to the 1-ZB form) are shown in Fig. 5 and Fig. 6. Structures 1-EB and 1-ZB are the lowest energy structures and are nearly identical in energy (a difference of 0.6 kcal mol−1), with the barrier to their interconversion being between 21–22 kcal mol−1. Fig. 7 shows a fifth structure (1-ZB′) that is the lowest energy structure related to these imines. It corresponds to a rotation around the Cim–Ci-Pr bond, with Hi-Pr pointed towards Npy (see Fig. 2 for atom labelling). While this structure is closely related to 1-ZB and is probably in equilibrium with 1-ZB, its structure has an important property that helps in the identification of the isomers observed in NMR experiments. We will represent an equilibrium mixture of 1-ZB and 1-ZB′ as 1-ZB/1-ZB′.
![]() | ||
| Fig. 7 Structure and relative energy of 1-ZB′. Geometry is optimized at the B3LYP/6-31G* level of theory and the relative energy is determined using single point MP2/6-311+G** calculations. The relative energy is referenced to 1-ZB. | ||
The computed barrier to interconversion of 1-EB and 1-ZB/1-ZB′ (21–22 kcal mol−1) shows reasonable agreement with the experimentally derived barrier of 15.3 kcal mol−1. It is relevant that the computed barrier is derived from gas-phase calculations and the experimental barrier is derived from NMR spectra acquired in methylene chloride-d2 solvent—this may account for some of the discrepancy in values.
The interconversion of atropisomers (1-EAversus1-EB, and 1-ZAversus1-ZB) was also probed. These relationships are defined as atropisomerism due to hindered rotation about the Cim–C2b single bond. It is anticipated that the energy required for interconversion of these atropisomers would be much lower than that required for interconversion of geometric isomers (passing through transition states 1-TSA and 1-TSB, respectively). This is indeed the case, where constrained optimization calculations reveal a minimal barrier to interconversion of 1-ZA and 1-ZB. Further, the high-energy 1-EA structure readily converts to the lower energy 1-EB structure without surpassing a measurable barrier. Vibrational energy considerations render 1-EA as a structure with more transition state character as opposed to that of a discrete intermediate. Despite this distinction, the 1-EA/1-EB interconversion remains relevant for complexation to a metal.
Fortunately, two such instances were observed for 1. The experimental 1H NMR chemical shift for H3b differs significantly between the two isomers (a 1.4 ppm difference at −60 °C, see Table 1), with the proton in the minor isomer having the downfield chemical shift. Structure 1-EB has a much larger chemical shift for H3b (> 1.1 ppm downfield) than all other energetically viable structures (Table 4), suggesting it corresponds to the minor isomer. While the pattern of shifts observed for the other pyridine protons appears to corroborate this conclusion, the differences between shifts between the major and minor isomers is not large enough to make the computational results meaningful (including H4b, which experimentally has a 0.42 ppm difference in chemical shift upon comparison of the major and minor isomers). However, the shift difference observed for Hi-Pr is relatively large (a 0.87 ppm difference at −60 °C), making a comparison between computational and experimental results more conclusive. The minor isomer has the upfield Hi-Pr chemical shift, which is also the case for 1-EB, thus corroborating that 1-EB is the minor isomer.
| H3b | H4b | H5b | H6b | Hi-Pr | |
|---|---|---|---|---|---|
| 1-ZB | 6.72 | 7.14 | 6.90 | 8.62 | 3.22 |
| 1-ZB′ | 6.86 | 7.15 | 6.93 | 8.66 | 4.33 |
| 1-ZA | 7.34 | 7.48 | 6.85 | 8.32 | 3.05 |
| 1-EB | 8.51 | 7.56 | 7.04 | 8.61 | 2.72 |
The most downfield predicted chemical shift for Hi-Pr is observed for the most stable species at the MP2 level, 1-ZB′. While this value seems too large compared to the experimental value observed for the major isomer (3.51 ppm at −60 °C), the actual value of the shift is likely a weighted average of the shifts for 1-ZB (3.22 ppm) and 1-ZB′ (4.33 ppm) since they are related through a simple rotation around the Cim–Ci-Pr bond and the structures are relatively close in energy (0.6 kcal mol−1 difference).
NMR spectra were recorded on a Varian Mercury 300 Plus spectrometer at 300 MHz (1H) and 75 MHz (13C) at temperatures over a 210–300 K range, depending on the purpose of the given experiment. All chemical shifts were referenced relative to the NMR solvent (either residual protio or 13C signals for the solvent peak(s)). Temperature calibrations were performed using the standard methanol sample and software provided by Varian (now Agilent). gCOSY and NOESY 2D-NMR were obtained using standard Varian (now Agilent) pulse sequences. NOESY relative distance 2D-NMR plots were obtained using the NOESY sequence on degassed samples at −40 °C and −60 °C with 1.2 s mixing times. EXSY 2D-NMR plots were obtained using the NOESY sequence with 100 to 200 ms mixing times on un-degassed samples, typically at 20–30 °C. The following abbreviations are used for NMR splitting patterns: cd (complex doublet), ct (complex triplet), cm (complex multiplet) and br s (broad singlet). Infrared spectra were recorded using a Perkin-Elmer Spectrum One FTIR System; samples were prepared by placing the compounds on a diamond attenuated total reflectance (ATR) plate in either solid or liquid form. Elemental analyses were performed at Atlantic Microlab, Inc. in Norcross, Georgia, USA.
O). 1H NMR (benzene-d6, 22 °C): δ 6.94 (3H, m, Ar), 5.71 (1H, br s, N–H), 2.06 (6H, s, CH3), 2.03 (1H, septet, J = 7.2 Hz, CH(CH3)2), 1.07 (6H, d, J = 7.2 Hz, CH(CH3)2). 13C{1H} NMR (benzene-d6, 22 °C): δ 174.3 (C
O), 135.8, 134.9, 128.2, 127.1, 35.5, 19.9, 18.4. Characterization data is consistent with prior findings.27
N). 1H NMR (benzene-d6, 22 °C): δ 6.97 (3H, m, Ar), 2.69 (1H, septet, J = 7.2 Hz, CH(CH3)2), 2.05 (6H, s, CH3), 1.14 (6H, d, J = 7.2 Hz, CH(CH3)2).
N). 1H NMR (methylene chloride-d2, −60 °C, major isomer): δ 8.57 (1H, ddd, J = 4.80 Hz, 1.75 Hz, 0.98 Hz, H6b), 7.39 (1H, m, ddd, J = 8.01 Hz, 7.59 Hz, 1.75 Hz, H4b), 7.17 (1H, ddd, J = 7.59 Hz, 4.80 Hz, 1.16 Hz, H5b), 6.86 (2H, d, J = 7.59 Hz, H3,5-a), 6.75 (1H, t, J = 7.59 Hz, H4a), 6.67 (1H, ddd, J = 8.01 Hz, 1.16 Hz, 0.98 Hz, H3b), 3.51 (1H, septet, J = 6.81 Hz, Hi-Pr), 1.91 (6H, s, Mea), 1.22 (6H, d, J = 6.81 Hz, Mei-Pr). 1H NMR (methylene chloride-d2, −60 °C, minor isomer): δ 8.60 (1H, ddd, J = 4.94 Hz, 1.76 Hz, 0.85 Hz, H6b), 8.04 (1H, m, ddd, J = 7.89 Hz, 7.38 Hz, 0.85 Hz, H3b), 7.81 (1H, ddd, J = 7.89 Hz, 7.76 Hz, 1.76 Hz, H4b), 7.36 (1H, ddd, J = 7.76 Hz, 4.94 Hz, 1.00 Hz, H5b), 7.05 (2H, d, J = 7.38 Hz, H3,5-a), 6.90 (1H, t, J = 7.38 Hz, H4a), 2.64 (1H, septet, J = 7.00 Hz, Hi-Pr), 2.02 (6H, s, Mea), 1.26 (6H, d, J = 7.00 Hz, Mei-Pr). See Fig. 2 for an atom labeling scheme for compound 1 and see the ESI† for additional 1H NMR data of 1 in methylene chloride-d2 at −40 °C and 20 °C, respectively.
Footnotes |
| † Electronic Supplementary Information (ESI) available: NMR spectral characterization data for 1. Cartesian coordinates (in Angstroms) and absolute energies for all structures reported. CIF file for the Ni(II) complex‡ of ligand 1. See DOI: 10.1039/c2ra20688a/ |
| ‡ Single crystals of ([2-{(2,6-Me2-C6H3)NC(i-Pr)}C5H4N]NiBr2)2 were obtained from slow evaporation of an anhydrous dichloromethane solution at 22 °C. A single crystal was mounted using Paratone oil onto a glass fiber and transferred to the cold gas stream of the diffractometer. Crystal data. C17H20Br2N2Ni, M = 470.88, monoclinic, a = 12.1735(8), b = 10.4451(7), c = 14.3322(9) Å, α = 90, β = 100.5480(10), γ = 90°, U = 1793.3(2) Å3, T = 100 K, space group P2(1)/n, Z = 4, 17630 reflections measured, 6803 unique (Rint = 0.0362) which were used in all calculations. The final wR(F2) was 0.0606 (all data). CCDC reference number 840092. |
| This journal is © The Royal Society of Chemistry 2012 |