Acid-Catalyzed Rearrangements of 3-Aryloxirane-2-Carboxamides: Novel DFT Mechanistic Insights.

Efficient synthesis of 3-arylquinolin-2(1H)-ones and N-(2-carboxyaryl)-oxalamides from protic acid-catalyzed rearrangements of 3-aryloxirane-2-carboxamides was achieved recently but not well understood. In contrast to the classical Meinwald rearrangement, extensive DFT calculations reveal that the proximal aryl and amide groups have strong synergetic effects to control the amide-aided and aryl-directed oxirane-opening and further rearrangement sequences. The ortho-nitro substituent of the proximal aryl is directly involved in a nucleophilic oxirane ring-opening, the amide C=O is an important proton shuttle for facile H-shifts, while the N-aryl may act as a potential ring-closing site via Friedel-Crafts alkylation. The mechanistic insights are useful for rational design of novel synthesis by changing the aryl and amide functional groups proximal to the oxirane ring.

Oxiranes (or epoxides) containing a saturated CCO three‐membered ring are one of the most versatile classes of organic compounds available to the synthetic chemist which are frequently used in atom‐economical rearrangement reactions.1 In classical (Lewis or protic) acid catalyzed Meinwald rearrangements (Scheme 1),2 mono‐functional n class="Chemical">oxiranes are converted into neutral carbonyl compounds (aldehyde or ketone) via acid‐induced C−O cleavage followed by 1,2‐shift of a hydride or alkyl group within the resultant carbenium intermediate.3 Recently, a series of intermolecular tandem Meinwald rearrangement reactions have also been developed.4

Scheme 1

Classical Meinwald rearrangement of monofunctional oxiranes (top) and acid‐catalyzed rearrangements of 3‐aryloxirane‐2‐amides 1  a and 1  b for the respective syntheses of 3‐phenylquinolin‐2(1H)‐one 2  a and N‐(2‐carboxyphenyl)‐oxalamide 3  b. Strong synergetic effects of the proximal amide and aryl functional groups are revealed by our extensive DFT calculations, which are useful for rational design of novel synthetic approaches.

Classical Meinwald rearrangement of monofunctional oxiranes (top) and acid‐catalyzed rearrangements of 3‐n class="Chemical">aryloxirane‐2‐amides 1  a and 1  b for the respective syntheses of 3‐phenylquinolin‐2(1H)‐one 2  a and N‐(2‐carboxyphenyl)‐oxalamide 3  b. Strong synergetic effects of the proximal amide and aryl functional groups are revealed by our extensive DFT calculations, which are useful for rational design of novel synthetic approaches.

Alternatively, functional groups can be introduced as substituent to the oxirane ring to realize novel conversions via either synergetic or sequential rearrangements.5 Recently, very interesting sulfuric acid (n class="Chemical">H2SO4) catalyzed rearrangements of 3‐aryloxirane‐2‐amides in acetic acid (CH3COOH) solution were achieved for very efficient one‐pot synthesis of 3‐arylquinolin‐2(1H)‐ones5c and N‐(2‐carboxyaryl)‐oxalamides5b (Scheme 1), tentatively assumed to be initialized by the classical Meinwald rearrangement (forming neutral carbonyl intermediates), and followed by intramolecular electrophilic Friedel‐Crafts alkylation of the amide N‐aryl or nucleophilic addition to the nitro group.5b, 5c Quinolin‐2‐ones and oxalamides are omnipresent in naturally occurring and synthetic compounds displaying a broad range of pharmacological activities and practical applications.6 It is still desirable to gain deep mechanistic insights into such novel synthetic approach, especially the role of the amide and 2‐nitroaryl functional groups proximal to the central oxirane ring, which can be very helpful for further rational design of more efficient and novel synthesis of new compounds.

In this theoretical work, extensive DFT calculations at the PW6B95‐D3+COSMO‐RS //TPSS‐D3+COSMO level in acetic acid solution (see below for computational details) are performed to explore the potential free‐energy paths for the protic n class="Chemical">H2SO4 catalyzed rearrangement reactions of two typical 3‐aryloxirane‐2‐carboxamides, i. e., trans‐N,3‐diphenyloxirane‐2‐carboxamide (1  a) and trans‐3‐(2‐nitrophenyl)oxirane‐2‐carboxamide (1  b), which selectively lead to quite different 3‐phenylquinolin‐2(1H)‐one (2  a) and N‐(2‐carboxyphenyl)‐oxalamide (3  b) products, respectively (Scheme 1). In contrast to the classical acid‐catalyzed Meinwald rearrangement of mono‐functional oxiranes leading to neutral carbonyl compound intermediates, our DFT calculations clearly show that the proximal amide and aryl functional groups do show strong synergetic effects that profoundly affect the complete protic acid‐catalyzed rearrangement process.

