Abstract

The determination of ring size can vary from simple to complex, but the process in difficult cases can be advantageously augmented by DFT calculation of NMR parameters such as the chemical shifts of (), (), and other nuclei as well as pertinent spin-spin (scalar) coupling constants, for example, those between protons (). Differentiation between 6- and 7-membered ring formation in the case of 3,4-dihydro-2H-3-hydroxymethyl-1,4-benzoxazine and 2,3,4,5-tetrahydro-1,5-benzoxazepine-3-ol was evaluated with a view to not only affecting 6- versus 7-membered ring differentiation generally for cases on hand, but also in the case of literature reports where the assigned structures may be in doubt. Thus, the main focus was on the usually reported NMR parameters of , , and and wherein the analysis was found to be highly successful, particularly for , and thus potentially amenable for broad application.

1. Introduction

Obviously the determination of ring size in the framework of organic compounds is a crucial and basic precept in structural chemistry, whether the compounds in question arose through synthetic processes or were obtained from natural sources. The distinction between, for example, 6- and 7-membered rings can, in particular cases, be surprisingly difficult even by the application of sophisticated NMR experiments [1]. Evermore so in heterocyclic systems where the availability of suitable hydrogen atoms within or proximal to the cyclic unit for analysis may be limited—whether concerning either solely 1H nuclei or 1H nuclei in conjunction with other nuclei, for example, 13C, the problems and challenges of discerning between 5- and 6-membered rings in heterocyclic systems with a deficient number of 1H nuclei by NMR have been considered [1] previously in depth using DFT calculations based on the consideration of various NMR parameters as depicted in Figure 1. For example, the chemical shifts of nuclei such as 1H (), 13C (), and 15N (); various spin-spin (scalar) coupling constants between nuclei such as 1H and 1H , 1H and 13C (), and 1H and 15N (); NOEs between nuclei such as 1H and 1H () and 1H and 13C (). The prediction of NMR parameters [2, 3] such as chemical shifts, but especially coupling constants, has come of age recently in terms of development and practical application [46], and especially with respect to conformational analysis [79].

We herein complement the previous study of 5- versus 6-membered ring formation with the consideration of 6- and 7-membered ring differentiation using as a basis for the study the pair of compounds, 3,4-dihydro-2H-3-hydroxymethyl-1,4-benzoxazine (1) and 2,3,4,5-tetrahydro-1,5-benzoxazepine-3-ol (2), presented in Figure 2. The chosen pair is notable and of interest for a number of reasons: the assigned structures of 1 have been erroneously reported (vide infra) which has come to light only very recently [12]; the compounds are important synthons [12]; they are interesting models for structural studies, in particular, conformational analysis [13]; furthermore, they and their analogs have had significant biological activities reported [12, 14]; finally, they belong to pharmacologically important compound classes [12].

In addition to solving this somewhat persistent problem as an example of the wider application of DFT-calculated NMR parameters such as , , and for differentiating between 6- and 7-membered rings, the question arises, what if workers do not have recourse to performing the NMR measurements? This may not necessarily be because workers do not have access to suitable instrumentation, but because they do not have the compounds on hand, either at present or in the foreseeable future. Workers may, akin to the case of García-Rubiño et al. [12], be contemplating a synthesis, or, alternatively, workers may not even be intent on physically procuring the compounds as they may be undergoing structure-activity relationships (SARs) or docking studies in the absence of physical measurements altogether and be hoping to use the ensuing results to explain past observations or predict future directions for work. Hence, in this examination we focus primarily on the basic NMR data which is generally available in the literature, namely, the chemical shifts of 1H () and 13C nuclei () and the 1H-1H spin-spin (scalar) coupling constants () even if the structural proofs of the compounds may be highly amenable, even facile, using modern NMR techniques on sophisticated, modern instrumentation. Thus, in order to save the time, effort, and cost of resynthesizing compounds or otherwise procuring them (not to mention avoiding confusion or even the reporting of erroneous results), the current study has been undertaken to examine the particular example case of 6- and 7-membered ring differentiation.

2. Computational Methods