Our DFT calculations show that in acetic acid solution, the neutral H‐bonded complex and the separated ion pair of n class="Chemical">HSO4 ‐ and CH3C(OH)2 + are 2.5 and 6.9 kcal/mol less stable than the solvated H2SO4 and CH3COOH molecules, respectively, and thus are thermodynamically unfavourable. Moreover, the proton transfers from H2SO4 to the nitro, oxirane and amide oxygen sites of 1  b are 19.0, 21.9 and 0.1 kcal/mol endergonic to form the corresponding protonated species along with the HSO4 ‐ anion. The selective and reversible protonation of 1  b at the amide C=O site by H2SO4 is almost barrierless. Similarly, the selective protonation of 1  a at the amide C=O site by H2SO4 is 1.1 kcal/mol endergonic over a rather low barrier of 2.7 kcal/mol, indicating a slightly reduced proton affinity of N‐aryl substituted amide group. Moreover, our DFT calculations suggest that most ionic and neutral species should exist as monomeric rather than hydrogen‐bonded complexes due to the high dissolving power of acetic acid as solvent, which makes our DFT mechanistic study in such polar and protic solvent simpler than expected.

The mechanism of H2SO4‐catalyzed rearrangement of 1  a is explored by our DFT calculations at first. As shown in Figure 1, the protonation of the n class="Chemical">amide C=O group of 1  a by H2SO4 is only 1.1 kcal/mol endergonic to form the separated ions of A + and HSO4 ‐ in solution. Due to the protonation‐enhanced C=N double bond character, trans‐to‐cis conversion of the amide group in A + is prevented by a high barrier of 29.7 kcal/mol to form Ac + that is −0.8 kcal/mol more stable; for comparison, similar conversion within neutral 1  a is 3.5 kcal/mol endergonic over a moderate barrier of 17.7 kcal/mol. Further proton transfer from the amide C=O to the oxirane oxygen within A may selectively cleave the phenyl‐connected C−O bond, which is 8.1 kcal/mol endergonic over a rather low barrier of only 11.0 kcal/mol (via transition structure TSB) to form the phenyl stabilized carbenium B. The recently proposed5c cleavage of the C−O bond proximal to the electron‐deficient amide group is prevented by a high barrier of 29.4 kcal/mol (via TSA), and thus is kinetically less favourable. The regioselective oxirane ring‐opening is thus strongly directed by the proximal π‐electron‐donating phenyl group,3 aided by the proximal amide C=O group as proton shuttle. The subsequent trans‐to‐cis conversion of the amide group of B (via TSC over a low barrier of 12.4 kcal/mol) can bring the amide N‐phenyl group more close to the benzylium‐like cation center, facilitating the intramolecular Friedel‐Crafts alkylation of the N‐phenyl at the ortho‐site possible (via TSD) to form the arenium cation D. Alternatively, the reactive B can be efficiently trapped by HSO4 ‐ to form the neutral complex Bs, followed by similar trans‐to‐cis conversion of the amide group (via TSCs) into Cs that may also lead to C after HSO4 ‐ elimination, which is kinetically more likely with an overall barrier of 21.9 kcal/mol and may avoid potential side reactions at reactive carbenium center. Subsequently, the new OH group of D may abstract the arenium proton from the same side of the C5N‐ring, which is −6.4 kcal/mol exergonic over a low barrier of 7.5 kcal/mol (via TSE) to form the complex E with a bound water (H2O) molecule. The new H2O molecule is then replaced by the adjacent phenyl group from the opposite ring side, which is −12.8 kcal/mol exergonic over a low barrier of 8.4 kcal/mol (via TSF) to form the site‐3‐carbon‐protonated product F. The final proton transfer from F to HSO4 − is −25.7 kcal/mol exergonic and almost barrierless to form the neutral 3‐phenyl‐substituted 2  a that can be easily protonated by H2SO4 at the C=O oxygen site. The relative trans‐configuration between the OH and Ph groups is crucial for the formation of 3‐phenyl‐substituted 2  a; 4‐phenylquinolin‐2(1H)‐one is expected for the cis‐configuration due to facile 1,2‐H‐shift instead. The strong synergetic effects of the proximal aryl and amide functional groups play a crucial role in directing the oxirane‐ring opening as well as the facile intramolecular Friedel‐Crafts alkylation. Such mechanism is quite different from the classical Meinwald rearrangement mechanism that involves a simple 1,2‐shift of a hydride or alkyl group after oxirane opening to form carbonyl compounds.

Figure 1

DFT computed free energy paths (in kcal/mol, at 298 K and 1 M concentration) for the H2SO4‐catalyzed rearrangement of trans‐N,3‐diphenyloxirane‐2‐carboxamide (1  a) to 3‐phenylquinolin‐2(1H)‐one (2  a). The free energies of neutral species (trapped or deprotonated by HSO4 ‐) are shown in red. The H, C, N and O atoms in ball‐stick models are shown as white, grey, blue and red balls. Partially breaking bonds are indicated by dashed lines.

DFT computed free energy paths (in kcal/mol, at 298 K and 1 M concentration) for the H2SO4‐catalyzed rearrangement of trans‐N,3‐diphenyln class="Chemical">oxirane‐2‐carboxamide (1  a) to 3‐phenylquinolin‐2(1H)‐one (2  a). The free energies of neutral species (trapped or deprotonated by HSO4 ‐) are shown in red. The H, C, N and O atoms in ball‐stick models are shown as white, grey, blue and red balls. Partially breaking bonds are indicated by dashed lines.

As shown in Figure 2, the protonation of the amide C=O group of 1  b by n class="Chemical">H2SO4 is only 0.1 kcal/mol endergonic to form the separated ions of G + and HSO4 ‐ in solution. The ortho‐nitro NO2 substituent to the proximal phenyl actually leads to a different mechanism of oxirane‐ring opening: the amide‐to‐oxirane proton transfer is now combined in a concerted way with the regioselective SN2 nucleophilic attack of one nitro oxygen toward the oxirane carbon connected to the phenyl, which is −6.9 kcal/mol exergonic over a low barrier of 17.6 kcal/mol (via transition structure TSH) to form the hetero‐cyclic intermediate H with a new C3NO five‐membered ring upon oxirane‐ring opening. In contrast, similar nitro attack at the amide‐connected oxirane carbon (via TSHa, see ESI) is prevented by a 11.1 kcal/mol higher barrier, while oxirane opening without nitro attack (via TSG) is kinetically 4.0 kcal/mol less favourable.

Figure 2

DFT computed free energy paths (in kcal/mol, at 298 K and 1 M concentration) for the H2SO4‐catalyzed rearrangement of trans‐3‐(2‐nitrophenyl)oxirane‐2‐carboxamide (1  b). The free energies of neutral species after deprotonation by HSO4 ‐ are shown in red. The H, C, N and O atoms in ball‐stick models are shown as white, grey, blue and red balls. Partially breaking bonds are indicated by dashed lines.

DFT computed free energy paths (in kcal/mol, at 298 K and 1 M concentration) for the H2SO4‐catalyzed rearrangement of trans‐3‐(2‐n class="Chemical">nitrophenyl)oxirane‐2‐carboxamide (1  b). The free energies of neutral species after deprotonation by HSO4 ‐ are shown in red. The H, C, N and O atoms in ball‐stick models are shown as white, grey, blue and red balls. Partially breaking bonds are indicated by dashed lines.

In the present case of the proximal 2‐nitrophenyl group, the extra proton on the n class="Chemical">C3NO ring of intermediate H can be easily transferred to the adjacent amide C=O group (via TSI) over a low barrier of 18.6 kcal/mol, which in turn induces a more facile cleavage of the endocyclic N−O bond (via TSJ) and subsequent nitroso rotation (via TSJc) to form the low‐lying intermediate Jc. In this way, one oxygen atom is transferred from the nitro group to oxidize the phenyl‐connected oxirane carbon. Further nucleophilic attack of the nitroso oxygen aided by the protonation of the new C=O group enable the facile 1,3‐shift of the cationic CHOHCONH2 + unit, which is −22.0 kcal/mol exergonic over a low free energy barrier of 10.9 kcal/mol (via TSK). The facile 1,2‐H‐shift along the C−N bond and concerted N−O bond cleavage (via TSL) eventually leads to the most stable form of the protonated product L, which is −55.5 kcal/mol exergonic over a barrier of 15.8 kcal/mol. The final deprotonation of L by the anion HSO4 − is −4.5 kcal/mol exergonic and almost barrierless to form the neutral product 3  b along with regenerated H2SO4 catalyst. Again, the aryl and amide functional groups proximal to the central oxirane ring lead to a new acid‐catalyzed rearrangement mechanism.