DFT quantum chemical calculations were performed using Gaussian09 [15] (version A.01) and analyzed using GaussView (version 3.07). Geometry optimization of the structures in the gas phase was performed using the M06-2X hybrid metadensity functional [16, 17] with the 6-31G(d) basis set in tandem with vibrational analysis and thermochemistry calculations at the same level of theory. It is worth noting that the use of a modest basis set for structural optimization affords computational efficacy where high structural fidelity is not required to obtain quality results for the subsequent calculation of NMR parameters [4]. Vibrational analyses, invoking the keyword freq  =  noraman, were conducted to confirm that optimized structures were true minima on the potential energy surface by not providing imaginary frequencies and to obtain the thermodynamic contributions at 298.15 K and 1 atm wherein frequencies were left unscaled. If necessitated by the presence of imaginary frequencies, structures were reoptimized (together also with vibrational analyses and thermochemistry calculations) using tight convergence criteria by invoking the keywords opt  =  tight and int  =  ultrafine.

Absolute chemical shieldings and spin-spin coupling constants were calculated for geometry-optimized structures in the gas phase using the GIAO method [18] and the B3LYP functional [19, 20] with the cc-pVTZ basis set. Chemical shifts relative to the reference compounds were determined by subtracting calculated shieldings from the calculated shieldings of the reference compounds: TMS for which both and = 0.00 ppm and CH3NO2 for which = 0.00 ppm and = 602.50 ppm. The chemical shifts of  1H, 13C, and 15N nuclei were calibrated following literature methodology [21] using the following relationships: The calculated values were not calibrated. Spin-spin (scalar) coupling constants such as and were calculated by invoking the keyword with the calculations incorporating all four contributing terms—the Fermi contact term, paramagnetic and diamagnetic spin-orbit and spin-dipolar terms—despite the increased computational cost of doing so. Coupling constant values obtained from the calculations were used as received and were not corrected.

3. Results and Discussion

3.1. Course of the Reaction

The example pair of compounds on which this analysis has been conducted was taken from the recent literature wherein the structure of 3,4-dihydro-2H-3-hydroxymethyl-1,4-benzoxazine (1) had been assigned to the products resulting from several different reported reactions [14, 22, 23]. In each case the product was the same despite the varying conditions and different research groups concerned. The correct structure of the resultant compound, in fact, has been shown very recently to be 2,3,4,5-tetrahydro-1,5-benzoxazepine-3-ol (2) by García-Rubiño et al. [12] who proved conclusively the structures of compounds 1 and 2 by both the application of NMR methods and synthetic procedures. The problem here is that structures were not originally assigned based on factual NMR measurements and workers inadvertently and unwittingly implied that multiplicity assignments of the 13C had been determined, presumably by the application of some multiplicity determining pulse sequence like the DEPT pulse sequence. It is unclear how the various workers settled on the assigned structures (assumed reaction course perhaps?), but subsequently it was demonstrated rather simply the proof [12] of the structures of 1 and 2 by readily available means to most chemists, namely, by DEPT and HMBC NMR experiments. Despite the ready and straightforward means to affect distinction, the erroneous structural assignments have, disparagingly, perpetuated [22, 23].

The reaction at hand, in various guises, is essentially composed of intramolecular nucleophilic attack by an amino group on a monosubstituted 1,2-epoxide as depicted in Scheme 1. By Baldwin’s [24] rules for ring closure, both routes are feasible processes, 6-exo-tet to lead to 1 and 7-exo-tet to lead to 2. The reaction is interesting as it compares 1° versus 2° attack to yield 7- or 6-membered rings, respectively. Thus the result is a confluence of thermodynamic constraints (a 6- or a 7-membered ring with the former generally preferred and borne out by calculations in this instance) versus kinetic factors (1° or 2° attack) versus proximity of the reacting centers. Thus, it is not necessarily a surprise that both compounds could be produced in comparable amounts depending on the reaction conditions applied [12]. Interestingly, the original source of the error may have come from work on a 6-aza analogue of 1, 4-acetyl-3-(hydroxymethyl)-3,4-dihydro-2H-pyrido[3,2-b]oxazine, prepared in a similar fashion to the preparation reported by Jiao et al. [14] and where too the workers [25] not only seem to have made the same inadvertent error regarding the multiplicities of the carbons, but have also made an unsound presumption regarding the course of the reaction. In that case, a 2° amide supposedly displaces an oxygen from a mono-substituted 1,2-epoxide group, not at the 1° carbon, but instead at the 2° carbon, to form the product. Unfortunately, there must now be grave doubt regarding the correctness of that structure [25] as well.