Interestingly, when trans‐3‐(2‐nitrophenyl)‐N‐phenyl‐n class="Chemical">oxirane‐2‐carboxamide (1  c) instead of 1  b was used as the substrate for which competitive rearrangements to both 3‐(2‐nitrophenyl)quinolin‐2(1H)‐one (2  c) and N‐(2‐carboxyphenyl)‐N‐phenyl‐oxalamide (3  c) become possible, only the oxalamide product 3  c was formed exclusively.5b Fully consistent with such experimental results, the DFT‐computed overall barriers for the catalytic rearrangements of very similar substrates 1  a and 1  c to the quinoline‐2‐one 2  a and the oxalamide 3  b are 21.9 and 18.6 kcal/mol, respectively, with the oxalamide product channel being kinetically 3 kcal/mol more favourable.

Our DFT mechanistic insights into the detailed role of various functional groups proximal to the central oxirane‐ring are useful for further rational design of novel synthesis of annulated compounds, for example, by replacing the ortho‐n class="Chemical">NO2 and amide NH2 groups with other unsaturated functional groups such as NO, aryl and acyl groups. Such idea is further supported by our test DFT calculations showing that similar nucleophilic oxirane opening reactions can also be induced by simple NO, CHO and CH=CH2 groups at the ortho‐site of the proximal phenyl group (see ESI). Interestingly, it was experimentally known that the replacement of the amide NH2 of 1  b with a phenyl group can lead to the same H2SO4‐catalyzed rearrangement into 2‐(2‐oxo‐2‐phenylacetamido)benzoic acid,5b though an additional O2‐elimination channel to form 2‐phenyl‐3‐hydroxyquinolin‐4‐one was also found very recently.7

In conclusion, detailed mechanisms of protic acid‐catalyzed rearrangement reactions of multifunctional 3‐n class="Chemical">aryloxirane‐2‐carboxamides for the synthesis of 3‐arylquinolin‐2(1H)‐ones and N‐(2‐carboxyaryl)‐oxalamides are revealed by state‐of‐the‐art DFT calculations. In contrast to monofunctional Meinwald rearrangement, the amide and aryl functional groups proximal to the central oxirane ring lead to novel rearrangement sequences, with the amide C=O, amide N‐aryl, aryl proximal to oxirane ring, and ortho‐nitro functional groups acting as proton shuttle, potential Friedel‐Crafts alkylation site, regioselective oxirane opening director, and nucleophile for SN2‐type oxirane opening, respectively. The mechanistic insights are useful for further rational design of novel synthesis by modifying the proximal functional groups.

Computational Methods

All DFT calculations are performed with the TURBOMOLE 7.3 suite of programs.8 The structures are fully optimized at the TPSS‐D3/def2‐TZVP+COSMO (HCOOH) level, which combines the TPSS meta‐GGA density functional9 with the BJ‐damped DFT−D3 dispersion correction10 and the def2‐TZVP basis set,11 using the Conductor‐like Screening Model (COSMO)12 for acetic acid solvent (dielectric constant ϵ=6.19 and diameter Rsolv=2.83 Å). The density‐fitting RI−J approach13 is used to accelerate the calculations. The optimized structures are characterized by frequency analysis (no imaginary frequency for true minima and only one imaginary frequency for transition states) to provide thermal free‐energy corrections (at 298.15 K and 1 atm) according to the modified ideal gas‐rigid rotor‐harmonic oscillator model.14

More accurate solvation free energies in acetic acid are computed with the COSMO‐RS model15 (parameter file: BP_TZVP_C30_1601.ctd) using the COSMOtherm package16 based on the TPSS‐D3 optimized structures, corrected by +1.89 kcal/mol to account for the 1 mol/L reference concentration in solution. To check the effects of the chosen DFT functional on the reaction energies and barriers, single‐point calculations at both TPSS‐D39 and hybrid‐meta‐GGA PW6B95‐D317 levels are performed using the larger def2‐QZVP11 basis set. Final reaction free energies (ΔG) are determined from the electronic single‐point energies plus TPSS‐D3 thermal corrections and COSMO‐RS solvation free energies. As expected, the overall results from both DFT functionals are in good mutual agreement (0.2±3.2 kcal/mol, mean difference between two functionals ± standard deviation) for reaction energies but somewhat higher reaction barriers (2.8±3.6 kcal/mol) are found at the PW6B95‐D3 level compared to the TPSS‐D3 results. In our discussion, the more reliable PW6B95‐D3+COSMO‐RS free energies (in kcal/mol, at 298.15 K and 1 mol/L concentration) are used unless specified otherwise. The applied DFT methods in combination with the large AO basis set provide usually accurate electronic energies leading to errors for chemical energies (including barriers) on the order of typically 1–2 kcal/mol. This has been tested thoroughly for the huge data base GMTKN5518 which is the common standard in the field of DFT benchmarking.

Conflict of interest

The authors declare no conflict of interest.

As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.

Supplementary

Click here for additional data file.