697156.sch.001

3.2. Conformational Analysis and Molecular Modeling

For compound 1, two half-chair conformations are possible for a 6-membered alicyclic fused to a benzene ring [8] with the –CH2OH group adopting either axial or equatorial orientations. Additionally, the NH proton can potentially also adopt either axial or equatorial orientations, and finally, rotation about the CH–CH2OH (i.e., C3–C4) and CH2–OH bonds is also possible (3 rotamers each). Altogether then, 2 × 2 × 3 × 3 = 36 conformations are conceivable. The former two mobilities were fully explored but for the latter two bond rotations, judicious selection based on hydrogen bonding opportunities and steric congestion considerations, as well as also taking into account the preference for a gauche disposition [26] between N-5 and O-4 (the “gauche effect”), was used to limit the number of structures examined. Altogether, 6 low-energy conformers () were located for 1 (Table 1); that is, 5 other conformers were within 2 kcal mol−1 of the global minimum energy conformer. Whilst the global minimum incorporated a hydrogen bond and a –CH2OH group that was equatorially oriented, of the other 5 low-energy conformers each with an axially oriented –CH2OH group, 2 conformers involved intramolecular hydrogen bonding of HO-4 whilst another 2 involved HN-5. The remaining conformer lacked intramolecular hydrogen bonding altogether. In 5 conformers, the NH proton was equatorially oriented whilst in the remaining conformer it was isoclinal. The two lowest-energy structures for compound 1 are presented in Figures 3 and 4.

Other high-energy structures were not included in the ensuing calculation of NMR parameters since their contribution is negligible (well less than 5%) and thus only conformers within 2 kcal mol−1 of the global minimum energy conformer were considered for further analysis, and in this instance, the global minimum energy conformer accounts for ca. 65% of the population. In line with expectations, intramolecularly hydrogen bonding structures were generally lower in energy, although surprisingly, one nonhydrogen bonded structure did rank amongst these low-energy structures. Even more surprising, it did not have a gauche interaction either. Of note, it is well known that intramolecularly hydrogen-bonded structures can persist in solution, even in CDCl3. Indeed, the observation of a 2H broad singlet—presumably the NH and OH protons—by Jiao et al. [14] is an indication of this for these compounds since if intermolecularly hydrogen-bonded structures only were present, then the NH and OH protons would likely be so broad as to be unobservable, as it is generally the case for CDCl3 solutions of compounds with labile protons.

7-Membered rings present a challenge in terms of conformational analysis which is generally substantially greater in comparison to 6-membered rings [13] due to additional conformational flexibility and a larger number of conformers that reside close in energy. For a 7-membered ring fused to an aromatic ring, chair (C) and twist-boat (TB) conformers are acknowledged [13] as being the preferred conformers, with the former much more so but this can be altered substantially depending on the particular case of a substituent at C-5 [13]. For 7-membered ring 2, the search protocol followed similarly to that for the 6-membered ring 1 and consisted of starting with one hand-built structure, for example, a chair, optimizing the structure and then using the result as a template for the construction of subsequent structures. Likewise, the first permutation considered was a change of the OH group to the other isoclinal/axial orientation by simply switching the OH and H moieties on C-3. Then, for each of these two structures, axial and equatorial orientations of the NH group were evaluated. Since the –CHOH segment is part of the ring, rotations about the C2–CHOH and CHOH–C4 bonds to yield various gauche interactions and hydrogen bonding possibilities were not possible without changing the ring conformation. Incidentally, for 2, two gauche interactions are possible, even simultaneously, from the disposition of the OH group to the N-5 and O-1 atoms in comparison to 1. The same process was repeated starting from a twist-boat conformation. Of these constructs, possibilities for hydrogen bonding (and alternatively, not) were examined by rotation about the C3–OH bond. Boat conformations were also examined but comparable low-energy structures were not located. By this means, a total of three low-energy conformations, one chair (), and two twist boats (), were found (Table 1) with the chair conformation, as expected, overwhelmingly dominant (>90% populated). Interestingly, in this conformation, a three-center hydrogen bond is formed and this conformation is further stabilized by two gauche interactions. The other two conformations each incorporated one hydrogen bond and one gauche interaction. Again, likewise to 1, in all these three conformations, the NH proton was equatorially oriented. The two lowest-energy structures for compound 2 are presented in Figures 5 and 6.

It must be said that, somewhat paradoxically, the 7-membered ring compound 2 was more acquiescent to conformational analysis than the 6-membered ring compound 1 in this case.

3.3. Analysis of the Calculated NMR Parameters

The of the labile protons HO-4 and HN-5 were not evaluated since they are usually highly variable in nature as they are very dependent on concentration, temperature, solvetic hydrogen bonding, and other environmental influences, as well as dynamic exchange effects yielding shifted or averaged chemical shifts. Moreover, since it is difficult to account for intermolecular hydrogen bonding, the CH–CH2OH and CH2–OH bond rotations in compound 1 and the CH2–OH bond rotation in compound 2 are not considered very reliable in terms of conformational analysis. Fortunately, the rotations were not found to be influential on the geometry of the –O–CH2–CH– fragment in 1 and thus, the –O–CH2–CH– spin system in 1 was the primary focus for the analysis whilst the full –O–CH2–CH–CH2– spin system was considered in 2. The calculated for all the non-labile protons together with the experimental [12] are presented in Table 2 wherein it can be seen that the calculated in relation to the observed are in fair agreement but not strongly indicative of the assigned structure. With both compounds on hand, a valued judgment could easily be made, for example, and and particularly with respect to their difference, but decidedly less if the experimental values of only one compound were available. For example, with only compound 1 available, an erroneous assignment could easily be made based on . This is not too surprising as are known to be problematic in relation to which are known to be reliable indicators [6, 9] of structure. Unexpectedly, and are much better predicted in 1 than in 2 despite the above concerns. The aromatic provide a reasonable indication and though they are not unequivocal, they nevertheless do insinuate correct assignment, and moreover, could do so even if only one compound was available.

The calculated for the carbon nuclei together with the experimental [12] are presented in Table 3 wherein it can be seen that the calculated are in excellent agreement with the observed and would convincingly enable the correct structural assignment to be made—even in the event of only one isomer being available. The , and results are particularly emphatic, but a similar assessment also holds for the aromatic carbons that are not in common with the heterocycle (i.e., , and ). Of particular note, the structural determination can readily be done even without the assignment of the carbons experimentally.

The calculated for the protons of the alicyclic ring together with their experimental analogous [12] are presented in Table 4. The calculated vicinal couplings and in relation to the observed values are in excellent agreement and would convincingly enable the correct structural assignment to be accomplished, again even if only one compound was available. The situation, however, with the and couplings is slightly more complicated, the values for 2 are in excellent agreement with the experimental values but those for 1, in concert with the aforementioned concerns, are not. Nevertheless, with due appreciation of the factors concerned, correct analysis can be effected if either both compounds were on hand or even if only one was available. The geminal couplings, and , were disappointing and do not provide any sort of indication of the correct structure. The long-range couplings were all evaluated as negligible, consistent with observations, except for in compound 2 which was calculated to be 1.5 Hz in magnitude. This is in fine agreement with the measured value of 1.4 Hz and arises from the ω-disposition of these two protons in the dominant chair conformer (90%).

García-Rubiño et al. [12] also reported that whilst correlations from the H-2s to C-9a were observed in the HMBC spectra for both compounds 1 and 2 (i.e., the must be sizeable, e.g., a few Hz at least in magnitude for each compound), only for 2 could correlations from protons of the heterocyclic ring be observed to C-5a, and these were from the H-4s, thus unequivocally proving its structure. This is because the in 2 extending over three bonds were expected to be sizeable whereas the in 1 extending over four bonds were not. From the calculations (Table 5), these surmises—which are considered in general as truisms—are well supported, but unexpectedly not by the population-weighted averaged (PWA) values which are small for all of , , , and in compound 2 and in compound 1 despite clear correlations being observed in each case. Evidently the conformer interconversion must be relatively slow on the NMR timescale to allow the development of antiphase magnetization in the HMBC spectrum during the initial delay time (typically ca. 60 ms if a of 8 Hz is desired for optimum observation) for creation of multiple quantum magnetization due to the large values present in the minor conformers. Thus, the correlations observed [12], with the exception of H-2β-C-9a, in fact arise from the minor conformers (ca. 24% and ca. 10% for 1 and 2, resp.). Enigmatically then, acquiring the HMBC spectrum at lower temperature could lead to a reduction in peak intensity due to a lowering of the content of the minor conformers whilst acquiring the HMBC spectrum at higher temperature could also lead to a reduction in peak intensity due to fast interconversion of the conformers leading to a PWA value being in effect. Of note though, large values were neither calculated for nor for in any of the conformations for compound 1.

Although García-Rubiño et al. [12] did not measure 15N NMR spectra, they are also likely to be indicative of the structures obtained based on the calculated values of (Table 5) where the chemical shift difference for between compounds 1 and 2 was found to be a sizeable 8.6 ppm. The 15N NMR spectra for compounds 1 and 2 are likely to be easily measureable, either directly using techniques such as DEPT and INEPT, [27], or indirectly using HSQC since N-5 bears a directly attached proton.

Similarly, though García-Rubiño et al. [12] did not measure 17O NMR spectra, differences were calculated for the two oxygen nuclei present (Table 5). For , the calculated chemical shift difference was actually only marginal, 8.3 ppm, given the breadth of the chemical shift range for 17O and the difficulty in measuring accurately due to the large breadth of the 17O NMR signals (of the order of several ppm depending on , etc.). However, for , the calculated chemical shift difference was a sizeable 23.4 ppm, in line with expectations [28, 29] for differences between primary and secondary alcohols. Moreover, the chemical shift difference between O-1 and O-3/4 contrasts nicely in the two compounds, 80.4 ppm for compound 1 and 48.7 ppm for compound 2. Thus, as a result of the chemical shift difference of the two oxygen nuclei present in the one molecule, 17O NMR may afford a greater degree of structural reliability over 15N NMR in this instance. Hence, despite the problems with acquiring 17O NMR data-low sensitivity arising from low γ, low isotopic abundance, and broad lines thus requiring relatively large sample amounts and/or long acquisition times in addition to the aforementioned accurate determination of —it could well be decisive and worth the effort.

4. Conclusions

In this work, the differentiation between 6- and 7-membered ring formation in the case of 3,4-dihydro-2H-3-hydroxymethyl-1,4-benzoxazine (1) and 2,3,4,5-tetrahydro-1,5-benzoxazepine-3-ol (2) was evaluated using DFT calculation of the NMR parameters , , and . The point was not only to confirm the correction of the erroneous structural assignments in the literature with respect to these particular compounds, but also with a view to affecting 6- versus 7-membered ring differentiation generally for cases on hand and to add to our methodology for differentiating between 5- and 6-membered ring formation. Since the distinction was readily facilitated in this instance using and as well as , the methodology can thus not only be applied generally, but can also be applied to other cases in the literature where assigned structures may be in doubt. Thus, the structural assignments can be either confirmed or dispelled a priori without wasting valuable time and resources. This is important in cases not only for workers contemplating or undergoing synthetic processes, but also others who may be involved with prospective compounds but not in physical possession of the compounds themselves and be using them for SAR or docking studies in the absence of physical measurements or for literature evaluations. In this particular evaluation, it was found that were not so reliable, the provided mixed results but were reliable if due care was applied, and finally, the were found to be extremely indicative of the correct structure at hand, even in the absence of experimental data for all candidate structures.

It is also worth bringing to attention several points regarding the obtainment of reliable calculations. Firstly, it is highly important that a systematic search of potential contributing conformational participants be undertaken. Secondly, that sufficient geometry optimization using a reasonable level of theory [4, 5] is utilized. In the present day with the computational power available to most practitioners, this can be a quite modest requirement. Thirdly, an appropriate basis set for the NMR parameter calculations [4, 5] is also utilized; a large number of conformations need not necessarily be calculated though, just the primary ones accounting for the majority of the populated states. And finally, appropriate treatment of the participating conformers by evaluating population-weighted averages [7, 9] is affected when appropriate.

Acknowledgment

Thanks are due to The CSC-IT Center for Science Ltd. for providing computational resources.